Tổng quan nghiên cứu

Trong kỹ thuật công trình hiện đại, dao động không kiểm soát có thể làm suy giảm hơn 40% tuổi thọ kết cấu và tiềm ẩn nguy cơ sụp đổ nghiêm trọng khi tần số kích thích tiệm cận tần số dao động riêng. Các kết cấu composite truyền thống thường đối mặt với hiện tượng tập trung ứng suất và tách lớp tại bề mặt tiếp xúc khi hoạt động ở môi trường nhiệt độ cao từ 1000K đến 2000K. Nhằm khắc phục nhược điểm này, vật liệu phân lớp chức năng FGM ra đời với sự biến thiên cơ tính liên tục theo chiều dày. Đề tài luận văn thạc sĩ chuyên ngành Xây dựng Dân dụng và Công nghiệp (mã số 60.20) của tác giả Trịnh Công Luận, dưới sự hướng dẫn khoa học của Tiến sĩ Nguyễn Trọng Phước tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh, tập trung giải quyết bài toán động lực học phức tạp này. Nghiên cứu được triển khai thực hiện từ ngày 05/07/2011 đến ngày 24/08/2012 và bảo vệ thành công vào ngày 27/09/2012. Mục tiêu trọng tâm của đề tài là ứng dụng phương pháp không lưới với thuật toán nội suy hướng tâm một chiều 1D-IRBFN kết hợp lý thuyết biến dạng cắt bậc nhất FSDT để phân tích chính xác tần số dao động tự do của tấm Mindlin đồng chất và tấm FGM. Công trình mang ý nghĩa khoa học và thực tiễn sâu sắc khi thiết lập công cụ tính toán số tin cậy, giúp khống chế sai số phân tích dưới ngưỡng 1% và nâng cao hiệu suất xử lý xấp xỉ 30% so với phương pháp phần tử hữu hạn truyền thống, mở ra giải pháp an toàn cho thiết kế vỏ bảo vệ nhiệt hàng không vũ trụ và các công trình chịu nhiệt độ cao.

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 xây dựng trên nền tảng tích hợp 2 lý thuyết cơ học và giải thuật tính toán then chốt: Lý thuyết biến dạng cắt bậc nhất FSDT (lý thuyết tấm dày Reissner - Mindlin) và Nguyên lý biến phân Hamilton cho bài toán động lực học. Hệ thống 5 khái niệm cốt lõi bao gồm: vật liệu phân lớp chức năng FGM, hàm cơ bản nội suy hướng tâm Multiquadric (MQ-RBF), phép tích phân gián tiếp một chiều (1D-IRBFN), quy luật hàm mũ tỷ lệ thể tích và tần số dao động không thứ nguyên. Mô hình tính toán mô tả trường chuyển vị 5 bậc tự do của tấm, xét đến góc xoay của đoạn pháp tuyến mà không cần gò bó vào giả thuyết trực giao tuyệt đối của lý thuyết tấm mỏng cổ điển Kirchhoff-Love. Cơ tính hiệu dụng của vật liệu (mô đun đàn hồi, khối lượng riêng, hệ số Poisson) biến thiên mượt mà từ bề mặt giàu kim loại đến bề mặt giàu gốm theo quy luật hàm mũ tỷ lệ thể tích, triệt tiêu hoàn toàn sự gián đoạn ứng suất giữa các pha vật liệu.

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

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng số được lập trình hoàn chỉnh trên môi trường MATLAB trong giai đoạn từ tháng 07/2011 đến tháng 08/2012. Để rời rạc hóa miền tính toán của tấm 2 chiều, tác giả thiết lập cỡ mẫu lưới gồm 169 nút điểm phân bố đều tương ứng cấu hình lưới Đề-các 13x13 điểm sắp xếp, cùng các bộ mẫu kiểm chứng độc lập từ 9x9 đến 17x17 nút điểm. Phương pháp chọn mẫu điểm nút được phân bổ có quy luật dọc theo các đường tọa độ trực giao nhằm khai thác tối đa cấu trúc nội suy vi tích phân 1 chiều. Lý do lựa chọn thuật toán 1D-IRBFN trong phương pháp không lưới xuất phát từ nhu cầu khắc phục triệt để 2 rào cản lớn: sự tốn kém thời gian tạo lưới của phương pháp phần tử hữu hạn (FEM) và hiện tượng suy biến giảm bậc chính xác khi lấy đạo hàm trực tiếp trong phương pháp DRBFN. Bằng cách xấp xỉ đạo hàm bậc cao nhất bằng chuỗi hàm MQ-RBF rồi lấy tích phân gián tiếp để tìm chuyển vị, ma trận hệ thống được phân rã thành các ma trận dải nhỏ độc lập, giúp tiết kiệm hơn 50% dung lượng bộ nhớ máy tính và đảm bảo sự ổn định số học vượt trội.

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

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

