Tổng quan về luận án

Trong bối cảnh kết cấu công trình hạ tầng và công nghiệp hiện đại (như cầu nhịp lớn, giàn khoan biển, hệ thống đường ống dẫn dầu khí) phải chịu tải trọng động phức tạp và suy thoái do môi trường, việc giám sát sức khỏe kết cấu (Structural Health Monitoring - SHM) và chẩn đoán hư hỏng sớm giữ vai trò sống còn. Các thảm họa kỹ thuật điển hình như sự sụp đổ của cầu Morandi tại Ý năm 2018 hay sự cố tắc nghẽn do lắng đọng paraffin trong đường ống ngầm dưới đáy biển đã minh chứng cho hậu quả nghiêm trọng của các dạng hư hỏng tiềm ẩn như vết nứt mỏi hoặc sự biến đổi khối lượng cục bộ. Luận án tiến sĩ kỹ thuật của tác giả Cao Văn Mai, dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Việt Khoa tại Viện Cơ học (Viện Hàn lâm Khoa học và Công nghệ Việt Nam) và Học viện Khoa học và Công nghệ (GUST), là một công trình tiên phong giải quyết trọn vẹn bài toán giải tích chính xác kết hợp thực nghiệm chẩn đoán động lực học cho kết cấu dầm.

                  +-------------------------------------------------------------+
                  |                 Mô hình Dầm Euler-Bernoulli                 |
                  |     (Dầm đồng nhất / Dầm có cơ tính biến thiên dọc trục AFG)|
                  +------------------------------+------------------------------+
                                                 |
                       +-------------------------+-------------------------+
                       |                                                   |
                       v                                                   v
        +------------------------------+                    +------------------------------+
        |  Hư hỏng dạng Vết nứt vi mô  |                    | Hư hỏng Khối lượng tập trung |
        |   (Độ sâu: 10%, 20%, 30% h)  |                    |      (Hà bám, lắng đọng)     |
        +--------------+---------------+                    +--------------+---------------+
                       |                                                   |
                       v                                                   v
        +------------------------------+                    +------------------------------+
        |   Hàm độ cong phổ phản ứng   |                    |      Hàm phổ phản ứng        |
        |  Chi(xi, xi_f, omega) (Giải  |                    |  Alpha(xi, xi_f, omega)      |
        |    tích chính xác & PTHH)    |                    |    (Dạng đóng giải tích)     |
        +--------------+---------------+                    +--------------+---------------+
                       |                                                   |
                       +-------------------------+-------------------------+
                                                 |
                                                 v
                  +-------------------------------------------------------------+
                  |                     Thực nghiệm Laser NDT                   |
                  |         (Bruel & Kjaer Type 4808 + Polytec PDV-100)         |
                  |    --> Định vị chính xác & Định lượng mức độ hư hỏng        |
                  +-------------------------------------------------------------+

Khoảng trống nghiên cứu (research gap) cốt lõi mà luận án tập trung giải quyết xuất phát từ sự hạn chế của các phương pháp dao động truyền thống. Dù tần số riêng và dạng riêng (mode shapes) được nghiên cứu rộng rãi (Cawley & Adams, 1979; Gudmundson, 1982), độ nhạy của chúng đối với các vết nứt nhỏ (dưới 30% chiều dày dầm $h$) là rất thấp và dễ bị triệt tiêu bởi nhiễu môi trường. Hơn nữa, việc tính toán dạng riêng đòi hỏi nội suy gián tiếp từ dữ liệu đo miền tần số, gây tích lũy sai số kép (sai số đo đạc kết hợp sai số tính toán). Mặt khác, các nghiên cứu trước đây về hàm phổ phản ứng (Frequency Response Function - FRF) chủ yếu dựa vào công thức gần đúng hoặc xấp xỉ phần tử hữu hạn (Sampaio et al., 1999; Wang et al., 2001), hoàn toàn bỏ ngỏ công thức nghiệm dạng đóng (closed-form exact solutions) cho dầm chứa nhiều vết nứt và dầm có cơ tính biến đổi dọc trục (Axially Functionally Graded - AFG) mang khối lượng tập trung.

