Tổng quan về luận án

Trong kỹ thuật cơ khí, giao thông và công trình xây dựng, việc chẩn đoán hư hỏng kết cấu nhằm duy trì khả năng làm việc liên tục và kéo dài tuổi thọ công trình đóng vai trò sống còn đối với an toàn công cộng và hiệu quả kinh tế. Trong các dạng khuyết tật thường gặp, vết nứt ngang cục bộ gây suy giảm độ cứng cục bộ nghiêm trọng, dẫn tới nguy cơ sụp đổ đột ngột nếu không được phát hiện kịp thời. Luận án tiến sĩ của tác giả Phạm Thị Ba Liên (thực hiện dưới sự hướng dẫn khoa học của GS. Nguyễn Tiến Khiêm và TS. Nguyễn Minh Tuấn tại Viện Cơ học, Viện Hàn lâm Khoa học và Công nghệ Việt Nam) là công trình tiên phong giải quyết bài toán ngược chẩn đoán đa vết nứt trong kết cấu thanh và dầm đàn hồi thông qua việc khai thác đồng thời tần số cộng hưởng (resonance frequencies) và tần số phản cộng hưởng (antiresonance frequencies).

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án xác định xuất phát từ giới hạn căn bản của các phương pháp chẩn đoán động lực học truyền thống. Tần số dao động riêng là đại lượng vô hướng toàn cục (global parameter), dễ đo đạc với độ chính xác cao nhưng lại kém nhạy cảm với các hư hỏng cục bộ và không mang lại nghiệm duy nhất khi kết cấu có điều kiện biên đối xứng (Morassi, 2001; Dilena & Morassi, 2004). Ngược lại, các đặc trưng hàm như dạng dao động riêng (mode shapes) hay hàm truyền đáp ứng tần số (Frequency Response Functions - FRF) rất nhạy với khuyết tật cục bộ nhưng lại phụ thuộc mật thiết vào vị trí đặt cảm biến, vị trí lực kích động, đồng thời rất khó đo đạc chính xác trong thực địa do ảnh hưởng của nhiễu và tương tác đa mode dao động. Luận án đã xác lập luận điểm đột phá: Tần số phản cộng hưởng—vốn là điểm không (zeros) của hàm đáp ứng tần số—chính là một đặc trưng số mang bản chất cục bộ, chứa đựng đồng thời tọa độ kích thích và vị trí đo, cho phép khắc phục triệt để nghịch lý giữa "tính dễ đo của đặc trưng số" và "tính nhạy cục bộ của đặc trưng hàm".

Các câu hỏi nghiên cứu và giả thuyết khoa học được xác định cụ thể:

  • RQ1: Làm thế nào để thiết lập phương trình tần số giải tích tường minh dạng đa thức đối với các tham số vết nứt trong dao động dọc trục của thanh và dao động uốn của dầm có nhiều vết nứt?
  • RQ2: Mở rộng công thức Rayleigh như thế nào để tính toán trực tiếp tần số riêng của thanh có nhiều vết nứt và xác định phạm vi tin cậy của các nghiệm xấp xỉ tiệm cận?
  • RQ3: Tần số phản cộng hưởng biến thiên như thế nào theo vị trí, chiều sâu vết nứt và làm sao xây dựng thuật toán chẩn đoán đa vết nứt ổn định từ hệ dữ liệu cộng hưởng - phản cộng hưởng?
  • H1: Phương trình tần số dạng đa thức bậc $n$ theo độ lớn vết nứt $\gamma_j$ cho phép xác định chính xác phổ tần số mà không cần giải hệ phương trình ma trận siêu việt phức tạp.
  • H2: Xấp xỉ tiệm cận bậc nhất của công thức Rayleigh đạt độ chính xác cao khi chiều sâu vết nứt dưới $20%$, trong khi xấp xỉ bậc hai và bậc ba duy trì độ tin cậy đến độ sâu $30% - 50%$.
  • H3: Việc kết hợp dữ liệu tần số phản cộng hưởng từ các điểm đo phân bố giải quyết triệt để tính chất không chỉnh (ill-posedness) của bài toán ngược, cho phép nhận dạng đồng thời vị trí và chiều sâu của đa vết nứt.

