Tổng quan về luận án

Trong kỷ nguyên phát triển vượt bậc của kỹ thuật xây dựng hiện đại, quy mô và độ phức tạp của các công trình cơ sở hạ tầng ngày càng gia tăng, tiêu biểu như siêu tháp Burj Khalifa cao 828 m (UAE), cầu treo Akashi-Kaikyo nhịp chính 3.911 m (Nhật Bản), hay tháp Tokyo Sky Tree cao 634 m. Đi cùng với quy mô đồ sộ là thách thức to lớn về độ an toàn kết cấu trước tác động cộng hưởng từ môi trường tự nhiên (động đất, bão gió, ăn mòn hóa học) và tải trọng khai thác lặp chu kỳ. Sự cố sụp đổ cầu vượt tại Canada (2006), cầu I-35W qua sông Mississippi tại Mỹ (2007), hay tại Việt Nam là sự cố tháp truyền hình Nam Định cao 180 m (2012) và cầu treo Chu Va 6 tại Lai Châu (2014) là những minh chứng điển hình cho hậu quả nghiêm trọng khi các khuyết tật tích tụ không được nhận diện kịp thời. Do đó, việc thiết lập hệ thống giám sát sức khỏe kết cấu (Structural Health Monitoring – SHM) theo thời gian thực đóng vai trò sống còn trong việc phòng ngừa thảm họa và tối ưu hóa chi phí vòng đời công trình.

Theo phân loại kinh điển của Rytter (1993), SHM bao gồm 4 cấp độ: Cấp độ 1 (Phát hiện sự xuất hiện hư hỏng), Cấp độ 2 (Khoanh vùng vị trí hư hỏng), Cấp độ 3 (Định lượng mức độ hư hỏng), và Cấp độ 4 (Dự báo tuổi thọ còn lại của kết cấu). Trong các cấu kiện chịu lực cơ bản, kết cấu tấm (plate structures) như sàn bê tông cốt thép (BTCT), bản mặt cầu, tường chắn chịu uốn chiếm tỷ trọng lớn nhưng công tác chẩn đoán lại đối mặt với nhiều rào cản kỹ thuật. Phương pháp năng lượng biến dạng dao động (Modal Strain Energy – MSE) kết hợp giải thuật di truyền (Genetic Algorithm – GA) được xem là công cụ hứa hẹn để đạt tới cấp độ 3. Tuy nhiên, các phương pháp hiện hành bộc lộ research gap lớn:

  1. Đa số thuật toán tính toán MSE truyền thống yêu cầu thu thập đủ 5 thành phần chuyển vị nút ($u, v, w, \theta_x, \theta_y$), trong đó thành phần chuyển vị xoay ($\theta_x, \theta_y$) hầu như không thể đo đạc chính xác trên công trình thực tế.
  2. Các công thức xấp xỉ sai phân hữu hạn dựa trên phần tử tứ giác 4 nút (Q4) chịu ảnh hưởng nặng nề bởi điều kiện biên và gây nhiễu loạn chẩn đoán tại các vị trí cận biên hoặc vùng có vết nứt nhỏ (< 10%).
  3. Khi bước khoanh vùng chẩn đoán dư (false positives), thuật toán GA truyền thống rơi vào bế tắc do không gian tìm kiếm đa chiều, dẫn đến hàm mục tiêu không hội tụ.

Luận án của NCS. Lê Thanh Cao (Trường Đại học Bách khoa – ĐHQG TP.HCM, 2023), dưới sự hướng dẫn khoa học của PGS. Hồ Đức Duy, mang tiêu đề "Phát triển phương pháp năng lượng biến dạng để chẩn đoán hư hỏng cho kết cấu tấm" (Chuyên ngành: Kỹ thuật Xây dựng, Mã số: 9580201) đã giải quyết triệt để khoảng trống nghiên cứu này. Luận án đặt ra hệ thống câu hỏi nghiên cứu (RQ) và giả thuyết khoa học (H) cốt lõi:

  • RQ1: Làm thế nào để tính toán chính xác chỉ số năng lượng biến dạng trên kết cấu tấm với điều kiện biên bất kỳ mà chỉ cần duy nhất dữ liệu chuyển vị thẳng theo phương đứng ($w$)?
  • RQ2: Cơ chế lọc dữ liệu dao động nào tối ưu hóa độ nhạy của thuật toán trước các kịch bản đa hư hỏng, hư hỏng nhỏ và hư hỏng cận biên?
  • RQ3: Làm cách nào để giải thuật tiến hóa tự động loại bỏ các phần tử cảnh báo sai từ bước khoanh vùng sơ bộ nhằm định lượng chính xác độ suy giảm độ cứng?
  • RQ4: Khả năng ứng dụng của quy trình cải tiến trong việc nhận diện vết nứt thực tế của tấm sàn BTCT làm việc ngoài giai đoạn đàn hồi dưới tác dụng tải trọng tĩnh tăng dần là gì?

Tương ứng với các câu hỏi trên là 4 giả thuyết khoa học:

  • H1: Việc tích phân số trường chuyển vị uốn thông qua hàm dạng phần tử đẳng tham số 9 nút Lagrange (Q9) sẽ triệt tiêu hoàn toàn sự phụ thuộc vào điều kiện biên và loại bỏ nhu cầu đo chuyển vị xoay.
  • H2: Việc chọn lọc các mode dao động có chỉ số tương quan Modal Assurance Criterion (MAC) cao nhất sẽ nâng cao độ chính xác khoanh vùng so với việc chọn ngẫu nhiên các mode bậc thấp.
  • H3: Thuật toán di truyền lặp (Iterative Genetic Algorithm – I-GA) đa pha cho phép đào thải tuần tự các phần tử không hư hỏng bị nhận diện nhầm, giúp hàm mục tiêu nhanh chóng hội tụ về nghiệm thực.
  • H4: Quy trình hai bước (Global-Local MSE kết hợp I-GA) có khả năng định vị chính xác vị trí nứt bê tông ngay tại cấp tải gây nứt ban đầu ($q_{crack} = 3,134\text{ kN/m}^2$) được kiểm chứng theo tiêu chuẩn ACI 318-02 và mô phỏng phần tử hữu hạn phi tuyến trên ANSYS.

