Tổng quan nghiên cứu

Trong bối cảnh công nghiệp cơ khí và hàng không hiện đại, nhu cầu giảm trọng lượng kết cấu nhưng vẫn đảm bảo khả năng chịu lực đóng vai trò then chốt. Hơn 60% các bộ phận máy bay và phương tiện vận tải thế hệ mới đang tận dụng vật liệu composite nhằm tối ưu hóa tải trọng bản thân. Kể từ khi xuất hiện lần đầu vào năm 1871 tại Mỹ, tấm các tông lõi lượn sóng đã trở thành giải pháp vật liệu dạng sandwich vượt trội cho ngành bao bì công nghiệp, xây dựng dân dụng và bảo vệ hàng hóa nhờ đặc tính nhẹ, chi phí thấp và khả năng tái chế cao.

Tuy nhiên, thách thức lớn nhất trong việc thiết kế kết cấu composite lõi sóng là việc dự đoán chính xác ứng xử cơ học phi tuyến khi chịu tải trọng phức tạp. Phần lớn các công bố khoa học trong và ngoài nước trước đây chỉ tập trung vào ứng xử đàn hồi tuyến tính đơn giản. Trong thực tế vận hành, các tác động ngoại lực như va đập mạnh, nén ép liên tục hoặc kéo uốn thường đẩy vật liệu vượt quá giới hạn đàn hồi, dẫn đến các biến dạng dẻo cục bộ nghiêm trọng làm phá hủy kết cấu.

Luận văn thạc sĩ kỹ thuật cơ khí của tác giả Nguyễn Thị Thu Phương tập trung nghiên cứu, tính toán và phát triển mô hình đồng nhất hóa hai chiều tương đương thay thế cho mô hình chi tiết ba chiều của tấm composite lõi lượn sóng chịu kéo và uốn. Mục tiêu trọng tâm là mô tả chính xác trạng thái biến dạng dẻo phi tuyến trên mẫu các tông sóng loại B, đồng thời tối ưu hóa thời gian tính toán mô phỏng số. Kết quả nghiên cứu có ý nghĩa thực tiễn to lớn khi giúp giảm hơn 94% thời gian phân tích trên hệ thống máy tính chuyên dụng, trong khi vẫn duy trì sai số dự đoán cơ học dưới ngưỡng 2% so với mô hình ba chiều nguyên bả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 được thiết lập dựa trên sự kết hợp chặt chẽ giữa ba nền tảng cơ học kết cấu: lý thuyết tấm biến dạng cắt bậc nhất Mindlin, mô hình dẻo tương đương IPE của Mäkelä và Östlund, cùng lý thuyết tấm nhiều lớp mở rộng cho cấu trúc sandwich rỗng. Vật liệu giấy cấu thành tấm sóng thể hiện rõ tính chất trực hướng với 9 hằng số đàn hồi độc lập cần được xác định trong không gian ba chiều. Theo các quan sát thực nghiệm cơ bản, độ giãn dài giới hạn đàn hồi của giấy chỉ đạt mức khoảng 0,2%, sau ngưỡng này vật liệu sẽ chuyển dịch sang miền phản ứng dẻo phi tuyến rõ rệt. Đồng thời, mô đun đàn hồi theo phương chiều dày thường thấp hơn xấp xỉ 200 lần so với phương máy chạy.

Bốn khái niệm cốt lõi định hình khung lý thuyết bao gồm: tính trực hướng vật liệu, tenxơ ứng suất lệch, mặt chảy dẻo phi tuyến và ma trận độ cứng tương đương. Bằng cách thiết lập hàm thế dẻo phi bậc hai và tích phân các quan hệ ứng suất - biến dạng qua chiều dày từng lớp thành phần, mô hình đã chuyển đổi trạng thái ứng suất phẳng phức tạp sang một hệ phương trình tiếp tuyến vi mô nhất quán, tạo nền tảng vững chắc cho quá trình đồng nhất hóa cấu trúc.

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

Quy trình nghiên cứu áp dụng phương pháp phân tích phần tử hữu hạn số kết hợp thuật toán đồng nhất hóa giải tích. Cỡ mẫu nghiên cứu được thiết lập bao gồm 532 phần tử vỏ bốn nút S4R cho mẫu thử kéo phẳng đơn lớp có chiều dày tiêu chuẩn 0,51 mm. Đối với thí nghiệm uốn ba điểm, quy mô mẫu nghiên cứu là 1330 phần tử tương đương và 1420 nút đối với mô hình đồng nhất hai chiều, đối chứng trực tiếp với 10990 phần tử và 10437 nút đối với mô hình chi tiết ba chiều trên mẫu có kích thước thực tế 140 mm chiều dài và 40 mm chiều rộng.

