Luận văn phân tích kết cấu tấm bằng phần tử ES-MITC3 - Nguyễn Duy Quang

Luận văn thạc sĩ phân tích kết cấu tấm bằng phần tử biến dạng trơn ES-MITC3. Tài liệu nghiên cứu cơ học kết cấu chi tiết và giải pháp tính toán hiệu quả.

• Tác giả: Nguyễn Duy Quang
• Người hướng dẫn: TS. Châu Đình Thành
• Trường học: Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh
• Chuyên ngành: Kỹ thuật xây dựng công trình dân dụng và công nghiệp
• Đề tài: Phân tích kết cấu tấm bằng phần tử biến dạng trơn ES-MITC3
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2016
• Địa điểm: Thành phố Hồ Chí Minh

Tóm tắt

I. Tổng quan về luận văn phân tích kết cấu tấm bằng phần tử biến dạng trơn ES MITC3

Phân tích kết cấu tấm bằng phương pháp phần tử hữu hạn (FEM) đã trở thành công cụ quan trọng trong kỹ thuật xây dựng. Tuy nhiên, phương pháp truyền thống gặp hạn chế như sai số do hiện tượng khóa cứng (shear locking) khi tấm mỏng. Kỹ thuật biến dạng trơn (ES-FEM) kết hợp phần tử MITC3 ra đời nhằm khắc phục nhược điểm này. ES-MITC3 sử dụng biến dạng trơn trên cạnh phần tử, giảm thiểu sai số do khóa cứng và cải thiện độ chính xác kết quả. Phương pháp này đặc biệt hiệu quả trong phân tích tĩnh, dao động và ổn định kết cấu tấm. Nghiên cứu này kế thừa các thành quả từ Liu [17] và Phùng Văn Phúc [18], đề xuất ứng dụng ES-MITC3 trong phân tích tấm vuông ngàm chịu tải phân bố đều.

1.1. Phương pháp phần tử hữu hạn truyền thống trong phân tích tấm

Phương pháp phần tử hữu hạn (FEM) truyền thống sử dụng các phần tử như MITC3, MITC4 để mô hình hóa tấm. Các phần tử này dựa trên lý thuyết biến dạng cắt bậc nhất (FSDT) để tránh hiện tượng khóa cứng. Tuy nhiên, khi tấm mỏng (t/L < 0.05), FEM truyền thống thường cho kết quả không chính xác do sai số tích phân số. Đặc biệt, phần tử MITC3 có độ chính xác thấp hơn MITC4 nhưng đơn giản hơn trong lập trình. Việc lựa chọn phần tử phụ thuộc vào yêu cầu độ chính xác và tài nguyên tính toán.

1.2. Sự phát triển của kỹ thuật biến dạng trơn ES FEM

Kỹ thuật biến dạng trơn (ES-FEM) cải tiến từ FEM bằng cách làm trơn biến dạng trên cạnh phần tử. Phương pháp này giảm thiểu sai số do khóa cứng và cải thiện độ hội tụ. ES-FEM kết hợp với phần tử MITC3 (ES-MITC3) cho kết quả vượt trội so với FEM truyền thống. Ưu điểm chính bao gồm: (1) độ chính xác cao hơn 20-30% so với FEM-Q4 cùng số nút, (2) ổn định trong phân tích động do không có năng lượng ảo, (3) đơn giản hơn trong triển khai lập trình. Các nghiên cứu trước đây đã chứng minh hiệu quả của ES-FEM trong phân tích tĩnh và dao động.

II. Phân tích vấn đề và hạn chế của phương pháp truyền thống

Phương pháp FEM truyền thống gặp nhiều hạn chế khi phân tích tấm mỏng. Hiện tượng khóa cứng (shear locking) xảy ra khi tấm mỏng (t/L < 0.05), dẫn đến sai số lớn trong kết quả độ võng và moment. Phần tử MITC3 có ưu điểm đơn giản nhưng độ chính xác thấp hơn MITC4. Ngoài ra, FEM truyền thống không ổn định trong phân tích động do xuất hiện năng lượng ảo. Nghiên cứu của Liu [17] chỉ ra rằng ES-FEM khắc phục được những hạn chế này nhờ kỹ thuật làm trơn biến dạng. Tuy nhiên, ứng dụng ES-MITC3 trong thực tế vẫn còn hạn chế do thiếu dữ liệu so sánh và chưa được chuẩn hóa rộng rãi.

2.1. Hiện tượng khóa cứng trong phân tích tấm mỏng

Khóa cứng xảy ra khi tấm mỏng (t/L < 0.05) do biến dạng cắt bị đánh giá quá cao. Trong FEM truyền thống, phần tử MITC3 sử dụng phép xấp xỉ tuyến tính cho biến dạng cắt, dẫn đến sai số tích phân số. Kết quả là độ võng và moment bị đánh giá thấp so với lý thuyết. Ví dụ, với tấm vuông ngàm chịu tải phân bố đều, phần tử MITC3 cho độ võng thấp hơn 15% so với lý thuyết khi t/L = 0.001. ES-MITC3 khắc phục vấn đề này bằng cách làm trơn biến dạng trên cạnh phần tử.

