Tổng quan nghiên cứu

Trong kỹ thuật công trình, khuyết tật nứt ngầm là nguyên nhân hàng đầu gây ra hơn 40% các sự cố suy giảm khả năng chịu tải và phá hủy đột ngột của kết cấu tấm. Việc nhận diện và đánh giá vết nứt ngay từ giai đoạn khởi phát, khi kích thước vết nứt mới chỉ chiếm từ 2% đến 5% kích thước hình học cấu kiện, đóng vai trò sống còn trong công tác duy tu và bảo trì công trình dân dụng và công nghiệp. Tuy nhiên, phương pháp phần tử hữu hạn truyền thống gặp rào cản kỹ thuật rất lớn khi mô tả vết nứt di động do bắt buộc phải tái tạo lưới liên tục, làm tăng thời gian tính toán lên gấp 3 đến 5 lần.

Luận văn thạc sĩ tập trung giải quyết bài toán cốt lõi: phân tích độ nhạy của chuyển vị tấm chịu uốn có vết nứt bằng cách kết hợp phương pháp phần tử hữu hạn mở rộng (XFEM) với công cụ biến đổi Wavelet liên tục. Mục tiêu cụ thể là xây dựng thuật toán số hóa hoàn chỉnh trên nền tảng MATLAB nhằm xác định trường chuyển vị tĩnh, tần số dao động riêng và đáp ứng động lực học của tấm, từ đó khai thác hệ số Wavelet để định vị chính xác vị trí cũng như kích thước vết nứt.

Nghiên cứu được triển khai tại Trường Đại học Bách Khoa – ĐHQG TP.HCM trong khoảng thời gian từ tháng 01/2014 đến tháng 11/2014. Phạm vi nghiên cứu khảo sát tấm phẳng Mindlin đẳng hướng dưới hai điều kiện biên chuẩn tắc là tựa đơn bốn cạnh và ngàm bốn cạnh, chịu tải trọng tĩnh hai điểm tập trung cùng các tải trọng động điều hòa với tần số góc từ 2π rad/s đến 20π rad/s và tải trọng xung hình sin. 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 chẩn đoán không phá hủy với độ tin cậy cao, giảm hơn 60% khối lượng chia lưới so với phương pháp số cổ điể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 xây dựng trên nền tảng ba trụ cột lý thuyết cơ học và toán học tính toán hiện đại:

Thứ nhất, lý thuyết tấm Mindlin-Reissner được áp dụng để mô tả trạng thái uốn của tấm phẳng đàn hồi đẳng hướng, có kể đến ảnh hưởng của biến dạng trượt cắt ngang với hệ số hiệu chỉnh k = 5/6 (tương đương khoảng 0.833). Trường chuyển vị tại mỗi nút bao gồm ba bậc tự do độc lập là độ võng uốn và hai góc xoay tiết diện. Để triệt tiêu hoàn toàn hiện tượng khóa cắt khi tỷ lệ bề dày trên nhịp giảm xuống dưới 0.01, mô hình phần tử hỗn hợp MITC4 với bốn nút đẳng tham số được tích hợp chặt chẽ.

Thứ hai, phương pháp phần tử hữu hạn mở rộng (XFEM) mở rộng không gian xấp xỉ tiêu chuẩn bằng cách bổ sung hàm bước nhảy Heaviside cho các nút bị vết nứt cắt ngang và bốn hàm tiệm cận cực địa phương cho các nút bao quanh đỉnh vết nứt. Kỹ thuật hàm đường mức (Level Set Method) được dùng để biểu diễn hình học vết nứt độc lập hoàn toàn với lưới phần tử.

Thứ ba, lý thuyết biến đổi sóng Wavelet liên tục (CWT) sử dụng họ sóng Symlets (điển hình là sym4 và sym8) làm công cụ phân tích đa độ phân giải. Biến đổi Wavelet cho phép chuyển đổi tín hiệu phân bố độ võng không gian thành các hệ số chi tiết, khuếch đại các điểm kỳ dị gián đoạn bậc cao tại vị trí hư hại.

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

Nguồn dữ liệu của luận văn được thu thập từ mô hình mô phỏng số thông qua chương trình tính toán độc lập viết trên ngôn ngữ MATLAB, sau đó kiểm chứng đối chiếu với phần mềm thương mại ANSYS 14.5 và các nghiệm giải tích lý thuyết.

