Tổng quan nghiên cứu

Trong kỹ thuật kết cấu hiện đại, sự suy giảm độ bền cơ học do hư hỏng cục bộ chiếm hơn 65% nguyên nhân dẫn đến các sự cố sụp đổ công trình bất ngờ. Để đảm bảo an toàn vận hành cho các cấu kiện dạng thanh như cọc chịu lực (pile) và thanh giàn không gian (truss), việc phát hiện sớm khuyết tật là yêu cầu sống còn. Vấn đề nghiên cứu trọng tâm là xử lý bài toán ngược chẩn đoán khuyết tật kết cấu: xác định chính xác số lượng, vị trí và kích thước vết nứt khi chỉ tiếp cận được các phản ứng dao động tổng thể.

Mục tiêu cụ thể của công trình là thiết lập hệ phương trình tần số dạng tường minh cho dao động dọc trục của thanh đàn hồi có số lượng vết nứt bất kỳ, đồng thời xây dựng thuật toán chẩn đoán tự động kết hợp phương pháp dò tìm vết nứt với kỹ thuật hiệu chỉnh toán học cao cấp. Phạm vi nghiên cứu tập trung vào các kết cấu thanh một chiều với các điều kiện biên kinh điển (ngàm - tự do, ngàm - ngàm, tự do - tự do) và liên kết đàn hồi tổng quát, được triển khai và hoàn thiện trong giai đoạn năm 2016 tại Đại học Công nghệ, Đại học Quốc gia Hà Nội.

Ý nghĩa khoa học và thực tiễn của nghiên cứu được lượng hóa rõ nét: phương pháp mới giúp giảm thiểu sai số định vị vết nứt xuống dưới 2% trên lưới quét 100 điểm chia, đồng thời duy trì độ ổn định tính toán ngay cả khi dữ liệu đo đạc thực nghiệm chứa mức độ nhiễu lên tới 10%. Kết quả này mở ra hướng tiếp cận chi phí thấp, tối ưu hóa công tác kiểm tra không phá hủy cho các công trình hạ tầng trọng yếu.

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 xây dựng trên nền tảng cơ học dao động liên tục và cơ học phá hủy đàn hồi tuyến tính. Lý thuyết dao động dọc trục của thanh thẳng được mô hình hóa thông qua nguyên lý D'Alembert và định luật Hooke, dẫn xuất từ phương trình truyền sóng một chiều với vận tốc truyền sóng âm đặc trưng phụ thuộc vào mô đun đàn hồi $E$ và khối lượng riêng $\rho$.

Khuyết tật vết nứt ngang được mô hình hóa bằng phần tử lò xo dọc trục tương đương theo mô hình cơ học của Chondros và Dimarogonas, trong đó độ cứng lò xo $K_j$ đại diện cho sự suy giảm độ cứng cục bộ và liên hệ trực tiếp với độ sâu tương đối của vết nứt $a/h$.

Khung phân tích tích hợp 4 khái niệm cốt lõi:

  1. Tần số dao động riêng ($\omega$) và hàm dạng dao động ($U(x)$).
  2. Hàm đáp ứng tần số (Frequency Response Function - FRF) thể hiện phản ứng động lực học của kết cấu dưới kích động điều hòa.
  3. Điểm nút tần số - vị trí đặc biệt mà vết nứt xuất hiện nhưng không làm thay đổi giá trị của một tần số riêng cụ thể.
  4. Ma trận truyền (Transfer Matrix Method - TMM) và kỹ thuật xấp xỉ tiệm cận tham số bé.

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

Nghiên cứu sử dụng nguồn dữ liệu mô phỏng số học chính xác cao được đối sánh trực tiếp với nghiệm giải tích và phương pháp ma trận truyền trên tập hợp 5 bậc tần số dao động riêng đầu tiên. Cỡ mẫu tính toán khảo sát bao gồm 100 điểm nút rời rạc phân bố đều dọc theo chiều dài chuẩn hóa $L = 1.0$ của thanh, kết hợp cùng các kịch bản thử nghiệm đa dạng gồm thanh có 1 vết nứt, 2 vết nứt và 3 vết nứt tại các tọa độ tương đối khác nhau ($e = 0.2; 0.5; 0.7$).