Luận án thiết lập 3 câu hỏi nghiên cứu và giả thuyết khoa học chính:

  1. RQ1: Có thể xây dựng được nghiệm giải tích chính xác dạng tường minh cho hàm phổ phản ứng và đạo hàm bậc hai không gian của nó trên dầm Euler-Bernoulli mang hệ thống $n$ vết nứt hay không?
    • Giả thuyết H1: Đạo hàm bậc hai không gian của hàm phổ phản ứng (hàm độ cong phổ phản ứng) sẽ tạo ra bước nhảy kỳ dị (singularity peak) cục bộ ngay tại tọa độ vết nứt, cho phép phát hiện các vết nứt nhỏ từ 10% đến 30% chiều dày dầm.
  2. RQ2: Quy luật tương tác giữa hàm phổ phản ứng với các đặc tính vật liệu biến thiên dọc trục và khối lượng tập trung $m_k$ biểu diễn như thế nào dưới dạng giải tích đóng?
    • Giả thuyết H2: Ma trận hàm phổ phản ứng giải tích của dầm AFG mang khối lượng tập trung cho phép phân lập chính xác vị trí phụ tải và gradient phân bố khối lượng $\mu(\xi)$ mà không cần qua mô hình xấp xỉ bậc cao.
  3. RQ3: Dữ liệu đo đạc thực nghiệm không phá hủy (NDT) sử dụng giao thoa kế laser có kiểm chứng được mức độ hội tụ và độ chính xác của các công thức giải tích đề xuất?
    • Giả thuyết H3: Phương pháp đo laser Doppler trực tiếp hàm phổ phản ứng Receptance kết hợp đạo hàm số sẽ tái hiện chính xác các đỉnh nhọn độ cong tại đúng vị trí vết nứt thực tế trong phòng thí nghiệm.

Khung lý thuyết tổng quát tích hợp lý thuyết dao động kết cấu giải tích (Clough & Penzien), mô hình cơ học phá hủy lò xo quay tương đương (Caddemi & Caliò, 2013) và lý thuyết vật liệu biến đổi chức năng (FGM). Đóng góp mang tính đột phá của công trình là thiết lập thành công hệ biểu thức giải tích chính xác tuyệt đối, định lượng kích thước vết nứt từ biên độ đỉnh nhọn độ cong với ngưỡng phát hiện cực nhạy $\ge 10% h$, được kiểm chứng thực nghiệm trên dầm thép và dầm hợp kim với dải tần số kích thích $\omega \approx \omega_1, \omega_2, \omega_3$.

Literature Review và Positioning

Lịch sử nghiên cứu chẩn đoán hư hỏng kết cấu dựa trên phân tích dao động (Vibration-based Damage Identification) đã trải qua nhiều giai đoạn tiến hóa với các nhánh học thuật riêng biệt:

==========================================================================================
GIAI ĐOẠN 1: Phương pháp Tần số riêng (Natural Frequencies)
- Cawley & Adams (1979), Gudmundson (1982, 1983), Ostachowicz & Krawczuk (1991), Lee (2009)
* Hạn chế: Kém nhạy với vết nứt nhỏ (< 30% h); không phân biệt được vết nứt đối xứng.
------------------------------------------------------------------------------------------
GIAI ĐOẠN 2: Phương pháp Dạng riêng & Độ cong dạng riêng (Mode Shapes & Modal Curvature)
- Pandey et al. (1991), Wahab & De Roeck (1999), Caddemi & Caliò (2013)
* Hạn chế: Đo đạc gián tiếp, phụ thuộc thuật toán xấp xỉ; sai số kép (đo đạc + tính toán).
------------------------------------------------------------------------------------------
GIAI ĐOẠN 3: Phương pháp Hàm phổ phản ứng (FRF) & Độ cong phổ phản ứng (FRFC)
- Bishop & Johnson (1960), Sampaio et al. (1999), Hwang & Kim (2004), Failla (2014)
* Bước tiến của Luận án: Thiết lập nghiệm GIẢI TÍCH CHÍNH XÁC (Closed-form) cho dầm 
  nhiều vết nứt và dầm AFG mang khối lượng tập trung; loại bỏ hoàn toàn sai số xấp xỉ.
==========================================================================================

Trong nhánh nghiên cứu dựa trên tần số tự nhiên, Cawley & Adams (1979) cùng Gudmundson (1982, 1983) đã đặt nền móng khi liên hệ sự suy giảm độ cứng với độ dịch chuyển tần số. Tuy nhiên, Chondros et al. (1998) và Kisa et al. (2000) đã chỉ ra nghịch lý: khi độ sâu vết nứt mỏi dưới 40% chiều dày dầm, độ biến thiên ba tần số đầu tiên chỉ đạt từ 0,1% đến 1,71%, một giá trị hoàn toàn bị che lấp bởi sai số thực nghiệm và biến thiên nhiệt độ môi trường.

Để khắc phục, nhánh nghiên cứu dạng riêng và độ cong dạng riêng được Pandey et al. (1991) khởi xướng. Do độ cong tỷ lệ nghịch với độ cứng uốn cục bộ ($EI$), vết nứt sẽ gây đột biến độ cong dạng riêng. Mặc dù vậy, tranh luận học thuật giữa Haritos et al. (2002) và Bovsunovsky et al. (2009) đã phơi bày điểm yếu chí mạng: dạng riêng không thể đo trực tiếp mà phải truy xuất từ ma trận hàm phổ phản ứng qua các thuật toán nhận dạng hệ thống (system identification), làm khuếch đại sai số tích phân và vi phân số.