Khung lý thuyết của luận án được xây dựng vững chắc dựa trên Lý thuyết tấm dày Mindlin-Reissner, Cơ học dao động kết cấu công trình, Lý thuyết Năng lượng biến dạng phần tử hữu hạn, Tiêu chuẩn chứng thực dao động MAC (Allemang & Brown, 1982), và Giải thuật tối ưu hóa tiến hóa (Holland, 1975; Goldberg, 1989). Phạm vi nghiên cứu bao gồm mô hình mô phỏng số chi tiết trên tấm nhôm, tấm bê tông (lưới phần tử từ $10\times10$, $20\times20$ đến $36\times24$) và mô hình tấm sàn BTCT chịu uốn phi tuyến thực tế.

Literature Review và Positioning

Phương pháp chẩn đoán hư hỏng dựa vào đặc trưng dao động (Vibration-Based Damage Identification – VBDI) đã trải qua hơn 4 thập kỷ phát triển mạnh mẽ kể từ những ứng dụng đầu tiên trong ngành hàng không vũ trụ và công trình cầu (Doebling et al., 1996; Sohn et al., 2003). Tổng quan các luồng nghiên cứu chính cho thấy sự phân hóa rõ nét giữa các trường phái phương pháp luận:

  1. Nhóm phương pháp dựa vào sự thay đổi tần số tự nhiên: Đại diện bởi Salawu (1997), Messina và cộng sự (1998) với các chỉ số DLAC (Damage Location Assurance Criterion) và MDLAC (Multiple Damage Location Assurance Criterion). Ưu điểm là việc đo đạc tần số có độ chính xác rất cao và chi phí thấp. Tuy nhiên, sự thay đổi tần số khi xuất hiện hư hỏng nhỏ thường dưới 1%, độ nhạy với hư hỏng cục bộ rất kém và chịu tác động lớn bởi biến thiên nhiệt độ môi trường (Farrar et al., 2001). Do đó, nhóm này không đủ năng lực khoanh vùng chính xác đa hư hỏng trong kết cấu hai chiều.
  2. Nhóm phương pháp dựa vào dạng dao động và độ cong dạng dao động: Pandey và cộng sự (1991) mở đầu hướng tiếp cận sử dụng đạo hàm bậc hai của dạng dao động (Mode Shape Curvature) để định vị vết nứt trên dầm, tiếp nối bởi Lieven & Ewins (1988) với chỉ số COMAC (Co-Ordinate Modal Assurance Criterion). Mặc dù độ cong dạng dao động nhạy hơn nhiều so với tần số, việc tính đạo hàm số bậc cao (bậc hai, bậc bốn như Gauthier et al., 2008) từ dữ liệu lưới đo rời rạc dẫn đến sự khuếch đại sai số đo lường (measurement noise amplification) và gặp sai số biên nghiêm trọng (boundary distortion).
  3. Nhóm phương pháp dựa vào Năng lượng biến dạng dao động (Modal Strain Energy - MSE): Được tiên phong bởi Stubbs và cộng sự (1995) trên kết cấu dầm và Cornwell, Doebling, Farrar (1999) mở rộng cho kết cấu tấm. Shi, Law & Zhang (1998, 2000) phát triển chỉ số thay đổi năng lượng biến dạng MSECR (Modal Strain Energy Change Ratio), chứng minh tính ổn định vượt trội trước nhiễu so với phương pháp độ cong và ma trận độ mềm (Flexibility Matrix). Sau đó, các biến thể như MSEBI (Jaishi & Ren, 2006), MSEEI (Seyedpoor, 2012) lần lượt ra đời.

Trong y văn tồn tại hai cuộc tranh luận học thuật sâu sắc:

  • Tranh luận 1: Phương pháp dựa vào mô hình (Model-based / FEM Updating) đối lập Phương pháp phi mô hình (Data-driven / Non-model based). Phái Data-driven (Sohn et al., 2003; Farrar & Worden, 2007) cho rằng mô hình số không thể phản ánh chính xác kết cấu thực tế quy mô lớn và đề xuất sử dụng chuỗi thời gian hoặc học máy thuần túy. Ngược lại, phái Model-based (Friswell & Mottershead, 1995) khẳng định chỉ có mô hình vật lý mới lượng hóa được chính xác mức độ tổn thương kết cấu (Cấp độ 3). Luận án của Lê Thanh Cao dung hòa hai trường phái này thông qua cấu trúc lai (hybrid): sử dụng MSE phi mô hình ở Bước 1 để sàng lọc không gian và FEMU tối ưu hóa bằng GA ở Bước 2.
  • Tranh luận 2: Cách tiếp cận xấp xỉ trường năng lượng biến dạng trong kết cấu tấm. Công trình quốc tế kinh điển của Cornwell và cộng sự (1999) hay Kim & Stubbs (2002) phân chia tấm thành các dải dầm 1D giao nhau và dùng sai phân trung tâm, dẫn đến sai số lớn tại các góc và biên ngàm/tựa. Nghiên cứu của Shi và cộng sự (2000) tiếp cận bằng tích ma trận độ cứng phần tử $[K_e]$ và véc-tơ chuyển vị đầy đủ ${u_e}$, nhưng lại vướng rào cản đo đạc thực tế 5 bậc tự do nút.

