Tổng quan về luận án

Trong bối cảnh công nghiệp hóa và đô thị hóa nhanh chóng, ngành xây dựng Việt Nam đang đối mặt với thách thức kép: sự cạn kiệt tài nguyên thiên nhiên (cát sông, đá tự nhiên) do khai thác quá mức và áp lực môi trường nặng nề từ chất thải rắn công nghiệp. Điển hình tại khu vực miền Trung, "theo báo cáo của Nhà máy nhiệt điện Vũng Áng - Hà Tĩnh, mỗi năm nhà máy phát sinh khoảng 1 triệu tấn tro và xỉ, và tại các mỏ đá có một lượng lớn là phế phẩm bột đá từ việc khai thác đá". Luận án tiến sĩ kỹ thuật chuyên ngành Cơ kỹ thuật (mã số 9 52 01 01) của tác giả Vương Lê Thắng, bảo vệ tại Trường Đại học Bách khoa – Đại học Đà Nẵng dưới sự hướng dẫn khoa học của PGS. Lê Cung và TS. Nguyễn Đình Sơn, mang tựa đề: "Nghiên cứu ứng dụng sóng siêu âm dự đoán cường độ chịu nén và vết nứt của bê tông sử dụng tro bay và bột đá". Đây là công trình tiên phong giải quyết bài toán cơ học tính toán và kiểm định không phá hủy (NDT - Non-Destructive Testing) cho dòng bê tông xanh sử dụng phụ phẩm công nghiệp.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các phương pháp kiểm tra chất lượng bê tông truyền thống dựa trên vận tốc xung siêu âm UPV (Ultrasonic Pulse Velocity) theo TCVN 9357:2012 thường chỉ xây dựng mối quan hệ hồi quy đơn biến $(f_c - UPV)$ với độ chính xác hạn chế và không thể lượng hóa đồng thời các yếu tố suy giảm biên độ sóng do môi trường tán xạ phức tạp. Hơn nữa, việc tích hợp đồng thời 20% tro bay thay thế xi măng và 20% bột đá thay thế cát tạo nên cấu trúc vi mô biến đổi sâu sắc, làm thay đổi vận tốc truyền sóng, độ suy giảm cơ học và ma trận cản kết cấu. Chưa có công trình nào trong nước và quốc tế xây dựng mô hình toàn diện dự đoán đồng thời cường độ chịu nén ($f_c$), các hệ số cản Rayleigh $(\alpha, \beta)$ và chiều sâu vết nứt mở vuông góc bề mặt ($D$) cho hệ vật liệu đặc thù này.

Nghiên cứu tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions) và kiểm chứng 3 giả thuyết (Hypotheses) tương ứng:

  • RQ1: Mối quan hệ định lượng giữa thành phần cấp phối, vận tốc xung UPV, tỉ lệ suy giảm biên độ sóng ($A_2/A_1$) và cường độ chịu nén 28 ngày ($f_c$) của bê tông sử dụng tro bay và bột đá được biểu diễn tối ưu bằng mô hình toán học nào? (Hypothesis H1: Mạng nơ-ron nhân tạo đa tầng vượt trội hơn mô hình hồi quy tuyến tính đa biến trong việc nắm bắt tương tác phi tuyến giữa các biến đầu vào).
  • RQ2: Làm thế nào để xác định chính xác ma trận cản trong bài toán mô phỏng lan truyền sóng siêu âm 2D bằng phương pháp phần tử hữu hạn (FEM) đối với bê tông chứa tro bay và bột đá? (Hypothesis H2: Mô hình cản Rayleigh với cặp hệ số $\alpha, \beta$ trích xuất từ dữ liệu thực nghiệm đáp ứng chính xác sự suy giảm biên độ sóng qua mẫu).
  • RQ3: Phương pháp xác định thời gian nhiễu xạ lan truyền sóng (Time of Flight Diffraction - TOFD) đạt độ chính xác và sai số như thế nào khi dự đoán chiều sâu vết nứt mở bề mặt trên cơ sở tích hợp giữa mô phỏng số Newmark và đo đạc thực nghiệm? (Hypothesis H3: Thời gian truyền sóng dọc nhiễu xạ tại đỉnh vết nứt cho phép định lượng chính xác chiều sâu nứt $D$ với sai lệch nhỏ hơn 10% so với kích thước thực tế).

