Tổng quan về luận án

Nghiên cứu về độ bền vững và khả năng chống sụp đổ của kết cấu công trình trước các biến cố thảm họa là một trong những nhánh trọng yếu của kỹ thuật công trình hiện đại. Luận án tiến sĩ kỹ thuật với đề tài "Nghiên cứu sụp đổ lũy tiến của khung bê tông cốt thép toàn khối chịu tác dụng nổ và đề xuất một số giải pháp kháng sập" do nghiên cứu sinh Phan Thành Trung thực hiệ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, Mã số: 9 58 02 06; dưới sự hướng dẫn khoa học của TSKH. Nguyễn Văn Hợi và TS. Nguyễn Quốc Bảo, bảo vệ năm 2021) là công trình nghiên cứu tiên phong tại Việt Nam tích hợp cơ học nổ, cơ học phá hủy phi tuyến và phân tích động lực học kết cấu không gian 3D nhằm giải quyết bài toán chống sập đổ dây chuyền cho công trình an ninh - quốc phòng và dân dụng trọng yếu.

                  +----------------------------------------------+
                  | TẢI TRỌNG NỔ (BLAST LOAD: TNT CONTACT/CLOSE) |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |  Hư hỏng cục bộ ban đầu (Cột A1, A3)         |
                  |  - Phá hoại cắt, vỡ nát bê tông, nén nổ      |
                  +----------------------+-----------------------+
                                         |
                                         v
                  +----------------------------------------------+
                  |  Cơ chế truyền tải thay thế (ALPM)           |
                  |  - Hình thành khớp dẻo liên tiếp             |
                  |  - Tái phân bố nội lực động học              |
                  +----------------------+-----------------------+
                                         |
                    +--------------------+--------------------+
                    |                                         |
                    v                                         v
+---------------------------------------+ +---------------------------------------+
|  KHÔNG GIA CƯỜNG (UNRETROFITTED)      | |  GIA CƯỜNG KHÁNG SẬP (RETROFITTED)    |
|  - Mất ổn định dây chuyền             | |  - Cột phụ (Auxiliary Columns)        |
|  - Phá hoại khớp dẻo                  | |  - Bọc thép cột (Steel Jacketing)     |
|  - SỤP ĐỔ LŨY TIẾN TOÀN BỘ            | |  - Tường chắn sóng (Blast Wall)       |
|    (0.005s - 0.900s)                  | |  -> TRIỆT TIÊU SỤP ĐỔ LŨY TIẾN       |
+---------------------------------------+ +---------------------------------------+

Khoảng trống học thuật (research gap) mà luận án tập trung giải quyết xuất phát từ sự bất cập của các tiêu chuẩn thiết kế quốc tế hiện hành như GSA (2003) và DoD UFC 4-023-03. Các tiêu chuẩn này chủ yếu áp dụng phương pháp đường dẫn tải thay thế (Alternate Load Path Method - ALPM) theo cách tiếp cận "độc lập với hiểm họa" (threat-independent), giả định loại bỏ tức thời một cấu kiện chịu lực chính mà bỏ qua bản chất vật lý của sóng nổ, tốc độ biến dạng siêu cao ($\dot{\varepsilon} > 10^2 - 10^3\text{ s}^{-1}$), hiệu ứng dội sóng và quá trình phá hủy nén nổ - phá vỡ dòn cục bộ trước khi chuyển sang giai đoạn động học sập đổ toàn thể. Tại Việt Nam, các nghiên cứu trước đây chủ yếu dừng lại ở khung phẳng đàn dẻo lý tưởng hoặc bài toán tĩnh tương đương (Алмазов và Cao Duy Khôi, 2010; Vũ Văn Hoàng, 2012), chưa mô tả được đầy đủ ứng xử phi tuyến vật liệu kết hợp nứt vỡ, phá hủy dưới tải nổ phức tạp.

Luận án thiết lập 3 câu hỏi nghiên cứu và 3 giả thuyết khoa học tương ứng:

  1. Câu hỏi nghiên cứu 1 (RQ1): Bộ tham số thực nghiệm cho mô hình vật liệu phi tuyến phá hủy động học nào phản ánh chính xác nhất ứng xử của bê tông và cốt thép chế tạo tại Việt Nam dưới tải nổ xung kích?
    • Giả thuyết 1 (H1): Mô hình cơ học phá hủy Holmquist-Johnson-Cook (HJC) đối với bê tông và Johnson-Cook (JC) đối với cốt thép, khi được xác định thông qua thí nghiệm nén 3 trục và hiệu chỉnh bằng thực nghiệm nổ hiện trường, sẽ mô tả chính xác cả 3 pha: đàn hồi, tăng bền áp lực và phá hủy suy giảm độ cứng.
  2. Câu hỏi nghiên cứu 2 (RQ2): Cơ chế phá hoại cục bộ ban đầu do nổ tiếp xúc và nổ gần chuyển tiếp sang quá trình sụp đổ lũy tiến của khung không gian bê tông cốt thép (BTCT) toàn khối diễn ra như thế nào theo thang thời gian mili-giây?
    • Giả thuyết 2 (H2): Sự phá hủy dẻo tại các khớp dẻo và đứt gãy tiết diện chịu lực ở các cột biên (A1) hoặc cột giữa (A3) sẽ gây ra phản ứng dây chuyền làm mất độ cứng tổng thể của khung không gian nhanh hơn đáng kể so với giả định thiết kế truyền thống do có sự kết hợp của hiệu ứng nổ phá vỡ cục bộ và gia tải động thứ cấp.
  3. Câu hỏi nghiên cứu 3 (RQ3): Các giải pháp kết cấu nào có khả năng ngăn chặn hoặc giảm thiểu tối đa nguy cơ sụp đổ lũy tiến cho khung không gian BTCT chịu tải nổ?
    • Giả thuyết 3 (H3): Việc can thiệp bằng giải pháp trực tiếp (bọc thép cột, bổ sung cột phụ) và giải pháp kiềm chế sự kiện (tường chắn sóng nổ) sẽ triệt tiêu được đường dẫn phá hủy dây chuyền, bảo tồn khả năng chịu lực dư của hệ kết cấu.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng cơ học môi trường liên tục phi tuyến (Nonlinear Continuum Mechanics), lý thuyết phá hủy động lực học (Dynamic Fracture Mechanics), phương trình trạng thái vật liệu nổ Jones-Wilkins-Lee (JWL) và phương pháp tích phân thời gian tường minh (Explicit Time Integration) trong môi trường phần mềm ABAQUS/Explicit.

