Tổng quan nghiên cứu

Trong hơn 20 năm qua, sự phát triển vượt bậc của khoa học vật liệu đã tạo ra các dạng kết cấu composite tiên tiến, nổi bật là vật liệu có cơ tính biến thiên (Functionally Graded Material - FGM). Khác với vật liệu phân lớp truyền thống dễ bị tách lớp do tập trung ứng suất tại mặt tiếp giáp, tấm FGM sở hữu thành phần cơ tính biến đổi liên tục từ gốm sang kim loại, giúp nâng cao khả năng chịu nhiệt độ cao lên tới hơn 1200 độ C và tăng độ bền chống ăn mòn bề mặt thêm 45%. Tuy nhiên, khi làm việc trong điều kiện tải trọng động và môi trường nền móng phức tạp, việc phân tích dao động và độ võng của tấm FGM gia cường gân vẫn là một thách thức kỹ thuật lớn đòi hỏi các phương pháp tính toán số có độ chính xác cao.

Mục tiêu cụ thể của nghiên cứu là xây dựng mô hình phần tử hữu hạn dựa trên lý thuyết biến dạng trượt bậc nhất để giải quyết bài toán tĩnh và bài toán dao động cưỡng bức của tấm chữ nhật FGM có gân gia cường đặt trên nền đàn hồi Winkler. Nghiên cứu được hoàn thành vào năm 2018 tại Khoa Toán - Cơ - Tin học, Trường Đại học Khoa học Tự nhiên thuộc Đại học Quốc gia Hà Nội. Công trình tập trung vào việc mô hình hóa các tấm kích thước tiêu chuẩn cạnh 0.4 mét, chiều dày biến thiên với các tỷ số kích thước hình học khác nhau. Ý nghĩa thực tiễn của luận văn thể hiện ở việc đưa ra giải pháp giảm độ võng kết cấu tới hơn 35% khi tối ưu hóa bố trí gân gia cường, đồng thời kiểm soát sai số tính toán dưới ngưỡng 2.5% so với nghiệm giải tích, cung cấp cơ sở khoa học tin cậy cho thiết kế kết cấu hàng không, tàu biển và xây dựng công trình hiện đạ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 ứng dụng đồng thời hai lý thuyết cơ học nền tảng là lý thuyết tấm Mindlin - Reissner (lý thuyết biến dạng trượt bậc nhất FSDT) và lý thuyết dầm Timoshenko kết hợp với mô hình nền đàn hồi Winkler. Khác với lý thuyết tấm mỏng cổ điển Kirchhoff bỏ qua biến dạng trượt ngang, lý thuyết Mindlin - Reissner xét đến ảnh hưởng của góc xoay mặt cắt với hệ số hiệu chỉnh trượt chuẩn là 5/6, cho phép tính toán chính xác ứng xử uốn của tấm có độ dày trung bình và tấm dày. Đồng thời, mô hình dầm Timoshenko được áp dụng để mô tả biến dạng trượt của hệ thống gân gia cường đối xứng.

Khung lý thuyết vận hành dựa trên 4 khái niệm trọng tâm:

  1. Quy luật biến thiên cơ tính lũy thừa (P-FGM): Mô đun đàn hồi E(z) và khối lượng riêng rho(z) biến đổi dọc theo chiều dày h theo hàm số mũ với chỉ số thể tích k chạy từ 0 đến 10.
  2. Quan hệ ứng suất - biến dạng đàn hồi theo định luật Hooke mở rộng cho vật liệu dị tính.
  3. Mô hình nền đàn hồi Winkler: Đất nền phản ứng lại chuyển vị thẳng đứng của tấm thông qua hệ số nền K với áp lực phản lực tỷ lệ thuận bậc nhất với độ võng w.
  4. Cân bằng động lực học kết cấu theo nguyên lý công ảo và phương trình Lagrange loại hai.

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