Khung lý thuyết của luận án được định hình từ 4 trụ cột: Động lực học đàn hồi môi trường liên tục (Elastodynamics), Lý thuyết Giảm chấn Rayleigh (Rayleigh Damping Theory), Lý thuyết Học máy (Artificial Neural Networks - ANN) và Cơ học sóng nhiễu xạ (TOFD Theory). Phạm vi thực nghiệm bao quát 72 cấp phối bê tông chế tạo từ vật liệu miền Trung với mẫu lập phương $15 \times 15 \times 15\text{ cm}$, khảo sát dải cấp độ bền chịu nén từ B10 đến B45 (tương đương mác 150 đến 500), xung siêu âm tần số danh định $54\text{ kHz}$ tại tuổi 28 ngày. Nghiên cứu mang ý nghĩa đột phá khi cung cấp bộ công cụ số và thực nghiệm toàn diện, hỗ trợ kỹ sư tối ưu hóa thiết kế cấp phối đạt cường độ thiết kế (ví dụ $300\text{ daN/cm}^2$) và kiểm soát khuyết tật kết cấu bê tông thân thiện môi trường.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về cơ học sóng siêu âm trong bê tông cho thấy sự phân hóa rõ rệt thành 3 hướng tiếp cận chính: mô phỏng số lan truyền sóng, dự đoán cường độ chịu nén $f_c$, và định lượng khuyết tật vết nứt.

Trong lĩnh vực mô phỏng số, hai trường phái tính toán chủ đạo là Phương pháp Sai phân hữu hạn (FDM/EFIT/FIT) và Phương pháp Phần tử hữu hạn (FEM/SEM/BFEM). Jean Virieux (1984, 1986) đã đặt nền móng cho FDM mô phỏng sóng dọc và sóng ngang trong môi trường đàn hồi thông qua sơ đồ rời rạc hai thành phần vận tốc - ứng suất. Dựa trên FDM, Fellinger (1995) và Frank Schubert (2004) phát triển kỹ thuật tích phân hữu hạn động lực học đàn hồi (EFIT) để mô phỏng sự lan truyền sóng trong môi trường dị tính gồm vữa xi măng và cốt liệu; Kazuyuki Nakahata et al. (2012, 2015) ứng dụng thành công kỹ thuật FIT 3D kết hợp chụp cắt lớp để mô phỏng bê tông 3 pha (vữa, cốt liệu, lỗ rỗng). Tuy nhiên, FDM bộc lộ nhược điểm chí mạng: khi xuất hiện vết nứt, việc xử lý biên so le ứng suất - vận tốc tại đỉnh nứt cực kỳ phức tạp (dù Asriana 2013 và Kazushi Kimoto 2016 đã nỗ lực đề xuất sơ đồ biên nhưng còn nhiều hạn chế tính toán). Ngược lại, FEM (Moser 1999, Brian R. Mace 2005, Magdalena Rucka 2017) giải quyết bài toán biên phức tạp và vết nứt linh hoạt hơn, nhưng bài toán nan giải là xác định ma trận cản $[C]$. Trong khi Nakahata (2012) bỏ qua ma trận cản, các nghiên cứu của Ramadas (2011) trên composite và Zhen Tian (2018) trên dầm bê tông đã khẳng định tính tất yếu của mô hình cản Rayleigh $[C] = \alpha [M] + \beta [K]$ để tái hiện chính xác sự tiêu tán năng lượng sóng.

Về dự đoán cường độ chịu nén, y văn tồn tại cuộc tranh luận giữa mô hình đơn biến và đa biến. Mô hình hồi quy phi tuyến một biến $f_c = a \cdot e^{b \cdot UPV}$ của Bogas et al. (2013) trên 84 cấp phối bê tông nhẹ cho thấy hệ số xác định bội $R^2$ chỉ đạt $0.61$ khi gộp chung mẫu; con số này chỉ tăng lên $0.83$ khi tách riêng từng loại cốt liệu (Arlita, Leca) và tỉ lệ nước/xi măng. Điều này chứng minh sóng siêu âm đơn lẻ không thể phản ánh trọn vẹn cường độ nếu thiếu các biến kiểm soát thành phần. Để khắc phục, các mô hình hồi quy đa biến ra đời: Wang et al. (2017) sử dụng tỉ lệ nước/chất kết dính ($w/b$), hàm lượng phế phẩm thủy tinh ($G$) và vận tốc sóng ($V_s$) để dự đoán $f'_c$ đạt $R^2 = 0.932$, sai số tuyệt đối trung bình $MAPE = 14.7%$; Mostafa Jalal (2017) tích hợp 6 tham số (xi măng, silica fume, zeolite, vụn cao su, tuổi, độ ẩm) đạt $R^2 = 0.896$; Ruoyu Jin (2017) dùng 10 thông số cho bê tông gạch tái chế bọc sợi FRP. Bước nhảy vọt tiếp theo là mô hình mạng nơ-ron nhân tạo (ANN): Łukasz Sadowski (2018) dùng mạng MLP 4 biến đầu vào đạt hệ số tương quan $R = 0.93$ (train) và $R = 0.91$ (test); Atici (2011) kết hợp tro bay, xỉ lò cao, UPV và trị số súng bật nảy RN đạt $R = 0.951$. Tại Việt Nam, các nghiên cứu của Lương Xuân Chiểu (2014) và Nguyễn Hồng Đức (2019) kết hợp UPV và súng bật nảy (hồi quy $R^2 = 0.74$), Hoàng Nhật Đức (2014) dự đoán bê tông tính năng cao ($R^2 = 0.83$), Đào Văn Đông (2020) dùng ANN cho bê tông Geopolymer ($R^2 = 0.851$), Lý Hải Bằng (2020) và Mai Thị Hải Vân (2021) ứng dụng ANN trên dữ liệu thứ cấp ($R \approx 0.9$). Tuy vậy, tất cả các công trình trên đều chưa thiết lập mô hình cho bê tông sử dụng đồng thời tro bay và bột đá tại miền Trung.