Quy mô và phạm vi nghiên cứu bao gồm: Thử nghiệm cơ tính vật liệu trong phòng thí nghiệm (nén đơn trục, kéo ép chẻ, nén 3 trục); 05 tổ hợp thử nghiệm nổ thực địa với lượng nổ TNT tiêu chuẩn đo đạc bằng hệ thống thu thập dữ liệu động NI SCXI-1000DC và cảm biến biến dạng chuyên dụng KC-60-120-A1-11; mô phỏng số chi tiết khung không gian BTCT toàn khối nhiều tầng nhiều nhịp chịu các kịch bản nổ tiếp xúc và nổ gần tại các vị trí cột xung yếu (cột góc A1, cột biên A3) trong khoảng thời gian từ $0,005\text{ s}$ đến $0,9\text{ s}$. Ý nghĩa khoa học và thực tiễn của luận án là cung cấp bộ thông số cơ học tin cậy cho vật liệu trong nước, hoàn thiện phương pháp luận phân tích kháng sập động lực học phi tuyến và đưa ra các giải pháp thiết kế định lượng cho công trình công sự, sở chỉ huy quân sự và hạ tầng phòng thủ dân sự.

Literature Review và Positioning

Lịch sử nghiên cứu hiện tượng sụp đổ lũy tiến (Progressive Collapse) được đánh dấu bởi các sự cố thảm họa mang tính bước ngoặt:

  • Vụ sập một góc 22 tầng của chung cư Ronan Point (Luân Đôn, 1968) do nổ khí gas tại tầng 18 làm bay tường bê tông đúc sẵn, kéo theo phản ứng dây chuyền dạng "trò chơi xếp hình" (house of cards) làm sập toàn bộ các tầng ngầm bên dưới.
  • Vụ đánh bom xe tải chở 1,8 tấn thuốc nổ TNT cách mặt phía bắc tòa nhà Alfred P. Murrah (Oklahoma, 1995) 4,9m, phá hủy 3 cột biên dẫn đến sự sụp đổ của gần một nửa tòa nhà.
  • Vụ tấn công khủng bố ngày 11/9/2001 vào tháp đôi Trung tâm Thương mại Thế giới (WTC, New York), nơi va chạm của máy bay Boeing 767 kết hợp cháy nhiên liệu diện rộng làm mất ổn định hệ giằng và cột chu vi, gây thảm họa sập đổ toàn phần làm 2.749 người thiệt mạng.
                                  TIẾN TRÌNH TIÊU CHUẨN KHÁNG SẬP
               
[1968: Thảm họa Ronan Point]  --->  [1970: Tiêu chuẩn Anh & NBC Canada đưa quy định SĐLT]
                                           |
                                           v
[1995: Tòa nhà Alfred P. Murrah] --->  [2000: Bộ tiêu chuẩn Quốc tế IBC & ASCE 7]
                                           |
                                           v
[2001: Thảm họa WTC 11/9]     --->  [2003: Hướng dẫn GSA & 2005: Tiêu chuẩn DoD UFC 4-023-03]
                                           |
                                           v
                              [2021: Luận án TS. Phan Thành Trung]
                              - Khung không gian 3D BTCT toàn khối
                              - Tải nổ thực tế (JWL + HJC + Explicit)
                              - Định lượng giải pháp gia cường kháng sập

Về mặt lý thuyết và quy chuẩn, sau sự cố Ronan Point, Anh (1970) và Canada (1970) là những quốc gia đầu tiên đưa yêu cầu kháng sập vào quy chuẩn xây dựng. Đến năm 1981, bộ tiêu chuẩn BOCA của Mỹ bắt đầu tích hợp các điều khoản này, và sau đó được hợp nhất trong Bộ tiêu chuẩn Xây dựng Quốc tế IBC (2000) cùng ASCE 7. Theo quy chuẩn ASCE 7 và IBC (2000):

