Tổng quan về luận án

Bối cảnh phòng thủ bờ biển và hải đảo của Việt Nam đặt ra yêu cầu cấp thiết về việc xây dựng hệ thống công trình ngầm (CTN) quân sự có khả năng kháng chấn và chịu tải trọng nổ cực hạn. Tại các đảo san hô xa bờ thuộc thềm lục địa phía Nam, hầu hết các công trình hiện hữu là công sự đào nông có cấp bảo vệ hạn chế. Môi trường trầm tích san hô bão hòa nước với cấu trúc dị hướng, xốp và dòn cùng tính chất xâm thực mặn khốc liệt của nước biển tạo nên thách thức đa trường vật lý phức tạp, khiến các giải pháp kết cấu bê tông cốt thép (BTCT) truyền thống nhanh chóng suy giảm phẩm cấp cơ lý do ăn mòn điện hóa.

Khoảng trống nghiên cứu then chốt (research gap) nằm ở chỗ các nghiên cứu trong nước trước đây (như Nguyễn Duy Túy [20], Nguyễn Trí Tá [13], Vũ Công Hoằng [5]) phần lớn tiếp cận theo phương pháp tải trọng truyền thống—không xét đến động học giãn nở của sản phẩm nổ và biến dạng phi tuyến lớn của đất nền trong thời gian cực ngắn ($t \le 0{,}1,\text{s}$). Đồng thời, các quan niệm nền theo mô hình đàn hồi tuyến tính đơn giản (như trong [11], [14]) hoàn toàn bất đối xứng và không phản ánh đúng cơ chế phá hủy dẻo nén vỡ của hạt san hô khi chịu áp lực sóng nén cao tần. Hơn nữa, việc ứng dụng vật liệu composite cốt sợi thủy tinh (GFRP) thay thế cốt thép trong công trình ngầm phòng thủ biển đảo chịu tác động nổ trực xạ và cận trường chưa từng được thiết lập cơ sở tính toán lý thuyết lẫn thực nghiệm toàn diện.

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

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô tả chính xác tương tác đa pha phi tuyến giữa sóng nổ bom đạn, môi trường địa chất cát san hô bão hòa nước và vỏ kết cấu ngầm trong không gian 3 chiều?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật ứng xử động lực học và khả năng kháng nổ của kết cấu bê tông cốt thanh polyme sợi thủy tinh (GFRP) có ưu việt hoặc tương đương với bê tông cốt thép truyền thống (BTCT) trong nền san hô hay không?
  • Giả thuyết 1 (H1): Phương pháp mô hình hóa tương tác đầy đủ (fully-coupled model) kết hợp thủy động lực học hạt làm mịn (SPH) và phần tử hữu hạn (FEM/Lagrange) sẽ giảm sai số áp lực sóng nổ xuống dưới 10% so với giả định tải trọng biên truyền thống.
  • Giả thuyết 2 (H2): Vật liệu GFRP với giới hạn đàn hồi cao và độ bền kéo lớn triệt tiêu được hiện tượng ăn mòn xâm thực mặn, đồng thời bảo đảm khả năng tiêu tán năng lượng nổ tương đương kết cấu thép dưới tải trọng xung kích cận tâm nổ.

Khung lý thuyết của luận án tích hợp hệ phương trình bảo toàn cơ học môi trường liên tục (khối lượng, động lượng, năng lượng), cơ chế biến dạng hạt xốp rời rạc qua mô hình bền Drucker-Prager mở rộng (MO Granular), phương trình trạng thái (Equation of State - EOS) Shock Hugoniot và Jones-Wilkins-Lee (JWL) cho sản phẩm nổ TNT. Luận án bao gồm 136 trang thuyết minh khoa học, 73 hình vẽ/đồ thị, 22 bảng tham số thực nghiệm, kế thừa và phát triển từ 80 tài liệu chuyên khảo trong nước và quốc tế.

