Tổng quan nghiên cứu

Trong kỹ thuật kết cấu hiện đại, kết cấu tấm và vỏ làm từ vật liệu composite nhiều lớp gia cường sườn đóng vai trò cốt lõi trong các ngành hàng không, hàng hải, ô tô và xây dựng công trình vượt nhịp lớn nhờ tỷ số độ cứng trên trọng lượng vượt trội. Tuy nhiên, việc mô phỏng ứng xử tĩnh của kết cấu này bằng phương pháp phần tử hữu hạn truyền thống (FEM) thường gặp hiện tượng khóa cắt (shear locking) nghiêm trọng khi chiều dày tấm giảm xuống dưới 1% kích thước cạnh, đồng thời ma trận độ cứng tính toán thường quá cứng so với ứng xử thực tế.

Nghiên cứu tập trung giải quyết bài toán phân tích tĩnh kết cấu tấm và vỏ composite nhiều lớp có sườn thông qua việc phát triển phần tử vỏ phẳng tam giác 3 nút kết hợp kỹ thuật phần tử hữu hạn trơn trên miền phần tử (CS-FEM) và kỹ thuật khử khóa cắt MITC3+ (CS-MITC3+). Mục tiêu cụ thể là thiết lập công thức toán học và thuật toán số tích hợp 6 bậc tự do trên mỗi nút đỉnh cùng 1 nút nổi ở trọng tâm, sau đó kết nối đồng bộ với phần tử dầm 2 nút gia cường.

Nghiên cứu được triển khai thực nghiệm số trên môi trường MATLAB tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh. Các khảo sát tính toán được thực hiện trên các mô hình tấm hình vuông kích thước 25,4 mm x 25,4 mm, chiều dày 0,254 mm, mô đun đàn hồi 117,215 x 10^3 N/mm2, chịu tải trọng phân bố đều 0,006895 N/mm2 và tải trọng tập trung. Kết quả đóng góp giải pháp tính toán giúp nâng cao độ chính xác chuyển vị, kiểm soát sai số hội tụ dưới 2% đối với tấm composite có sườn so với các lời giải giải tích chuẩn.

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 xây dựng trên nền tảng lý thuyết biến dạng cắt bậc nhất (FSDT - First-order Shear Deformation Theory) của Reissner và Mindlin cho tấm/vỏ và dầm. Khác với lý thuyết tấm cổ điển Kirchhoff-Love vốn bỏ qua biến dạng cắt ngang (gamma_xz = gamma_yz = 0), FSDT giả thiết đoạn thẳng vuông góc với mặt trung bình trước khi biến dạng vẫn thẳng nhưng không còn vuông góc sau khi biến dạng, cho phép tính toán hiệu quả cho cả kết cấu dày và mỏng với hệ số hiệu chỉnh cắt tiêu chuẩn 5/6.

Mô hình cấu tạo vật liệu sử dụng lý thuyết lớp tương đương (Equivalent Single Layer) để thiết lập quan hệ ứng suất - biến dạng của composite nhiều lớp trong hệ tọa độ cục bộ và toàn cục thông qua ma trận độ cứng giảm chuyển đổi. Nghiên cứu tích hợp 3 khái niệm cơ học cốt lõi: trường chuyển vị kết hợp bậc tự do xoay vuông góc mặt phẳng trung bình (drilling rotation) theo công thức Allman với 6 bậc tự do tại mỗi nút, kỹ thuật nội suy ten-xơ hỗn hợp MITC3+ nhằm loại bỏ năng lượng cắt ký sinh, và kỹ thuật làm trơn biến dạng trên miền phần tử con (Cell-based Smoothing) giúp ma trận độ cứng trở nên linh hoạt hơn.

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