Khung lý thuyết của công trình tích hợp cơ học phá hủy đàn hồi tuyến tính (Linear Elastic Fracture Mechanics), lý thuyết dao động kết cấu một chiều Euler-Bernoulli và lý thuyết bài toán ngược không chỉnh Ti-khô-nôv (Tikhonov Regularization). Phạm vi nghiên cứu tập trung vào kết cấu thanh đàn hồi chịu dao động dọc trục và dầm chịu uốn thuần túy chứa từ 1 đến 19 vết nứt ngang hở, với độ sâu vết nứt khảo sát từ $10%$ đến $50%$ chiều cao tiết diện dưới đa dạng điều kiện biên (ngàm - ngàm, ngàm - tự do, tự do - tự do, gối tựa đơn). Ý nghĩa khoa học của nghiên cứu thể hiện qua việc khép lại khoảng trống giải tích tồn tại hàng thập kỷ trong việc tìm mối liên hệ tường minh giữa tham số khuyết tật và đặc trưng động lực học phản cộng hưởng.

Literature Review và Positioning

Lịch sử chẩn đoán hư hỏng kết cấu dựa trên dao động (Vibration-Based Structural Health Monitoring) đã trải qua nhiều giai đoạn phát triển quan trọng, được tổng kết toàn diện qua các công trình tổng quan kinh điển của Doebling và cộng sự (1996), Sohn và cộng sự (2001), Fan và Qiao (2011), và gần nhất là Hou và cộng sự (2019). Ngay từ năm 1978, Adams và các cộng sự đã đặt nền móng cho việc chẩn đoán một vết nứt trong thanh bằng mô hình lò xo dọc trục tương đương thông qua phương trình tần số dạng $d_0(\omega) + \gamma d_1(\omega, e) = 0$. Tuy nhiên, khi mở rộng cho thanh có $n$ vết nứt, mãi đến năm 2004, Ruotolo và Surace mới thiết lập được phương trình tần số ở dạng định thức $\det[A(\omega, e_1, \dots, \gamma_n)] = 0$ với ma trận cấp 4 chứa các phần tử siêu việt cực kỳ phức tạp. Trong kết cấu dầm chịu uốn, các nghiên cứu của Shifrin và Bratus (1999), Morassi (2001), và Khiêm và các cộng sự (2002, 2011) đã cố gắng giải quyết bài toán thuận và nghịch, song việc chẩn đoán đa vết nứt bằng tần số riêng vẫn vấp phải rào cản toán học do hệ phương trình siêu việt không có dạng hiển, đòi hỏi phải tuyến tính hóa với giả thiết vết nứt rất nhỏ ($\gamma \ll 1$).

Một cuộc tranh luận học thuật kéo dài giữa hai trường phái chẩn đoán: Trường phái dựa trên triệu chứng (symptom-based/data-driven) sử dụng xử lý tín hiệu hiện đại (Wavelet, Hilbert-Huang, Học máy) có ưu điểm phát hiện bất thường nhanh nhưng không định lượng được kích thước vật lý của khuyết tật và rất nhạy cảm với nhiễu môi trường; trong khi Trường phái dựa trên mô hình (model-based) cho phép xác định chính xác tọa độ hình học và độ sâu vết nứt nhưng lại phụ thuộc vào việc xây dựng quan hệ giải tích chính xác giữa hư hỏng và phản ứng động lực học. Luận án của Phạm Thị Ba Liên định vị vững chắc trong trường phái tiếp cận dựa trên mô hình, giải quyết triệt để sự mâu thuẫn này bằng cách cung cấp các công thức giải tích tường minh.

Nghiên cứu về tần số phản cộng hưởng bắt đầu từ các bài toán điều khiển dao động (Miu, 1993) và cập nhật mô hình (Friswell & Mottershead, 1995). Tiếp đó, Wang và cộng sự (2001) đã phân tích có hệ thống tần số phản cộng hưởng của dầm công-xôn và chỉ ra rằng khi vị trí điểm đo hoặc điểm đặt lực trùng với điểm nút của một dạng riêng, tần số phản cộng hưởng sẽ trùng với tần số riêng tương ứng. Bamnios và các cộng sự (2002) chứng minh vị trí điểm đặt lực trong phép đo trở kháng cơ học (mechanical impedance) là nguồn dữ liệu quý giá để định vị vết nứt. Dilena và Morassi (2004) cùng Rubio và cộng sự (2006) khẳng định việc tích hợp tần số phản cộng hưởng giúp khắc phục tính đối xứng của nghiệm chẩn đoán. Đáng chú ý nhất, Meruane và Heylen (2010) đã kết luận rằng:

"sử dụng tần số phản cộng hưởng cùng với tần số riêng để chẩn đoán hư hỏng kết cấu hiệu quả hơn việc sử dụng tần số riêng cùng với dạng dao động riêng."

Luận án tiến sĩ này đã tiến xa hơn các công bố quốc tế nêu trên bằng cách: (1) Thay thế các phương trình định thức ma trận cồng kềnh của Ruotolo và Surace (2004) bằng phương trình đại số đa thức tường minh đối với các tham số vết nứt; (2) Mở rộng định lý Rayleigh cổ điển (Fernandez-Saez et al., 1999) cho thanh có vô số vết nứt dọc trục; (3) Thiết lập thuật toán giải bài toán ngược kết hợp khai triển kỳ dị ma trận (SVD) và phương pháp điều chỉnh Ti-khô-nôv cho phép nhận dạng vết nứt với độ chính xác cao ngay cả khi số điểm đo thực nghiệm bị giới hạn.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án tạo ra bước tiến quan trọng cho cơ học kết cấu và động lực học công trình thông qua việc mở rộng các lý thuyết cơ bản sau:

  1. Mở rộng lý thuyết Adams - Ruotolo - Surace: Luận án đã chứng minh dạng tổng quát của phương trình đặc trưng tần số dao động dọc trục của thanh có $n$ vết nứt dưới dạng đa thức tường minh bậc $n$ theo các tham số vết nứt $\gamma_j$: $$D(\omega) = d_0(\omega) + \sum_{j=1}^n \gamma_j d_1(\omega, e_j) + \sum_{j=2}^n \sum_{k=1}^{j-1} \gamma_j \gamma_k d_2(\omega, e_j, e_k) + \dots + (\gamma_1 \dots \gamma_n) d_n(\omega, e_1, \dots, e_n) = 0$$ Cấu trúc đa thức này giúp loại bỏ hoàn toàn nhu cầu tính toán các định thức ma trận phức tạp, cung cấp nền tảng giải tích hoàn chỉnh để tính toán độ nhạy dao động mà không cần giả định vết nứt vô cùng bé.

  2. Thiết lập công thức Rayleigh mở rộng cho thanh có nhiều vết nứt: Mở rộng thương số Rayleigh từ dầm uốn sang dao động dọc trục của thanh có khuyết tật gián đoạn: $$\omega^2 / \omega_0^2 = \frac{1 - \sum_{j=1}^n \gamma_j \phi_0^2(e_j)}{1 - 2\sum_{j=1}^n \gamma_j \phi_0^2(e_j) + \omega_0^2 \sum_{j,k=1}^n a_{jk} \gamma_j \gamma_k \phi_0'(e_j) \phi_0'(e_k)}$$ Khai triển chuỗi Taylor công thức này thiết lập mối liên hệ trực tiếp giữa tần số riêng và các bậc xấp xỉ tiệm cận, chứng minh xấp xỉ bậc nhất hoàn toàn tương thích với các kết quả cổ điển của Morassi (2001) và Shifrin (1999).

  3. Xác lập mô hình giải tích cho tần số phản cộng hưởng: Luận án chứng minh rằng tần số phản cộng hưởng $\omega_a$ của kết cấu thanh và dầm có vết nứt chính là nghiệm của định thức ma trận hàm đáp ứng tần số $\det[H(\omega)] = 0$, tương đương với tần số dao động riêng của một hệ kết cấu phụ bị bổ sung các liên kết ngàm ảo tại điểm đặt lực hoặc điểm đo phản ứng.

graph TD
    A["Mô hình Thanh/Dầm có Đa Vết nứt"] --> B["Cơ học Phá hủy: Lò xo tương đương (Kj, γj)"]
    B --> C["Hàm dạng dao động tổng quát Φ(x)"]
    C --> D["Phương trình đặc trưng đa thức tường minh"]
    C --> E["Công thức Rayleigh mở rộng"]
    D --> F["Tần số cộng hưởng (Resonance)"]
    D --> G["Tần số phản cộng hưởng (Antiresonance)"]
    F & G --> H["Hệ phương trình phi tuyến chẩn đoán"]
    H --> I["Bài toán ngược: Điều chỉnh Ti-khô-nôv + SVD"]
    I --> J["Nhận dạng Tọa độ (ej) và Chiều sâu (aj/h)"]

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp chặt chẽ ba trụ cột lý thuyết: Cơ học biến dạng đàn hồi môi trường liên tục, Cơ học phá hủy đàn hồi tuyến tính (mô hình hóa vết nứt bằng lò xo độ cứng $K_j = EA / (L\gamma_j)$ với tham số không thứ nguyên $\gamma_j = 2(1-\nu^2)(h/L)f(a_j/h)$), và Lý thuyết bài toán ngược phi tuyến.