Nguồn dữ liệu tính toán được thiết lập từ các thông số vật liệu thực tế đã được công bố quốc tế, bao gồm pha gốm Silicon Nitride với mô đun đàn hồi 380 GPa, khối lượng riêng 3800 kg/m3 và pha kim loại Nhôm với mô đun đàn hồi 70 GPa, khối lượng riêng 2702 kg/m3, hệ số Poisson cố định ở mức 0.3. Timeline nghiên cứu được thực hiện trong thời gian 24 tháng với quy trình chuẩn hóa từ lập công thức giải tích đến lập trình tính toán số.

Phương pháp phân tích chủ đạo là phương pháp phần tử hữu hạn (FEM) kết hợp thuật toán tích phân từng bước Newmark trong miền thời gian. Cỡ mẫu tính toán sử dụng lưới chia 16 phần tử tấm chữ nhật 8 nút, mỗi nút gồm 3 bậc tự do (độ võng w0 và 2 góc xoay theta x, theta y) tạo nên véc-tơ chuyển vị 24 thành phần cho mỗi phần tử. Hệ thống gân gia cường được rời rạc hóa bằng phần tử dầm 3 nút với 6 bậc tự do. Lý do lựa chọn mô hình này là khả năng mô phỏng chính xác trường chuyển vị phi tuyến, kết hợp kỹ thuật tích phân số 2 điểm Gauss nhằm triệt tiêu triệt để hiện tượng khóa trượt (shear locking) thường gặp trong phần tử tấm bậc thấp. Toàn bộ thuật toán được lập trình tự động trên nền tảng phần mềm toán học Maple.

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

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

Nghiên cứu đã chỉ ra 4 phát hiện quan trọng về ứng xử tĩnh và dao động của tấm FGM gia cường gân:

  1. Ảnh hưởng của chỉ số thể tích k: Khi chỉ số k tăng từ 0 (tấm 100% gốm) lên k = 1, độ võng tĩnh cực đại tại tâm tấm tăng mạnh hơn 42% do tỷ lệ kim loại tăng làm giảm độ cứng tổng thể của tấm. Tuy nhiên, khi k tiếp tục tăng từ 1 lên 10, tốc độ tăng độ võng chậm lại đáng kể, chỉ tăng thêm khoảng 12% đến 15%.
  2. Ảnh hưởng của tỷ số kích thước hình học a/h: Khi tỷ số a/h tăng từ 10 lên 30 (tấm mỏng dần), độ võng tĩnh và biên độ dao động tại tâm tấm tăng gấp 3.2 lần do mô men kháng uốn của mặt cắt giảm theo hàm bậc ba của chiều dày.
  3. Vai trò của nền đàn hồi Winkler: Với hệ số nền K tăng từ 0 lên 0.5x10^8 N/m3, độ võng cực đại của tấm chữ nhật giảm đến 48.6%, đồng thời biên độ dao động cưỡng bức ở tần số kích thích 1650 rad/s giảm hơn 50% so với khi không đặt trên nền.
  4. Hiệu quả của cấu trúc gân gia cường: Tấm được bố trí 2 gân vuông góc chữ thập tại tâm cho độ võng nhỏ hơn 32% so với tấm chỉ có 2 gân song song và giảm hơn 56% biên độ dao động so với tấm trơn không có gân gia cường.

Thảo luận kết quả

Nguyên nhân chính của các hiện tượng trên bắt nguồn từ sự tương phản cơ tính giữa hai pha vật liệu: gốm có mô đun đàn hồi gấp 5.42 lần kim loại, đóng vai trò tạo khung chịu lực chính, trong khi kim loại duy trì tính dẻo dai. Khi gân gia cường được bố trí dạng chữ thập, mô men quán tính uốn theo cả hai phương trục x và y đều tăng vọt, giúp phân bổ tải trọng cục bộ thành tải trọng phân bố rộng khắp mặt tấm.

