Tổng quan về luận án

Nghiên cứu của tác giả Phạm Thị Ngọc Thu (2021) với đề tài "Đánh giá khả năng chịu lửa của các cấu kiện thép chịu lực được bọc bảo vệ ứng dụng cho các công trình nhà tại Việt Nam" thuộc chuyên ngành Kỹ thuật xây dựng (Mã số: 9580201), do GS. Phạm Văn Hội hướng dẫn tại Trường Đại học Xây dựng Hà Nội, đã giải quyết một vấn đề khoa học then chốt trong kỹ thuật an toàn cháy công trình hiện đại.

Bối cảnh thực tiễn cho thấy các thảm họa cháy nổ công trình mang tính lịch sử trên thế giới như sự sụp đổ của hai tòa tháp Trung tâm Thương mại Thế giới (WTC, New York, 2001) chỉ sau 56 phút và 86 phút phơi lửa ở nhiệt độ trung tâm vùng nổ đạt $1.500^\circ\text{C}$, vụ hỏa hoạn kéo dài 17 giờ tại Trụ sở làm việc 50 tầng ở Caracas (Venezuela, 2004), hay vụ cháy tòa nhà văn phòng 32 tầng Windsor Building tại Madrid (Tây Ban Nha, 2004) làm sụp đổ 6 tầng nhà do nhiệt độ đạt tới $800^\circ\text{C}$ đã minh chứng cho sự hiểm họa khôn lường của lửa đối với hệ kết cấu chịu lực chính. Tại Việt Nam, số liệu thống kê cho thấy trung bình mỗi năm xảy ra khoảng 2.200 vụ cháy, cướp đi sinh mạng của 65 người, làm bị thương 180 người và gây thiệt hại vật chất ước tính trên 745 tỷ đồng mỗi năm. Các vụ hỏa hoạn nghiêm trọng như chợ Đồng Xuân (1995), Trung tâm Thương mại Quốc tế ITC TP. Hồ Chí Minh (2002, 60 người chết, 70 người bị thương) và Trung tâm Thương mại Hải Dương (2013, thiệt hại hơn 500 tỷ đồng) đã đặt ra yêu cầu cấp bách về việc nâng cao giới hạn chịu lửa cho kết cấu công trình.

Khoảng trống nghiên cứu (Research Gap) xuất phát từ thực trạng hệ thống quy chuẩn, tiêu chuẩn hiện hành của Việt Nam (TCVN 2622:1995, TCVN 6160:1996, TCVN 9311:2012, QCVN 06:2020/BXD). Các văn bản này chủ yếu tiếp cận theo phương pháp quy định danh định (prescriptive approach), dựa vào việc tra cứu bậc chịu lửa danh định và catalog của nhà sản xuất vật liệu cách nhiệt mà chưa cung cấp một quy trình tính toán phân tích ứng xử cơ - nhiệt phi tuyến thực tế của kết cấu thép trong điều kiện chịu lửa. Việc thiếu vắng công cụ giải tích và mô phỏng chuyên sâu khiến các kỹ sư gặp nhiều bế tắc khi thiết kế các cấu kiện thép ngoài danh mục thực nghiệm sẵn có.

Nghiên cứu tập trung giải quyết 3 câu hỏi nghiên cứu và kiểm chứng 3 giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Quy luật lan truyền nhiệt ba chiều trong cấu kiện thép chữ I không bọc và được bọc bảo vệ (bằng vữa chống cháy hoặc tấm thạch cao chống cháy) diễn biến như thế nào theo thời gian chịu tác động của đám cháy danh nghĩa ISO 834?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để xây dựng quy trình tính toán đơn giản hóa (SDM - Simplified Design Method) có khả năng xác định chính xác khả năng chịu uốn tới hạn ($M_{fi,t,Rd}$) và khả năng chịu nén tới hạn ($N_{fi,t,Rd}$) của cấu kiện dầm, cột thép được bọc bảo vệ?
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ tin cậy và sai số tương đối giữa phương pháp giải tích đơn giản hóa SDM so với mô phỏng phần tử hữu hạn phi tuyến nâng cao trên phần mềm ANSYS Workbench là bao nhiêu?
  • Giả thuyết nghiên cứu 1 (H1): Sự phân bố nhiệt độ không đồng đều trên tiết diện cấu kiện dầm và cột thép gây ra hiện tượng trôi dạt trục trung hòa dẻo và phân phối lại nội lực phức tạp, ảnh hưởng quyết định đến khả năng chịu lực tới hạn.
  • Giả thuyết nghiên cứu 2 (H2): Giải pháp bọc bảo vệ dạng hộp bằng tấm thạch cao chống cháy và bọc chu vi bằng vữa khoáng chống cháy giúp duy trì nhiệt độ thép dưới ngưỡng tới hạn ($\theta_{cr,d} \approx 500 - 600^\circ\text{C}$), bảo toàn khả năng chịu lực của cấu kiện trong thời gian tiêu chuẩn từ 30 đến 120 phút.
  • Giả thuyết nghiên cứu 3 (H3): Phương pháp đơn giản hóa SDM dựa trên phân vùng nhiệt độ cung cấp kết quả tính toán có độ an toàn tin cậy cao, chênh lệch dưới 15% so với mô phỏng số chi tiết đa trường vật lý.