Tính độc đáo của khung phân tích thể hiện ở việc chuyển hóa bài toán tìm điểm không của trường hàm sóng phức tạp thành bài toán giải đại số tuyến tính có điều kiện điều chỉnh. Luận án xác định rõ các điều kiện biên tổng quát tại hai đầu biên $x=0$ và $x=1$ thông qua các tham số $\alpha_0, \beta_0, \alpha_1, \beta_1$, bao quát toàn bộ các liên kết thực tế: ngàm - ngàm ($\alpha_0=\alpha_1=1, \beta_0=\beta_1=0$), ngàm - tự do ($\alpha_0=1, \beta_0=0, \alpha_1=0, \beta_1=1$), tự do - tự do ($\alpha_0=\alpha_1=0, \beta_0=\beta_1=1$), cũng như các liên kết đàn hồi biên thông qua lò xo dọc trục $K_0, K_1$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo thế giới quan thực chứng (positivism) và cơ học duy lý, kết hợp giải tích toán học thuần túy với mô phỏng số học nâng cao. Thiết kế nghiên cứu đa cấp độ (multi-level design) bao gồm:

  • Cấp độ 1: Khảo sát bài toán thuận dao động dọc trục của thanh đồng nhất chứa từ 1 đến 9 vết nứt nhằm kiểm chứng độ hội tụ và sai số của các nghiệm giải tích so với nghiệm chính xác.
  • Cấp độ 2: Khảo sát bài toán thuận dao động uốn của dầm Euler-Bernoulli chứa từ 1 đến 19 vết nứt, phân tích tính chất phân bố của các điểm nút tần số (frequency nodes).
  • Cấp độ 3: Khảo sát bài toán ngược chẩn đoán đa vết nứt bằng cách thiết lập ma trận nhạy cảm và giải thuật nhận dạng số trên nền tảng MATLAB.
flowchart LR
    subgraph ThietKe [Thiết kế nghiên cứu]
        M1["Mô hình giải tích chính xác"] --> M2["Xấp xỉ tiệm cận bậc 1, 2, 3"]
        M2 --> M3["Mô phỏng số MATLAB"]
    end
    subgraph ThuNghiem [Kịch bản thử nghiệm]
        T1["Thanh: 4, 8, 9 vết nứt"]
        T2["Dầm: 1, 2, 3, 19 vết nứt"]
        T3["Độ sâu nứt: 10% - 50%"]
    end
    subgraph ChanDoan [Giải thuật chẩn đoán ngược]
        C1["Hệ điểm đo: 10 - 20 cảm biến"]
        C2["Ma trận kỳ dị SVD"]
        C3["Điều chỉnh Ti-khô-nôv"]
    end
    ThietKe --> ThuNghiem --> ChanDoan

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ tính chặt chẽ học thuật cao nhất:

  1. Rời rạc hóa tương đương: Vết nứt ngang không lan truyền được mô hình hóa qua ma trận mềm dẻo cục bộ, liên kết hai đoạn thanh/dầm liên tục bằng các điều kiện tương thích ứng lực và chuyển vị gián đoạn: $$u(e_j + 0) - u(e_j - 0) = \gamma_j u'(e_j), \quad u'(e_j + 0) = u'(e_j - 0)$$
  2. Thiết lập hàm dạng tổng quát: Xây dựng nghiệm dạng truy hồi $\Phi_{j+1}(x) = \Phi_j(x) + \gamma_j \Phi_j(e_j) \cos[\lambda(x - e_j)]$ thỏa mãn đồng thời phương trình vi phân dao động và điều kiện biên trái.
  3. Giải phương trình siêu việt: Ứng dụng phương pháp dây cung (secant method) lập trình trong môi trường MATLAB để trích xuất nghiệm tần số riêng chính xác và tần số phản cộng hưởng với sai số hội tụ khống chế dưới $10^{-8}$.
  4. Khử tính không chỉnh của bài toán ngược: Do số lượng vết nứt thực tế $n$ chưa biết trước và số lượng cảm biến đo $M$ luôn hữu hạn, bài toán ngược dẫn đến hệ phương trình đại số tuyến tính kỳ dị $A \mu = b$. Luận án đã tích hợp kỹ thuật phân tích khai triển giá trị kỳ dị (SVD) kết hợp chuẩn hóa Ti-khô-nôv (Tikhonov Regularization) để lọc bỏ thành phần nhiễu và tìm nghiệm vector vết nứt ổn định nhất.

