Tổng quan về luận án

Bối cảnh an ninh toàn cầu và sự gia tăng các mối đe dọa từ những vụ nổ có chủ đích nhắm vào cơ sở hạ tầng dân sinh lẫn quân sự (tiêu biểu như vụ đánh bom Tòa nhà Liên bang Alfred P. Murrah tại Oklahoma năm 1995, vụ Tháp Khobar năm 1996, và biến cố 11/9/2001) đã đặt ra yêu cầu cấp bách về việc nâng cao khả năng kháng chấn và chống nổ cho các công trình xây dựng. Trong các dạng tải trọng phát sinh từ vụ nổ, áp lực sóng xung kích lan truyền trong không khí – đặc biệt là áp lực phản xạ – được xác định là nguyên nhân cốt lõi gây phá hoại kết cấu trên diện rộng. Đối với các công trình hiện hữu chưa được tính toán kháng nổ ngay từ khâu thiết kế ban đầu, việc can thiệp kết cấu đòi hỏi các giải pháp bảo vệ tăng cường có khả năng hấp thụ năng lượng vượt trội, thi công linh hoạt mà không làm gia tăng đáng kể tĩnh tải công trình.

Nghiên cứu của nghiên cứu sinh Nguyễn Công Nghị, dưới sự hướng dẫn khoa học của Cố GS.TS Vũ Đình Lợi và TS Lê Anh Tuấn tại Học viện Kỹ thuật Quân sự (Chuyên ngành: Kỹ thuật Xây dựng Công trình Đặc biệt), tập trung giải quyết khoảng trống nghiên cứu then chốt: Thiết lập cơ chế hấp thụ năng lượng động học và chuẩn hóa mô hình số tương tác đa vật liệu cho kết cấu bảo vệ dạng tấm ốp sử dụng vật liệu có cấu trúc lỗ rỗng (như xốp Extruded Polystyrene - XPS và nhôm bọt Aluminium Foam).

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

  • Câu hỏi nghiên cứu 1 (RQ1): Cơ chế tiêu tán năng lượng của vật liệu có cấu trúc lỗ rỗng diễn ra như thế nào dưới tác động của sóng xung kích nổ cực hạn?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để xây dựng và chuẩn hóa mô hình phần tử hữu hạn phi tuyến trong LS-DYNA nhằm mô tả chính xác tương tác va chạm - nổ giữa khối thuốc nổ, môi trường không khí, lớp ốp xốp và kết cấu bê tông cốt thép (BTCT)?
  • Câu hỏi nghiên cứu 3 (RQ3): Cấu tạo phân lớp (sandwich), chiều dày tấm dàn lực và kích thước phân đoạn module nào sẽ tối ưu hóa hiệu quả giảm áp lực và gia tốc truyền vào kết cấu?
  • Giả thuyết 1 (H1): Ứng xử nén dẻo với thềm chảy phẳng dài (plateau stress) của vật liệu xốp cho phép kéo dài thời gian duy trì pha nén $t_d$, làm giảm áp lực đỉnh $P_s$ truyền vào bề mặt bê tông xuống dưới ngưỡng ứng suất kéo phá hoại.
  • Giả thuyết 2 (H2): Mô hình số tích hợp phương pháp Euler-Lagrange tùy biến (ALE) kết hợp thẻ vật liệu xốp nén vỡ phi đàn hồi (MAT_063/MAT_057) có khả năng dự báo phản ứng động học của kết cấu bảo vệ với sai số dưới 15% so với thực nghiệm nổ hiện trường.

Khung lý thuyết của luận án tích hợp Cơ học môi trường liên tục, Cơ học vật liệu có cấu trúc lỗ rỗng (Cellular Solids Mechanics) và Lý thuyết thủy động lực học sóng xung kích (Rankine-Hugoniot & JWL Equation of State). Phạm vi nghiên cứu bao gồm các thí nghiệm nén tĩnh đơn trục mẫu xốp XPS/nhôm bọt, các chuỗi thực nghiệm nổ hiện trường (khối lượng thuốc nổ TNT tại các khoảng cách $60\text{ cm}$, $80\text{ cm}$, $100\text{ cm}$), và mô phỏng số phi tuyến kết cấu tấm, cột BTCT $100\text{ mm}$ chịu nổ gần.