2.2. So sánh hiệu suất giữa phần tử MITC3 và ES MITC3

Nghiên cứu so sánh hiệu suất giữa MITC3 và ES-MITC3 cho thấy ES-MITC3 vượt trội trong hầu hết các trường hợp. Với tấm vuông ngàm chịu tải phân bố đều, ES-MITC3 cho độ võng chính xác hơn 25% so với MITC3 khi t/L = 0.001. Đối với moment, ES-MITC3 cũng cho kết quả tốt hơn, đặc biệt trong phân tích tấm dày (t/L = 0.1). Ưu điểm này do ES-MITC3 giảm thiểu sai số do khóa cứng và cải thiện độ hội tụ. Tuy nhiên, ES-MITC3 đòi hỏi tài nguyên tính toán cao hơn do thuật toán phức tạp.

III. Phương pháp nghiên cứu và đề xuất giải pháp ES MITC3

Đề xuất phương pháp ES-MITC3 dựa trên kỹ thuật biến dạng trơn (ES-FEM) kết hợp phần tử MITC3. Phương pháp này sử dụng biến dạng trơn trên cạnh phần tử để giảm sai số do khóa cứng. Thuật toán bao gồm ba bước: (1) chia lưới tấm thành các phần tử tam giác 3-node, (2) tính toán biến dạng trơn trên cạnh phần tử, (3) tích hợp kết quả để thu được độ võng và moment. Phương pháp được triển khai bằng ngôn ngữ lập trình Python và kiểm chứng qua các ví dụ số. Kết quả cho thấy ES-MITC3 cho độ chính xác cao hơn 20% so với FEM truyền thống.

3.1. Cơ sở lý thuyết của phương pháp ES MITC3

ES-MITC3 dựa trên lý thuyết biến dạng cắt bậc nhất (FSDT) kết hợp kỹ thuật biến dạng trơn (ES-FEM). Biến dạng trơn được tính toán bằng cách làm mịn biến dạng trên cạnh phần tử, giảm thiểu sai số do khóa cứng. Phương pháp này kế thừa từ Liu [17], sử dụng phần tử tam giác 3-node (MITC3) để đơn giản hóa mô hình. Thuật toán ES-MITC3 bao gồm các bước: chia lưới, tính biến dạng trơn, tích hợp kết quả. Phương pháp này ổn định trong phân tích động do không có năng lượng ảo.

3.2. Quy trình triển khai thuật toán ES MITC3

Quy trình triển khai ES-MITC3 gồm bốn bước chính: (1) chia lưới tấm thành các phần tử tam giác 3-node, (2) tính toán ma trận độ cứng phần tử bằng phương pháp biến dạng trơn, (3) lắp ghép ma trận toàn cầu, (4) giải hệ phương trình tuyến tính để thu được độ võng và moment. Phương pháp này được triển khai bằng ngôn ngữ Python, sử dụng thư viện NumPy cho tính toán ma trận. Kết quả được kiểm chứng qua các ví dụ số, so sánh với lý thuyết và FEM truyền thống.

IV. Kết quả nghiên cứu và ứng dụng thực tiễn của ES MITC3

Nghiên cứu tiến hành phân tích tấm vuông ngàm chịu tải phân bố đều bằng phương pháp ES-MITC3. Kết quả cho thấy ES-MITC3 cho độ võng và moment chính xác hơn 20% so với FEM truyền thống. Đặc biệt, ES-MITC3 ổn định trong phân tích động do không có năng lượng ảo. Phương pháp này cũng đơn giản hơn trong triển khai lập trình. Ứng dụng ES-MITC3 trong thực tiễn giúp cải thiện độ chính xác kết cấu tấm, đặc biệt trong các dự án xây dựng đòi hỏi độ chính xác cao.

4.1. Kết quả phân tích tấm vuông ngàm bằng ES MITC3

Phân tích tấm vuông ngàm chịu tải phân bố đều bằng ES-MITC3 cho thấy độ võng và moment đạt độ chính xác cao. Ví dụ, với tấm mỏng (t/L = 0.001), ES-MITC3 cho độ võng chính xác hơn 25% so với MITC3. Đối với tấm dày (t/L = 0.1), ES-MITC3 cũng cho kết quả tốt hơn, đặc biệt trong phân tích moment. Kết quả này chứng minh hiệu quả của ES-MITC3 trong cải thiện độ chính xác kết cấu tấm.

4.2. Ứng dụng ES MITC3 trong kỹ thuật xây dựng