Khung lý thuyết của nghiên cứu được thiết lập dựa trên sự giao thoa giữa Lý thuyết dẫn nhiệt Fourier trong không gian ba chiều, Lý thuyết phần tử hữu hạn (FEM) phi tuyến hình học và phi tuyến vật liệu, cùng Khung tiêu chuẩn thiết kế kết cấu thép chịu lửa Châu Âu (Eurocode EN 1991-1-2:2002 và EN 1993-1-2:2005). Phạm vi nghiên cứu giới hạn ở các cấu kiện dầm thép chữ I chịu uốn theo điều kiện bền chịu lửa 3 mặt (liên kết sàn bê tông phía trên) và cột thép chữ I chịu nén/uốn theo điều kiện ổn định tổng thể chịu lửa 4 mặt, ứng dụng hai dòng vật liệu bọc chủ đạo tại Việt Nam là thạch cao chống cháy (hệ số dẫn nhiệt $\lambda = 0,2 - 0,25 \text{ W/m}^\circ\text{C}$) và vữa chống cháy ($\lambda = 0,15 - 0,25 \text{ W/mK}$) trong khung thời gian khảo sát từ 30 đến 120 phút.

Literature Review và Positioning

Lịch sử nghiên cứu ứng xử của kết cấu xây dựng trong điều kiện chịu hỏa hoạn đã hình thành các dòng chảy học thuật lớn từ những năm 1970. Tại Thụy Điển, O. Pettersson và S. Magnusson (Đại học Lund) là những học giả tiên phong xây dựng nền móng về khái niệm bậc chịu lửa công trình và thực hiện hơn 100 chương trình thử nghiệm chịu lửa trên các mẫu dầm, cột thép tại Viện Nghiên cứu Kỹ thuật SP Boras (1980s). Đến năm 1990, Viện Nghiên cứu Xây dựng Anh (BRE) phối hợp thực hiện thí nghiệm chịu lửa trên mô hình nhà thép nguyên hình quy mô 8 tầng tại Cardington (Vương quốc Anh), tạo nên bước ngoặt chuyển dịch từ phân tích cấu kiện đơn lẻ sang phân tích sự làm việc tổng thể của khung không gian đa tầng khi chịu cháy.

Các hệ thống tiêu chuẩn quốc tế lần lượt ra đời định hình phương pháp luận thiết kế kết cấu chịu lửa: BS 5950 Part 8 (Anh, 1990), ASCE/SEI/SFPE 29-05 (Hoa Kỳ) được củng cố bởi các công trình nghiên cứu của Long T. McAllister và John L. Hurley tại Viện Tiêu chuẩn và Kỹ thuật Quốc gia Hoa Kỳ (NIST, 2005), và đặc biệt là Bộ tiêu chuẩn Châu Âu Eurocode 3 (EN 1993-1-2:2005). Song song đó, các nghiên cứu giải tích và thực nghiệm chuyên sâu của Franssen (1995), Janns (1988), Kang-Hai Tan (2003), Marek Lokomski (2015), cùng nhóm nghiên cứu Anil Agarwal và Lisa Choe (2014, 2018) đã làm sáng tỏ hiện tượng mất ổn định tổng thể, uốn dọc và mất ổn định cục bộ của cột thép chịu nén dưới trường nhiệt phi tuyến. Về mặt mô phỏng số, Schleich (1980s) và Phan (2000s) đã ứng dụng phương pháp phần tử hữu hạn thông qua các phần mềm SAFIR, ABAQUS, ANSYS và FDS (Fire Dynamics Simulator) để mô phỏng tương tác cháy - nhiệt - kết cấu.

Trong y văn học thuật, tồn tại cuộc tranh luận kéo dài giữa hai trường phái:

  • Trường phái thiết kế truyền thống (Prescriptive Design): Dựa vào việc gán bậc chịu lửa cố định và tra bảng chiều dày bọc bảo vệ theo các kích thước danh định từ thử nghiệm đốt lò tiêu chuẩn ISO 834. Trường phái này bị chỉ trích là thiếu linh hoạt, gây lãng phí vật liệu hoặc đánh giá sai nguy cơ sụp đổ khi sơ đồ kết cấu thực tế khác biệt với mẫu đốt lò.
  • Trường phái thiết kế theo hiệu năng (Performance-Based Structural Fire Engineering): Đòi hỏi phân tích chính xác trường nhiệt độ biến thiên theo thời gian và không gian, tính toán trực tiếp sự suy giảm các thông số cơ lý của vật liệu thép (mô đun đàn hồi $E_{a,\theta}$, giới hạn chảy $f_{y,\theta}$) và trạng thái ứng suất - biến dạng của hệ kết cấu.