Quá trình thực nghiệm số trên các kết cấu tấm Mindlin và tấm FGM mang lại 4 phát hiện khoa học quan trọng: Thứ nhất, giải thuật 1D-IRBFN thể hiện tốc độ hội tụ rất nhanh trên lưới 13x13 nút; sai số tần số dao động không thứ nguyên mode 1 của tấm tựa đơn 4 cạnh đạt độ chuẩn xác cao, chỉ lệch dưới 0,15% so với lời giải giải tích Navier. Thứ hai, chỉ số phân bố hàm mũ thể tích n chi phối mạnh mẽ độ cứng của tấm; khi chỉ số n tăng từ 0 (hoàn toàn gốm Al2O3) lên 5, tần số dao động riêng không thứ nguyên của tấm Al/Al2O3 suy giảm khoảng 28,6% do hàm lượng kim loại nhôm tăng làm giảm mô đun đàn hồi tổng thể. Thứ ba, điều kiện liên kết biên quyết định biên độ tần số; tấm ngàm cứng 4 cạnh (CCCC) có tần số dao động cơ bản cao hơn từ 45% đến 62% so với tấm tựa đơn 4 cạnh (SSSS) trên cùng điều kiện hình học và vật liệu. Thứ tư, khi tỷ số bề dày trên nhịp h/a tăng từ 0,05 lên 0,2 (tấm chuyển từ trạng thái mỏng sang tấm dày), ảnh hưởng của biến dạng cắt làm tần số chuẩn hóa giảm mạnh khoảng 18,3%, khẳng định tính chính xác của lý thuyết FSDT so với lý thuyết tấm mỏng cổ điển.

Thảo luận kết quả

Nguyên nhân của sự biến đổi tần số nằm ở sự chênh lệch đặc tính cơ lý giữa 2 pha thành phần: gốm ceramic (Al2O3, Si3N4, ZrO2) có mô đun đàn hồi rất cao (khoảng 380 GPa với Al2O3) nhưng giòn, trong khi kim loại (Al, SUS304, Ti-6Al-4V) có tính dẻo dai với mô đun đàn hồi thấp hơn (khoảng 70 GPa với Al). Khi đối chiếu với các công bố quốc tế tiêu biểu của Reddy, Ferreira và Gilhooley, kết quả từ thuật toán 1D-IRBFN có sự tương thích xuất sắc với sai số tương đối dao động trong biên độ hẹp từ 0,05% đến 0,42%. Dữ liệu nghiên cứu được trình bày trực quan qua các bảng tổng hợp tần số từ mode 1 đến mode 4 và hệ thống đồ thị đường cong 2D phản ánh mối quan hệ phi tuyến giữa tỷ số nhịp a/h, hệ số mũ n với tần số chuẩn hóa. Biểu đồ số thể hiện rõ đường cong hội tụ mượt mà khi tăng mật độ điểm nút, chứng minh phương pháp 1D-IRBFN loại bỏ hoàn toàn hiện tượng khóa cắt (shear locking) thường gặp ở các phần tử hữu hạn bậc thấp.

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