Đối với bài toán định lượng chiều sâu vết nứt mở, 4 kỹ thuật NDT cạnh tranh gồm:

  1. Impact-Echo (Carino & Sansalone 1984, Çam 2013, Hsiao 2008): Chiều sâu tính theo $T = c/(2 f_T)$, nhưng tần số đỉnh $f_T$ biến thiên phi quy luật tùy hình học mẫu.
  2. Truyền sóng bề mặt - Surface Wave Transmission (Popovic 2000, Foo Wei Lee 2015, Seong-Hoon Kee 2011): Đo tỉ số truyền $d_{BC}(f)$, song đòi hỏi điều kiện nghiêm ngặt $h/\lambda \ge 1.5$ và không gian bề mặt lớn, khó áp dụng cho nứt sâu.
  3. Siêu âm khuếch tán - Diffusion Method (Ramamoorthy 2004, Becker 2003, Matthias Seher 2013, Eunjong Ahn 2018): Đo thời gian trễ đỉnh năng lượng (Lag time), hiệu quả với vi nứt phân tán nhưng nhạy cảm với biên kết cấu.
  4. Xác định thời gian nhiễu xạ lan truyền - TOFD (Sansalone 1998, Liou 2010, TCVN 9357:2012): Khai thác sóng dọc sơ cấp đến đỉnh nứt sinh sóng thứ cấp, cho phép xác định trực tiếp chiều sâu $D = \sqrt{(c_p t / 2)^2 - H^2}$.

Luận án định vị nghiên cứu chính xác tại giao điểm chưa được khai phá: phát triển mô hình phần tử hữu hạn 2D tích hợp ma trận cản Rayleigh thực nghiệm để mô phỏng lan truyền sóng, đồng thời ứng dụng ANN và TOFD để giải quyết trọn vẹn bài toán kép (cường độ - vết nứt) cho bê tông tro bay - bột đá.

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

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

Luận án đóng góp quan trọng vào lý thuyết động lực học môi trường đàn hồi dị tính tương đương vĩ mô. Bê tông thông thường vốn là vật liệu compozit đa pha phức tạp; khi bổ sung 20% tro bay (thay thế xi măng) và 20% bột đá (thay thế cát), hệ vi cấu trúc matrix xi măng - cốt liệu mịn biến đổi sâu sắc về độ đặc chắc, diện tích tiếp xúc pha và sự phân bố lỗ rỗng mao dẫn.

  1. Mở rộng Lý thuyết Giảm chấn Rayleigh (Rayleigh Damping Formulation): Nghiên cứu chứng minh rằng ma trận cản $[C]$ trong phương trình dao động kết cấu: $$[M]{\ddot{Q}(t)} + [C]{\dot{Q}(t)} + [K]{Q(t)} = {F(t)}$$ không thể giả định cố định hoặc bỏ qua như các mô hình đơn giản hóa (ví dụ Nakahata 2012). Cặp hệ số cản khối lượng $\alpha$ $(\text{rad/s})$ và cản độ cứng $\beta$ $(\text{s/rad})$ trong biểu thức $[C] = \alpha [M] + \beta [K]$ phụ thuộc phi tuyến chặt chẽ vào tỉ lệ nước/chất kết dính, hàm lượng tro bay và bột đá. Luận án đã lượng hóa và thiết lập thành công quy luật biến thiên của $\alpha$ và $\beta$ ứng với 72 cấp phối thực nghiệm dưới tác động của xung siêu âm tần số cao ($54\text{ kHz}$).
  2. Chuyển dịch mô hình dự đoán (Paradigm Shift in Strength Prediction): Chuyển dịch từ mô hình xác định luận thuần túy (dựa trên quan hệ tuyến tính UPV - $f_c$ của TCVN 9357:2012) sang mô hình phi tuyến đa tham số tích hợp xử lý tín hiệu. Sự kết hợp giữa vận tốc truyền pha UPV và tỉ lệ suy giảm biên độ sóng $A_2/A_1$ tạo nên bộ mô tả trạng thái đàn hồi - tiêu tán hoàn chỉnh cho vật liệu.

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

Khung phân tích của luận án là sự hợp nhất có hệ thống giữa 3 trụ cột khoa học: Cơ học tính toán môi trường liên tục (Continuum Mechanics), Kỹ thuật xử lý tín hiệu sóng siêu âm (Ultrasonic Signal Processing), và Trí tuệ tính toán (Computational Intelligence / Neural Networks).

   +-----------------------------------------------------------------------+
   |                       KHUNG PHÂN TÍCH TOÀN DIỆN                       |
   +-----------------------------------------------------------------------+
                                       |
        +------------------------------+------------------------------+
        |                                                             |
        v                                                             v