"Nhà cửa và các kết cấu khác nhau được thiết kế sao cho để nếu xảy ra phá hoại cục bộ, thì toàn bộ hệ thống kết cấu còn lại vẫn bền vững và không bị phá hoại đến mức không cân xứng so với những thiệt hại cục bộ ban đầu."

Ủy ban An ninh Liên ngành Mỹ (ISC, 2001) cũng quy định:

"Ở mức tối thiểu, tất cả các bộ phận mới cần được thiết kế đảm bảo sao cho nếu tổn thất của một cột trong tầng một tại chu vi tòa nhà sẽ không xảy ra sự sụp đổ lũy tiến."

Trong y văn thế giới, tồn tại hai luồng quan điểm và tranh luận học thuật sâu sắc:

  1. Luồng quan điểm thiết kế gián tiếp và phương pháp đường dẫn tải thay thế đơn giản hóa: Marjanishvili & Agnew (2006), Powell (2005), Izzuddin et al. (2007) đề xuất các phương pháp phân tích tĩnh phi tuyến (pushover) hoặc động phi tuyến tương đương dựa trên khái niệm bảo toàn năng lượng và khớp dẻo tập trung. Nhóm nghiên cứu này ủng hộ việc đơn giản hóa bài toán bằng cách loại bỏ trực tiếp cột để tính toán phân bố lại nội lực nhằm phục vụ thiết kế nhanh. Tuy nhiên, nhược điểm là bỏ qua hoàn toàn tải trọng bất thường ban đầu và hiệu ứng phá hoại dẻo kèm vỡ nát vật liệu ở vùng lân cận.
  2. Luồng quan điểm mô phỏng số chi tiết đa tỷ lệ và phân tích cơ học phá hủy nổ: Kwasniewski (2009) thực hiện mô phỏng khung thép hỗn hợp 8 tầng Cardington với 1,08 triệu phần tử, mất 19 ngày tính toán trên 60 bộ vi xử lý song song. Sadek et al. (2008) tại NIST và Yu et al. (2006) sử dụng LS-DYNA mô hình hóa hệ sàn liên hợp với mô hình bê tông phi tuyến nhằm làm rõ vai trò của liên kết kháng uốn và hiệu ứng màng (catenary action). Kaewkulchai & Williamson (2004) chỉ ra rằng phân tích động phi tuyến có xét đến lực va chạm của cấu kiện rơi là bắt buộc vì phương pháp tĩnh không thể dự báo chính xác nhu cầu biến dạng cực hạn.

Tại Việt Nam, các công trình của Ngô Đức Tuấn (Đại học Melbourne) và Phạm Xuân Đạt (Đại học Xây dựng Hà Nội) chủ yếu tập trung vào thử nghiệm cấu kiện chịu nổ đơn lẻ hoặc ALPM thuần túy. Luận án của Phan Thành Trung định vị ở điểm giao thoa đột phá: không xem việc loại bỏ cấu kiện là một giả định trừu tượng mà mô hình hóa trực tiếp quá trình kích nổ lượng nổ TNT, tương tác sóng nổ - kết cấu, phá hủy nứt vỡ cục bộ cấu kiện BTCT và động học sụp đổ lũy tiến tiếp theo trên mô hình khung không gian 3D hoàn chỉnh.

Đó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 và hoàn thiện lý thuyết ứng xử phi tuyến của kết cấu bê tông cốt thép chịu tải trọng đặc biệt thông qua các đóng góp lý thuyết cốt lõi:

  • Hiệu chỉnh và xác lập bộ tham số mô hình phá hủy kết hợp Holmquist-Johnson-Cook (HJC, 1993) cho bê tông Việt Nam: Mô hình HJC nguyên bản được xây dựng cho các loại bê tông tiêu chuẩn quân sự phương Tây. Luận án đã mở rộng mô hình này cho bê tông sử dụng xi măng PC30 chế tạo trong nước bằng cách thiết lập thực nghiệm toàn diện 21 tham số vật liệu bao gồm phương trình trạng thái (EOS), hàm mặt dẻo tăng cứng áp lực và quy luật tích lũy biến phá hủy.
  • Chuyển dịch hệ hình (Paradigm shift) từ lý thuyết khớp dẻo truyền thống sang cơ học phá hoại dẻo - suy giảm độ cứng động học: Trong tính toán thiết kế truyền thống (Baker, 1956), khớp dẻo sau khi xuất hiện được coi là duy trì khả năng chịu lực uốn không đổi. Luận án chứng minh rằng dưới tác dụng nổ, khi biến dạng dẻo tương đương đạt giá trị cực hạn $\varepsilon_p \ge \varepsilon_{f}$, khớp dẻo bị phá hủy hoàn toàn (đứt gãy tiết diện, tách phần tử khỏi nút khung), dẫn đến sự suy giảm độ cứng đột ngột và kích hoạt phản ứng sụp đổ dây chuyền với tốc độ phá hoại nhanh gấp nhiều lần dự báo của lý thuyết dẻo cổ điển.
                           MÔ HÌNH TOÁN HỌC VẬT LIỆU HJC
                           
         Ứng suất chuẩn hóa (sigma*)               Phương trình trạng thái P(mu)
                    ^                                       ^
           Smax ----|--------+                              |        P_lock (Giai đoạn 3: Nén chặt)
                    |       /                               |         /
                    |      / (Mặt dẻo HJC)                  |        /
                    |     /                                 |  P_crush (Giai đoạn 2: Nứt vỡ)
                    |    /                                  |    /
                    |   /                                   |   / (Giai đoạn 1: Đàn hồi)
                    |  /                                    |  /
                    +--------------------> P*               +--------------------> mu
                       (Áp lực chuẩn hóa)                      (Độ biến dạng thể tích)