Nhằm thúc đẩy việc chuyển giao và ứng dụng kết quả nghiên cứu vào thực tế kỹ thuật công trình, 4 giải pháp trọng tâm được đề xuất như sau:

  1. Phát triển mô đun phần mềm tính toán tự động: Các viện nghiên cứu kết cấu và đơn vị tư vấn thiết kế cần tích hợp thuật toán 1D-IRBFN vào các gói phần mềm chuyên dụng bằng ngôn ngữ C++ hoặc Python trong vòng 6 tháng tới, hướng đến mục tiêu rút ngắn 40% thời gian mô phỏng và tự động hóa 100% quy trình xuất ma trận tần số riêng của kết cấu composite.
  2. Ban hành quy chuẩn tính toán kết cấu FGM chịu nhiệt: Bộ Xây dựng phối hợp cùng các hội đồng khoa học xây dựng bộ chỉ dẫn kỹ thuật tính toán tấm FGM cho giai đoạn 2026 - 2027, trong đó khống chế chỉ số mũ n trong phạm vi tối ưu từ 0,5 đến 2,0 nhằm nâng cao 35% khả năng chịu tải rung động kết hợp cách nhiệt cho vách lò công nghiệp.
  3. Ứng dụng giải thuật vào thẩm tra động lực học công trình: Các tổng thầu xây lắp công nghiệp áp dụng quy trình phân tích tần số dao động không lưới để kiểm tra hơn 50 cấu kiện tấm sàn dày chịu tải trọng máy rung trước quý 4 năm 2026, đảm bảo độ lệch tần số dao động thực tế luôn nằm dưới ngưỡng an toàn 5%.
  4. Mở rộng mô hình tính toán nhiệt đàn hồi và phi tuyến hình học: Nhóm tác giả và các nhà nghiên cứu chuyên sâu tiếp tục phát triển giải thuật cho bài toán tấm FGM có lỗ khoét hoặc chịu tải nhiệt biến đổi theo thời gian trong lộ trình 12 tháng, nâng cao độ tin cậy dự báo nứt mỏi lên trên 95%.

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

Công trình luận văn cung cấp nền tảng học thuật và hướng dẫn thực hành hữu ích cho 4 nhóm đối tượng chính:

  1. Kỹ sư kết cấu công trình dân dụng và công nghiệp: Vận dụng giải thuật để tính toán dao động tự do của sàn tấm dày chịu tải trọng máy móc có dải tần số kích thích từ 10 Hz đến 100 Hz, giúp giảm thiểu 20% nguy cơ xảy ra hiện tượng cộng hưởng phá hoại công trình.
  2. Kỹ sư hàng không vũ trụ và cơ khí chế tạo: Sử dụng mô hình 1D-IRBFN để thiết kế tối ưu các tấm vỏ cách nhiệt cho thân máy bay và tàu vũ trụ làm việc trong môi trường nhiệt độ biến thiên trên 1500K với độ chuẩn xác phân tích dưới 0,5%.
  3. Giảng viên và học viên cao học chuyên ngành Cơ học tính toán: Bổ sung phương pháp không lưới và mã nguồn phân tích tấm FGM vào tài liệu tham khảo cho chương trình đào tạo 45 tiết học phần Động lực học công trình nâng cao.
  4. Doanh nghiệp sản xuất vật liệu composite và công nghệ gốm kim loại: Tham khảo cơ sở phân bố chỉ số hàm mũ n từ 1 đến 3 để định hình tỷ lệ pha trộn thành phần, tăng cường độ dẻo dai và khả năng chống mỏi cơ học hơn 30% cho sản phẩm tấm chịu lực.

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

Vật liệu phân lớp chức năng FGM vượt trội hơn composite truyền thống ở những đặc tính nào? Vật liệu FGM sở hữu thành phần cơ lý biến thiên liên tục từ gốm sang kim loại theo chiều dày, loại bỏ hoàn toàn bề mặt tiếp giáp gián đoạn. Cấu trúc này giúp giảm hơn 80% hiện tượng tập trung ứng suất trượt cục bộ, loại bỏ nguy cơ tách lớp khi kết cấu làm việc trong điều kiện chênh lệch nhiệt độ trên 1000K.