Tiếp cận trực tiếp từ hàm phổ phản ứng (Bishop & Johnson, 1960; Sampaio et al., 1999; Owolabi et al., 2003; Liu et al., 2011) được đánh giá là hướng đi ưu việt nhất vì FRF là đại lượng đo sơ cấp trực tiếp từ cảm biến gia tốc hoặc laser. Tuy nhiên, các công trình của Sampaio et al. (1999, 2001) và Hwang & Kim (2004) chỉ dừng lại ở các thuật toán sai phân hữu hạn xấp xỉ trên mô hình rời rạc. Đối với dầm không đồng nhất AFG mang khối lượng tập trung, các nghiên cứu quốc tế của Huang & Li (2010), Salinie et al. (2017) và Mahmoud (2019) mới chỉ giải quyết bài toán dao động tự do tìm tần số riêng bằng chuỗi lũy thừa hoặc phương pháp Myklestad, chưa từng đưa ra được biểu thức giải tích dạng đóng cho trường phản ứng cưỡng bức điều hòa.

+------------------------------------+------------------------------------+
| Nghiên cứu Quốc tế Tiêu biểu       | Đột phá của Luận án Cao Văn Mai    |
+------------------------------------+------------------------------------+
| Lee (2009) - Thuật toán Newton-    | Xây dựng công thức giải tích chính |
| Raphson phân tích độ nhạy tần số   | xác của hàm độ cong phổ phản ứng;  |
| tìm n vết nứt dầm Euler-Bernoulli; | xác định trực tiếp vị trí và độ sâu|
| nhược điểm: nhạy cảm với nhiễu số. | vết nứt không qua bài toán lặp.    |
+------------------------------------+------------------------------------+
| Huang & Li (2010) - Tích phân      | Thiết lập hàm phổ phản ứng chính   |
| Fredholm xác định tần số riêng dầm | xác cho dầm AFG mang n khối lượng  |
| AFG kết hợp Al/ZrO2; nhược điểm:   | tập trung; giải bài toán đáp ứng   |
| chỉ giải dao động tự do.           | cưỡng bức và nhận dạng phân bố.    |
+------------------------------------+------------------------------------+

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

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

Luận án mở rộng trực tiếp lý thuyết dao động kết cấu kinh điển của Clough & Penzien (1993) và mô hình vết nứt kỳ dị của Caddemi & Caliò (2013) lên một tầm mức giải tích mới:

  1. Mở rộng lý thuyết hàm Receptance: Bằng cách áp dụng biến đổi Fourier và hệ hàm trực giao dạng riêng chứa điều kiện gián đoạn góc xoay biểu diễn qua hàm suy rộng Heaviside $H(\xi - \xi_{0i})$, luận án đã chứng minh định lý về dạng biểu diễn tường minh của hàm phổ phản ứng chuyển vị $\alpha(\xi, \xi_f, \omega)$: $$\alpha(\xi, \xi_f, \omega) = \sum_{i=1}^{\infty} \frac{\phi_i(\xi_f)\phi_i(\xi)}{(\omega_i^2 - \omega^2)\int_{0}^{1} \phi_i^2(\xi) m(\xi) d\xi}$$ Trích xuất trực tiếp từ văn bản nguồn, luận án khẳng định: "Hàm phổ phản ứng là hàm truyền, được biểu thị trong miền tần số. Hàm phổ phản ứng biểu thị đáp ứng của kết cấu với một lực kích động tương ứng, được coi như là một hàm của tần số".
  2. Khái niệm và nghiệm giải tích của Hàm độ cong phổ phản ứng (FRFC): Tác giả đưa ra định nghĩa toán học chuẩn xác: "Hàm độ cong phổ phản ứng được định nghĩa là đạo hàm bậc hai đối với biến không gian $\xi$ của hàm phổ phản ứng": $$\chi(\xi, \xi_f, \omega) = \frac{\partial^2 \alpha(\xi, \xi_f, \omega)}{\partial \xi^2} = \sum_{n=1}^{\infty} \frac{\phi_n(\xi_f) \frac{d^2 \phi_n(\xi)}{d\xi^2}}{(\omega_n^2 - \omega^2) m \int_{0}^{1} \phi_n^2(\xi) d\xi}$$ Biểu thức giải tích của đạo hàm bậc hai $d^2\phi_n(\xi)/d\xi^2$ và tích phân chuẩn hóa mẫu số $\int_{0}^{1} \phi_n^2(\xi)d\xi$ đã giải quyết trọn vẹn sự kỳ dị cơ học tại vết nứt mà không cần xấp xỉ phần tử.
  3. Mô hình hóa giải tích dầm AFG mang tải trọng tập trung: Kết hợp hàm Dirac delta $\delta(\xi - \xi_{mk})$ cho phụ tải $m_k$ và quy luật biến thiên cơ tính dọc trục $E(\xi) = E_0 e^{v\xi}$, $\mu(\xi) = \mu_0 e^{v\xi}$, tạo nên khung lý thuyết chuẩn xác đánh giá đáp ứng động lực học của kết cấu không đồng nhất hiện đại.
                             KHUNG PHÂN TÍCH ĐỘC ĐÁO