Mô hình bền HJC được xác định bởi phương trình ứng suất tương đương chuẩn hóa: $$\sigma^* = \left[ A(1 - D) + B P^{N} \right] \left[ 1 + C \ln(\dot{\varepsilon}^) \right] \le S_{max}$$

Biến phá hủy tích lũy $D$ ($0 \le D \le 1$) phát triển theo công thức: $$D = \sum \frac{\Delta \varepsilon_p + \Delta \mu_p}{\varepsilon_p^f + \mu_p^f}$$ trong đó biến dạng phá hủy dẻo dưới áp suất $P^$ được biểu diễn qua: $$\varepsilon_p^f + \mu_p^f = D_1 (P^ + T^*)^{D_2} \ge e_{fmin}$$

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

Luận án đã thiết lập một khung phân tích đa vật lý (Multiphysics Framework) tích hợp chặt chẽ 3 nền tảng lý thuyết:

  1. Lý thuyết thủy động lực học sóng nổ khí quyển: Sử dụng phương trình trạng thái Jones-Wilkins-Lee (JWL) để mô tả động học giải phóng năng lượng của lượng nổ TNT trong không gian: $$P = A_{tn} \left(1 - \frac{\omega}{R_1 V}\right) e^{-R_1 V} + B_{tn} \left(1 - \frac{\omega}{R_2 V}\right) e^{-R_2 V} + \frac{\omega E_{tn}}{V}$$
  2. Cơ học phá hủy môi trường liên tục (Continuum Damage Mechanics): Mô hình hóa sự suy giảm mô đun đàn hồi hiệu dụng $E_{eff} = E_0 (1 - D)$ và tính dòn, tính nén nổ của bê tông kết hợp với mô hình Johnson-Cook mô tả tính chảy dẻo, hóa bền biến dạng và phá hủy cắt của cốt thép.
  3. Lý thuyết phần tử hữu hạn phi tuyến động học: Tích phân phương trình chuyển động từng bước theo sơ đồ sai phân trung tâm tường minh (Explicit Central Difference): $$\mathbf{M} \ddot{\mathbf{u}}{(t)} = \mathbf{F}^{ext}{(t)} - \mathbf{I}^{int}_{(t)}$$ với điều kiện ổn định Courant-Friedrichs-Lewy: $\Delta t \le \min \left( \frac{L_e}{c_d} \right)$.

Điều kiện biên và phạm vi ứng dụng của khung phân tích được xác lập rõ ràng cho kết cấu khung không gian BTCT toàn khối chịu tải trọng nổ tiếp xúc và nổ gần, loại trừ các ảnh hưởng nhiệt do hỏa hoạn kéo dài và tải trọng va chạm thứ cấp quy mô lớn của mảnh văng vĩ mô từ các công trình lân cận.

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

+-----------------------------------------------------------------------------------+
|                           QUY TRÌNH NGHIÊN CỨU TOÀN DIỆN                          |
+-----------------------------------------------------------------------------------+
|                                                                                   |
|  [PHA 1: PHÒNG THÍ NGHIỆM]                                                        |
|  - Thí nghiệm nén đơn trục, kéo ép chẻ (Mẫu trụ & lập phương xi măng PC30)        |
|  - Thí nghiệm nén 3 trục (Triaxial Cell) ở các cấp áp lực buồng khác nhau         |
|  -> Xác định 21 tham số mô hình HJC: fc, T, A, B, N, C, Smax, D1, D2, efmin, K...|
|                                         |                                         |
|                                         v                                         |
|  [PHA 2: THỰC NGHIỆM NỔ HIỆN TRƯỜNG]                                              |
|  - 05 cấu kiện cột/dầm BTCT thực tế                                               |
|  - Thuốc nổ TNT + Kíp nổ điện số 8 + Máy điểm hỏa FD200                           |
|  - Đo động: Hệ thống NI SCXI-1000DC + Cảm biến strain gauge KC-60-120-A1-11      |
|  -> Thu thập lịch sử biến dạng LE33, đo kích thước phễu nổ, chấn sụp, vỡ nát      |
|                                         |                                         |
|                                         v                                         |
|  [PHA 3: MÔ PHỎNG HIỆU CHỈNH SỐ TRÊN ABAQUS/EXPLICIT]                             |
|  - Mô hình Lagrange phần tử 3D Solid C3D8R (bê tông) + Truss T3D2 (thép)          |
|  - Tương tác tiếp xúc nổ & nổ gần                                                 |
|  - So sánh đối chứng kết quả nứt vỡ & biến dạng (Sai số < 10%)                    |
|                                         |                                         |
|                                         v                                         |
|  [PHA 4: PHÂN TÍCH KHUNG KHÔNG GIAN VÀ ĐỀ XUẤT GIẢI PHÁP]                         |
|  - Kịch bản nổ cột biên A1, nổ cột giữa A3 (Tiếp xúc & Nổ gần)                    |
|  - Đánh giá quá trình sụp đổ lũy tiến (0.005s - 0.9s)                             |
|  - Khảo sát 3 giải pháp: Cột phụ, Bọc thép, Tường chắn sóng                       |
+-----------------------------------------------------------------------------------+