Luận án định vị chính xác tại điểm nghẽn này: Đề xuất công thức năng lượng biến dạng dựa trên phần tử đẳng tham số 9 nút (Q9) độc lập hoàn toàn với điều kiện biên và chỉ dùng 1 chuyển vị đứng $w$, đồng thời giải quyết triệt để hiện tượng phân kỳ của GA thông qua giải thuật di truyền lặp đa pha.

Đó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 khung lý thuyết cơ học dao động công trình và nhận dạng hệ thống qua các đóng góp nền tảng:

      LÝ THUYẾT TẤM DÀY MINDLIN-REISSNER (5 DOFs)
      GIẢM BẬC TỰ DO: CHỈ SỬ DỤNG CHUYỂN VỊ ĐỨNG w(x,y)
      PHẦN TỬ ĐẲNG THAM SỐ 9 NÚT (Q9)
      • 9 hàm Lagrange: N_i(ξ, η)
      • Đạo hàm tọa độ tự nhiên + Ma trận Jacobi [J]
      XẤP XỈ NĂNG LƯỢNG BIẾN DẠNG MODAL (MSE)
      • Khử triệt để hiệu ứng biên (boundary-invariant)
      • Khử hoàn toàn sai số sai phân trung tâm
      CHUẨN HÓA CHỈ SỐ HƯ HỎNG (nMSEBI) & BỘ ĐÁNH GIÁ (A, B, C, D)
  1. Mở rộng lý thuyết phần tử hữu hạn trong không gian năng lượng dao động: Thiết lập mối liên hệ giải tích giữa trường chuyển vị uốn $w(x,y)$ và năng lượng biến dạng phần tử $MSE_e$ thông qua toán tử biến dạng uốn $[B_b^e]$ trên lưới 9 nút Lagrange. Định lý xấp xỉ mới chứng minh rằng bằng cách chuyển đổi đạo hàm không gian sang tọa độ tự nhiên $(\xi, \eta)$ qua ma trận chuyển đổi Jacobi $[J]$, giá trị $MSE_e$ không còn bị chi phối bởi các ràng buộc biên hình học (ngàm, gối tựa đơn hay tự do).
  2. Khung lý thuyết tối ưu hóa tìm kiếm thích nghi (Adaptive Search Space Reduction Theory): Luận án chứng minh định lý về sự suy giảm không gian trạng thái trong giải thuật di truyền lặp (I-GA). Khi không gian tìm kiếm của bài toán nghịch đảo có số chiều ban đầu là $N_{var}$, việc loại bỏ từng phần tử $j$ có mức độ hư hỏng $\alpha_j \to 0$ sau mỗi pha $k$ sẽ giảm số chiều xuống $N_{var} - k$, biến bài toán tối ưu phi tuyến không lồi (non-convex) thành chuỗi các bài toán con hội tụ đơn điệu về nghiệm toàn cục.

Hệ thống giả thuyết và mệnh đề lý thuyết (Theoretical Propositions):

  • Proposition 1: Năng lượng biến dạng uốn $MSE_e^k$ của phần tử thứ $e$ tại dạng dao động thứ $k$ biểu diễn chính xác qua tích phân Gauss-Legendre $3\times3$ điểm trên miền tự nhiên $[-1, 1] \times [-1, 1]$ mà không cần giả định chuyển vị xoay nút: $$MSE_e^k = \frac{1}{2} \int_{-1}^{1} \int_{-1}^{1} {\kappa^k}^T [D_b] {\kappa^k} |J(\xi, \eta)| , d\xi , d\eta$$ trong đó ${\kappa^k} = [B_b^e] {w_e^k}$ là véc-tơ độ cong uốn chính, $[D_b]$ là ma trận độ cứng chống uốn của tấm.
  • Proposition 2: Mức độ tương quan dao động $MAC(i, j)$ đóng vai trò hàm trọng số lọc nhiễu tự nhiên; các mode có $MAC \ge 0.90$ mang thông tin suy giảm độ cứng thuần khiết nhất, loại bỏ hiện tượng ghép mode (mode-coupling) do tổn thương cục bộ.

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 giữa 3 trụ cột lý thuyết: (1) Động lực học công trình phi tuyến (Mindlin-Reissner Plate Dynamics), (2) Lý thuyết Độ nhạy Năng lượng biến dạng (Sensitivity-based MSE), và (3) Lý thuyết Tiến hóa Sinh học thích nghi (Adaptive Evolutionary Optimization).

Quy trình vận hành theo cơ chế phân cấp chặt chẽ:

  • Tiêu chuẩn ngưỡng kép (Dual-threshold Criterion): Ở bước tổng thể (Global Step), ngưỡng hư hỏng $\beta_{th1} = 20% \times \max(\beta)$ được áp dụng nhằm tối đa hóa độ nhạy, chấp nhận cảnh báo dư để không bỏ sót bất kỳ vết nứt tiềm tàng nào. Ở bước cục bộ (Local Step), ngưỡng $\beta_{th2} = 30% \times \max(\beta_{local})$ được áp dụng độc lập trên từng vùng phụ đã cô lập, loại bỏ ngay $40-60%$ các phần tử nhiễu xung quanh.
  • Điều kiện biên áp dụng (Boundary Conditions Invariance): Áp dụng đồng nhất cho mọi dạng liên kết: 4 cạnh ngàm (CCCC), 4 cạnh tựa đơn (SSSS), 4 cạnh tự do (FFFF) và các liên kết hỗn hợp (CSCS, CFCF).

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

Thiết kế nghiên cứu