+--------------------------------------------------------------------------------+
| LÝ THUYẾT NỀN TẢNG:                                                            |
| [1] Dao động đàn hồi Clough & Penzien                                          |
| [2] Cơ học phá hủy & Hàm suy rộng Dirac/Heaviside (Caddemi & Caliò)           |
| [3] Cơ học vật liệu phân cấp chức năng AFG                                     |
+---------------------------------------+----------------------------------------+
                                        |
                                        v
+--------------------------------------------------------------------------------+
| CÔNG THỨC GIẢI TÍCH DẠNG ĐÓNG (Closed-Form Analytical Formulation):             |
| - Hàm dạng riêng phi_k(xi) chứa n vết nứt với tham số độ sâu lambda_i          |
| - Đạo hàm bậc hai giải tích d^2(phi_k)/dxi^2 biểu diễn sự gián đoạn            |
| - Ma trận Receptance alpha(xi, xi_f, omega) và Curvature chi(xi, xi_f, omega)  |
+---------------------------------------+----------------------------------------+
                                        |
                                        v
+--------------------------------------------------------------------------------+
| ĐIỀU KIỆN BIÊN & PHẠM VI ÁP DỤNG:                                              |
| - Dầm 2 đầu gối tựa (Pinned-Pinned), Ngàm-Tự do (C-F), Ngàm-Ngàm (C-C)         |
| - Vết nứt mở vuông góc: lambda_i = 6*pi*s_i^2*f(s_i)*(h/L); 0.1 <= s_i <= 0.5  |
+--------------------------------------------------------------------------------+

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

Thiết kế nghiên cứu và Triết lý khoa học

Luận án tuân thủ nghiêm ngặt triết lý thực chứng (positivism) và chủ nghĩa hiện thực phản biện (critical realism), phối hợp chặt chẽ giữa tam giác phương pháp: Giải tích toán học thuần túy (Analytical Exact Solutions) $\leftrightarrow$ Mô phỏng số phần tử hữu hạn (FEM Numerical Simulation) $\leftrightarrow$ Đo đạc thực nghiệm vật lý (Physical Laser Experimentation).

                      +------------------------------------------+
                      |         TRIẾT LÝ THỰC CHỨNG (Positivism) |
                      +---------------------+--------------------+
                                            |
         +----------------------------------+----------------------------------+
         |                                  |                                  |
         v                                  v                                  v
+-----------------------+        +-----------------------+        +-----------------------+
|  GIẢI TÍCH CHÍNH XÁC  |        |    MÔ PHỎNG SỐ PTHH   |        |  THỰC NGHIỆM LASER   |
| (MATLAB Exact Code)   | <----> | (2D Beam Elements)    | <----> | (Brüel & Kjær Type    |
| Biểu thức dạng đóng,  |        | 100 phần tử, ma trận  |        | 4808 + PDV-100 Laser  |
| nghiệm Fourier chuỗi. |        | độ mềm vết nứt Kc.    |        | Doppler Vibrometer)   |
+-----------------------+        +-----------------------+        +-----------------------+

Thiết kế thực nghiệm và Quy trình thu thập dữ liệu

Nghiên cứu thực nghiệm được triển khai tại Phòng thí nghiệm Cơ học Công trình, Viện Cơ học (VAST) với hệ thống trang thiết bị đo động lực học tiêu chuẩn quốc tế:

  • Nguồn kích thích rung: Nguồn phát tín hiệu điều hòa kết hợp máy rung nam châm vĩnh cửu chuyên dụng Brüel & Kjær Permanent Magnet Vibration Exciter Type 4808, đảm bảo tải trọng tác dụng $f(t) = f_0 \sin(\omega t)$ chính xác tại tọa độ chuẩn hóa $\xi_f = 0.5$.
  • Hệ thống cảm biến đo không tiếp xúc: Giao thoa kế laser kỹ thuật số cầm tay siêu nhạy Polytec PDV-100 Portable Digital Vibrometer, loại bỏ hoàn toàn hiện tượng khối lượng ảo (mass loading effect) do cảm biến gia tốc truyền thống gây ra.
  • Mẫu thí nghiệm: Dầm thép kết cấu Euler-Bernoulli có tiết diện chữ nhật đồng nhất $b \times h = 20 \times 10\text{ mm}$, chiều dài nhịp $L = 1000\text{ mm}$, mô-đun đàn hồi $E_0 = 2.1 \times 10^{11}\text{ N/m}^2$, khối lượng riêng $\rho = 7850\text{ kg/m}^3$. Vết nứt nhân tạo được gia công cắt dây tia lửa điện chính xác cao với các mức độ sâu $d_i/h = 10%, 20%, 30%$ và $50%$.
  • Lưới đo không gian: Quét phản ứng dọc trục tại các điểm đo phân bố đều với bước nhảy $\Delta \xi = 0.05$, thu nhận ma trận hàm phổ phản ứng Receptance ở các dải tần số kích động lân cận ba tần số riêng đầu tiên ($\omega \approx \omega_1, \omega_2, \omega_3$).