Khi so sánh với kết quả giải tích của các công trình cơ học kinh điển trên thế giới, sai số tương đối của mô hình phần tử hữu hạn trong nghiên cứu chỉ dao động từ 1.8% đến 2.4%, chứng minh độ tin cậy vượt trội của thuật toán. Trong các báo cáo kỹ thuật, dữ liệu này có thể được trình bày trực quan qua bảng so sánh độ võng cực đại theo các điều kiện biên bốn cạnh tựa đơn (SSSS) hoặc hai cạnh ngàm hai cạnh tựa đơn (CCSS), kết hợp với biểu đồ đường cong dao động w(t) theo tỷ số thời gian t/T để hiển thị rõ pha triệt tiêu biên độ dao động khi có lực cản của nền móng.

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

Dựa trên kết quả mô phỏng số, 4 khuyến nghị kỹ thuật có tính ứng dụng cao được đề xuất:

  1. Tối ưu hóa thành phần vật liệu chế tạo: Khống chế chỉ số thể tích k trong khoảng 0.5 đến 1.0 khi gia công tấm FGM nhằm tiết kiệm 25% khối lượng pha gốm đắt tiền mà vẫn duy trì được trên 85% độ cứng chống uốn tổng thể. Thời gian áp dụng từ năm 2026 đến 2027, do các kỹ sư vật liệu tại các viện nghiên cứu luyện kim và composite đảm trách.
  2. Tiêu chuẩn hóa thiết kế gân gia cường dạng chữ thập: Bắt buộc áp dụng sơ đồ bố trí tối thiểu 2 gân trực giao tại khu vực chịu mô men uốn lớn nhất trên các vách ngăn tàu thủy và thân vỏ máy bay để tăng 35% khả năng chịu tải va đập. Mục tiêu hoàn thiện quy chuẩn thiết kế trong vòng 12 tháng, giao cho các phòng thiết kế cơ khí chính xác thực hiện.
  3. Ứng dụng nền đàn hồi giảm chấn công trình: Tích hợp các lớp lót đệm đàn hồi composite có hệ số đàn hồi K tối thiểu 0.5x10^8 N/m3 dưới các bệ máy và sàn chịu tải rung để giảm thiểu 40% biên độ dao động cộng hưởng trước quý 4 năm 2027, do các kỹ sư kết cấu xây dựng dân dụng và công nghiệp triển khai.
  4. Nâng cấp thuật toán phần tử hữu hạn đa trường: Mở rộng chương trình tính toán số trên Maple sang các kết cấu phức tạp như vỏ trụ cong và vỏ nón FGM tích hợp ảnh hưởng của trường nhiệt độ trong 18 tháng tới, do nhóm nghiên cứu cơ học vật rắn tại các trường đại học chủ trì.

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

  1. Kỹ sư kết cấu hàng không và vũ trụ: Tận dụng các công thức phân tích dao động tấm FGM để tối ưu hóa thiết kế cánh máy bay siêu âm, vỏ tên lửa chịu tải trọng rung động cưỡng bức ở dải tần số trên 1500 rad/s.
  2. Học viên cao học và nghiên cứu sinh ngành Cơ học vật rắn: Sử dụng phương pháp xây dựng ma trận độ cứng phần tử 8 nút và thuật toán Newmark làm tài liệu tham khảo mẫu cho việc phát triển các mô hình tính toán số nâng cao.
  3. Kỹ sư xây dựng công trình ngầm và thủy công: Khai thác mô hình tấm FGM trên nền Winkler để tính toán gia cố bản đáy hầm, sàn đáy bể chứa chịu áp lực thủy tĩnh và phản lực nền móng biến thiên.
  4. Chuyên gia phát triển phần mềm mô phỏng CAE: Sử dụng thuật toán tích phân 2 điểm Gauss và phương pháp khử khóa trượt trong luận văn để lập trình các module mở rộng cho các phần mềm thương mại như Abaqus, Ansys.

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

  1. Tại sao vật liệu FGM lại vượt trội hơn vật liệu composite nhiều lớp? Vật liệu FGM có thành phần cơ tính thay đổi trơn liên tục từ gốm sang kim loại, giúp loại bỏ hoàn toàn các mặt tiếp giáp đột ngột. Điều này triệt tiêu hiện tượng tập trung ứng suất trượt cục bộ, ngăn ngừa nguy cơ tách lớp vốn làm giảm hơn 50% độ bền của kết cấu composite phân lớp truyền thống.

  2. Lý do luận văn lựa chọn lý thuyết Mindlin - Reissner thay cho lý thuyết Kirchhoff? Lý thuyết Kirchhoff cổ điển bỏ qua hoàn toàn biến dạng trượt ngang, chỉ phù hợp cho tấm rất mỏng với tỷ số a/h lớn hơn 50. Lý thuyết Mindlin - Reissner đưa vào góc xoay bổ sung và hệ số hiệu chỉnh 5/6, giúp tính toán chính xác độ võng cho cả tấm có độ dày trung bình với sai số dưới 2.5%.

  3. Thuật toán Newmark đóng vai trò gì trong việc giải bài toán dao động? Thuật toán Newmark là phương pháp tích phân số từng bước trực tiếp trong miền thời gian. Thuật toán thiết lập mối liên hệ giữa véc-tơ chuyển vị, vận tốc và gia tốc tại các bước thời gian kế tiếp, cho phép xác định chính xác đáp ứng động lực học của tấm dưới tải trọng tuần hoàn sin mà không làm suy giảm độ ổn định số.

  4. Nền đàn hồi Winkler tác động như thế nào đến độ võng của kết cấu tấm? Nền đàn hồi Winkler cung cấp một trường phản lực phân bố liên tục tỷ lệ thuận với độ võng. Khi hệ số nền K đạt 0.5x10^8 N/m3, lực kháng uốn của toàn hệ thống tăng lên, giúp hấp thụ năng lượng dao động và làm suy giảm độ võng tĩnh cực đại tới hơn 48%.

  5. Kỹ thuật tích phân 2 điểm Gauss giải quyết bài toán kỹ thuật nào trong phần tử dầm? Khi tính toán ma trận độ cứng của phần tử dầm Timoshenko, hiện tượng khóa trượt thường làm cho phần tử trở nên quá cứng giả tạo. Kỹ thuật đổi biến và tích phân số 2 điểm Gauss giúp làm mềm ma trận độ cứng trượt, khôi phục độ chính xác và đảm bảo tốc độ hội tụ cao cho lời giải.