ES-MITC3 có tiềm năng ứng dụng rộng rãi trong kỹ thuật xây dựng, đặc biệt trong phân tích tấm sàn, tấm mái, và tấm vách. Phương pháp này giúp cải thiện độ chính xác kết cấu, giảm thiểu sai số do khóa cứng. Ngoài ra, ES-MITC3 ổn định trong phân tích động, phù hợp cho các công trình chịu tải trọng động như cầu, nhà cao tầng. Ứng dụng ES-MITC3 trong thực tiễn đòi hỏi tài nguyên tính toán cao hơn, nhưng đem lại lợi ích lâu dài về độ chính xác.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình dân dụng và công nghiệp, kết cấu tấm chiếm hơn 60% khối lượng chịu lực của hệ thống sàn, vách cứng, tường chắn và móng bè. Việc tính toán chính xác ứng xử cơ học của tấm có ý nghĩa quyết định đến mức độ an toàn và tối ưu hóa chi phí vật liệu của dự án. Nghiên cứu của tác giả Nguyễn Duy Quang, thực hiện năm 2016 tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh, tập trung giải quyết triệt để vấn đề "khóa cắt" (shear locking) trong phân tích tĩnh kết cấu tấm bằng phương pháp phần tử hữu hạn.

Lý thuyết tấm cổ điển Kirchhoff-Love bỏ qua hoàn toàn biến dạng trượt, dẫn đến sai số lớn khi áp dụng cho tấm dày có tỉ lệ bề dày trên nhịp $h/a \ge 1/20$. Ngược lại, lý thuyết biến dạng cắt bậc nhất Reissner-Mindlin giải quyết được tấm dày nhưng lại gặp hiện tượng khóa cắt nghiêm trọng khi chuyển sang vùng tấm mỏng ($h/a < 1/20$), khiến độ cứng của kết cấu bị phóng đại quá mức và kết quả chuyển vị bị triệt tiêu sai thực tế. Mục tiêu cốt lõi của luận văn là xây dựng và kiểm chứng phần tử tấm tam giác 3 nút mới, ký hiệu là ES-MITC3, kết hợp kỹ thuật làm trơn biến dạng trên cạnh (ES-FEM) với kỹ thuật nội suy các thành phần tensor hỗn hợp (MITC3).

Nghiên cứu được triển khai trong giai đoạn từ tháng 09/2014 đến tháng 04/2016, khảo sát trên 7 bài toán chuẩn đối sánh kinh điển với dải độ mảnh biến thiên từ $L/t = 10$ đến $L/t = 1000$. Kết quả tính toán cho thấy phần tử ES-MITC3 triệt tiêu hoàn toàn hiện tượng khóa cắt, đưa sai số chuyển vị võng và mô men uốn tại tâm tấm xuống mức cực thấp, chỉ từ 0,00% đến 0,96%, cải thiện từ 15% đến 25% độ chính xác so với phần tử MITC3 nguyên bản trên cùng mật độ chia lưới.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu dựa trên nền tảng lý thuyết tấm biến dạng cắt bậc nhất (First-order Shear Deformation Theory - FSDT), hay còn gọi là lý thuyết tấm Reissner-Mindlin. Lý thuyết này giả định đoạn thẳng vuông góc với mặt trung bình của tấm trước khi biến dạng vẫn duy trì dạng thẳng nhưng không nhất thiết phải vuông góc với mặt trung bình sau khi biến dạng. Hệ số hiệu chỉnh cắt được thiết lập ở giá trị chuẩn $5/6$ (xấp xỉ 0,833) nhằm hiệu chỉnh sự phân bố phi tuyến của ứng suất tiếp qua chiều dày tấm.

Khung phân tích tích hợp hai mô hình tính toán tiên tiến bao gồm:

1. Kỹ thuật nội suy biến dạng cắt trượt MITC3: Sử dụng các điểm liên kết (tying points) trên từng cạnh của phần tử tam giác 3 nút để giả định biến dạng cắt ngang là hằng số, ngăn ngừa suy biến ma trận độ cứng cắt $C_s$ khi chiều dày tấm tiến về 0.
2. Phương pháp làm trơn biến dạng dựa trên cạnh (Edge-based Smoothed Finite Element Method - ES-FEM): Phân chia miền tính toán của tấm thành các ô làm trơn $\Omega^{(k)}$ gắn liền với từng cạnh biên phần tử, giúp tính toán ma trận độ cứng $K^{(k)}$ thông qua tích phân biến dạng trơn thay cho đạo hàm hàm dạng thông thường.

Ba khái niệm học thuật trọng tâm chi phối mô hình gồm: trường chuyển vị 3 bậc tự do trên mỗi nút ($w, \theta_x, \theta_y$), độ cứng trụ uốn của tấm $D = \frac{Eh^3}{12(1-\nu^2)}$, và hiện tượng khóa cắt xảy ra khi năng lượng biến dạng cắt lấn át năng lượng uốn.

Phương pháp nghiên cứu

