Tổng quan nghiên cứu

Vật liệu composite nhiều lớp đang trở thành giải pháp vật liệu tiên tiến then chốt trong kỹ thuật kết cấu hiện đại, đóng góp tới hơn 50% khối lượng cấu trúc trên các dòng tàu hỏa cao tốc, máy bay thương mại thế hệ mới và các công trình dân dụng chịu tải trọng khắc nghiệt. Nhờ tỷ số độ bền trên khối lượng vượt trội và khả năng tùy biến cơ tính theo phương dệt sợi, kết cấu dầm composite ngày càng được ứng dụng rộng rãi. Tuy nhiên, khi vận hành trong môi trường đồng thời chịu lực nén dọc trục, nhiệt độ tăng cao lên tới 3000 độ C và độ ẩm bão hòa, kết cấu rất dễ bị mất ổn định cục bộ và suy giảm nghiêm trọng tần số dao động riêng.

Vấn đề nghiên cứu trọng tâm của luận văn là xây dựng mô hình toán học giải tích chính xác nhằm phân tích đáp ứng dao động tự do và tải trọng mất ổn định tới hạn của dầm composite nhiều lớp dưới tác động đồng thời của ba trường lực: cơ học, nhiệt độ và độ ẩm. Mục tiêu cụ thể là thiết lập hệ phương trình năng lượng tổng quát, áp dụng phương pháp Ritz kết hợp hàm dạng đa thức bậc cao để giải quyết bài toán ứng xử phi tuyến đối với đa dạng điều kiện liên kết biên.

Nghiên cứu được hoàn thành vào năm 2017 tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh, khảo sát toàn diện 5 nhóm vật liệu composite trực hướng với tỷ số mô đun đàn hồi kéo nén dao động từ 10 đến 50 và tỷ lệ kích thước chiều dài trên chiều cao từ 5 đến 50. Ý nghĩa thực tiễn của công trình thể hiện qua việc cung cấp bộ công cụ tính toán giải tích chuẩn xác, giúp tối ưu hóa kết cấu dầm composite, giảm thiểu sai số thiết kế xuống dưới 1% so với nghiệm phần tử hữu hạn 3D và nâng cao độ an toàn chịu lực trong môi trường biến đổi nhiệt ẩm khắc nghiệt.

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 ứng dụng lý thuyết dầm biến dạng cắt bậc cao nhằm vượt qua các giới hạn cố hữu của lý thuyết dầm cổ điển Euler-Bernoulli và lý thuyết biến dạng cắt bậc nhất Timoshenko. Lý thuyết biến dạng cắt bậc cao giả thiết trường chuyển vị dọc trục biến thiên phi tuyến bậc ba theo chiều dày dầm, qua đó thỏa mãn điều kiện ứng suất cắt bằng 0 tại mặt trên và mặt dưới của tiết diện mà không cần sử dụng hệ số hiệu chỉnh cắt.

Mô hình nghiên cứu xem xét dầm composite trực hướng có chiều dài L, chiều rộng b và chiều cao h, bao gồm nhiều lớp sợi xếp chồng với các góc xoay sợi khác nhau. Trạng thái cơ học của từng lớp được mô tả qua ma trận độ cứng giảm trong hệ tọa độ tổng thể với 9 hằng số đàn hồi độc lập. Hệ phương trình động lực học tổng quát được xây dựng dựa trên nguyên lý cực trị thế năng kết hợp phương trình vi phân Lagrange. Các khái niệm cốt lõi được định nghĩa chặt chẽ bao gồm: trường chuyển vị phi tuyến, ma trận độ cứng liên kết uốn kéo dãn, ten-xơ biến dạng cắt ngang, công ngoại lực do ứng suất nhiệt và ứng suất ẩm sinh ra.

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

Nguồn dữ liệu của luận văn được tổng hợp và chuẩn hóa từ hệ thống dữ liệu thực nghiệm và kết quả công bố của các nhà nghiên cứu đầu ngành trên thế giới. Phương pháp phân tích chủ đạo là phương pháp biến phân Ritz sử dụng hàm dạng đa thức trực giao kết hợp kỹ thuật nhân tử Lagrange để xử lý nghiêm ngặt các điều kiện biên hình học phức tạp như tựa đơn, khớp-khớp, ngàm-ngàm, ngàm-khớp và ngàm-tự do.