Phương pháp chọn mẫu tuân thủ nghiêm ngặt quy cách chế tạo tấm sóng công nghiệp loại B, lấy mẫu kiểm tra trực giao theo ba hướng chịu lực điển hình: phương máy chạy, phương ngang và phương nghiêng 45 độ nhằm ghi nhận đầy đủ hiệu ứng dị hướng. Lý do lựa chọn phương pháp phân tích tích hợp chương trình con UGENS trên Abaqus/Standard và VUMAT trên Abaqus/Explicit xuất phát từ nhu cầu giải quyết triệt để nút thắt cổ chai về tài nguyên tính toán. Phương pháp này cho phép cập nhật ma trận tiếp tuyến cục bộ tại 3 điểm tích phân Gauss trên mỗi lớp mà không cần xây dựng toàn bộ hình học ba chiều phức tạp. Toàn bộ chu trình nghiên cứu được triển khai tuần tự qua 3 giai đoạn: xây dựng thuật toán giải tích, lập trình chương trình con mô phỏng, và kiểm chứng chéo với dữ liệu thực nghiệm kéo và uốn.

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

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

Nghiên cứu đã mang lại 4 phát hiện quan trọng có giá trị khoa học và ứng dụng thực tiễn cao:

Thứ nhất, mô hình đồng nhất hóa hai chiều đạt độ chính xác tương đồng ấn tượng với mô hình ba chiều hoàn chỉnh khi chịu kéo. Sai lệch tương đối cực đại giữa hai mô hình trên các đường cong lực kéo và độ giãn dài theo cả ba phương máy chạy, phương ngang và phương nghiêng 45 độ đều duy trì ở mức dưới 2%.

Thứ hai, tốc độ giải số được nâng cao đột biến. Trong bài toán kéo đơn giản, thời gian tính toán CPU của mô hình đồng nhất hóa giảm từ 2 đến 4 lần so với mô hình ba chiều. Đặc biệt, trong bài toán mô phỏng uốn ba điểm, thời gian CPU của mô hình ba chiều là 3541 giây, trong khi mô hình đồng nhất hai chiều chỉ tiêu tốn 181 giây, tức rút ngắn thời gian tính toán nhanh hơn 19,56 lần, tương đương mức giảm 94,89% chi phí tài nguyên máy tính.

Thứ ba, cấu trúc lưới phần tử hữu hạn được tối ưu hóa vượt trội. Mô hình đồng nhất hai chiều chỉ yêu cầu 1330 phần tử và 1420 nút, giúp cắt giảm 87,9% số lượng phần tử và 86,4% số lượng nút so với cấu trúc 10990 phần tử của mô hình ba chiều.

Thứ tư, thuật toán dẻo tương đương IPE phản ánh xuất sắc trạng thái dẻo phi tuyến của vật liệu giấy các tông sóng, mô tả chính xác điểm uốn chuyển tiếp và vùng biến dạng dẻo tương thích hoàn toàn với các số liệu thực nghiệm uốn độc lập.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự cắt giảm mạnh mẽ thời gian CPU mà không làm suy giảm độ chính xác nằm ở cơ chế tích phân đại số độ cứng trực tiếp qua chiều dày tấm. Thay vì phải giải quyết hàng chục nghìn phương trình cân bằng phi tuyến tiếp xúc giữa các vách sóng ba chiều rỗng, mô hình hai chiều tương đương đại diện cho chu kỳ sóng đơn vị bằng một mặt phẳng trung bình có xét đến biến dạng cắt ngang theo lý thuyết Mindlin. So với các công bố trước đây vốn chỉ dừng lại ở phạm vi đàn hồi tuyến tính, nghiên cứu này đã tạo ra bước tiến quan trọng khi mở rộng thành công sang miền phi tuyến dẻo nhờ ma trận tiếp tuyến trong chương trình con UGENS.