Nguồn dữ liệu của đề tài được tạo lập và kiểm chứng từ các bài toán chuẩn (benchmarks) trong cơ học tính toán kết cấu tấm, vỏ phẳng, vỏ trụ và vỏ cầu. Phương pháp phân tích trung tâm là phương pháp phần tử hữu hạn trơn CS-FEM phát triển trên phần tử tam giác 3 nút MITC3+. Mỗi phần tử tam giác được chia thành 3 miền tam giác con liên kết từ 3 đỉnh đến nút nổi tại trọng tâm. Biến dạng màng và biến dạng uốn được làm trơn trên từng miền tam giác con, chuyển tích phân diện tích sang tích phân đường trên các cạnh biên với 2 điểm Gauss. Biến dạng cắt ngoài mặt phẳng được khử khóa cắt thông qua 6 điểm buộc (A, B, C, D, E, F) với thông số định vị d = 0,0001.

Mô hình dầm gia cường được rời rạc bằng phần tử dầm 2 nút với 3 bậc tự do mỗi nút (chuyển vị dọc trục, chuyển vị đứng và góc xoay uốn), sau đó chuyển đổi và ghép nối vào hệ tọa độ tổng thể 6 bậc tự do để làm việc đồng thời với tấm/vỏ. Cỡ mẫu khảo sát gồm 8 bài toán kết cấu phức hợp với mật độ lưới phần tử chia mịn dần theo các cấp độ 8x8x2, 12x12x2, 16x16x2, 20x20x2 và 24x24x2 phần tử tam giác. Toàn bộ thuật toán được lập trình tự động hóa bằng ngôn ngữ MATLAB trong thời gian thực hiện đề tài.

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ố trên các mô hình tấm và vỏ khẳng định tính ổn định và khả năng khử khóa cắt vượt trội của phần tử CS-MITC3+:

Thứ nhất, đối với tấm hình vuông đồng nhất có 1 sườn giữa chịu tải trọng phân bố đều 0,006895 N/mm2, độ võng tại tâm tấm hội tụ ổn định khi số phần tử tăng từ 8x8x2 lên 24x24x2. Giá trị độ võng không thứ nguyên đạt mức xấp xỉ 0,0132 mm, sai lệch so với nghiệm tham khảo của các tác giả quốc tế dưới 1,5%.

Thứ hai, ở bài toán tấm chữ nhật và tấm vuông composite nhiều lớp gia cường 1 hoặc 2 sườn chịu tải phân bố đều và tải tập trung, phần tử đề xuất phản ánh chính xác ảnh hưởng của góc sắp xếp lớp sợi và độ cứng sườn. Với điều kiện biên 4 cạnh tựa đơn hoặc liên kết ngàm hỗn hợp, độ võng cực đại giảm từ 35% đến 58% khi bố trí hệ sườn trực giao so với tấm không có sườn.

Thứ ba, đối với kết cấu vỏ trụ đồng nhất ngàm một cạnh chịu tải tập trung và vỏ cầu đồng nhất có 2 sườn giao nhau, chuyển vị hướng tâm và chuyển vị đứng cho kết quả tương đồng cao với lời giải phần tử hữu hạn thương mại, tốc độ hội tụ nhanh hơn phần tử tam giác chuẩn khoảng 12%.

Thứ tư, khi phân tích vỏ composite nhiều lớp ngàm một cạnh có 3 sườn chịu tải tập trung, kết quả chuyển vị hướng tâm và chuyển vị tiếp tuyến của phần tử CS-MITC3+ xuất hiện độ lệch đáng kể (khoảng 15% đến 22%) so với các công bố tham khảo trước đó, chỉ ra tính phi tuyến hình học và hiệu ứng tương tác uốn - xoắn phức tạp ở vỏ cong composite.

Thảo luận kết quả

Khả năng hội tụ chính xác của CS-MITC3+ bắt nguồn từ việc kết hợp phép làm trơn biến dạng miền (CS-FEM) giúp làm mềm ma trận độ cứng của hệ phần tử tam giác 3 nút – vốn thường quá cứng trong FEM truyền thống. Việc chuyển đổi tích phân đạo hàm hàm dạng sang tích phân đường trên biên các tam giác con giúp giảm thiểu sự nhạy cảm với độ méo của lưới phần tử.