Cỡ mẫu nghiên cứu bao gồm 8 bài toán mô phỏng kết cấu chuyên sâu với hệ thống lưới phần tử từ 400 đến 1600 phần tử tứ giác Q4. Phương pháp chọn mẫu theo quy trình quét tham số tất định có kiểm soát, khảo sát toàn diện sự biến thiên của ba nhóm biến số then chốt: chiều dài vết nứt không thứ nguyên thay đổi từ 0.1 đến 0.6, tọa độ vị trí vết nứt thay đổi từ 0.1 đến 0.9 theo chiều rộng tấm, và tỷ lệ chiều dày tấm thay đổi từ 0.01 đến 0.2.

Lý do lựa chọn phương pháp phân tích này là vì XFEM loại bỏ triệt để nhu cầu làm mịn và tái cấu trúc lưới quanh đỉnh nứt, trong khi phương pháp tích phân từng bước Newmark dạng chuyển vị (với hai tham số gia tốc trung bình alpha = 0.25 và delta = 0.5) bảo đảm tính ổn định vô điều kiện khi giải bài toán động. Toàn bộ quy trình mô phỏng, phân tích dữ liệu và đánh giá độ nhạy được thực hiện liên tục trong tiến trình 10 tháng.

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

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

Quá trình mô phỏng số và phân tích tín hiệu đã mang lại các phát hiện cốt lõi sau:

Thứ nhất, thuật toán XFEM tự phát triển đạt độ chính xác vượt trội khi so sánh với ANSYS 14.5 và lời giải giải tích. Độ lệch tương đối về chuyển vị uốn cực đại chỉ nằm trong khoảng 0.15% đến 1.25% đối với tấm tựa đơn và dưới 1.10% đối với tấm ngàm bốn cạnh. Tần số dao động riêng bậc một đến bậc bốn có sai số so với lý thuyết không vượt quá 0.95%.

Thứ hai, chiều dài vết nứt ảnh hưởng trực tiếp đến độ suy giảm độ cứng của tấm. Khi tỷ lệ chiều dài vết nứt tăng từ 0.1 lên 0.5, độ võng tĩnh cực đại tại tâm tấm tăng từ 14.8% đến 28.6% tùy thuộc vào điều kiện biên. Đồng thời, biên độ hệ số chi tiết Wavelet tại vị trí vết nứt tăng vọt hơn 3.5 lần so với trạng thái tấm lành lặn, tạo thành đỉnh sóng nhận dạng vô cùng rõ nét.

Thứ ba, trong bài toán động học với tải trọng điều hòa ở tần số góc 2π rad/s và 10π rad/s, biến đổi Wavelet phát hiện chính xác tọa độ vết nứt với sai số không gian dưới 0.8% kích thước hình học tấm. Ngay cả dưới tác dụng của tải trọng xung ngắn hạn, tín hiệu Wavelet vẫn giữ được độ nhạy phân giải vị trí vết nứt rất cao.

Thứ tư, điều kiện biên ngàm giúp giảm độ võng động tổng thể từ 42% đến 58% so với điều kiện biên tựa đơn, nhưng sự gián đoạn của đạo hàm trường chuyển vị qua vết nứt vẫn được toán tử Wavelet tách lọc và khuếch đại thành công mà không bị triệt tiêu bởi độ cứng biên.

Thảo luận kết quả

Sự xuất hiện của vết nứt gây ra gián đoạn trường ứng suất và làm suy giảm mô men quán tính cục bộ của tiết diện. Mặc dù độ võng tổng thể thay đổi rất nhỏ (chỉ vài phần trăm milimét trong các kịch bản nứt sớm), nhưng đạo hàm bậc một và bậc hai của độ võng theo tọa độ không gian lại bị đứt đoạn mạnh mẽ tại đúng đường nứt. Biến đổi Wavelet với đặc tính định xứ không gian - tần số đã đóng vai trò như một bộ lọc thông cao tinh vi, biến các gián đoạn vi mô này thành các đỉnh cực trị sắc nhọn trên đồ thị hệ số chi tiết.

Kết quả này hoàn toàn tương đồng và mở rộng các kết luận trước đây của các nhà khoa học trong ngành về độ nhạy độ võng dầm, đồng thời khắc phục triệt để hiện tượng khóa cắt mà các nghiên cứu dùng phần tử cổ điển thường gặp phải.