Nghiên cứu sử dụng phương pháp số mô phỏng kết cấu dựa trên mã nguồn tự lập trình hoàn toàn trên môi trường Matlab. Nguồn dữ liệu kiểm chứng bao gồm 7 bài toán chuẩn đối sánh (benchmark problems) có lời giải giải tích chính xác hoặc nghiệm số tin cậy cao từ các nhà cơ học nổi tiếng thế giới. Cỡ mẫu khảo sát bao gồm hệ thống lưới phần tử tam giác được làm mịn liên tục từ 2x2, 4x4, 8x8, 10x10, 12x12, 16x16 đến 32x32 phần tử (tương đương quy mô từ 8 đến 2048 phần tử cho toàn bộ kết cấu tấm).

Phương pháp chọn mẫu dựa trên tính đa dạng hình học và điều kiện biên đại diện cho công trình thực tế:

• Tấm hình vuông ngàm và tựa đơn 4 cạnh chịu tải trọng phân bố đều.
• Tấm hình thoi Razzaque và tấm hình thoi Morley chịu tải phân bố đều.
• Tấm hình thoi ngàm một cạnh chịu tải phân bố đều với các góc nghiêng 20 độ, 40 độ và 60 độ.
• Tấm hình tròn ngàm chịu tải trọng phân bố đều với các tỉ lệ độ dày $t/R = 0,02$ và $t/R = 0,2$.

Lý do lựa chọn phương pháp phân tích số ES-MITC3 là vì phương pháp này cho phép tự động hóa chia lưới với phần tử tam giác 3 nút đơn giản, thích ứng hoàn hảo với các hình học phức tạp mà không đòi hỏi thêm các bậc tự do bổ sung hay kỹ thuật hiệu chỉnh ma trận phức tạp. Toàn bộ quá trình tính toán và kiểm tra nghiệm thu được tiến hành liên tục trong suốt 20 tháng nghiên cứu.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng số và phân tích hội tụ trên 7 nhóm mô hình hình học đã mang lại 4 phát hiện quan trọng:

1. Triệt tiêu hoàn toàn hiện tượng khóa cắt trên tấm mỏng: Khi khảo sát tấm vuông mỏng ngàm 4 cạnh có tỉ lệ $t/L = 0,001$, giá trị độ võng chuẩn hóa tại tâm $w_c/(qL^4/100D)$ đạt 0,1267 ở mức chia lưới 16x16. Độ lệch so với lời giải chính xác chỉ là 0,39%, vượt trội so với mức sai số trên 15% của phần tử MITC3 tiêu chuẩn khi chưa làm trơn.
2. Độ chính xác tuyệt đối cho kết cấu tấm dày: Đối với bài toán tấm dày có tỉ lệ $t/L = 0,1$, phần tử ES-MITC3 cho giá trị độ võng chuẩn hóa đạt 0,1499, trùng khớp 100% với lời giải giải tích chính xác (độ lệch 0,00%). Đồng thời, mô men uốn chuẩn hóa $M_c/(qL^2/10)$ đạt 0,2310 với độ lệch cực thấp là 0,48%.
3. Tính miễn nhiễm cao với lưới phần tử bị méo (distorted mesh): Trong điều kiện chia lưới tam giác méo dạng ngẫu nhiên ở cấp lưới 16x16, phần tử ES-MITC3 vẫn kiểm soát độ lệch võng của tấm mỏng ở mức 0,87% và độ lệch võng tấm dày ở mức 0,00%, trong khi mô men uốn chỉ lệch 0,14%.
4. Sự đồng nhất về hướng chia lưới và đối xứng hình học: Việc phân chia lưới theo hướng chéo trái hay chéo phải đều đem lại kết quả chuyển vị và nội lực đồng nhất 100%. Sơ đồ tính toán trên 1/4 tấm đối xứng cho thấy giá trị mô men chuẩn hóa hội tụ nhanh hơn từ 10% đến 15% so với phân tích toàn phần.

Thảo luận kết quả

Nguyên nhân căn bản giúp phần tử ES-MITC3 vượt trội so với các phần tử FEM truyền thống nằm ở cơ chế làm mềm độ cứng (softening effect) của kỹ thuật làm trơn trên cạnh. Trong phương pháp phần tử hữu hạn tiêu chuẩn, mô hình luôn có xu hướng quá cứng (over-stiff) do các ràng buộc chuyển vị liên tục. Kỹ thuật ES-FEM tạo ra các ô làm trơn trung bình trên từng cạnh, làm giảm độ cứng giả tạo của ma trận $K$, từ đó nâng cao độ chính xác của trường ứng suất và mô men uốn.

So sánh với các nghiên cứu cùng thời kỳ như ES-DSG3, NS-DSG hay MIN3, phần tử ES-MITC3 thể hiện tốc độ hội tụ nhanh hơn từ 18% đến 22% trên cùng mật độ nút. Dữ liệu từ các bảng tính có thể được biểu diễn trực quan qua biểu đồ đường cong hội tụ năng lượng và biểu đồ quan hệ chuyển vị - tỉ lệ $L/t$. Trên đồ thị biểu diễn dải độ mảnh từ 10 đến 1000, đường cong của ES-MITC3 duy trì dạng nằm ngang ổn định, không xuất hiện hiện tượng gãy khúc hay suy giảm chuyển vị như phần tử tam giác thông thường. Điều này chứng minh tính ứng dụng bền vững của thuật toán cho mọi loại kết cấu từ vỏ mỏng kim loại đến sàn bê tông dày.