flowchart TD
    A["Tải trọng sóng xung kích (TNT, JWL EOS)"] --> B["Tấm dàn lực mặt trước (Steel Plate)"]
    B --> C["Lớp lõi xốp lỗ rỗng (XPS / Nhôm bọt)"]
    C -- "Biến dạng nén dẻo thềm phẳng (Plateau Phase)" --> D["Hấp thụ năng lượng thể tích Wvol"]
    D --> E["Giảm đỉnh áp lực & Kéo dài thời gian tác dụng"]
    E --> F["Kết cấu BTCT được bảo vệ (Giảm nứt, chấn sụp & gia tốc)"]

Literature Review và Positioning

Các nghiên cứu giảm thiểu tác hại nổ cho kết cấu công trình trên thế giới và trong nước hiện chia thành ba dòng tiếp cận chính:

  1. Dòng nghiên cứu gia cường nội tại vật liệu: Nâng cao độ dẻo và khả năng tiêu tán năng lượng của bê tông bằng bê tông cốt sợi thép (SFRC - Nam et al., 2006), bê tông cốt sợi polyetylen (PFRC - Maalej et al., 2005), và bê tông cốt sợi thép cường độ siêu cao (UHSFRC - Song et al., 2007; Yi et al., 2012 với năng lượng đứt gãy đạt $20.000 - 40.000\text{ J/m}^2$). Dù nâng cao đáng kể khả năng chống chấn sụp (spalling), phương pháp này chỉ áp dụng được cho công trình xây mới.
  2. Dòng nghiên cứu gia cố bề mặt bằng vật liệu composite dẻo: Dán tấm polyme sợi carbon (CFRP), sợi thủy tinh (GFRP) và sợi aramid (AFRP) (Vũ Đình Lợi và cộng sự, 2010; Muszynski & Purcell, 2003; Razaqpur et al., 2007; Ohkubo et al., 2008), hoặc phủ màng lai CFRP/polyurea (CPU) (Ha et al., 2011). Điểm hạn chế căn bản của giải pháp này là kết cấu BTCT vẫn phải tiếp nhận toàn bộ áp lực nổ trực tiếp và truyền tải xung lực xuống móng, kèm theo yêu cầu thi công phức tạp với thời gian đông kết bảo dưỡng kéo dài trên 72 giờ.
  3. Dòng nghiên cứu lớp ốp hy sinh hấp thụ năng lượng bằng vật liệu xốp: Khai thác tính năng nén ép thể tích lớn của nhôm bọt (Aluminium foam), xốp tổ ong Nomex và Expandable Polystyrene (EPS) (Ashby et al., 2000; Mukai et al., 2006; Ackland, 2008; Schenker et al., 2005, 2008; Wu et al., 2010).

Trong y văn tồn tại một cuộc tranh luận học thuật lớn giữa hai quan điểm:

  • Quan điểm 1 (Gia tăng độ cứng tuyệt đối): Nhóm ủng hộ tăng cường độ và bề dày tấm chắn thép/composite để phản xạ sóng xung kích, cho rằng vật liệu xốp dễ bị nén bẹp hoàn toàn (densification) dẫn đến hiện tượng "truyền dẫn sóng xung kích tăng cường" (shock enhancement) làm phá hoại kết cấu nặng nề hơn (Karagiozova et al., 2009).
  • Quan điểm 2 (Phân tán năng lượng dẻo): Nhóm các nhà nghiên cứu đứng đầu bởi Ashby (2000), Hanssen (2002) và Schenker (2008) chứng minh rằng nếu thiết kế ứng suất thềm phẳng của lớp xốp phù hợp với cấp tải trọng nổ, năng lượng sóng nổ sẽ chuyển hóa thành công biến dạng dẻo mà không gây truyền áp lực đỉnh sang kết cấu.

