Tổng quan nghiên cứu

Trong kỹ thuật kết cấu công trình dân dụng và công nghiệp, các khuyết tật nứt tế vi chiếm hơn 70% nguyên nhân dẫn đến suy giảm khả năng chịu lực dài hạn và đe dọa an toàn công trình. Việc phát hiện sớm các vết nứt này là bài toán cấp bách, nhưng các phương pháp phần tử hữu hạn truyền thống (FEM) thường gặp rào cản lớn khi mô phỏng kích thích tần số cao trong dải Kilohertz (từ 10 kHz đến hàng trăm kHz) do đòi hỏi chia lưới cực mịn làm bùng nổ khối lượng tính toán.

Luận văn thạc sĩ kỹ thuật công trình của tác giả Nguyễn Thanh Tú, thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh năm 2015 dưới sự hướng dẫn của PGS. Nguyễn Thị Hiền Lương và TS. Bùi Quốc Tính, tập trung giải quyết bài toán truyền sóng trong dầm Euler-Bernoulli có vết nứt ngang thông qua phương pháp phần tử hữu hạn phổ dựa trên Wavelet (Wavelet Spectral Finite Element Method – WSFEM). Mục tiêu cốt lõi của nghiên cứu là xây dựng công cụ số có độ chính xác cao để giải bài toán ngược xác định vị trí và chiều sâu vết nứt, phục vụ trực tiếp cho công tác chẩn đoán và giám sát sức khỏe kết cấu (SHM).

Nghiên cứu khảo sát toàn diện trên kết cấu dầm console và dầm tựa đơn trong các cửa sổ thời gian từ 512 microgiây đến 2048 microgiây. Ý nghĩa khoa học và thực tiễn của đề tài thể hiện ở việc nâng cao hiệu suất phân tích động lực học kết cấu, kiểm soát sai số nhận diện vị trí vết nứt ở mức dưới 1%, đồng thời tiết kiệm hơn 80% tài nguyên tính toán so với các phương pháp phần tử hữu hạn 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 tích hợp của ba khung lý thuyết chuyên sâu:

Thứ nhất là lý thuyết biến đổi Wavelet trực giao Daubechies (db), cụ thể là các bậc Daubechies N = 4, 6, 12 và 22 tương ứng với số moment triệt tiêu từ 2 đến 11. Các hàm tỉ lệ (father wavelet) và hàm Wavelet mẹ (mother wavelet) được sử dụng để rời rạc hóa các biến trạng thái trên miền thời gian, cho phép biến đổi và cục bộ hóa tín hiệu hiệu quả mà không làm mất thông tin thời gian như biến đổi Fourier truyền thống.

Thứ hai là lý thuyết dầm Euler-Bernoulli mở rộng kết hợp cơ học rạn nứt đàn hồi tuyến tính (LEFM). Vết nứt ngang trong dầm được mô hình hóa bằng liên kết lò xo xoay tương đương với độ mềm cục bộ phát sinh được tính toán trực tiếp từ hệ số cường độ ứng suất và hàm mật độ giải phóng năng lượng biến dạng.

Thứ ba là lý thuyết phân tích phổ và động lực học sóng đàn hồi. Các phương trình đạo hàm riêng (PDE) mô tả chuyển vị dọc trục u(x,t) và chuyển vị ngang uốn w(x,t) được phân tích thành hệ phương trình vi phân thường (ODE), từ đó xác định các thông số sóng cốt lõi bao gồm số sóng (wavenumber), vận tốc pha, vận tốc nhóm và các chế độ lan truyền hoặc tắt dần tương ứng với tần số cắt.

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

Nguồn dữ liệu của luận văn được thu thập từ quá trình mô phỏng số chi tiết trên 7 bài toán khảo sát đặc thù, bao gồm dầm không nứt, dầm có vết nứt biến thiên độ sâu từ 10% đến 50% chiều cao tiết diện, và dầm có vết nứt biến thiên vị trí tại 30% và 50% chiều dài dầm. Cỡ mẫu khảo sát gồm hơn 20 kịch bản tương tác sóng với các bước thời gian lấy mẫu delta t lần lượt là 1 microgiây, 2 microgiây và 4 microgiây.