Tại Việt Nam, các nghiên cứu của GS. Phạm Văn Hội (2000s) về khung nhà thép chịu lửa, TS. Chu Thị Bình (2010) về thiết kế an toàn cháy bằng phần mềm SAFIR, NCS. Nguyễn Đức Việt về cột thép bọc thạch cao chịu nén đúng tâm, NCS. Trương Quang Vinh (2015) về kết cấu liên hợp thép - bê tông, và NCS. Hoàng Anh Giang (2018) về sàn bê tông rỗng ứng lực trước đã đặt những viên gạch đầu tiên. Tuy nhiên, chưa có công trình nào tích hợp hoàn chỉnh quy trình từ mô hình truyền nhiệt 3D sang phương pháp tính toán dẻo đơn giản hóa SDM kết hợp kiểm chứng mô phỏng số 3D cho cả dầm và cột thép chữ I bọc vữa và thạch cao trong điều kiện xây dựng Việt Nam.

Luận án này định vị chính xác tại điểm giao thoa giữa lý thuyết Eurocode tiên tiến và thực tiễn kỹ thuật Việt Nam. Nghiên cứu thúc đẩy lĩnh vực tiến lên thông qua việc so sánh với hai nghiên cứu quốc tế điển hình: nghiên cứu của Franssen (1995) về hiện tượng oằn cột thép và nghiên cứu của Anil Agarwal & Lisa Choe (2014) về tác động của gradient nhiệt không đều, chứng minh rằng mô hình giải tích đề xuất đạt độ tương thích cao nhưng giảm thiểu đáng kể chi phí tính toán.

Đó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à phát triển Lý thuyết Kỹ thuật Kết cấu Chịu lửa (Structural Fire Engineering Theory) thông qua việc hệ thống hóa và cụ thể hóa các mô hình suy giảm đặc tính nhiệt - cơ của vật liệu thép kết cấu theo nhiệt độ $\theta$. Nghiên cứu áp dụng và kiểm chứng mô hình vật liệu phi tuyến theo tiêu chuẩn Eurocode EN 1993-1-2:2005:

  • Hệ số suy giảm mô đun đàn hồi: $$k_{E,\theta} = \frac{E_{a,\theta}}{E_a}$$
  • Hệ số suy giảm giới hạn chảy: $$k_{y,\theta} = \frac{f_{y,\theta}}{f_y}$$
  • Hệ số suy giảm giới hạn tỷ lệ: $$k_{p,\theta} = \frac{f_{p,\theta}}{f_p}$$
                Biểu đồ quan hệ Ứng suất - Biến dạng của Thép ở Nhiệt độ cao
    Ứng suất (σ)
         ^
  fy,θ --+------------------- Giai đoạn chảy dẻo (fp,θ < σ ≤ fy,θ)
         |                  / 
  fp,θ --+                 / 
         |                /  Giai đoạn đàn hồi tuyến tính (σ = Ea,θ * ε)
         |               /
         |              /
         +-------------+----------------------------------------> Biến dạng (ε)
         0           εp,θ   εy,θ       εt,θ                   εu,θ

Khung lý thuyết của luận án đưa ra 4 mệnh đề khoa học được chứng minh bằng thực nghiệm số:

  • Mệnh đề 1 (P1): Sự suy giảm khả năng chịu lực của dầm thép liên kết bản sàn bê tông trong đám cháy chịu sự chi phối chủ yếu bởi sự chênh lệch nhiệt độ giữa cánh trên và cánh dưới, làm dịch chuyển trục trung hòa dẻo $z_{pl,\theta}$ về phía bản cánh có nhiệt độ thấp hơn.
  • Mệnh đề 2 (P2): Sự mất ổn định tổng thể của cột thép chữ I chịu nén - uốn đồng thời trong đám cháy xảy ra do sự tương tác phi tuyến giữa độ mảnh $\lambda_{\theta,max}$, độ lệch tâm tải trọng $e$, và sự phân bố nhiệt độ đối xứng hoặc bất đối xứng trên 4 mặt cột.
  • Mệnh đề 3 (P3): Cơ chế truyền nhiệt qua lớp bảo vệ thạch cao là quá trình thu nhiệt phân rã phân tử nước kết tinh (21% khối lượng nước bị canxi hóa ở nhiệt độ dưới $1.200^\circ\text{C}$), tạo thành lớp đệm hơi nước làm chậm quá trình tăng nhiệt của lõi thép.
  • Mệnh đề 4 (P4): Trạng thái ứng xử dẻo của kết cấu trong điều kiện cháy được mô tả chính xác thông qua việc phân tích quan hệ phi tuyến đa đoạn (Multilinear Isotropic Hardening - MISO) giữa tensor ứng suất nhiệt ${\sigma}$ và tensor biến dạng tổng ${\varepsilon_{tt}} = {\varepsilon_{el}} + {\varepsilon_{pl}} + {\varepsilon_{th}} + {\varepsilon_{cr}}$.

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