+-------------------------------+                            +-------------------------------+
|     TRỤ CỘT CƠ HỌC SỐ         |                            |     TRỤ CỘT THỰC NGHIỆM & AI  |
| - PTHH đàn hồi 2D             |                            | - 72 Cấp phối Tro bay/Bột đá  |
| - Ma trận cản Rayleigh        |                            | - Đo UPV, A2/A1, E, rho, fc   |
|   [C] = alpha[M] + beta[K]    |                            | - Tích hợp mạng MLP-ANN       |
| - Tích phân Newmark trực tiếp |                            | - Hồi quy phi tuyến đa biến   |
+-------------------------------+                            +-------------------------------+
        |                                                             |
        +------------------------------+------------------------------+
                                       |
                                       v
        +-------------------------------------------------------------+
        |                 KẾT QUẢ ĐẦU RA ĐỘT PHÁ                      |
        | 1. Bộ thông số Rayleigh (alpha, beta) chuẩn xác             |
        | 2. Mô hình ANN dự đoán fc & tối ưu cấp phối (R2 > 0.95)     |
        | 3. Kiểm chứng TOFD định lượng chiều sâu nứt D (sai số < 7%) |
        +-------------------------------------------------------------+

Điều kiện biên và giả thiết giới hạn (Boundary Conditions) được xác lập chặt chẽ:

  • Bê tông được giả thiết là môi trường đàn hồi, đồng nhất và đẳng hướng ở cấp độ vĩ mô (với các thông số khối lượng riêng $\rho$, mô-đun đàn hồi $E$, hệ số Poisson $\nu = 0.2$).
  • Khảo sát lan truyền sóng phẳng 2D trong điều kiện biến dạng phẳng (Plane Strain) hoặc ứng suất phẳng (Plane Stress).
  • Vết nứt khảo sát là vết nứt mở vuông góc bề mặt, không chứa chất lỏng hay tạp chất chèn kẹt, đỉnh nứt có bán kính cong xấp xỉ không.

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 triết lý Thực chứng luận (Positivism) kết hợp phương pháp luận suy diễn thực nghiệm - mô phỏng số khép kín (Deductive Computational & Empirical Triangulation). Thiết kế nghiên cứu đa tầng bao gồm:

  • Tầng 1 - Đặc trưng hóa vật liệu: Xác định chỉ tiêu cơ lý của từng thành phần (cát, đá dăm, xi măng PC50 Sông Gianh, tro bay nhiệt điện Vũng Áng, bột đá mỏ Phước Tường).
  • Tầng 2 - Quy hoạch thực nghiệm đa thông số: Thiết kế ma trận 72 cấp phối bê tông, thay đổi có hệ thống hàm lượng nước ($N$), xi măng ($X$), tro bay ($TB = 20% X$), bột đá ($BĐ = 20% C$), cát ($C$), đá ($Đ$), kiểm soát cấp độ bền từ B10 đến B45.
  • Tầng 3 - Đo đạc đặc tính siêu âm & cơ học: Mẫu lập phương cạnh $15\text{ cm}$ chuẩn được dưỡng hộ 28 ngày theo tiêu chuẩn TCVN 3105:1993, tiến hành đo UPV, suy giảm biên độ $A_2/A_1$, nén xác định $f_c$ theo TCVN 3118:1993, đo khối lượng riêng $\rho$ và mô-đun đàn hồi $E$.
  • Tầng 4 - Mô phỏng số học 2D FEM: Xây dựng thuật toán giải hệ phương trình động lực học trên nền tảng phần mềm Matlab R2019.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ nghiêm ngặt các tiêu chuẩn kiểm soát chất lượng và khử nhiễu:

  1. Đo vận tốc xung siêu âm (UPV): Sử dụng thiết bị đo siêu âm chuyên dụng với đầu dò tần số $54\text{ kHz}$, đường kính tiếp xúc chuẩn, sử dụng chất tiếp âm chuyên dụng (gel/mỡ bôi trơn) để triệt tiêu lớp không khí tiếp xúc. Vận tốc truyền sóng dọc được tính qua: $$UPV = \frac{L}{\Delta t}$$ với $L = 0.15\text{ m}$ và $\Delta t$ là thời gian truyền xung qua mẫu.
  2. Xác định tỉ lệ suy giảm biên độ sóng ($A_2/A_1$): Ghi nhận tín hiệu xung điện áp tại đầu phát ($A_1$) và biên độ xung điện áp nhận được sau khi truyền qua chiều dày khối bê tông $15\text{ cm}$ ($A_2$). Tỉ lệ $A_2/A_1$ phản ánh trực tiếp mức độ tiêu tán và tán xạ năng lượng của môi trường vật liệu.
  3. Thuật toán xác định hệ số cản Rayleigh $(\alpha, \beta)$: Thiết lập quy trình lặp tối ưu hóa sai số giữa dạng sóng điện áp thực nghiệm và dạng sóng chuyển vị mô phỏng. Bằng cách phát xung siêu âm $54\text{ kHz}$ qua khối bê tông $15\text{ cm}$, hệ số suy giảm thực nghiệm $k_w$ được tính toán. Qua vòng lặp thuật toán số, xác định giá trị gia số $\Delta k_R$, từ đó hội tụ về nghiệm $(\alpha, \beta)$ chính xác cho từng cấp phối.
  4. Mô phỏng lan truyền sóng bằng tích phân số Newmark: Hệ phương trình vi phân chuyển động bậc hai (1.4) được giải bằng phương pháp Newmark trực tiếp với các tham số vô điều kiện ổn định $\beta_N = 0.25$ và $\gamma_N = 0.5$. Kích thước lưới phần tử $\Delta L$ và bước thời gian tích phân $\Delta t$ thỏa mãn triệt để điều kiện ổn định Courant-Friedrichs-Lewy (CFL): $$\Delta L \le \frac{\lambda_{\min}}{10} = \frac{c_{\min}}{10 f_{\max}}, \quad \Delta t \le \Delta t_{cr} = \frac{\Delta L}{c_p}$$ đảm bảo sóng không bị méo dạng số (numerical dispersion) hay suy hao nhân tạo (numerical dissipation).