Data và phân tích

Thông số hình học và cơ lý của mô hình khảo sát bao gồm: Mô đun đàn hồi $E = 200\text{ GPa}$, khối lượng riêng $\rho = 7850\text{ kg/m}^3$, hệ số Poisson $\nu = 0.3$, chiều dài thanh/dầm $L = 1.0\text{ m}$, tiết diện hình chữ nhật $b \times h = 0.02\text{ m} \times 0.02\text{ m}$.

Các thử nghiệm chẩn đoán ngược được khảo sát với cấu hình dày đặc:

  • Số lượng điểm đo phân bố: Khảo sát lưới đo $M \in {10, 13, 14, 15, 16, 18, 20}$ điểm đo dọc theo chiều dài dầm.
  • Số lượng vết nứt giả thiết: Thiết lập mạng lưới vết nứt giả định $N \in {15, 20}$ vị trí.
  • Kịch bản kiểm tra độ sâu: Khảo sát các vết nứt có chiều sâu $a/h$ từ $10%$ đến $50%$. Trong kịch bản dầm công-xôn chứa 3 vết nứt phức tạp, các thông số thực được kiểm chứng tại $e_1 = 0.3$ ($a_1/h = 30%$), $e_2 = 0.6$ ($a_2/h = 20%$) và $e_3 = 0.9$ ($a_3/h = 50%$).

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Xác lập ranh giới áp dụng của các phương trình xấp xỉ tiệm cận: Kết quả tính toán số so sánh giữa phương trình chính xác (EX) và các xấp xỉ bậc một (A1), bậc hai (A2), bậc ba (A3) trên thanh chứa 9 vết nứt (ngàm - ngàm), 8 vết nứt (ngàm - tự do) và 4 vết nứt (tự do - tự do) cho thấy: Xấp xỉ bậc nhất chỉ đảm bảo độ chính xác khi chiều sâu vết nứt $a/h \le 20%$. Xấp xỉ bậc hai và bậc ba duy trì độ tin cậy vượt trội với chiều sâu vết nứt lên tới $30% - 40%$. Đặc biệt đối với mode dao động thứ nhất của thanh hai đầu tự do, xấp xỉ bậc ba cho kết quả tiệm cận trùng khít với nghiệm chính xác ngay cả khi chiều sâu vết nứt đạt tới $50%$.
Bậc xấp xỉ Phạm vi chiều sâu nứt ($a/h$) Sai số tương đối tối đa Ứng dụng khuyến nghị
Bậc nhất (A1) $< 20%$ $< 1.5%$ Vết nứt chân chim giai đoạn sớm
Bậc hai (A2) $20% - 35%$ $< 2.0%$ Vết nứt phát triển trung bình
Bậc ba (A3) $30% - 45%$ ($50%$ cho Mode 1) $< 1.8%$ Vết nứt nghiêm trọng / Đa khuyết tật
Chính xác (EX) Mọi chiều sâu ($10% - 50%$) $0.00%$ (Chuẩn đối chiếu) Chẩn đoán chuẩn mực & Đánh giá an toàn
  1. Khám phá quy luật phân bố điểm nút tần số (Frequency Nodes): Luận án đã chỉ ra rằng trên kết cấu tồn tại các vị trí không gian đặc biệt mà khi vết nứt xuất hiện tại đó, một tần số riêng cụ thể sẽ hoàn toàn không bị thay đổi ($d_1(\omega_0, x) = 0$). Ví dụ, đối với thanh hai đầu ngàm, mode 1 có điểm nút tại vị trí giữa nhịp $e = 0.5$; mode 2 có các điểm nút tại $e = 0.25$ và $e = 0.75$. Phát hiện này giải thích nguyên nhân các phương pháp chẩn đoán đơn tần số truyền thống thường bỏ sót khuyết tật nghiêm trọng nếu vết nứt vô tình rơi vào điểm nút.