Phương pháp chọn mẫu áp dụng kỹ thuật lưới quét đồng nhất dọc trục (Crack Scanning Method), chia thanh thành 100 đoạn bằng nhau để đảm bảo không bỏ sót bất kỳ vùng khuyết tật tiềm năng nào trong không gian khảo sát.

Lý do lựa chọn phương pháp giải tích kết hợp phương pháp điều chỉnh Tikhonov và phép khai triển giá trị kỳ dị (Singular Value Decomposition - SVD) là vì bài toán ngược chẩn đoán hư hại luôn thuộc dạng bài toán phi chỉnh (ill-posed): số lượng phương trình đo đạc thực tế ($m = 5$) nhỏ hơn rất nhiều so với số ẩn số cần tìm ($n = 100$ điểm lưới quét), đồng thời dữ liệu đo luôn bị ảnh hưởng bởi sai số thiết bị. Kỹ thuật điều chỉnh Tikhonov cho phép loại trừ tính mất ổn định của ma trận hệ số, mang lại nghiệm hội tụ đơn nhất và ổn định. Toàn bộ quy trình tính toán được hoàn thành trong chu kỳ nghiên cứu 24 tháng với 3 giai đoạn phân tích chặt chẽ.

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

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

Thứ nhất, nghiên cứu đã thiết lập thành công phương trình tần số giải tích dạng tường minh cho thanh có số lượng vết nứt tùy ý với điều kiện biên tổng quát. Kết quả kiểm chứng trên 5 trị riêng đầu tiên cho thấy phương pháp giải tích đạt độ trùng khớp tuyệt đối, sai số dưới 0.01% khi so sánh với phương pháp ma trận truyền (TMM).

Thứ hai, sự biến thiên của tần số riêng phụ thuộc phi tuyến vào độ sâu và vị trí vết nứt. Khi độ sâu vết nứt tăng từ $a/h = 0.1$ lên $a/h = 0.5$, tần số dao động riêng thứ nhất của thanh ngàm - tự do giảm mạnh tới 18.2%, trong khi các bậc tần số cao hơn có mức suy giảm dao động từ 5.4% đến 14.6% tùy thuộc vào dạng mode dao động.

Thứ ba, nghiên cứu đã lập bảng xác định chính xác các điểm nút tần số: đối với thanh hai đầu ngàm, vết nứt tại vị trí chính giữa $e = 0.5$ hoàn toàn không làm thay đổi tần số riêng bậc 1 và bậc 3, nhưng làm suy giảm đáng kể tần số riêng bậc 2 khoảng 12.8%.

Thứ tư, khi áp dụng giải thuật dò quét kết hợp điều chỉnh Tikhonov trên lưới 100 điểm với 5 tần số đầu tiên, thuật toán đã cô lập chính xác tọa độ của 2 vết nứt đồng thời tại $e_1 = 0.3$ và $e_2 = 0.7$, ước lượng độ sâu vết nứt đạt độ tin cậy trên 95.8%, triệt tiêu hoàn toàn các nghiệm giả hoặc giá trị âm không hợp lý.

Thảo luận kết quả

Cơ chế vật lý của hiện tượng suy giảm tần số bắt nguồn từ sự tập trung ứng suất và suy giảm độ cứng cục bộ tại vị trí vết nứt. Tại các vùng có biên độ biến dạng dọc trục cực đại, ảnh hưởng của vết nứt lên tần số là lớn nhất; ngược lại, tại các điểm nút dao động (nơi biến dạng bằng 0), sự hiện diện của vết nứt hầu như không gây ra bất kỳ biến thiên tần số nào.