Khung phân tích độc đáo của luận án tích hợp 3 lý thuyết nền tảng:

  1. Lý thuyết truyền nhiệt phi tuyến đa trường (Fourier Heat Conduction combined with Stefan-Boltzmann Radiation and Convection): Phân tích sự trao đổi nhiệt giữa luồng khí nóng của đám cháy và bề mặt cấu kiện.
  2. Lý thuyết cơ học kết cấu phi tuyến nhiệt - cơ: Xác định trạng thái nội lực tới hạn $M_{fi,t,Rd}, N_{fi,t,Rd}$ dựa trên trường nhiệt độ phân bố cục bộ.
  3. Lý thuyết tổ hợp tải trọng tai nạn cháy (EN 1991-1-2:2002): Xác định tải trọng tính toán tác dụng lên kết cấu trong điều kiện chịu lửa thông qua hệ số giảm tải $\eta_{fi}$:

$$q_{fi} = \gamma_{GA} G_k + \psi_{1,1} Q_{k,1} + \sum_{i \ge 2} \psi_{1,i} Q_{k,i} + A_d$$

$$\eta_{fi} = \frac{\gamma_{GA} G_k + \psi_{1,1} Q_{k,1}}{\gamma_G G_k + \gamma_{Q,1} Q_{k,1}}$$

Trong đó:

  • $\gamma_{GA} = 1,0$ là hệ số vượt tải của tải trọng thường xuyên trong điều kiện cháy.
  • $\gamma_G = 1,35$ và $\gamma_{Q,1} = 1,5$ là các hệ số vượt tải ở điều kiện làm việc bình thường.
  • $\psi_{1,1} \in [0,5; 0,9]$ là hệ số tổ hợp hoạt tải chính.
  • $A_d$ là tác động gián tiếp do biến dạng nhiệt gây ra trong hệ kết cấu siêu tĩnh.

Điều kiện biên và phạm vi áp dụng (Boundary Conditions) được xác định rõ ràng: cấu kiện không có giảm yếu tiết diện dọc chiều dài; điều kiện liên kết hai đầu dầm/cột được duy trì không đổi trong suốt quá trình chịu lửa; đường cong nhiệt độ đám cháy danh nghĩa tuân thủ tiêu chuẩn ISO 834-1:

$$T = 345 \log_{10}(8t + 1) + 20$$

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường nhận thức luận Thực chứng (Positivism) kết hợp Phân tích suy diễn toán học (Deductive Mathematical Modeling). Thiết kế nghiên cứu đa phương pháp (Mixed Methods / Multi-method Analytical Design) kết hợp chặt chẽ giữa 3 trục kỹ thuật:

  1. Phát triển thuật toán phần tử hữu hạn truyền nhiệt 3D (Thuật toán DT3D lập trình trên nền tảng MATLAB).
  2. Xây dựng quy trình tính toán kỹ thuật đơn giản hóa SDM (Simple Design Method).
  3. Mô phỏng số chi tiết đa trường vật lý tương hỗ (Thermal - Structural Coupled FEA) trên phần mềm chuyên dụng ANSYS Workbench.
                              QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU
                              
   +------------------------------------------------------------------------+
   |                       MÔ HÌNH ĐÁM CHÁY DANH NGHĨA                     |
   |              Đường cong nhiệt độ - thời gian tiêu chuẩn ISO 834        |
   +-----------------------------------+------------------------------------+
                                       |
                                       v
   +-----------------------------------+------------------------------------+
   |                          BÀI TOÁN TRUYỀN NHIỆT                         |
   |   - Phương trình Fourier 3D vi phân đạo hàm riêng                      |
   |   - Thuật toán DT3D (MATLAB) & Thermal Analysis (ANSYS Workbench)      |
   |   - Kết quả: Trường nhiệt độ T(x,y,z,t) trên toàn bộ cấu kiện         |
   +-----------------+------------------------------------+-----------------+
                     |                                    |
                     v                                    v
   +-----------------+------------------+   +-------------+------------------+
   |   PHƯƠNG PHÁP ĐƠN GIẢN HÓA (SDM)   |   |     MÔ PHỎNG SỐ NÂNG CAO (FEM)   |
   | - Phân chia tiết diện thành n phần |   | - ANSYS Structural Analysis    |
   | - Suy giảm cơ lý theo EN 1993-1-2  |   | - Mô hình vật liệu dẻo MISO    |
   | - Tính M_fi,t,Rd & N_fi,t,Rd       |   | - Phân tích Eigenvalue Buckling|
   +-----------------+------------------+   +-------------+------------------+
                     |                                    |
                     +-----------------+------------------+
                                       |
                                       v
   +-----------------------------------+------------------------------------+
   |            ĐỐI CHIẾU, ĐÁNH GIÁ VÀ KIỂM CHỨNG KẾT QUẢ                   |
   |   - So sánh với thí nghiệm thực tế (SP Boras, Thụy Điển)               |
   |   - Đánh giá sai số tương đối, xác lập phạm vi ứng dụng cho TCVN      |
   +------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình giải bài toán truyền nhiệt ba chiều dựa trên phương trình vi phân truyền dẫn nhiệt:

$$\frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} + \frac{\partial^2 T}{\partial z^2} = \frac{\rho C}{\lambda} \frac{\partial T}{\partial t}$$

Trong đó:

  • $T(x,y,z,t)$ là nhiệt độ tại điểm có tọa độ $(x,y,z)$ ở thời điểm $t$.
  • $\rho, C, \lambda$ lần lượt là khối lượng thể tích, nhiệt dung riêng và hệ số dẫn nhiệt của vật liệu (biến thiên phi tuyến theo nhiệt độ).

Thuật toán DT3D sử dụng phần tử hữu hạn 3 chiều 6 mặt (phần tử lục diện 8 nút), xây dựng ma trận truyền nhiệt phần tử $[k^e]$ và ma trận truyền nhiệt toàn hệ thống $[K]$:

$$[K] {T} + [C] \left{ \frac{\partial T}{\partial t} \right} = {Q}$$

Quá trình rời rạc hóa thời gian sử dụng sơ đồ sai phân hữu hạn ẩn Crank-Nicolson nhằm đảm bảo tính ổn định tuyệt đối của nghiệm số:

$$\left( \frac{1}{\Delta t} [C] + \theta_{diff} [K] \right) {T}{n+1} = \left( \frac{1}{\Delta t} [C] - (1 - \theta{diff}) [K] \right) {T}n + {Q}{n+\theta_{diff}}$$

Tính hợp thức và độ tin cậy của thuật toán DT3D và mô hình mô phỏng số ANSYS được kiểm chứng tam giác hóa (Triangulation) thông qua việc so sánh với các dữ liệu thực nghiệm đốt lò chuẩn quốc tế của SP Boras (Thụy Điển) trên dầm thép không bọc và cột thép chèn gạch. Các điểm đo nhiệt độ trên bản cánh ($F_1, F_3, F_4, F_6$) và bản bụng ($W_1, W_3, W_4, W_6$) đều cho thấy sự trùng khít xuất sắc giữa đường cong mô phỏng và thực nghiệm với hệ số tương quan $R^2 > 0,96$.

Data và phân tích

Phân tích số được thực hiện trên tập hợp các cấu kiện dầm và cột thép chữ I tiêu chuẩn với các thông số hình học và cơ lý xác thực:

  • Dầm thép: Chiều dài $L = 8\text{m}$, tiết diện I định hình, liên kết hai đầu ngàm hoặc khớp, chịu tải trọng phân bố đều $q_{fi}$, chịu lửa 3 mặt (cánh trên được bảo vệ bởi sàn bê tông), bọc vữa dày $t_{bv} = 15\text{mm}, 25\text{mm}$ hoặc bọc thạch cao dạng hộp dày $t_{tc} = 9\text{mm}, 13\text{mm}, 16\text{mm}$.
  • Cột thép: Chiều cao $H = 3,6\text{m}$, chịu lực nén đúng tâm ($e = 0$) và lệch tâm ($e = 200\text{mm}$), chịu lửa 4 mặt xung quanh, bọc vữa chu vi hoặc thạch cao dạng hộp.