gantt
    title Tọa độ điểm nút tần số riêng (Frequency Nodes) trên thanh
    dateFormat X
    axisFormat %s
    section Hai đầu ngàm
    Mode 1 (e = 0.50)           :active, 50, 51
    Mode 2 (e = 0.25; 0.75)      :active, 25, 26
    Mode 2 (e = 0.75)           :active, 75, 76
    section Hai đầu tự do
    Mode 1 (Không có điểm nút)  :crit, 0, 1
    Mode 2 (e = 0.50)           :active, 50, 51
    section Ngàm - Tự do
    Mode 1 (Không có điểm nút)  :crit, 0, 1
    Mode 2 (e = 0.33)           :active, 33, 34

  1. Tính ưu việt của tần số phản cộng hưởng trong nhận dạng đa vết nứt: Khi tăng số lượng điểm đo từ 10 lên 14, 16, 18 và 20 điểm trên dầm công-xôn, thuật toán chẩn đoán hội tụ chính xác tuyệt đối về tọa độ thực và độ sâu thực của các vết nứt ($e_1=0.3, a_1/h=30%; e_2=0.6, a_2/h=20%; e_3=0.9, a_3/h=50%$). Khi số điểm đo $M \ge 16$, các đỉnh giả định (ghost peaks) bị triệt tiêu hoàn toàn.

  2. Hiện tượng tương đương động học giữa phản cộng hưởng và dạng riêng: Luận án chứng minh bằng toán học rằng việc thu thập phổ phản cộng hưởng tại các vị trí phân bố thực chất cung cấp thông tin không gian tương đương với việc đo ma trận dạng dao động riêng đầy đủ nhưng với độ tin cậy tín hiệu vượt trội, loại bỏ hoàn toàn các lỗi tích lũy do đạo hàm số học dạng dao động.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp công cụ toán học hoàn chỉnh cho bài toán động lực học của hệ thanh/dầm gián đoạn, mở đường cho việc phát triển các mô hình giải tích cho hệ khung không gian, dầm nhiều nhịp và kết cấu tấm vỏ.
  • Về mặt phương pháp luận: Thiết lập quy trình chẩn đoán bài toán ngược vững chắc dựa trên sự kết hợp giữa phân tích giải tích chính xác, xấp xỉ tiệm cận Taylor, và thuật toán tối ưu hóa điều chỉnh Tikhonov.
  • Về mặt thực tiễn kỹ thuật: Cho phép kỹ sư kiểm định cầu đường và công trình biển phát hiện chính xác vị trí và độ sâu của nhiều vết nứt ẩn giấu bên trong kết cấu chỉ với một số lượng hữu hạn các cảm biến gia tốc thông thường kết hợp búa gõ kích động.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn học thuật:

  1. Giả thiết vết nứt hở không lan truyền: Mô hình giả định vết nứt luôn luôn mở (open crack) và không thay đổi trạng thái đóng mở trong quá trình dao động (breathing crack), do đó chưa phản ánh tính chất phi tuyến hình học khi kết cấu chịu tải trọng uốn đổi chiều biên độ lớn.
  2. Giới hạn mô hình dầm Euler-Bernoulli: Nghiên cứu chủ yếu áp dụng lý thuyết dầm mảnh Euler-Bernoulli bỏ qua biến dạng trượt và quán tính quay, điều này có thể dẫn đến sai số khi áp dụng cho dầm sâu (dầm Timoshenko) hoặc ở các bậc tần số dao động rất cao.
  3. Môi trường không cản trong mô hình giải tích: Để thiết lập phương trình đa thức hiển, ma trận cản được xem là nhỏ hoặc bị triệt tiêu, chưa khảo sát sâu ảnh hưởng của cản phi tỷ lệ trong môi trường thực tế.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng chính:

  • Mở rộng lý thuyết tần số phản cộng hưởng cho mô hình vết nứt thở phi tuyến (breathing cracks) sử dụng phương pháp cân bằng điều hòa.
  • Phát triển phương trình phản cộng hưởng cho dầm Timoshenko và kết cấu tấm 2D vật liệu biến tính chức năng (Functionally Graded Materials - FGM).
  • Tích hợp trí tuệ nhân tạo (Mạng nơ-ron tích chập 1D, Physics-Informed Neural Networks - PINNs) để tự động hóa quá trình giải bài toán ngược từ dữ liệu phổ phản cộng hưởng thời gian thực.
  • Thử nghiệm thực nghiệm quy mô phòng thí nghiệm và hiện trường trên dầm cầu bê tông cốt thép dự ứng lực sử dụng thiết bị đo rung laser Doppler (SLDV).