+---------------------------------------------------------------------------------------+
|                KHUNG TƯƠNG TÁC ĐẦY ĐỦ ĐA TRƯỜNG VẬT LÝ (AUTODYN 3D)                   |
|                                                                                       |
|   +-------------------+    Sóng nổ cận tâm     +----------------------------------+   |
|   |  Túi thuốc TNT    | =====================> | Cát san hô bão hòa nước (SPH)    |   |
|   |  (EOS: JWL)       |    (SPH Particles)     | (EOS đa pha, MO Granular)        |   |
|   +-------------------+                        +----------------------------------+   |
|                                                                 ||                    |
|                                                      Truyền sóng áp lực               |
|                                                      (Coupling SPH-FEM)               |
|                                                                 \/                    |
|   +------------------------------------+       +----------------------------------+   |
|   |  Vỏ hầm Bê tông cốt GFRP / BTCT   | <==== | Cát san hô xa tâm nổ (Lagrange)  |   |
|   |  (RHT Concrete, Elastic-Fail GFRP) |       | (Transmitting Boundaries)        |   |
|   +------------------------------------+       +----------------------------------+   |
+---------------------------------------------------------------------------------------+

Literature Review và Positioning

Các nhánh nghiên cứu chuyên sâu về công trình ngầm chịu tải trọng nổ trên nền địa chất san hô phân kỳ thành ba trường phái học thuật chính:

Trường phái thứ nhất tập trung vào phương pháp tải trọng tĩnh quy đổi và tải trọng động đặt trên biên. Các quy chuẩn của Bộ Quốc phòng Hoa Kỳ như UFC 3-340-02 và công cụ SBEDs (Single Degree of Freedom Blast Effects Design Spreadsheet) [56], [67] mô hình hóa sóng nén dưới dạng xung lực tức thời $p_e = i \cdot \omega$ hoặc tải trọng ngắn hạn $p_e = K_d \cdot K_f \cdot p_m(z)$. Tại Việt Nam, Nguyễn Duy Túy [20] khảo sát lan truyền sóng nổ trong đất bão hòa nước, Nguyễn Trí Tá [13] ứng dụng mô hình mũ đàn dẻo 2D, và Vũ Công Hoằng [5] áp dụng ANSYS giải bài toán công sự nửa nổi nửa chìm. Tuy nhiên, trường phái này tách rời nguồn nổ khỏi môi trường tương tác, dẫn đến sự đánh giá sai lệch biên độ sóng nén phản xạ khi tỷ lệ khoảng cách nổ rút gọn giảm sâu.

Trường phái thứ hai nghiên cứu cơ học nền san hô nhiệt đới. Các công trình thực địa của Nguyễn Hoa Thịnh [18], Hoàng Xuân Lượng [9], và Nguyễn Thái Chung [2] đã cung cấp cơ sở dữ liệu nền tảng về địa mạo phân tầng tại quần đảo Trường Sa, xác lập các chỉ tiêu cơ lý: dung trọng tự nhiên $\gamma = 1{,}973,\text{g/cm}^3$, hệ số rỗng $e_0 = 0{,}656$, góc ma sát trong $\varphi = 30^\circ54'$, cường độ nén đơn trục san hô cành bão hòa nước đạt $5,\text{kg/cm}^2$ ($0{,}49,\text{MPa}$) và san hô tảng cứng dao động từ $25$ đến $85,\text{kg/cm}^2$ ($2{,}45 - 8{,}34,\text{MPa}$). Điểm tranh luận cốt lõi nằm ở mô hình hóa: trong khi các nghiên cứu tĩnh xem san hô là môi trường đàn hồi liên tục [11], [14], các thử nghiệm va chấn thanh nén chia tách Hopkinson (SHPB) [57] chứng minh rằng dưới tốc độ biến dạng cao ($\dot{\varepsilon} > 10^2,\text{s}^{-1}$), cát san hô trải qua quá trình nghiền vỡ cấu trúc hạt (grain crushing), làm gia tăng áp lực nước lỗ rỗng tức thời và tái sắp xếp pha hạt.

Trường phái thứ ba khảo sát kết cấu composite FRP trong công trình phòng thủ. Các tiêu chuẩn quốc tế ACI 440.1R-06, ACI 440.5-08 của Viện Bê tông Hoa Kỳ và hệ thống tiêu chuẩn Nga GOST 31938-2012, CTO HOCTPOЙ 2.9-2103 đã quy chuẩn hóa thanh composite phi kim loại. Các nghiên cứu quốc tế của Deitz et al. [36] và Razaqpur et al. [29] cùng công trình trong nước của Nguyễn Ngọc Long [3], Nguyễn Đức Năng [1] chủ yếu tập trung vào dầm bản chịu uốn, xoắn tĩnh hoặc mỏi ăn mòn.

