Tổng quan về luận án

Sự toàn vẹn kết cấu của các công trình cơ khí, xây dựng dân dụng, cầu đường và hàng không luôn phải đối mặt với nguy cơ suy giảm khả năng chịu tải do tải trọng lặp, tác động môi trường và thảm họa thiên nhiên (như động đất). Dưới tác động của tải trọng mỏi, các vết nứt tế vi hình thành, tích tụ và phát triển thành hư hỏng cục bộ đe dọa sự ổn định tổng thể của toàn hệ thống kết cấu. Luận án tiến sĩ kỹ thuật cơ khí và cơ kỹ thuật của tác giả Nguyễn Văn Quang với đề tài "Phát triển và ứng dụng các phương pháp phân tích tín hiệu trong chẩn đoán vết nứt kết cấu hệ thanh" (Chuyên ngành Cơ kỹ thuật, Mã số: 9520101, thực hiện tại Học viện Khoa học và Công nghệ - Viện Hàn lâm Khoa học và Công nghệ Việt Nam dưới sự hướng dẫn của PGS.TS. Nguyễn Việt Khoa) định vị là một công trình tiên phong trong lĩnh vực Giám sát Sức khỏe Kết cấu (Structural Health Monitoring - SHM) và Đánh giá Không phá hủy dựa trên Dao động (Vibration-based Non-Destructive Testing).

Khoảng trống nghiên cứu (research gap) cốt lõi mà luận án giải quyết xuất phát từ các giới hạn căn bản của y văn quốc tế:

  1. Các phương pháp chẩn đoán truyền thống dựa trên sự thay đổi tần số dao động riêng (natural frequency) thường chỉ nhạy khi tỷ lệ độ sâu vết nứt vượt quá ngưỡng nghiêm trọng ($a/h \ge 40% - 50%$), trong khi các vết nứt nhỏ ban đầu gây ra độ suy giảm tần số không đáng kể ($< 1% - 2%$), dễ bị triệt tiêu bởi nhiễu môi trường (Chondros et al., 1998; Kisa et al., 1998).
  2. Các phương pháp nhận dạng dựa trên dạng riêng (mode shapes) đòi hỏi mật độ mạng lưới cảm biến dày đặc, quy trình đo đạc tốn kém, và khó áp dụng cho các liên kết giàn không gian hoặc kết cấu phức tạp bị che khuất (Pandey et al., 1991; Haritos et al., 2004).
  3. Các kỹ thuật phân tích Wavelet kinh điển chủ yếu khai thác đáp ứng tĩnh hoặc dạng riêng theo tọa độ không gian để định vị vết nứt tĩnh, hoàn toàn thiếu vắng khả năng phát hiện thời điểm xuất hiện tức thời của vết nứt đột ngột trong quá trình kích động động lực học thời gian thực (như động đất).
  4. Chưa có giải pháp tích hợp hiện tượng tương tác giữa khối lượng tập trung di động và hàm tương quan vị trí - tần số (Frequency-Mass Location - FML) để khuếch đại độ nhạy vết nứt trong hệ dầm kép.
  5. Việc xác định ma trận độ cứng thực nghiệm từ trước đến nay đều phải thông qua bước trung gian là trích xuất dạng dao động riêng, dẫn đến sai số lan truyền lớn; thiếu một phương pháp tái thiết trực tiếp chỉ số phân bố độ cứng phần tử (Element Stiffness Distribution - ESD) từ các hàm đáp ứng tần số (FRF).

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để nhận dạng đồng thời sự tồn tại, vị trí và thời điểm xuất hiện của vết nứt đột ngột trong kết cấu dầm chịu tải trọng động đất mà không làm mất thông tin thời gian?
  • Câu hỏi nghiên cứu 2 (RQ2): Sự tích hợp của khối lượng tập trung di động trên hệ dầm kép tương tác đàn hồi có thể khuếch đại độ nhạy của phép biến đổi Wavelet để xác định vết nứt nông như thế nào?
  • Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để xây dựng ma trận độ cứng và tái cấu trúc phân bố chỉ số suy giảm độ cứng cục bộ trực tiếp từ tín hiệu miền tần số mà không cần khôi phục dạng dao động riêng?
  • Giả thuyết khoa học (H1): Phổ năng lượng Wavelet liên tục (CWT Power Scalogram) kết hợp với đường đặc trưng tần số tức thời (Instantaneous Frequency - IF) sẽ tạo ra bước nhảy năng lượng cục bộ $\Delta f$ tại thời điểm nứt đột ngột, với biên độ đỉnh tỷ lệ thuận trực tiếp với độ sâu vết nứt $a/h$.
  • Giả thuyết khoa học (H2): Phép biến đổi Wavelet trên hàm tương quan tần số - vị trí khối lượng (FML) của hệ dầm kép tạo ra điểm kỳ dị cực đại tại đúng tọa độ vết nứt $L_c$.
  • Giả thuyết khoa học (H3): Hàm phân bố chỉ số độ cứng phần tử trong miền tần số cho phép định vị và định lượng chính xác độ suy giảm độ cứng của từng thanh trong kết cấu khung/giàn không gian 3D ngay cả trong điều kiện tín hiệu bị ô nhiễm bởi nhiễu trắng Gaussian.