Về phương diện trực quan hóa, toàn bộ hệ dữ liệu nghiên cứu được biểu diễn thông qua các biểu đồ so sánh đường cong ứng suất - biến dạng uốn, đồ thị quan hệ lực tác dụng và độ võng tâm, cùng biểu đồ cột thể hiện mức tiêu hao thời gian CPU. Một bảng dữ liệu đối chuẩn chi tiết đã xác nhận độ võng tâm mẫu uốn kích thước 140 mm nhân 40 mm có sai lệch cực đại chỉ xấp xỉ 3,5% giữa thực nghiệm và mô phỏng số. Điều này chứng minh mô hình tương đương hoàn toàn đủ độ tin cậy để thay thế các thử nghiệm uốn ba điểm thực tế trên giá đỡ bán kính 10 mm.

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

Dựa trên những đóng góp khoa học của luận văn, 4 giải pháp ứng dụng và mở rộng được khuyến nghị triển khai cụ thể:

Một là, mở rộng thư viện phần tử đồng nhất hóa sang các kết cấu lõi composite dạng sóng tam giác, sóng hình thang và cấu trúc tổ ong sandwich. Tích hợp và chuẩn hóa các thuật toán tiếp tuyến phi tuyến nhằm đạt độ chính xác sai số dưới 3% trên toàn bộ miền biến dạng dẻo trong khung thời gian 6 tháng, do nhóm nghiên cứu cơ học vật liệu và tính toán kết cấu tại các viện cơ học kỹ thuật đảm nhận.

Hai là, phát triển gói công cụ tự động hóa trên nền tảng phần mềm kỹ thuật số thương mại. Lập trình giao diện tự động xuất ma trận vật liệu cho chương trình con UGENS nhằm cắt giảm 80% thời gian thiết lập tiền xử lý hình học cho kỹ sư thiết kế trong vòng 9 tháng, do nhóm kỹ sư phần mềm tính toán mô phỏng công nghiệp chủ trì.

Ba là, tiến hành kiểm nghiệm mô hình dưới các điều kiện tải trọng động lực học phức tạp. Thực hiện mô phỏng va đập vận tốc cao và rung động ngẫu nhiên nhằm đảm bảo sai lệch kết quả gia tốc không vượt quá 5% so với đo đạc cảm biến thực tế trong thời hạn 12 tháng, do phòng thí nghiệm cơ học thực nghiệm và kiểm định kết cấu hàng không thực hiện.

Bốn là, thương mại hóa quy trình thiết kế bao bì bảo vệ thông minh phục vụ logistics xanh. Áp dụng phương pháp đồng nhất hóa vào dây chuyền sản xuất hộp chịu lực xuất khẩu để tiết kiệm tối thiểu 15% lượng giấy tiêu hao và giảm 30% chi phí mẫu thử vật lý hàng năm trong lộ trình 18 tháng, do bộ phận nghiên cứu và phát triển sản phẩm tại các doanh nghiệp sản xuất bao bì công nghiệp triển khai.

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

Luận văn thạc sĩ này mang lại giá trị tham khảo chuyên môn sâu sắc cho 4 nhóm đối tượng trọng tâm:

Thứ nhất, các kỹ sư tính toán kết cấu và chuyên gia mô phỏng cơ khí trong ngành công nghiệp ô tô và hàng không. Tài liệu cung cấp quy trình số hóa tấm composite lõi sóng, giúp tiết kiệm hơn 90% thời gian chạy phần mềm phân tích va chạm và tối ưu hóa độ bền phương tiện giao thông.

Thứ hai, giảng viên, nhà nghiên cứu và học viên cao học chuyên ngành Cơ học kỹ thuật và Khoa học vật liệu. Bản luận văn đóng vai trò tài liệu tham khảo giá trị về phương pháp biến dạng cắt Mindlin kết hợp mô hình dẻo IPE để phát triển các đề tài mở rộng về kết cấu sandwich cho hơn 10 dạng hình học lõi khác nhau.

Thứ ba, các giám đốc kỹ thuật và kỹ sư R&D trong ngành công nghiệp sản xuất bao bì giấy sóng. Luận văn cung cấp công cụ dự báo chính xác độ bền uốn và nén thùng hộp khi xếp chồng qua hàng trăm chu kỳ bảo quản kho vận thực tế mà không cần tốn kém chi phí thử phá hủy mẫu.

Thứ tư, các chuyên gia thiết kế nội thất xanh và kiến trúc ứng dụng vật liệu siêu nhẹ. Kết quả tính toán giúp tạo lập các module nội thất chịu lực từ bìa sóng có khả năng nâng đỡ an toàn tải trọng trên 100 kg mà vẫn đảm bảo tính thẩm mỹ và độ bền cơ lý.

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