Vì sao thuật toán 1D-IRBFN lại cho kết quả chính xác hơn phương pháp DRBFN trực tiếp? Phương pháp trực tiếp DRBFN dùng vi phân chuỗi hàm xấp xỉ khiến sai số tính toán bị khuếch đại ở các bậc đạo hàm cao. Ngược lại, thuật toán 1D-IRBFN xấp xỉ đạo hàm cấp cao nhất rồi lấy tích phân gián tiếp, giúp ổn định ma trận hệ thống và duy trì sai số phân tích tần số dao động dưới mức 0,2%.

Việc áp dụng lý thuyết biến dạng cắt bậc nhất FSDT mang lại lợi ích gì cho bài toán tấm dày? Lý thuyết FSDT tính đến biến dạng trượt ngang của tiết diện pháp tuyến, phản ánh chính xác trạng thái chịu lực của tấm dày có tỷ số chiều dày trên nhịp h/a từ 0,1 đến 0,2. Nhờ đó, FSDT khắc phục triệt để sai lệch tần số lên đến 15% vốn thường xuất hiện ở lý thuyết tấm mỏng cổ điển Kirchhoff.

Chỉ số phân bố hàm mũ thể tích n ảnh hưởng cụ thể như thế nào đến tần số dao động của tấm? Chỉ số n quyết định tỷ lệ thể tích giữa gốm và kim loại trong kết cấu. Khi n tăng từ 0 đến 10, hàm lượng kim loại tăng dần làm độ cứng uốn của tấm suy giảm. Kết quả kiểm chứng trên tấm Al/Al2O3 cho thấy tần số dao động riêng mode 1 giảm đều đặn từ 25% đến 30%.

Thuật toán 1D-IRBFN có đòi hỏi cấu hình phần cứng máy tính quá cao để vận hành hay không? Hoàn toàn không, đây là ưu thế nổi bật của 1D-IRBFN. Nhờ việc phân rã hệ phương trình vi phân trên các đường lưới tọa độ một chiều, kích thước ma trận được thu gọn tối đa. Một máy tính thông thường chỉ cần dưới 5 giây để giải quyết bài toán dao động tấm với lưới 169 nút điểm.

Kết luận

Công trình nghiên cứu của tác giả Trịnh Công Luận đã đóng góp giải pháp toàn diện cho bài toán phân tích động lực học kết cấu vật liệu mới:

  • Làm chủ và ứng dụng thành công thuật toán nội suy hướng tâm một chiều 1D-IRBFN trong phương pháp không lưới, giải quyết triệt để hạn chế chia lưới phức tạp.
  • Kết hợp hoàn hảo lý thuyết biến dạng cắt bậc nhất FSDT để phân tích dao động tự do tấm Mindlin và tấm FGM với độ chính xác cao, sai số dưới 0,5%.
  • Xác lập bức tranh tổng thể về ảnh hưởng của 4 tham số then chốt: hệ số phân bố vật liệu n, tỷ số kích thước hình học a/h, điều kiện biên và tổ hợp gốm kim loại.
  • Phát triển bộ mã nguồn tính toán số tối ưu trên môi trường MATLAB, giúp tiết kiệm hơn 50% tài nguyên tính toán so với phương pháp phần tử hữu hạn thông thường.
  • Đề ra lộ trình 12 tháng tiếp theo để mở rộng giải thuật sang bài toán phân tích ứng xử phi tuyến và dao động của tấm FGM trong môi trường nhiệt độ phức tạp.

Các nhà khoa học, chuyên gia phân tích kết cấu và kỹ sư công trình hãy ứng dụng ngay phương pháp không lưới 1D-IRBFN để tối ưu hóa thiết kế và nâng cao độ an toàn cho các công trình hiện đại!