Mô hình phân tích kết cấu (Structural Analysis) trên ANSYS Workbench tích hợp mô hình biến dạng dẻo đẳng hướng đa tuyến MISO (Multilinear Isotropic Hardening) phản ánh chính xác biểu đồ $\sigma - \varepsilon$ của thép ở các mức nhiệt độ từ $20^\circ\text{C}$ đến $1.000^\circ\text{C}$. Đối với phân tích cột, kỹ thuật phân tích mất ổn định Eigenvalue Buckling kết hợp phân tích phi tuyến hình học Large Deflection được sử dụng để bắt chính xác tải trọng tới hạn $N_{fi,t,Rd}$ và các dạng mốt mất ổn định uốn dọc, oằn vặn bên ($M_{b,fi,t,Rd}$).

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

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

  1. Quy luật gradient nhiệt độ không đồng đều trên tiết diện dầm: Khi dầm thép chữ I chịu tác động của đám cháy từ 3 phía dưới, cánh trên tiếp xúc với bản sàn bê tông có tốc độ tăng nhiệt chậm hơn đáng kể so với cánh dưới và bản bụng. Tại thời điểm $t = 60\text{ phút}$, chênh lệch nhiệt độ giữa cánh trên và cánh dưới có thể lên tới $200 - 320^\circ\text{C}$. Điều này dẫn đến sự suy giảm không đồng đều về cường độ dẻo và gây ra hiện tượng dịch chuyển vị trí trục trung hòa dẻo của tiết diện lên phía trên tới 15 - 22% chiều cao tiết diện dầm.
  2. Hiệu năng vượt trội của giải pháp bọc vữa chu vi và thạch cao dạng hộp: Đối với dầm $L=8\text{m}$ bọc vữa chống cháy dày $25\text{mm}$, sau 120 phút chịu lửa tiêu chuẩn ISO 834, nhiệt độ lớn nhất trong cốt thép chỉ đạt $412^\circ\text{C}$ (thấp hơn nhiều so với nhiệt độ tới hạn nguy hiểm $\theta_{cr} \approx 550 - 600^\circ\text{C}$), đảm bảo dầm hoàn toàn không bị phá hoại dẻo và giữ độ võng trong giới hạn cho phép. Tương tự, tấm thạch cao chống cháy ốp hộp dày $16\text{mm}$ duy trì nhiệt độ cột thép dưới $480^\circ\text{C}$ sau 90 phút.
  3. Mối quan hệ phi tuyến giữa độ lệch tâm tải trọng và khả năng chịu lửa của cột: Với cột thép chịu nén đúng tâm ($e = 0$), sự suy giảm khả năng chịu lực diễn ra chậm trong 30 phút đầu nhưng giảm đột ngột khi nhiệt độ vượt ngưỡng $400^\circ\text{C}$ do mất ổn định uốn dọc tổng thể ($\chi_{fi}$). Khi có độ lệch tâm ($e = 200\text{mm}$), mô men uốn phụ do biến dạng nhiệt gây ra làm giảm khả năng chịu lực tới hạn thêm 28 - 35% so với trường hợp nén đúng tâm ở cùng thời điểm khảo sát.
  4. Mức độ tương thích và sai số an toàn của phương pháp SDM so với ANSYS: Kết quả tính toán khả năng chịu lực theo quy trình đơn giản hóa SDM luôn cho giá trị nhỏ hơn (thiên về an toàn) so với kết quả mô phỏng số 3D chi tiết trên ANSYS Workbench từ 5,8% đến 12,4%. Sai số này nằm hoàn toàn trong phạm vi cho phép của kỹ thuật công trình, khẳng định phương pháp SDM hoàn toàn đủ độ tin cậy để đưa vào ứng dụng thực tế.
+----------------------------------------------------------------------------------------------------+
|               BẢNG SO SÁNH KHẢ NĂNG CHỊU LỰC DẦM VÀ CỘT THÉP THEO SDM VÀ MÔ PHỎNG ANSYS             |
+----------------------+--------------------+---------------+---------------+---------------+--------+
| Cấu kiện khảo sát    | Giải pháp bọc      | Thời gian (t) | Khả năng SDM  | Mô phỏng ANSYS| Sai số |
+----------------------+--------------------+---------------+---------------+---------------+--------+
| Dầm I ngàm L=8m      | Vữa chu vi 25mm    | 30 phút       | 285.4 kNm     | 302.1 kNm     | -5.5%  |
| Dầm I ngàm L=8m      | Vữa chu vi 25mm    | 60 phút       | 242.8 kNm     | 261.5 kNm     | -7.1%  |
| Dầm I ngàm L=8m      | Vữa chu vi 25mm    | 120 phút      | 168.2 kNm     | 188.0 kNm     | -10.5% |
| Cột I H=3.6m (e=0)   | Thạch cao hộp 16mm | 30 phút       | 1420 kN       | 1510 kN       | -5.9%  |
| Cột I H=3.6m (e=0)   | Thạch cao hộp 16mm | 60 phút       | 1180 kN       | 1290 kN       | -8.5%  |
| Cột I H=3.6m (e=0)   | Thạch cao hộp 16mm | 120 phút      | 650 kN        | 740 kN        | -12.1% |
| Cột I H=3.6m (e=200) | Thạch cao hộp 16mm | 60 phút       | 740 kN        | 815 kN        | -9.2%  |
+----------------------+--------------------+---------------+---------------+---------------+--------+

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở khoa học hoàn chỉnh để mở rộng lý thuyết tính toán trạng thái giới hạn chịu lực của kết cấu thép trong môi trường nhiệt độ cao phi tuyến, làm rõ cơ chế tương tác giữa gradient nhiệt và ứng xử dẻo của tiết diện chữ I.
  • Về mặt phương pháp luận: Thiết lập thành công chuỗi công cụ tính toán phân tầng: từ thuật toán chuyên sâu DT3D (nghiên cứu hàn lâm) đến mô hình mô phỏng ANSYS (kiểm tra dự án phức tạp) và quy trình đơn giản hóa SDM (thiết kế thực hành hàng ngày).
  • Về mặt thực tiễn kỹ sư: Giúp các kỹ sư kết cấu tại Việt Nam chủ động tính toán, lựa chọn chính xác chiều dày và hình thức bọc bảo vệ (vữa khoáng hoặc thạch cao) cho từng vị trí cấu kiện cụ thể, thay thế hoàn toàn việc phỏng đoán hay tra bảng thiếu căn cứ.
  • Về mặt chính sách và quy chuẩn: Đóng vai trò là tài liệu tham khảo cốt lõi để Bộ Xây dựng và Cục Cảnh sát PCCC & CNCH xây dựng Tiêu chuẩn Thiết kế Kết cấu Thép Chịu Lửa của Việt Nam, tương thích với Eurocode 3 và QCVN 06:2020/BXD.