Luận án tuân thủ chặt chẽ thế giới quan thực chứng (Positivism) và chủ nghĩa hiện thực phản biện (Critical Realism) trong khoa học kỹ thuật, kết hợp giữa chứng minh giải tích toán học, mô phỏng số học chi tiết (Numerical Experiments) và kiểm chứng mô hình vật lý phi tuyến. Thiết kế nghiên cứu đa tầng được tổ chức theo cấp độ phức tạp tăng dần:

                                  THIẾT KẾ NGHIÊN CỨU ĐA TẦNG
 [TẦNG 1: MÔ PHỎNG SỐ TẤM NHÔM]          [TẦNG 2: MÔ PHỎNG SỐ TẤM BÊ TÔNG]        [TẦNG 3: MÔ HÌNH SÀN BTCT PHI TUYẾN]
 • Kích thước: 0.6 x 0.4 x 0.002 m       • Kích thước: 1.0 x 1.0 x 0.05 m         • Kích thước: 1.5 x 1.5 x 0.06 m
 • Lưới: 36 x 24 phần tử (864 elements)  • Lưới: 10 x 10 và 20 x 20 phần tử       • Phần tử SOLID65 (Bê tông) & LINK180 (Thép)
 • Khảo sát biên: SSSS, CCCC, FFFF       • Khảo sát: Đa hư hỏng, hư hỏng nhỏ (10%)• Gia tải tĩnh tăng dần đến phá hoại
 • So sánh: Phần tử Q4 vs Phần tử Q9     • Đánh giá: 4 Hàm mục tiêu OF1-OF4       • Chẩn đoán vết nứt tại q1, q2, q3
  • Tầng 1 - Kiểm chứng tính bất biến của phần tử Q9: Khảo sát trên tấm hợp kim nhôm ($E = 70\text{ GPa}$, $\rho = 2700\text{ kg/m}^3$, $\nu = 0.33$, kích thước $0.6\text{ m} \times 0.4\text{ m} \times 0.002\text{ m}$), lưới $36 \times 24$ phần tử. Kiểm tra khả năng vượt qua rào cản điều kiện biên của Q9 so với Q4 trong các kịch bản biên ngàm, tựa, tự do.
  • Tầng 2 - Khảo sát đa kịch bản phức tạp trên tấm bê tông: Kích thước $1.0\text{ m} \times 1.0\text{ m} \times 0.05\text{ m}$ ($E = 30\text{ GPa}$, $\rho = 2500\text{ kg/m}^3$, $\nu = 0.2$), lưới $10 \times 10$ và $20 \times 20$. Kiểm thử kịch bản đa hư hỏng, hư hỏng nhỏ mức độ $10-20%$, hư hỏng liền kề và hư hỏng cận biên.
  • Tầng 3 - Kiểm chứng trên mô hình sàn BTCT làm việc sau đàn hồi: Tấm sàn BTCT bốn cạnh tựa đơn kích thước $1.5\text{ m} \times 1.5\text{ m} \times 0.06\text{ m}$, cốt thép $\phi 6@100$, mô phỏng phi tuyến trên ANSYS với mô hình vật lý Willam-Warnke cho bê tông và mô hình song tuyến tính (bilinear kinematic hardening) cho cốt thép.

Quy trình nghiên cứu rigorous

1. Xây dựng toán tử phần tử đẳng tham số 9 nút (Q9)

Chuyển vị uốn $w(\xi, \eta)$ trong miền tọa độ tự nhiên $[-1, 1] \times [-1, 1]$ được nội suy qua 9 nút bằng hàm dạng bậc hai $N_i(\xi, \eta)$ ($i = 1, \dots, 9$): $$w(\xi, \eta) = \sum_{i=1}^{9} N_i(\xi, \eta) w_i = {N}^T {w_e}$$

Ma trận chuyển đổi đạo hàm không gian được xác định thông qua ma trận Jacobi $[J]$ kích thước $2 \times 2$: $$[J] = \begin{bmatrix} \frac{\partial x}{\partial \xi} & \frac{\partial y}{\partial \xi} \ \frac{\partial x}{\partial \eta} & \frac{\partial y}{\partial \eta} \end{bmatrix} = \begin{bmatrix} \sum \frac{\partial N_i}{\partial \xi} x_i & \sum \frac{\partial N_i}{\partial \xi} y_i \ \sum \frac{\partial N_i}{\partial \eta} x_i & \sum \frac{\partial N_i}{\partial \eta} y_i \end{bmatrix}$$

Toán tử độ cong đạo hàm bậc hai được xác định chính xác tại 9 điểm nút mà không cần xấp xỉ sai phân hữu hạn: $${\kappa} = \begin{Bmatrix} -\frac{\partial^2 w}{\partial x^2} \ -\frac{\partial^2 w}{\partial y^2} \ -2\frac{\partial^2 w}{\partial x \partial y} \end{Bmatrix} = [B_b^e] {w_e}$$

Năng lượng biến dạng $MSE_e^k$ được tích phân số chính xác bằng phương pháp cầu phương Gauss-Legendre bậc $3 \times 3$: $$MSE_e^k = \frac{1}{2} \sum_{m=1}^{3} \sum_{n=1}^{3} W_m W_n \left( {\kappa^k}^T [D_b] {\kappa^k} |J| \right)_{(\xi_m, \eta_n)}$$

2. Chỉ số năng lượng biến dạng chuẩn hóa (nMSEBI)

Tỷ số năng lượng biến dạng phần tử $F_{k,ij}$ cho mode thứ $k$ tại phần tử $(i,j)$ và chỉ số hư hỏng $\beta_{ij}$ qua $M$ dạng dao động được chọn lọc: $$F_{k,ij} = \frac{MSE_{k,ij}}{MSE_k^{total}}, \quad \beta_{ij} = \frac{\sum_{k=1}^M F_{k,ij}^}{\sum_{k=1}^M F_{k,ij}}$$ Chuẩn hóa thống kê z-score để loại bỏ nhiễu nền: $$nMSEBI_{ij} = \frac{\beta_{ij} - \bar{\beta}}{\sigma_\beta}$$ trong đó dấu $$ ký hiệu trạng thái kết cấu có hư hỏng, $\bar{\beta}$ và $\sigma_\beta$ là giá trị trung bình và độ lệch chuẩn của toàn bộ các phần tử.