Tiêu chí so sánh Nghiên cứu của Darina Fišerová (2006) [37] Nghiên cứu của Z. L. Wang et al. (2005) [72-74] Luận án của Trịnh Trung Tiến (2020)
Môi trường nền Đất sét/cát 3 pha tiêu chuẩn lục địa Nền đất bùn trầm tích phù sa Cát san hô lẫn cành vụn bão hòa nước biển
Mô hình vật liệu nền Mô hình hóa đa pha chuẩn tích hợp AutoDyn Mô hình đàn dẻo kết hợp Drucker-Prager trong LS-DYNA Xây dựng thuật toán chuyên dụng MCORAL gán vào User Subroutine AutoDyn
Vật liệu cốt kết cấu Thép xây dựng truyền thống Thép cường độ cao Thanh cốt sợi thủy tinh GFRP chịu nổ cận trường
Phương pháp giải Euler-Lagrange Coupling Phần tử hữu hạn phi tuyến FEM Khớp nối đa phương pháp SPH - Lagrange 3D
Kiểm chứng thực tế Thí nghiệm trong phòng thí nghiệm (nổ buồng kín) Dữ liệu mô phỏng tham chiếu Thử nghiệm nổ hiện trường kết hợp đo biến dạng động và gia tốc vi sai

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

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

Luận án thách thức giả định truyền thống về tính liên tục tuyến tính của định luật Hooke và các giới hạn đơn giản của mô hình Mohr-Coulomb cổ điển khi áp dụng cho tải trọng sóng xung kích. Trong không gian ứng suất chính 3D, mặt chảy Mohr-Coulomb xuất hiện các điểm kỳ dị góc cạnh ($J_1 \sin\varphi + \sqrt{J_2} \sin\theta + \frac{1}{3}\sqrt{J_2}\cos\theta\sin\varphi - c\cos\varphi = 0$), gây mất ổn định nghiêm trọng trong giải thuật số khi biến dạng thể tích tăng đột biến.

              ^ Trục ứng suất chính / Mật độ
              |
              |             /  Mặt chảy MO Granular: \sigma_y = f(p) + g(\rho)
              |            /   (Tính đến tái bền nén chặt và hóa mềm)
              |           /
              |          /       /  Mô hình Drucker-Prager (Nón tròn nhẵn)
              |         /       /
              |        /       /       /  Mô hình Mohr-Coulomb (Lục giác có điểm kỳ dị)
              |       /       /       /
              |      /       /       /
              +------------------------------------> Áp suất thủy tĩnh (p)

Luận án đã mở rộng lý thuyết dẻo cho vật liệu hạt xốp bằng cách tích hợp mô hình bền MO Granular, xác lập hàm chảy ứng suất: $$\sigma_y = \sigma_p(p) + \sigma_\rho(\rho)$$ Trong đó $\sigma_p(p)$ là thành phần giới hạn chảy phụ thuộc phi tuyến vào áp suất thủy tĩnh $p$, và $\sigma_\rho(\rho)$ là thành phần tái bền phụ thuộc vào sự gia tăng mật độ môi trường $\rho$ trong quá trình sóng nén đi qua. Luận án xây dựng thành công bộ thông số phương trình trạng thái dạng Shock Hugoniot: $$U_s = c_0 + s \cdot u_p$$ với $c_0$ là tốc độ âm thanh ban đầu trong môi trường san hô bão hòa nước và $s$ là hệ số nén động thể tích, mô tả chính xác sự chuyển tiếp trạng thái từ nén dẻo đến nén chặt tuyệt đối.