Limitations và Future Research

Nghiên cứu thẳng thắn nhìn nhận 4 giới hạn khoa học cần được tiếp tục hoàn thiện:

  1. Giả định về sự bám dính toàn vẹn của lớp bọc bảo vệ: Mô hình tính toán giả định lớp vữa chống cháy và tấm thạch cao duy trì liên kết hoàn hảo với bề mặt thép trong suốt thời gian chịu lửa, chưa xét đến hiện tượng nứt vỡ cục bộ, co ngót nhiệt hoặc bong tróc cơ học (delamination) do gradient nhiệt và áp suất hơi nước tích tụ.
  2. Điều kiện biên liên kết lý tưởng: Các cấu kiện dầm và cột được khảo sát dưới dạng liên kết ngàm hoặc khớp lý tưởng không đổi. Trong thực tế đám cháy, độ cứng của nút liên kết (dầm - cột, bu lông, đường hàn) suy giảm phi tuyến theo nhiệt độ, gây ra sự suy giảm bậc siêu tĩnh của toàn khung.
  3. Mô hình đám cháy danh nghĩa đơn nhất: Nghiên cứu tập trung vào đường cong cháy chuẩn ISO 834 (đám cháy cellulose), chưa khảo sát đám cháy nhiên liệu hydrocarbon có tốc độ gia nhiệt cực nhanh hoặc mô hình đám cháy di động (traveling fire) trong không gian rộng lớn không có khoang cháy.
  4. Bỏ qua giai đoạn làm nguội (Cooling Phase): Chưa đánh giá ứng xử cơ học và biến dạng dẻo dư của kết cấu trong giai đoạn đám cháy suy tàn và làm nguội bằng nước chữa cháy - giai đoạn thường gây ra hiện tượng co ngắn đột ngột và đứt gãy liên kết bu lông.

Chương trình nghiên cứu tương lai (Future Research Agenda) mở ra 4 hướng tiếp cận:

  • Nghiên cứu cơ chế phá hủy liên kết nhiệt - cơ giữa lớp vật liệu bọc bảo vệ và bề mặt thép dưới tác động của gradient nhiệt lớn.
  • Mô phỏng ứng xử chịu lửa của các dạng nút liên kết dầm - cột thép (khung không gian 3D) có xét đến biến dạng nhiệt toàn cục.
  • Tích hợp mô hình cháy thực tế CFD bằng phần mềm FDS với phân tích kết cấu phi tuyến ANSYS (FDS-ANSYS Co-simulation).
  • Đánh giá khả năng chịu lực còn lại (Residual Load-Bearing Capacity) và phương án gia cường kết cấu thép sau khi xảy ra hỏa hoạn.

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

  • Tác động học thuật (Academic Impact): Luận án đóng góp nguồn dữ liệu giải tích và mô phỏng số có giá trị tham khảo cao, dự kiến thu hút nhiều trích dẫn khoa học trong các tạp chí chuyên ngành Kỹ thuật Xây dựng, Kết cấu và An toàn Cháy tại Việt Nam và khu vực.
  • Chuyển đổi công nghiệp xây dựng (Industry Transformation): Việc tối ưu hóa quy trình tính toán chiều dày bọc bảo vệ giúp các doanh nghiệp sản xuất nhà thép tiền chế, nhà thầu xây dựng công nghiệp (như tại Tổ hợp Hóa dầu Long Sơn, Nhà máy Samsung, LG Display) tiết kiệm từ 15% đến 25% chi phí vật liệu chống cháy trong khi vẫn bảo đảm tuyệt đối bậc chịu lửa theo thiết kế.
  • Ảnh hưởng chính sách nhà nước (Policy Influence): Cung cấp luận cứ khoa học thực chứng phục vụ công tác thẩm duyệt thiết kế PCCC của cơ quan quản lý nhà nước, tháo gỡ các vướng mắc kéo dài trong việc áp dụng Quy chuẩn QCVN 06 đối với các công trình công nghiệp và dân dụng kết cấu thép.
  • Lợi ích xã hội và môi trường (Societal & Environmental Benefits): Giảm thiểu nguy cơ sụp đổ công trình khi xảy ra hỏa hoạn, kéo dài thời gian thoát hiểm cho con người và tạo điều kiện an toàn cho lực lượng cứu hỏa tác chiến, giảm thiểu tối đa thiệt hại nhân mạng và của cải xã hội.

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

  • Nghiên cứu sinh và Giảng viên ngành Xây dựng: Tiếp cận thuật toán DT3D nguồn mở, phương pháp lập trình phần tử hữu hạn truyền nhiệt 3D trên MATLAB và kỹ thuật mô phỏng phi tuyến trên ANSYS Workbench.
  • Kỹ sư Kết cấu và Tư vấn Thiết kế: Sở hữu quy trình tính toán đơn giản hóa SDM rõ ràng từng bước, có thể lập bảng tính Excel để xác định nhanh khả năng chịu lửa của dầm, cột thép bọc bảo vệ.
  • Chủ đầu tư và Tổng thầu thi công: Lựa chọn được giải pháp bọc cách nhiệt tối ưu kinh tế - kỹ thuật (vữa chống cháy hay thạch cao chống cháy) phù hợp với công năng và kiến trúc công trình.
  • Cơ quan Thẩm duyệt và Ban hành Tiêu chuẩn: Có bộ số liệu đối chiếu chuẩn mực để xây dựng các tiêu chuẩn thiết kế và quy trình nghiệm thu an toàn cháy cho kết cấu thép tại Việt Nam.

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 thiết lập thành công mô hình giải tích dự báo sự dịch chuyển trục trung hòa dẻo của tiết diện dầm chữ I dưới trường nhiệt độ không đồng đều 3 chiều, mở rộng Lý thuyết Kỹ thuật Kết cấu Chịu lửa (Eurocode EN 1993-1-2:2005) trong điều kiện biên bọc cách nhiệt cục bộ tại Việt Nam.

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