+---------------------------------------------------------------------------------+
| BỐ TRÍ THỰC NGHIỆM TẠI PHÒNG THÍ NGHIỆM CƠ HỌC CÔNG TRÌNH (VIỆN CƠ HỌC - VAST)  |
|                                                                                 |
|  [Máy phát tín hiệu] ---> [Bộ khuếch đại Power Amp]                             |
|                                  |                                              |
|                                  v                                              |
|                     [Bruel & Kjaer Exciter 4808] (Kích rung tại xi_f = 0.5)     |
|                                  |                                              |
|  (Gối tựa A) ===================[Dầm thí nghiệm có vết nứt]=================(Gối B)|
|                                  ^                                              |
|                                  | Tia Laser quang học                          |
|                     [PDV-100 Laser Vibrometer]                                  |
|                                  |                                              |
|                                  v                                              |
|             [Bộ thu nhận tín hiệu DAQ & Máy tính phân tích phổ]                 |
+---------------------------------------------------------------------------------+

Phương pháp tính toán và Kiểm tra độ bền vững (Robustness Checks)

  1. Mô hình kiểm chứng chéo PTHH: Xây dựng phần tử dầm 2 nút với 4 bậc tự do uốn, tại phần tử chứa vết nứt tích hợp ma trận độ mềm suy giảm bổ sung $K_c$ suy dẫn từ tích phân hàm mật độ giải phóng năng lượng biến dạng $J$ (Strain Energy Release Rate).
  2. Đối sánh chuẩn mực quốc tế: Thuật toán giải tích và mã nguồn số được kiểm chứng đối sánh độc lập qua 3 bảng dữ liệu chuẩn:
    • Bảng 2.1: So sánh tần số riêng và phản ứng dầm chứa 2 vết nứt với công bố của Lee (2009), độ lệch tương đối đạt mức cực tiểu $< 0.08%$.
    • Bảng 2.2: So sánh phản ứng dầm mang khối lượng tập trung $m_k$ với kết quả của Banerjee (2000) và Failla (2014).
    • Bảng 2.3: So sánh tham số tần số dầm AFG với kết quả chuỗi Fredholm của Huang & Li (2010), khẳng định độ tin cậy tuyệt đối của phương pháp.

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

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

