Tổng quan nghiên cứu

Vật liệu composite nhiều lớp ngày càng đóng vai trò then chốt trong các ngành công nghiệp kỹ thuật cao như hàng không vũ trụ, hàng hải và xây dựng công trình hiện đại nhờ tỷ số độ cứng trên trọng lượng vượt trội hơn từ 30% đến 50% so với thép kết cấu truyền thống. Tuy nhiên, việc mô phỏng và phân tích chính xác ứng xử cơ học của các kết cấu tấm composite gặp nhiều thách thức lớn, đặc biệt là hiện tượng khóa cắt (shear locking) xảy ra khi tỷ số bề dày trên chiều dài tấm (t/L) giảm xuống mức từ 0,01 đến 0,0001. Bên cạnh đó, các phương pháp phần tử hữu hạn truyền thống thường đòi hỏi lưới phần tử tứ giác chất lượng cao, gây khó khăn khi xử lý các miền hình học phức tạp và làm giảm tốc độ hội tụ của lời giải số.

Đề tài luận văn thạc sĩ chuyên ngành Kỹ thuật Xây dựng và Cơ học Ứng dụng được thực hiện bởi tác giả Nguyễn Hòa dưới sự hướng dẫn của Tiến sĩ Châu Đình Thành tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh, hoàn thành vào tháng 10 năm 2015. Mục tiêu trọng tâm của nghiên cứu là xây dựng và phát triển phần tử tam giác 3 nút kết hợp kỹ thuật trơn hóa biến dạng dựa trên cạnh và nội suy hỗn hợp thành phần ten xơ, tạo nên mô hình tính toán số ES-MITC3. Luận văn ứng dụng mô hình này để phân tích toàn diện ứng xử tĩnh và dao động tự do của tấm composite nhiều lớp dựa trên lý thuyết biến dạng cắt bậc cao của Reddy với 7 bậc tự do trên mỗi nút. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu sai số tính toán dao động xuống dưới 2%, tối ưu hóa thời gian xử lý ma trận độ cứng và cung cấp công cụ tính toán số tin cậy cho kỹ sư thiết kế kết cấu composite mỏng và dày vừa.

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 được thiết lập dựa trên hệ thống cơ sở lý thuyết cơ học kết cấu vững chắc, kết hợp hai nền tảng lý thuyết chủ đạo:

Thứ nhất là Lý thuyết biến dạng cắt bậc cao (TSDT hay HSDT) do Reddy phát triển năm 1997. Khác với lý thuyết tấm cổ điển (CLPT) dựa trên giả thuyết Kirchhoff chỉ phù hợp cho tấm mỏng hay lý thuyết biến dạng cắt bậc nhất (FSDT) đòi hỏi hệ số hiệu chỉnh cắt k = 5/6, lý thuyết của Reddy giả thiết trường chuyển vị theo hàm bậc 3 của tọa độ bề dày z. Mô hình này thiết lập 7 bậc tự do độc lập cho mỗi nút gồm 3 chuyển vị và 4 thành phần xoay, độ cong, đảm bảo điều kiện biên ứng suất cắt bằng 0 tại mặt trên z = h/2 và mặt dưới z = -h/2 của tấm mà không cần bất kỳ hệ số hiệu chỉnh nhân tạo nào.

Thứ hai là Phương pháp phần tử hữu hạn trơn dựa trên cạnh (ES-FEM) do Liu và các cộng sự đề xuất năm 2008, kết hợp kỹ thuật nội suy hỗn hợp ten xơ biến dạng (MITC3) của Bathe. Mô hình tính toán tích hợp các khái niệm cốt lõi bao gồm miền trơn cục bộ liên kết với từng cạnh phần tử, điểm buộc hình học (tying points) và ma trận độ cứng trơn hóa. Kỹ thuật này giúp chuyển đổi phép tính tích phân từ diện tích phần tử sang các miền trơn biên, loại bỏ hoàn toàn hiện tượng suy biến dạng mode và triệt tiêu độ cứng ảo do khóa cắt khi tỷ lệ chiều dày tấm t/L đạt ngưỡng 10 mũ trừ 4.

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

Nghiên cứu sử dụng phương pháp mô phỏng số thông qua việc lập trình toàn bộ thuật toán phần tử ES-MITC3 trên môi trường phần mềm MATLAB, được thực hiện trong giai đoạn từ đầu năm 2015 đến tháng 10 năm 2015.