Data và phân tích

Phân tích dữ liệu thực nghiệm được triển khai song song qua hai mô hình dự báo tiên tiến:

  1. Mô hình Hồi quy tuyến tính đa biến (Multivariable Linear Regression): Xây dựng 3 dạng phương trình dự báo cường độ chịu nén $f_c$ tuổi 28 ngày:
  • Mô hình 1 (Theo cấp phối): $f_c = a_0 + a_1 X + a_2 TB + a_3 C + a_4 BĐ + a_5 Đ + a_6 N$
  • Mô hình 2 (Theo đặc tính siêu âm & vật lý): $f_c = b_0 + b_1 UPV + b_2 (A_2/A_1) + b_3 \rho + b_4 E$
  • Mô hình 3 (Tích hợp toàn diện): Kết hợp cả thành phần cấp phối và các tham số sóng siêu âm. Đánh giá mô hình qua hệ số xác định hiệu chỉnh $R^2_{\text{adj}}$, tổng bình phương sai số phần dư $SSE$, và phân tích biểu đồ phân dư (Residual plots) để đảm bảo giả định phân phối chuẩn không bị vi phạm.
  1. Mô hình Mạng nơ-ron nhân tạo (ANN): Cấu trúc mạng nhiều lớp truyền thẳng (Feedforward Multilayer Perceptron - MLP) huấn luyện bằng thuật toán Levenberg-Marquardt backpropagation:
  • Lớp đầu vào (Input layer): Số nơ-ron tương ứng với các biến thành phần (Mô hình 1: 6 đầu vào; Mô hình 2: 4 đầu vào; Mô hình 3: 10 đầu vào).
  • Lớp ẩn (Hidden layer): Khảo sát tối ưu từ 3 đến 12 nơ-ron, hàm kích hoạt phi tuyến sigmoid/tansig.
  • Lớp đầu ra (Output layer): 1 nơ-ron biểu diễn cường độ nén $f_c$ hoặc 2 nơ-ron dự đoán đồng thời hệ số cản Rayleigh $(\alpha, \beta)$.
  • Bộ dữ liệu 72 cấp phối được phân chia ngẫu nhiên chuẩn: 70% Training, 15% Validation, 15% Testing.

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

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