Hệ thống mệnh đề lý thuyết được chứng minh:

  • Mệnh đề 1 (P1): Trong môi trường cát san hô bão hòa nước, năng lượng sóng nén bị tiêu tán mạnh bởi quá trình vỡ vụn khung xương san hô cacbonat hữu cơ, làm suy giảm biên độ đỉnh áp suất nhanh hơn 18-24% so với môi trường đất cát thạch anh lục địa cùng độ ẩm.
  • Mệnh đề 2 (P2): Liên kết tiếp xúc giữa vỏ bê tông và nền san hô là liên kết cơ học một chiều có ma sát trượt thấp; việc mô hình hóa bằng phần tử trượt tiếp xúc Lagrange-Lagrange loại bỏ hoàn toàn các ứng suất kéo ảo phát sinh tại mặt sau kết cấu ngầm.

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

Khung phân tích của luận án kết hợp đồng thời ba cấu trúc lý thuyết:

  1. Lý thuyết thủy động lực học nổ: Áp dụng phương trình trạng thái Jones-Wilkins-Lee (JWL) cho thuốc nổ TNT: $$p = A\left(1 - \frac{\omega}{R_1 V}\right)e^{-R_1 V} + B\left(1 - \frac{\omega}{R_2 V}\right)e^{-R_2 V} + \frac{\omega E}{V}$$
  2. Lý thuyết dẻo môi trường rời rạc đa pha: Mô hình hóa hạt san hô, nước lỗ rỗng và bọt khí vi mô thông qua áp suất giả nhớt $q$ và tenxơ lệch ứng suất $s_{ij}$.
  3. Lý thuyết tương tác động lực học phi tuyến kết cấu composite: Xây dựng mô hình phá hủy ứng suất giới hạn (directional failure criteria) cho thanh cốt GFRP trực giao kết hợp mô hình bê tông RHT/Drucker-Prager.

Điều kiện biên hạn định (boundary conditions): Khung phân tích được giới hạn trong phạm vi nền cát san hô lẫn cành vụn bão hòa nước ngập mặn; tải trọng bom đạn nổ phá có đương lượng từ 100g đến 16,77kg TNT; khoảng cách nổ cận trường ($R < 2{,}0,\text{m}$) và nổ mặt nóc công trình ngầm đặt nông chiều sâu chôn $H \le 3{,}0,\text{m}$.

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ủ thế giới quan khoa học thực chứng (positivism) kết hợp chủ nghĩa hiện thực phản biện (critical realism). Thiết kế nghiên cứu hỗn hợp đa tầng (multi-level computational & empirical design) được thiết lập chặt chẽ:

+-----------------------------------------------------------------------------------------+
|                                    QUY TRÌNH NGHIÊN CỨU                                 |
|                                                                                         |
|   +--------------------------+     +--------------------------+     +---------------+   |
|   | 1. Khảo sát địa chất     |     | 2. Thí nghiệm nổ sóng    |     | 3. Lập trình  |   |
|   |    Cơ lý san hô bão hòa  | ==> |    Đo độ nén động        | ==> |    Mô đun số  |   |
|   |    (Hoàng Xuân Lượng [9])|     |    (Cảm biến áp điện)    |     |    `MCORAL`   |   |
|   +--------------------------+     +--------------------------+     +---------------+   |
|                                                                             ||          |
|                                                                             \/          |
|   +--------------------------+     +--------------------------+     +---------------+   |
|   | 6. Kiểm chứng sai số     | <== | 5. Nổ thử nghiệm dã ngoại| <== | 4. Mô phỏng   |   |
|   |    Đánh giá độ tin cậy   |     |    Đo biến dạng/gia tốc  |     |    AutoDyn 3D |   |
|   |    (Dưới 9.5% sai lệch)  |     |    (NI SCXI, PCB 353B33) |     |    (SPH-FEM)  |   |
|   +--------------------------+     +--------------------------+     +---------------+   |
+-----------------------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Chiến lược thu thập và xử lý dữ liệu được chuẩn hóa qua hệ thống thiết bị đo động lực học tối tân:

  • Thiết bị thu thập tín hiệu: Hệ thống thu nhận dữ liệu tốc độ cao NI SCXI-1000DC của National Instruments với tần số lấy mẫu $100,\text{kHz}$.
  • Cảm biến đo lường: Đầu đo gia tốc áp điện công nghiệp PCB Piezotronics mã hiệu 353B33 (dải đo $\pm 5000,\text{g}$) và đầu đo biến dạng động chuyên dụng PL-60-11 (chiều dài lá đo $60,\text{mm}$, điện trở $120,\Omega$) gắn trực tiếp tại các tiết diện xung yếu (giữa sườn, nóc vòm, bản đáy hầm).
  • Quy trình tam giác hóa (Triangulation): Dữ liệu được đối chiếu 3 chiều giữa giải tích phương trình trạng thái số $\leftrightarrow$ mô phỏng phần tử rời rạc AutoDyn 3D $\leftrightarrow$ kết quả hiện trường thực nghiệm tại bãi thử Tiểu đoàn 93 - Bộ Tư lệnh Công binh.
  • Độ tin cậy của phép đo: Độ biến thiên tín hiệu giữa các lần đo nổ lặp lại đạt hệ số tin cậy $R^2 > 0{,}94$, loại trừ nhiễu điện từ xung nổ bằng cáp bọc kim ba lớp và nối đất độc lập.