Về cỡ mẫu và cấu hình khảo sát, tác giả tiến hành thử nghiệm trên hệ thống 5 cấp độ phân chia lưới phần tử gồm 2x2, 4x4, 8x8, 16x16 và 32x32 phần tử trên cạnh tấm, tương ứng với số lượng phần tử tam giác lên tới 2048 phần tử và khảo sát 6 dạng mode dao động riêng đầu tiên. Phương pháp chọn mẫu tập trung vào các trường hợp kết cấu chuẩn mực kinh điển trong cơ học kết cấu để đối sánh, bao gồm: tấm vuông đồng nhất đẳng hướng với mô đun Young E = 1092000, tấm composite 4 lớp trực hướng (0/90/90/0) chịu tải phân bố đều và tải hình sin, tấm composite 16 lớp đối xứng với góc nghiêng sợi 45 độ, và tấm composite 3 lớp (0/90/0) chịu liên kết ngàm 4 cạnh.

Lý do lựa chọn phương pháp phân tích ES-MITC3 kết hợp lý thuyết Reddy xuất phát từ tính ưu việt của phần tử tam giác 3 nút trong việc tự động hóa chia lưới cho các kết cấu phức tạp. Việc trơn hóa biến dạng trên cạnh giúp làm mềm ma trận độ cứng của mô hình phần tử hữu hạn liên tục, tăng tốc độ hội tụ năng lượng và đảm bảo kết quả ứng suất tiếp xúc giữa các lớp sợi đạt độ chính xác tương đương lời giải đàn hồi 3D nhưng tiết kiệm tới hơn 40% chi phí tính toán.

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

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

Quá trình phân tích số trên các mô hình tấm composite đã mang lại 4 phát hiện quan trọng:

Thứ nhất, phần tử ES-MITC3 thể hiện khả năng khử hiện tượng khóa cắt vượt trội khi tỷ số độ dày tấm t/L giảm liên tục từ 10 mũ trừ 1 xuống 10 mũ trừ 4. Ở tỷ số t/L = 0,0001, độ võng không thứ nguyên tại tâm tấm của mô hình đạt giá trị xấp xỉ chính xác lời giải giải tích với độ lệch dưới 0,05%, trong khi các phần tử tam giác cổ điển thường bị cứng hóa và sai lệch trên 85%.

Thứ hai, tốc độ hội tụ của ES-MITC3 nhanh hơn đáng kể so với các phần tử tam giác 3 nút cùng thế hệ như MIN3, DSG3 và ESDSG3. Khi tăng mật độ lưới từ 2x2 lên 32x32 trên cạnh biên, độ võng tấm hội tụ ổn định và giảm sai số dư thừa hơn 15% so với mô hình FSDT có hệ số hiệu chỉnh k = 5/6.

Thứ ba, đối với tấm composite 4 lớp trực hướng (0/90/90/0) có tỷ số mô đun E1/E2 = 25 chịu tải trọng hình sin và phân bố đều, giá trị ứng suất tiếp và ứng suất cắt trong mặt phẳng hoàn toàn tương thích với lời giải giải tích đàn hồi 3 chiều 3D-FEM của tác giả Noor và phương pháp RPF-PS-Layerwise, với sai số ứng suất pháp chỉ dao động từ 1,2% đến 2,8%.

Thứ tư, khi phân tích dao động tự do trên tấm composite 3 lớp (0/90/0) có tỷ lệ mô đun đàn hồi E1/E2 = 40, tần số dao động không thứ nguyên ở mode 1 tại tỷ số L/t = 5 đạt 4,738, tiệm cận chặt chẽ với kết quả của tác giả Liew bằng phương pháp p-Ritz năm 1996 và phương pháp không lưới MLSDQ năm 2003 với mức chênh lệch dưới 1,5%.

Thảo luận kết quả

Sự chính xác và ổn định của phương pháp ES-MITC3 bắt nguồn từ việc kết hợp kỹ thuật trơn hóa cạnh ES-FEM giúp khắc phục hiện tượng quá cứng của phần tử hữu hạn tiêu chuẩn, kết hợp với các điểm buộc tying points theo cơ chế MITC để khử triệt để năng lượng cắt ký sinh. So sánh với các nghiên cứu của Liew (1996), Zhen và Wanji (2006) hay Ferreira và Fasshauer, mô hình ES-MITC3 cho thấy sự vượt trội về độ ổn định vì không làm xuất hiện các dạng mode dao động giả (spurious modes) trong toàn bộ 6 mode dao động khảo sát.

Trong các báo cáo kỹ thuật và công bố khoa học, toàn bộ kết quả phân tích nên được trực quan hóa thông qua hệ thống bảng so sánh đa chiều và biểu đồ đường logarit. Cụ thể, dữ liệu độ võng không thứ nguyên theo sự thay đổi log(L/t) từ 1 đến 4 được trình bày trực quan trên biểu đồ phân tán dạng đường để làm nổi bật ưu thế không bị khóa cắt của ES-MITC3 so với MIN3 và DSG3. Đồng thời, biểu đồ phân bố ứng suất cắt ngang sigma_xz và sigma_yz theo chiều cao tọa độ z từ âm h/2 đến dương h/2 được thể hiện dưới dạng đường cong bậc hai liên tục qua các lớp vật liệu. Cách trình bày này giúp các kỹ sư kết cấu dễ dàng nhận biết sự biến thiên ứng suất tiếp xúc giữa các lớp sợi 0 độ và 90 độ, từ đó phát hiện sớm nguy cơ phá hủy tách lớp (delamination) trong các kết cấu composite chịu uốn.

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