Phương pháp chọn mẫu dựa trên nguyên tắc biến thiên đơn tham số nhằm cô lập và đánh giá chính xác tác động của từng yếu tố: bậc Wavelet Daubechies (từ N=4 đến N=22), chiều sâu vết nứt và cửa sổ thời gian quan sát (từ 512 đến 2048 microgiây). Tác giả lựa chọn phương pháp phân tích WSFEM thay vì phương pháp phần tử hữu hạn thông thường vì WSFEM có khả năng thiết lập ma trận độ cứng động lực học chính xác từ nghiệm giải tích, giải quyết triệt để bài toán kích thích tần số cao mà không bị giới hạn bởi hiện tượng tán sắc số hay kích thước phần tử.

Toàn bộ quy trình tính toán được thực hiện trong thời gian 6 tháng (từ tháng 7/2014 đến tháng 12/2014), sử dụng thuật toán phân tích giá trị đơn (SVD) và phân tích trị riêng để giải hệ phương trình truyền sóng và đối chuẩn trực tiếp với phương pháp phần tử hữu hạn phổ Laplace (LSFEM).

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 và thực nghiệm số mang lại 4 phát hiện quan trọng:

Một là, bậc của hàm tỉ lệ Wavelet Daubechies đóng vai trò quyết định đến độ chính xác của phản ứng động lực học. Khi sử dụng Daubechies bậc cao N = 22 với cửa sổ thời gian 1024 microgiây, đáp ứng vận tốc sóng dọc và sóng ngang triệt tiêu hoàn toàn các dao động sóng giả ở biên, cải thiện độ chính xác hơn 99,5% so với bậc N = 6.

Hai là, quan hệ giữa biên độ sóng phản xạ và độ sâu vết nứt tuân theo quy luật biến thiên rõ rệt. Khi tỷ lệ độ sâu vết nứt hd/h tăng từ 20% lên 50%, biên độ xung sóng phản xạ ghi nhận tại đầu tự do của dầm console tăng từ 35% lên hơn 60% so với biên độ sóng tới ban đầu.

Ba là, phương pháp WSFEM giải quyết thành công bài toán ngược nhận diện hư hại. Với bước dịch chuyển lấy mẫu delta x = 4 x 10^-7 m, sai số xác định vị trí vết nứt tại tọa độ Ld = 0.3L và Ld = 0.5L đều đạt mức dưới 0,5%, chứng minh độ tin cậy vượt trội trong việc phát hiện khuyết tật hình học.

Bốn là, phân tích phổ làm rõ sự khác biệt bản chất giữa sóng dọc và sóng ngang. Sóng dọc lan truyền với số sóng tuyến tính theo tần số nên là sóng không tán sắc với vận tốc không đổi; ngược lại, sóng uốn ngang thể hiện tính tán sắc mạnh mẽ, phụ thuộc trực tiếp vào tần số cắt và hình học mặt cắt dầm.

Thảo luận kết quả

Về mặt cơ chế vật lý, sự hiện diện của vết nứt làm gia tăng độ mềm cục bộ, gây ra sự gián đoạn trở kháng cơ học trong dầm Euler-Bernoulli. Khi sóng tới gặp mặt cắt nứt, một phần năng lượng bị phản xạ ngược lại và phần còn lại tiếp tục truyền qua, tạo nên các xung tín hiệu đặc trưng trên biểu đồ miền thời gian.

Các kết quả nghiên cứu được trực quan hóa thông qua hệ thống biểu đồ đáp ứng vận tốc sóng dọc và sóng ngang theo thời gian từ 0 đến 2048 microgiây, kết hợp các bảng đối chiếu sai số định vị. Dữ liệu thể hiện rõ ràng các đỉnh xung sóng tới và sóng phản xạ, giúp kỹ sư dễ dàng trích xuất khoảng thời gian di chuyển của sóng để tính toán khoảng cách đến khuyết tật.

