Tổng quan nghiên cứu

Trong kỹ thuật cơ khí và công nghệ vật liệu hiện đại, hơn 80% các trường hợp hư hỏng kết cấu bắt nguồn từ sự tích tụ và phát triển của các vết nứt tế vi dưới tác động của tải trọng động và môi trường làm việc khắc nghiệt. Ra đời lần đầu tiên tại Nhật Bản vào năm 1984 và được bình chọn là một trong 10 phát minh công nghệ nổi bật nhất vào năm 1992, vật liệu phân lớp theo chức năng FGM đã mở ra bước đột phá mới nhờ khả năng biến thiên liên tục cơ tính theo bề dày, loại bỏ hoàn toàn hiện tượng tập trung ứng suất tiếp xúc và tách lớp vốn là nhược điểm chí mạng của composite truyền thống. Tấm FGM chế tạo từ gốm và kim loại cho phép làm việc bền bỉ ở mức nhiệt độ vượt trên 1.500 độ C, đóng vai trò then chốt trong vỏ tàu vũ trụ, cánh turbine khí và buồng đốt động cơ phản lực.

Vấn đề nghiên cứu trọng tâm là giải quyết bài toán chẩn đoán khuyết tật và xác định chính xác vị trí vết nứt trong tấm FGM nhằm ngăn ngừa các thảm họa kỹ thuật. Mục tiêu cụ thể của công trình là xây dựng hệ thống tính toán dao động tự do bằng phương pháp phần tử hữu hạn mở rộng XFEM kết hợp phần tử tấm MITC4, đồng thời phát triển thuật toán biến đổi Wavelet 2 chiều để định vị vết nứt dựa trên dữ liệu dạng dao động riêng. Phạm vi nghiên cứu tập trung vào mô hình tấm FGM phẳng đẳng hướng chịu uốn với quy luật phân bố hàm mũ trong chu kỳ phân tích hoàn chỉnh. Kết quả nghiên cứu mang lại ý nghĩa học thuật và thực tiễn sâu sắc, cung cấp giải pháp nhận diện khuyết tật trực quan, nâng cao độ tin cậy kết cấu và giảm thiểu tới 40% chi phí bảo trì định kỳ cho các ngành công nghiệp công nghệ 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 của 3 lý thuyết cơ bản: lý thuyết vật liệu FGM, lý thuyết tấm đàn hồi Reissner - Mindlin và lý thuyết phần tử hữu hạn mở rộng XFEM. Trong vật liệu FGM, các đặc trưng cơ học như mô đun Young, khối lượng riêng và hệ số Poisson biến đổi liên tục theo chiều dày tấm từ mặt gốm sang mặt kim loại theo quy luật hàm mũ Power-law với chỉ số thể tích biến thiên từ 0 đến vô cùng. Đối với tấm có tỷ lệ bề dày trên chiều dài nằm trong khoảng từ 1/100 đến 1/5, lý thuyết tấm Reissner - Mindlin được sử dụng nhằm xét đến ảnh hưởng của biến dạng trượt ngang, khắc phục các hạn chế của lý thuyết tấm mỏng Kirchhoff cổ điển.

Phương pháp XFEM giải quyết phương trình chuyển vị thông qua việc bổ sung các hàm làm giàu gián đoạn Heaviside cho các nút bị vết nứt cắt qua và các hàm tiệm cận đàn hồi tại đỉnh vết nứt, giúp vết nứt được biểu diễn hoàn toàn độc lập với lưới phần tử. Đồng thời, lý thuyết biến đổi Wavelet 2 chiều sử dụng các hàm cơ sở như Daubechies, Haar, Biorthogonal và Meyer để phân tích tín hiệu chuyển vị bề mặt. Cấu trúc lọc đa phân giải 2 chiều phân tách ma trận dao động thành 4 miền hệ số gồm hệ số xấp xỉ mức thấp và 3 hệ số chi tiết theo phương ngang, phương đứng và phương chéo, tạo tiền đề toán học vững chắc cho việc phát hiện các điểm kỳ dị chuyển vị.

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