Thiết kế nghiên cứu

Luận án áp dụng thế giới quan bản thể luận hiện thực phản biện (Critical Realism) và nhận thức luận thực chứng thực nghiệm (Positivist-Empirical Paradigm). Nghiên cứu kết hợp đa phương pháp (Mixed-Method Quantitative Approach) bao gồm: Thực nghiệm vật liệu vi mô $\rightarrow$ Thực nghiệm hiện trường trung mô $\rightarrow$ Mô phỏng số động lực học phi tuyến vĩ mô. Thiết kế nghiên cứu đa cấp độ (Multi-level Design) cho phép kiểm chứng độc lập các phương trình trạng thái trước khi đưa vào bài toán kết cấu phức tạp.

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm và mô phỏng được chuẩn hóa nghiêm ngặt:

  1. Thí nghiệm trong phòng: Chế tạo các bộ mẫu bê tông mác chuẩn sử dụng xi măng PC30. Thực hiện thí nghiệm nén đơn trục trên máy nén thủy lực xác định cường độ nén $f_c$ và mô đun đàn hồi $E$; thí nghiệm kéo ép chẻ xác định độ bền kéo đơn trục $T$; thí nghiệm nén ba trục (Triaxial Compression Test) tại các mức áp lực hông khác nhau để xác định đường cong bao phá hủy, hệ số làm cứng áp lực chuẩn hóa $B$ và số mũ tăng cứng $N$.
  2. Thực nghiệm nổ ngoài bãi thử: Chế tạo 05 mẫu cấu kiện cột BTCT kích thước thực tế kèm gối tựa chuyên dụng. Bố trí lượng nổ TNT tiếp xúc bề mặt và nổ gần ở cự ly xác định. Hệ thống đo đạc động lực học gồm máy đo NI SCXI-1000DC tốc độ cao kết hợp đầu đo biến dạng điện trở KC-60-120-A1-11 gắn trực tiếp tại các vị trí xung yếu trên bề mặt bê tông và cốt thép; kích nổ chuẩn xác bằng kíp điện số 8 và máy điểm hỏa FD200.
  3. Hiệu chỉnh mô hình số: Thiết lập mô hình phần tử hữu hạn 3D tương đương trên ABAQUS. Kết quả thực nghiệm và số được đối chứng chéo về hình thái nứt vỡ, độ sâu chấn sụp, chiều dài phễu phá hủy và đáp ứng biến dạng thời gian thực.
Thông số hình học & Phá hủy Thí nghiệm thực tế (Mẫu 2 & 3) Mô phỏng số ABAQUS Mức độ sai lệch (%)
Nổ tiếp xúc: Chiều dài vùng phá hủy $32,5\text{ cm}$ $30,8\text{ cm}$ $5,23%$
Nổ tiếp xúc: Chiều sâu phá hoại nén nổ $12,0\text{ cm}$ $11,4\text{ cm}$ $5,00%$
Nổ gần: Chiều dài vùng chấn sụp mặt trước $45,0\text{ cm}$ $42,6\text{ cm}$ $5,33%$
Nổ gần: Chiều cao vùng vỡ nát mặt sau $28,0\text{ cm}$ $26,1\text{ cm}$ $6,78%$
Đỉnh biến dạng dọc trục ($LE_{33}$) tại Điểm 1 $1450,\mu\varepsilon$ $1380,\mu\varepsilon$ $4,83%$

Độ sai lệch giữa mô phỏng số và thực nghiệm hiện trường đều dưới $7%$, khẳng định độ tin cậy và tính hợp lệ (construct validity and internal validity) của mô hình tính toán.

Data và phân tích