Khi so sánh với các công bố khoa học quốc tế sử dụng phương pháp LSFEM và FEM 2D, kết quả từ WSFEM cho độ tương thích hình dạng sóng đạt trên 98%. Tuy nhiên, phương pháp WSFEM chỉ cần 1 phần tử phổ cho một đoạn dầm liên tục thay vì hàng nghìn phần tử như FEM 2D, giúp giảm hơn 80% thời gian tính toán và loại bỏ hoàn toàn hiện tượng phản xạ sóng nhân tạo tại biên mô hình.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 khuyến nghị và giải pháp hành động cụ thể:

Thứ nhất, tích hợp thuật toán WSFEM vào các hệ thống quan trắc sức khỏe kết cấu (SHM) tự động. Các doanh nghiệp công nghệ xây dựng và đơn vị quản lý hạ tầng cần xây dựng mô đun xử lý tín hiệu cảm biến áp điện PZT trong dải tần 50 kHz đến 200 kHz, đặt mục tiêu phát hiện sớm các vết nứt có tỷ lệ hd/h từ 10% trong vòng 3 đến 6 tháng triển khai.

Thứ hai, chuẩn hóa quy trình đo đạc sóng đàn hồi tần số cao tại hiện trường. Các trung tâm kiểm định chất lượng công trình cần nâng cấp thiết bị thu nhận dữ liệu đạt tần số lấy mẫu tối thiểu 1 MHz (khoảng thời gian 1 microgiây) nhằm đảm bảo sai số định vị hư hại thực tế luôn duy trì dưới mức 1% trong lộ trình 6 đến 12 tháng tới.

Thứ ba, mở rộng phạm vi ứng dụng của mô hình phổ Wavelet cho kết cấu hai chiều và vật liệu tiên tiến. Các viện nghiên cứu chuyên ngành cần tiếp tục xây dựng công thức phần tử tấm mỏng và tấm composite nhiều lớp có khuyết tật tách lớp, hướng tới mục tiêu nâng tỷ lệ chẩn đoán chính xác đa vết nứt phức tạp lên trên 95% trong giai đoạn 2026–2028.

Thứ tư, kết hợp phương pháp WSFEM với các thuật toán học máy (Machine Learning) và giải thuật di truyền (GA). Các nhóm nghiên cứu liên ngành cần số hóa ma trận đáp ứng động lực học để huấn luyện mạng nơ-ron nhân tạo, tự động hóa bài toán ngược chẩn đoán đa khuyết tật với thời gian phản hồi dưới 5 giây.

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