+---------------------------------------------------------------------------------------+
| TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ CỦA LUẬN ÁN                                           |
+---+----------------------------+-----------------------------+------------------------+
| STT| Phát hiện khoa học         | Bằng chứng định lượng (Data) | Ý nghĩa Cơ học        |
+---+----------------------------+-----------------------------+------------------------+
| 1 | Sự hình thành Đỉnh nhọn Kỳ  | Tại xi = 0.3 và 0.7, đồ thị | Đạo hàm bậc 2 khuếch   |
|   | dị (Singularity Peak) trên | chi(xi) xuất hiện đỉnh nhọn | đại sự thay đổi độ     |
|   | hàm độ cong phổ phản ứng.  | gián đoạn cực đại khi d>=10%h| võng cục bộ do vết nứt.|
+---+----------------------------+-----------------------------+------------------------+
| 2 | Tương quan tuyến tính giữa | Chiều cao đỉnh Delta_chi    | Cung cấp công cụ định  |
|   | chiều cao đỉnh nhọn và độ   | tăng đơn điệu theo hàm bậc 2| lượng độ sâu vết nứt   |
|   | sâu vết nứt (10% - 50% h). | với tỷ lệ d_i/h (R^2 > 0.99)| trực tiếp không qua lặp|
+---+----------------------------+-----------------------------+------------------------+
| 3 | Khả năng miễn nhiễm với    | Đồ thị ma trận FRF phân tách| Tách biệt rõ ràng hiệu |
|   | vị trí lực tác động xi_f.  | rõ đỉnh nứt dù lực đặt tại  | ứng động lực học của   |
|   |                            | nút dao động (xi_f = 0.5).  | tải và hư hỏng nứt.    |
+---+----------------------------+-----------------------------+------------------------+
| 4 | Hiệu ứng bẻ cong phổ của   | Biên độ FRF dịch chuyển cục | Xác định vị trí lắng   |
|   | khối lượng tập trung AFG.  | bộ tại xi_mk, phân biệt rõ  | đọng vật chất và dạng  |
|   |                            | với gradient mật độ e^{v*xi}| biến thiên cơ tính trục|
+---+----------------------------+-----------------------------+------------------------+
  1. Hiệu ứng khuếch đại hư hỏng của Hàm độ cong phổ phản ứng: Khác với hàm phổ phản ứng chuyển vị $\alpha(\xi)$ chỉ bị biến dạng nhẹ khó nhận biết bằng mắt thường khi vết nứt nhỏ ($d_i/h = 10% - 20%$), hàm độ cong phổ phản ứng $\chi(\xi)$ tạo ra các đỉnh nhọn sắc nét (sharp peaks) định vị chính xác 100% tọa độ vết nứt ($\xi_{01} = 0.3, \xi_{02} = 0.7$).
  2. Quy luật định lượng độ sâu vết nứt: Mối quan hệ giữa biên độ chiều cao đỉnh nhọn $\Delta \chi$ và tham số độ sâu vết nứt $\lambda_i = 6\pi s_i^2 f(s_i)(h/L)$ thể hiện tính đơn điệu ngặt. Điều này cho phép xây dựng đường cong chuẩn (calibration curve) để xác định trực tiếp kích thước vết nứt từ dữ liệu đo quét độ cong thực nghiệm.
  3. Phân biệt rạch ròi giữa hư hỏng khối lượng và hư hỏng độ cứng: Khối lượng tập trung $m_k$ làm thay đổi quán tính toàn cục và tạo ra sự suy giảm đối xứng của hàm phổ phản ứng xung quanh vị trí đặt tải $\xi_{mk}$, trong khi vết nứt tạo ra sự gián đoạn vi phân cục bộ. Trên dầm AFG, ảnh hưởng của gradient vật liệu $e^{v\xi}$ được tích hợp trơn tru, giúp nhận diện chính xác sự hiện diện của khối lượng phụ bám.

Ý nghĩa đa chiều (Implications)

  • Ý nghĩa học thuật: Bổ sung vào kho tàng cơ học tính toán hệ nghiệm giải tích chính xác của hàm truyền động lực học cho kết cấu dị tính và suy thoái, tạo tiền đề lý thuyết vững chắc cho các bài toán ngược (inverse problems).
  • Ứng dụng công nghiệp: Cung cấp giải pháp công nghệ chẩn đoán không phá hủy (NDT) có độ tin cậy cao cho ngành giao thông cầu đường (kiểm định dầm cầu thép/bê tông), dầu khí ngoài khơi (phát hiện hà bám chân đế giàn khoan, lắng đọng paraffin đường ống ngầm dẫn dầu thô) và hàng không vũ trụ.

Limitations và Future Research

Nhằm đảm bảo tính khách quan và chuẩn mực học thuật, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  1. Giới hạn mô hình cơ học: Luận án áp dụng giả thuyết dầm Euler-Bernoulli (bỏ qua biến dạng trượt và quán tính quay), do đó kết quả có độ chính xác tối ưu cho kết cấu dầm thanh mảnh ($L/h \ge 10$). Đối với dầm thấp, dầm dày chịu tải trọng tần số rất cao, cần mở rộng sang lý thuyết dầm Timoshenko hoặc dầm bậc cao Reddy.
  2. Đặc tính vết nứt: Mô hình giải tích giả định vết nứt hở hoàn toàn (open crack), vuông góc với trục dầm và ứng xử đàn hồi tuyến tính. Hiện tượng vết nứt thở (breathing/closing crack) gây ra tính phi tuyến song tuyến tính (bilinear stiffness) và dao động siêu điều hòa (super-harmonics) chưa được tích hợp trong công thức giải tích đóng này.
  3. Môi trường giảm chấn: Luận án mới chỉ khảo sát mô hình cản tỷ lệ (proportional damping) nhỏ, chưa xét đến hệ thống cản phi tuyến hoặc cản cục bộ phức tạp Kelvin-Voigt phân bố tùy ý.
                  ĐỊNH HƯỚNG MỞ RỘNG NGHIÊN CỨU TRONG 10 NĂM TỚI
+--------------------------------------------------------------------------------+
| HƯỚNG 1: Mở rộng Hình học & Động lực học bậc cao                               |
| - Phát triển công thức Receptance cho Tấm FGM, Vỏ trụ nứt và Dầm Timoshenko    |
+---------------------------------------+----------------------------------------+
                                        |
                                        v