Luận án của NCS. Nguyễn Công Nghị định vị chính xác vào khoảng trống này: Giải quyết mâu thuẫn giữa hiện tượng nén chặt sớm và hiệu quả hấp thụ năng lượng thông qua cấu trúc bảo vệ đa lớp (Multilayer Sandwich) kết hợp tấm dàn lực mặt trước (Front faceplate) và cơ chế phân mảnh module. So sánh với hai công trình quốc tế kinh điển:

  • So với nghiên cứu của Schenker và cộng sự (2008) thử nghiệm nổ quy mô $1000\text{ kg}$ TNT tại Israel/Đức trên tường BTCT phủ nhôm bọt ghi nhận mức giảm gia tốc $50%$, luận án đã mở rộng xây dựng quy trình chuẩn hóa số chi tiết từ vật liệu xốp mật độ thấp (XPS) đến nhôm bọt, tối ưu hóa sự tương thích giữa độ dày tấm thép dàn lực và chiều dày lớp xốp.
  • So với công trình của H. Ackland (2008) tại Úc so sánh thực nghiệm nhôm bọt, Nomex, Skydex và EPS nhưng bị giới hạn bởi độ dày mẫu không đồng nhất ($14\text{ mm}$ đến $25\text{ mm}$), luận án đã đồng nhất hóa các điều kiện biên và giải quyết triệt để bài toán mô hình hóa vật liệu vĩ mô (macro-scale constitutive modeling) thay cho mô hình phân vùng Voronoi 3D vốn không khả thi về chi phí tính toán trên kết cấu công trình lớn.

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

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

Luận án kế thừa và mở rộng thuyết biến dạng vật liệu cấu trúc lỗ rỗng của Gibson & Ashby (1997) sang miền tải trọng động lực học phi tuyến cao do nổ gây ra. Dưới tác dụng của sóng xung kích nổ, ứng xử của lớp xốp được mô hình hóa qua ba giai đoạn đặc thù:

$$\sigma(\varepsilon) = \begin{cases} E_f \cdot \varepsilon & \text{khi } \varepsilon < \varepsilon_y \text{ (Đàn hồi tuyến tính)} \ \sigma_{pl} \approx \text{const} & \text{khi } \varepsilon_y \le \varepsilon \le \varepsilon_D \text{ (Thềm chảy dẻo sụp đổ)} \ \sigma_{pl} + E_D \left(\frac{\varepsilon - \varepsilon_D}{1 - \varepsilon}\right) & \text{khi } \varepsilon > \varepsilon_D \text{ (Nén chặt - Densification)} \end{cases}$$

Năng lượng tiêu tán trên một đơn vị thể tích $W_{vol}$ được xác định bằng diện tích dưới đường cong ứng suất - biến dạng:

$$W_{vol} = \int_0^{\varepsilon_D} \sigma(\varepsilon) , d\varepsilon \approx \sigma_{pl} \cdot \varepsilon_D$$

  Ứng suất σ
      0              ε_y          ε_D

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

  • Mệnh đề 1 (P1): Lớp xốp đóng vai trò là bộ lọc cơ học (mechanical low-pass filter), làm giảm ứng suất truyền vào kết cấu nền xuống mức xấp xỉ ứng suất thềm $\sigma_{pl}$, triệt tiêu hoàn toàn các gai áp lực phản xạ đỉnh cực ngắn của sóng nổ.
  • Mệnh đề 2 (P2): Tấm dàn lực mặt trước đóng vai trò phân phối động lượng va chạm cục bộ của sóng nổ thành chuyển động tịnh tiến đồng đều của khối lượng $M_s$, ép đều lên bề mặt xốp và tối đa hóa thể tích vật liệu xốp tham gia biến dạng dẻo.
  • Mệnh đề 3 (P3): Bước chuyển hệ hình (Paradigm shift) chuyển từ cơ chế "kháng cự thụ động bằng độ cứng và khối lượng kết cấu" sang cơ chế "tiêu tán chủ động qua biến dạng thể tích của lớp vật liệu hy sinh".

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 thủy động lực học sóng nổ khí quyển: Phương trình trạng thái Jones-Wilkins-Lee (JWL) cho sản phẩm nổ TNT và Linear Polynomial EOS cho không khí lý tưởng.
  2. Cơ học vật liệu xốp phi tuyến: Mô hình vật liệu xốp nén vỡ đẳng hướng MAT_CRUSHABLE_FOAM (MAT_063) và MAT_LOW_DENSITY_FOAM (MAT_057) với tham số dỡ tải động (HU, SHAPE).
  3. Lý thuyết hư hại động lực học bê tông cốt thép: Mô hình vật liệu bê tông Winfrith (MAT_084) và Karagozian & Case Concrete (KCC/MAT_072R3) kết hợp hiệu ứng tốc độ biến dạng (Strain rate effects).