Đề xuất và khuyến nghị

Nhằm chuyển hóa các kết quả nghiên cứu lý thuyết vào thực tiễn kỹ thuật xây dựng và công nghệ mô phỏng, 4 giải pháp hành động cụ thể được đề xuất như sau:

1. Tích hợp thuật toán phần tử ES-MITC3 vào các gói phần mềm phân tích kết cấu thương mại (như ETABS, SAP2000 hoặc các module mã nguồn mở trong OpenSees) trong lộ trình 12 tháng, hướng tới mục tiêu giảm từ 20% đến 30% thời gian tính toán lặp khi phân tích sàn bê tông cốt thép khẩu độ lớn.
2. Xây dựng quy trình tự động hóa chia lưới tam giác thích nghi (adaptive triangular meshing) cho các kết cấu sàn có lỗ mở kỹ thuật hoặc hình học vát góc phức tạp, đảm bảo độ lệch chuyển vị kiểm soát dưới 1,0% trong giai đoạn thiết kế kỹ thuật từ quý 1 đến quý 4 năm tới do các kỹ sư kết cấu chủ trì.
3. Mở rộng khung tính toán ES-MITC3 sang bài toán phân tích động lực học, dao động tự do và tấm vật liệu composite phân lớp chức năng (FGM), triển khai bởi các nhóm nghiên cứu tại các trường đại học kỹ thuật trong thời gian 24 tháng với chỉ tiêu công bố ít nhất 2 bài báo quốc tế uy tín.
4. Tổ chức các chương trình tập huấn chuyển giao công nghệ phân tích phần tử hữu hạn trơn cho các công ty tư vấn thiết kế xây dựng, với mục tiêu đào tạo từ 50 đến 100 kỹ sư kết cấu hàng năm nhằm tối ưu hóa kích thước cấu kiện và tiết kiệm từ 5% đến 8% chi phí bê tông, cốt thép.

Đối tượng nên tham khảo luận văn

Tài liệu luận văn là nguồn tham khảo chuyên sâu và hữu ích cho 4 nhóm đối tượng trọng tâm:

1. Kỹ sư kết cấu công trình dân dụng và công nghiệp: Ứng dụng thuật toán để kiểm tra độ võng và phân phối mô men của sàn không dầm, sàn ứng suất trước và móng bè nhà cao tầng, đảm bảo độ tin cậy vượt mức 99% theo tiêu chuẩn thiết kế.
2. Học viên cao học và nghiên cứu sinh ngành Cơ học tính toán và Kỹ thuật xây dựng: Sử dụng hơn 60 trang công thức toán học chi tiết về FSDT, ES-FEM và MITC3 làm tài liệu nền tảng để phát triển các đề tài luận án chuyên sâu về phần tử tấm, vỏ.
3. Chuyên gia phát triển phần mềm CAE/CAD: Tiếp cận giải thuật lập trình phần tử tam giác 3 nút không tăng bậc tự do toàn cục để tối ưu hóa nhân tính toán (solver), giúp tiết kiệm từ 15% đến 25% dung lượng bộ nhớ RAM khi giải hệ phương trình lớn.
4. Giảng viên và nhà nghiên cứu tại các trường đại học kỹ thuật: Bổ sung kiến thức về các phương pháp phần tử hữu hạn nâng cao vào bài giảng chuyên đề cho sinh viên năm cuối và học viên sau đại học.

Câu hỏi thường gặp

1. Kỹ thuật ES-MITC3 giải quyết hiện tượng khóa cắt như thế nào?
   Thuật toán kết hợp phép nội suy biến dạng cắt trượt MITC3 trên các điểm liên kết cạnh với kỹ thuật tích phân làm trơn ES-FEM trên từng miền biên. Sự phối hợp này triệt tiêu năng lượng cắt ảo khi chiều dày tấm tiến về 0, giúp mô hình duy trì độ võng chuẩn xác 0,1267 trên dải tỉ lệ mỏng $t/L = 0,001$.

2. Phần tử ES-MITC3 có làm tăng số lượng bậc tự do toàn cục không?
   Không. Phần tử ES-MITC3 giữ nguyên cấu trúc phần tử tam giác 3 nút cơ bản với đúng 3 bậc tự do trên mỗi nút bao gồm chuyển vị thẳng đứng $w$ và hai góc xoay $\theta_x, \theta_y$. Do đó, số lượng bậc tự do toàn cục hoàn toàn tương đương với phần tử FEM cổ điển, không làm tăng thời gian giải hệ phương trình $Ku = F$.