Trong phân tích số trên ABAQUS/Explicit:

  • Vật liệu bê tông được rời rạc hóa bằng phần tử khối 8 nút suy giảm tích phân (C3D8R) với kỹ thuật kiểm soát chế độ biến dạng ảo (Hourglass Control). Kích thước lưới phần tử được tinh chỉnh cục bộ tại vùng nổ ($10\text{ mm} - 15\text{ mm}$) và mở rộng dần ra xa ($25\text{ mm} - 50\text{ mm}$) để đảm bảo thu bắt chính xác gradient ứng suất xung kích.
  • Cốt thép được mô hình hóa bằng phần tử thanh chống/kéo 2 nút không gian (T3D2), liên kết hoàn toàn (Embedded Region Constraint) với bê tông khối.
  • Thuốc nổ TNT và không khí được mô hình hóa bằng lưới Euler hoặc tải trọng áp suất mặt nổ định sẵn (Abaqus Blast Loading Interface).
  • Quá trình phân tích động phi tuyến kiểm soát chặt chẽ biến dạng logarit $LE_{33}$, ứng suất Von Mises và biến phá hủy $D$ tại các phần tử nút dầm - cột then chốt (như các phần tử PT 739, PT 793, PT 847).

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

               ĐÁP ỨNG BIẾN DẠNG ĐỘNG LE33 VÀ CHUYỂN VỊ NÚT KHUNG THEO THỜI GIAN
               
   Biến dạng LE33 (x10^-3)
      ^
 3.5 -|                                     /---\ (Không gia cường: Mất ổn định & Sập đổ)
 3.0 -|                                    /     \
 2.5 -|                      /------------/       \====> Đứt gãy hoàn toàn (t = 0.85s)
 2.0 -|                     /
 1.5 -|          /---------/ (Hình thành khớp dẻo liên tiếp)
 1.0 -|         /
 0.5 -|--------/ (Có bọc thép / Cột phụ: Dao động ổn định tắt dần <= 0.0012)
      +--------+---------+---------+---------+---------+---------+---------> Thời gian (s)
      0      0.15      0.30      0.45      0.60      0.75      0.90

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

  1. Cơ chế lan truyền phá hủy hai pha dưới tác dụng nổ tiếp xúc: Khi khối thuốc nổ TNT phát nổ tiếp xúc tại chân cột A1 hoặc cột A3, quá trình phá hủy diễn ra cực kỳ nhanh: trong pha sóng xung kích ($0 - 0,007\text{ s}$), sóng nén nổ cường độ hàng chục GPa gây vỡ vụn bê tông cục bộ và đứt gãy cốt thép tại chỗ. Ngay sau đó ($0,007\text{ s} - 0,45\text{ s}$), cột bị loại bỏ hoàn toàn khả năng chịu lực, làm phát sinh phản lực động và chuyển vị thẳng đứng đột ngột tại nút khung phía trên, kích hoạt chuyển đổi cơ chế truyền lực uốn sang cơ chế lực màng/lực căng kéo dây chuyền trên các dầm liền kề.
  2. Quy luật sụp đổ lũy tiến theo thời gian của khung không gian (Kịch bản Cột A1 và A3):
    • Tại cột góc A1: Sự mất liên kết góc gây mất ổn định xoắn không gian 2 phương. Từ $t = 0,55\text{ s}$ đến $t = 0,90\text{ s}$, các phần tử dầm sàn tầng 1 và tầng 2 lần lượt xuất hiện khớp dẻo thứ cấp tại các tiết diện gối tựa; biến dạng $LE_{33}$ tại các phần tử 739, 793 vượt ngưỡng phá hủy dẻo dẫn đến sự sụp đổ dây chuyền toàn bộ góc biên công trình.
    • Tại cột biên giữa A3: Nhịp dầm phía trên cột A3 chịu độ võng tăng gấp 4,2 lần tải trọng tĩnh. Tại thời điểm $t = 0,85\text{ s}$, toàn bộ 3 nhịp dầm trục ngang bị kéo đứt gãy khỏi nút khung, gây sụp đổ lũy tiến 4 gian sàn liền kề.
  3. Hiện tượng sóng phản xạ kéo khi chịu tải nổ gần: Trong trường hợp nổ gần ngoài mặt nhà, sóng nổ tác dụng vào mặt trước cột tạo ra ứng suất nén; khi truyền đến mặt sau (mặt tự do), sóng bị phản xạ thành sóng kéo. Do cường độ chịu kéo của bê tông chỉ bằng khoảng $8 - 10%$ cường độ nén ($T \approx 1,85\text{ MPa}$ so với $f_c \approx 23,5\text{ MPa}$), hiện tượng bong tróc mảng lớn (spalling) xảy ra dữ dội ở mặt sau cột, làm suy giảm $65%$ diện tích tiết diện làm việc của cột trước khi tải trọng thẳng đứng kịp phân bố lại.
  4. Hiệu quả định lượng vượt trội của các giải pháp kháng sập:
    • Giải pháp bổ sung cột phụ: Khi bố trí cột phụ lân cận cột A3, độ võng động cực đại của dầm sàn giảm từ $385\text{ mm}$ xuống còn $42\text{ mm}$ (giảm $89,1%$), triệt tiêu hoàn toàn khả năng hình thành khớp dẻo phá hoại dây chuyền.
    • Giải pháp bọc thép cột: Ống thép bọc ngoài tạo hiệu ứng bó hông (confinement effect) 3 trục làm tăng độ bền nén bê tông lõi lên $2,4$ lần, ngăn chặn hiện tượng vỡ nát bê tông dưới nổ tiếp xúc; biến dạng dẻo $LE_{33}$ duy trì dưới ngưỡng $0,0015$, kết cấu khung hoàn toàn ổn định sau nổ ($t = 0,24\text{ s}$).
    • Giải pháp tường chắn sóng nổ: Tường chắn làm chệch hướng dòng khí nổ áp suất cao, giảm áp lực sóng nổ đỉnh tác dụng trực tiếp lên khung công trình từ $48%$ đến $62%$, chuyển phá hoại từ mức sụp đổ cục bộ sang trạng thái hư hại nhẹ có thể sửa chữa.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm và cơ sở toán học khẳng định tính tất yếu của việc tích hợp cơ học nổ xung kích vào mô hình phân tích sụp đổ lũy tiến kết cấu, thay thế cho các giả định tĩnh/động tuyến tính đơn giản hóa đang tồn tại trong ASCE 7 và GSA.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình thực nghiệm xác định bộ 21 thông số HJC và phương pháp mô phỏng số ABAQUS/Explicit có thể chuyển giao áp dụng cho nhiều dạng kết cấu công trình ngầm, công trình chịu va chạm tốc độ cao và công sự quốc phòng.
  • Về ứng dụng thực tiễn: Cung cấp hướng dẫn kỹ thuật cụ thể về kích thước, chiều dày bản thép bọc cột và quy cách bố trí cột phụ dự phòng cho các kỹ sư kết cấu khi thiết kế các tòa nhà trung tâm hành chính, đại sứ quán và sở chỉ huy quân sự.
  • Về chính sách và quy chuẩn: Đóng góp cơ sở khoa học trực tiếp cho Bộ Xây dựng và Bộ Quốc phòng trong việc xây dựng Tiêu chuẩn Quốc gia (TCVN) chuyên biệt về thiết kế công trình kháng nổ và kháng sập lũy tiến tại 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 4 hạn chế khoa học:

  1. Chưa xét đến tác động đồng thời của nhiệt độ cháy (Fire-Blast Interaction): Trong các vụ nổ thực tế (như sự cố WTC 2001), tải trọng nổ thường kéo theo cháy lớn làm suy giảm nhanh giới hạn chảy của thép và cường độ bê tông theo thời gian nhiệt độ.
  2. Giản lược hóa tải trọng va chạm của các cấu kiện bị phá hủy rơi xuống: Mô hình tập trung vào sự hình thành và đứt gãy khớp dẻo của khung không gian nhưng chưa mô phỏng chi tiết động lượng va đập của các tấm sàn bê tông rơi tự do xuống các tầng sàn bên dưới.
  3. Mới dừng lại ở kết cấu bê tông cốt thép toàn khối thông thường: Chưa khảo sát ứng xử kháng sập của các loại vật liệu tiên tiến mới như bê tông cường độ siêu cao (UHPC), bê tông tự lèn gia cường sợi thép (SFRC) hay kết cấu lắp ghép có khớp nối bán cứng.
  4. Quy mô nổ hiện trường giới hạn ở cấp cấu kiện: Do điều kiện an toàn bãi thử và trang thiết bị, thử nghiệm nổ thực địa mới tiến hành trên 05 cấu kiện cột/dầm đơn lẻ chứ chưa thể kích nổ trên mô hình khung không gian 3D toàn phần quy mô 1:1.

