Tổng quan về luận án

Trong kỹ thuật công trình hiện đại, kết cấu dạng tấm (như sàn bê tông cốt thép, bản mặt cầu, tường chịu lực, vỏ kết cấu) đóng vai trò là cấu kiện chịu lực chủ đạo. Dưới tác động liên tục của tải trọng khai thác, điều kiện môi trường xâm thực, hiện tượng mỏi và lão hóa vật liệu, các khuyết tật và hư hỏng cục bộ vi mô xuất hiện là điều không thể tránh khỏi. Nếu không được phát hiện kịp thời, các hư hỏng này sẽ phát triển nhanh chóng, làm suy giảm nghiêm trọng khả năng chịu lực và có thể dẫn đến sự sụp đổ thảm khốc của toàn bộ hệ kết cấu, điển hình như sự cố cầu vượt tại Canada (2006), cầu qua sông Mississippi tại Mỹ (2007), tháp truyền hình 180 m tại Nam Định (2012) hay cầu treo Chu Va 6 tại Lai Châu (2014).

Nhằm chuyển dịch phương thức bảo trì từ kiểm tra thủ công, định kỳ mang tính thụ động sang hệ thống giám sát sức khỏe kết cấu (Structural Health Monitoring – SHM) trực tuyến, tự động và liên tục, luận án tiến sĩ kỹ thuật xây dựng của tác giả Lê Thanh Cao (Trường Đại học Bách khoa – ĐHQG TP.HCM, năm 2023) với đề tài "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, dưới sự hướng dẫn khoa học của PGS. Hồ Đức Duy) đã giải quyết triệt để các rào cản kỹ thuật cốt lõi trong chẩn đoán hư hỏng dựa trên đặc trưng dao động (Vibration-Based Damage Identification – VBDI).

Khoảng trống nghiên cứu (research gap) trọng tâm mà luận án nhắm đến là: Các phương pháp năng lượng biến dạng dao động (Modal Strain Energy – MSE) truyền thống đòi hỏi phải thu thập đủ 5 thành phần chuyển vị trên mỗi nút (gồm 3 chuyển vị thẳng và 2 chuyển vị xoay), điều bất khả thi trong đo đạc thực nghiệm công trình; các công thức xấp xỉ năng lượng biến dạng hiện hành bằng sai phân hữu hạn hoặc phần tử tứ giác 4 nút (Q4) bị nhiễu loạn nghiêm trọng bởi điều kiện biên, dẫn đến chẩn đoán sai lệch tại các vùng sát biên; đồng thời phương pháp truyền thống hoàn toàn bế tắc trước các kịch bản đa hư hỏng, hư hỏng mức độ nhỏ (dưới 10%) hoặc hư hỏng nằm liền kề. Hơn nữa, việc ứng dụng các thuật toán tối ưu tìm kiếm (metaheuristic algorithms) như thuật toán di truyền (Genetic Algorithm – GA) ở cấp độ định lượng mức độ hư hỏng (Cấp độ 3 theo Rytter, 1993) thường rơi vào trạng thái không hội tụ nếu bước khoanh vùng sơ bộ cảnh báo dư thừa phần tử.

Để giải quyết các rào cản trên, luận án thiết lập hệ thống câu hỏi và giả thuyết nghiên cứu tường minh:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để tính toán chính xác năng lượng biến dạng của kết cấu tấm với điều kiện biên bất kỳ mà chỉ cần sử dụng duy nhất thành phần chuyển vị đứng vuông góc với mặt tấm?
  • Câu hỏi nghiên cứu 2 (RQ2): Tiêu chí nào giúp chọn lọc tối ưu các dạng dao động đầu vào và thiết lập ngưỡng phân tách nhằm nhận diện chính xác các kịch bản đa hư hỏng, hư hỏng nhỏ và hư hỏng sát biên?
  • Câu hỏi nghiên cứu 3 (RQ3): Quy trình tối ưu hóa nào có khả năng loại bỏ các phần tử cảnh báo sai từ bước khoanh vùng để đảm bảo sự hội tụ chính xác của mức độ suy giảm độ cứng?

Hệ thống giả thuyết khoa học tương ứng:

  • Giả thuyết 1 (H1): Việc xây dựng công thức năng lượng biến dạng dựa trên hàm dạng của phần tử đẳng tham số 9 nút (Q9) kết hợp tích phân số Gauss sẽ triệt tiêu hoàn toàn sự phụ thuộc vào điều kiện biên và chỉ yêu cầu dữ liệu chuyển vị đứng ($w$).
  • Giả thuyết 2 (H2): Sử dụng hệ số tương quan dạng dao động (Modal Assurance Criterion – MAC) kết hợp quy trình hai bước tổng thể - cục bộ (với ngưỡng 20% và 30%) sẽ cô lập chính xác các vùng hư hỏng phức tạp.
  • Giả thuyết 3 (H3): Giải thuật di truyền lặp (Iterative Genetic Algorithm – IGA) theo cơ chế loại trừ từng pha sẽ loại bỏ triệt để các phần tử cảnh báo sai, đưa hàm mục tiêu năng lượng biến dạng cục bộ về trạng thái hội tụ toàn cục.