Về mặt biểu diễn dữ liệu, toàn bộ diễn biến động lực học của kết cấu được mô tả sinh động qua các đồ thị quan hệ chuyển vị - thời gian với bước thời gian delta t = 0.001 giây, các biểu đồ phân bố hệ số Wavelet hai chiều thể hiện rõ đỉnh xung định vị, cùng hệ thống bảng số liệu đối chuẩn đa cột so sánh trực tiếp kết quả giữa MATLAB, ANSYS 14.5 và giải tích.

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

Dựa trên các kết quả đạt được, luận văn đưa ra bốn nhóm giải pháp mang tính ứng dụng cao:

  1. Chuẩn hóa quy trình chẩn đoán khuyết tật kết cấu dạng tấm: Các cơ quan kiểm định xây dựng và viện chuyên ngành cần ứng dụng tích hợp thuật toán XFEM và biến đổi Wavelet vào quy trình kiểm tra công trình định kỳ. Mục tiêu đặt ra là phát hiện sớm các vết nứt bề mặt có chiều sâu dưới 1.0 mm, giúp rút ngắn 40% thời gian quan trắc so với các phương pháp thử tải phá hủy, triển khai áp dụng từ quý 1 đến quý 4 năm tới.

  2. Tự động hóa module xử lý tín hiệu quan trắc: Nhóm nghiên cứu kết cấu và các kỹ sư lập trình cần đóng gói bộ công cụ MATLAB/Python thành phần mềm thương mại có giao diện trực quan. Hệ thống cần hướng đến mục tiêu tiếp nhận dữ liệu cảm biến đo độ võng hiện trường và trích xuất vị trí hư hỏng theo thời gian thực với độ trễ xử lý dưới 30 giây cho diện tích sàn 50 m2, hoàn thiện trong vòng 12 tháng.

  3. Mở rộng mô hình tính toán cho vật liệu tiên tiến: Các phòng thí nghiệm trọng điểm và trường đại học kỹ thuật nên tiếp tục phát triển giải thuật cho tấm vật liệu composite biến tính FGM và kết cấu tấm có đa vết nứt tương tác phức tạp. Mục tiêu nâng cao độ chính xác dự báo tuổi thọ mỏi lên trên 95% trong lộ trình nghiên cứu 24 tháng.

  4. Tổ chức đào tạo chuyển giao công nghệ: Hội đồng khoa học ngành xây dựng cần chủ trì tổ chức từ 3 đến 5 khóa tập huấn chuyên sâu hàng năm về cơ học phá hủy tính toán và xử lý tín hiệu số, đào tạo hơn 500 kỹ sư kết cấu thành thạo công nghệ chẩn đoán hiện đại từ năm 2026.

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

Nội dung và kết quả của luận văn mang lại giá trị thiết thực cho bốn nhóm đối tượng chính:

Thứ nhất, học viên cao học và nghiên cứu sinh ngành Xây dựng công trình, Cơ học kỹ thuật và Kỹ thuật công trình biển. Luận văn cung cấp toàn bộ cơ sở toán học của phần tử MITC4, kỹ thuật làm giàu nút trong XFEM và thuật toán Wavelet rời rạc, làm tài liệu tham khảo mẫu mực cho các đề tài mô phỏng rạn nứt kết cấu.

Thứ hai, kỹ sư kết cấu tại các công ty tư vấn thiết kế và thẩm tra công trình. Nghiên cứu giúp kỹ sư nắm rõ quy luật suy giảm khả năng chịu uốn và tập trung ứng suất quanh vết nứt, phục vụ tính toán gia cường an toàn cho sàn bê tông nhịp lớn, bản mặt cầu hoặc vỏ thép công trình chịu trên 100.000 chu kỳ tải trọng động.

Thứ ba, chuyên gia kiểm định và quan trắc sức khỏe công trình (SHM). Tài liệu là cẩm nang ứng dụng biến đổi sóng Wavelet để xử lý tín hiệu đo đạc từ hiện trường, giúp định vị chính xác vị trí khuyết tật ngầm trên các nhịp dầm cầu dài từ 15 m đến 30 m mà không cần phá vỡ cấu trúc.

Thứ tư, giảng viên đại học chuyên trách các bộ môn Cơ học kết cấu, Phương pháp tính và Phương pháp phần tử hữu hạn nâng cao, sử dụng làm học liệu bài giảng và định hướng nghiên cứu chuyên sâu cho sinh viên.

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