graph LR
    subgraph Multi_Physics ["Khung phân tích đa vật lý tương tác"]
        A["1. Nổ khí quyển<br>(JWL EOS + ALE Formulation)"] --> B["2. Tương tác tấm ốp<br>(Contact Penalty + Plate Dynamics)"]
        B --> C["3. Lõi xốp biến dạng dẻo<br>(MAT_063 / MAT_057 Energy Absorption)"]
        C --> D["4. Kết cấu BTCT chịu tải<br>(MAT_084 Winfrith Damage & Cracking)"]
    end

Điều kiện biên xác lập rõ: Tải trọng nổ tác dụng trong không gian tự do; khoảng cách tỷ lệ $Z = R / \sqrt[3]{W}$ từ vùng nổ gần đến trung bình; cấu trúc xốp không phục hồi biến dạng sau khi phá hủy.

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 thực chứng (Positivism) kết hợp thực chứng hậu hiện đại (Post-positivism), sử dụng phương pháp luận kết hợp chặt chẽ giữa Thực nghiệm vật lý hiện trường và Mô phỏng số động lực học phi tuyến cao. Thiết kế nghiên cứu đa cấp độ (Multi-level experimental-computational framework) được triển khai qua ba tầng:

  • Cấp độ 1: Thí nghiệm nén tĩnh xác định quan hệ $\sigma - \varepsilon$ của mẫu vật liệu xốp XPS và nhôm bọt nguyên mẫu.
  • Cấp độ 2: Thực nghiệm nổ hiện trường trên tấm BTCT trần và tấm BTCT có kết cấu bảo vệ tại Thao trường Học viện Kỹ thuật Quân sự.
  • Cấp độ 3: Mô hình hóa số phần tử hữu hạn 3D toàn phần trong phần mềm chuyên dụng LS-DYNA để phân tích sâu các biến thể kết cấu (chiều dày, số lớp, phân đoạn module).
Thí nghiệm nén tĩnh mẫu xốp (XPS / Nhôm bọt)
     Hiệu chuẩn thông số vật liệu (MAT_063 / MAT_057)
Thực nghiệm nổ hiện trường (TN_01, TN_02: Đo Áp lực, Gia tốc, Chuyển vị)
      Xác minh chéo (Cross-validation) với mô hình số LS-DYNA ALE
Mô phỏng khảo sát mở rộng cấu tạo (Phân lớp Sandwich, Phân đoạn Module, Tấm dàn lực)

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu bảo đảm tính nghiêm ngặt học thuật cao:

  1. Thu thập dữ liệu nén tĩnh: Mẫu xốp XPS kích thước tiêu chuẩn được thí nghiệm nén đơn trục trên máy nén thủy lực tự động điều khiển biến dạng, ghi nhận chính xác mô đun đàn hồi $E_f$, ứng suất thềm chảy dẻo $\sigma_{pl}$ và điểm bắt đầu nén chặt $\varepsilon_D$.
  2. Thực nghiệm nổ hiện trường:
    • Sử dụng thuốc nổ TNT đúc chuẩn dạng trụ/khối, kích nổ bằng kíp điện số 8.
    • Bố trí các cảm biến áp lực áp điện (Piezoelectric pressure transducers) tại các khoảng cách $60\text{ cm}$, $80\text{ cm}$ và $100\text{ cm}$ để đo áp lực sóng xung kích tự do và áp lực phản xạ.
    • Bố trí cảm biến đo gia tốc trở kháng (Piezoresistive accelerometers) và đầu đo chuyển vị quang học/LVDT ở mặt sau của tấm BTCT kích thước $1000 \times 1000 \times 100\text{ mm}$.
  3. Kỹ thuật tam giác giác hóa (Triangulation): Kết hợp đối chiếu đồng thời giữa lý thuyết giải tích sóng nổ (quan hệ Rankine-Hugoniot), số liệu đo thực nghiệm hiện trường và kết quả xuất từ bộ giải LS-DYNA Explicit Solver.
  4. Kiểm định độ tin cậy và giá trị: Kiểm tra độ nhạy lưới (Mesh convergence study) cho miền không khí và khối TNT với các tỷ lệ lưới từ $1/10$ đến $1/25$ (TNT1, TNT2, TNT3, TNT4). Đánh giá hệ số tỷ lệ bước thời gian (Time step scale factor - TSSFAC) biến thiên từ $0.1$ đến $0.6$ (TSF1 - TSF6) để đảm bảo điều kiện ổn định Courant-Friedrichs-Lewy (CFL).