Dữ liệu hội tụ được thể hiện trực quan qua đồ thị quan hệ giữa độ võng tại tâm và mật độ chia lưới (từ 128 phần tử lên 1.152 phần tử tam giác), minh chứng đường cong chuyển vị tiệm cận trơn tru về nghiệm chính xác. Bên cạnh đó, các bảng so sánh chuyển vị không thứ nguyên cho thấy phần tử mới loại bỏ triệt để hiện tượng khóa cắt ngay cả khi tỷ lệ bề dày trên nhịp giảm xuống mức 1/1000.

Sự khác biệt số liệu ở bài toán vỏ cong composite nhiều lớp có sườn được giải thích do lý thuyết FSDT tuyến tính chưa bao quát hết biến dạng trượt bậc cao giữa các lớp khi vỏ chịu tải tập trung tại điểm đặt lực lệch tâm, mở ra định hướng hiệu chỉnh ma trận ghép nối vỏ - dầm cong.

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

Dựa trên các phát hiện nghiên cứu, các giải pháp và khuyến nghị kỹ thuật được đề xuất như sau:

Thứ nhất, tối ưu hóa góc định hướng sợi composite và kích thước mặt cắt sườn gia cường. Doanh nghiệp thiết kế kết cấu nên áp dụng cấu hình góc dệt sợi 0/90 độ kết hợp +/-45 độ với chiều cao sườn tối thiểu gấp 10 lần chiều dày vỏ để nâng cao khả năng chịu uốn lên 40%, áp dụng ngay trong quy trình tính toán thiết kế vỏ tàu thủy và thân vỏ xe cơ giới.

Thứ hai, mở rộng thuật toán phần tử CS-MITC3+ kết hợp lý thuyết biến dạng cắt bậc cao (HSDT) và lý thuyết từng lớp (Layerwise Theory). Nhóm nghiên cứu tại các trường đại học cần triển khai nâng cấp mã nguồn trong khung thời gian 6 đến 12 tháng nhằm giải quyết triệt để sai số chuyển vị trên vỏ cong composite, hướng tới mục tiêu giảm sai lệch tính toán xuống dưới 3%.

Thứ ba, phát triển module tính toán phân tích phi tuyến hình học, dao động tự do và ổn định động lực học. Các kỹ sư CAE cần mở rộng chương trình tính toán hiện tại để đánh giá tải trọng tới hạn gây mất ổn định mất nếp cục bộ của sườn trong vòng 9 tháng tới.

Thứ tư, thương mại hóa và chuyển đổi mã nguồn MATLAB sang các thư viện tính toán hiệu năng cao như C++ hoặc Python kết nối API vào các phần mềm thương mại. Chủ thể thực hiện là các viện nghiên cứu cơ học kết cấu phối hợp cùng doanh nghiệp phần mềm trong lộ trình 18 tháng.

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

Nội dung luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cho 4 nhóm đối tượng chính:

Nhóm kỹ sư thiết kế kết cấu công trình và giao thông: Tiếp cận quy trình tính toán chính xác ứng xử chịu lực của sàn bản phẳng, vách mỏng composite có sườn tăng cứng, giúp tiết kiệm từ 15% đến 20% khối lượng vật liệu trong các dự án thực tế.

Nhóm lập trình viên và chuyên gia phát triển phần mềm mô phỏng CAE: Sử dụng trọn bộ công thức toán học ma trận độ cứng phần tử CS-MITC3+, kỹ thuật nén bậc tự do tại nút nổi và ma trận chuyển đổi tọa độ toàn cục để xây dựng các phần tử mới cho phần mềm phần tử hữu hạn.

Các nhà nghiên cứu và học viên cao học ngành Cơ kỹ thuật, Kỹ thuật Xây dựng: Sử dụng luận văn như tài liệu tham khảo chuyên sâu về phương pháp phần tử hữu hạn trơn (S-FEM), kỹ thuật khử khóa cắt MITC và cơ học composite nhiều lớp.

Nhóm chuyên gia thẩm tra và kiểm chuẩn kết cấu hàng không, hàng hải: Ứng dụng các bảng số liệu chuẩn và mô hình phân tích để đối chuẩn độ an toàn của vỏ tàu thủy composite hoặc cánh vỏ máy bay chịu tải trọng tĩnh phân bố và tập trung.

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