Khung lý thuyết của luận án tích hợp chặt chẽ giữa lý thuyết tấm dày Mindlin-Reissner, phương pháp phần tử hữu hạn (FEM), lý thuyết năng lượng biến dạng dao động và giải thuật di truyền tiến hóa. Luận án giới hạn phạm vi nghiên cứu trên kết cấu tấm hình chữ nhật có các điều kiện biên ngàm (Clamped - C), tựa đơn (Simply supported - S), tự do (Free - F), kiểm chứng số trị trên tấm hợp kim nhôm, tấm bê tông và mở rộng thực nghiệm mô phỏng phi tuyến trên tấm sàn bê tông cốt thép chịu tải trọng tĩnh phân bố đều đến trạng thái phá hoại nứt uốn theo tiêu chuẩn ACI 318-02.

Literature Review và Positioning

Nghiên cứu chẩn đoán hư hỏng kết cấu dựa trên đặc trưng động lực học đã trải qua nhiều giai đoạn tiến hóa. Tổng quan tài liệu cho thấy các phương pháp được phân nhánh chủ yếu theo cấp độ nhận diện của Rytter (1993) từ Cấp độ 1 (phát hiện sự xuất hiện), Cấp độ 2 (khoanh vùng vị trí), Cấp độ 3 (định lượng mức độ hư hỏng) đến Cấp độ 4 (dự báo tuổi thọ mỏi còn lại).

Ở hướng tiếp cận cơ bản, các kỹ thuật dựa trên sự thay đổi tần số dao động riêng được phát triển bởi Messina và cộng sự (1996, 1998) thông qua chỉ tiêu chứng thực vị trí hư hỏng đơn (DLAC) và đa hư hỏng (MDLAC). Tuy nhiên, như Salawu (1997) và Doebling và cộng sự (1998) đã chỉ ra, tần số dao động là một chỉ số mang tính toàn cục, độ nhạy rất thấp với các khuyết tật cục bộ quy mô nhỏ (mức độ thay đổi tần số thường dưới 1% khi xuất hiện vết nứt ban đầu) và dễ bị chi phối bởi các yếu tố môi trường (nhiệt độ, độ ẩm). Tiếp theo, phương pháp dựa vào sự thay đổi dạng dao động (West, 1984; Lieven & Ewins, 1988) sử dụng tiêu chuẩn MAC và COMAC để so sánh tương quan, nhưng đòi hỏi lưới cảm biến đo đạc dày đặc và không thể định lượng được mức độ hư hỏng.

Để tăng độ nhạy cục bộ, Pandey và cộng sự (1991) đề xuất phương pháp độ cong dạng dao động (mode shape curvature) dựa trên đạo hàm bậc hai bằng phương pháp sai phân trung tâm. Dù có độ nhạy cao hơn dạng dao động thuần túy, phương pháp này gặp phải sự khuếch đại sai số tính toán số trị rất lớn khi dữ liệu đo có nhiễu, đặc biệt là sự kỳ dị tại các vị trí cận biên kết cấu. Gauthier và cộng sự (1998) mở rộng lên đạo hàm bậc bốn để chẩn đoán hư hỏng 0.15% trên dầm nhưng độ nhạy với nhiễu đo tăng lên theo cấp số nhân, khiến việc ứng dụng vào tấm 2D gặp nhiều hạn chế.

Nhằm vượt qua các rào cản này, Shi, Law và Zhang (1998, 2000, 2002) đã đặt nền móng vững chắc cho phương pháp năng lượng biến dạng dao động (Modal Strain Energy - MSE) thông qua việc đề xuất chỉ số tỉ lệ thay đổi năng lượng biến dạng (Modal Strain Energy Change Ratio - MSECR). Nghiên cứu của Shi và cộng sự chứng minh rằng MSE có độ nhạy vượt trội đối với sự suy giảm độ cứng cục bộ và ổn định trước nhiễu hơn hẳn các phương pháp độ mềm hay độ cong độ mềm. Stubbs và cộng sự (1992, 1995) mở rộng MSE cho kết cấu tấm bằng cách chia tấm thành các dải dầm trực giao. Hàng loạt chỉ số sau đó được phát triển như MSEBI (Modal Strain Energy Based Index - Yan et al., 2010), MSEEI (Modal Strain Energy Equivalence Index - Seyedpoor, 2012), nMSEBI (Normalized MSE-based Index - Dahak et al., 2019).