Từ các hạn chế trên, hướng nghiên cứu tiếp theo (Future Research Agenda) được xác lập gồm:

  • Mở rộng bài toán phân tích sụp đổ lũy tiến chịu tác động đa hiểm họa liên hoàn: Nổ $\rightarrow$ Cháy nhiệt độ cao $\rightarrow$ Sụp đổ dây chuyền (Blast-Fire-Collapse Multi-hazard Coupling).
  • Ứng dụng vật liệu composite polyme cốt sợi (CFRP/AFRP) và ống thép nhồi bê tông cường độ siêu cao (CFST/UHPC) để tăng cường sức kháng nổ cục bộ.
  • Xây dựng mô hình phân tích ngẫu nhiên (Probabilistic Progressive Collapse Analysis) có xét đến tính bất định của vị trí nổ, lượng nổ TNT và khuyết tật hình học của vật liệu.

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

  • Tác động học thuật (Academic Impact): Luận án là tài liệu tham khảo nền tảng trong lĩnh vực công trình chịu tải trọng đặc biệt tại Việt Nam, mở ra hướng ứng dụng cơ học nổ tính toán (Computational Blast Mechanics) trong các chương trình đào tạo tiến sĩ và thạc sĩ chuyên ngành Xây dựng Công trình Đặc biệt và Công trình Quốc phòng.
  • Chuyển đổi ngành công nghiệp xây dựng: Cung cấp giải pháp công nghệ bọc thép kết hợp cột phụ giúp các doanh nghiệp tư vấn thiết kế và nhà thầu xây dựng tối ưu hóa chi phí gia cường bảo vệ công trình trọng điểm mà không cần tăng tiết diện kết cấu quá mức.
  • Ảnh hưởng đến an ninh - quốc phòng: Trực tiếp phục vụ nhiệm vụ thiết kế, thẩm tra và nâng cấp độ an toàn cho hệ thống sở chỉ huy các cấp, công trình chiến đấu, kho tàng vũ khí đạn dược và công trình phòng thủ bờ biển trước các loại vũ khí công nghệ cao.
  • Giá trị xã hội và quốc tế: Nâng cao năng lực chống chịu thảm họa của các đô thị hiện đại, giảm thiểu thương vong và thiệt hại kinh tế khi xảy ra các sự cố nổ công nghiệp hoặc khủng bố, tiếp cận trình độ nghiên cứu của các viện nghiên cứu quân sự và công nghệ hàng đầu thế giới (như US Army Corps of Engineers, ERDC).

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

  • Nghiên cứu sinh và học giả: Tiếp cận quy trình thực nghiệm nén 3 trục chuẩn mực và phương pháp số hóa bộ tham số HJC/JC trong ABAQUS/Explicit để phát triển các đề tài nghiên cứu chuyên sâu về cơ học phá hủy động.
  • Các viện nghiên cứu và trường đại học: Sử dụng làm tài liệu giảng dạy chuyên đề cao học Kỹ thuật Công trình Đặc biệt, Kỹ thuật Công trình Quốc phòng, Động lực học Công trình.
  • Kỹ sư kết cấu và chuyên gia R&D công trình: Sở hữu công cụ tính toán định lượng và cẩm nang giải pháp gia cố (bọc thép cột, bố trí cột phụ, thiết kế tường chắn sóng) có thể áp dụng trực tiếp vào dự án thực tế.
  • Cơ quan quản lý nhà nước và Bộ Quốc phòng: Có luận cứ khoa học thực chứng để ban hành quy chuẩn kỹ thuật quốc gia về đánh giá nguy cơ và thiết kế phòng chống sụp đổ lũy tiến cho các công trình cấp đặc biệt.

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 mô hình phá hủy cơ học liên tục Holmquist-Johnson-Cook (HJC) để xác lập tường minh hệ thống 21 tham số động học cho bê tông Việt Nam mác PC30, kết nối thành công giữa phương trình trạng thái phi tuyến thủy động lực học ($P - \mu$) với quy luật tích lũy biến phá hủy dẻo ($D$) dưới tốc độ biến dạng cực cao ($\dot{\varepsilon} > 10^2\text{ s}^{-1}$).

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây là gì?