Tác động và ảnh hưởng

Nghiên cứu của tác giả Phạm Thị Ba Liên tạo ra tác động sâu rộng trên cả ba phương diện: Học thuật, Công nghiệp và Quản lý kết cấu:

  • Tác động học thuật: Định hình một hướng đi mới trong cộng đồng Chẩn đoán Kỹ thuật và Giám sát Sức khỏe Kết cấu (SHM) tại Việt Nam và quốc tế, mở rộng năng lực giải quyết bài toán ngược từ việc chỉ dùng tần số riêng sang trường dữ liệu phản cộng hưởng đa kênh.
  • Chuyển đổi công nghiệp: Ứng dụng trực tiếp trong quy trình kiểm định định kỳ các dầm cầu vượt giao thông, cọc móng công trình ngầm, cần cẩu tháp cảng biển và turbine gió. Phương pháp giúp chuyển đổi từ bảo trì thụ động sang bảo trì dự đoán (predictive maintenance), giảm thiểu rủi ro sự cố công trình lên tới $70%$.
  • Lợi ích kinh tế - xã hội: Giúp tiết kiệm hàng tỷ đồng chi phí kiểm định phá hủy hoặc dừng khai thác công trình không cần thiết, kéo dài tuổi thọ khai thác an toàn của cơ sở hạ tầng giao thông quốc gia.

Đối tượng hưởng lợi

mindmap
  root((Hệ sinh thái hưởng lợi))
    Nghiên cứu sinh & Giới học thuật
      Mô hình giải tích đa thức mới
      Hướng nghiên cứu FGM và Timoshenko
      Khung lý thuyết bài toán ngược
    Doanh nghiệp R&D & Kiểm định
      Phần mềm chẩn đoán hư hỏng tự động
      Quy trình đo đạc phản cộng hưởng
      Tối ưu hóa số lượng cảm biến
    Cơ quan quản lý & Bộ GTVT
      Quy chuẩn kỹ thuật kiểm định cầu
      Hệ thống SHM thời gian thực
      Đảm bảo an toàn hạ tầng quốc gia
  1. Nghiên cứu sinh và Giới học thuật Cơ học: Tiếp cận hệ phương trình tần số giải tích tường minh và công thức Rayleigh mở rộng để phát triển các đề tài chuyên sâu về dao động phi tuyến và cơ học phá hủy.
  2. Kỹ sư R&D và Các đơn vị Kiểm định Công trình (GTVT, Xây dựng): Sở hữu thuật toán số tin cậy đã kiểm chứng trên MATLAB để tích hợp vào các thiết bị phân tích dao động chuyên dụng và phần mềm chẩn đoán tự động.
  3. Cơ quan Quản lý Cơ sở Hạ tầng Giao thông: Có cơ sở khoa học vững chắc để ban hành các quy chuẩn, tiêu chuẩn kỹ thuật mới về giám sát sức khỏe kết cấu nhịp cầu và công trình cao tầng.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo và có tính đột phá nhất của luận án là gì?
Đóng góp độc đáo nhất là việc thiết lập thành công phương trình tần số dao động dọc trục của thanh có $n$ vết nứt dưới dạng đa thức tường minh bậc $n$ theo các tham số vết nứt $\gamma_j$ và mở rộng công thức Rayleigh cho hệ thanh có khuyết tật gián đoạn. Công trình đã chuyển hóa một bài toán ma trận siêu việt phức tạp từng tồn tại từ nghiên cứu của Ruotolo và Surace (2004) thành một công thức đại số tường minh dễ tính toán và phân tích độ nhạy.