Data và phân tích

Giải pháp tính toán AutoDyn 3D sử dụng thuật toán tích phân tường minh sai phân trung tâm (explicit Leapfrog time integration): $$x_i^{n+1} = x_i^n + \dot{x}_i^{n+1/2} \Delta t^{n+1/2}$$ $$\dot{x}_i^{n+1/2} = \dot{x}_i^{n-1/2} + \ddot{x}_i^n \Delta t^n$$ với gia tốc nút $\ddot{x}_i = \frac{F_i}{m} + b_i$.

Để triệt tiêu hiện tượng méo lưới (mesh distortion) khi thuốc nổ và đất sát tâm nổ chịu áp lực hàng chục Gigapascal, tác giả chia tách không gian tính toán thành các phân vùng hình học tối ưu:

  1. Vùng thuốc nổ và đất sát tâm nổ: Mô hình hóa bằng phương pháp hạt không lưới SPH (Smoothed Particle Hydrodynamics) với bán kính làm mịn biến thiên $h$.
  2. Vùng đất xa và kết cấu hầm ngầm: Rời rạc hóa theo phương pháp Lagrange 3D, gán liên kết tiếp xúc phi tuyến SPH-FEM (Sliding Interface).
  3. Môi trường không khí tự do: Mô hình hóa lưới Euler trực giao.

Bước thời gian tích phân được kiểm soát ngặt nghèo bởi tiêu chuẩn hội tụ Courant-Friedrichs-Lewy (CFL): $$\Delta t \le \frac{k \cdot l}{c_a}$$ với $l$ là kích thước phần tử nhỏ nhất ($10,\text{mm}$), hệ số ổn định chọn $k = 0{,}66$, và $c_a$ là tốc độ lan truyền sóng âm trong san hô bão hòa ($D \approx 1500 - 1850,\text{m/s}$).

Chương trình phần mềm tự biên soạn MCORAL (được viết bằng mã nguồn tích hợp) tự động chuyển đổi các chỉ tiêu địa kỹ thuật tĩnh của san hô thành ma trận hằng số phương trình trạng thái Hugoniot $p(\rho, e)$, mô đun cắt biến thiên $G(\rho)$ và tiêu chí phá hủy kéo $P_{\text{hydrostatic}} \le -P_{\text{cutoff}}$.

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

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

Thứ nhất, kiểm chứng thành công sự tương thích tuyệt đối của mô hình vật liệu cát san hô bão hòa nước thông qua thực nghiệm nổ định chuẩn. Dưới lượng nổ $100,\text{g}$ và $200,\text{g}$ TNT, đỉnh áp lực sóng nén đo được tại 3 điểm đo so với kết quả mô phỏng AutoDyn 3D có độ lệch cực đại chỉ $7{,}2% - 9{,}1%$.

Áp lực sóng nén đỉnh (kPa)
  |                                                  
  |                               [Thí nghiệm thực tế]
  |                               [Mô phỏng AutoDyn 3D]
  |         ______
  |        |      |
  |        |      |
  |        |  ||  |        ______
  |        |  ||  |       |      |
  |        |  ||  |       |  ||  |        ______
  |        |  ||  |       |  ||  |       |      |
  |        |  ||  |       |  ||  |       |  ||  |
  +--------+------+-------+------+-------+------+--------> Vị trí cảm biến
           Điểm đo 1      Điểm đo 2      Điểm đo 3
           (Sai số 7.2%)  (Sai số 8.4%)  (Sai số 9.1%)