So với các nghiên cứu kinh điển trước đây vốn chỉ giải quyết trường hợp một vết nứt đơn lẻ hoặc yêu cầu dạng dao động đầy đủ, phương pháp đề xuất vượt trội nhờ khả năng nhận dạng đồng thời nhiều vết nứt chỉ bằng một số lượng hạn chế các tần số riêng đo được.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị đường cong phân tán 2D biểu diễn quan hệ giữa tỷ số tần số $\omega/\omega_0$ theo vị trí vết nứt $e \in [0, 1]$, kết hợp cùng bảng ma trận độ nhạy so sánh 5 trị riêng giải tích với sai số chuẩn hóa. Các biểu đồ biên độ đáp ứng tần số FRF thể hiện rõ nét các đỉnh cộng hưởng bị dịch chuyển về phía tần số thấp và sự biến đổi rõ rệt của các tần số phản cộng hưởng khi vết nứt xuất hiện và phát triển.

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

  1. Tích hợp thuật toán dò quét vết nứt và phương pháp điều chỉnh Tikhonov vào các hệ thống quan trắc sức khỏe công trình (Structural Health Monitoring - SHM) theo thời gian thực; hướng tới mục tiêu giảm sai số nhận dạng vị trí khuyết tật xuống dưới 2.5%, triển khai thử nghiệm trong vòng 6 tháng bởi các kỹ sư kết cấu và đơn vị quản lý vận hành cầu đường.

  2. Tiêu chuẩn hóa quy trình thu thập dữ liệu dao động hiện trường bằng cách đo đạc tối thiểu 5 mode tần số riêng đầu tiên kết hợp cảm biến gia tốc dải tần rộng; hướng tới chỉ tiêu lọc bỏ trên 90% nhiễu trắng môi trường, thực hiện trong khung thời gian 3 tháng bởi các trung tâm kiểm định chất lượng công trình xây dựng.

  3. Xây dựng và hoàn thiện bộ công cụ phần mềm tính toán tự động hóa giải thuật ma trận SVD trên các lưới quét mật độ cao từ 100 đến 500 điểm phân chia; đặt mục tiêu rút ngắn thời gian xử lý dữ liệu chẩn đoán xuống dưới 5 giây cho mỗi chu kỳ kiểm tra, hoàn thành trong vòng 12 tháng bởi các nhóm nghiên cứu công nghệ tính toán chuyên sâu.

  4. Mở rộng mô hình lý thuyết từ kết cấu thanh một chiều sang hệ kết cấu dầm chịu uốn - xoắn phức hợp và kết cấu dàn thép không gian đa liên kết đàn hồi; hướng tới nâng cao phạm vi bao phủ chẩn đoán lên 80% các dạng kết cấu chịu lực thực tế, lộ trình thực hiện trong 18 tháng bởi các viện nghiên cứu chuyên ngành cơ học công trình.

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

  • Kỹ sư kết cấu và chuyên gia giám định công trình: Nắm vững phương pháp kiểm tra không phá hủy (NDT) dựa trên dao động, ứng dụng trực tiếp vào việc đánh giá độ toàn vẹn của cọc móng sâu, thanh giàn thép và kết cấu cầu cảng.
  • Giảng viên và nghiên cứu viên ngành Cơ kỹ thuật, Kỹ thuật Công trình: Khai thác hệ phương trình tần số dạng tường minh và mô hình lò xo tương đương làm tài liệu tham khảo giảng dạy và phát triển các bài toán dao động nâng cao.
  • Học viên cao học và sinh viên chuyên ngành Cơ học tính toán, Xây dựng: Tiếp cận phương pháp luận mẫu mực về giải quyết bài toán ngược phi chỉnh bằng công cụ toán học hiện đại như khai triển SVD và hiệu chỉnh Tikhonov.
  • Các công ty công nghệ phát triển giải pháp IoT và cảm biến thông minh: Ứng dụng thuật toán tối ưu hóa vào việc đóng gói các thuật toán chẩn đoán tự động trên chip vi xử lý của thiết bị giám sát công trình liên tục.

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

Phương pháp tần số riêng có thể phát hiện được vết nứt ở giai đoạn rất sớm không? Có, phương pháp hoàn toàn có thể phát hiện vết nứt khi độ sâu tương đối mới đạt mức khoảng 5% đến 10% chiều cao tiết diện. Nhờ vào thuật toán điều chỉnh Tikhonov kết hợp với việc khảo sát đồng thời 5 bậc tần số riêng đầu tiên, các biến thiên rất nhỏ của tần số (dưới 1%) vẫn được nhận dạng chính xác mà không bị chìm trong nhiễu đo đạc.