Cỡ mẫu khảo sát bao gồm 10 bài toán mô phỏng số chi tiết với chuỗi thử nghiệm số lượng đa thức m biến thiên từ 6 đến 18 trên 5 hệ vật liệu composite cốt sợi thủy tinh và carbon tiêu chuẩn. Phương pháp chọn mẫu tập trung vào việc quét toàn diện các thông số hình học và cơ lý tính nhạy cảm. Lý do lựa chọn phương pháp Ritz kết hợp hàm đa thức là khả năng thiết lập ma trận phần tử giải tích tự động trên phần mềm tính toán kỹ thuật, đảm bảo độ ổn định số cao, tốc độ hội tụ nhanh vượt trội và không gặp hiện tượng khóa cắt thường thấy trong phân tích số thông thường.

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

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

Thứ nhất, phân tích số chứng minh tính hội tụ tối ưu của chuỗi hàm đa thức Ritz. Với số số hạng đa thức m = 12, tải trọng ổn định tới hạn chuẩn hóa của dầm đã hoàn toàn ổn định, và khi nâng lên m = 14, tần số dao động tự nhiên vô thứ nguyên của dầm tựa đơn hướng sợi đối xứng đạt giá trị hội tụ chính xác 9.6032 tại tỷ lệ L/h = 5 và tỷ số E1/E2 = 40.

Thứ hai, cấu hình sắp xếp lớp sợi quyết định trực tiếp đến độ cứng tổng thể của dầm. Kết cấu dầm có lớp sợi đối xứng (0 độ / 90 độ / 0 độ) thể hiện khả năng chịu lực nén tới hạn và tần số dao động cao hơn từ 128% đến 315% so với cấu hình bất đối xứng (0 độ / 90 độ) trong cùng điều kiện hình học và tải trọng. Khi góc xoay hướng sợi tăng dần từ 0 độ lên 90 độ, tần số dao động vô thứ nguyên giảm mạnh liên tục từ 4.6201 xuống còn 0.2617 đối với điều kiện biên ngàm-ngàm.

Thứ ba, sự gia tăng của tỷ số mô đun đàn hồi E1/E2 từ 10 lên 40 làm tăng tải trọng tới hạn không thứ nguyên từ mức 4.945 lên 8.906 đối với dầm tựa đơn. Khi tỷ lệ kích thước L/h tăng từ 5 lên 50, tần số dao động cơ bản của dầm đối xứng tăng thêm khoảng 23.4% do giảm thiểu ảnh hưởng của hiệu ứng biến dạng cắt ngang.

Thứ tư, tác động của trường nhiệt độ và độ ẩm làm suy giảm nghiêm trọng khả năng làm việc của kết cấu. Khi nhiệt độ môi trường tăng cao, lực nén nhiệt dọc trục phát triển nhanh chóng, khiến tần số dao động riêng của dầm sụt giảm tuyến tính về 0 tại điểm nhiệt độ tới hạn, dẫn đến trạng thái mất ổn định nhiệt hoàn toàn.

Thảo luận kết quả

Nguyên nhân căn bản của sự chênh lệch đáp ứng giữa các cấu hình sợi nằm ở ma trận độ cứng uốn kéo dãn B. Ở dầm composite bất đối xứng, hiện tượng ghép nối cơ học giữa kéo dãn và uốn làm suy giảm mô men kháng uốn hiệu dụng, dẫn đến suy giảm tải trọng tới hạn. Khi góc sợi dịch chuyển về 90 độ, phương sợi chính vuông góc với trục dọc dầm làm dầm mất đi phần lớn độ cứng uốn dọc trục.

Khi so sánh với các nghiên cứu giải tích của Vo và Thai, Murthy cùng cộng sự, kết quả tính toán của luận văn cho thấy sự trùng khớp rất cao với sai số trung bình dưới 0.5%, chứng minh độ tin cậy tuyệt đối của thuật toán Ritz đa thức. Trên các biểu đồ quan hệ giữa nhiệt độ và tần số dao động riêng, các đường cong suy giảm dốc đứng thể hiện rõ cơ chế dầm bị nén trước do giãn nở nhiệt bị cản trở bởi các liên kết ngàm cứng, làm sáng tỏ bản chất của bài toán dao động phi tuyến trong môi trường cơ nhiệt ẩm đồng thời.

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