3. Bộ 4 chỉ số đánh giá độ tin cậy chẩn đoán (Diagnostic Metrics)

Để định lượng hóa độ chính xác một cách khách quan, luận án thiết lập bộ chỉ số:

  • Chỉ số A (Độ chính xác vùng hư hỏng - True Positive Rate): $$A = \frac{A_{dam}^{predict, in}}{A_{dam}^{real}} \times 100%$$
  • Chỉ số B (Độ chính xác vùng không hư hỏng - True Negative Rate): $$B = \frac{A_{undam}^{predict}}{A_{undam}^{real}} \times 100%$$
  • Chỉ số C (Độ chính xác tổng hợp - Overall Accuracy Index): $$C = w_{dam} A + w_{undam} B, \quad \left(w_{dam} = \frac{A_{dam}^{real}}{A_{gross}}, , w_{undam} = \frac{A_{undam}^{real}}{A_{gross}}\right)$$
  • Chỉ số D (Mức độ chẩn đoán sai vùng hư hỏng - False Alarm Severity Index): $$D = \frac{A_{dam}^{predict, out}}{A_{dam}^{real}} \times 100%$$

Data và phân tích

1. Khảo sát 4 hàm mục tiêu định lượng trong GA

Luận án tiến hành so sánh 4 cấu trúc hàm mục tiêu (Objective Functions - OF) để tìm ra công thức hội tụ tối ưu nhất:

  • $OF_1$ (Dựa vào tần số): $OF_1({\alpha}) = \sum_{i=1}^M \left(\frac{f_i^d - f_i^s({\alpha})}{f_i^d}\right)^2$
  • $OF_2$ (Dựa vào dạng dao động): $OF_2({\alpha}) = \sum_{i=1}^M \left(1 - \sqrt{MAC(\phi_i^d, \phi_i^s({\alpha}))}\right)^2$
  • $OF_3$ (Dựa vào MSE toàn tấm): $OF_3({\alpha}) = \sum_{i=1}^M \sum_{e=1}^{N_e} \left(\frac{MSE_{e,i}^d - MSE_{e,i}^s({\alpha})}{MSE_{e,i}^d}\right)^2$
  • $OF_4$ (Dựa vào MSE vùng cục bộ): $OF_4({\alpha}) = \sum_{i=1}^M \sum_{e \in \Omega_{local}} \left(\frac{MSE_{e,i}^d - MSE_{e,i}^s({\alpha})}{MSE_{e,i}^d}\right)^2$

2. Thuật toán di truyền lặp đa pha (Iterative Genetic Algorithm - I-GA)

Thông số GA tiêu chuẩn: Kích thước quần thể (Population size) = 100 cá thể; Xác suất lai ghép (Crossover probability) $P_c = 0.85$; Xác suất đột biến (Mutation probability) $P_m = 0.01$; Số thế hệ tối đa mỗi pha = 200 thế hệ.

Quy trình loại trừ của I-GA:

  1. Pha 1: Chạy GA với tất cả $k$ phần tử nghi ngờ từ Bước 1. Sau khi kết thúc, xác định phần tử $e_{min}$ có mức độ hư hỏng chẩn đoán nhỏ nhất $\alpha_{min} < 5%$.
  2. Pha $m$: Cố định độ cứng phần tử $e_{min}$ về nguyên vẹn ($\alpha = 0$), giảm số biến tối ưu xuống $k - (m - 1)$, tái khởi tạo quần thể và chạy GA.
  3. Điều kiện dừng: Thuật toán kết thúc khi tất cả các biến còn lại đều có $\alpha_j > 10%$ và giá trị hàm mục tiêu $OF_4 \le 10^{-4}$.

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

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

Các kết quả mô phỏng số và phân tích thực nghiệm đã mang lại 5 phát hiện khoa học mang tính đột phá:

                                    5 PHÁT HIỆN KHOA HỌC THEN CHỐT
[PHÁT HIỆN 1]      [PHÁT HIỆN 2]       [PHÁT HIỆN 3]       [PHÁT HIỆN 4]       [PHÁT HIỆN 5]
Q9 vượt trội Q4:   Lọc mode bằng MAC:  Quy trình 2 bước:   I-GA đa pha:        Bắt vết nứt sớm:
Khử sạch nhiễu     Chỉ số C tăng       C = 96.8%           Định lượng sai số   Định vị nứt tại
biên, C = 98.2%    12-18% so với ngẫu  Tiết kiệm 70%       < 1.2%              q_crack = 3.134
vs Q4 (C = 74.5%)  nhiên các mode đầu  dữ liệu đo đạc      Khử 100% false pos  kN/m2 theo ACI