Tại sao việc chẩn đoán vết nứt lại được xếp vào nhóm bài toán ngược phi chỉnh? Bài toán mang tính chất phi chỉnh vì số lượng thông số cần tìm (vị trí và độ sâu của hàng trăm điểm tiềm năng) vượt xa số lượng dữ liệu đo đạc thực tế (thường chỉ đo được 3 đến 5 tần số). Khi ma trận hệ số không vuông kết hợp với sai số đo đạc thực tế, nghiệm của phương trình đại số rất dễ bị phân kỳ hoặc không ổn định nếu không có công cụ điều chỉnh.

Điểm nút tần số là gì và xử lý thế nào nếu vết nứt nằm đúng điểm nút? Điểm nút tần số là vị trí trên thanh mà khi vết nứt xuất hiện tại đó, tần số riêng của một bậc dao động nhất định không hề bị thay đổi. Để tránh bỏ sót khuyết tật, nghiên cứu chỉ ra rằng cần đo đạc phối hợp ít nhất 3 bậc dao động khác nhau, bởi vì điểm nút của bậc này sẽ trùng với vùng biến dạng cực đại của bậc dao động khác.

Điều kiện biên của thanh ảnh hưởng như thế nào đến độ nhạy chẩn đoán? Điều kiện biên quyết định trực tiếp quy luật phân bố ứng suất dọc trục. Ví dụ, thanh ngàm - tự do có độ nhạy nứt cao nhất tại vùng sát đầu ngàm (vị trí $x = 0$) và giảm dần về đầu tự do ($x = 1.0$), trong khi thanh ngàm hai đầu lại có độ nhạy cao ở cả hai đầu liên kết và suy giảm cục bộ tại điểm chính giữa.

Cần tối thiểu bao nhiêu tần số dao động riêng để xác định chính xác một vết nứt? Trong điều kiện lý tưởng không có nhiễu, chỉ cần đo đạc chính xác 2 tần số dao động riêng là đủ để xác định tọa độ và độ sâu của một vết nứt đơn lẻ. Tuy nhiên, trong điều kiện thực nghiệm có sai số đo đạc, việc sử dụng từ 3 đến 5 tần số riêng sẽ đảm bảo kết quả nhận dạng đạt độ tin cậy trên 95%.

Kết luận

  • Thiết lập thành công hệ phương trình tần số dạng tường minh cho dao động dọc trục của thanh đàn hồi mang số lượng vết nứt bất kỳ với điều kiện liên kết tổng quát.
  • Xây dựng hoàn chỉnh thuật toán dò quét khuyết tật kết hợp phương pháp điều chỉnh Tikhonov và kỹ thuật SVD, khắc phục triệt để tính chất phi chỉnh của bài toán ngược.
  • Xác lập bản đồ phân bố các điểm nút tần số trên các kết cấu thanh ngàm - tự do, ngàm - ngàm và tự do - tự do, cung cấp cơ sở nhận dạng sơ bộ vết nứt nhanh chóng.
  • Kiểm chứng số học trên lưới 100 điểm phân chia khẳng định độ chính xác vượt trội với sai số định vị dưới 2% ngay cả trong điều kiện dữ liệu đo chứa nhiễu.
  • Đóng góp giải pháp khoa học thực tiễn giúp nâng cao hiệu quả kiểm định không phá hủy và giảm thiểu rủi ro sự cố kết cấu trong xây dựng và cơ khí.
  • Kế hoạch phát triển tiếp theo tập trung vào việc thử nghiệm thực nghiệm trên mô hình cọc khoan nhồi thực tế và mở rộng giải thuật cho hệ khung không gian trong lộ trình 12 tháng tới.
  • Hãy tham khảo toàn văn luận văn để ứng dụng ngay các công thức giải tích chuẩn xác và giải thuật tối ưu vào các dự án nghiên cứu cũng như công tác kiểm định công trình thực tế!