Quá trình phân tích thực nghiệm và mô phỏng số mang lại 4 phát hiện khoa học mang tính đột phá:

  1. Quy luật tiêu tán sóng trong bê tông tro bay - bột đá: Việc thay thế đồng thời 20% tro bay cho xi măng và 20% bột đá cho cát không làm suy giảm đột biến vận tốc truyền sóng mà làm thay đổi đặc tính tiêu tán tần số cao. Cụ thể, vận tốc xung UPV dao động trong khoảng $3800\text{ m/s}$ đến $4750\text{ m/s}$ tương ứng với dải cường độ chịu nén 28 ngày từ $12\text{ MPa}$ đến $48\text{ MPa}$ (cấp độ bền B10 - B45). Tỉ lệ suy giảm biên độ sóng $A_2/A_1$ phản ánh độ đặc chắc và mức độ điền đầy lỗ rỗng của bột đá mịn tốt hơn so với thông số UPV đơn lẻ.
  2. Lượng hóa chính xác hệ số cản Rayleigh $(\alpha, \beta)$: Lần đầu tiên, một bộ cơ sở dữ liệu hoàn chỉnh gồm cặp hệ số cản Rayleigh cho 72 cấp phối bê tông tro bay - bột đá được thiết lập. Kết quả chứng minh rằng việc giả định $\alpha = 0$ hoặc cố định $\beta$ trong các mô phỏng trước đây dẫn đến sai lệch biên độ sóng nhận được lên tới 35-50%. Khi áp dụng đúng cặp $(\alpha, \beta)$ thực nghiệm vào mô hình FEM Newmark trên Matlab, đồ thị chuyển vị - thời gian tại điểm thu sóng đạt sự trùng khớp cao với tín hiệu điện áp máy hiện sóng đo được trên mẫu thực tế.
  3. Sự vượt trội của mô hình học máy ANN trong dự báo cường độ và tối ưu hóa cấp phối:
    • Mô hình hồi quy tuyến tính đa biến đạt hệ số xác định $R^2 \approx 0.82 - 0.88$.
    • Mô hình mạng nơ-ron nhân tạo (ANN) cấu trúc tối ưu (ví dụ mạng 10-8-1 cho Mô hình 3) đạt hệ số tương quan $R > 0.95$ trên cả tập huấn luyện và tập kiểm tra, với sai số tương đối $MAPE < 5%$.
    • Từ mô hình ANN đã huấn luyện, nghiên cứu đã xây dựng biểu đồ đường bao (contour plot) và xác định chính xác cấp phối tối ưu để đạt cường độ chịu nén chuẩn $300\text{ daN/cm}^2$ ($30\text{ MPa}$), giúp tiết kiệm đáng kể lượng xi măng tiêu thụ.
  4. Định lượng chính xác chiều sâu vết nứt bằng phương pháp TOFD tích hợp: Mô phỏng số 2D và đo đạc thực nghiệm trên các mẫu bê tông có tạo vết nứt mở vuông góc bề mặt với các chiều sâu nhân tạo khác nhau ($D = 20\text{ mm}, 40\text{ mm}, 60\text{ mm}, 80\text{ mm}$) chứng minh: thời gian lan truyền xung sóng dọc nhiễu xạ $t$ tuân theo chính xác quy luật hình học: $$D = \sqrt{\left(\frac{c_p t}{2}\right)^2 - H^2}$$ Sai số giữa kết quả dự đoán từ mô phỏng số và thực nghiệm so với kích thước vết nứt thực tế luôn duy trì dưới $7.2%$, khẳng định tính khả thi tuyệt đối của việc áp dụng TOFD trong kiểm định vết nứt công trình bê tông tro bay - bột đá.

Implications đa chiều

  • Về mặt Lý thuyết: Bổ sung vào kho tàng cơ học vật liệu và động lực học công trình hệ thống tham số cản nhớt Rayleigh cho bê tông xanh, mở đường cho việc mô phỏng chính xác các bài toán tải trọng động (va chạm, động đất, nổ) trên kết cấu bê tông tái chế.
  • Về mặt Phương pháp luận: Thiết lập một quy trình chuẩn hóa kết hợp giữa đo đạc NDT đa thông số (UPV, $A_2/A_1$), thuật toán tối ưu hóa số Newmark, và công cụ trí tuệ nhân tạo ANN. Khung phương pháp này hoàn toàn có thể chuyển giao để nghiên cứu các loại bê tông bền vững khác (bê tông geopolymer, bê tông cốt sợi tái chế, bê tông xỉ lò cao).
  • Về mặt Thực tiễn & Kỹ thuật: Cung cấp cho các kỹ sư xây dựng, các trạm trộn bê tông thương phẩm tại miền Trung giải pháp thiết kế cấp phối nhanh chóng, chính xác; đồng thời cung cấp phương pháp kiểm tra không phá hủy hiện trường để nghiệm thu cường độ và chẩn đoán vết nứt mà không cần khoan lấy mẫu phá hủy kết cấu.
  • Về mặt Chính sách & Môi trường: Hiện thực hóa chủ trương của Chính phủ về kinh tế tuần hoàn và phát triển vật liệu xây không nung/vật liệu xanh, giải quyết bài toán tiêu thụ 1 triệu tấn tro xỉ nhiệt điện Vũng Áng và hàng trăm nghìn tấn bột đá phế phẩm mỗi năm, giảm thiểu tàn phá tài nguyên cát sông tự nhiên.

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu (Boundary Conditions & Limitations):

  1. Giới hạn chiều không gian mô hình: Luận án chỉ tập trung xây dựng chương trình mô phỏng số cho bài toán phẳng hai chiều (2D). Trong thực tế kết cấu, trường sóng siêu âm phân tán theo không gian ba chiều (3D) với sự xuất hiện của sóng phản xạ từ các mặt bên phức tạp.
  2. Giả thiết đồng nhất hóa vĩ mô: Mô hình số giả định bê tông là môi trường đàn hồi đồng nhất và đẳng hướng, bỏ qua tương tác tán xạ vi mô tại ranh giới hạt cốt liệu thô (đá dăm) và pha vữa xi măng - tro bay.
  3. Đặc điểm hình học khuyết tật: Mới chỉ khảo sát thực nghiệm và mô phỏng trên vết nứt mở đơn lẻ, phương vuông góc với bề mặt. Các dạng khuyết tật phức tạp như vết nứt xiên góc, hệ thống nứt chân chim phân nhánh, hoặc vết nứt ngập nước chưa được khảo sát.
  4. Phạm vi nguồn gốc vật liệu: Dữ liệu thực nghiệm tập trung vào tro bay Nhà máy nhiệt điện Vũng Áng và bột đá mỏ Phước Tường (Đà Nẵng); đặc tính hóa lý của tro xỉ từ các nhà máy nhiệt điện khác có thể có sự khác biệt về độ mịn và hàm lượng mất khi nung (LOI).