Data và phân tích

Mô hình số được thiết lập trong LS-DYNA với các thiết lập thẻ phần tử và mô hình vật liệu cao cấp:

  • Môi trường nổ & không khí: Miền không khí và lượng nổ TNT được rời rạc hóa bằng phần tử khối Euler 8 nút, gán thuật toán *CONTROL_ALE. Không khí sử dụng thẻ *MAT_NULL kết hợp phương trình trạng thái *EOS_LINEAR_POLYNOMIAL:

$$P = C_0 + C_1\mu + C_2\mu^2 + C_3\mu^3 + (C_4 + C_5\mu + C_6\mu^2)E_0$$

  • Thuốc nổ TNT: Thẻ *MAT_HIGH_EXPLOSIVE_BURN kết hợp phương trình trạng thái *EOS_JWL:

$$P = A_{jwl}\left(1 - \frac{\omega_{jwl}}{R_1 V_N}\right)e^{-R_1 V_N} + B_{jwl}\left(1 - \frac{\omega_{jwl}}{R_2 V_N}\right)e^{-R_2 V_N} + \frac{\omega_{jwl} E_0}{V_N}$$

  • Vật liệu xốp: Sử dụng thẻ *MAT_CRUSHABLE_FOAM (MAT_063) nhập trực tiếp đường cong $\sigma - \varepsilon$ thu được từ thực nghiệm nén tĩnh, hệ số Poisson $\nu \approx 0$.
  • Bê tông và cốt thép: Bê tông sử dụng phần tử Solid 8 nút gán thẻ *MAT_WINFRITH_CONCRETE (MAT_084) có khả năng mô tả sự hình thành vết nứt và phá hủy chấn sụp; cốt thép sử dụng phần tử Beam Hughes-Liu gán thẻ *MAT_PLASTIC_KINEMATIC (MAT_003) liên kết đồng nút (Constrained Lagrange in Solid) với mạng bê tông.
  • Tương tác biên: Sử dụng thẻ *CONSTRAINED_LAGRANGE_IN_SOLID cho tương tác chất lưu - kết cấu (FSI) giữa môi trường khí nổ ALE và tấm chắn thép mặt trước, kết hợp thẻ *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE giữa các lớp sandwich.

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

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

  1. Cơ chế kéo dài xung lực và triệt tiêu đỉnh áp lực: Trái ngược với quan niệm tấm thép đặc phản xạ xung lượng lớn, sự xuất hiện của lớp xốp lõi giúp kéo dài thời gian tác dụng của pha nén $t_d$, biến đổi xung áp lực đỉnh cực hạn ngắn hạn thành dạng tải phân bố đều có biên độ thấp xấp xỉ ứng suất chảy dẻo $\sigma_{pl}$.
  2. Hiệu ứng bảo vệ phi tuyến chống hiện tượng chấn sụp (Spalling elimination): Thực nghiệm hiện trường (TN_01 và TN_02) chứng minh: tấm BTCT trần chịu nổ gần bị nứt vỡ nghiêm trọng hình mai rùa, bong tróc mảng lớn ở mặt đáy do sóng ứng suất kéo phản xạ. Ngược lại, mẫu BTCT có tấm ốp bảo vệ bằng xốp hoàn toàn không xuất hiện vết nứt phá hoại mặt sau, kết cấu giữ nguyên khả năng chịu lực đàn hồi.
  3. Hiệu quả suy giảm gia tốc động học vượt bậc: Dữ liệu gia tốc ghi nhận tại tâm mặt sau kết cấu giảm tới hơn $40% - 50%$ so với mẫu đối chứng không có lớp xốp bảo vệ, phù hợp và thẩm định chính xác các quan sát của Schenker và cộng sự (2007) trên nhôm bọt Alulight.
  4. Quy luật tương tác tấm dàn lực - lõi xốp: Chiều dày tấm thép mặt trước $\delta_{td}$ đóng vai trò quyết định. Nếu tấm dàn lực quá mỏng, sóng nổ gây đâm thủng cục bộ; nếu tấm quá dày, khối lượng quán tính lớn làm tăng động năng đập vào lớp xốp. Tồn tại một tỷ lệ tối ưu giữa độ cứng tấm dàn lực và chiều dày hấp thụ của xốp $h_{blast}$.
  5. Đột phá cấu tạo phân đoạn module (Structure segmentation): Việc phân nhỏ kết cấu bảo vệ thành các module kích thước $20 \times 20\text{ cm}$, $25 \times 25\text{ cm}$, $33 \times 33\text{ cm}$ thay vì một tấm phẳng liền khối kích thước lớn giúp triệt tiêu sóng uốn truyền dọc bề mặt, giảm hiện tượng tập trung ứng suất tại biên và nâng cao hiệu quả tiêu tán năng lượng của toàn hệ thống thêm $15% - 22%$.
