Tổng quan nghiên cứu

Trong kỹ thuật xây dựng và kết cấu công trình hiện đại, các cấu kiện tấm mỏng chiếm hơn 60% tỷ trọng kết cấu trong nhiều hạng mục trọng điểm như sàn thép, bản mặt cầu, vỏ tàu thủy và kết cấu công trình công nghiệp. Tuy nhiên, sự suy giảm cơ học do nứt mỏi bề mặt tiềm ẩn nguy cơ phá hủy nghiêm trọng, điển hình như sự cố sụp đổ cầu vượt de la Concorde tại Quebec, Canada năm 2006 gây tổn thất lớn về người và tài sản do không phát hiện kịp thời các khuyết tật tích tụ. Vấn đề cốt lõi đặt ra là các phương pháp quan trắc sức khỏe kết cấu (SHM - Structural Health Monitoring) truyền thống dựa vào sự suy giảm tần số dao động chỉ nhạy với hư hỏng tổng thể hoặc các vết nứt xuyên tấm lớn, hoàn toàn hạn chế trong việc định vị chính xác các vết nứt bề mặt nhỏ nông.

Mục tiêu trọng tâm của luận văn thạc sĩ do học viên Lê Thanh Cao thực hiện dưới sự hướng dẫn khoa học của Tiến sĩ Hồ Đức Duy tại Trường Đại học Bách Khoa - Đại học Quốc gia TP.HCM là phát triển và mở rộng phương pháp năng lượng biến dạng (Strain Energy Method - SEM) hai chiều (2-D) nhằm định vị chính xác vị trí và ước lượng kích thước vết nứt bề mặt trên kết cấu tấm mỏng. Phạm vi nghiên cứu tập trung vào các kết cấu tấm đẳng hướng hình vuông và hình chữ nhật kích thước tiêu chuẩn 240 x 240 x 2 mm dưới ba điều kiện biên cơ bản: tự do hoàn toàn, liên kết khớp và liên kết ngàm, được triển khai thực hiện từ ngày 10/02/2014 đến ngày 14/11/2014. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp thuật toán kiểm định không phá hủy với chi phí tính toán tối ưu, có khả năng nhận dạng các vùng suy giảm độ cứng cục bộ từ 5% đến 10% và giảm thiểu tới hơn 40% chi phí đo đạc so với phương pháp siêu âm rà quét bề mặt thủ công.

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 dựa trên nền tảng Lý thuyết tấm mỏng Kirchhoff với giả thiết mặt trung bình không chịu biến dạng trượt ngang và Lý thuyết dao động động lực học công trình cho hệ nhiều bậc tự do không cản ($M\ddot{v} + Kv = 0$). Năng lượng biến dạng uốn đàn hồi của tấm Kirchhoff hai chiều được mô tả thông qua hàm tích phân vi phân năng lượng biến dạng danh nghĩa của D. Young (1956).

Khung lý thuyết của luận văn tích hợp 4 khái niệm và mô hình then chốt:

  1. Phương pháp năng lượng biến dạng hai chiều (SEM 2-D): Mở rộng từ mô hình gốc của Stubbs et al. (1995) và Cornwell et al. (1997), tính toán tỷ số năng lượng biến dạng phân đoạn giữa trạng thái nguyên vẹn ($F_{k,ij}$) và trạng thái hư hỏng ($F^*{k,ij}$) để xác định chỉ số hư hỏng $\beta{ij}$ cùng chỉ số chuẩn hóa $Z_{ij}$.
  2. Phương pháp sai phân trung tâm bậc hai (CDM - Central Differential Method): Ứng dụng khai triển chuỗi Taylor bậc 4 nhằm xấp xỉ hóa đạo hàm bậc hai của hàm dạng dao động theo các biến tọa độ $x, y$ với sai số rút gọn bậc $O(h^2)$.
  3. Tiêu chuẩn chứng thực dao động (MAC - Modal Assurance Criterion): Thước đo đối sánh véc-tơ dạng dao động trước và sau hư hỏng trong phạm vi từ 0 đến 1 nhằm kiểm soát mức độ tương quan dao động.
  4. Ngưỡng hư hỏng ($Z_0$) và Khả năng chẩn đoán ($KNCD$): Bộ công cụ định lượng đánh giá tỷ lệ nhận dạng chiều dài vết nứt theo ba mức phân loại: Kém ($KNCD < 50%$), Trung bình ($50% \le KNCD \le 80%$) và Tốt ($KNCD > 80%$).

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