Thứ hai, phát hiện hiện tượng giảm tải chênh lệch ứng suất cục bộ ở vỏ hầm bê tông cốt sợi GFRP. Tại thời điểm $t = 1{,}5 - 2{,}0,\text{ms}$ khi chịu lượng nổ $100,\text{g}$ TNT trên nóc vòm, ứng suất tương đương Von-Mises trong kết cấu cốt GFRP phân bố đồng đều hơn kết cấu BTCT truyền thống $14{,}6%$. Do mô đun đàn hồi của thanh GFRP ($E_f \approx 45 - 55,\text{GPa}$) mềm dẻo hơn đáng kể so với cốt thép ($E_s = 200,\text{GPa}$), kết cấu hấp thụ biến dạng uốn đàn hồi lớn hơn mà không xuất hiện tập trung ứng suất cắt gây nứt vỡ giòn chân vòm.

Thứ ba, kiểm chứng phản ứng chịu tải cực hạn dưới vụ nổ công phá lớn đương lượng $16{,}77,\text{kg}$ TNT cách sườn công trình $2{,}0,\text{m}$ tại thời điểm $9{,}0,\text{ms}$. Lực dọc trục cực đại phát sinh trong cốt thép đạt $128,\text{kN}$, trong khi cốt GFRP chỉ chịu lực dọc trục đỉnh $82,\text{kN}$ nhờ khả năng dãn dài đàn hồi tích lũy, giúp tiết diện bê tông xung quanh giảm thiểu diện tích phá hoại nứt toác bề mặt tiếp xúc đất.

Thứ tư, kết quả so sánh dạng sóng gia tốc và biến dạng động tại hiện trường bãi thử: Dưới tác dụng nổ $100,\text{g}$ TNT cách thành sườn hầm $0{,}5,\text{m}$, biên độ gia tốc cực đại theo phương X đo được tại điểm 1 là $a_{\max} = 1420,\text{m/s}^2$ (mô phỏng: $1310,\text{m/s}^2$, sai số $7{,}7%$), biến dạng động tương đối $\varepsilon_{\max} = 485,\mu\varepsilon$ (mô phỏng: $445,\mu\varepsilon$, sai số $8{,}2%$).

Implications đa chiều

  • Về mặt lý thuyết: Thiết lập mô hình giải tích toàn diện mô tả cơ chế lan truyền sóng kích động trong môi trường rạn san hô bão hòa nước ngập mặn—bổ khuyết vào khoảng trống của cơ học địa kỹ thuật quân sự thế giới đối với địa hình karst/san hô biển đảo.
  • Về phương pháp luận: Chuẩn hóa quy trình mô phỏng liên hợp SPH-Lagrange 3D cho các bài toán nổ phá cận kết cấu ngầm, cung cấp bộ công cụ thuật toán mở MCORAL dùng chung cho cộng đồng nghiên cứu động lực học công trình ngầm.
  • Về ứng dụng thực tiễn: Khẳng định tính khả thi vượt trội của việc thay thế toàn bộ cốt thép bằng cốt sợi thủy tinh GFRP trong xây dựng công sự biển đảo. Giải pháp này loại bỏ nguy cơ ăn mòn cốt kim loại do ion $\text{Cl}^-$, nâng cao tuổi thọ công trình từ $15 - 20$ năm lên trên $50$ năm trong môi trường khắc nghiệt mà vẫn bảo đảm hệ số an toàn kháng nổ.
  • Về chính sách quốc phòng: Cung cấp cơ sở khoa học để Cục Tác chiến, Bộ Tư lệnh Công binh và Bộ Quốc phòng ban hành sổ tay hướng dẫn tính toán thiết kế công sự ngầm kiên cố cấp chiến lược trên các đảo tiền tiêu.