Cấu hình khảo sát Áp lực đỉnh truyền vào BTCT Độ võng dư (Chuyển vị đỉnh) Tỷ lệ giảm gia tốc đỉnh Tình trạng nứt mặt sau BTCT
BTCT trần (Đối chứng) $100%$ (Cực đại) $100%$ (Chuẩn) $0%$ Phá hoại chấn sụp diện rộng, bê tông vỡ mảnh
BTCT + Tấm thép dàn lực đơn thuần $85 - 90%$ $82%$ $10 - 15%$ Xuất hiện vết nứt uốn và cắt chéo rõ rệt
BTCT + Cấu trúc cơ bản (Thép - Xốp) $35 - 45%$ $42%$ $40 - 48%$ Không chấn sụp, chỉ có vi nứt chân chim
BTCT + Sandwich nhiều lớp (Thép-Xốp-Thép-Xốp) $25 - 30%$ $31%$ $52 - 58%$ Bề mặt bê tông nguyên vẹn hoàn toàn
BTCT + Module phân đoạn $25 \times 25\text{ cm}$ (Nhôm bọt) $18 - 22%$ $24%$ $60 - 65%$ Hoàn toàn đàn hồi, không tổn thương

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình định lượng xác lập mối quan hệ giữa cấu trúc vi mô vật liệu xốp, mật độ tương đối $\rho_f / \rho_s$, ứng suất thềm phẳng $\sigma_{pl}$ và khả năng suy giảm sóng xung kích trong động lực học công trình.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình mô phỏng FSI đa vật chất trong LS-DYNA kết hợp mô hình phần tử xốp vĩ mô tương đương, cho phép các kỹ sư tính toán kết cấu chống nổ quy mô lớn mà không cần giải mã lưới vi mô Voronoi 3D tốn kém tài nguyên tính toán.
  • Về mặt ứng dụng thực tiễn: Tạo ra giải pháp tấm ốp bảo vệ module hóa: khối lượng nhẹ, khả năng cơ động cao, lắp dựng tháo dỡ nhanh chóng, sẵn sàng nâng cấp khả năng phòng thủ cho các trụ sở cơ quan trọng yếu, đại sứ quán, trạm chỉ huy quân sự và kết cấu cột trụ cầu hầm đô thị mà không làm gián đoạn công năng sử dụng của công trình.
  • Về mặt chính sách và tiêu chuẩn: Cung cấp cơ sở khoa học để biên soạn các hướng dẫn kỹ thuật thiết kế gia cố công trình chống khủng bố và nổ đô thị tại Việt Nam, tiệm cận các tiêu chuẩn quốc tế như DoD UFC 3-340-02 (Mỹ).

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  1. Giới hạn hiệu ứng nhiệt và cháy: Nghiên cứu chỉ tập trung vào thành phần áp lực cơ học của sóng xung kích nổ, chưa tích hợp trường nhiệt độ cực cao từ quả cầu lửa (fireball) và tác động cháy của sản phẩm nổ tác động lên vật liệu xốp polyme.
  2. Tính chất vật liệu một lần (Sacrificial nature): Các loại xốp XPS và nhôm bọt nghiên cứu là vật liệu hấp thụ năng lượng không phục hồi (irreversible crushing); sau khi chịu một vụ nổ làm xốp bị nén bẹp hoàn toàn đến vùng nén chặt $\varepsilon_D$, lớp bảo vệ mất đi khả năng kháng nổ cho các vụ nổ thứ cấp liên tiếp.
  3. Phạm vi trường nổ: Thí nghiệm hiện trường bị giới hạn ở khối lượng nổ TNT cục bộ tương ứng với năng lực an toàn thao trường, chưa kiểm chứng với các vụ nổ cực lớn quy mô hàng tấn TNT ngoài đô thị.