Nghiên cứu đã xây dựng thuật toán DT3D độc lập trên MATLAB giải bài toán truyền nhiệt 3D phi tuyến và tích hợp trực tiếp dữ liệu trường nhiệt vào phương pháp đơn giản hóa SDM. Khi so sánh với các nghiên cứu của Franssen (1995, Bỉ) và Anil Agarwal (2014, Hoa Kỳ), phương pháp của luận án cho phép các kỹ sư thực hành tính toán nhanh khả năng chịu lực của cấu kiện bọc bảo vệ mà không đòi hỏi phần mềm thương mại đắt tiền, với độ lệch kiểm chứng dưới 12,4%.

3. Phát hiện bất ngờ nhất từ dữ liệu nghiên cứu là gì?

Phát hiện bất ngờ nhất là hiện tượng: đối với dầm thép bọc vữa 3 mặt liên kết sàn bê tông, dù nhiệt độ cánh dưới đạt tới $400^\circ\text{C}$ sau 60 phút, khả năng chịu mô men uốn của dầm vẫn duy trì trên 80% giá trị thiết kế ban đầu nhờ sự bảo vệ nhiệt của bản sàn phía trên và sự tham gia chịu nén dẻo của vùng thép cánh trên còn nguội.

4. Luận án có cung cấp đầy đủ quy trình và giao thức để tái lặp nghiên cứu (Replication Protocol) không?

Có. Luận án công bố chi tiết toàn bộ hệ phương trình vi phân truyền nhiệt, ma trận độ cứng phần tử, sơ đồ khối thuật toán DT3D, mã nguồn chương trình MATLAB (tại Phụ lục 1), cùng các thông số khai báo mô hình vật liệu dẻo MISO và chia lưới trên ANSYS Workbench, cho phép cộng đồng học thuật tái lập và kiểm chứng hoàn toàn kết quả.

5. Lộ trình phát triển nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Lộ trình 10 năm định hướng: (1) Nghiên cứu ứng xử nhiệt - cơ của liên kết nút khung thép; (2) Tích hợp mô hình cháy tương tác FDS - ANSYS cho các công trình nhịp lớn; (3) Ứng dụng trí tuệ nhân tạo (Machine Learning) để dự báo nhanh thời gian sụp đổ kết cấu chịu lửa; (4) Xây dựng bộ Tiêu chuẩn Thiết kế Kết cấu Thép Chịu Lửa Quốc gia của Việt Nam.

Kết luận

  1. Nghiên cứu đã hệ thống hóa và làm sáng tỏ cơ sở khoa học về quá trình truyền nhiệt và ứng xử cơ lý của kết cấu thép bọc bảo vệ trong điều kiện hỏa hoạn.
  2. Xây dựng thành công thuật toán và chương trình máy tính DT3D (MATLAB) giải bài toán truyền nhiệt ba chiều với độ chính xác cao được kiểm chứng qua thực nghiệm quốc tế.
  3. Đề xuất quy trình thiết kế đơn giản hóa SDM giúp xác định nhanh chóng và an toàn khả năng chịu lực của dầm và cột thép chữ I bọc vữa và thạch cao chống cháy.
  4. Thực hiện mô phỏng số phi tuyến chuyên sâu trên ANSYS Workbench, làm rõ cơ chế biến dạng, suy giảm độ cứng và mất ổn định tổng thể của cấu kiện dầm, cột.
  5. Xác lập phạm vi ứng dụng và chứng minh tính ưu việt về hiệu năng chịu lửa từ 30 đến 120 phút của các giải pháp bọc bảo vệ phổ biến tại thị trường Việt Nam.
  6. Đóng góp trực tiếp vào sự chuyển dịch từ phương pháp thiết kế danh định sang thiết kế hiệu năng, mở ra 3 nhánh nghiên cứu chuyên sâu về liên kết chịu lửa, mô phỏng cháy FDS-FEM tương hỗ và tiêu chuẩn hóa thiết kế an toàn cháy cho ngành xây dựng Việt Nam.