Khung lý thuyết của luận án xây dựng trên nền tảng: Cơ học phá hủy đàn hồi tuyến tính (Linear Elastic Fracture Mechanics - LEFM), Lý thuyết dầm Euler-Bernoulli, Nguyên lý công ảo và định lý Castigliano, Phương pháp Phần tử Hữu hạn (FEM) cho dầm 2D/3D chứa vết nứt, và Phép biến đổi Wavelet liên tục (Continuous Wavelet Transform - CWT). Phạm vi thực nghiệm và mô phỏng số bao gồm: dầm 2D, hệ dầm kép mang khối lượng tập trung, hệ khung phẳng và mô hình kết cấu giàn cao tầng không gian 3D chịu kích động xung và động đất trên bàn rung, sử dụng hệ thống thu thập và phân tích dữ liệu đa kênh chuyên dụng Brüel & Kjær PULSE tại Viện Cơ học (VAST).

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy lịch sử phát triển của các phương pháp chẩn đoán hư hỏng dựa trên tham số động lực học kết cấu chia thành ba dòng nghiên cứu chính:

Dòng nghiên cứu 1: Nhận dạng vết nứt dựa trên tần số riêng và mô hình suy giảm độ cứng cục bộ. Chondros và đồng nghiệp (1998) [1] đã phát triển lý thuyết dao động cho dầm Euler-Bernoulli có vết nứt mở, mô hình hóa vết nứt như một lò xo xoay cục bộ dựa trên trường biến dạng cơ học phá hủy, chứng minh tần số tự nhiên giảm khi độ sâu vết nứt tăng. Lee và đồng nghiệp (2000) [2] kết hợp ma trận độ mềm suy ra từ cơ học phá hủy với phương pháp Armon’s Rank-ordering và phương trình Gudmundson để tìm vị trí và độ sâu nứt, tuy nhiên đã kết luận rằng: "tần số riêng thay đổi nhỏ khi độ sâu vết nứt là nhỏ và chỉ thay đổi đáng kể khi độ sâu vết nứt lớn đến 40% độ cao của dầm". Kisa và đồng nghiệp (1998) [6] sử dụng FEM kết hợp tổng hợp dạng dao động thành phần, Saez và cộng sự (2002) [7] mô hình hóa khớp đàn hồi, đều ghi nhận sự bất khả thi trong việc phát hiện vết nứt nông ($a/h < 30%$) nếu chỉ dùng phổ Fourier truyền thống.

Dòng nghiên cứu 2: Phân tích vết nứt đóng - mở (breathing cracks) và ứng xử phi tuyến song tuyến tính. Carlson (1967) [19] và Gudmundson (1983) [20] chứng minh rằng dưới tác động của vết nứt "thở", sự thay đổi tần số dao động là cực nhỏ khi độ sâu $a/h < 50%$ do hai mép nứt tiếp xúc nhau trong chu kỳ nén làm hồi phục độ cứng. Chondros và cộng sự (2001) [22, 23] chỉ ra tính chất song tuyến tính (bilinear stiffness) của kết cấu có vết nứt mỏi, trong đó: "khi độ sâu vết nứt là 40%, sự thay đổi tần số thấp nhất cho vết nứt mở là 1.9%, còn đối với vết nứt đóng - mở chỉ là 0.5%". Cheng và đồng nghiệp (1999) [24] khẳng định với $a/h = 30%$, độ dịch chuyển tần số của vết nứt thở chỉ đạt $1%$. Bovsunovsky (2004) [46, 47], Ruotolo (1997) [48], Rivola (1998) [27] tiếp tục mở rộng các phân tích phi tuyến miền thời gian và hàm đáp ứng tần số bậc cao, tuy nhiên các công trình này đều bế tắc trong việc thiết lập thuật toán định lượng tự động tọa độ và kích thước nứt nếu không có dữ liệu thực nghiệm chuẩn ban đầu.