Các định hướng nghiên cứu tương lai (10-Year Research Agenda):

  • Phát triển các dòng vật liệu xốp cơ tính biến thiên (Functionally Graded Foams) và cấu trúc siêu vật liệu Auxetic có khả năng tự phục hồi hoặc tiêu tán đa xung nổ.
  • Tích hợp bài toán liên hợp Nhiệt - Nổ - Va đập mảnh văng (Thermo-mechanical Blast-Fragment coupled analysis).
  • Mở rộng thuật toán tối ưu hóa hình học đa mục tiêu cho lớp ốp sandwich bảo vệ móng ngầm và kết cấu vỏ hầm chịu nổ trong lòng đất.

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

  • Tác động học thuật (Academic Impact): Đặt nền móng học thuật cho phân ngành kết cấu bảo vệ bằng vật liệu xốp tại Việt Nam; cung cấp bộ dữ liệu thực nghiệm nổ hiện trường quý giá cho cộng đồng nghiên cứu động lực học công trình công sự và công trình phòng thủ. Dự báo mang lại tỷ lệ trích dẫn cao trong các tạp chí chuyên ngành Kết cấu, Vật liệu & Kỹ thuật Quốc phòng.
  • Chuyển đổi công nghiệp & An ninh quốc phòng (Industry & Defense Transformation): Chuyển giao công nghệ chế tạo module ốp chống nổ dạng tiền chế cho các nhà máy cơ khí quốc phòng và doanh nghiệp vật liệu xây dựng đặc chủng. Ứng dụng trực tiếp bảo vệ các tháp kiểm soát không lưu, ngân hàng dữ liệu số, giàn khoan dầu khí ngoài khơi và công trình biên giới hải đảo.
  • Lợi ích xã hội (Societal Benefits): Giảm thiểu tối đa thương vong do mảnh bê tông văng sập trong các tình huống tai nạn nổ công nghiệp (kho hóa chất, trạm khí gas) và các vụ tấn công khủng bố, 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 & Nhà khoa học chuyên sâu: Khai thác quy trình chuẩn hóa mô hình LS-DYNA ALE, bộ tham số vật liệu MAT_063/MAT_084 và hướng nghiên cứu mở rộng về siêu vật liệu hấp thụ xung lực.
  • Kỹ sư thiết kế kết cấu & Tư vấn xây dựng: Tiếp cận cẩm nang tính toán chi tiết về kích thước tấm dàn lực, tỷ lệ chiều dày lõi xốp và phương pháp phân đoạn module cho các dự án công trình đặc biệt.
  • Cơ quan quản lý Nhà nước & Bộ Quốc phòng: Sử dụng kết quả nghiên cứu làm cơ sở khoa học để ban hành định mức, quy chuẩn kỹ thuật xây dựng công sự chiến đấu, hầm trú ẩn và quy chuẩn phòng chống khủng bố đô thị.

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 độc đáo nhất là việc mở rộng Thuyết cơ học vật liệu cấu trúc lỗ rỗng của Gibson & Ashby (1997) sang bài toán động lực học va chạm nổ tốc độ cao, chứng minh rằng sự kết hợp giữa tấm dàn lực mặt trước và lõi xốp biến đổi trạng thái ứng suất phá hoại của sóng xung kích nổ từ dạng tải trọng động xung kích đỉnh nhọn thành dạng tải trọng nén tĩnh tương đương có chặn trên bởi ứng suất thềm phẳng $\sigma_{pl}$, từ đó vô hiệu hóa cơ chế tạo sóng kéo gây chấn sụp mặt sau của bê tông.

2. Tính đổi mới trong phương pháp nghiên cứu khi so sánh với các công trình quốc tế?

So với nghiên cứu của Schenker et al. (2008) chỉ tập trung vào nhôm bọt tấm phẳng liền khối và Karagiozova et al. (2009) dừng lại ở phân tích vi mô mẫu nhỏ, luận án đã đổi mới bằng cách: (1) Thiết lập thành công phương pháp mô hình hóa vĩ mô tương đương trong LS-DYNA thông qua thực nghiệm nén chuẩn hóa; (2) Phát triển giải pháp cấu tạo phân đoạn module hóa (segmentation) giúp giảm hiện tượng phản xạ sóng uốn và tập trung ứng suất mép biên; (3) Xác minh chéo toàn diện giữa thực nghiệm đo áp lực/gia tốc nổ hiện trường và mô phỏng FSI phi tuyến.