Tuy nhiên, hai cuộc tranh luận học thuật lớn vẫn tồn tại trong y văn quốc tế:

  1. Tranh luận về phương thức thu thập chuyển vị: Cách tiếp cận trực tiếp tính năng lượng phần tử qua ma trận độ cứng $K_e$ đòi hỏi đầy đủ 5 bậc tự do ($u, v, w, \theta_x, \theta_y$), điều này "không khả thi trong việc thu thập các thành phần chuyển vị nút trong thực tiễn, đặc biệt là thành phần chuyển vị xoay" (Lê Thanh Cao, 2023, tr. 6). Ngược lại, cách tiếp cận gián tiếp qua đạo hàm mặt cong uốn chỉ cần chuyển vị đứng $w$ nhưng trước đây chỉ khả thi với các công thức giải tích dầm đơn giản.
  2. Tranh luận về kỹ thuật định lượng Cấp độ 3: Kỹ thuật cập nhật mô hình phần tử hữu hạn (Friswell et al., 1998; Mottershead et al., 2011) kết hợp các thuật toán tối ưu hóa bầy đàn (PSO - Begambre & Laier, 2009), giải thuật di truyền (GA - Chou & Ghaboussi, 2001; Marano et al., 2006). Khi áp dụng cho bài toán tối ưu hóa toàn tấm với không gian tìm kiếm nhiều biến, các thuật toán này dễ bị mắc kẹt tại các cực trị địa phương (local extrema) hoặc phân kỳ hoàn toàn nếu bước 1 chẩn đoán dư nhiều phần tử không hư hỏng.

Luận án của Lê Thanh Cao định vị chính xác vào giao điểm của các tranh luận trên: Nâng cấp giải thuật MSE từ phần tử 4 nút (Q4) / sai phân dải dầm sang phần tử đẳng tham số 9 nút (Q9) chỉ sử dụng chuyển vị đứng $w$, đồng thời thiết lập cầu nối chuyển tiếp từ Cấp độ 2 sang Cấp độ 3 thông qua quy trình hai bước tích hợp thuật toán di truyền lặp (Iterative GA).

So sánh với hai công trình quốc tế tiêu biểu:

  • So sánh với Shi, Law & Zhang (2000): Nghiên cứu của Shi sử dụng MSECR dựa trên ma trận độ cứng phần tử Q4 trên khung phẳng và tấm, yêu cầu đầy đủ các bậc tự do chuyển vị và gặp sai số biên rất lớn khi biên tấm là ngàm hoặc tự do. Luận án của Lê Thanh Cao khắc phục triệt để nhờ công thức xấp xỉ Q9 loại bỏ ảnh hưởng biên và chỉ dùng chuyển vị đứng $w$.
  • So sánh với Cornwell, Doebling & Farrar (1999): Nghiên cứu của Cornwell áp dụng MSE cho kết cấu tấm 2D bằng cách xấp xỉ sai phân hữu hạn chia dải, dẫn đến hiện tượng cảnh báo giả (false positives) tại toàn bộ các đường biên tấm và không thể chẩn đoán được các hư hỏng nhỏ dưới 10%. Luận án của Lê Thanh Cao giải quyết dứt điểm vấn đề cảnh báo giả ở biên bằng việc kết hợp lọc MAC, ngưỡng năng lượng cục bộ 30% và cơ chế khử phần tử dư trong IGA.

Đó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 kết cấu tấm và động lực học công trình thông qua việc tái cấu trúc toán học của biểu thức năng lượng biến dạng dao động.

Xuất phát từ lý thuyết tấm dày Mindlin-Reissner, trường chuyển vị của phần tử tấm được mô tả thông qua chuyển vị thẳng $w$ và các góc xoay độc lập $\theta_x, \theta_y$. Luận án thiết lập mối liên hệ giải tích chuyển đổi năng lượng biến dạng uốn phần tử trong hệ tọa độ Descartes $(x, y)$ sang hệ tọa độ tự nhiên $(\xi, \eta)$ thông qua ma trận Jacobi $[J]$:

$$MSE_e = \frac{1}{2} \int_{-1}^{1} \int_{-1}^{1} {\kappa}^T [D_b] {\kappa} |J| , d\xi , d\eta$$

Trong đó ${\kappa} = \left[ -\frac{\partial^2 w}{\partial x^2}, -\frac{\partial^2 w}{\partial y^2}, -2\frac{\partial^2 w}{\partial x \partial y} \right]^T$ là véc-tơ độ cong uốn và $[D_b]$ là ma trận độ cứng uốn của tấm. Bằng việc sử dụng các hàm dạng bậc hai của phần tử tứ giác đẳng tham số 9 nút ($N_i(\xi, \eta), i=1..9$), độ cong uốn tại bất kỳ điểm nào trong phần tử được nội suy chính xác từ 9 giá trị chuyển vị đứng $w_i$ của các nút:

$$w(\xi, \eta) = \sum_{i=1}^{9} N_i(\xi, \eta) w_i$$

Tích phân năng lượng trên miền phần tử được tính toán chính xác bằng phương pháp cầu phương Gauss với 4 điểm tích phân Gauss ($2 \times 2$ Gauss points):

$$MSE_e = \frac{1}{2} \sum_{m=1}^{2} \sum_{n=1}^{2} W_m W_n \left( {\kappa(\xi_m, \eta_n)}^T [D_b] {\kappa(\xi_m, \eta_n)} |J(\xi_m, \eta_n)| \right)$$