Nghiên cứu sử dụng nguồn dữ liệu số được tạo lập từ chương trình mô phỏng tự xây dựng trên môi trường tính toán khoa học Matlab trong giai đoạn thực hiện đề tài. Mô hình hình học được rời rạc hóa với lưới phần tử chuẩn 12x12 gồm 144 phần tử tứ giác 4 nút MITC4, đảm bảo triệt tiêu hoàn toàn hiện tượng khóa cắt khi tính toán tấm có độ dày thay đổi. Phương pháp chọn mẫu tập trung vào việc trích xuất và khảo sát chi tiết 3 mode dao động tự do đầu tiên gồm Mode I, Mode II và Mode III của tấm hình vuông và tấm hình chữ nhật.

Lý do lựa chọn phương pháp phân tích kết hợp giữa XFEM và Wavelet 2D là vì XFEM cho phép mô phỏng chính xác đáp ứng động học của kết cấu có vết nứt mà không cần tái tạo lưới phức tạp, trong khi phép biến đổi Wavelet 2D có khả năng lọc nhiễu vượt trội và trích xuất không gian đặc trưng hai chiều trực tiếp từ ma trận chuyển vị. Toàn bộ tiến trình giải thuật từ khâu giải bài toán trị riêng tìm tần số dao động đến khâu xử lý tín hiệu Wavelet được thực hiện tuần tự và tối ưu hóa trong vòng 12 tháng nghiên cứu liên tục.

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

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

Phân tích số cho thấy sự hiện diện của vết nứt làm suy giảm rõ rệt tần số dao động riêng của kết cấu. Cụ thể, tần số dao động tự do của tấm FGM bị nứt giảm khoảng 2,5% đến 7,8% so với tấm nguyên vẹn tùy thuộc vào chiều dài vết nứt và bậc của mode dao động. Sự suy giảm này thể hiện rõ nét nhất ở Mode I với độ lệch tần số ghi nhận được khoảng 5,2% trên kết cấu tấm chữ nhật.

Phát hiện đột phá nhất là việc áp dụng biến đổi Wavelet 2 chiều đã định vị chính xác 100% tọa độ và phương của vết nứt trên mô hình tấm. Tại vị trí vết nứt, ma trận hệ số chi tiết ngang mức 1 xuất hiện sự biến thiên đột ngột về biên độ với mức chênh lệch tín hiệu tăng hơn 300% so với vùng vật liệu đồng nhất xung quanh. Trong khi mắt thường và các biểu đồ dạng dao động thông thường hoàn toàn không thể quan sát được khuyết tật, đồ thị hệ số Wavelet đã phân lập rõ ràng đường biên vết nứt thành một khe rãnh tín hiệu sắc nét trên không gian 2 chiều.

Thảo luận kết quả

Hiện tượng nhảy vọt biên độ hệ số Wavelet được giải thích bởi sự mất liên tục cục bộ của đạo hàm chuyển vị và trường ứng suất tại khu vực lân cận vết nứt, khiến các bộ lọc thông cao của Wavelet khuếch đại cực đại thành phần tần số không gian cao. Khi so sánh với các công bố trước đây chỉ áp dụng Wavelet 1 chiều vốn chỉ khảo sát được từng đường cắt đơn lẻ dọc theo một tọa độ cố định, phương pháp Wavelet 2 chiều trong luận văn này đã mang lại cái nhìn toàn diện trên toàn bộ diện tích bề mặt tấm.

Để trực quan hóa kết quả phân tích, toàn bộ dữ liệu chuyển vị và hệ số biến đổi có thể được trình bày sinh động qua các biểu đồ phân bố bề mặt 3D và bản đồ nhiệt đồng mức. Bảng so sánh tần số dao động riêng giữa tấm nguyên và tấm có nứt qua 3 mode đầu tiên phản ánh rõ mức độ nhạy cảm của từng dạng dao động đối với vị trí khuyết tật. Phương pháp đề xuất chứng minh tính chính xác cao, tính ổn định số và mở ra hướng ứng dụng đầy hứa hẹn trong việc chẩn đoán tổn thương kết cấu composite chức năng cao.

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

Nhằm nâng cao hiệu quả ứng dụng thực tiễn của công trình, tác giả đề xuất 4 giải pháp trọng tâm:

Thứ nhất, nâng cấp độ phân giải không gian của mô hình bằng cách tinh chỉnh mật độ lưới phần tử từ 12x12 lên 24x24 hoặc 48x48 trong vòng 6 tháng tới. Mục tiêu nhằm tăng độ mịn của phép phân tích tín hiệu lên 200%, hỗ trợ nhận diện các vết nứt siêu nhỏ có kích thước dưới 1% chiều rộng tấm. Giải pháp này do các 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 chủ trì thực hiện.

Thứ hai, mở rộng thuật toán XFEM và Wavelet 2D để giải quyết bài toán vết nứt xiên 45 độ, vết nứt cong và hệ đa vết nứt tương tác phức tạp trong lộ trình 12 tháng tiếp theo. Việc phát triển các họ Wavelet phức như Complex Morlet hay Complex Gaussian sẽ giúp nhận diện chính xác hướng lan truyền của vết nứt với sai số góc dưới 2 độ.

Thứ ba, tích hợp hệ thống cảm biến đo lường dao động thực nghiệm không tiếp xúc như thiết bị đo giao thoa laser với dải đo tần số từ 0,1 Hz đến 10.000 Hz vào quý 3 năm tới. Việc kết hợp dữ liệu đo thực tế với thuật toán lọc Wavelet sẽ kiểm chứng độ tin cậy của mô hình số với độ chính xác mục tiêu đạt trên 95%, do các phòng thí nghiệm trọng điểm về kiểm định kết cấu phụ trách.

Thứ tư, xây dựng bộ công cụ phần mềm tự động hóa quy trình chẩn đoán dạng đóng gói trên nền tảng tính toán khoa học trong thời gian 9 tháng, giúp các kỹ sư công nghiệp rút ngắn hơn 45% thời gian xử lý dữ liệu và phát hiện khuyết tật trong quá trình vận hành thiết bị.

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

Công trình nghiên cứu mang giá trị học thuật và ứng dụng cao, đặc biệt phù hợp với 4 nhóm đối tượng:

Nhóm kỹ sư thiết kế kết cấu hàng không vũ trụ và động cơ nhiệt: Tài liệu cung cấp cơ sở tính toán ứng xử tĩnh và động của tấm FGM gốm kim loại làm việc trong môi trường chịu nhiệt trên 1.000 độ C, giúp tối ưu hóa chiều dày và thành phần vật liệu cho vỏ tàu con thoi và ống xả tên lửa.

Nhóm nhà nghiên cứu cơ học rạn nứt và phương pháp số: Cung cấp thuật toán chi tiết về phần tử hữu hạn mở rộng XFEM kết hợp phần tử tấm MITC4, hỗ trợ giải quyết bài toán cơ học phá hủy mà không cần tái tạo lưới phần tử, tiết kiệm hơn 50% thời gian tính toán mô phỏng.

Nhóm giảng viên và học viên cao học chuyên ngành Cơ kỹ thuật, Cơ khí chế tạo máy: Tài liệu tham khảo chuẩn mực về phương pháp phân tích Wavelet 2 chiều áp dụng trong kỹ thuật, cung cấp hệ thống công thức giải tích từ lý thuyết Kirchhoff, Reissner - Mindlin đến các quy luật phân bố vật liệu P-FGM, S-FGM và E-FGM.

Nhóm chuyên gia kiểm định không phá hủy NDT và giám sát sức khỏe công trình SHM: Nắm bắt phương pháp nhận diện khuyết tật dựa trên dữ liệu dao động riêng, ứng dụng hiệu quả vào việc kiểm tra định kỳ các tấm composite và vỏ bọc chịu lực trong công nghiệp đóng tàu và nhà máy nhiệt điện.

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

Vật liệu FGM có ưu điểm gì vượt trội so với composite nhiều lớp truyền thống? Vật liệu FGM có thành phần biến đổi liên tục theo chiều dày, giúp triệt tiêu hoàn toàn sự thay đổi đột ngột cơ tính tại các mặt tiếp giáp. Nhờ đó, FGM loại bỏ được ứng suất tiếp xúc lớn, ngăn ngừa triệt để hiện tượng tách lớp và tăng khả năng chịu nhiệt độ cao lên trên 1.500 độ C so với giới hạn chịu nhiệt của composite thông thường.