Phương pháp đồng nhất hóa hai chiều giúp tiết kiệm thời gian tính toán bao nhiêu so với mô hình ba chiều truyền thống? Trong bài toán uốn ba điểm, mô hình đồng nhất hai chiều hoàn thành mô phỏng chỉ trong 181 giây so với 3541 giây của mô hình ba chiều nguyên bản. Phương pháp này giúp rút ngắn thời gian tính toán nhanh hơn 19,56 lần, tương đương mức tiết giảm đến 94,89% chi phí tài nguyên tính toán trên hệ thống máy tính chuyên dụng.

Vì sao bắt buộc phải nghiên cứu ứng xử dẻo phi tuyến thay vì dừng lại ở mô hình đàn hồi tuyến tính? Giấy và bìa cứng có miền làm việc đàn hồi rất hẹp với giới hạn độ giãn dài chỉ khoảng 0,2%. Khi chịu tải trọng uốn, kéo lớn hoặc va chạm trong thực tế, vật liệu sẽ lập tức rơi vào miền biến dạng dẻo phi tuyến, khiến mô hình thuần đàn hồi trước đây hoàn toàn mất độ chính xác.

Chương trình con nào của phần mềm Abaqus được sử dụng để tích hợp mô hình dẻo tương đương IPE? Nghiên cứu ứng dụng chương trình con VUMAT trên nền Abaqus/Explicit cho bài toán kéo tấm phẳng đơn lớp, đồng thời lập trình chương trình con người dùng UGENS trên Abaqus/Standard để giải quyết liên kết ma trận tiếp tuyến độ cứng 6x6 của mô hình đồng nhất hóa hai chiều tương đương.

Độ chính xác của mô hình đồng nhất hai chiều có được bảo chứng qua các kiểm định thực nghiệm không? Mô hình đạt độ chuẩn xác vượt trội khi đối chuẩn với thực nghiệm uốn ba điểm và mô phỏng ba chiều chi tiết, ghi nhận mức sai lệch lực kéo theo các hướng máy chạy, phương ngang và góc 45 độ dưới 2%, cùng sai số biến dạng uốn thực tế nhỏ hơn 3,5%.

Mô hình đồng nhất hóa này có thể áp dụng cho các cấu trúc lõi sandwich composite nào khác? Phương pháp toán học hoàn toàn có thể mở rộng cho các tấm composite lõi sóng hình tam giác, hình thang hoặc cấu trúc lõi tổ ong chịu tải uốn xoắn phức tạp, bằng cách thay đổi hàm số hình học mô tả đường sinh chu kỳ lõi và cập nhật các hằng số vật liệu trực hướng tương ứng.

Kết luận

  • Xây dựng thành công mô hình đồng nhất hóa hai chiều tương đương cho tấm composite lõi lượn sóng chịu kéo và uốn có xét đến ứng xử biến dạng dẻo phi tuyến.
  • Ứng dụng xuất sắc mô hình dẻo IPE và tích hợp trọn vẹn thông qua các chương trình con UGENS và VUMAT trong môi trường tính toán Abaqus.
  • Đạt độ tương thích và chuẩn xác cao với sai số dự báo dưới 2% ở bài toán kéo ba hướng và dưới 3,5% ở thí nghiệm uốn ba điểm.
  • Giảm tải vượt bậc 94,89% thời gian CPU máy tính từ 3541 giây xuống 181 giây và cắt giảm 87,9% tổng số lượng phần tử hữu hạn.
  • Đặt nền móng vững chắc cho việc mô phỏng tối ưu hóa các kết cấu tấm sandwich siêu nhẹ chịu tải va chạm phức tạp trong tương lai.

Đóng góp chính của luận văn là đã giải quyết triệt để bài toán dung hòa giữa độ chính xác vật lý phi tuyến và hiệu năng giải số cho kết cấu tấm composite lõi sóng. Kế hoạch tiếp theo trong vòng 12 đến 24 tháng tới là mở rộng kiểm nghiệm cho các bài toán va chạm động học vận tốc cao và rung động ngẫu nhiên. Quý độc giả, kỹ sư và các đơn vị sản xuất quan tâm có thể tiếp cận trọn vẹn công trình nghiên cứu và liên hệ trực tiếp với nhóm tác giả để nhận chuyển giao mô hình số tiên tiến này vào thực tiễn kỹ thuật.