Dựa trên những phát hiện và hạn chế được chỉ ra trong nghiên cứu, 4 khuyến nghị và giải pháp chiến lược được đề xuất như sau:

Thứ nhất, tích hợp lý thuyết biến dạng cắt lớp thông minh (Layer-wise theory) vào thuật toán phần tử ES-MITC3 hiện hữu. Nhóm nghiên cứu cơ học tính toán tại các trường đại học kỹ thuật cần chủ trì thực hiện giải pháp này trong lộ trình 12 tháng, nhằm nâng cao độ chính xác dự báo ứng suất cắt ngang qua bề mặt tiếp xúc giữa các lớp vật liệu, hướng tới giảm thiểu sai số cục bộ từ mức 8% hiện tại xuống dưới 2%.

Thứ hai, phát triển mô-đun phần mềm tính toán tự động hóa hoàn toàn trên nền tảng MATLAB và Python. Các kỹ sư phần mềm kết cấu và chuyên gia phân tích CAE cần triển khai trong vòng 6 tháng để đóng gói thuật toán ES-MITC3 thành một hộp công cụ (toolbox) thương mại hóa hoặc mã nguồn mở, giúp giảm thời gian thiết lập bài toán từ 4 giờ xuống còn dưới 15 phút cho các kỹ sư công trình.

Thứ ba, mở rộng phạm vi ứng dụng phần tử ES-MITC3 cho bài toán phân tích phi tuyến hình học và phân tích phá hủy vật liệu composite chịu tải trọng động lực học cao. Các viện nghiên cứu hàng không và đóng tàu cần tiến hành thử nghiệm trong thời gian 18 tháng, hướng tới mục tiêu tối ưu hóa trọng lượng kết cấu thân vỏ tàu biển và cánh máy bay giảm từ 15% đến 25% mà vẫn đảm bảo hệ số an toàn chịu lực trên 1,5.

Thứ tư, chuẩn hóa quy trình kiểm định và đối chiếu số liệu mô phỏng số với dữ liệu thực nghiệm đo đạc biến dạng thực tế. Ban quản lý chất lượng và các phòng thí nghiệm trọng điểm công trình cần phối hợp thực hiện trong thời hạn 9 tháng với tối thiểu 50 mẫu thử nghiệm uốn tĩnh và rung động thực tế, nhằm xác lập ngưỡng tin cậy đạt trên 98% trước khi áp dụng vào sản xuất đại trà.

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

Luận văn thạc sĩ là nguồn tài liệu học thuật và kỹ thuật giá trị cao cho 4 nhóm đối tượng trọng tâm sau:

Nhóm thứ nhất là học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Cơ học tính toán, Kỹ thuật Xây dựng và Cơ kỹ thuật. Tài liệu cung cấp trọn vẹn cơ sở toán học 7 bậc tự do và kỹ thuật lập trình phần tử ES-MITC3, phục vụ việc phát triển các đề tài nghiên cứu chuyên sâu và công bố bài báo khoa học quốc tế.

Nhóm thứ hai là kỹ sư thiết kế kết cấu trong ngành công nghiệp hàng không vũ trụ và đóng tàu cao tốc. Tài liệu mang lại giải pháp mô phỏng chính xác độ võng và 6 mode dao động của tấm composite mỏng với tỷ số t/L đạt 10 mũ trừ 4, giúp tối ưu hóa cấu trúc tấm composite 16 lớp giảm trọng lượng khung vỏ tới 20%.

Nhóm thứ ba là chuyên gia phát triển phần mềm mô phỏng CAE (Computer-Aided Engineering). Nhóm này có thể ứng dụng trực tiếp các ma trận độ cứng trơn hóa và thuật toán nội suy ten xơ để tích hợp tính năng phân tích tấm mỏng vào các phần mềm thương mại, giúp tăng tốc độ giải ma trận lên 30%.

Nhóm thứ tư là các đơn vị tư vấn thẩm tra và quản lý chất lượng dự án công trình dân dụng sử dụng vật liệu mới. Tài liệu cung cấp cơ sở dữ liệu số chuẩn xác với hơn 5 bảng tra cứu tần số dao động và biểu đồ phân bố ứng suất, làm căn cứ thẩm định an toàn chịu lực cho các kết cấu sàn composite nhịp lớn.

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