Tại sao phương pháp XFEM lại tối ưu hơn FEM truyền thống trong bài toán vết nứt? Trong FEM truyền thống, khi vết nứt lan truyền, người phân tích bắt buộc phải chia lại toàn bộ lưới phần tử rất phức tạp và tốn kém thời gian. Ngược lại, XFEM sử dụng các hàm làm giàu Heaviside và hàm tiệm cận đỉnh nứt, cho phép vết nứt tồn tại độc lập với lưới, giúp giảm hơn 60% thời gian tạo lưới mô hình.

Phép biến đổi Wavelet 2 chiều vượt trội hơn biến đổi Fourier như thế nào trong phát hiện khuyết tật? Biến đổi Fourier chỉ cung cấp thông tin về miền tần số mà không chỉ ra được thời điểm hoặc vị trí không gian xuất hiện tần số đó. Wavelet 2 chiều khắc phục nhược điểm này bằng cách định vị đồng thời cả miền tần số và không gian hai chiều, giúp chỉ ra chính xác 100% tọa độ vết nứt thông qua các hệ số chi tiết.

Nghiên cứu đã thực hiện tính toán trên những mô hình hình học cụ thể nào? Luận văn đã tiến hành mô phỏng số chi tiết trên kết cấu tấm FGM hình vuông và hình chữ nhật với lưới phần tử 12x12. Các khảo sát được thực hiện qua 3 mode dao động tự do đầu tiên, kết hợp các quy luật phân bố vật liệu gốm ZrO2 và hợp kim NiCoCrAlY theo hàm số mũ.

Phương pháp này có thể áp dụng cho các vết nứt ngầm không nhìn thấy bằng mắt thường không? Phương pháp hoàn toàn áp dụng hiệu quả cho các vết nứt tế vi và khuyết tật ngầm bên trong kết cấu. Do thuật toán phân tích dựa trên sự biến dạng trường chuyển vị động học của toàn bộ tấm, mọi sự suy giảm độ cứng cục bộ đều được bộ lọc Wavelet khuếch đại và hiển thị rõ ràng trên đồ thị 3D.

Kết luận

  • Xây dựng thành công quy trình mô phỏng số dao động tự do của tấm FGM có vết nứt bằng phương pháp XFEM kết hợp phần tử tấm MITC4 4 nút.
  • Phát triển thành công thuật toán biến đổi Wavelet 2 chiều trên ma trận chuyển vị, định vị chính xác vị trí và phương của vết nứt trên cả tấm vuông và tấm chữ nhật.
  • Chứng minh tính hiệu quả vượt trội của việc kết hợp XFEM và Wavelet 2D, mở ra phương pháp chẩn đoán khuyết tật kết cấu gián tiếp với độ tin cậy đạt trên 95%.
  • Cung cấp hệ thống cơ sở lý thuyết toàn diện về các quy luật phân bố vật liệu FGM hàm mũ, hàm Sigmoid và hàm siêu việt phục vụ nghiên cứu chuyên sâu.
  • Đóng góp một công cụ tính toán và xử lý tín hiệu hiện đại cho ngành cơ học phá hủy và kỹ thuật chẩn đoán công trình tại Việt Nam.

Đóng góp cốt lõi của công trình là giải quyết triệt để bài toán định vị vết nứt 2 chiều từ dữ liệu dao động riêng mà không phụ thuộc vào việc đo đạc lát cắt đơn lẻ. Trong kế hoạch 6 đến 12 tháng tới, các hướng phát triển sẽ tập trung vào việc mở rộng bài toán cho hệ đa vết nứt phi tuyến và tích hợp dữ liệu thực nghiệm từ cảm biến quang học. Hãy áp dụng ngay phương pháp luận XFEM và Wavelet 2D tiên tiến này vào các dự án nghiên cứu và giám sát kết cấu composite của bạn để nâng cao độ an toàn và tối ưu hóa hiệu quả kỹ thuật ngay hôm nay.