Limitations và Future Research

  1. Giới hạn về mô hình hóa địa chất: Luận án mới chỉ khảo sát chi tiết lớp "cát san hô lẫn cành vụn bão hòa nước" đồng nhất, chưa xét đến tính bất liên tục của các hang hốc karst ngầm quy mô lớn và sự phân tầng dị hướng phức hợp (xen kẹp san hô khối tảng đặc sít).
  2. Giới hạn về đặc tính phá hủy dài hạn: Thử nghiệm hiện trường chỉ ghi nhận đáp ứng động lực học tức thời ($t < 20,\text{ms}$), chưa khảo sát biến dạng dẻo dư tích lũy và suy giảm độ dính bám giữa thanh GFRP và bê tông sau nhiều chu kỳ nổ lặp lại.
  3. Chương trình nghiên cứu 5-10 năm tiếp theo:
    • Mở rộng thuật toán MCORAL cho mô hình đất san hô 3 pha hoàn chỉnh (hạt rắn - nước mặn - bọt khí không bão hòa).
    • Nghiên cứu ứng dụng bê tông Geopolymer tro bay/xỉ lò cao kết hợp cốt sợi Bazan (BFRP) và Carbon (CFRP) trong môi trường nổ nhiệt độ cao.
    • Ứng dụng trí tuệ nhân tạo (PINNs - Physics-Informed Neural Networks) để dự báo nhanh trường áp lực sóng nổ 3D theo thời gian thực phục vụ chỉ huy công binh chiến trường.

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

Nghiên cứu tạo bước ngoặt trong kỹ thuật xây dựng công trình quốc phòng ven biển và hải đảo. Về mặt học thuật, công trình đóng góp bộ tham số nén động Hugoniot của san hô vùng biển nhiệt đới Đông Nam Á, dự kiến trở thành điểm quy chiếu cho các bài báo khoa học quốc tế chuyên ngành (ISI/Scopus nhóm Q1 về Computers and Geotechnics, International Journal of Impact Engineering).

Về kinh tế - xã hội, việc chuyển đổi sang kết cấu bê tông cốt sợi GFRP tiết kiệm hàng trăm tỷ đồng chi phí sửa chữa, duy tu bảo dưỡng định kỳ các công trình phòng thủ biển đảo, đồng thời mở ra tiềm năng ứng dụng rộng rãi cho các kết cấu dân dụng như hầm cáp điện, bể ngầm chứa nước ngọt và móng cọc cảng biển trên các quần đảo xa bờ.

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

  • Nghiên cứu sinh / Học viên Cao học ngành Xây dựng Công trình đặc biệt: Tiếp cận tường minh phương pháp mô hình hóa tương tác SPH-FEM trong AutoDyn 3D và mã nguồn thuật toán MCORAL.
  • Giảng viên, Nhà khoa học Cơ kỹ thuật & Địa kỹ thuật: Kế thừa bộ cơ sở dữ liệu thực nghiệm về chỉ tiêu cơ lý động của san hô bão hòa nước dưới tải trọng nổ.
  • Kỹ sư thiết kế, Đơn vị tư vấn công binh quân sự: Sở hữu công cụ tính toán, lựa chọn khẩu độ và chiều dày vỏ hầm bê tông cốt GFRP tối ưu theo cấp tải trọng bom đạn thiết kế.
  • Cơ quan chỉ huy chiến lược (Bộ Quốc phòng, Quân chủng Hải quân): Có luận cứ khoa học chuẩn xác để quy hoạch, nâng cấp hệ thống phòng thủ đảo đảo vững chắc bảo vệ chủ quyền thềm lục địa.

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ì?
Luận án đã mở rộng mô hình chảy dẻo nén chặt MO Granular và phương trình trạng thái Shock Hugoniot cho cát san hô lẫn cành vụn bão hòa nước biển, giải quyết triệt để tính chất phi tuyến dẻo và cơ chế nghiền nát hạt xốp dưới tải trọng nổ cực hạn mà các mô hình Mohr-Coulomb hay Drucker-Prager tiêu chuẩn không mô tả được.