Nội dung luận văn mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng chính:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Xây dựng, Cơ kỹ thuật và Cơ học kết cấu: Tài liệu cung cấp cơ sở toán học chặt chẽ về phân tích đa giải Wavelet Daubechies bậc N=22 và phương pháp thành lập ma trận độ cứng động lực học để phát triển các đề tài nghiên cứu chuyên sâu.

  2. Kỹ sư kiểm định và chẩn đoán công trình cầu đường, nhà cao tầng: Tài liệu đóng vai trò là cẩm nang ứng dụng phương pháp sóng siêu âm không phá hủy (NDT) để xác định chính xác vị trí và độ sâu nứt ngầm trong dầm bê tông cốt thép và dầm thép.

  3. Các chuyên gia và doanh nghiệp phát triển thiết bị giám sát kết cấu (SHM): Công trình cung cấp thuật toán tối ưu giúp giảm tải 80% bộ nhớ tính toán trên các thiết bị xử lý biên (edge computing) kết nối cảm biến rung động.

  4. Giảng viên và nhà nghiên cứu tại các trường đại học kỹ thuật: Luận văn là tài liệu tham khảo chất lượng cao gồm 7 chương lý thuyết và ví dụ số hoàn chỉnh, phục vụ trực tiếp cho công tác giảng dạy các môn Động lực học công trình nâng cao và Cơ học rạn nứt.

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

  1. Phương pháp WSFEM vượt trội hơn phương pháp FEM truyền thống ở những điểm nào? WSFEM giải quyết bài toán truyền sóng tần số cao từ 10 kHz đến hàng trăm kHz mà không cần chia lưới dày đặc. Trong khi FEM đòi hỏi hàng nghìn phần tử khiến thời gian tính toán kéo dài, WSFEM chỉ cần 1 phần tử cho mỗi đoạn dầm đồng nhất, giúp giảm hơn 80% chi phí tính toán và loại bỏ triệt để hiện tượng tán sắc số.

  2. Vì sao nghiên cứu lại ưu tiên sử dụng hàm Wavelet Daubechies bậc N = 22? Wavelet Daubechies bậc N = 22 sở hữu 11 moment triệt tiêu, cho phép xấp xỉ chính xác các hàm đa thức bậc cao và xử lý mượt mà tín hiệu có biến đổi đột ngột. So với bậc N = 6, bậc N = 22 trong cửa sổ thời gian 1024 microgiây giúp triệt tiêu hoàn toàn sóng phản xạ biên giả, đảm bảo độ chính xác trên 99,5%.

  3. Làm thế nào để xác định chính xác vị trí và độ sâu vết nứt từ tín hiệu sóng? Vị trí vết nứt được xác định dựa vào thời gian đến của xung sóng phản xạ nhân với vận tốc nhóm của sóng. Chiều sâu vết nứt trong khoảng 20% đến 50% chiều cao dầm được đánh giá thông qua mức độ suy giảm biên độ xung sóng tới và sự gia tăng biên độ của xung sóng phản xạ trên biểu đồ miền thời gian.

  4. Mô hình lò xo tương đương mô tả vết nứt trong dầm có độ tin cậy ra sao? Mô hình lò xo tương đương được thiết lập dựa trên ma trận độ mềm cục bộ từ lý thuyết cơ học rạn nứt đàn hồi tuyến tính (LEFM). Mô hình này phản ánh chính xác sự suy giảm độ cứng uốn cục bộ tại tiết diện nứt, cho kết quả truyền sóng trùng khớp với mô hình phần tử hữu hạn 2D với sai số dưới 1,5%.

  5. Phương pháp WSFEM có khả năng nhận diện dầm có nhiều vết nứt cùng lúc không? Có. Trong chương 6 của luận văn, tác giả đã thử nghiệm thành công bài toán dầm có 2 vết nứt đồng thời. Bằng cách phân tích chuỗi xung sóng phản xạ liên tiếp trên biểu đồ miền thời gian từ 0 đến 2048 microgiây, thuật toán WSFEM bóc tách chính xác từng điểm dị thường với sai số vị trí dưới 1%.

Kết luận

  • Luận văn xây dựng thành công phương pháp phần tử hữu hạn phổ Wavelet (WSFEM) để giải bài toán truyền sóng tần số cao trong dầm Euler-Bernoulli có vết nứt ngang.
  • Ma trận độ cứng động lực học được thiết lập chính xác thông qua mô hình lò xo tương đương dựa trên độ mềm cục bộ và cơ học rạn nứt đàn hồi tuyến tính.
  • Thuật toán giải quyết hiệu quả bài toán ngược, xác định chính xác vị trí nứt (Ld = 0.3L, 0.5L) và chiều sâu nứt (hd/h = 20% - 50%) với sai số dưới 1%.
  • Phương pháp giúp tiết kiệm hơn 80% thời gian tính toán và dung lượng bộ nhớ so với các mô hình phần tử hữu hạn truyền thống trong dải thời gian 512 đến 2048 microgiây.
  • Cung cấp cơ sở dữ liệu số hóa tin cậy và hệ thống 7 bài toán kiểm chứng từ dầm đơn giản đến dầm đa vết nứt phức tạp.

Công trình tạo ra bước tiến quan trọng trong việc ứng dụng công cụ toán học Wavelet vào lĩnh vực kỹ thuật xây dựng và giám sát công trình. Trong giai đoạn 2026–2028, các kỹ sư và nhà quản lý công trình cần nhanh chóng ứng dụng phương pháp WSFEM vào thực tế kiểm định để nâng cao độ an toàn và kéo dài tuổi thọ khai thác cho các kết cấu hạ tầng trọng điểm.