So với nghiên cứu của Sadek et al. (2008) tại NIST và Kwasniewski (2009) vốn chỉ áp dụng phương pháp ALPM loại bỏ cột tức thời bằng thuật toán số thuần túy, luận án đã kết hợp trực tiếp mô hình sóng nổ thực (JWL EOS) với thử nghiệm nổ thực địa đo bằng hệ thống thu thập động NI SCXI-1000DC, cho phép mô hình hóa trọn vẹn cả pha phá hủy cục bộ nén nổ ban đầu và pha sụp đổ lũy tiến tổng thể trên khung không gian 3D.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm và mô phỏng là gì?

Phát hiện bất ngờ nhất là hiện tượng phá hủy do sóng phản xạ kéo (spalling) ở mặt sau của cấu kiện khi chịu nổ gần còn nguy hiểm hơn sóng nén ở mặt trước, làm mất hơn $60%$ tiết diện làm việc của cột chỉ trong vòng chưa đầy $0,016\text{ s}$, khiến cột bị mất ổn định chịu lực đứng trước cả khi khớp dẻo uốn hình thành đầy đủ.

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

Có. Luận án cung cấp chi tiết trong Chương 2, Chương 3 và Phần Phụ lục toàn bộ cấp phối bê tông, sơ đồ đấu nối cảm biến KC-60-120-A1-11, thông số máy đo NI SCXI, bảng 21 tham số vật liệu HJC và mã nguồn thiết lập bài toán trong phần mềm ABAQUS/Explicit, cho phép các nhà nghiên cứu độc lập tái lập hoàn toàn kết quả mô phỏng và thực nghiệm.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm tập trung vào 3 trọng tâm: (1) Phát triển mô hình tích hợp Đa hiểm họa Nổ - Cháy - Va đập cơ học; (2) Ứng dụng vật liệu composite tự phục hồi và bê tông siêu tính năng UHPC trong kết cấu kháng sập; (3) Số hóa quy chuẩn kháng sập lũy tiến thành bộ tiêu chuẩn thiết kế bắt buộc cho công trình đặc biệt tại Việt Nam.

Kết luận

Luận án tiến sĩ của tác giả Phan Thành Trung đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra với 6 đóng góp khoa học và thực tiễn nổi bật:

  1. Xác lập thành công hệ thống tham số thực nghiệm toàn diện của mô hình vật liệu bê tông HJC và cốt thép JC phù hợp với điều kiện chế tạo tại Việt Nam, được kiểm chứng với độ tin cậy cao qua thí nghiệm nén 3 trục và thử nghiệm nổ hiện trường (sai số $< 7%$).
  2. Làm sáng tỏ bản chất cơ học và thang thời gian mili-giây của quá trình chuyển tiếp từ phá hoại cục bộ sang sụp đổ lũy tiến trong khung không gian BTCT toàn khối dưới tác dụng của tải trọng nổ tiếp xúc và nổ gần.
  3. Chỉ ra quy luật phân bố lại nội lực động học, sự hình thành và phá hủy của các khớp dẻo thứ cấp dẫn đến sự mất ổn định không gian khi mất các vị trí cột xung yếu (cột góc A1, cột giữa A3).
  4. Đề xuất và định lượng hóa hiệu quả của 3 giải pháp công nghệ kháng sập: bổ sung cột phụ (giảm võng $89,1%$), bọc thép cột (tăng khả năng chịu lực nén nổ gấp 2,4 lần) và xây dựng tường chắn sóng nổ (giảm áp lực nổ tác dụng từ $48 - 62%$).
  5. Đóng góp bước tiến quan trọng về phương pháp luận nghiên cứu: xây dựng quy trình kết hợp chặt chẽ giữa cơ sở lý thuyết, thử nghiệm hiện trường đo động tốc độ cao và mô phỏng số 3D tiên tiến trên ABAQUS/Explicit.
  6. Đặt nền móng lý luận và thực nghiệm vững chắc cho hệ thống tiêu chuẩn thiết kế kháng chấn, kháng nổ và chống sụp đổ lũy tiến của ngành xây dựng công trình đặc biệt và công trình quốc phòng tại Việt Nam trong giai đoạn phát triển mới.