2. Điểm mới về phương pháp luận so với các nghiên cứu quốc tế trước đây là gì?
Khác với Morassi (2001) chỉ dùng tần số riêng đơn thuần (dẫn đến nghiệm không duy nhất) hay Meruane và Heylen (2010) sử dụng thuật toán di truyền tính toán nặng nề, luận án đã tích hợp trực tiếp phương trình tần số phản cộng hưởng giải tích vào thuật toán điều chỉnh Ti-khô-nôv kết hợp khai triển giá trị kỳ dị (SVD). Điều này giúp nghiệm của bài toán ngược hội tụ đơn trị, ổn định và xác định chính xác đồng thời cả vị trí lẫn độ sâu của đa vết nứt mà không bị rơi vào cực trị địa phương.

3. Phát hiện nào trong dữ liệu mô phỏng số gây bất ngờ nhất?
Phát hiện bất ngờ nhất là tính chính xác vượt trội của xấp xỉ tiệm cận bậc ba (A3) của công thức Rayleigh: Đối với mode dao động thứ nhất của thanh tự do hai đầu, công thức gần đúng bậc ba có thể tính toán chính xác tần số riêng tương đương phương trình giải tích chính xác ngay cả khi vết nứt có chiều sâu lên tới $50%$ tiết diện thanh. Ngoài ra, việc bổ sung các điểm đo phản cộng hưởng giúp triệt tiêu hoàn toàn các đỉnh nứt giả (ghost cracks) khi số lượng cảm biến $M \ge 16$.

4. Luận án có cung cấp quy trình tái lập (replication protocol) đầy đủ không?
Luận án cung cấp đầy đủ và chi tiết toàn bộ các ma trận hệ số, điều kiện biên tổng quát ($\alpha_0, \beta_0, \alpha_1, \beta_1$), hàm đáp ứng tần số giải tích $H(\omega)$, sơ đồ khối thuật toán chẩn đoán và các tham số vật liệu mô phỏng trong môi trường MATLAB, cho phép bất kỳ nhà nghiên cứu nào cũng có thể tái lập hoàn toàn các kết quả số.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
Lộ trình 10 năm tập trung vào việc mở rộng hệ phương trình phản cộng hưởng sang hệ kết cấu 2D/3D phức tạp (tấm composite FGM, vỏ trụ), tích hợp cảm biến thông minh không dây (MEMS/áp điện Piezoelectric), phát triển thuật toán chẩn đoán thời gian thực dựa trên học sâu (Deep Learning) và thương mại hóa phần mềm SHM chuyên dụng cho hệ thống hạ tầng giao thông quốc gia.

Kết luận

Luận án tiến sĩ của tác giả Phạm Thị Ba Liên là một công trình khoa học mẫu mực, kết hợp hài hòa giữa toán học giải tích đỉnh cao và kỹ thuật tính toán hiện đại để giải quyết trọn vẹn bài toán chẩn đoán vết nứt trong kết cấu công trình.

Các thành tựu cốt lõi của công trình bao gồm:

  1. Thiết lập thành công hệ phương trình đặc trưng tần số cộng hưởng và phản cộng hưởng dạng đa thức tường minh cho kết cấu thanh (dao động dọc trục) và dầm (dao động uốn) chứa nhiều vết nứt dưới điều kiện biên tổng quát.
  2. Xây dựng và chứng minh công thức Rayleigh mở rộng cho thanh có nhiều vết nứt, đồng thời định lượng hóa chính xác phạm vi tin cậy của các nghiệm xấp xỉ tiệm cận bậc nhất ($a/h < 20%$), bậc hai và bậc ba ($a/h$ từ $30% - 50%$).
  3. Khám phá bản chất cơ học và quy luật toán học của các điểm nút tần số (frequency nodes), làm sáng tỏ nguyên nhân bỏ sót khuyết tật của các phương pháp dao động cổ điển.
  4. Phát triển thành công thuật toán chẩn đoán ngược đa vết nứt ổn định bằng cách kết hợp dữ liệu tần số phản cộng hưởng phân bố với phương pháp điều chỉnh Ti-khô-nôv và khai triển kỳ dị SVD.
  5. Thử nghiệm số học toàn diện trên môi trường MATLAB chứng minh khả năng nhận dạng chính xác đồng thời vị trí và chiều sâu của hệ thống đa vết nứt phức tạp (lên đến 19-20 vết nứt).
  6. Mở ra nhánh nghiên cứu đột phá về khai thác đặc trưng động lực học phản cộng hưởng trong lĩnh vực Giám sát Sức khỏe Kết cấu (SHM), đặt nền tảng vững chắc cho các ứng dụng kiểm định thông minh trong kỷ nguyên Công nghiệp 4.0.