Thứ nhất, kỹ sư kết cấu cần áp dụng ngay quy trình bố trí lớp sợi đối xứng (0 độ / 90 độ / 0 độ) hoặc ưu tiên góc dệt sợi dọc trục 0 độ cho các cấu kiện dầm chịu lực uốn nén chính. Giải pháp này giúp gia tăng hơn 150% khả năng chịu tải mất ổn định cục bộ mà không làm phát sinh thêm trọng lượng vật liệu.

Thứ hai, các đơn vị thiết kế phải thiết lập hệ số an toàn nhiệt ẩm riêng biệt trong hồ sơ kỹ thuật. Cần kiểm soát chặt chẽ nhiệt độ làm việc dưới ngưỡng 70% nhiệt độ tới hạn tính toán và bổ sung lớp phủ ngăn ẩm bề mặt để triệt tiêu lực nén dọc trục ký sinh, duy trì tần số dao động riêng luôn cao hơn 15% so với tần số kích thích của môi trường xung quanh.

Thứ ba, các kỹ sư R&D trong ngành hàng không và xây dựng nên tích hợp thuật toán Ritz đa thức vào phần mềm tính toán sơ bộ giai đoạn đầu dự án. Việc này cho phép rút ngắn thời gian phân tích tối ưu hóa đa thông số từ vài tuần xuống còn vài giờ với độ chính xác tương đương mô hình phần tử hữu hạn 3D phức tạp.

Thứ tư, nhà sản xuất composite cần chuẩn hóa quy trình kiểm định chất lượng mô đun đàn hồi dọc trục E1 và hệ số giãn nở nhiệt ẩm của từng lô sợi trước khi gia công hàng loạt. Mục tiêu duy trì sai số cơ tính vật liệu dưới 3% trong suốt vòng đời dự án kéo dài từ 20 đến 50 năm.

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

Nghiên cứu này mang lại giá trị học thuật và ứng dụng thực tiễn to lớn cho bốn nhóm đối tượng chuyên môn sau:

Nhóm thứ nhất là các nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Xây dựng, Cơ kỹ thuật và Cơ học Vật liệu. Luận văn cung cấp toàn bộ quy trình giải tích mẫu mực từ thiết lập thế năng Lagrange đến kỹ thuật lập trình ma trận Ritz trên môi trường tính toán số.

Nhóm thứ hai là các kỹ sư kết cấu công trình dân dụng, giao thông và cầu cảng. Tài liệu là cẩm nang hữu ích để tính toán chính xác khả năng chịu lực của dầm xà gồ, sàn composite làm việc tại các vùng ven biển có độ ẩm bão hòa và bức xạ nhiệt cao.

Nhóm thứ ba là các chuyên gia thiết kế trong ngành công nghiệp đóng tàu, hàng không vũ trụ và quân sự. Các bài toán số trong luận văn đóng vai trò là dữ liệu tham chiếu chuẩn để thiết kế khung vỏ tàu thuyền composite, cánh máy bay và ống phóng khí tài.

Nhóm thứ tư là các nhà quản lý kỹ thuật và chuyên gia thẩm tra dự án. Luận văn cung cấp các tiêu chí định lượng và ngưỡng giới hạn an toàn nhiệt ẩm rõ ràng, hỗ trợ hiệu quả công tác nghiệm thu thiết kế kết cấu vật liệu mới.

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

Lý thuyết biến dạng cắt bậc cao HOBT có ưu điểm vượt trội gì so với lý thuyết Euler-Bernoulli và Timoshenko? Lý thuyết dầm biến dạng cắt bậc cao mô tả chính xác sự phân bố phi tuyến parabol của ứng suất cắt dọc theo chiều cao dầm composite. Do đó, phương pháp thỏa mãn hoàn toàn điều kiện biên ứng suất cắt bằng 0 tại bề mặt tự do mà không cần dùng hệ số hiệu chỉnh cắt nhân tạo, hạn chế tối đa sai số lên tới 25% khi tính dầm dày.