Dòng nghiên cứu 3: Khai thác biến đổi Wavelet (CWT, DWT, SWT) trong chẩn đoán không phá hủy. Để vượt qua giới hạn mất thông tin thời gian của biến đổi Fourier tích phân (FFT), Ovanesova và Suarez (2004) [59] đã áp dụng wavelet bior6.8 trên đường cong chuyển vị tĩnh để định vị đỉnh dị thường. Gentile và Messina (2003) [67] phối hợp đạo hàm, tích chập và làm trơn trong CWT để nhận dạng vết nứt từ dạng riêng uốn. Loutridis và cộng sự (2004, 2005) [68, 82] khai thác số mũ Lipschitz (Lipschitz exponent) và hàm Mexican Hat để đánh giá mức độ suy giảm độ cứng từ dạng riêng. Rucka và Wilde (2006) [64], Wu và Wang (2011) [104] sử dụng xử lý ảnh và CWT chuyển vị tĩnh để phát hiện vết nứt dầm. Tuy nhiên, toàn bộ các nghiên cứu này đòi hỏi phải đo đạc dạng riêng không gian với độ phân giải cực cao, rất nhạy cảm với nhiễu thực nghiệm, và hoàn toàn không xử lý được bài toán vết nứt phát sinh đột ngột trong biến cố động học.

Vị thế khoa học của luận án được xác lập thông qua việc khắc phục triệt để các hạn chế trên: Tác giả phát triển phương pháp phổ năng lượng Wavelet liên tục để nắm bắt đồng thời tọa độ thời gian và biên độ biến thiên tần số tức thời khi xuất hiện vết nứt tức thời; thiết lập mô hình hệ dầm kép tích hợp khối lượng tập trung di động tạo hiệu ứng cộng hưởng cục bộ; và kiến tạo thuật toán phân bố chỉ số độ cứng phần tử trực tiếp từ đáp ứng tần số (FRF), loại bỏ hoàn toàn sự phụ thuộc vào dạng dao động riêng.

Đó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 nền tảng Cơ học phá hủy đàn hồi tuyến tính của Irwin-Griffith và Lý thuyết dao động kết cấu dầm Euler-Bernoulli thông qua việc xây dựng hệ thống ma trận độ cứng giải tích chính xác cho phần tử hữu hạn dầm phẳng 2D và dầm không gian 3D chứa vết nứt mở ngang.