3. Phát hiện nào gây bất ngờ nhất về mặt cơ học công trình?

Phát hiện bất ngờ nhất là: Gia tăng chiều dày tấm thép chắn mặt trước không đồng nghĩa với việc gia tăng khả năng bảo vệ. Khi tấm thép mặt trước quá dày, xung lượng nổ biến toàn bộ tấm thép thành một "đầu đạn phẳng" có quán tính cực lớn, ép lớp xốp nhanh chóng vượt qua ngưỡng thềm dẻo $\varepsilon_D$ rơi vào vùng nén chặt (densification), gây truyền dẫn sóng xung kích tăng cường (shock enhancement) làm kết cấu BTCT bên dưới chịu lực ép và chuyển vị còn lớn hơn trường hợp tấm thép mỏng tối ưu.

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

Có. Toàn bộ thông số vật liệu thí nghiệm (mô đun đàn hồi, cường độ chảy, đường cong $\sigma - \varepsilon$), các hệ số phương trình trạng thái JWL ($A_{jwl}, B_{jwl}, R_1, R_2, \omega_{jwl}$), bảng thông số thẻ vật liệu LS-DYNA (MAT_063, MAT_057, MAT_084, MAT_003), kích thước hình học, điều kiện biên và quy trình chia lưới ALE/Lagrange đều được lập bảng chi tiết trong Chương 2, Chương 3 và Chương 4 của luận án, cho phép các nhà nghiên cứu độc lập tái lập mô phỏng chính xác.

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

Chương trình nghiên cứu 10 năm tập trung vào ba giai đoạn:

  • Giai đoạn 1 (1-3 năm): Tối ưu hóa cấu trúc sandwich đa tầng FGM (Functionally Graded Materials) kết hợp nhôm bọt mật độ biến thiên.
  • Giai đoạn 2 (4-6 năm): Nghiên cứu vật liệu xốp chịu tải trọng liên hợp Cháy nổ - Va đập mảnh văng đạn đạo (Blast-Thermal-Impact combined loads).
  • Giai đoạn 3 (7-10 năm): Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) dự báo nhanh phản ứng phá hủy kết cấu và tự động hóa thiết kế hệ thống vỏ bọc thông minh cho toàn bộ công trình hạ tầng đô thị trọng yếu.

Kết luận

Công trình nghiên cứu của NCS. Nguyễn Công Nghị đã giải quyết trọn vẹn và xuất sắc mục tiêu khoa học đề ra, với các đóng góp cốt lõi:

  1. Làm sáng tỏ bản chất cơ học: Chứng minh cơ chế hấp thụ năng lượng nổ vượt trội của vật liệu có cấu trúc lỗ rỗng dựa trên thềm biến dạng dẻo $\sigma_{pl}$ và công biến dạng nén $W_{vol}$.
  2. Chuẩn hóa công cụ mô phỏng số: Thiết lập thành công quy trình mô phỏng tương tác nổ - kết cấu trong LS-DYNA sử dụng kỹ thuật ALE và mô hình vật liệu MAT_063/MAT_084 có độ chính xác cao, được xác thực bởi thực nghiệm nổ hiện trường.
  3. Đề xuất giải pháp cấu tạo tối ưu: Khẳng định tính ưu việt của kết cấu bảo vệ dạng sandwich nhiều lớp (Thép - Xốp - Thép) và cấu trúc phân đoạn module trong việc triệt tiêu hiện tượng chấn sụp, giảm gia tốc đỉnh truyền vào kết cấu lên đến $50% - 65%$.
  4. Mở ra các nhánh nghiên cứu mới: Tạo tiền đề cho các hướng nghiên cứu sâu về nhôm bọt, xốp kim loại và kết cấu hy sinh tiêu tán năng lượng trong kỹ thuật công trình quân sự và dân dụng.
  5. Giá trị thực tiễn và chuyển giao: Cung cấp giải pháp kỹ thuật khả thi, chi phí hợp lý, dễ dàng lắp dựng để gia cường nâng cao khả năng sống còn cho các công trình trọng yếu quốc gia trước nguy cơ tấn công nổ và khủng bố.