3. Tốc độ hội tụ của ES-MITC3 có vượt trội hơn phần tử tứ giác 4 nút không?
   Có. Kết quả số đối chứng cho thấy ES-MITC3 đạt tốc độ hội tụ nhanh hơn từ 18% đến 22% so với phần tử MITC3 và tương đương, thậm chí tốt hơn phần tử tứ giác MITC4 trên cùng mật độ nút phân bố, đặc biệt là khi lưới phần tử chịu biến dạng méo.

4. Phần tử ES-MITC3 hoạt động ra sao khi lưới tam giác bị méo dạng?
   Phần tử thể hiện độ ổn định cơ học xuất sắc. Ngay cả khi các nút lưới bị xê dịch ngẫu nhiên tạo ra các tam giác méo, sai số chuẩn hóa của độ võng chỉ dao động trong khoảng từ 0,00% đến 0,87%, đảm bảo kết quả tính toán không bị suy giảm chất lượng trong thực tế chia lưới phức tạp.

5. Luận văn đã kiểm chứng thuật toán trên bao nhiêu dạng hình học kết cấu tấm?
   Nghiên cứu đã thực hiện kiểm chứng toàn diện trên 7 bài toán mẫu gồm bài toán kiểm tra đơn biến, tấm vuông ngàm, tấm vuông tựa đơn, tấm hình thoi Razzaque, tấm hình thoi Morley, tấm thoi ngàm một cạnh với góc nghiêng từ 20 đến 60 độ, và tấm tròn ngàm với các tỷ lệ độ dày khác nhau.

Kết luận

• Xây dựng thành công phần tử tấm tam giác 3 nút ES-MITC3 tích hợp lý thuyết biến dạng cắt bậc nhất và kỹ thuật làm trơn trên cạnh.
• Khử bỏ hoàn toàn hiện tượng khóa cắt, xử lý hiệu quả cho cả kết cấu tấm mỏng ($t/L = 0,001$) và tấm dày ($t/L = 0,1$).
• Đạt độ chính xác vượt trội với sai số chuyển vị võng dưới 0,39% và sai số mô men uốn dưới 0,48% trên lưới chuẩn 16x16.
• Duy trì tính ổn định cơ học cao trên hệ thống lưới tam giác méo dạng với sai số kiểm soát nghiêm ngặt dưới 0,87%.
• Cung cấp mã nguồn thuật toán Matlab chuẩn xác, tạo tiền đề vững chắc cho việc mở rộng sang các bài toán động lực học và vật liệu composite phi tuyến trong 1 đến 2 năm tới.

Các kỹ sư và nhà nghiên cứu kết cấu nên khai thác và ứng dụng ngay giải thuật ES-MITC3 để nâng cao độ chính xác và hiệu quả kinh tế trong thiết kế công trình xây dựng hiện đại.

Trích đoạn nội dung tài liệu

Chương 1 TỔNG QUAN 1. Kết cấu tấm là một trong những kết cấu phổ biến trong tự nhiên. Ngoài ra, còn có vô số kết cấu tấm do chính con người tạo ra. Do có đặc trưng mỏng, nhẹ, khả năng chịu uốn, vượt nhịp lớn, nên tấm đã được ứng dụng rất nhiều vào các lĩnh vực xây dựng (sàn, tường, tôn, vách), cơ khí (tấm cách nhiệt, tua bin, lò phản ứng), hàng không (thân vỏ máy bay),… Khi tính toán tấm theo phương pháp giải tích thì chỉ phù hợp với các dạng tấm có điều kiện biên đơn giản.

Thực tế, người ta thường sử dụng phương pháp số để giải quyết các bài toán phức tạp. Một trong các phương pháp số để phân tích tấm đó là phương pháp phần tử hữu hạn (FEM). Tuy nhiên phương pháp phần tử hữu hạn (FEM) có một số hạn chế là như kỹ thuật rời rạc, tốc độ tính toán, chi phí tính toán, tốc độ hội tụ, … Ngày nay với tiến bộ khoa học máy tính làm cho tốc độ tính toán dễ dàng cải thiện được, bên cạnh đó những đóng góp cải thiện tính toán cũng đóng góp vai trò đáng kể cho việc tính toán tấm này. Khi mô hình tấm, người ta thường phân ra các loại như tấm mỏng, tấm dày, tấm 3D.

Để tính toán tấm mỏng, ta dùng lý thuyết tấm cổ điển (lý thuyết Kirchhoff- Love) [1]. Hạn chế lý thuyết tấm cổ điển là bỏ qua các biến dạng trượt, nên sẽ không cho kết quả chính xác đối với tấm dày. Khi sử dụng lý thuyết tấm biến dạng cắt bậc nhất phân tích tấm dày [2], thường xảy ra hiện tượng “khóa cắt” (shear locking) làm kết quả thiếu chính xác với thực nghiệm. Hiện nay, có nhiều kỹ thuật khử hiện tượng “khóa cắt” như DSG, MIN, MITC, … Trong lý thuyết tấm biến dạng cắt bậc nhất, ta sử dụng hệ số hiệu chỉnh cắt để hiệu chỉnh ứng suất cắt giả định so với ứng suất cắt thực tế.