Nguồn dữ liệu phân tích được trích xuất hoàn toàn từ mô hình mô phỏng số phần tử hữu hạn (FEM) trên phần mềm ANSYS Mechanical APDL phiên bản 14.0, sau đó được lập trình xử lý thuật toán tự động trên môi trường MATLAB kết hợp Microsoft Excel.

Về cỡ mẫu và kỹ thuật chia lưới: Tấm kết cấu kích thước 240 x 240 x 2 mm được chia lưới mịn với kích thước phần tử 1 mm, tạo ra cỡ mẫu gồm 58.081 nút trên mặt phẳng trung bình để trích xuất trường chuyển vị dạng dao động.

Về phương pháp chọn phần tử mô hình: Nghiên cứu áp dụng quy trình mô hình hai giai đoạn: sử dụng phần tử vỏ 4 nút SHELL181 (6 bậc tự do mỗi nút) để chia lưới bề mặt kiểm soát hình học, sau đó đùn tạo khối thành phần tử thể tích 8 nút SOLID185 (3 bậc tự do tịnh tiến mỗi nút) nhằm mô phỏng vết nứt không gian 3D thực tế với kích thước chiều dài 40 mm, chiều rộng 0.8 mm và chiều sâu 1 mm (chiếm đúng 50% chiều dày tấm).

Lý do lựa chọn phương pháp phân tích: Việc ứng dụng phương pháp sai phân trung tâm kết hợp với kỹ thuật ngoại suy đa thức bậc ba spline hai chiều (thông qua hàm interp2) tại các điểm biên ảo cho phép khắc phục hiện tượng mất dữ liệu vi phân ở biên tấm. Phương pháp này không yêu cầu chuẩn hóa véc-tơ dạng dao động theo khối lượng, tạo điều kiện thuận lợi khi ứng dụng dữ liệu đo đạc thực tế từ kích thích môi trường. Toàn bộ quy trình khảo sát 6 bài toán số được thực hiện chặt chẽ trong khoảng thời gian nghiên cứu 9 tháng năm 2014.

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

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

  1. Khả năng định vị vết nứt chính xác vượt trội ở ngưỡng tối ưu: Trên lưới phân tích 48 x 48 phần tử (gồm 2.304 phần tử con), với ngưỡng hư hỏng chọn lọc $Z_0 = 50% \beta_{max}$, phương pháp SEM phát hiện chính xác 100% vị trí vết nứt dài 40 mm nằm tại trung tâm tấm ở cả ba điều kiện biên: tự do hoàn toàn, liên kết khớp và liên kết ngàm. Khi hạ ngưỡng xuống mức $30% \beta_{max}$, tỷ lệ chẩn đoán sai lệch tăng lên khoảng 18% do xuất hiện các cảnh báo giả tại vùng góc biên.
  2. Mức độ nhạy cảm của chỉ số năng lượng so với tần số dao động: Kết quả tính toán tiêu chuẩn MAC cho 6 mode dao động đầu tiên luôn duy trì ở mức rất cao, từ 0.9985 đến 0.9999. Trong khi tần số dao động riêng chỉ suy giảm rất nhỏ từ 0.05% đến 0.42% (không đủ để kết luận hư hỏng bằng phương pháp đo tần số), thì chỉ số hư hỏng chuẩn hóa $Z_{ij}$ tại vị trí có vết nứt lại tăng vọt từ 2.5 đến 4.8 lần so với giá trị trung bình của tấm nguyên vẹn.
  3. Hiệu quả nhận dạng các kịch bản đa vết nứt phức tạp: Trong bài toán tấm vuông biên khớp có hai vết nứt vuông góc và bài toán tấm chữ nhật biên ngàm có hai vết nứt song song, thuật toán SEM vẫn định dạng chính xác đồng thời cả hai vết nứt với chỉ số khả năng chẩn đoán $KNCD$ đạt từ 85% đến 92.5% tại ngưỡng $50% \beta_{max}$.
  4. Tối ưu hóa độ chính xác bằng thuật toán SEM cục bộ hai bước: Khi áp dụng phương pháp SEM cục bộ trên lưới siêu mịn 240 x 240 phần tử (tương đương 57.600 phần tử con), chỉ số $KNCD$ đạt mức hoàn hảo 100%, định vị vết nứt với sai số hình học nhỏ hơn 0.5 mm.

Thảo luận kết quả

Nguyên nhân cốt lõi giúp SEM 2-D đạt độ chính xác cao là do thuật toán tập trung vào đạo hàm bậc hai của hàm dạng chuyển vị ($\partial^2 \phi / \partial x^2$ và $\partial^2 \phi / \partial y^2$). Đạo hàm bậc hai đại diện cho độ cong uốn, vốn biến thiên đột ngột tại đúng vị trí suy giảm mômen quán tính tiết diện do nứt bề mặt.