Chương trình nghiên cứu tiếp theo (Future Research Agenda) bao gồm 5 định hướng:

  1. Phát triển thuật toán mô phỏng phần tử hữu hạn 3D song song (Parallel 3D FEM/SEM) xét đến cấu trúc vi mô 3 pha ngẫu nhiên (vữa, cốt liệu thô, vùng chuyển tiếp ITZ).
  2. Nghiên cứu ảnh hưởng của độ ẩm, độ bão hòa nước trong vết nứt và nhiệt độ môi trường đến vận tốc sóng UPV và hệ số cản Rayleigh.
  3. Ứng dụng kỹ thuật học sâu (Deep Learning / Convolutional Neural Networks - CNN) xử lý trực tiếp toàn bộ ảnh dạng sóng thời gian (A-scan/B-scan) để tự động nhận dạng vị trí và biên dạng vết nứt xiên.
  4. Mở rộng bộ dữ liệu thực nghiệm cho các cấp độ bền siêu cao (HPC/UHPC) sử dụng tro bay hạt siêu mịn và nano-silica.
  5. Nghiên cứu thực nghiệm trên kết cấu bê tông cốt thép chịu tải uốn thực tế để đánh giá ảnh hưởng của mật độ cốt thép đến độ chính xác của phương pháp TOFD.

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

Luận án tạo ra tác động lan tỏa mạnh mẽ trên nhiều phương diện:

  • Ảnh hưởng học thuật: Cung cấp bộ dữ liệu thực nghiệm chuẩn mực 72 cấp phối và mã nguồn mô phỏng Matlab tin cậy, đóng vai trò là tài liệu tham khảo chất lượng cao cho các nghiên cứu sinh, học viên cao học ngành Cơ kỹ thuật, Kỹ thuật Xây dựng và Cơ học Vật liệu.
  • Chuyển đổi công nghiệp vật liệu xây dựng: Thúc đẩy các doanh nghiệp sản xuất bê tông thương phẩm tự tin ứng dụng tro bay và bột đá ở tỉ lệ thay thế cao (20%), chuẩn hóa quy trình NDT giúp rút ngắn thời gian nghiệm thu công trình và giảm chi phí khoan kiểm định phá hủy.
  • Giá trị xã hội và môi trường: Đóng góp trực tiếp vào mục tiêu giảm phát thải khí nhà kính (giảm $20%$ lượng clinker xi măng tương đương giảm hàng trăm nghìn tấn $CO_2$ mỗi năm) và giải phóng diện tích đất chứa bãi thải tro xỉ tại khu vực Bắc Trung Bộ.

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

  • Nghiên cứu sinh & Nhà khoa học chuyên ngành Cơ kỹ thuật / Vật liệu: Thụ hưởng khung lý thuyết giảm chấn Rayleigh hoàn chỉnh, phương pháp mô phỏng sóng Newmark 2D và phương pháp luận tích hợp ANN trong xử lý dữ liệu cơ học.
  • Giảng viên và Viện nghiên cứu: Sử dụng luận án như một case study mẫu mực về nghiên cứu kết hợp giữa cơ học lý thuyết, mô phỏng số trên Matlab và kiểm chứng thực nghiệm nghiêm ngặt.
  • Kỹ sư R&D trạm trộn bê tông & Doanh nghiệp thi công: Ứng dụng trực tiếp bảng cấp phối tối ưu để chế tạo bê tông đạt cường độ mác $200 - 500$ từ nguồn phế thải sẵn có với giá thành hạ.
  • Cơ quan kiểm định chất lượng công trình (Sở Xây dựng, Trung tâm Kiểm định): Tiếp nhận quy trình đo đạc TOFD và tiêu chuẩn suy giảm biên độ siêu âm để đánh giá chính xác độ an toàn kết cấu và chiều sâu nứt cầu, đập thủy điện, công trình dân dụng.

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 mở rộng Lý thuyết Giảm chấn Rayleigh trong động lực học môi trường đàn hồi dị tính. Thay vì coi hệ số cản là hằng số phi vật lý, luận án đã thiết lập quy trình thực nghiệm - số hóa trích xuất chính xác cặp hệ số cản khối lượng $\alpha$ và cản độ cứng $\beta$ cho 72 cấp phối bê tông tro bay - bột đá ở tần số siêu âm $54\text{ kHz}$, hoàn thiện ma trận $[C] = \alpha [M] + \beta [K]$ trong bài toán mô phỏng lan truyền sóng bằng FEM.

2. Đột phá về phương pháp luận nghiên cứu so với các công trình quốc tế tiêu biểu?