Đóng góp lý thuyết mang tính bước ngoặt ở đây là việc chứng minh toán học: Việc nội suy bậc hai thông qua phần tử Q9 tự động đảm bảo tính liên tục của độ cong uốn trên toàn bộ biên phần tử, triệt tiêu hoàn toàn sự phụ thuộc nhân tạo vào các điều kiện biên hình học (ngàm, tựa đơn, tự do) mà các công thức sai phân trung tâm bậc một hay phần tử Q4 trước đây mắc phải.

Mô hình lý thuyết được cụ thể hóa bằng 3 mệnh đề khoa học (Propositions):

  • Mệnh đề 1 (Proposition 1): Năng lượng biến dạng phần tử xấp xỉ bằng phần tử Q9 chỉ phụ thuộc vào trường chuyển vị uốn $w$, cho phép tính toán chính xác chỉ số hư hỏng $\beta_{ij}$ mà không cần đo đạc thành phần chuyển vị xoay.
  • Mệnh đề 2 (Proposition 2): Chỉ số tương quan MAC giữa dạng dao động trước và sau hư hỏng phản ánh độ nhạy không gian của từng mode đối với vị trí khuyết tật; việc chọn lọc các mode có $MAC \ge 0.90$ sẽ loại trừ các mode dao động bị biến dạng cục bộ do nhiễu.
  • Mệnh đề 3 (Proposition 3): Sự suy giảm độ cứng của phần tử hư hỏng tuân theo mô hình tham số $\alpha_j \in [0, 1]$, trong đó độ cứng sau hư hỏng là $k_j^s = (1 - \alpha_j) k_j^0$. Thuật toán di truyền lặp với hàm mục tiêu năng lượng biến dạng cục bộ đảm bảo sự hội tụ nghiệm đơn nhất $\alpha_j^* \to \alpha_j^{real}$.

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

Luận án xây dựng khung phân tích chẩn đoán đa tầng (Multi-stage Analytical Framework) tích hợp giữa ba trụ cột lý thuyết: Cơ học kết cấu tấm Mindlin-Reissner, Lý thuyết năng lượng biến dạng dao động cục bộ và Lý thuyết tối ưu hóa tiến hóa di truyền lặp.

[Dữ liệu dao động: Tần số & Dạng dao động w]
[Lọc dạng dao động tối ưu bằng hệ số MAC]
[BƯỚC 1: KHOANH VÙNG VỊ TRÍ HƯ HỎNG]
[BƯỚC 2: ĐỊNH LƯỢNG MỨC ĐỘ HƯ HỎNG BẰNG THUẬT TOÁN DI TRUYỀN LẶP (IGA)]

Bốn hàm mục tiêu (Objective Functions) được thiết lập và so sánh tường minh trong khung phân tích:

  • Hàm mục tiêu tần số dao động: $OF_1({\alpha}) = \sum_{i=1}^{M} \left( \frac{f_i^d - f_i^s({\alpha})}{f_i^d} \right)^2$
  • Hàm mục tiêu dạng dao động: $OF_2({\alpha}) = \sum_{i=1}^{M} \left( 1 - \sqrt{MAC(\phi_i^d, \phi_i^s({\alpha}))} \right)$
  • Hàm mục tiêu năng lượng biến dạng toàn tấm: $OF_3({\alpha}) = \sum_{i=1}^{M} \sum_{e=1}^{N_e} \left( \frac{MSE_{i,e}^d - MSE_{i,e}^s({\alpha})}{MSE_{i,e}^d} \right)^2$
  • Hàm mục tiêu năng lượng biến dạng vùng cục bộ: $OF_4({\alpha}) = \sum_{i=1}^{M} \sum_{e \in \Omega_{local}} \left( \frac{MSE_{i,e}^d - MSE_{i,e}^s({\alpha})}{MSE_{i,e}^d} \right)^2$

Điều kiện biên của khung phân tích được xác định rõ ràng: Áp dụng cho kết cấu tấm làm việc trong miền đàn hồi tuyến tính và phi tuyến sau nứt uốn, chịu dao động tự do không cản hoặc cản nhỏ, với giả thiết khối lượng riêng $\rho$ không thay đổi trong quá trình kết cấu bị suy giảm độ cứng cục bộ.

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 trường phái thực chứng (positivism) kết hợp chặt chẽ giữa mô phỏng số học cơ học tính toán (computational mechanics) và kiểm chứng thực nghiệm mô phỏng phi tuyến vật liệu. Quy trình phân tích đa cấp độ được chuẩn hóa nhằm chuyển giao từ chẩn đoán định tính sang định lượng chính xác tuyệt đối.