So sánh với nghiên cứu thực nghiệm của Hu và các cộng sự (2008) vốn chỉ áp dụng phương pháp cầu phương vi phân (DQM) cho tấm nhôm biên tự do với một vết nứt đơn lẻ, nghiên cứu này đã chứng minh giải thuật sai phân trung tâm (CDM) mở rộng hoàn toàn thích ứng và duy trì độ tin cậy tuyệt đối trên các điều kiện biên ngàm và khớp phức tạp.

Dữ liệu chẩn đoán được biểu diễn trực quan thông qua đồ thị không gian 3 chiều của chỉ số $Z_{ij}$ và biểu đồ phân bố bề mặt. Trên các đồ thị này, các đỉnh nhọn vượt qua mặt phẳng ngưỡng $Z_0$ tạo thành các vệt thẳng hàng, thể hiện sự tương thích hoàn toàn với tọa độ hình học thực tế của vết nứt trong ma trận index dữ liệu.

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

  1. Tiêu chuẩn hóa quy trình chẩn đoán kết cấu hai cấp độ: Các đơn vị tư vấn và ban quản lý công trình cần áp dụng quy trình quét sơ bộ bằng lưới thưa 24 x 24 trong vòng 15 phút để khoanh vùng hư hỏng, sau đó tự động tái lập lưới mịn 240 x 240 tại vùng cục bộ nhằm xác định chính xác chiều dài vết nứt với độ tin cậy đạt 98% đến 100%, hoàn thiện hướng dẫn kỹ thuật trước quý 4 năm 2025.
  2. Tích hợp giải thuật SEM vào hệ thống quan trắc cảm biến không dây (WSSN): Các trung tâm quản lý cầu đường và công trình biển cần triển khai mạng lưới gia tốc kế đa kênh để đo đạc dạng dao động định kỳ 6 tháng một lần, hướng tới mục tiêu phát hiện sớm các vết nứt bề mặt có độ sâu từ 0.5 mm trở lên trên các kết cấu dầm bản thép.
  3. Thiết lập ngưỡng hư hỏng $Z_0$ linh hoạt theo điều kiện biên: Nhóm kỹ sư phân tích dữ liệu nên cài đặt ngưỡng mặc định $Z_0 = 50% \beta_{max}$ cho các kết cấu tấm liên kết ngàm hoặc khớp, và điều chỉnh lên mức $60% \beta_{max}$ đối với kết cấu có biên tự do nhằm triệt tiêu 100% các điểm nhiễu giả biên trong giai đoạn bảo trì giai đoạn 2025-2026.
  4. Mở rộng thử nghiệm trên vật liệu composite và kết cấu tấm có gân: Các viện nghiên cứu chuyên ngành xây dựng và hàng không cần chủ trì thực hiện các nghiên cứu thực nghiệm mở rộng trên tấm vật liệu phân lớp dị hướng trong giai đoạn 2025-2027 nhằm gia tăng khả năng ứng dụng thực tiễn trong ngành hàng không vũ trụ và đóng tàu cao tốc.

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

  1. Kỹ sư kết cấu và chuyên viên kiểm định công trình: Nắm vững quy trình chẩn đoán hư hỏng không phá hủy để ứng dụng trực tiếp vào công tác bảo trì định kỳ bản mặt cầu thép, sàn liên hợp và vỏ bồn bể công nghiệp, giúp tiết kiệm khoảng 35% chi phí vận hành bảo dưỡng.
  2. Học viên cao học và nghiên cứu sinh ngành Kỹ thuật Xây dựng và Cơ học: Tiếp cận tài liệu tham khảo chi tiết về kỹ thuật mô phỏng phần tử hữu hạn nâng cao trên ANSYS APDL, kỹ thuật xử lý dữ liệu ma trận chuyển vị trên MATLAB và phương pháp tính đạo hàm sai phân trung tâm.
  3. Giảng viên và nhà khoa học chuyên ngành Động lực học công trình: Sử dụng toàn bộ khung lý thuyết mở rộng của SEM 2-D và tiêu chuẩn MAC làm học liệu giảng dạy chuyên đề cao học về giám sát sức khỏe kết cấu (SHM) và đánh giá an toàn công trình.
  4. Doanh nghiệp công nghệ phát triển giải pháp IoT giám sát công trình: Khai thác giải thuật trích xuất đặc trưng dao động để lập trình nhúng vào các thiết bị đo lường thông minh và hệ thống cảm biến WSSN phục vụ cảnh báo sớm sự cố kết cấu trong thời gian thực.

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

  1. Phương pháp năng lượng biến dạng (SEM) có ưu thế gì vượt trội so với phương pháp theo dõi tần số tự nhiên? Phương pháp tần số tự nhiên chỉ phản ánh sự suy giảm độ cứng tổng thể của toàn bộ kết cấu và rất kém nhạy với các khuyết tật nông (tần số riêng thường chỉ suy giảm dưới 0.5%). Ngược lại, SEM khai thác độ cong dạng dao động thông qua đạo hàm bậc hai, giúp định vị chính xác các vết nứt bề mặt cục bộ chỉ sâu 1 mm trên tấm dày 2 mm.

  2. Vì sao luận văn sử dụng phương pháp sai phân trung tâm (CDM) thay vì cầu phương vi phân (DQM)? Phương pháp sai phân trung tâm (CDM) có cấu trúc giải thuật đơn giản, tốc độ xử lý ma trận nhanh và dễ tự động hóa trên MATLAB khi phân tích hệ thống lưới điểm lớn. Kết hợp cùng hàm ngoại suy đa thức bậc ba spline tại các nút cận biên ảo, CDM triệt tiêu sai số biên và mang lại độ chính xác xấp xỉ bậc $O(h^2)$ trên toàn bộ miền tấm.

  3. Dữ liệu dạng dao động đo đạc thực tế có bắt buộc phải chuẩn hóa theo khối lượng không? Thuật toán SEM 2-D chỉ sử dụng tỷ số năng lượng biến dạng phân đoạn giữa trạng thái ban đầu và trạng thái hư hỏng, do đó các dạng dao động không cần chuẩn hóa theo khối lượng. Đặc tính này mang lại lợi thế ứng dụng thực tế rất lớn khi kỹ sư chỉ cần sử dụng các kỹ thuật đo kích thích môi trường xung quanh mà không cần đo lực kích thích đầu vào.

  4. Ngưỡng hư hỏng $Z_0$ nào được khuyến nghị áp dụng để đạt kết quả chẩn đoán tối ưu? Dựa trên kết quả khảo sát số của 6 bài toán với đầy đủ các điều kiện biên tự do, khớp và ngàm, ngưỡng $Z_0 = 50% \beta_{max}$ mang lại hiệu quả tối ưu nhất. Tại ngưỡng này, chỉ số nhận dạng chiều dài vết nứt ($KNCD$) đạt mức từ 85% đến 100% mà hoàn toàn không bị ảnh hưởng bởi hiện tượng báo động giả tại các vùng biên.

  5. Việc phối hợp hai loại phần tử SHELL181 và SOLID185 trong mô hình ANSYS nhằm mục đích gì? Phần tử SHELL181 được sử dụng ở bước tạo mặt ban đầu để chia lưới bề mặt phẳng đồng đều và kiểm soát nút chính xác. Sau đó, cấu trúc được đùn thành phần tử khối 3D SOLID185 với 8 nút nhằm mô phỏng đúng hình học không gian của khe nứt (sâu 1 mm, rộng 0.8 mm), giúp phản ánh chân thực trạng thái tập trung ứng suất và biến dạng cục bộ.