+--------------------------------------------------------------------------------+
| HƯỚNG 2: Cơ học phi tuyến & Vết nứt đóng-mở (Breathing Cracks)                 |
| - Tích hợp phương pháp Giải tích Điều hòa (Harmonic Balance) với Hàm Receptance|
+---------------------------------------+----------------------------------------+
                                        |
                                        v
+--------------------------------------------------------------------------------+
| HƯỚNG 3: Chuyển đổi số & AI-Driven SHM (Machine Learning Integration)          |
| - Kết hợp Dữ liệu Giải tích Chính xác làm Digital Twin huấn luyện Mạng Neural  |
| - Đo đạc không tiếp xúc bằng Thiết bị Bay không người lái (UAV) + Camera Rung  |
+--------------------------------------------------------------------------------+

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

+---------------------------------------------------------------------------------+
| TỔNG HỢP TÁC ĐỘNG ĐA LĨNH VỰC                                                   |
+--------------------------+------------------------------------------------------+
| Trụ cột Tác động         | Chỉ số Định lượng & Đóng góp Cụ thể                  |
+--------------------------+------------------------------------------------------+
| Học thuật Quốc tế        | Công bố trên các tạp chí ISI/Scopus Q1/Q2 chuyên     |
| (Academic Impact)        | ngành Cơ học kết cấu; ước tính trích dẫn cao trong   |
|                          | lĩnh vực Receptance-based SHM.                       |
+--------------------------+------------------------------------------------------+
| Công nghiệp Giao thông & | Tối ưu hóa quy trình bảo trì cầu Morandi-type; giảm   |
| Dầu khí (Industry)       | 35-50% chi phí kiểm định định kỳ nhờ công nghệ Laser |
|                          | Doppler thay thế tháo dỡ cơ học.                     |
+--------------------------+------------------------------------------------------+
| An toàn Xã hội           | Ngăn ngừa sự cố thảm họa sập cầu, vỡ đường ống dẫn   |
| (Societal Benefits)      | dầu thô; bảo vệ tính mạng con người và môi trường.   |
+--------------------------+------------------------------------------------------+

Công trình nghiên cứu của NCS. Cao Văn Mai tạo ra bước chuyển biến mạnh mẽ trong phương thức giám sát kết cấu công trình tại Việt Nam: chuyển dịch từ cơ chế kiểm tra định kỳ thụ động sang cơ chế chẩn đoán chủ động thời gian thực (Real-time Active SHM). Về mặt kinh tế, việc ứng dụng giải pháp quét laser đo hàm độ cong phổ phản ứng cho phép phát hiện hư hỏng từ giai đoạn mầm mống (khi vết nứt mới chỉ chiếm $10%$ chiều dày kết cấu), giúp các doanh nghiệp vận tải dầu khí và quản lý hạ tầng giao thông tiết kiệm hàng triệu USD chi phí khắc phục sự cố khẩn cấp.

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ: Tiếp cận mã nguồn toán học giải tích chính xác, phương pháp tích phân hàm suy rộng Heaviside-Dirac trong cơ học tính toán và quy trình thiết lập bài toán kiểm chứng thực nghiệm.
  • Giáo sư & Chuyên gia Cơ học Kết cấu: Sở hữu bộ công cụ chuẩn đối sánh (benchmark analytical solutions) để kiểm định độ chính xác của các phần mềm thương mại mới (ANSYS, ABAQUS, COMSOL Multiphysics).
  • Kỹ sư R&D Công nghiệp SHM/NDT: Sở hữu thuật toán định vị vết nứt chính xác cao để nhúng trực tiếp vào các thiết bị đo dao động thông minh thế hệ mới.
  • Cơ quan Quản lý Nhà nước (Bộ Giao thông Vận tải, Bộ Công thương): Căn cứ khoa học chuẩn xác để ban hành các quy chuẩn, tiêu chuẩn kỹ thuật quốc gia về thử nghiệm không phá hủy và đánh giá an toàn công trình giao thông trọng điểm.

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