2. Đổi mới phương pháp luận so với các công trình quốc tế tiền nhiệm?
So với nghiên cứu của Darina Fišerová [37] (chỉ xét đất sét/cát 3 pha thông thường) và Z.L. Wang [72] (mô hình hóa bán tương tác trong LS-DYNA), luận án đã thiết lập quy trình tích hợp toàn phần SPH - Lagrange 3D trong AutoDyn, tự lập trình module vật liệu chuyên biệt MCORAL và đối chiếu trực tiếp với thử nghiệm hiện trường đo đa kênh bằng đầu đo áp điện tốc độ cao.

3. Phát hiện bất ngờ nhất có số liệu thực nghiệm chứng minh?
Mặc dù thanh cốt sợi GFRP có mô đun đàn hồi chỉ bằng $\approx 25%$ cốt thép, kết cấu vỏ hầm bê tông cốt GFRP không bị sập đổ giòn mà còn làm giảm lực dọc trục đỉnh truyền vào cốt đến $35{,}9%$ ($82,\text{kN}$ so với $128,\text{kN}$ của thép dưới lượng nổ $16{,}77,\text{kg}$ TNT). Độ linh hoạt đàn hồi của cốt GFRP giúp tiêu tán và san bằng phổ ứng suất Von-Mises trên toàn tiết diện vòm hầm.

4. Khả năng tái lập quy trình nghiên cứu (Replication Protocol)?
Luận án cung cấp chi tiết toàn bộ các phương trình vi phân chuyển đổi, thuật toán rời rạc hóa, cấu trúc mã nguồn phụ lục MCORAL, bảng thông số ngõ vào vật liệu (khối lượng riêng, $c_0$, $s$, $G$, góc ma sát nội $\varphi$, lực dính $c$) cùng sơ đồ đấu nối cảm biến NI SCXI-1000DC, bảo đảm khả năng tái lập độc lập 100% trên các hệ thống phần mềm động lực học tiêu chuẩn.

5. Lộ trình phát triển học thuật 10 năm tới?
Hình thành trường phái nghiên cứu vật liệu composite tiên tiến chịu tải trọng xung kích cực hạn trong môi trường biển đảo Việt Nam, phát triển phần mềm mô phỏng chuyên dụng đa trường vật lý tương tác đất san hô - kết cấu ngầm - sóng nổ phục vụ chiến lược phòng thủ quốc gia.

Kết luận

  1. Xây dựng và chuẩn hóa thành công phương pháp sử dụng mô hình tương tác đầy đủ (fully-coupled model) trên nền tảng AutoDyn 3D, giải quyết đồng thời bài toán nổ, lan truyền sóng xung kích và tương tác động lực học giữa nền cát san hô bão hòa nước và công trình ngầm trong không gian 3 chiều.
  2. Thiết lập bộ tham số cơ lý động lực học và phát triển thành công chương trình MCORAL xây dựng mô hình vật liệu (EOS Hugoniot, mô hình bền MO Granular, mô hình phá hủy) cho cát san hô lẫn cành vụn bão hòa nước, đạt độ sai lệch dưới $9{,}1%$ so với thực nghiệm nổ định chuẩn.
  3. Chứng minh một cách khoa học tính vượt trội và khả năng làm việc an toàn của kết cấu bê tông cốt sợi thủy tinh (GFRP) dưới tác dụng nổ sóng cận trường, mở ra giải pháp đột phá thay thế kết cấu BTCT dễ bị ăn mòn tại các đảo tiền tiêu.
  4. Hoàn thành xuất sắc thử nghiệm dã ngoại nổ mô hình công trình ngầm tại hiện trường với hệ thống đo động NI SCXI-1000DC, cung cấp bộ dữ liệu thực chứng tin cậy về gia tốc ($1420,\text{m/s}^2$) và biến dạng động ($485,\mu\varepsilon$).
  5. Mở ra 3 hướng nghiên cứu mới: Cơ học động lực học san hô đa pha chứa khí; Kết cấu bảo vệ composite lai đa lớp chịu đạn xuyên phá nổ sâu; Hệ thống công trình ngầm thông minh tự triệt tiêu sóng chấn động biển đảo.
  6. Đóng góp sản phẩm cụ thể, có giá trị ứng dụng trực tiếp phục vụ nhiệm vụ chiến lược củng cố quốc phòng, bảo vệ vững chắc chủ quyền biển, đảo và thềm lục địa thiêng liêng của Tổ quốc.