Trả lời: So với nghiên cứu của Bogas et al. (2013) chỉ dùng hồi quy đơn biến UPV ($R^2 = 0.61$), Wang et al. (2017) dùng hồi quy đa biến ($R^2 = 0.932$, $MAPE = 14.7%$), và Nakahata et al. (2012) bỏ qua ma trận cản, luận án đã đột phá khi: (1) Tích hợp đồng thời tham số động lực học sóng ($UPV$, tỉ lệ suy giảm biên độ $A_2/A_1$) với thành phần vật liệu; (2) Sử dụng mạng ANN Levenberg-Marquardt đạt $R > 0.95$, $MAPE < 5%$; (3) Khép kín chu trình giữa mô phỏng số Newmark có cản Rayleigh và đo đạc sóng thực tế.

3. Phát hiện nào bất ngờ và phản trực giác nhất từ số liệu thực nghiệm?

Trả lời: Phát hiện bất ngờ nhất là: việc thay thế đồng thời 20% tro bay và 20% bột đá không làm suy giảm đơn điệu vận tốc sóng dọc $c_p$ như các dự đoán định tính trước đó, mà tạo ra sự tái cấu trúc biên độ suy giảm $A_2/A_1$. Các hạt bột đá siêu mịn đóng vai trò chèn lấp lỗ rỗng (micro-filler effect), làm tăng mô-đun đàn hồi cục bộ, bù đắp cho sự suy giảm tốc độ thủy hóa ban đầu của tro bay.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh không?

Trả lời: Hoàn toàn có. Luận án cung cấp chi tiết: toàn bộ 72 bảng thành phần cấp phối, quy trình trộn và dưỡng hộ 28 ngày, thông số thiết bị đo siêu âm $54\text{ kHz}$, mã nguồn thuật toán tích phân Newmark trên phần mềm Matlab R2019, sơ đồ lặp xác định hệ số cản $k_R$, và cấu trúc chi tiết từng lớp mạng ANN (input, hidden nodes, transfer functions).

5. Chương trình nghị sự nghiên cứu 10 năm (10-Year Research Agenda) được phác thảo thế nào?

Trả lời: Lộ trình 10 năm gồm: (Giai đoạn 1: 1-3 năm) Nâng cấp mô hình FEM 3D song song và tích hợp phân bố vi cấu trúc cốt liệu thực tế; (Giai đoạn 2: 3-5 năm) Phát triển thuật toán AI tự động quét khuyết tật vết nứt 3D trên hiện trường từ mảng đầu dò siêu âm (Phased Array Ultrasound); (Giai đoạn 3: 5-10 năm) Xây dựng hệ thống bản sao số (Digital Twin) theo dõi sức khỏe kết cấu thời gian thực (Structural Health Monitoring - SHM) cho các công trình hạ tầng trọng điểm sử dụng bê tông xanh.

Kết luận

Luận án tiến sĩ của tác giả Vương Lê Thắng đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp mang tính cột mốc:

  1. Xây dựng thành công chương trình mô phỏng số 2D phần tử hữu hạn bằng phần mềm Matlab R2019 mô tả chân thực quá trình lan truyền sóng siêu âm trong bê tông, tích hợp đầy đủ ma trận cản Rayleigh thông qua thuật toán tích phân Newmark trực tiếp.
  2. Thiết lập bộ cơ sở dữ liệu thực nghiệm toàn diện gồm 72 cấp phối bê tông sử dụng 20% tro bay và 20% bột đá miền Trung, bao quát dải cấp độ bền B10 đến B45 với đầy đủ các thông số: UPV, tỉ lệ suy giảm biên độ $A_2/A_1$, khối lượng riêng $\rho$, mô-đun đàn hồi $E$, và cường độ nén $f_c$ tuổi 28 ngày.
  3. Xác định chính xác bảng giá trị hệ số cản Rayleigh $(\alpha, \beta)$ cho 72 cấp phối và phát triển thành công mô hình mạng ANN dự đoán tức thời cặp hệ số cản này cho cấp phối bất kỳ.
  4. Phát triển mô hình mạng nơ-ron nhân tạo (ANN) dự báo cường độ chịu nén 28 ngày với độ chính xác vượt bậc ($R > 0.95$), thiết lập biểu đồ đường bao cho phép định lượng cấp phối tối ưu đạt mác thiết kế ($300\text{ daN/cm}^2$).
  5. Làm chủ và kiểm chứng thành công phương pháp xác định thời gian nhiễu xạ lan truyền sóng (TOFD) trong việc định lượng chiều sâu vết nứt mở vuông góc bề mặt, với độ sai lệch giữa mô phỏng số và thực nghiệm dưới $7.2%$.
  6. Mở ra hướng đi mang tính bước ngoặt cho công nghệ bê tông bền vững tại Việt Nam, kết nối hài hòa giữa cơ học lý thuyết, mô phỏng số học, trí tuệ nhân tạo và bài toán bảo vệ môi trường thế kỷ 21.