Tại sao số số hạng đa thức m = 14 được chọn làm chuẩn cho toàn bộ tính toán? Qua khảo sát thực nghiệm số với m chạy từ 6 đến 18, kết quả cho thấy tải trọng ổn định tới hạn hội tụ hoàn toàn từ m = 12 và tần số dao động hội tụ tuyệt đối tại m = 14 với độ biến thiên kết quả nhỏ hơn 0.001%. Cỡ mẫu m = 14 đảm bảo tối ưu hóa giữa độ chính xác và tài nguyên máy tính tính toán.

Hướng xếp sợi ảnh hưởng như thế nào đến tần số dao động tự nhiên của dầm composite? Góc sợi càng nhỏ so với trục dầm thì độ cứng kháng uốn càng cao, dẫn đến tần số dao động càng lớn. Khi góc xoay sợi tăng từ 0 độ lên 90 độ, tần số dao động vô thứ nguyên của dầm ngàm-ngàm giảm mạnh từ 4.6201 xuống 0.2617, tương đương mức suy giảm độ cứng uốn hơn 94%.

Nhiệt độ và độ ẩm tác động đồng thời lên dầm composite theo cơ chế vật lý nào? Sự gia tăng nhiệt độ và hấp thụ độ ẩm làm giãn nở mạng lưới phân tử polymer, suy giảm trực tiếp mô đun đàn hồi của vật liệu. Đồng thời, liên kết biên ngàm cản trở sự giãn nở, sinh ra lực nén dọc trục nội tại rất lớn, làm triệt tiêu độ cứng uốn và kéo tần số dao động về 0 tại điểm mất ổn định.

Điều kiện biên liên kết nào mang lại khả năng chịu lực và chống mất ổn định tốt nhất? Điều kiện biên ngàm hai đầu C-C luôn cho tần số dao động tự nhiên và tải trọng ổn định tới hạn cao nhất nhờ khả năng kiềm chế hoàn toàn chuyển vị xoay và chuyển vị thẳng tại hai đầu dầm, theo sau lần lượt là ngàm-khớp C-H, khớp-khớp H-H, tựa đơn S-S và thấp nhất là ngàm-tự do C-F.

Kết luận

Luận văn đã giải quyết trọn vẹn bài toán cơ học phức tạp về đáp ứng dao động và ổn định của dầm composite nhiều lớp dưới tác động đồng thời của trường tải trọng cơ, nhiệt độ và độ ẩm.

Năm đóng góp khoa học và kỹ thuật cốt lõi của công trình bao gồm:

  • Xây dựng thành công mô hình giải tích chính xác dựa trên lý thuyết biến dạng cắt bậc cao kết hợp phương trình Lagrange tổng quát.
  • Ứng dụng xuất sắc phương pháp Ritz và hàm dạng đa thức với kỹ thuật nhân tử Lagrange, giải quyết triệt để bài toán biên phức tạp.
  • Xác lập chuẩn hội tụ tin cậy tại m = 14 với độ chính xác cao và sai số dưới 0.5% so với các công bố quốc tế uy tín.
  • Định lượng hóa toàn diện quy luật suy giảm tần số và tải trọng tới hạn theo góc xoay sợi, tỷ lệ hình học L/h và hệ số đàn hồi E1/E2.
  • Làm sáng tỏ cơ chế mất ổn định nhiệt ẩm đồng thời, thiết lập cơ sở dữ liệu quan trọng cho thiết kế dầm composite làm việc trong môi trường khắc nghiệt.

Kế hoạch nghiên cứu tiếp theo cần mở rộng mô hình cho dầm composite có tiết diện thay đổi, xét đến tính chất phi tuyến hình học lớn và phân bố nhiệt ẩm không đều trong giai đoạn 2024-2026. Hãy áp dụng ngay các công thức và giải pháp tính toán từ luận văn này vào quy trình thiết kế kết cấu composite của bạn để nâng cao tối đa hiệu quả chịu lực và độ an toàn công trình!