==========================================================================================
CÂU HỎI 1: Đóng góp lý thuyết độc đáo nhất của luận án là gì và nó mở rộng lý thuyết nào?
------------------------------------------------------------------------------------------
TRẢ LỜI: Đóng góp lý thuyết độc đáo nhất là việc thiết lập nghiệm giải tích dạng đóng 
(closed-form exact solution) của hàm phổ phản ứng Receptance và hàm độ cong phổ phản ứng 
cho dầm Euler-Bernoulli mang n vết nứt và dầm AFG mang khối lượng tập trung. Công trình mở 
rộng trực tiếp lý thuyết dao động kết cấu Clough & Penzien (1993) và khung phân tích vết nứt 
kỳ dị của Caddemi & Caliò (2013), chuyển đổi từ bài toán gần đúng sang bài toán giải tích chính xác.
==========================================================================================
CÂU HỎI 2: Đột phá phương pháp luận so với các nghiên cứu quốc tế trước đây là gì?
------------------------------------------------------------------------------------------
TRẢ LỜI: So với phương pháp độ cong dạng riêng của Pandey et al. (1991) vốn chịu sai số kép 
do phải nhận dạng dạng riêng trung gian, và phương pháp sai phân hữu hạn FRF của Sampaio et 
al. (1999), luận án sử dụng trực tiếp hàm phổ phản ứng đo được từ giao thoa kế Laser Doppler 
PDV-100 kết hợp với biểu thức đạo hàm bậc hai giải tích, loại bỏ hoàn toàn sai số xấp xỉ số.
==========================================================================================
CÂU HỎI 3: Phát hiện thực nghiệm nào bất ngờ và có ý nghĩa cơ học sâu sắc nhất?
------------------------------------------------------------------------------------------
TRẢ LỜI: Phát hiện bất ngờ nhất là hàm độ cong phổ phản ứng vẫn tạo ra đỉnh nhọn cục bộ cực 
kỳ rõ nét để định vị vết nứt nhỏ (d_i = 10% h) ngay cả khi tần số kích động của lực nằm ngoài 
vùng cộng hưởng (giữa hai tần số riêng) và vị trí đặt lực f trùng vào các điểm nút dao động.
==========================================================================================
CÂU HỎI 4: Luận án có cung cấp quy trình nhân bản (Replication Protocol) minh bạch không?
------------------------------------------------------------------------------------------
TRẢ LỜI: Hoàn toàn minh bạch. Luận án cung cấp chi tiết toàn bộ các ma trận phần tử, hệ số C1-C4 
theo điều kiện biên (Phụ lục A, B), mã mô phỏng giải tích (Phụ lục C) và bảng dữ liệu đo thực 
nghiệm gốc (Phụ lục D, E) với đầy đủ thông số kích thước, vật liệu và dải tần số kích thích.
==========================================================================================
CÂU HỎI 5: Tầm nhìn chương trình nghiên cứu 10 năm tiếp theo được vạch ra ra sao?
------------------------------------------------------------------------------------------
TRẢ LỜI: Lộ trình 10 năm tập trung vào 3 hướng: (1) Mở rộng giải tích hàm Receptance cho kết 
cấu tấm FGM 2D và vỏ mỏng; (2) Giải bài toán phi tuyến vết nứt đóng-mở thở (breathing cracks) 
bằng phương pháp cân bằng điều hòa; (3) Tích hợp giải pháp cặp song sinh số (Digital Twin) với 
mạng học sâu (Deep Learning) phục vụ giám sát tự động bằng thiết bị không người lái.
==========================================================================================

Kết luận

Luận án tiến sĩ của tác giả Cao Văn Mai là một công trình khoa học xuất sắc, mẫu mực về tính chuẩn xác lý thuyết và tính ứng dụng thực tiễn, đóng góp toàn diện cho ngành Cơ học Kỹ thuật và Kỹ thuật Kết cấu:

  1. Xây dựng thành công công thức chính xác dạng đóng của hàm phổ phản ứng chuyển vị cho dầm Euler-Bernoulli đồng nhất chứa hệ thống nhiều vết nứt với các điều kiện biên ngàm, gối tựa tổng quát.
  2. Đề xuất và thiết lập nghiệm giải tích của Hàm độ cong phổ phản ứng, chứng minh bằng toán học và thực nghiệm khả năng khuếch đại vượt trội để nhận diện các vết nứt vi mô ($\ge 10% h$).
  3. Thiết lập biểu thức giải tích chính xác hàm phổ phản ứng cho dầm có cơ tính biến đổi dọc trục (AFG) mang khối lượng tập trung, mở đường cho việc chẩn đoán phân bố vật liệu và phụ tải bám.
  4. Xây dựng hệ thống mô phỏng số PTHH đối chứng độc lập, khẳng định sự hội tụ và chính xác tuyệt đối của các thuật toán giải tích đã phát triển.
  5. Thực hiện thành công hệ thực nghiệm kiểm chứng không phá hủy NDT với thiết bị kích rung Brüel & Kjær Type 4808 và giao thoa kế Laser Doppler PDV-100, xác thực hoàn toàn tính khả thi công nghệ.
  6. Mở ra các nhánh nghiên cứu mới trong lĩnh vực chẩn đoán sức khỏe kết cấu thông minh, đóng góp thiết thực cho công cuộc bảo đảm an toàn công trình hạ tầng trọng điểm quốc gia và quốc tế.