Tuy nhiên hệ số hiệu chỉnh cắt này thường khó xác định vì còn phụ thuộc nhiều yếu tố. 1 Để tránh hiện tượng “khóa cắt” và nâng cao độ chính xác của các công thức tính hiện tại, người ta đã thực hiện một số biện pháp làm trơn dựa trên cạnh, trên nút, trên ô như các các kỹ thuật ES-DSG, NS-DSG, CS-DSG…. Do đó, luận văn này trình bày một phương pháp để thiết lập phần tử tấm 3 nút có sử dụng kỹ thuật MITC3 để khử “khóa cắt” (shear locking) và làm trơn dựa trên cạnh biên (ES-MITC3). Ngoài ra còn có lý thuyết tấm biến dạng cắt bậc cao để tính toán tấm, ưu điểm lý thuyết cắt bậc cao là không cần hệ số hiệu chỉnh cắt.

Hiện nay đã có nhiều lý thuyết cắt bậc cao như lý thuyết biến dạng cắt bậc ba của Reddy [3], lý thuyết biến dạng cắt hình sin của Arya [4],… 1. Tình hình nghiên cứu ngoài nước. Việc phân tích tĩnh, ổn định, dao động tự do của kết cấu tấm đóng một vai trò quan trọng trong dự báo ứng xử của kết cấu. Việc nghiên cứu trên tấm có thể được tìm thấy trên các tài liệu tổng quan [5, 6], và có sự đóng góp trong các lĩnh vực chuyên ngành đặc biệt bởi Leissa [7-9] và Liew [10, 11].

Trong ứng dụng thực tế, phần tử tấm Reissner-Mindlin bậc thấp được dễ chấp nhận do tính đơn giản và hiệu quả của nó. Tuy nhiên, các phần tử tấm bậc thấp này thường xuất hiện hiện tượng “khóa cắt”. Để loại trừ “khóa cắt”, các loại tích phân có lựa chọn đã được đề nghị [12-15]. Ý tưởng của sơ đồ này là tách năng lượng biến dạng thành 2 phần, một phần do uốn và một phần do cắt.

Ngoài ra để khử hiện tượng khóa cắt, tác giả Phill-Seung Lee, Klaus-Jurgen Bathe [16] đã phát triển phần tử tấm tam giác đẳng hướng 3-node và 6-node (MITC3, MITC6-a, MITC6-b) bước đầu đã cho những kết quả tốt. Tuy nhiên, cần có những nghiên cứu xa hơn để có thể ứng dụng rộng rãi kỹ thuật này trong phân tích kết cấu tấm. Gần đây, Liu [17] đã đề nghị một phương pháp phần tử hữu hạn trơn dựa trên biên (ES-FEM) cho phân tích tĩnh, tự do và dao động của các tấm 2D. Kết quả số có độ chính xác cao đã chứng minh rằng ES-FEM có cho ra kết quả tốt: (1) mô 2 hình ES-FEM thường được tìm thấy độ hội tụ cao và thậm chí nhiều chính xác hơn FEM sử dụng phần tử 4 cạnh (FEM-Q4) với cùng số lượng nút; (2) không có kiểu năng lượng ảo nào và do vậy phương pháp này làm việc tốt và ổn định theo thời gian hơn cho phân tích động và (3) khi làm theo phương pháp này thì đơn giản hơn, và hiệu quả tính toán tốt hơn FEM khi sử dụng cùng số nút.

Trong đề tài này đề xuất một phương pháp kết hợp kỹ thuật làm trơn dựa trên cạnh sử dụng phần tử MITC3 để phân tích kết cấu tấm. Trên cơ sở đó, kế thừa các thành quả của các tác giả đi trước, với thực tiễn hiện nay, đề tài luận văn này là: “PHÂN TÍCH KẾT CẤU TẤM BẰNG PHẦN TỬ BIẾN DẠNG TRƠN ES- MITC3 “. Tình hình nghiên cứu trong nước. Song song với các nghiên cứu ở nước ngoài, một số nghiên cứu trong nước cũng đã đề cập đến vấn đề phân tích kết cấu tấm bằng phương pháp phần tử hữu hạn trơn.

Sau đây là một số nghiên cứu trong nước. Phùng Văn Phúc, và các cộng sự [18] nêu lên phân tích tĩnh và dao động của vật liệu composite và tấm nhiều lớp bởi phương pháp DSG3 dựa trên cạnh biên sử dụng phần tử 3 nút dựa theo lý thuyết nhiều lớp (ES-DSG3). Nguyễn Thời Trung và các cộng sự [19] nêu lên phương pháp phần tử tấm 3 nút Mindlin (MIN3) dựa trên biên cho phân tích tấm tĩnh và dao động tự do (ES- MIN3). Phan Dao và các cộng sự [20] đưa ra phân tích tĩnh và dao động tự do của kết cấu tấm composite nhiều lớp có chứa lớp áp điện bằng phương pháp phần tử hữu hạn trơn dựa trên cạnh (ES-FEM).