========================================================================================
BẢNG SO SÁNH HIỆU QUẢ CHẨN ĐOÁN VỊ TRÍ HƯ HỎNG GIỮA PHẦN TỬ Q4 VÀ Q9 (BIÊN TỰ DO - FFFF)
========================================================================================
========================================================================================
  1. Phần tử Q9 loại bỏ triệt để hiện tượng méo tín hiệu tại biên: Trong điều kiện biên tự do (FFFF) và ngàm (CCCC), phương pháp Q4 truyền thống bị nhiễu biên nặng nề dẫn đến chỉ số báo sai $D = 26.2%$ và độ chính xác tổng thể chỉ đạt $C = 74.5%$. Ngược lại, phần tử Q9 nâng chỉ số $B$ lên $98.1%$, kéo chỉ số $D$ xuống mức danh nghĩa $1.9%$, đưa độ chính xác tổng thể $C$ đạt $98.2%$. Kết quả này chứng minh xấp xỉ bậc cao triệt tiêu hoàn toàn sự phụ thuộc vào điều kiện biên.
  2. Cơ chế lọc mode bằng tương quan MAC tăng đột biến độ nhạy chẩn đoán: Việc lựa chọn các mode dao động có hệ số tương quan $MAC \ge 0.92$ (thay vì lấy ngẫu nhiên 4-6 mode đầu tiên) giúp loại bỏ hiện tượng che mờ hư hỏng (damage masking). Đối với kịch bản đa hư hỏng phức tạp với mức độ suy giảm chỉ $10%$, việc lọc mode giúp chỉ số $C$ tăng thêm từ $12.4%$ đến $18.6%$.
  3. Quy trình hai bước Global-Local tiết kiệm 70% dữ liệu cảm biến: Bước tổng thể với ngưỡng $20%$ khoanh đúng $100%$ các vùng nghi ngờ ($A = 100%$). Bước cục bộ với ngưỡng $30%$ thu hẹp chính xác tọa độ phần tử nứt, giảm số lượng biến tối ưu cho Bước 2 từ 100 phần tử xuống còn 4-8 phần tử, giảm thiểu $70%$ khối lượng đo đạc cảm biến trên hiện trường.
  4. Hàm mục tiêu $OF_4$ kết hợp I-GA định lượng chính xác tuyệt đối mức suy giảm độ cứng: So sánh giữa 4 hàm mục tiêu cho thấy $OF_1$ (tần số) và $OF_2$ (dạng dao động) hoàn toàn không hội tụ khi có trên 3 phần tử hư hỏng. $OF_3$ (MSE toàn tấm) hội tụ chậm và sai số trên $15%$. Riêng $OF_4$ (MSE cục bộ) kết hợp I-GA đạt độ hội tụ sai số mức độ hư hỏng $\Delta \alpha < 1.2%$. Toàn bộ các phần tử cảnh báo sai ban đầu đều bị loại bỏ tuần tự về $\alpha = 0%$ sau 2-3 pha lặp.
  5. Khả năng chẩn đoán sớm vết nứt bê tông cốt thép phi tuyến: Trên mô hình tấm sàn BTCT chịu tải trọng phân bố đều, tại cấp tải gây nứt $q_1 = q_{crack} = 3,134\text{ kN/m}^2$ (khi độ võng giữa nhịp mới chỉ đạt $0.42\text{ mm}$), thuật toán đã khoanh vùng chính xác các phần tử nứt tại thớ dưới giữa sàn, hoàn toàn trùng khớp với vùng nứt phần tử SOLID65 trên ANSYS và tính toán theo ACI 318-02. Khi tải trọng tăng lên các cấp $q_2 = 6.0\text{ kN/m}^2$ và $q_3 = 9.0\text{ kN/m}^2$, thuật toán tiếp tục theo dõi chính xác sự lan truyền vết nứt hướng về 4 góc bản.

Implications đa chiều

  • Về mặt Lý thuyết: Đặt nền móng cho việc mở rộng cơ học năng lượng biến dạng sang các kết cấu hai chiều phức tạp có tính dị hướng (anisotropic), composite nhiều lớp hoặc kết cấu vỏ mỏng cong.
  • Về mặt Phương pháp luận: Đưa ra chuẩn mực đánh giá mới cho các thuật toán SHM thông qua bộ chỉ số định lượng $(A, B, C, D)$, chấm dứt việc đánh giá trực quan định tính vốn thiếu tính nhất quán trong các nghiên cứu trước đây.
  • Về mặt Ứng dụng thực tiễn: Quy trình chỉ yêu cầu cảm biến gia tốc đo chuyển vị đứng $w$ (như gia tốc kế đơn trục Piezoelectric hoặc cảm biến LVDT/Laser Doppler), giảm thiểu $80%$ chi phí đầu tư phần cứng so với các hệ thống đo chuyển vị đa trục đắt tiền.
  • Về mặt Tiêu chuẩn & Chính sách: Cung cấp cơ sở khoa học để các cơ quan quản lý nhà nước (Bộ Xây dựng, Bộ Giao thông Vận tải) xây dựng quy trình kiểm định công trình bằng dao động động lực học, hỗ trợ ban hành quy chuẩn bảo trì định kỳ cho các công trình hạ tầng trọng điểm quốc gia.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, tác giả luận án thẳng thắn chỉ ra các giới hạn nghiên cứu (boundary conditions & limitations):

  1. Hình học kết cấu: Nghiên cứu mới tập trung kiểm chứng trên kết cấu tấm chữ nhật phẳng có chiều dày không đổi; chưa mở rộng cho các tấm có hình học phức hợp (tam giác, đa giác, tấm có lỗ khoét lớn) hoặc tấm có sườn tăng cường (stiffened plates).
  2. Điều kiện môi trường: Các mô hình phân tích số giả định môi trường nhiệt độ và độ ẩm không đổi, chưa tích hợp mô hình bù trừ sai lệch do nhiệt độ (Thermal Compensation Model) vào tiêu chuẩn MAC.
  3. Mô hình tải trọng động: Thử nghiệm trên tấm sàn BTCT mới dừng lại ở việc trích xuất dạng dao động tự do sau các cấp gia tải tĩnh, chưa khảo sát ứng xử dao động cưỡng bức dưới tải trọng xe chạy (moving load) hoặc kích động động đất ngẫu nhiên.

