Phân tích dao động tấm dày bằng phương pháp phần tử tự do galerkin (EFGM)

Đối với các bài toán cơ học môi trường liên tục, việc chia lưới là cả một vấn đề, ta có thể gặp khó khăn khi đưa ra một mắt lưới thích hợp do cấu trúc hình học phức tạp, đôi khi quá trình xử lý có thể yêu cầu chia lại lưới thường xuyên. Thời gian cần thiết để ta tạo các mắt lưới hoặc chia lại lưới thường gấp nhiều lần so với thời gian cần thiết hình thành và giải quyết các hệ phương trình. Để tránh hoặc giảm công việc tạo lưới khó khăn này, người ta đã đề xuất một số phương pháp, thường được gọi chung là phương pháp không lưới, chẳng hạn: phương pháp SPH (SPHM), phương pháp RKP (RKPM), phương pháp EFG (EFGM), phương pháp PUFE (PUFEM)…
Nội dung trích xuất từ tài liệu:
Phân tích dao động tấm dày bằng phương pháp phần tử tự do galerkin (EFGM) TẠP CHÍ KHOA HỌC, Đại học Huế, Tập 74B, Số 5 (2012), 75-80 PHÂN TÍCH DAO ĐỘNG TẤM DÀY BẰNG PHƯƠNG PHÁP PHẦN TỬ TỰ DO GALERKIN (EFGM) Tôn Thất Hoàng Lân Khoa Xây Dựng, Trường Đại học Kiến Trúc Tp. HCM Tóm tắt. Đối với các bài toán cơ học môi trường liên tục, việc chia lưới là cả một vấn đề, ta có thể gặp khó khăn khi đưa ra một mắt lưới thích hợp do cấu trúc hình học phức tạp, đôi khi quá trình xử lý có thể yêu cầu chia lại lưới thường xuyên. Thời gian cần thiết để ta tạo các mắt lưới hoặc chia lại lưới thường gấp nhiều lần so với thời gian cần thiết hình thành và giải quyết các hệ phương trình. Để tránh hoặc giảm công việc tạo lưới khó khăn này, người ta đã đề xuất một số phương pháp, thường được gọi chung là phương pháp không lưới, chẳng hạn: phương pháp SPH (SPHM), phương pháp RKP (RKPM), phương pháp EFG (EFGM), phương pháp PUFE (PUFEM)… Phương pháp không lưới với các cách tiếp cận khác nhau của nó đã được áp dụng cho một loạt các vấn đề kỹ thuật. Chẳng hạn, G.R.Liu sử dụng phương pháp EFGM để phân tích vỏ và tấm mỏng. Việc áp dụng phương pháp không lưới phần tử tự do Galerkin (EFGM) để nghiên cứu dao động tấm dày chữ nhật đồng chất là chủ đề bài báo. Từ khoá: phương pháp không lưới, phương pháp EFGM, dao động, bình phương tối thiểu động, tấm dày. 1. Hàm xấp xỉ bình phương tối thiểu động (MLS) Hàm MLS đã được phát triển bởi Lancaster và Salkauskas để xấp xỉ đường cong và bề mặt. Xem xét miền Ω có chứa một tập hợp các nút phân tán xi(1 ≤ i ≤ n), tương ứng giả định các giá trị wi. Xấp xỉ MLS của hàm liên tục w trên Ω gọi là wh(x) được cho bởi: m w h (x)   P (x)α (x)  p T (x)αx  (1) i i1 i trong đó p(x) là một tổ hợp m hàm độc lập tuyến tính, p T (x)   p (x)  1 p (x) .... p (x)  2 m  (2) và α(x) là một tổ hợp các thông số chưa xác định,  α T (x)  α 0 (x) α1 (x) α 2 (x) .... α m (x)  (3) Các thông số α(x) được tìm thấy tại điểm x bất kỳ bằng cách cực tiểu: 75 76 Phân tích dao động tấm dày bằng phương pháp phần tử tự do Galerkin (EFGM) n(x) 2 J(x)    (x  x ). w h (x )  w i i i i1 i   (4) n(x) 2 J(x)    (x  x ). p T (x ). (x) - w i i i i1 i   (5)  i (x  x i ) là hàm trọng số có giá trị khác 0 trên miền ảnh hưởng của nút xi. Chỉ có các nút xi mà miền ảnh hưởng chứa điểm x sẽ xuất hiện trong công thức trên. Kích thước miền ảnh hưởng của mỗi nút và cách lựa chọn hàm trọng số là yếu tố quyết định sự gần đúng của MLS. Cực tiểu J (x) để biết các thông số α(x) α(x)  A 1 (x).B(x).w (6)  B  1 (x  x1).p(x1 )  2 (x  x 2 ).p(x 2 )... n (x  x n ).p(x n ) n(x) A    (x - x ).p T (x ).p(x ) i i i i1 i  (7) (8)  Ta thay kết quả và suy ra: w T  w1 w 2 ...... w n n w h (x)   φ (x).w  Φ(x).w i i 1 i  (9) (10) m Với: φ (x)   p (x)(A 1(x)B(x))  p T A 1B i ji i j1 j (11) Φ(x)  φ1 (x)  (12) φ 2 (x) ... φ n (x)   2. Lý thuyết biến dạng cắt bậc nhất (Lý thuyết tấm dày Mindlin-Reissner) Đối với tấm dày, pháp tuyến của mặt trung hòa vẫn thẳng nhưng không nhất thiết phải vuông góc với mặt phẳng này trong suốt quá trình biến dạng. Lý thuyết của Mindlin-Reissner được sử dụng để tính toán biến dạng và ứng suất trong tấm dày trong khi lý thuyết Kirchhoff-Love được áp dụng cho các tấm mỏng. Lý thuyết của Mindlin giả định rằng có một sự biến đổi tuyến tính của chuyển vị thông qua độ dày tấm nhưng độ dày tấm không thay đổi trong suốt quá trình biến dạng. Điều này ngụ ý rằng, ứng suất pháp thông qua chiều dày tấm bị bỏ qua, trạng thái ứng suất phẳng xảy ra. Mặt khác, lý thuyết của Reissner giả định rằng ứng suất uốn là tuyến tính trong khi ứng suất cắt là bậc hai trên suốt chiều dày của tấm. Điều này dẫn đến một trạng thái khác mà ở đó chuyển vị trên suốt chiều dày tấm không nhất thiết phải tuyến tính và độ dày tấm có thể thay đổi trong quá trình biến dạng. Lý thuyết tấm dày Mindlin-Reissner là sự kết hợp lý thuyết của hai người và từ đây thường được gọi là lý thuyết biến dạng cắt bậc nhất. Trường chuyển vị được thể hiện theo tài liệu [10]. TÔN THẤT HOÀNG LÂN 77 u (x)  u o (x , x )  x  α α 1 2 3 α (13) u (x)  w(x , x ) 3 1 2 (14) α  1,2 (15) trong đó x1, x2 , x3, 1, 2 và w lần lượt được diễn đạt như hình vẽ và u o là α chuyển vị trong mặt phẳng trung hòa theo 2 phương 1, 2. Tấm có chiều dày h thì độ cứng uốn D  Eh 3 12(1  ν 2 ) (16) Hình 1. Hệ trục tọa độ và các thông số của tấm dày  Phương trình vi phân chủ đạo: D 2 2 w   q(x, t)  2hρw (17)  Đối với dao động tự do D 2 2 w  2hρw (18) 3. Phương pháp EFG và ví dụ số xét dao động của tấm hình chữ nhật Phương pháp EFG là một phương pháp kết hợp giữa xấp xỉ bình phương tối thiểu động MLS và dạng yếu Galerkin. Ở đây dạng yếu Galerkin thể hiện bởi phương trình T T T T  dΩ  0  δε σdΩ   δλ wdΓ w   λ δwdΓ w   ρδw w Ω Ω Ω Ω trong đó λ T  λ u (19) λ v  là toán tử Lagrange Xét dao động của tấm dày đồng chất hình chữ nhật có kích thước 4x4m, dày 0.4m. Nghiệm chính xác của tần số dao động tự do theo tài liệu [10] sẽ được so sánh với nghiệm có được từ phương pháp EFG và ta thể hiện sự sai số giữa chúng theo các bảng 1, 2, 3. 78 Phân tích dao động tấm dày bằng phương pháp phần tử tự do Galerkin (EFGM) Hình 2. Thể hiện 3 mode dao động đầu tiên Ngôn ngữ Matlab đã được sử dụng để lập trình tính toán dao động tấm dày dựa theo phương pháp không lưới phần tử tự do Galerkin (EFGM). Bảng 1. Sai số giữa nghiệm gần đúng và nghiệm chính xác của mode 1 Mode 1 55 1010 1515 2020  EFG   chínhxác  100  chínhxác 82% 61% 1.1% 0.2% Bảng 2. Sai số giữa nghiệm gần đúng và nghiệm chính xác của mode 2 Mode 2 55 1010 1515 2020  EFG   chínhxác  100  chínhxác 51% 47% 0.4% 0.05% TÔN THẤT HOÀNG LÂN 79 Bảng 3. Sai số giữa nghiệm gần đúng và nghiệm chính xác của mode 3 Mode 3 55 1010 1515 2020  EFG   chínhxác  100  chínhxác 8.7% 7.9% 0.33% 0.04% 4. Kết luận Bài báo này đã áp dụng phư ...