Mô hình năng lượng biến dạng bổ sung do vết nứt gây ra được thiết lập từ mật độ giải phóng năng lượng biến dạng $J_s$: $$J_s = \frac{1}{E'} \left[ \left(\sum K_{Ii}\right)^2 + \left(\sum K_{IIi}\right)^2 + (1+\nu)\left(\sum K_{IIIi}\right)^2 \right]$$ Trong đó $E' = E$ đối với trạng thái ứng suất phẳng và $E' = E/(1-\nu^2)$ đối với biến dạng phẳng; $K_I, K_{II}, K_{III}$ là các hệ số cường độ ứng suất tương ứng với ba kiểu nứt cơ bản: Kiểu I (mở), Kiểu II (trượt dọc) và Kiểu III (xé rách).

Đối với phần tử dầm 2D chịu uốn và cắt, năng lượng biến dạng bổ sung $W^{(1)}$ do vết nứt có độ sâu $a$ trên bề rộng dầm $b$ được tích phân chính xác: $$W^{(1)} = b \int_0^a \frac{(K_{IM} + K_{IP})^2 + K_{IIP}^2}{E'} da$$ Với các hàm hiệu chỉnh hình học không thứ nguyên $F_I(s)$ và $F_{II}(s)$ ($s = a/h$): $$F_I(s) = \sqrt{\frac{\tan(\pi s / 2)}{\pi s / 2}} \frac{0.923 + 0.199(1 - \sin(\pi s / 2))^4}{\cos(\pi s / 2)}$$ $$F_{II}(s) = \frac{1.122 - 0.561s + 0.085s^2 + 0.18s^3}{\sqrt{1 - s}}$$

Từ định lý Castigliano, ma trận độ mềm bổ sung của phần tử nứt $c^{(1)}{ij} = \frac{\partial^2 W^{(1)}}{\partial P_i \partial P_j}$ được ghép nối trực tiếp với ma trận độ mềm phần tử nguyên vẹn $c^{(0)}{ij}$, từ đó xác lập ma trận độ cứng chính xác của phần tử chứa vết nứt thông qua ma trận chuyển đổi cân bằng $T$: $$K_c = T^T \left( c^{(0)} + c^{(1)} \right)^{-1} T$$

   P1, M1                                                     P2, M2

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

Khung phân tích của luận án tích hợp liên ngành ba trụ cột lý thuyết: (1) Cơ học kết cấu tính toán (Computational Structural Mechanics), (2) Lý thuyết xử lý tín hiệu số nâng cao miền Thời gian - Tần số (Time-Frequency Signal Processing), và (3) Lý thuyết bài toán ngược trong chẩn đoán kỹ thuật (Inverse Problems in SHM).

  1. Kỹ thuật Phổ Năng lượng Wavelet Liên tục (CWT Power Scalogram): Khai thác hàm Wavelet mẹ Mexican Hat $\psi(t) = (1 - t^2)e^{-t^2/2}$ và Morlet để phân tích tín hiệu gia tốc động lực học. Tần số tức thời $IF(t)$ được trích xuất dọc theo các đường cực đại (ridge lines) của phổ wavelet, phản ánh chính xác sự suy giảm độ cứng nhảy bậc tại thời điểm vết nứt hình thành.
  2. Lý thuyết Tương quan Vị trí Khối lượng - Tần số (FML) trên Hệ Dầm Kép: Khi một khối lượng tập trung $M_c$ di chuyển dọc theo hệ dầm kép liên kết đàn hồi Winkler $K$, hàm tần số tự nhiên $f_1(x_{mass})$ xuất hiện điểm gián đoạn đạo hàm bậc nhất tại đúng tọa độ vết nứt $L_c$. Việc áp dụng biến đổi Wavelet không gian lên chuỗi dữ liệu FML đóng vai trò như một bộ lọc thông cao khuếch đại điểm kỳ dị này thành đỉnh nhọn đột biến (peak), cho phép định vị vết nứt ngay cả khi tỷ số nứt cực nhỏ ($a/h = 5%$).
  3. Phương pháp Phân bố Độ cứng Phần tử trong Miền Tần số (Frequency-Domain ESD Index): Thiết lập thuật toán tái tạo trực tiếp độ cứng từng thanh từ ma trận hàm truyền đáp ứng tần số $H(\omega)$: $${X(\omega)} = \left[ -\omega^2 [M] + i\omega [C] + [K] \right]^{-1} {F(\omega)}$$ Định nghĩa chỉ số biến thiên độ cứng phần tử $\Delta k_e$ thông qua sự sai lệch giữa năng lượng biến dạng phần tử đo đạc và mô hình gốc mà không cần trung gian chuyển đổi sang vector riêng.

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

Thiết kế nghiên cứu

Luận án vận hành theo hệ hình thực chứng nghiêm ngặt (strict positivist paradigm) với lập trường nhận thức luận khách quan (objective epistemology), kết hợp đa tầng giữa: (1) Mô phỏng giải tích và số trị phần tử hữu hạn (FEM/ANSYS/MATLAB), và (2) Thực nghiệm kiểm chứng vật lý trên mô hình thực tại phòng thí nghiệm chuyên sâu.

Thiết kế nghiên cứu đa cấp độ bao gồm 4 cấu hình kết cấu tiêu biểu:

  • Cấu hình 1: Dầm công xôn và dầm đơn giản 2D chịu gia tốc nền động đất El Centro / kích động sóng điều hòa.
  • Cấu hình 2: Hệ hai dầm song song (hệ dầm kép) liên kết bởi lớp đệm đàn hồi liên tục $K = 10^5 \text{ N/m}^2$, mang khối lượng tập trung $M_c = 0.1 M_{beam}$ di động.
  • Cấu hình 3: Khung phẳng chịu tải trọng ngang và tải trọng điều hòa.
  • Cấu hình 4: Kết cấu giàn cao tầng không gian 3D gồm nhiều tầng thanh giằng phức tạp.

Quy trình nghiên cứu rigorous

  1. Quy trình thu thập dữ liệu số và mô phỏng: Phương trình vi phân dao động tổng thể $[M]{\ddot{u}} + [C]{\dot{u}} + [K]{u} = {F(t)}$ được giải bằng phương pháp tích phân từng bước Newmark-$\beta$ và phương pháp lặp cân bằng lực. Ma trận cản được thiết lập theo mô hình cản tỷ lệ Rayleigh $[C] = \alpha[M] + \beta[K]$.
  2. Quy trình thực nghiệm vật lý: Thí nghiệm kiểm chứng được triển khai tại Phòng thí nghiệm Động lực học Công trình - Viện Cơ học (VAST).
    • Mô hình dầm thép có kích thước hình học chính xác, được gá đặt trên bàn rung điều khiển tần số để mô phỏng sự cố nứt gãy tức thời thông qua cơ cấu nhả chốt cơ học tự động.
    • Mô hình giàn thép không gian 3D được kích thích bằng búa gõ động lực học (Impact Hammer) chuyên dụng có gắn đầu đo lực.
    • Tín hiệu gia tốc phản hồi được thu thập bằng hệ thống gia tốc kế áp điện độ nhạy cao kết nối với máy phân tích phổ đa kênh Brüel & Kjær PULSE (Đan Mạch). Dữ liệu được số hóa với tần số lấy mẫu $f_s = 2048 \text{ Hz}$ nhằm tránh hiện tượng méo chồng phổ (aliasing).

Data và phân tích

  • Khảo sát ma trận dữ liệu phong phú với các cấp độ nứt biến thiên từ $a/h = 5%, 10%, 20%, 30%, 40%, 50%$ tại các tọa độ vết nứt $L_c/L = 0.2; 0.4; 0.5; 0.6; 0.8$.
  • Kiểm tra độ bền vững của thuật toán (Robustness Check) bằng cách bổ sung nhiễu trắng Gaussian nhân tạo (Additive White Gaussian Noise) ở các mức $1%, 3%, 5%$ và $10%$ vào tín hiệu đáp ứng gia tốc mô phỏng.
  • Phân tích tương quan hồi quy phi tuyến giữa độ suy giảm tần số tức thời $\Delta f$, chiều cao đỉnh phổ Wavelet $dh$, và chỉ số phân bố độ cứng phần tử với độ sâu thực tế của vết nứt.

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

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

  1. Khám phá quy luật tương quan tuyến tính giữa bước nhảy tần số tức thời $\Delta f$ và bình phương độ sâu vết nứt ($s^2$): Trong bài toán dầm chịu tác động của động đất, phân tích phổ Wavelet liên tục cho thấy tần số tức thời $IF(t)$ giữ giá trị không đổi trước thời điểm nứt $t_c$, sau đó sụt giảm đột ngột tại đúng thời điểm xuất hiện vết nứt $t = t_c$ và tiếp tục dao động ở mức tần số mới thấp hơn. Bước nhảy tần số $\Delta f = f_{before} - f_{after}$ tỷ lệ đơn điệu và phi tuyến bậc hai với độ sâu vết nứt: $$\Delta f \sim \gamma \cdot \left(\frac{a}{h}\right)^2$$ Điều này cung cấp một công thức tường minh để trích xuất trực tiếp kích thước vết nứt từ tín hiệu đo gia tốc đơn điểm mà không cần biết trước tải trọng kích động.

  2. Cơ chế khuếch đại điểm kỳ dị FML trên hệ dầm kép: Khi khối lượng tập trung $M_c$ di chuyển qua vị trí vết nứt $L_c$, tần số dao động riêng thứ nhất $f_1$ của hệ dầm kép xuất hiện sự biến dạng cục bộ dạng gãy khúc. Sau khi áp dụng biến đổi Wavelet (hàm Mexican Hat) lên đường FML, một đỉnh nhọn đột biến $dh$ xuất hiện với tọa độ đỉnh trùng khớp $100%$ với vị trí $L_c$. Chiều cao đỉnh $dh$ tăng mạnh và tuyến tính theo tỷ lệ độ sâu vết nứt từ $a/h = 10%$ đến $50%$.

  3. Khả năng định vị vết nứt chính xác tuyệt đối của phương pháp Phân bố Độ cứng Phần tử (ESD): Trong mô hình giàn cao tầng không gian 3D, khi thanh số #17 bị nứt ngầm với các độ sâu khác nhau ($10% - 50%$), biểu đồ tái cấu trúc chỉ số độ cứng phần tử cho thấy một đỉnh vượt ngưỡng đột biến tại đúng phần tử #17. Ngay cả khi dữ liệu đo bị ô nhiễm bởi $5%$ nhiễu ngẫu nhiên, đỉnh chỉ số $\Delta k_{17}$ vẫn vượt trội hoàn toàn so với nền nhiễu của các phần tử lân cận, loại trừ hoàn toàn các báo động giả (false positives).

  4. Bằng chứng thực nghiệm xác thực độ tin cậy của thuật toán: Kết quả đo đạc thực nghiệm trên hệ thống PULSE tại Viện Cơ học hoàn toàn khớp với mô phỏng giải tích:

  • Sai số xác định thời điểm xảy ra vết nứt: $e_t < 0.005 \text{ giây}$.
  • Sai số định vị phần tử nứt trong kết cấu giàn không gian: $0%$ (nhận dạng chính xác $100%$ thanh #17 bị tổn thương).
  • Sai số định lượng độ sâu vết nứt: $e_a < 3.8%$.
Thông số chẩn đoán Kết quả mô phỏng FEM Dữ liệu thực nghiệm (PULSE/VAST) Độ sai lệch tương đối
Vị trí thanh nứt (Giàn không gian) Phần tử #17 Phần tử #17 0.00% (Tuyệt đối)
Thời điểm nứt dầm $t_c$ 2.500 s 2.512 s 0.48%
Độ suy giảm tần số $\Delta f$ ($a/h=30%$) 1.45 Hz 1.41 Hz 2.76%
Biên độ đỉnh chỉ số $dh$ ($a/h=40%$) 0.082 0.079 3.66% 
  Biên độ Chỉ số Độ cứng (dh)
             #1    #3  #5    #7  #9   #11 #13   #15  #17 #19 #21  (Số thứ tự phần tử giàn)

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cầu nối toán học hoàn chỉnh giữa Cơ học phá hủy phi tuyến và Lý thuyết biểu diễn tín hiệu Wavelet nhiều mức, chứng minh rằng các biến đổi vi mô cục bộ hoàn toàn có thể giải mã được từ các tín hiệu đo động lực học vĩ mô.
  • Về mặt phương pháp luận: Thiết lập một quy trình chẩn đoán không xâm lấn trực tiếp từ dữ liệu FRF, giúp cộng đồng khoa học loại bỏ hoàn toàn các thuật toán nhận dạng mode shape cồng kềnh, giảm thiểu chi phí tính toán đến $60%$.
  • Về mặt thực tiễn kỹ thuật: Cho phép chế tạo các mô-đun giám sát thông minh gắn tại các công trình trọng điểm (cầu dây văng, dàn khoan ngoài khơi, kết cấu tháp truyền hình, turbine gió), giúp cảnh báo sớm sự cố đứt gãy ngay trong thời gian thực khi xảy ra bão cát, va chạm tàu bè hoặc động đất.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu mang tính biên điều kiện (boundary conditions):

  1. Giả định hình học vết nứt: Mô hình giải tích tập trung vào vết nứt mở phẳng 1 cạnh, có phương vuông góc với trục thanh, chiều sâu không đổi dọc theo bề rộng, chưa xét đến các vết nứt xiên không gian phức tạp hoặc vết nứt phân nhánh (branching cracks).
  2. Bản chất vật liệu: Nghiên cứu thực hiện trên vật liệu kim loại đồng nhất đẳng hướng tuyến tính (thép kết cấu), chưa mở rộng sang vật liệu composite nhiều lớp có hiện tượng tách lớp (delamination) hoặc bê tông cốt thép có tính phi tuyến vật liệu cao.
  3. Điều kiện môi trường: Các kiểm chứng thực nghiệm được thực hiện trong môi trường phòng thí nghiệm với nhiệt độ và độ ẩm được kiểm soát, chưa xét đến biến thiên nhiệt độ môi trường hàng ngày - yếu tố có thể gây dịch chuyển tần số tự nhiên tương đương mức suy giảm do vết nứt gây ra.
  4. Mô hình tiếp xúc vết nứt thở: Chưa tích hợp đầy đủ ma sát coulomb giữa hai bề mặt mép nứt trong trạng thái đóng - trượt ma sát phức tạp.

Chương trình nghiên cứu tiếp nối trong tương lai (Future Research Agenda):

  • Phát triển mô hình giải tích FEM cho vết nứt xiên 3D kết hợp vết nứt mỏi dạng thở có xét đến lực ma sát bề mặt tiếp xúc.
  • Ứng dụng thuật toán phân bố độ cứng phần tử (ESD) cho kết cấu vật liệu composite dạng tấm/vỏ và vật liệu biến tính chức năng (FGM - Functionally Graded Materials).
  • Tích hợp mạng học sâu (Deep Learning/Convolutional Neural Networks) với đầu vào là ảnh phổ năng lượng Wavelet Scalogram để tự động hóa quy trình phân loại hư hỏng thời gian thực.
  • Thử nghiệm hệ thống giám sát không dây quy mô lớn trên các công trình cầu nhịp lớn và tháp gió ngoài khơi thực tế chịu tác động của tải trọng ngẫu nhiên do gió và sóng biển.

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

Luận án của tác giả Nguyễn Văn Quang tạo ra những tác động khoa học và ứng dụng thực tiễn sâu rộng:

  • Ảnh hưởng học thuật: Đóng góp các bài báo chất lượng cao trên các tạp chí chuyên ngành Cơ học và Kỹ thuật kết cấu uy tín; mở ra hướng nghiên cứu kết hợp Wavelet với hàm truyền FRF trực tiếp trong cộng đồng SHM tại Việt Nam và khu vực.
  • Chuyển giao công nghiệp: Cung cấp công cụ tính toán cốt lõi cho các viện nghiên cứu thiết kế cơ khí, hàng không vũ trụ và các tổng công ty tư vấn xây dựng công trình giao thông thủy lợi trong việc đánh giá tuổi thọ còn lại (Remaining Useful Life - RUL) của kết cấu giàn thép, khung dầm chịu tải trọng nặng.
  • Chính sách và tiêu chuẩn an toàn: Cung cấp luận cứ khoa học vững chắc phục vụ việc xây dựng các tiêu chuẩn kỹ thuật quốc gia về quan trắc sức khỏe công trình hạ tầng trọng điểm, giảm thiểu nguy cơ thảm họa sụp đổ công trình dân dụng trong vùng rủi ro địa chấn cao.

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

  • Nghiên cứu sinh và Giới học thuật (Academic Researchers): Thụ hưởng khung lý thuyết hoàn chỉnh về phần tử hữu hạn dầm 2D/3D chứa vết nứt và thuật toán phân tích Wavelet liên tục chuyên sâu, có thể kế thừa trực tiếp để phát triển cho các bài toán kết cấu phức tạp hơn.
  • Kỹ sư R&D và Giám sát Công trình (Industry Practitioners): Nhận được thuật toán phân bố độ cứng phần tử (ESD) tinh gọn, cho phép triển khai trực tiếp vào phần mềm nhúng của các thiết bị quan trắc dao động thương mại với số lượng cảm biến tối thiểu.
  • Nhà quản lý hạ tầng và cơ quan hoạch định chính sách: Sở hữu giải pháp công nghệ nội địa hóa với chi phí thấp, giúp tối ưu hóa ngân sách bảo trì định kỳ các cây cầu thép, cảng biển và dàn khoan dầu khí quốc gia.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Động lực học Cơ học Phá hủy thông qua việc thiết lập giải tích ma trận độ mềm bổ sung $c^{(1)}$ cho phần tử dầm 3D tổng quát chịu đồng thời uốn 2 trục, cắt 2 phương, kéo-nén và xoắn thuần túy, kết hợp xây dựng hàm phân bố chỉ số độ cứng phần tử (ESD) trong miền tần số. Công trình mở rộng trực tiếp lý thuyết dầm nứt của Chondros et al. (1998) và Zheng et al. (2004) từ hệ phẳng 2D sang không gian 3D phức tạp.

2. Sự đổi mới về phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm?

So với nghiên cứu của Ovanesova (2004) và Gentile (2003) vốn đòi hỏi đo đạc toàn bộ dạng riêng không gian, phương pháp của luận án:

  1. Cho phép nhận dạng vết nứt từ tín hiệu đo gia tốc đơn điểm thông qua phổ Wavelet liên tục;
  2. Tái tạo trực tiếp ma trận độ cứng từ hàm FRF mà không cần trích xuất dạng dao động riêng, triệt tiêu sai số làm mịn và tiết kiệm $80%$ số lượng cảm biến so với phương pháp của Haritos (2004) hay Verboven (2002).

3. Phát hiện gây bất ngờ nhất có số liệu thực nghiệm chứng minh là gì?

Phát hiện bất ngờ nhất là hiện tượng: Tần số dao động tức thời $IF(t)$ trong miền Wavelet không trôi dần đều mà xuất hiện một bước nhảy gián đoạn dốc đứng (vertical drop) tại đúng mili-giây xảy ra nứt đột ngột trong biến cố động đất mô phỏng, với độ lệch thời gian thực nghiệm so với lý thuyết chỉ $0.012 \text{ giây}$ ($< 0.5%$), chứng minh tính khả thi của việc phát hiện "thời khắc vàng" xảy ra phá hủy kết cấu.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án cung cấp tường minh:

  • Toàn bộ công thức giải tích tích phân ma trận độ mềm phần tử dầm $K_c$;
  • Các phương trình giải tích sóng Wavelet Mexican Hat và Morlet;
  • Sơ đồ khối thuật toán phân bố độ cứng phần tử;
  • Thông số chi tiết của thiết bị đo đạc thực nghiệm (máy PULSE, cảm biến gia tốc áp điện, búa gõ kích động, tần số lấy mẫu $2048 \text{ Hz}$).

5. Định hướng lộ trình nghiên cứu 10 năm được vạch ra như thế nào?

Lộ trình 10 năm tập trung vào:

  1. Mở rộng thuật toán ESD cho vật liệu dị hướng Functionally Graded Materials (FGM) và kết cấu vỏ mỏng hàng không;
  2. Xây dựng nền tảng đám mây Digital Twin kết hợp IoT thu thập dữ liệu rung động thời gian thực từ cầu đường thực tế;
  3. Tích hợp AI/Machine Learning để tự động hóa hoàn toàn quy trình ra quyết định cảnh báo nguy hiểm trong quản lý hạ tầng đô thị thông minh.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Văn Quang đã hoàn thành xuất sắc các mục tiêu khoa học đề ra, đóng góp 5 giá trị đột phá then chốt:

  1. Thiết lập thành công mô hình phần tử hữu hạn 2D và 3D giải tích cho dầm chứa vết nứt dựa trên lý thuyết cơ học phá hủy đàn hồi tuyến tính chính xác.
  2. Tiên phong ứng dụng phổ năng lượng Wavelet liên tục để phát hiện đồng thời vị trí, độ sâu và thời điểm xuất hiện của vết nứt đột ngột trong quá trình kích động động đất.
  3. Khám phá và khai thác hiệu quả quy luật tương tác Vị trí khối lượng - Tần số (FML) trên hệ dầm kép để khuếch đại độ nhạy nhận dạng các vết nứt nông ($a/h \ge 5%$).
  4. Kiến tạo phương pháp Phân bố độ cứng phần tử trong miền tần số (ESD), cho phép định lượng chính xác suy giảm độ cứng cục bộ của kết cấu giàn không gian đa tầng trực tiếp từ hàm FRF mà không cần trích xuất dạng riêng.
  5. Kiểm chứng thực nghiệm thành công toàn diện các thuật toán đề xuất trên hệ thống thiết bị đo động lực học hiện đại Brüel & Kjær PULSE tại Viện Cơ học (VAST).

Công trình đánh dấu bước tiến quan trọng của chuyên ngành Cơ kỹ thuật Việt Nam trong việc làm chủ các công nghệ chẩn đoán kỹ thuật không phá hủy tiên tiến, khẳng định năng lực hội nhập học thuật quốc tế đỉnh cao và mang lại giá trị thực tiễn to lớn cho sự an toàn của các công trình công nghiệp và dân dụng trọng điểm quốc gia.