Phần tử CS-MITC3+ có ưu điểm gì vượt trội so với phần tử tam giác 3 nút tiêu chuẩn?
Phần tử CS-MITC3+ giải quyết triệt để hiện tượng khóa cắt nhờ kỹ thuật nội suy biến dạng cắt MITC3+ qua 6 điểm buộc. Đồng thời, kỹ thuật làm trơn trên 3 miền tam giác con giúp ma trận độ cứng linh hoạt hơn, cho kết quả hội tụ chuẩn xác với số lượng phần tử ít hơn khoảng 20% so với phần tử tam giác thông thường.

Tại sao cần bổ sung bậc tự do xoay quanh trục pháp tuyến (drilling rotation) vào phần tử?
Trong bài toán tấm vỏ không gian, việc bổ sung bậc tự do thứ 6 tại mỗi nút theo công thức Allman giúp việc ghép nối ma trận độ cứng giữa các phần tử nghiêng bất kỳ và giữa phần tử tấm với phần tử dầm trong hệ tọa độ 3 chiều diễn ra thuận lợi, tránh xuất hiện ma trận suy biến.

Kỹ thuật nén bậc tự do (static condensation) đóng vai trò gì trong mô hình?
Phần tử sử dụng 1 nút nổi ở trọng tâm có 2 bậc tự do xoay nhằm làm giàu hàm dạng uốn. Kỹ thuật nén bậc tự do giúp khử các biến số của nút nổi ở cấp độ phần tử trước khi lắp ghép ma trận tổng thể, giữ nguyên kích thước ma trận 18x18, tối ưu hóa bộ nhớ và tốc độ giải hệ phương trình.

Nguyên nhân dẫn đến sai số chuyển vị ở bài toán vỏ cong composite nhiều lớp có sườn là gì?
Sai số xuất phát từ giả thiết biến dạng cắt bậc nhất FSDT chưa mô tả trọn vẹn sự phân bố phi tuyến của ứng suất cắt qua chiều dày của vỏ cong nhiều lớp không đối xứng, kết hợp với hiệu ứng lệch tâm hình học phức tạp giữa sườn dầm và mặt trung bình vỏ.

Luận văn sử dụng công cụ nào để mô phỏng và kiểm chứng kết quả?
Tác giả tự xây dựng toàn bộ mã nguồn tính toán bằng ngôn ngữ MATLAB, thiết lập thuật toán từ bước tạo lưới, tính tích phân Gauss trên cạnh tam giác con, lắp ghép ma trận tổng thể đến giải hệ phương trình và xuất biểu đồ chuyển vị so sánh với các nghiên cứu uy tín.

Kết luận

Luận văn đã hoàn thành các mục tiêu nghiên cứu trọng tâm thông qua các kết quả cụ thể:

  • Thiết lập thành công công thức phần tử vỏ tam giác phẳng trơn 3 nút CS-MITC3+ có 18 bậc tự do kết hợp bậc tự do xoay Allman và kỹ thuật khử khóa cắt.
  • Xây dựng mô hình tương tác hoàn chỉnh giữa phần tử tấm/vỏ CS-MITC3+ và phần tử dầm composite 2 nút gia cường theo lý thuyết biến dạng cắt bậc nhất FSDT.
  • Lập trình hoàn chỉnh hệ thống tính toán số trên MATLAB, giải quyết chính xác 8 bài toán phân tích tĩnh kết cấu tấm và vỏ có sườn.
  • Đạt độ chính xác cao và tốc độ hội tụ vượt trội cho các kết cấu tấm phẳng đồng nhất và tấm composite nhiều lớp có sườn với sai số dưới 2%.
  • Chỉ ra giới hạn tính toán của mô hình FSDT tuyến tính trên vỏ cong composite có sườn, tạo tiền đề vững chắc cho việc phát triển lý thuyết bậc cao HSDT trong giai đoạn 2026-2027.

Các đơn vị nghiên cứu và kỹ sư kết cấu quan tâm có thể ứng dụng trực tiếp thuật toán CS-MITC3+ vào các bài toán tính toán tối ưu hóa kết cấu tấm composite gia cường nhằm nâng cao hiệu quả kỹ thuật và giảm giá thành công trình.