Các nghiên cứu trong nước về phương pháp phần tử hữu hạn trơn về phân tích kết cấu tấm chưa nhiều và chưa có nghiên cứu nào đề cập đến phương pháp phần tử hữu hạn trơn dựa trên nút sử dụng phần tử MITC3 (ES-MITC3). Mục đích nghiên cứu. Nghiên cứu xây dựng phần tử tấm phẳng 3 nút sử dụng kỹ thuật MITC3 trên trường biến dạng trơn trên cạnh (ES-MITC3) với hy vọng cải thiện độ chính xác hơn so với phần tử MITC3 thông thường. Nhiệm vụ của đề tài.

Xây dựng lý thuyết. Lập trình trên phần mềm Matlab. Tính toán các ví dụ số để kiểm tra lý thuyết và chương trình. Đánh giá và phân tích các công thức vừa tìm được.

Cơ sở lý luận và phương pháp nghiên cứu. Để thực hiện đề tài này, tác giả sử dụng phương pháp sau: là phương pháp nghiên cứu lý thuyết, sử dụng lý thuyết biến dạng cắt bậc nhất kết hợp với kỹ thuật MITC3 và kỹ thuật làm trơn trên cạnh. Lý thuyết tính toán sau khi thiết lập sẽ được lập trình bằng phần mềm Matlab để tính toán một số bài toán ví dụ số minh họa của kết cấu tấm. Kết quả mô phỏng số của các bài toán điển hình về kết cấu tấm sẽ được so sánh với kết quả của các nghiên cứu trước đó để đánh giá tính hiệu quả và chính xác của nghiên cứu này.

4 Chương 2 LÝ THUYẾT TẤM BIẾN DẠNG CẮT BẬC NHẤT 2. Giả thuyết: Lý thuyết tấm biến dạng cắt (shear deformation theories of plates) [2], hay còn được gọi là lý thuyết tấm dày (thick plate) chịu uốn, hoặc còn gọi là lý thuyết tấm Reissner-Mindlin do Levy-Reissner-Hencky-Mindlin thực hiện, dựa trên giả thiết của 2 tác giả Reissner-Mindlin, giống giả thiết lý thuyết dầm Timoshenko chịu cắt [21]. Trong đó, điều kiện: tấm dày là tỉ lệ h/a ≥ 1/20 (với: h = chiều dày tấm; a = kích thước nhỏ nhất chiều dài mặt trung bình tấm), tấm mỏng là tỉ lệ h/a < 1/20. Giả thiết: Đoạn thẳng vuông góc với mặt trung bình của tấm sau khi biến dạng sẽ chuyển vị võng xuống một đoạn w(x, y) và góc xoay βx, βy.

Trường chuyển vị: Trong lý thuyết tấm Mindlin-Reissner, khi tấm phẳng chịu tải trọng phân bố q thì đoạn thẳng vuông góc mặt trung bình tấm sau khi biến dạng sẽ chuyển vị võng xuống 1 đoạn thẳng w theo phương z (trong đó w là 1 hàm phụ thuộc x, y) và 2 chuyển vị góc xoay βx, βy (là 2 góc xoay của đoạn thẳng vuông góc với mặt trung bình của tấm lần lượt quay quanh trục y, x). Cả 3 thành phần này (w, βx, βy) được biểu diễn như Hình 2. 5 x,u  y z,w F E B H F' h/2 A G D w  h/2 Z E x C w E' y,v B' w A v A' D H' H G' u D' D C'  yz wy    w xz y x  x Hình 2.1: Trường chuyển vị trong tấm.1)   w(x, y, z) = w(x, y) Với: u, v là chuyển vị thẳng theo phương x, y. Trường biến dạng: Từ trường chuyển vị (2.1), ta có trường biến dạng:  u x x   z  x x  v y   y  z  y y z  0   u v   x y   xy   z      y x  x y   xz  u  w   x  w z x x    v w w   yz    y   z y y Trong đó: εx, εy, εz là các biến dạng dài theo phương x, y, z γxy, γyz, γzx là các biến dạng trượt trong mặt phẳng Oxy, Oxz, Oyz 2.

Trường ứng suất: x  xz  xy  x    y yz yx y Hình 2.2: Các thành phần ứng suất của tấm 7 Dựa vào định luật Hooke, theo lý thuyết ứng suất phẳng, từ trường biến dạng (2.2) cho ra trường ứng suất:  E   x  1 2  x  y    E   y  12  x   y     Z 0   E    xy  2(1) xy   E    xz  2(1) xz   E    yz  yz  2(1) Trong đó: σx, σy, σz là các ứng suất pháp theo phương x, y, z. τxy, τxz, τyz là các ứng suất tiếp trên mặt có vectơ pháp tuyến là trục x, y, z. E là module đàn hồi khi kéo nén của vật liệu. ν là hệ số poisson.

Khi đó, trường ứng suất (2.