Tổng quan về luận án

Sự phát triển nhanh chóng của kết cấu hạ tầng giao thông và công trình xây dựng dân dụng - công nghiệp đặt ra yêu cầu cấp thiết về việc kiểm soát chất lượng vật liệu cũng như đánh giá suy thoái kết cấu theo thời gian thực. Trong cơ học công trình và kỹ thuật vật liệu xây dựng, bê tông là vật liệu nhân tạo được ứng dụng phổ biến nhất, với chỉ tiêu chất lượng cốt lõi được quy định theo TCVN 4453:1995 là cường độ chịu nén ($f_c, f'_c$, hay cấp độ bền B). Tuy nhiên, việc khai thác tài nguyên thiên nhiên quá mức để sản xuất bê tông truyền thống (đặc biệt là cát sông tự nhiên và đá dăm) đang gây ra những hệ lụy môi trường nghiêm trọng. Tại khu vực miền Trung Việt Nam, sự tích tụ nguồn phế thải công nghiệp quy mô lớn—điển hình là khoảng 1 triệu tấn tro xỉ phát sinh hàng năm từ Nhà máy Nhiệt điện Vũng Áng (Hà Tĩnh) cùng khối lượng khổng lồ phế phẩm bột đá từ quá trình khai thác tại Mỏ đá Phước Tường (Đà Nẵng)—đã thúc đẩy xu hướng phát triển bê tông xanh bền vững thông qua giải pháp thay thế một phần xi măng và cát tự nhiên.

Khoảng trống học thuật (Research Gap) then chốt hiện nay nằm ở sự thiếu hụt các mô hình toán học - cơ học và công cụ kiểm tra không phá hủy (Non-Destructive Testing - NDT) có khả năng định lượng chính xác sự tương quan giữa cấu trúc vi mô, đặc tính lan truyền sóng đàn hồi và các thông số cơ lý của bê tông tái chế đa thành phần. Các phương pháp chuẩn hóa hiện hành trong nước, điển hình là TCVN 9357:2012, chủ yếu dựa trên mô hình hồi quy một biến phi tuyến tính giữa vận tốc xung siêu âm (Ultrasonic Pulse Velocity - UPV) và cường độ chịu nén. Cách tiếp cận này bộc lộ sai số lớn do bỏ qua sự suy giảm biên độ sóng, ma trận cản nhớt - đàn hồi và độ phức tạp của cấp phối. Đồng thời, việc chẩn đoán chiều sâu các khuyết tật nứt bề mặt trong môi trường bê tông dị tính chứa tro bay và bột đá vẫn chưa được mô hình hóa số một cách tường minh kết hợp kiểm chứng thực nghiệm.

Luận án của nghiên cứu sinh Vương Lê Thắng thực hiện 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 đã giải quyết triệt để bài toán khoa học này. Nghiên cứu đặt ra ba câu hỏi nghiên cứu (Research Questions) và giả thuyết khoa học tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa chính xác quá trình lan truyền sóng siêu âm trong bê tông dị tính có xét đến ma trận cản suy giảm năng lượng động học?
    • Giả thuyết 1 (H1): Mô hình giảm chấn Rayleigh kết hợp phương pháp phần tử hữu hạn (FEM) và tích phân số Newmark có thể mô phỏng trung thực trường chuyển vị và sự tiêu tán biên độ sóng qua môi trường bê tông biến tính.
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ tương quan và tính ưu việt của mô hình mạng nơ-ron nhân tạo (ANN) so với mô hình hồi quy tuyến tính đa biến (MLR) trong việc dự đoán cường độ chịu nén $f'_c$ và các hệ số cản động lực học từ dữ liệu sóng siêu âm là gì?
    • Giả thuyết 2 (H2): Mạng nơ-ron nhân tạo đa tầng có khả năng ánh xạ phi tuyến chính xác các tương tác phức tạp giữa cấp phối bê tông, vận tốc UPV và tỉ lệ suy giảm biên độ $A_2/A_1$, vượt trội hoàn toàn so với mô hình hồi quy truyền thống.
  • Câu hỏi nghiên cứu 3 (RQ3): Phương pháp xác định thời gian nhiễu xạ lan truyền (Time of Flight Diffraction - TOFD) đạt độ chính xác ra sao khi chẩn đoán kích thước chiều sâu vết nứt mở vuông góc bề mặt trong bê tông sử dụng phế phẩm công nghiệp?
    • Giả thuyết 3 (H3): Tín hiệu sóng nhiễu xạ từ đỉnh vết nứt thu thập bằng cảm biến siêu âm kết hợp mô phỏng số cho phép tính toán chiều sâu vết nứt với sai số kỹ thuật tiệm cận thực nghiệm dưới 5%.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của Cơ học động lực học đàn hồi (Elastodynamics), Lý thuyết giảm chấn Rayleigh (Rayleigh Damping Theory), Phương pháp phân tích số Phần tử hữu hạn (FEM), và Lý thuyết học máy kết nối (Connectionist Neural Networks). Nghiên cứu có phạm vi thực nghiệm toàn diện trên 72 cấp phối bê tông khác nhau đạt cấp độ bền chịu nén từ B10 đến B45 (tương đương mác M150 đến M500), sử dụng mẫu lập phương tiêu chuẩn $15 \times 15 \times 15\text{ cm}$ dưỡng hộ chuẩn 28 ngày, với hàm lượng thay thế cố định 20% tro bay cho xi măng PC50 Sông Gianh và 20% bột đá cho cát mịn, sử dụng nguồn phát sóng siêu âm tần số trung tâm 54 kHz.

Literature Review và Positioning

Nghiên cứu mô phỏng động lực học sóng siêu âm và đánh giá không phá hủy kết cấu bê tông đã trải qua nhiều giai đoạn phát triển với những cuộc tranh luận học thuật sâu sắc giữa các trường phái tiếp cận lý thuyết và thực nghiệm trên thế giới.

Về phương pháp mô phỏng số truyền sóng:

  • Trường phái Sai phân hữu hạn (FDM) khởi xướng bởi Jean Virieux (1984, 1986) đã đặt nền móng cho việc giải hệ phương trình vận tốc - ứng suất bậc một trong môi trường đàn hồi hai chiều. Mở rộng từ FDM, Frank Schubert (2004) cùng Kazuyuki Nakahata và cộng sự (2003, 2005) đã phát triển Kỹ thuật tích phân hữu hạn động lực học đàn hồi (EFIT/FIT) để mô phỏng 3D sóng siêu âm qua cấu trúc bê tông ba pha (vữa xi măng, cốt liệu thô, lỗ rỗng vi mô). Tuy nhiên, nhược điểm chí mạng của FDM/EFIT nằm ở sự phức tạp khi xử lý điều kiện biên gián đoạn hình học tại các vị trí vết nứt mở, như phân tích của Asriana (2018) và Kazushi Kimoto (2014).
  • Ngược lại, trường phái Phần tử hữu hạn (FEM) được Friedrich Moser (1999) và Brian R. Mace (2005) chứng minh có tính linh hoạt vượt trội trong việc mô hình hóa các hình học phức tạp. Magdalena Rucka (2010) đã triển khai Phương pháp phần tử phổ (SEM) để theo dõi suy thoái dầm bê tông cốt thép có vết nứt. Mặc dù vậy, đa số các công trình FEM/SEM trước đây thường giả định bỏ qua ma trận cản $[\mathbf{C}]$ hoặc xem bê tông là môi trường đàn hồi thuần túy không tiêu tán năng lượng, dẫn đến sai lệch lớn về biên độ phản hồi so với thực tế.

Về dự đoán cường độ chịu nén không phá hủy:

  • Tranh luận học thuật tồn tại giữa mô hình một biến và mô hình đa biến. Nghiên cứu thực nghiệm kinh điển của Bogas và cộng sự (2013) trên 84 cấp phối bê tông nhẹ chỉ ra rằng mô hình hồi quy đơn biến giữa UPV và $f'_c$ đạt hệ số xác định rất thấp ($R^2 = 0,61$), và chỉ cải thiện lên $R^2 = 0,83$ khi phân tách đơn lẻ từng loại cốt liệu và tỷ lệ nước/chất kết dính ($w/b$).
  • Wang và cộng sự (2017) đã ứng dụng hồi quy đa biến trên bê tông chứa phế phẩm thủy tinh với 3 thông số đầu vào ($w/b$, % thủy tinh, UPV), đạt $R^2 = 0,932$ và sai số phần trăm tuyệt đối trung bình $\text{MAPE} = 14,7%$. Mostafa Jalal (2013) đưa 6 thông số vào mô hình hồi quy đạt $R^2 = 0,896$.
  • Xu hướng tích hợp trí tuệ nhân tạo ghi nhận bước tiến nhảy vọt: Łukasz Sadowski (2015) sử dụng mạng ANN (cấu trúc 4-9-1) dự đoán bê tông biến tính bụi khoáng đạt hệ số tương quan $R = 0,93$ trên tập huấn luyện; Atici (2011) kết hợp tro bay, xỉ lò cao và trị số súng bật nảy trong mạng ANN đạt $R = 0,951$.
  • Tại Việt Nam, các nghiên cứu của Lương Xuân Chiểu (2015) và Nguyễn Hồng Đức (2018) trên bê tông geopolymer chỉ dừng lại ở đồ thị tương quan kết hợp súng bật nảy và siêu âm với độ chính xác hạn chế ($R^2 = 0,74$, $\text{SSE} = 209,8$). Hoàng Nhật Đức (2014) và Đào Văn Đông (2020) đã bước đầu thử nghiệm hồi quy đa biến và ANN ($R^2 = 0,851$) nhưng chưa có nghiên cứu nào trong nước khảo sát đồng thời bê tông sử dụng đồng thời tro bay Vũng Áng và bột đá Phước Tường kết hợp tham số suy giảm biên độ sóng siêu âm.

Về đánh giá vết nứt: Bốn nhánh phương pháp chính gồm Tác động tiếng vang (Impact-Echo - Sansalone & Carino, 1984; Çam, 2008), Lan truyền sóng bề mặt (Surface Wave Transmission - Popovic, 2000; Foo Wei Lee, 2014), Siêu âm khuếch tán (Ultrasonic Diffusion - Ramamoorthy, 2004; Becker, 2003; Seher, 2013), và Thời gian nhiễu xạ lan truyền (Time of Flight Diffraction - TOFD - Sansalone, 1998; Liou et al., 2004). Trong đó, phương pháp TOFD thể hiện tính ưu việt tuyệt đối về cơ sở toán học tường minh, không bị phụ thuộc vào tần số đỉnh bất định như Impact-Echo hay điều kiện tỷ lệ bước sóng hạn chế ($h/\lambda \ge 1,5$) của sóng bề mặt.

Luận án này định vị tại giao điểm tiên phong: kết hợp mô hình giảm chấn Rayleigh thực nghiệm vào thuật toán FEM giải sóng 2D bằng tích phân Newmark trên nền MATLAB R2019, đồng thời xây dựng hệ thống mô hình mạng nơ-ron nhân tạo đa biến kết hợp tỉ lệ suy giảm sóng ($A_2/A_1$) để dự đoán toàn diện đặc tính cơ học và kích thước khuyết tật.

Đó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 động lực học kết cấu công trình thông qua việc tích hợp tường minh tensor tiêu tán năng lượng vào phương trình vi phân chuyển động sóng tổng quát:

$$[\mathbf{M}]\ddot{\mathbf{Q}}(t) + [\mathbf{C}]\dot{\mathbf{Q}}(t) + [\mathbf{K}]\mathbf{Q}(t) = \mathbf{F}(t)$$

Trong đó:

  • $[\mathbf{M}]$, $[\mathbf{C}]$, $[\mathbf{K}]$ lần lượt là ma trận khối lượng, ma trận cản nhớt, và ma trận độ cứng tổng thể của hệ kết cấu rời rạc.
  • $\mathbf{Q}(t)$, $\dot{\mathbf{Q}}(t)$, $\ddot{\mathbf{Q}}(t)$ là các vector chuyển vị, vận tốc và gia tốc nút tại thời điểm $t$.
  • $\mathbf{F}(t)$ là vector tải trọng kích động động lực học ngoài.

Đóng góp lý thuyết đột phá thể hiện qua việc chứng minh mô hình giảm chấn tỷ số Rayleigh:

$$[\mathbf{C}] = \alpha[\mathbf{M}] + \beta[\mathbf{K}]$$

Trong đó $\alpha$ (rad/s) là hệ số cản khối lượng và $\beta$ (s/rad) là hệ số cản độ cứng, hoàn toàn có khả năng biểu diễn chính xác cơ chế tán xạ và hấp thụ năng lượng sóng đàn hồi trong môi trường bê tông xanh chứa 20% tro bay và 20% bột đá. Thay vì áp dụng các hằng số cản lý thuyết áp đặt từ vật liệu thép hay bê tông truyền thống, luận án đã thiết lập quy luật biến thiên thực nghiệm của $\alpha$ và $\beta$ theo 72 cấp phối cụ thể.

Mô hình lý thuyết dự đoán cường độ chịu nén được chuẩn hóa từ hàm truyền thống đơn biến sang không gian trạng thái đa chiều $f'_c = \Phi(\mathbf{X})$, với vector đầu vào $\mathbf{X} = [\text{Xi măng}, \text{Cát}, \text{Đá}, \text{Bột đá}, \text{Tro bay}, \text{Nước}, \text{UPV}, A_2/A_1, \rho, E]^T$, tạo bước chuyển dịch từ mô tả hiện tượng học thuần túy sang mô hình hóa cơ chế truyền tải sóng đa tham số.

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

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết:

  1. Lý thuyết Động lực học Đàn hồi Môi trường Liên tục (Elastodynamics): Thiết lập phương trình cân bằng vi phân Navier-Cauchy với quan hệ ứng suất - biến dạng theo Định luật Hooke tổng quát $\sigma_{ij} = C_{ijkl}\varepsilon_{kl}$.
  2. Cơ học Môi trường Tiêu tán Năng lượng (Dissipative Continuum Mechanics): Tích hợp suy giảm biên độ sóng do tán xạ cấu trúc hạt và suy hao nhớt ma sát nội tại thông qua hệ số suy giảm thực nghiệm $k_w$ và quy đổi số gia $\Delta k_R$.
  3. Lý thuyết Học máy Kết nối Phi tuyến (Connectionist Learning Theory): Mô hình hóa mạng Perceptron đa tầng (MLP) với thuật toán lan truyền ngược (Backpropagation) để giải mã mối quan hệ phi tuyến cao giữa vi cấu trúc vật liệu và đặc trưng lan truyền âm.

Điều kiện biên và giới hạn áp dụng (Boundary Conditions) được xác lập nghiêm ngặt: môi trường bê tông được xem là vật liệu đàn hồi đồng nhất và đẳng hướng tương đương ở quy mô vĩ mô (macroscopic scale); bài toán mô phỏng khảo sát trong hệ tọa độ 2D phẳng; sóng kích hoạt có tần số trung tâm $54\text{ kHz}$; vết nứt khảo sát là vết nứt hở, định hướng vuông góc với bề mặt tiếp nhận sóng.

+-----------------------------------------------------------------------------------+
|                            KHUNG PHÂN TÍCH TÍCH HỢP                               |
+-----------------------------------------------------------------------------------+
|  [Động lực học Đàn hồi]  +  [Mô hình Cản Rayleigh]  +  [Mạng Nơ-ron Nhân tạo ANN] |
|   (Navier-Cauchy & FEM)       ([C] = a[M] + b[K])       (Dự đoán f'c, a, b & TOFD)|
+-----------------------------------------------------------------------------------+
                                         |
                                         v
+-----------------------------------------------------------------------------------+
|                               QUY TRÌNH THỰC THI                                  |
|  1. 72 Cấp phối thực nghiệm (20% Tro bay Vũng Áng + 20% Bột đá Phước Tường)      |
|  2. Thu thập dữ liệu: UPV, Suy giảm A2/A1, Khối lượng riêng, Mô-đun đàn hồi, f'c  |
|  3. Giải số thuật toán Newmark trên MATLAB R2019 (Kiểm soát bước thời gian Dt)   |
|  4. Định lượng kích thước vết nứt qua Thời gian Nhiễu xạ Lan truyền (TOFD)        |
+-----------------------------------------------------------------------------------+

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ thế giới quan thực chứng (Positivism) kết hợp phương pháp luận số trị và thực nghiệm lượng hóa chính xác. Thiết kế đa tầng tích hợp giữa:

  • Tầng 1 - Mô phỏng số cao cấp: Xây dựng mã nguồn chuyên dụng trên nền tảng MATLAB R2019 giải phương trình dao động đàn hồi bằng phương pháp phần tử hữu hạn 2D. Thuật toán tích phân Newmark được thiết lập với các tham số trung bình không đổi $\gamma_N = 0,5$ và $\beta_N = 0,25$, đảm bảo tính ổn định vô điều kiện về mặt toán học.
  • Tầng 2 - Thực nghiệm vật lý: Khảo sát toàn diện 72 tổ hợp cấp phối bê tông biến tính tại phòng thí nghiệm chuyên sâu, tạo lập cơ sở dữ liệu mẫu lớn phục vụ huấn luyện và kiểm định mô hình học máy.

Tiêu chí ổn định bước thời gian giải tích được kiểm soát chặt chẽ theo điều kiện Courant-Friedrichs-Lewy (CFL):

$$\Delta t \le \Delta t_{cr} = \frac{\Delta L}{c_p}$$

Trong đó $\Delta L$ là kích thước phần tử nhỏ nhất của lưới phần tử hữu hạn tam giác/chữ nhật, $c_p = \sqrt{\frac{\lambda + 2\mu}{\rho}} = \sqrt{\frac{E(1-\nu)}{\rho(1+\nu)(1-2\nu)}}$ là vận tốc truyền sóng dọc lý thuyết trong bê tông, với $\lambda, \mu$ là các hệ số Lamé, $E$ là mô-đun đàn hồi, $\nu$ là hệ số Poisson, và $\rho$ là khối lượng riêng của vật liệu.

Quy trình nghiên cứu rigorous

Quy trình chế tạo mẫu và thực nghiệm được tiêu chuẩn hóa theo chuỗi các bước:

  1. Lựa chọn và kiểm chuẩn vật liệu:
    • Xi măng Portland PC50 Sông Gianh với khối lượng thể tích $3,12\text{ g/cm}^3$, giới hạn bền nén 28 ngày đạt $52,5\text{ MPa}$.
    • Cát vàng sông tự nhiên qua sàng chuẩn hạt.
    • Đá dăm Mỏ đá Phước Tường (Đà Nẵng) cỡ hạt $5-20\text{ mm}$.
    • Bột đá phế phẩm sàng qua mắt lưới nhỏ thay thế 20% khối lượng cát.
    • Tro bay Nhà máy nhiệt điện Vũng Áng (Hà Tĩnh) đạt độ mịn cao thay thế 20% lượng xi măng.
  2. Chế tạo mẫu: Đúc 72 tổ hợp cấp phối với các tỷ lệ nước/chất kết dính, cốt liệu biến thiên có kiểm soát; dưỡng hộ mẫu lập phương $150 \times 150 \times 150\text{ mm}$ trong môi trường ẩm tiêu chuẩn nhiệt độ $27 \pm 2^\circ\text{C}$ đủ 28 ngày.
  3. Đo đạc tham số sóng và cơ lý:
    • Sử dụng thiết bị đo siêu âm chuyên dụng phát xung tần số $54\text{ kHz}$ để xác định vận tốc xung siêu âm xuyên mẫu ($UPV$).
    • Thu nhận biên độ sóng phát $A_1$ và biên độ sóng nhận $A_2$ để tính tỉ lệ suy giảm biên độ $A_2/A_1$.
    • Cân mẫu xác định khối lượng riêng $\rho$ và kiểm tra mô-đun đàn hồi tĩnh $E$.
    • Thí nghiệm nén phá hủy trên máy nén thủy lực tự động để ghi nhận cường độ chịu nén thực tế $f'_c$ tại tuổi 28 ngày.
  4. Xác định hệ số cản Rayleigh thực nghiệm: Áp dụng thuật toán lặp xác định số gia suy giảm sóng $\Delta k_R$ và hệ số cản thực nghiệm $k_w$ qua sự sai lệch dạng sóng thực nghiệm và mô phỏng trên mẫu chuẩn, từ đó giải ngược hệ phương trình tìm cặp hệ số $\alpha$ và $\beta$.
  5. Thực nghiệm đo chiều sâu vết nứt TOFD: Chế tạo các khối bê tông có tạo rãnh nứt nhân tạo vuông góc bề mặt với độ sâu kiểm soát ($D = 10, 20, 30, 40, 50\text{ mm}$). Bố trí đầu phát và đầu thu đối xứng qua vết nứt tại hai cự ly $H_1 = 150\text{ mm}$ và $H_2 = 300\text{ mm}$. Chiều sâu vết nứt thực nghiệm được tính toán độc lập theo biểu thức chuẩn hóa:

$$D = \sqrt{\frac{c_p^2 t^2}{4} - H^2} = \frac{c_p}{2}\sqrt{t_2^2 - t_1^2}$$

Data và phân tích

Phân tích thống kê và mô hình hóa dữ liệu được thực thi qua hai phương pháp chủ đạo:

  • Mô hình Hồi quy Tuyến tính Đa biến (MLR): Thiết lập phương trình hồi quy tuyến tính dạng ma trận $\mathbf{Y} = \mathbf{X}\mathbf{\Theta} + \mathbf{E}$, giải bằng phương pháp bình phương tối thiểu (OLS).
  • Mạng Nơ-ron Nhân tạo (ANN): Thiết lập kiến trúc mạng truyền thẳng Feedforward đa tầng với hàm kích hoạt phi tuyến (Sigmoid/Tanh) tại lớp ẩn và hàm tuyến tính (Purelin) tại lớp đầu ra. Dữ liệu từ 72 cấp phối được phân chia ngẫu nhiên thành 3 tập dữ liệu độc lập: 70% tập Huấn luyện (Training), 15% tập Kiểm chuẩn (Validation), và 15% tập Kiểm tra (Testing). Quá trình tối ưu trọng số liên kết sử dụng thuật toán Levenberg-Marquardt (trainlm).

Độ tin cậy và sự phù hợp của mô hình được đánh giá thông qua các chỉ số thống kê chuẩn mực: Hệ số xác định bội ($R^2$), Sai số bình phương trung bình (MSE), Sai số tuyệt đối trung bình (MAE), và Sai số phần trăm tuyệt đối trung bình (MAPE):

$$\text{MAPE} = \frac{1}{n}\sum_{i=1}^n \left|\frac{y_i - \hat{y}_i}{y_i}\right| \times 100%$$

$$R^2 = 1 - \frac{\sum_{i=1}^n (y_i - \hat{y}i)^2}{\sum{i=1}^n (y_i - \bar{y})^2}$$

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

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

  1. Hiệu năng vượt trội của mô hình mạng nơ-ron nhân tạo đa biến: Kết quả phân tích đối chiếu chỉ ra rằng mô hình hồi quy tuyến tính đa biến (MLR) chỉ giải thích được khoảng 75% đến 83% sự biến thiên của cường độ chịu nén ($R^2 = 0,75 - 0,83$), trong khi mạng ANN đa tầng đạt hệ số tương quan vượt bậc ($R^2 > 0,93 - 0,95$, $R_{\text{test}} \approx 0,91 - 0,94$) trên toàn bộ dải cấp độ bền từ B10 đến B45.
  2. Vai trò định lượng của tham số suy giảm biên độ sóng ($A_2/A_1$): Lần đầu tiên, biến số tỉ lệ suy giảm biên độ sóng siêu âm $A_2/A_1$ được chứng minh là thông số có trọng số độ nhạy cao trong mô hình học máy. Việc bổ sung $A_2/A_1$ cùng với $UPV$ giúp giảm sai số dự đoán cường độ nén hơn 35% so với mô hình chỉ dựa thuần túy vào vận tốc truyền sóng.
  3. Xác lập bộ thông số cản Rayleigh ($\alpha, \beta$) chuẩn hóa cho bê tông tro bay - bột đá: Luận án đã giải mã và công bố bảng giá trị hệ số cản Rayleigh thực nghiệm cho 72 cấp phối bê tông biến tính miền Trung. Kết quả cho thấy hệ số cản khối lượng $\alpha$ dao động trong khoảng $150 - 450\text{ rad/s}$ và hệ số cản độ cứng $\beta$ biến thiên trong khoảng $1,2 \times 10^{-6} - 4,8 \times 10^{-6}\text{ s/rad}$, phản ánh rõ nét mức độ cản dị tính do cấu trúc hạt bột đá và phản ứng puzolan của tro bay tạo nên.
  4. Sự đồng quy chính xác giữa mô phỏng FEM và thực nghiệm TOFD: Mô phỏng số sự lan truyền sóng dọc nhiễu xạ qua đỉnh vết nứt cho dạng sóng chuyển vị nút tại điểm thu trùng khớp hoàn toàn về thời điểm xuất hiện xung đầu (first-arrival time) so với tín hiệu đo thực tế trên máy hiện sóng số. Sai số tuyệt đối trong việc dự đoán chiều sâu vết nứt mở từ $10\text{ mm}$ đến $50\text{ mm}$ chỉ dao động trong khoảng $1,2 - 2,8\text{ mm}$ (sai số tương đối $< 6,5%$).
  5. Thiết lập biểu đồ bao (Contour Plot) phục vụ tối ưu cấp phối ngược: Nghiên cứu đã xây dựng thành công công cụ đồ họa không gian cho phép kỹ sư tra cứu ngược tỷ lệ cấp phối tro bay và bột đá tối ưu nhằm đạt đúng cường độ thiết kế mục tiêu $300\text{ daN/cm}^2$ (cấp độ bền B22.5 / mác 300) mà vẫn tận dụng tối đa 20% phế thải công nghiệp.

Implications đa chiều

  • Đóng góp học thuật lý thuyết: Mở rộng phạm vi ứng dụng của lý thuyết đàn hồi - tiêu tán năng lượng Rayleigh sang lĩnh vực vật liệu xây dựng composite phức hợp; cung cấp cơ sở dữ liệu thực nghiệm chuẩn mực đầu tiên về tham số âm học của bê tông sử dụng phế phẩm miền Trung.
  • Đột phá phương pháp luận: Chuẩn hóa quy trình kết hợp mô phỏng số FEM Newmark với phương pháp kiểm tra không phá hủy TOFD, cung cấp công cụ chẩn đoán khuyết tật nứt với độ tin cậy vượt trội so với các công thức kinh nghiệm trong TCVN 9357:2012.
  • Ứng dụng thực tiễn kỹ thuật: Cung cấp cho các trạm trộn bê tông thương phẩm và nhà thầu xây dựng quy trình kiểm soát chất lượng nhanh, cho phép nghiệm thu cường độ bê tông trực tiếp tại hiện trường mà không cần khoan lấy mẫu lõi phá hủy kết cấu.
  • Chính sách và Môi trường: Đóng góp luận cứ khoa học vững chắc hỗ trợ các cơ quan quản lý ban hành tiêu chuẩn tái sử dụng tro xỉ nhiệt điện và bột đá mỏ, góp phần thực hiện hóa mục tiêu Net Zero Carbon và kinh tế tuần hoàn trong ngành công nghiệp sản xuất vật liệu xây dựng Việt Nam.

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 các giới hạn nghiên cứu (Boundary Conditions):

  1. Giả thiết môi trường mô phỏng 2D: Mô hình số phần tử hữu hạn mới dừng lại ở bài toán biến dạng phẳng/ứng suất phẳng 2D, chưa xét đến hiệu ứng tán xạ không gian 3D phức tạp tại các góc cạnh và sự phản xạ sóng từ cốt thép chịu lực bên trong kết cấu.
  2. Hình học vết nứt cố định: Nghiên cứu chỉ tập trung vào vết nứt hở thẳng đứng, định hướng vuông góc với bề mặt. Thực tế các công trình thường xuất hiện vết nứt xiên, vết nứt nhánh ngầm hoặc vết nứt bị lấp đầy bởi nước và tạp chất.
  3. Cố định tỷ lệ thay thế vật liệu phế phẩm: Thực nghiệm chỉ khảo sát cố định ở mức 20% tro bay thay thế xi măng và 20% bột đá thay thế cát, chưa xây dựng ma trận thay đổi đa cấp tỷ lệ (ví dụ: 10%, 30%, 40%) để tìm điểm uốn tới hạn của cấu trúc vật liệu.
  4. Tần số kích hoạt đơn lẻ: Thí nghiệm và mô phỏng thực hiện tại tần số trung tâm $54\text{ kHz}$, chưa mở rộng khảo sát hiện tượng phân tán sóng (wave dispersion) trên dải tần rộng (Broadband Ultrasonics từ $20\text{ kHz}$ đến $500\text{ kHz}$).

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

  • Phát triển mô hình mô phỏng số 3D đa pha Meso-scale mô tả vi cấu trúc hạt cốt liệu thực tế kết hợp kỹ thuật chụp cắt lớp vi tính công nghiệp (Industrial CT-Scan).
  • Mở rộng thuật toán TOFD và mạng ANN để nhận dạng tự động vết nứt xiên bất kỳ và mạng lưới nứt vi mô phân tán trong bê tông cốt thép.
  • Nghiên cứu ảnh hưởng của môi trường ăn mòn ven biển miền Trung (ion Cl-, độ ẩm, chu kỳ khô - ướt) đến các hệ số cản Rayleigh của bê tông tro bay theo thời gian phục vụ dài hạn (1 năm, 5 năm, 10 năm).

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

Công trình nghiên cứu 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 nguồn trích dẫn khoa học phong phú trong các tạp chí chuyên ngành kết cấu (như NDT & E International, Construction and Building Materials, Cement and Concrete Composites) nhờ cung cấp trọn vẹn bộ dữ liệu 72 cấp phối và tham số giảm chấn Rayleigh chi tiết.
  • Tác động công nghiệp: Tạo tiền đề chuyển giao công nghệ cho các doanh nghiệp sản xuất bê tông và các đơn vị kiểm định chất lượng công trình (QA/QC) tại khu vực miền Trung và cả nước, rút ngắn thời gian thử nghiệm cấp phối từ 28 ngày xuống chẩn đoán tức thời qua sóng âm.
  • Tác động xã hội và môi trường: Giải quyết bài toán tiêu thụ hàng triệu tấn tro xỉ nhiệt điện Vũng Áng và phế phẩm mỏ đá Đà Nẵng, giảm thiểu khai thác cát sông trái phép, bảo vệ lòng sông và giảm lượng phát thải khí nhà kính $\text{CO}_2$ từ công nghiệp clinker xi măng.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận tích hợp mô phỏng số động lực học đàn hồi với kỹ thuật học máy hiện đại; khai thác bộ dữ liệu mở về đặc tính âm học bê tông biến tính.
  • Chuyên gia R&D và Kỹ sư Vật liệu: Sở hữu công cụ biểu đồ bao tối ưu hóa cấp phối tro bay - bột đá đạt mác thiết kế từ B10 đến B45 với chi phí sản xuất thấp nhất.
  • Kỹ sư Giám sát và Kiểm định Kết cấu: Ứng dụng quy trình đo đạc TOFD độ chính xác cao để đánh giá chiều sâu vết nứt công trình mà không làm suy giảm khả năng chịu lực của kết cấu hiện hữu.
  • Cơ quan Hoạch định Chính sách (Bộ Xây dựng, Sở KH&CN): Tiếp nhận cơ sở khoa học để ban hành và cập nhật các định mức kinh tế - kỹ thuật, tiêu chuẩn thiết kế bê tông sử dụng phế phẩm công nghiệp bền vữ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?

Đóng góp lý thuyết độc đáo nhất là việc tích hợp thành công mô hình tiêu tán năng lượng Rayleigh ($[\mathbf{C}] = \alpha[\mathbf{M}] + \beta[\mathbf{K}]$) vào phương trình vi phân chuyển động lan truyền sóng siêu âm 2D cho bê tông sử dụng phế phẩm công nghiệp, đồng thời xác lập quy trình giải ngược thực nghiệm để tìm chính xác cặp tham số $\alpha, \beta$. Nghiên cứu đã mở rộng Lý thuyết Động lực học Đàn hồi Môi trường Liên tục của Navier-Cauchy và phương pháp giải phần tử hữu hạn động lực học của Bathe & Wilson sang miền vật liệu composite hạt dị tính tiêu tán năng lượng phức tạp.

2. Đột phá về phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế tiêu biểu?

So với nghiên cứu của Bogas và cộng sự (2013) vốn chỉ sử dụng hồi quy đơn biến UPV với độ chính xác rất thấp ($R^2 = 0,61$), luận án đã đột phá khi tích hợp tỉ lệ suy giảm biên độ $A_2/A_1$ cùng mạng ANN đa biến đạt $R^2 > 0,93$. So với công trình mô phỏng FEM của Nakahata và cộng sự (2005) bỏ qua ma trận cản $[\mathbf{C}]$ dẫn đến bảo toàn năng lượng giả tạo, luận án đã đưa ma trận cản Rayleigh thực nghiệm vào thuật toán Newmark trên MATLAB R2019, tái tạo chính xác cả pha thời gian và độ suy giảm biên độ của sóng thực tế.

3. Phát hiện bất ngờ nhất (Most Surprising Finding) từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là tham số tỉ lệ suy giảm biên độ sóng siêu âm ($A_2/A_1$) có độ tương quan nghịch rất nhạy với cấu trúc lỗ rỗng vi mô do bột đá tạo ra, đóng vai trò như một "chỉ thị âm học" giúp nâng cao độ chính xác dự đoán cường độ nén $f'_c$ vượt bậc khi kết hợp với $UPV$. Ngoài ra, dù thay thế đến 40% vật liệu truyền thống bằng phế phẩm (20% tro bay + 20% bột đá), bê tông vẫn đạt cấp độ bền cao lên đến B45 (mác 500) khi được tối ưu hóa tỷ lệ cấp phối thích hợp.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án trình bày chi tiết toàn bộ quy trình: chỉ tiêu cơ lý cát, đá Phước Tường, tro bay Vũng Áng, xi măng PC50 Sông Gianh; bảng thành phần định lượng cho 72 cấp phối; sơ đồ kết nối thiết bị phát - thu xung siêu âm $54\text{ kHz}$; mã nguồn thuật toán Newmark giải hệ phương trình FEM trên MATLAB R2019; cấu trúc các lớp mạng ANN cùng bộ trọng số khởi tạo; và sơ đồ bố trí cảm biến đo thời gian lan truyền sóng nhiễu xạ TOFD tại hai khoảng cách $H = 150\text{ mm}$ và $300\text{ mm}$.

5. Định hướng nghiên cứu 10 năm tới (10-Year Research Agenda) được vạch ra như thế nào?

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

  • Giai đoạn 1 (1-3 năm): Nâng cấp mô hình số lên mô phỏng 3D vi cấu trúc thực (Meso-scale) kết hợp học sâu (Deep Learning/CNN) tự động xử lý ảnh tín hiệu sóng siêu âm thời gian thực.
  • Giai đoạn 2 (4-6 năm): Phát triển thiết bị phần cứng NDT cầm tay thông minh tích hợp sẵn chip vi xử lý nạp thuật toán ANN của luận án để chẩn đoán tại hiện trường.
  • Giai đoạn 3 (7-10 năm): Mở rộng nghiên cứu sang hệ thống giám sát sức khỏe kết cấu liên tục (Structural Health Monitoring - SHM) sử dụng mạng lưới cảm biến áp điện không dây (Piezoelectric Sensors) theo dõi sự hình thành vết nứt vi mô trên các công trình cầu cảng, đập thủy điện chịu tác động ăn mòn môi trường biển.

Kết luận

Luận án tiến sĩ của tác giả Vương Lê Thắng là một công trình khoa học công phu, mẫu mực, kết hợp nhuần nhuyễn giữa cơ sở lý thuyết động lực học kết cấu, kỹ thuật phân tích số trị tiên tiến và thực nghiệm kiểm chứng quy mô lớn.

Sáu đóng góp cụ thể, nổi bật của luận án bao gồm:

  1. Xây dựng hoàn chỉnh chương trình mô phỏng số 2D phần tử hữu hạn trên MATLAB R2019 giải bài toán lan truyền sóng siêu âm trong bê tông có xét đến ma trận cản Rayleigh bằng thuật toán tích phân Newmark.
  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 biến tính (20% tro bay Vũng Áng, 20% bột đá Phước Tường) đạt cấp độ bền từ B10 đến B45 với đầy đủ các chỉ tiêu âm học ($UPV, A_2/A_1$) và cơ lý ($E, \rho, f'_c$).
  3. Xác định chính xác bộ hệ số cản Rayleigh thực nghiệm ($\alpha, \beta$) cho 72 cấp phối bê tông phế phẩm, làm sáng tỏ cơ chế suy giảm năng lượng sóng trong môi trường vật liệu xanh.
  4. Phát triển hệ thống mô hình mạng nơ-ron nhân tạo (ANN) dự đoán cường độ chịu nén và hệ số cản Rayleigh với độ chính xác vượt trội ($R^2 > 0,93$), vượt trội hoàn toàn so với mô hình hồi quy tuyến tính cổ điển.
  5. Kiểm chứng thành công phương pháp thời gian nhiễu xạ lan truyền (TOFD) trong việc định lượng chiều sâu vết nứt mở vuông góc bề mặt, đạt độ sai lệch kỹ thuật dưới 6,5% giữa mô phỏng số và thực nghiệm.
  6. Thiết lập biểu đồ bao công nghệ hỗ trợ thiết kế cấp phối tối ưu ngược, giải quyết trực tiếp bài toán tận dụng nguồn phế thải công nghiệp tại miền Trung, mang lại giá trị to lớn về mặt kinh tế, kỹ thuật và bảo vệ môi trường.