Kết luận

  • Luận văn đã mở rộng và hoàn thiện thành công phương pháp năng lượng biến dạng (SEM 2-D) kết hợp giải thuật sai phân trung tâm để chẩn đoán chính xác vết nứt bề mặt trên kết cấu tấm mỏng dưới ba điều kiện biên: tự do, khớp và ngàm.
  • Thiết lập hoàn chỉnh quy trình liên thông tự động từ mô hình hóa 3D trên ANSYS APDL (với hơn 58.000 nút) đến thuật toán phân tích ma trận chỉ số hư hỏng trên MATLAB.
  • Xác lập bộ tiêu chuẩn đánh giá định lượng gồm ngưỡng hư hỏng $Z_0 = 50% \beta_{max}$ và chỉ số $KNCD$, đạt tỷ lệ nhận dạng chiều dài vết nứt từ 85% đến 100% trên các lưới phân tích 48 x 48 và 240 x 240.
  • Đóng góp giải pháp khoa học có tính khả thi cao, mở đường cho việc tích hợp thuật toán vào các hệ sinh thái mạng cảm biến thông minh không dây (WSSN) phục vụ giám sát công trình trong giai đoạn 2025-2030.
  • Các kỹ sư và chuyên gia công trình hãy ứng dụng ngay giải thuật SEM 2-D để nâng cao năng lực quan trắc, tối ưu hóa chi phí bảo trì và đảm bảo an toàn tuyệt đối cho các hệ thống kết cấu tấm hiện đại.