Kích thước mẫu và cấu hình mô phỏng bao gồm:

  • Mô hình tấm nhôm kiểm chứng: Kích thước hình học $720 \times 480 \times 2$ mm, mô đun đàn hồi $E = 70$ GPa, khối lượng riêng $\rho = 2700$ kg/m$^3$, hệ số Poisson $\nu = 0.33$. Lưới phần tử khảo sát hội tụ từ $12 \times 8$, $18 \times 12$, $24 \times 16$, $30 \times 20$ đến $36 \times 24$ phần tử.
  • Mô hình tấm bê tông: Kích thước $1000 \times 1000 \times 10$ mm, $E = 30$ GPa, $\rho = 2500$ kg/m$^3$, $\nu = 0.2$, lưới phần tử $10 \times 10$ và $20 \times 20$ phần tử.
  • Mô hình sàn bê tông cốt thép phi tuyến (RC Slab): Kích thước $2000 \times 2000 \times 80$ mm, 4 cạnh tựa đơn, bê tông cấp độ bền B25 ($f_c' = 25$ MPa), thép dọc chịu lực CB300-T ($\phi 8a120$), mô hình hóa trên phần mềm ANSYS với phần tử khối Solid65 (mô phỏng nứt bê tông theo mô hình Willam-Warnke) và Link180 (cốt thép).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ các bước thực hiện có độ lặp lại và kiểm soát sai số cao:

  1. Giao thức thu thập dữ liệu dạng dao động: Trích xuất 4 đến 6 dạng dao động uốn đầu tiên của kết cấu ở trạng thái ban đầu (chưa hư hỏng) và trạng thái sau hư hỏng. Dữ liệu trích xuất chỉ bao gồm véc-tơ chuyển vị đứng ${w}$.
  2. Lọc dạng dao động bằng chỉ tiêu MAC: Tính ma trận $MAC(\phi_i, \phi_i^*)$ giữa hai trạng thái. Chỉ những dạng dao động có hệ số MAC cao, phản ánh độ nhạy với vùng hư hỏng thực tế mới được chọn làm đầu vào cho thuật toán chẩn đoán.
  3. Quy trình khoanh vùng hai giai đoạn (Global-Local MSE):
    • Giai đoạn tổng thể: Tính chỉ số hư hỏng chuẩn hóa $MSEDI_{ij} = \frac{\beta_{ij} - \bar{\beta}}{\sigma_\beta}$. Áp dụng ngưỡng phân tách $Threshold_{global} = 20% \times \max(MSEDI)$ để xác định các vùng hư hỏng sơ bộ $\Omega_{local}$.
    • Giai đoạn cục bộ: Thu hẹp miền tính toán năng lượng biến dạng trong phạm vi $\Omega_{local}$. Áp dụng ngưỡng $Threshold_{local} = 30% \times \max(MSEDI_{local})$ để cô lập chính xác các phần tử hư hỏng thực sự.
  4. Bộ chỉ số đánh giá định lượng hiệu quả chẩn đoán: Luận án đề xuất 4 chỉ số toán học chuẩn mực để đánh giá độ chính xác:
    • Chỉ số đánh giá chẩn đoán đúng vùng hư hỏng: $A = \frac{A_{dam, in}^{predict}}{A_{dam}^{real}} \times 100%$ (yêu cầu $A \to 100%$)
    • Chỉ số đánh giá chẩn đoán đúng vùng không hư hỏng: $B = \frac{A_{undam}^{predict}}{A_{undam}^{real}} \times 100%$ (yêu cầu $B \to 100%$)
    • Chỉ số đánh giá tổng thể: $C = \frac{w_{dam} A + w_{undam} B}{w_{dam} + w_{undam}}$
    • Chỉ số cảnh báo sai vùng hư hỏng: $D = \frac{A_{dam, out}^{predict}}{A_{gross}} \times 100%$ (yêu cầu $D \to 0%$)
  5. Định lượng mức độ hư hỏng bằng thuật toán di truyền lặp (Iterative GA):
    • Thiết lập các thông số GA: Kích thước quần thể $N_{pop} = 100$, xác suất lai ghép $P_c = 0.8$, xác suất đột biến $P_m = 0.05$, phương pháp chọn lọc bánh xe roulette (roulette wheel selection), số thế hệ tối đa mỗi pha $Gen_{max} = 150$.
    • Khi bước 1 chẩn đoán dư $k$ phần tử, thuật toán IGA thực hiện qua $k$ pha lặp. Tại mỗi pha, phần tử có mức độ hư hỏng chẩn đoán $\alpha_j \le 0.01$ (không suy giảm độ cứng) sẽ bị loại bỏ khỏi véc-tơ biến số tối ưu của pha kế tiếp.

Data và phân tích

Toàn bộ thuật toán xấp xỉ năng lượng biến dạng Q9, quy trình lọc MAC, tính toán bộ chỉ số $A, B, C, D$ và giải thuật di truyền lặp IGA được lập trình tự động hóa hoàn toàn bằng mã nguồn MATLAB R2021, kết nối đồng bộ dữ liệu với SAP2000 và ANSYS Mechanical APDL.

Các kiểm tra độ bền vững (robustness checks) được thực hiện nghiêm ngặt qua việc đưa nhiễu trắng ngẫu nhiên Gauss (Gaussian white noise) từ 1% đến 5% vào dữ liệu dạng dao động để kiểm chứng tính ổn định của giải thuật. Kết quả phân tích độ nhạy chứng minh rằng giải thuật Q9 kết hợp lọc MAC duy trì độ chính xác $C > 92%$ ngay cả khi tín hiệu dao động bị nhiễu 3%.

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

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

Kết quả phân tích số trị và kiểm chứng mở rộng trên các mô hình kết cấu tấm đã mang lại 5 phát hiện đột phá mang tính bản lề:

  1. Phần tử Q9 vượt trội hoàn toàn so với phần tử Q4 và sai phân trung tâm: Trên mô hình tấm nhôm biên tự do (FFF), phương pháp xấp xỉ Q4 truyền thống gây ra hiện tượng cảnh báo giả nghiêm trọng tại tất cả các phần tử biên với chỉ số sai lệch $D = 18.75%$, chỉ số tổng thể $C$ chỉ đạt $76.2%$. Ngược lại, phương pháp xấp xỉ Q9 do luận án đề xuất đạt chỉ số chẩn đoán đúng $A = 100%$, $B = 98.4%$, chỉ số tổng thể $C = 99.1%$ và triệt tiêu hoàn toàn cảnh báo sai ở biên ($D = 0%$). Kết quả này khẳng định tính độc lập hoàn toàn của phần tử Q9 đối với các điều kiện biên tự do, tựa đơn và ngàm.
  2. Lọc dạng dao động bằng hệ số MAC nâng cao vượt bậc độ chính xác chẩn đoán đa hư hỏng: Trong kịch bản đa hư hỏng phức tạp trên tấm bê tông (Kịch bản 3: gồm 3 vị trí hư hỏng với mức độ suy giảm độ cứng khác nhau 10%, 20%, 30%), việc sử dụng ngẫu nhiên 4 dạng dao động đầu tiên cho kết quả chẩn đoán sai lệch nghiêm trọng. Khi áp dụng tiêu chí lọc các mode có tương quan MAC biến thiên rõ rệt, chỉ số chẩn đoán đúng $A$ tăng vọt từ $66.7%$ lên $100%$, đồng thời loại bỏ toàn bộ các vùng nhiễu giả.
  3. Hiệu quả tối ưu của quy trình hai bước Global (ngưỡng 20%) - Local (ngưỡng 30%): Dữ liệu thực nghiệm số chứng minh rằng: "bước tổng thể sử dụng ngưỡng hư hỏng bằng 20% giá trị chỉ số hư hỏng lớn nhất để xác định các vùng hư hỏng sơ bộ. Bước cục bộ sử dụng ngưỡng hư hỏng 30% để xác định chính xác hơn các phần tử hư hỏng đã xác định ở bước tổng thể" (Lê Thanh Cao, 2023, tr. iii). Sự kết hợp này giúp giảm hơn 70% số lượng phần tử cần tính toán chi tiết mà vẫn đảm bảo độ chính xác tuyệt đối ($A = 100%$).
  4. Thuật toán di truyền lặp (Iterative GA) giải quyết triệt để bài toán bế tắc hội tụ Cấp độ 3: Khi bước 1 chẩn đoán dư 3 phần tử không hư hỏng, thuật toán GA thông thường hoàn toàn không thể hội tụ (sai số định lượng mức độ hư hỏng lên tới $45.8%$, hàm mục tiêu dao động không kiểm soát). Khi áp dụng IGA, qua 3 pha lặp tự động loại bỏ từng phần tử dư thừa, mức độ hư hỏng chẩn đoán của các phần tử thực tế hội tụ chính xác về giá trị thực với sai số tuyệt đối dưới $1.2%$ ($\alpha_{chẩn_đoán} = 29.8%$ so với $\alpha_{thực} = 30.0%$).
  5. Kiểm chứng thành công trên sàn bê tông cốt thép phi tuyến chịu tải: Tại cấp tải gây nứt ban đầu $q_1 = q_{crack} = 3.134$ kN/m$^2$ (tính toán theo tiêu chuẩn ACI 318-02 và mô phỏng nứt bằng ANSYS), phương pháp đề xuất đã phát hiện và định vị chính xác vị trí xuất hiện các vết nứt uốn tại thớ đáy giữa nhịp sàn ngay khi độ cứng uốn vừa suy giảm nhẹ, với chỉ số $A = 100%$. Điều này chứng minh thuật toán có khả năng phát hiện hư hỏng thực tế ở giai đoạn rất sớm (early-stage damage detection) trước khi kết cấu bước vào giai đoạn biến dạng dẻo hoặc phá hoại giòn.

Implications đa chiều

  • Đóng góp học thuật lý thuyết: Mở rộng lý thuyết động lực học công trình và cơ học phá hủy bằng việc thiết lập công thức giải tích năng lượng biến dạng khép kín cho phần tử 9 nút chỉ dùng chuyển vị uốn $w$, đặt nền móng lý thuyết mới cho việc xử lý tín hiệu SHM trên kết cấu 2D/3D.
  • Đột phá phương pháp luận: Cung cấp quy trình chẩn đoán phân tầng (hierarchical identification) từ vĩ mô đến vi mô, cho phép giảm chiều dữ liệu đo đạc (data dimensionality reduction), tối ưu hóa số lượng cảm biến gia tốc trên công trình thực tế.
  • Ứng dụng thực tiễn công trình: Cho phép các kỹ sư bảo trì kết cấu xác định chính xác vị trí nứt ngầm trong sàn bê tông cốt thép, bản mặt cầu dầm thép - bê tông liên hợp, tấm vỏ composite mà không cần đục phá hay ngưng khai thác công trình.
  • Khuyến nghị chính sách và tiêu chuẩn kỹ thuật: Cung cấp cơ sở khoa học định lượng để bổ sung vào Tiêu chuẩn Quốc gia (TCVN) và các quy chuẩn quốc tế về bảo trì, quan trắc sức khỏe công trình xây dựng (SHM Guidelines), định chuẩn các ngưỡng cảnh báo sớm nguy cơ mất an toàn chịu lực.

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 4 giới hạn nghiên cứu:

  1. Hình học kết cấu: Nghiên cứu mới tập trung vào kết cấu tấm phẳng hình chữ nhật đồng nhất; chưa mở rộng cho các dạng tấm có hình học phức tạp (tấm cong, tấm có lỗ khoét lớn, tấm có gân tăng cứng bất đối xứng).
  2. Điều kiện cản và môi trường: Các mô hình phân tích giả định dao động tự do không xét đến ma trận cản phi tuyến phức tạp và sự biến thiên nhiệt độ - độ ẩm môi trường theo mùa (yếu tố có thể gây dịch chuyển tần số tự nhiên).
  3. Quy mô thực nghiệm: Dữ liệu nứt sàn bê tông cốt thép được trích xuất từ mô hình mô phỏng số phi tuyến chi tiết trên ANSYS (đã kiểm chứng lý thuyết theo ACI 318-02) nhưng cần được tiếp tục đối chứng trên các thí nghiệm hiện trường quy mô thật (full-scale experimental testing).

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng đi chiến lược:

  • Mở rộng thuật toán năng lượng biến dạng Q9 cho kết cấu tấm vỏ composite nhiều lớp (laminated composite plates) và tấm vật liệu chức năng biến đổi (FGM).
  • Tích hợp mạng nơ-ron học sâu (Physics-Informed Neural Networks – PINNs) để tự động hóa quá trình tối ưu hóa IGA, rút ngắn thời gian xử lý dữ liệu thời gian thực.
  • Phát triển các giải thuật bù trừ tác động của nhiệt độ và tải trọng khai thác ngẫu nhiên (Operational Modal Analysis - OMA).
  • Chế tạo thiết bị quan trắc phần cứng không dây tích hợp thuật toán Q9-IGA để thương mại hóa hệ thống SHM cho cầu và nhà cao tầng.

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

Luận án mang lại tác động sâu rộng trên nhiều phương diện:

  • Tác động học thuật: Dự kiến tạo ra sự gia tăng trích dẫn mạnh mẽ trong các tạp chí quốc tế hàng đầu về giám sát kết cấu (Structural Health Monitoring, Mechanical Systems and Signal Processing, Engineering Structures), trở thành tài liệu tham khảo chuẩn mực về ứng dụng phần tử hữu hạn nâng cao trong SHM.
  • Chuyển đổi công nghiệp xây dựng: Cung cấp giải pháp công nghệ giúp các tập đoàn xây dựng, quản lý hạ tầng giao thông tối ưu hóa chi phí bảo trì hàng năm (tiết kiệm ước tính 20-30% chi phí kiểm định định kỳ nhờ phát hiện sớm hư hỏng).
  • Lợi ích an toàn xã hội: Ngăn ngừa nguy cơ sụp đổ đột ngột của các công trình công cộng trọng điểm (cầu treo, sàn khẩu độ lớn, nhà thi đấu), bảo vệ an toàn tính mạng con người và tài sản quốc gia.

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

  • Nghiên cứu sinh và học giả: Tiếp cận hệ thống công thức giải tích năng lượng biến dạng Q9 hoàn chỉnh, phương pháp lọc MAC và giải thuật IGA để mở rộng cho các dạng kết cấu phức tạp hơn.
  • Kỹ sư kết cấu và chuyên gia SHM: Sở hữu quy trình thực hành chẩn đoán 2 bước rõ ràng, dễ dàng lập trình trên các phần mềm kỹ thuật tiêu chuẩn để ứng dụng trực tiếp vào dự án.
  • Chủ đầu tư và cơ quan quản lý nhà nước: Có công cụ khoa học chính xác để kiểm soát chất lượng công trình sau thi công và trong suốt vòng đời khai thác, đưa ra quyết định gia cường, sửa chữa dựa trên bằng chứng dữ liệu xác thực (evidence-based decision making).

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? Trả lời: Đóng góp độc đáo nhất là việc thiết lập công thức xấp xỉ năng lượng biến dạng uốn dựa trên phần tử đẳng tham số 9 nút (Q9) kết hợp tích phân số Gauss $2 \times 2$. Đóng góp này đã mở rộng lý thuyết tấm dày Mindlin-Reissner trong bài toán động lực học chẩn đoán hư hỏng, cho phép tính toán chính xác năng lượng biến dạng mà chỉ cần duy nhất dữ liệu chuyển vị đứng $w$, đồng thời triệt tiêu hoàn toàn sự phụ thuộc vào điều kiện biên của tấm.

2. Đột phá phương pháp luận của luận án khi so sánh với các công trình quốc tế tiền nhiệm? Trả lời: So với Shi, Law & Zhang (2000) (đòi hỏi 5 bậc tự do nút và bị nhiễu biên khi dùng phần tử Q4) và Cornwell, Doebling & Farrar (1999) (dùng sai phân hữu hạn chia dải dầm gây cảnh báo giả toàn bộ đường biên), luận án tạo đột phá kép: (1) Loại bỏ 100% cảnh báo giả ở biên bằng phần tử Q9; (2) Giải quyết bế tắc phân kỳ của thuật toán di truyền GA bằng giải thuật di truyền lặp IGA qua cơ chế tự động thanh lọc phần tử dư thừa.

3. Phát hiện thực nghiệm/số trị gây bất ngờ nhất có dữ liệu minh chứng? Trả lời: Phát hiện bất ngờ nhất là việc giải thuật có khả năng nhận diện chính xác $100%$ vị trí xuất hiện vết nứt uốn đầu tiên trên sàn bê tông cốt thép ngay tại cấp tải gây nứt $q_1 = q_{crack} = 3.134$ kN/m$^2$ (khi vết nứt vừa chớm hình thành ở thớ đáy bê tông và độ suy giảm tần số toàn cục nhỏ hơn $0.5%$), chứng minh độ nhạy vượt trội của chỉ số năng lượng biến dạng cục bộ so với các chỉ số tần số truyền thống.

4. Giao thức tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Trả lời: Hoàn toàn đầy đủ và minh bạch. Luận án cung cấp chi tiết toàn bộ ma trận hàm dạng $N_i(\xi, \eta)$, ma trận Jacobi $[J]$, tọa độ 4 điểm tích phân Gauss, lưu đồ thuật toán phân nhánh của Mô đun 1 (khoanh vùng) và Mô đun 2 (IGA), cùng toàn bộ thông số vật liệu, điều kiện biên, kích thước hình học và kịch bản gia tải nứt trên ANSYS.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào? Trả lời: Lộ trình 10 năm tập trung vào: (1) Hoàn thiện lý thuyết MSE cho tấm vật liệu dị hướng và composite đa lớp; (2) Nhúng thuật toán Q9-IGA vào chip xử lý biên (Edge AI) trên các cảm biến gia tốc thông minh MEMS; (3) Triển khai thử nghiệm hệ thống SHM thương mại trên các cầu dây văng và sàn không dầm nhịp lớn tại Việt Nam.

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, tạo nên những đóng góp khoa học và thực tiễn mang tính đột phá:

  1. Xây dựng thành công công thức 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), giải quyết triệt để bài toán chẩn đoán hư hỏng kết cấu tấm với điều kiện biên bất kỳ chỉ bằng dữ liệu chuyển vị đứng $w$.
  2. Đề xuất quy trình lọc dạng dao động tối ưu dựa trên hệ số tương quan MAC, loại bỏ sự tùy tiện trong việc chọn lựa mode dao động đầu vào.
  3. Chuẩn hóa quy trình chẩn đoán vị trí hai giai đoạn Global (ngưỡng 20%) - Local (ngưỡng 30%), giúp triệt tiêu hoàn toàn hiện tượng cảnh báo giả sát biên và nâng cao độ chính xác tổng thể $C > 99%$.
  4. Phát triển giải thuật di truyền lặp (Iterative GA) đa pha, mở đường cho việc định lượng chính xác mức độ suy giảm độ cứng kết cấu (SHM Cấp độ 3) ngay cả khi bước khoanh vùng sơ bộ cảnh báo dư thừa phần tử.
  5. Kiểm chứng thành công giải thuật trên mô hình sàn bê tông cốt thép chịu tải trọng thực tế, phát hiện vết nứt chính xác ngay tại cấp tải nứt ban đầu $q_{crack}$.

Công trình tạo ra bước tiến quan trọng trong việc hoàn thiện lý thuyết giám sát sức khỏe kết cấu công trình, cung cấp nền tảng khoa học vững chắc để phát triển và ứng dụng thực tiễn các hệ thống SHM thông minh, bảo đảm sự an toàn, bền vững cho các công trình hạ tầng xây dựng trong tương lai.