Chương trình nghiên cứu tương lai (10-Year Research Agenda):

  • Giai đoạn 1 (1-3 năm): Mở rộng thuật toán Q9 cho kết cấu vỏ cong (shells) và tấm composite nhiều lớp chịu hiện tượng tách lớp (delamination). Tích hợp thuật toán lọc Kalman mở rộng (EKF) để xử lý nhiễu ngẫu nhiên trong điều kiện đo đạc thực địa.
  • Giai đoạn 2 (3-5 năm): Kết hợp phương pháp MSE cục bộ với Mạng nơ-ron tích chập (CNN) và Mạng nơ-ron vật lý (Physics-Informed Neural Networks - PINNs) nhằm tăng tốc độ tối ưu hóa ở Bước 2 lên mức thời gian thực (real-time).
  • Giai đoạn 3 (5-10 năm): Phát triển phần cứng hệ thống cảm biến không dây thông minh (Wireless Smart Sensor Networks - WSSNs) tích hợp chip xử lý biên (Edge Computing) nhúng sẵn thuật toán I-GA để thương mại hóa giải pháp SHM trên các công trình cầu vượt, hầm ngầm và nhà cao tầng.

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

  • Tác động Học thuật: Luận án đã công bố nhiều công trình trên các tạp chí khoa học uy tín thuộc danh mục ISI/Scopus và hội nghị chuyên ngành quốc tế. Phương pháp luận Q9 MSE và I-GA mở ra hướng trích dẫn đột phá trong cộng đồng nghiên cứu cơ học tính toán và SHM quốc tế, ước tính tạo ra tiềm năng trích dẫn lớn trong các nghiên cứu về kết cấu bê tông và tấm vỏ.
  • Chuyển đổi Ngành Kỹ thuật Xây dựng: Cung cấp công cụ chẩn đoán không phá hủy cực kỳ chính xác cho các tổng công ty tư vấn kiểm định, trung tâm giám định chất lượng xây dựng, cho phép phát hiện sớm vết nứt trong sàn cầu, dầm bản mặt cầu bê tông trước khi xuất hiện dấu hiệu suy thoái nghiêm trọng trên bề mặt.
  • Tác động Kinh tế - Xã hội: Giúp chuyển đổi chiến lược bảo trì từ "bảo trì khắc phục" (reactive maintenance - tốn kém và nguy hiểm) sang "bảo trì dự đoán" (predictive maintenance). Tiết kiệm ước tính $30-45%$ ngân sách bảo dưỡng hàng năm cho các đơn vị quản lý vận hành kết cấu hạ tầng giao thông và dân dụng.
  • Tầm ảnh hưởng Toàn cầu: Giải quyết trọn vẹn thách thức loại bỏ chuyển vị xoay nút trong chẩn đoán kết cấu 2D – một vấn đề học thuật kéo dài nhiều năm trên các diễn đàn kết cấu quốc tế như Journal of Sound and Vibration, Mechanical Systems and Signal Processing, Engineering Structures.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận khung lý thuyết toán học hoàn chỉnh về phép biến đổi đẳng tham số trong tính toán năng lượng dao động, cùng kho dữ liệu mở về các kịch bản kiểm chứng mô phỏng số.
  • Giảng viên & Chuyên gia nghiên cứu cao cấp: Ứng dụng mô hình I-GA đa pha và hệ thống chỉ số $(A, B, C, D)$ vào các đề tài nghiên cứu chuyên sâu về giám sát an toàn công trình biển, tháp gió ngoài khơi và kết cấu hàng không.
  • Kỹ sư R&D & Doanh nghiệp Công nghệ SHM: Sở hữu lõi thuật toán tối ưu hóa chuyển vị đứng $w$ để lập trình trực tiếp vào firmware của các thiết bị cảm biến rung động IOT, thương mại hóa sản phẩm máy đo chẩn đoán khuyết tật kết cấu.
  • Cơ quan Quản lý Cơ sở hạ tầng: Có được công cụ đánh giá sức khỏe công trình minh bạch, chính xác để đưa ra quyết định cấp phép tiếp tục khai thác hoặc hạn chế tải trọng đối với các công trình đã khai thác lâu năm.

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à nó mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập công thức tính toán năng lượng biến dạng dao động (MSE) dựa trên phần tử đẳng tham số 9 nút Lagrange (Q9). Nghiên cứu đã mở rộng Lý thuyết Tấm dày Mindlin-Reissner trong cơ học kết cấu bằng cách chứng minh giải tích rằng trường độ cong uốn ${\kappa} = [B_b^e] {w_e}$ hoàn toàn có thể nội suy chính xác bậc hai thông qua ma trận Jacobi $[J]$ chỉ từ 1 thành phần chuyển vị đứng $w$. Điều này xóa bỏ hoàn toàn giả định bắt buộc phải có 5 thành phần chuyển vị nút $(u, v, w, \theta_x, \theta_y)$ của lý thuyết phần tử hữu hạn cổ điển khi tính năng lượng uốn, đồng thời giải phóng phương pháp MSE khỏi sự phụ thuộc vào điều kiện biên của kết cấu tấm.

2. Phương pháp luận của luận án có điểm gì đổi mới khi so sánh với các nghiên cứu quốc tế tiêu biểu?

So với nghiên cứu kinh điển của Shi, Law & Zhang (1998, 2000) (sử dụng chỉ số MSECR dựa trên ma trận độ cứng đòi hỏi chuyển vị xoay) và Cornwell, Doebling & Farrar (1999) (chia dải dầm 1D và dùng sai phân trung tâm):

  • Luận án của Lê Thanh Cao không dùng phép vi phân sai phân trung tâm (vốn gây sai số biên cục bộ lên đến $30-40%$), mà dùng phép tích phân số Gauss $3 \times 3$ trên phần tử Q9, giảm sai số nhận diện biên xuống dưới $1.9%$.
  • Khác với quy trình tối ưu hóa một bước của Friswell và cộng sự (1998) vốn bị bế tắc khi số lượng biến lớn, luận án phát triển thuật toán di truyền lặp đa pha (I-GA) kết hợp hàm mục tiêu MSE cục bộ ($OF_4$), tự động đào thải các phần tử cảnh báo sai để đạt độ chính xác định lượng tổn thất độ cứng gần như tuyệt đối ($\Delta \alpha < 1.2%$).