Kết luận

  • Xây dựng thành công mô hình phần tử hữu hạn 8 nút theo lý thuyết biến dạng trượt bậc nhất để phân tích chính xác ứng xử tĩnh và dao động của tấm chữ nhật FGM.
  • Thiết lập hoàn chỉnh ma trận độ cứng, ma trận khối lượng và véc-tơ tải trọng cho hệ tấm FGM tích hợp gân gia cường dầm Timoshenko đặt trên nền đàn hồi Winkler.
  • Khảo sát toàn diện ảnh hưởng của 6 tham số then chốt, chứng minh gân chữ thập và nền đàn hồi giúp giảm hơn 50% biên độ dao động cưỡng bức.
  • Kiểm chứng thuật toán trên phần mềm Maple với độ chính xác cao, sai số so với nghiệm giải tích chỉ từ 1.8% đến 2.4%.
  • Định hướng phát triển mở rộng mô hình tính toán cho kết cấu vỏ cong và vật liệu gia cường ống nano carbon trong giai đoạn 2026 - 2028.

Các đơn vị nghiên cứu và doanh nghiệp sản xuất kết cấu composite tiên tiến hãy liên hệ ngay với nhóm tác giả tại Khoa Toán - Cơ - Tin học, Trường Đại học Khoa học Tự nhiên - ĐHQGHN để hợp tác ứng dụng và chuyển giao công nghệ mô phỏng vật liệu mới.