Phương pháp XFEM có ưu thế vượt trội gì so với phương pháp phần tử hữu hạn truyền thống? Phương pháp XFEM giải quyết nhược điểm phải chia lại lưới của FEM truyền thống khi vết nứt xuất hiện. Nhờ các hàm làm giàu Heaviside và hàm tiệm cận đỉnh nứt, XFEM mô phỏng chính xác sự gián đoạn mà không cần lưới phần tử trùng với đường nứt, giúp giảm hơn 60% khối lượng chia lưới và hạn chế sai số biến dạng phần tử.

Biến đổi Wavelet định vị vết nứt trong kết cấu dựa trên cơ chế nào? Biến đổi Wavelet hoạt động như một kính hiển vi toán học phân tích tín hiệu trên cả miền không gian và tần số. Khi trường chuyển vị của tấm đi qua vị trí có vết nứt, dù độ võng chỉ biến thiên rất nhỏ khoảng 0.02 mm, hệ số chi tiết Wavelet họ Symlets vẫn tạo ra các bước nhảy cực trị rõ rệt, xác định chính xác tọa độ nứt.

Tại sao luận văn lại kết hợp lý thuyết Mindlin với phần tử MITC4? Lý thuyết Mindlin-Reissner xét đến biến dạng trượt cắt với hệ số k = 5/6, áp dụng tốt cho cả tấm mỏng lẫn tấm dày với tỷ lệ h/a từ 0.01 đến 0.2. Việc ứng dụng phần tử MITC4 giúp loại bỏ triệt để hiện tượng khóa cắt, bảo đảm nghiệm số hội tụ ổn định với sai số dưới 1.0% khi chiều dày tấm giảm dần.

Thuật toán Newmark giải bài toán dao động của tấm nứt như thế nào? Luận văn sử dụng phương pháp tích phân từng bước Newmark dạng chuyển vị với bước thời gian delta t = 0.001 giây và tham số gia tốc trung bình alpha = 0.25, delta = 0.5. Thuật toán này bảo đảm tính ổn định vô điều kiện, tái hiện chính xác đáp ứng động lực học của tấm dưới các tải trọng điều hòa có tần số từ 2π đến 20π rad/s.

Độ tin cậy của mô hình tính toán trong luận văn được kiểm chứng qua đâu? Các kết quả số lập trình trên MATLAB được đối chiếu toàn diện với phần mềm mô phỏng ANSYS 14.5 và các lời giải giải tích tiêu chuẩn. Mức độ sai lệch về tần số dao động riêng dưới 1.15% và chuyển vị uốn dưới 1.25%, khẳng định tính chính xác tuyệt đối và khả năng ứng dụng thực tế cao của đề tài.

Kết luận

  • Xây dựng thành công thuật toán số tích hợp XFEM và phần tử tấm MITC4, loại bỏ hiện tượng khóa cắt trong bài toán uốn tấm Mindlin có vết nứt.
  • Tính toán chính xác trường chuyển vị tĩnh, tần số dao động riêng và đáp ứng động lực học dưới tải trọng điều hòa (2π đến 20π rad/s) cùng tải trọng xung.
  • Ứng dụng xuất sắc công cụ biến đổi Wavelet họ Symlets để trích xuất độ nhạy chuyển vị, định vị chính xác vị trí vết nứt với sai số không gian dưới 0.8%.
  • Kiểm chứng chặt chẽ mô hình số với phần mềm ANSYS 14.5 và nghiệm giải tích, đạt độ sai lệch kiểm định tin cậy dưới 1.25%.
  • Đề xuất khung giải pháp toàn diện và lộ trình 12 đến 24 tháng nhằm ứng dụng xử lý tín hiệu Wavelet vào hệ thống quan trắc sức khỏe công trình thực tế.

Luận văn thạc sĩ là một công trình nghiên cứu công phu, đặt nền móng vững chắc cho việc ứng dụng các phương pháp số hiện đại và xử lý tín hiệu nâng cao trong chẩn đoán kỹ thuật công trình. Quý kỹ sư, nhà khoa học và các đơn vị kiểm định hãy chủ động khai thác các giải thuật tính toán này để nâng cao độ an toàn và chất lượng cho các công trình xây dựng ngay hôm nay.