Câu hỏi 1: Hiện tượng khóa cắt (shear locking) là gì và phần tử ES-MITC3 giải quyết vấn đề này như thế nào? Trả lời: Hiện tượng khóa cắt xảy ra khi kết cấu tấm rất mỏng, khiến năng lượng biến dạng cắt bị phóng đại giả tạo làm tấm cứng bất thường. Phần tử ES-MITC3 giải quyết triệt để vấn đề này bằng cách kết hợp kỹ thuật nội suy biến dạng ten xơ tại các điểm buộc hình học và làm trơn biến dạng trên cạnh, giúp duy trì độ chính xác cao ngay cả khi tỷ số t/L giảm xuống mức 0,0001.

Câu hỏi 2: Tại sao luận văn lại sử dụng phần tử tam giác 3 nút thay vì phần tử tứ giác 4 nút hoặc 8 nút? Trả lời: Phần tử tam giác 3 nút có khả năng thích ứng linh hoạt vượt trội trong việc chia lưới tự động cho các miền hình học phức tạp và biên cong. Mặc dù phần tử tam giác tiêu chuẩn thường có độ cứng lớn, sự kết hợp với phương pháp trơn hóa ES-FEM giúp tăng tốc độ hội tụ hơn 15% mà không cần bổ sung nút nội suy bậc cao tốn kém.

Câu hỏi 3: Ưu điểm nổi bật của lý thuyết biến dạng cắt bậc cao Reddy so với lý thuyết Mindlin là gì? Trả lời: Lý thuyết Reddy sử dụng hàm chuyển vị bậc 3 theo bề dày, cho phép phân bố ứng suất cắt dạng parabol và tự động triệt tiêu ứng suất cắt tại hai mặt biên tự do trên và dưới tấm. Điều này loại bỏ hoàn toàn sự phụ thuộc vào hệ số hiệu chỉnh cắt k = 5/6 vốn chỉ mang tính giả định của lý thuyết Mindlin.

Câu hỏi 4: Phương pháp ES-MITC3 có thể áp dụng cho các cấu hình tấm composite nhiều lớp nào? Trả lời: Phương pháp đã được kiểm chứng thành công trên nhiều cấu hình phức tạp, bao gồm tấm trực hướng 4 lớp (0/90/90/0), tấm 3 lớp đối xứng và tấm đa hướng 16 lớp ((45/90/-45/0)2)sym. Kết quả tính toán ứng suất và độ võng đều đạt độ tin cậy cao khi so sánh với lời giải giải tích 3 chiều của Noor.

Câu hỏi 5: Phần mềm nào được sử dụng để tính toán và làm thế nào để tái lập kết quả nghiên cứu? Trả lời: Tác giả xây dựng toàn bộ thuật toán phân tích tĩnh và dao động tự do trên nền tảng MATLAB. Người đọc có thể tái lập kết quả bằng cách thiết lập ma trận độ cứng ES-MITC3 với 7 bậc tự do cho mỗi nút và chạy thử nghiệm trên các cấu hình lưới từ 2x2 đến 32x32 phần tử đã được mô tả chi tiết.

Kết luận

Luận văn thạc sĩ của tác giả Nguyễn Hòa đã giải quyết thành công bài toán cơ học kết cấu phức tạp đối với tấm composite nhiều lớp, thể hiện qua 5 nội dung trọng tâm:

  • Phát triển thành công phần tử tam giác 3 nút ES-MITC3 tích hợp 7 bậc tự do trên mỗi nút, loại bỏ hoàn toàn hiện tượng khóa cắt ở tỷ số t/L đạt 10 mũ trừ 4.
  • Ứng dụng xuất sắc lý thuyết biến dạng cắt bậc cao của Reddy giúp mô phỏng chính xác ứng suất cắt mà không cần dùng hệ số hiệu chỉnh cắt.
  • Kiểm chứng độ chính xác và độ hội tụ của mô hình số qua 5 cấp độ lưới từ 2x2 đến 32x32 phần tử và 6 mode dao động riêng với sai số dưới 1,5%.
  • Đóng góp quy trình tính toán số tối ưu trên MATLAB giúp giảm hơn 30% khối lượng tính toán so với mô hình đàn hồi 3 chiều truyền thống.
  • Định hình lộ trình nghiên cứu tiếp theo trong 12 tháng tới để kết hợp phần tử ES-MITC3 với lý thuyết lớp thông minh Layer-wise nhằm nâng cao hơn nữa độ chính xác ứng suất tiếp xúc.

Hãy tải ngay toàn văn luận văn và liên hệ nhóm nghiên cứu để ứng dụng công nghệ mô phỏng tiên tiến này vào các dự án kết cấu composite thực tế của bạn.