3. Phát hiện bất ngờ nhất trong quá trình phân tích dữ liệu là gì?

Phát hiện bất ngờ và có giá trị thực tiễn cao nhất là "Hiệu ứng nghịch đảo tương quan mode (Mode Correlation Inversion Effect)": Việc sử dụng các dạng dao động cơ bản đầu tiên (Mode 1, Mode 2) không phải lúc nào cũng đem lại kết quả chẩn đoán tốt nhất, do một số mode bậc thấp có các đường nút (nodal lines) đi trùng qua vị trí vết nứt khiến độ nhạy năng lượng biến dạng bằng 0. Thay vào đó, việc sàng lọc các mode có chỉ số $MAC \ge 0.92$ giữa hai trạng thái đã tự động chọn ra các mode dao động bậc cao hơn nhưng có năng lượng uốn tập trung cực đại tại vùng nứt, giúp phát hiện các vết nứt siêu nhỏ ($10%$) mà phương pháp truyền thống hoàn toàn bỏ sót.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không?

Có, luận án cung cấp chi tiết toàn bộ quy trình tái lập:

  • Thông số hình học và vật liệu: Kích thước tấm, mô đun đàn hồi $E$, hệ số Poisson $\nu$, khối lượng riêng $\rho$, ma trận độ cứng $[D_b]$.
  • Mã nguồn giải thuật: Cung cấp cấu trúc hàm dạng 9 nút $N_i(\xi, \eta)$, tọa độ 9 điểm Gauss, ma trận Jacobi, các hàm mục tiêu $OF_1 - OF_4$, cùng thông số GA ($P_c = 0.85$, $P_m = 0.01$, Population = 100).
  • Quy trình mô phỏng ANSYS: Chi tiết các bước gán phần tử SOLID65, LINK180, mô hình vật liệu bê tông vỡ nứt Willam-Warnke và quy trình tăng tải trọng tĩnh phục vụ trích xuất tần số và mode shape.

5. Chương trình nghiên cứu 10 năm được phác thảo như thế nào?

Chương trình nghiên cứu 10 năm tập trung vào 3 trọng tâm: (1) Mở rộng toán tử Q9 cho kết cấu tấm cong, composite và tích hợp mô hình bù trừ nhiệt độ môi trường; (2) Tích hợp học sâu (PINNs, CNN) để giải bài toán định lượng mức độ hư hỏng thời gian thực thay cho GA; (3) Chế tạo thiết bị cảm biến thông minh không dây (Edge-AI WSSN) ứng dụng trực tiếp thuật toán vào công tác kiểm định các công trình cầu và tháp cao tầng quy mô lớn.

Kết luận

Luận án tiến sĩ của tác giả Lê Thanh Cao đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra, mang lại 6 đóng góp cốt lõi cho chuyên ngành Kỹ thuật Xây dựng và lĩnh vực Giám sát Sức khỏe Kết cấu:

  1. Xây dựng thành công phương pháp xấp xỉ năng lượng biến dạng dao động dựa trên phần tử đẳng tham số 9 nút (Q9), chỉ sử dụng duy nhất thành phần chuyển vị đứng $w$, loại bỏ hoàn toàn ảnh hưởng của điều kiện biên và nâng độ chính xác tổng thể chẩn đoán lên $98.2%$.
  2. Thiết lập tiêu chuẩn chọn lọc mode dao động dựa trên hệ số tương quan MAC, giải quyết triệt để bài toán nhận diện đa hư hỏng, hư hỏng nhỏ ($10%$) và hư hỏng cận biên.
  3. Đề xuất quy trình chẩn đoán 2 bước tối ưu (Global-Local MSE), kết hợp bộ ngưỡng hư hỏng kép ($20%$ cho bước tổng thể và $30%$ cho bước cục bộ), giúp tinh giảm $70%$ khối lượng dữ liệu đo đạc.
  4. Phát triển thuật toán di truyền lặp đa pha (I-GA) với hàm mục tiêu MSE cục bộ ($OF_4$), loại bỏ hoàn toàn các phần tử cảnh báo sai và định lượng chính xác mức suy giảm độ cứng với sai số $\Delta \alpha < 1.2%$.
  5. Đề xuất bộ 4 chỉ số định lượng $(A, B, C, D)$, chuẩn hóa phương pháp đánh giá hiệu năng của các thuật toán SHM trên kết cấu hai chiều.
  6. Kiểm chứng thành công quy trình trên mô hình tấm sàn BTCT chịu uốn phi tuyến, định vị chính xác vị trí xuất hiện vết nứt đầu tiên tại cấp tải gây nứt $q_{crack} = 3,134\text{ kN/m}^2$ theo tiêu chuẩn ACI 318-02 và ANSYS.

Nghiên cứu đánh dấu bước chuyển dịch quan trọng trong chuyên ngành Kỹ thuật Xây dựng: đưa phương pháp năng lượng biến dạng từ các phân tích dầm 1D lý thuyết thành công cụ chẩn đoán thực tế cho kết cấu tấm 2D phức tạp. Kết quả của luận án mở ra 3 hướng nghiên cứu mới về chẩn đoán kết cấu composite, tích hợp trí tuệ nhân tạo (PINNs) và phát triển thiết bị SHM thông minh, đóng góp nền tảng khoa học vững chắc bảo đảm sự an toàn, bền vững cho các công trình hạ tầng trong tương lai.