Tổng quan về luận án
Bảo đảm an toàn sinh mạng và duy trì tính toàn vẹn của kết cấu công trình chịu tác động của hỏa hoạn là một trong những thách thức kỹ thuật cốt lõi của ngành kỹ thuật xây dựng hiện đại. Các thảm họa cháy 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, Hoa Kỳ, 2001) sau 56 phút và 86 phút khi nhiệt độ vùng nổ đạt xấp xỉ $1500^\circ\text{C}$ làm bung tróc hoàn toàn lớp bảo vệ chống cháy, vụ cháy kéo dài 17 giờ tại Tòa nhà Trụ sở 50 tầng ở Caracas (Venezuela, 2004), hay vụ hỏa hoạn tại Tòa nhà Windsor 32 tầng (Madrid, Tây Ban Nha) ở mức nhiệt $800^\circ\text{C}$ làm sập đổ 6 tầng sàn và phá hủy hệ thống cột-vách chịu lực chính, đã minh chứng cho tính chất hiểm họa khốc liệt của lửa đối với công trình xây dựng. Tại Việt Nam, thống kê trung bình hàng năm ghi nhận khoảng 2.200 vụ cháy với thiệt hại vật chất ước tính trên 745 tỷ đồng, điển hình là các vụ cháy nghiêm trọng tại Trung tâm Thương mại Quốc tế ITC TP. Hồ Chí Minh (2002) và Trung tâm Thương mại Hải Dương (2013).
Đối với kết cấu thép, mặc dù sở hữu tỷ lệ cường độ trên trọng lượng vượt trội, vật liệu thép lại có độ nhạy nhiệt rất cao: giới hạn chảy và mô đun đàn hồi suy giảm nghiêm trọng khi nhiệt độ vượt quá ngưỡng tới hạn $400^\circ\text{C} - 600^\circ\text{C}$. Nhằm nâng cao bậc chịu lửa, việc bọc bảo vệ kết cấu bằng các vật liệu cách nhiệt nhẹ như tấm thạch cao chống cháy và vữa chống cháy trở thành giải pháp kỹ thuật tối ưu, vừa giảm thiểu tải trọng tĩnh tác dụng lên móng so với bọc bê tông truyền thống, vừa tối ưu hóa chi phí so với giải pháp sơn chống cháy trương nở.
+------------------------------------------------------------------------------------+
| KỊCH BẢN CHÁY TỔNG THỂ |
| Mô hình cháy danh nghĩa ISO 834: T = 345*log10(8t + 1) + 20 |
+------------------------------------------+-----------------------------------------+
|
v
+------------------------------------------------------------------------------------+
| TRUYỀN NHIỆT KHÔNG GIAN BA CHIỀU |
| Phương trình Fourier 3D: d²T/dx² + d²T/dy² + d²T/dz² = (rho*C / lambda)*(dT/dt) |
| Thuật toán PTHH MATLAB DT3D & Mô phỏng ANSYS Workbench (Thermal Analysis) |
+------------------------------------------+-----------------------------------------+
|
v
+------------------------------------------------------------------------------------+
| ỨNG XỬ CƠ - NHIỆT VÀ KHẢ NĂNG CHỊU LỰC |
| Suy giảm đặc tính: fy,theta = ky,theta * fy | Ea,theta = kE,theta * Ea |
| * Phương pháp đơn giản hóa (SMD): Mfi,t,Rd (Dầm uốn) | Nfi,t,Rd (Cột nén) |
| * Phương pháp mô phỏng số (ANSYS): Mô hình MISO & Phân tích Eigenvalue Buckling |
+------------------------------------------------------------------------------------+
Khoảng trống nghiên cứu (Research Gap) lớn nhất tại Việt Nam hiện nay nằm ở sự thiếu hụt các quy chuẩn tính toán trực tiếp khả năng chịu lực của kết cấu thép bọc bảo vệ trong điều kiện nhiệt độ cao. Các tiêu chuẩn hiện hành như QCVN 06:2020/BXD hay TCVN 9311 (ISO 834) chỉ dừng lại ở việc quy định thời gian chịu lửa danh định và đánh giá thông qua phương pháp tra bảng catalog từ nhà sản xuất vật liệu cách nhiệt, thiếu vắng cơ sở khoa học để phân tích ứng xử phi tuyến cơ - nhiệt của cấu kiện thép chịu lực dưới các kịch bản cháy cụ thể. Luận án tiến sĩ này tập trung giải quyết khoảng trống học thuật và thực tiễn nêu trên thông qua các câu hỏi và giả thuyết nghiên cứu được xác lập tường minh:
- Câu hỏi nghiên cứu 1 (RQ1): Quy luật phân bố trường nhiệt độ ba chiều theo thời gian trong cấu kiện dầm và cột thép chữ I được bọc bảo vệ bằng tấm thạch cao (dạng hộp/chu vi) và vữa chống cháy (theo chu vi) diễn ra như thế nào dưới tác động của đường cong cháy tiêu chuẩn?
- Câu hỏi nghiên cứu 2 (RQ2): Khả năng chịu lực tới hạn về uốn ($M_{fi,t,Rd}$) của dầm và ổn định tổng thể ($N_{fi,t,Rd}$) của cột thay đổi ra sao theo các mốc thời gian chịu lửa ($t = 30, 60, 90, 120, 150$ phút) khi áp dụng phương pháp giải tích đơn giản hóa SMD (Simplified Method)?
- Câu hỏi nghiên cứu 3 (RQ3): Mức độ tương quan và độ tin cậy giữa phương pháp tính toán giải tích đơn giản hóa SMD và mô phỏng số phần tử hữu hạn phi tuyến nâng cao (ANSYS Workbench) là bao nhiêu, và phạm vi ứng dụng hợp lý cho từng phương pháp trong thực hành kỹ thuật xây dựng tại Việt Nam được xác định như thế nào?
Các giả thuyết nghiên cứu tương ứng:
- Giả thuyết 1 (H1): Thuật toán truyền nhiệt 3 chiều DT3D giải trên nền tảng phần tử hữu hạn có khả năng mô tả chính xác sự suy giảm gradient nhiệt độ qua lớp bọc cách nhiệt và trường nhiệt phi đối xứng trên tiết diện thép chữ I với độ sai lệch dưới $5%$ so với thực nghiệm quốc tế.
- Giả thuyết 2 (H2): Cơ chế canxi hóa giải phóng $21%$ khối lượng nước kết tinh của thạch cao chống cháy và hệ số dẫn nhiệt thấp ($\lambda = 0,15 - 0,25\text{ W/mK}$) của vữa chống cháy tạo ra hiệu ứng làm chậm quá trình tăng nhiệt độ lõi thép, giúp cấu kiện duy trì khả năng chịu lực vượt trên các tổ hợp tải trọng tai nạn cháy theo EN 1991-1-2:2002.
- Giả thuyết 3 (H3): Phương pháp đơn giản hóa SMD cung cấp kết quả đánh giá khả năng chịu lực thiên về an toàn ($8% - 14%$) so với mô phỏng số phi tuyến ANSYS, hoàn toàn đủ độ tin cậy để ứng dụng trong thiết kế kỹ thuật công trình.
Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa lý thuyết truyền nhiệt phi tuyến Fourier, mô hình ứng xử cơ - nhiệt suy giảm theo Eurocode 3 (EN 1993-1-2:2005) và nguyên lý tổ hợp tải trọng tác động bất thường EN 1991-1-2:2002. Phạm vi nghiên cứu bao quát các cấu kiện dầm thép chữ I (chịu uốn theo điều kiện bền dưới tác động cháy 3 mặt) và cột thép chữ I (chịu nén đúng tâm và lệch tâm theo điều kiện ổn định tổng thể dưới tác động cháy 4 mặt), sử dụng vật liệu thạch cao chống cháy dày 9mm, 13mm, 16mm và vữa chống cháy dày từ 14mm đến 50mm, hướng tới ứng dụng cho các công trình công nghiệp và dân dụng quy mô lớn tại Việt Nam (như Tổ hợp Samsung Bắc Ninh/Thái Nguyên, LG Display Hải Phòng, Tổ hợp Hóa dầu Long Sơn).
Literature Review và Positioning
Nghiên cứu về ứng xử của kết cấu xây dựng trong điều kiện hỏa hoạn đã trải qua hơn năm thập kỷ phát triển mạnh mẽ với các trường phái học thuật tiêu biểu trên thế giới:
Trường phái Bắc Âu được khởi xướng từ những năm 1970 bởi O. Pettersson và S. Magnusson tại Đại học Lund (Thụy Điển), đặt nền móng cho khái niệm bậc chịu lửa của công trình dựa trên năng lượng đám cháy và đường cong nhiệt độ - thời gian. Đến thập niên 1980, Đại học Lund cùng Phòng thí nghiệm Quốc gia SP (Boras, Thụy Điển) dưới sự dẫn dắt của Y. Anderberg đã thực hiện hơn 100 mẫu thí nghiệm quy mô lớn nhằm khảo sát khả năng chịu uốn và nén của dầm, cột thép có bọc và không bọc bảo vệ, chỉ ra sự phụ thuộc chặt chẽ giữa độ mảnh ($\lambda$), độ lệch tâm ban đầu ($e$) và sự phân bố nhiệt độ không đều đối với thời gian sụp đổ kết cấu.
Trường phái Anh và Tây Âu phát triển đột phá vào thập niên 1990 với dự án thử nghiệm đốt cháy toàn phần tòa nhà thép 8 tầng nguyên hình tại Viện Nghiên cứu Xây dựng (BRE, Cardington, Vương quốc Anh). Kết quả từ BRE Cardington đã chứng minh hiện tượng màng chịu kéo của sàn bê tông liên hợp và sự tái phân phối nội lực trong hệ khung không gian giúp kết cấu duy trì khả năng chịu lực vượt xa dự báo của các mô hình cấu kiện đơn lẻ. Tiêu chuẩn BS 5950 Part 8 (Anh, 1990) và sau đó là Bộ tiêu chuẩn Eurocode (EN 1991-1-2, EN 1993-1-2, EN 1994-1-2) được hoàn thiện bởi các nhà khoa học như Franssen (phát triển phần mềm SAFIR tại Đại học Liège, Bỉ) và Schleich (Luxembourg), chuẩn hóa phương pháp giải tích đơn giản hóa và mô hình số ngoài đàn hồi cho kết cấu thép và liên hợp thép - bê tông.
Trường phái Bắc Mỹ sau biến cố WTC 2001 đã tập trung nghiên cứu sâu tại Viện Tiêu chuẩn và Kỹ thuật Quốc gia Hoa Kỳ (NIST) với các công trình của Long T. McAllister, John L. Hurley, cùng các nghiên cứu thực nghiệm và mô phỏng của Anil Agarwal và Lisa Choe (2014) về cột thép chịu tải trọng lệch trục trong điều kiện cháy không đồng đều. Các kết quả này được tích hợp vào hệ tiêu chuẩn ASCE/SEI/SFPE 29-05.
+------------------------------------------------------------------------------------+
| LỊCH SỬ PHÁT TRIỂN & ĐỊNH VỊ LUẬN ÁN |
+------------------------------------------------------------------------------------+
| 1970s - Thụy Điển (Pettersson, Magnusson, Anderberg): Bậc chịu lửa, 100 mẫu thử SP |
| 1990s - Anh & Bỉ (BRE Cardington 8 tầng, Franssen/SAFIR): Khung 3D, Eurocode EC3 |
| 2000s - Hoa Kỳ (NIST WTC, ASCE 29-05, Agarwal & Choe): FDS CFD, Cháy bất đối xứng |
| Hiện nay - Việt Nam (Luận án NCS): Tích hợp MATLAB DT3D + SMD + ANSYS cho dầm/cột |
| thép chữ I bọc thạch cao và vữa chống cháy theo QCVN 06:2020 / EN 1993-1-2 |
+------------------------------------------------------------------------------------+
Trong bức tranh tổng quan học thuật, tồn tại hai luồng quan điểm và tranh luận kỹ thuật chính:
- Tranh luận về mô hình phân bố nhiệt: Luồng quan điểm truyền thống (áp dụng trong Eurocode đơn giản hóa) giả định nhiệt độ phân bố đều trên toàn bộ tiết diện ngang cấu kiện để đơn giản hóa quá trình tính toán. Ngược lại, các nhà nghiên cứu thực nghiệm hiện đại (như Agarwal & Choe, 2014) chỉ ra rằng trong cấu kiện thực tế—đặc biệt là dầm đỡ sàn bê tông hay cột tiếp xúc nhiệt không đối xứng—gradient nhiệt độ theo phương đứng và phương ngang gây ra sự uốn nhiệt (thermal bowing) và làm lệch trục trung hòa dẻo, đòi hỏi phải giải bài toán truyền nhiệt vi phân 2D/3D.
- Tranh luận về phương pháp tiếp cận thiết kế: Một trường phái ưu tiên phương pháp tra bảng và giải tích bán thực nghiệm (Prescriptive/Simplified approach) nhằm tối ưu hóa tính tiện dụng cho kỹ sư thiết kế; trong khi trường phái số hóa (Performance-based engineering) yêu cầu mô phỏng phần tử hữu hạn phi tuyến đa vật lý (Coupled Thermal-Stress FEA) để nắm bắt chính xác hiện tượng mất ổn định cục bộ và hiệu ứng dẻo hóa biến dạng lớn.
Tại Việt Nam, các nghiên cứu trước đây của GS.TS. Phạm Văn Hội (Đại học Xây dựng), TS. Chu Thị Bình (Đại học Kiến trúc Hà Nội), NCS. Nguyễn Đức Việt (Đại học PCCC, 2018), NCS. Trương Quang Vinh (2019) và NCS. Hoàng Anh Giang (IBST) chủ yếu tập trung vào kết cấu bê tông, panel sàn rỗng ứng lực trước, kết cấu liên hợp thép - bê tông hoặc khảo sát thực nghiệm đơn lẻ trên cột thép chèn gạch/bọc thạch cao chịu nén đúng tâm.
Vị trí học thuật của luận án được định vị là công trình tiên phong nghiên cứu toàn diện và có hệ thống cả ba mắt xích của bài toán chịu lửa: từ mô hình cháy, truyền nhiệt 3 chiều (thuật toán DT3D), đến quy trình giải tích đơn giản hóa SMD và mô phỏng số phi tuyến trên ANSYS Workbench cho cấu kiện dầm (chịu uốn 3 mặt lửa) và cột (chịu nén đúng tâm, lệch tâm 4 mặt lửa) được bảo vệ bằng vật liệu thạch cao và vữa chống cháy tại Việt Nam.
Đó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 cơ học kết cấu công trình chịu tác động nhiệt độ cao thông qua việc tích hợp sâu các quy luật suy giảm vật liệu theo Eurocode 3 (EN 1993-1-2:2005) vào điều kiện hình học và vật liệu bọc đặc thù:
+------------------------------------------------------------------------------------+
| HỆ THỐNG MỆNH ĐỀ LÝ THUYẾT VÀ QUY LUẬT CƠ - NHIỆT (P1-P4) |
+------------------------------------------------------------------------------------+
| [P1] Trục trung hòa dẻo dầm dịch chuyển lên trên do gradient nhiệt độ cánh-bụng |
| [P2] Khả năng chịu uốn Mfi,t,Rd suy giảm phi tuyến theo hệ số ky,theta và kp,theta|
| [P3] Tải trọng tới hạn cột Nfi,t,Rd chịu chi phối bởi hệ số uốn dọc phi tuyến |
| chi_fi = f(lambda_theta, ky,theta, kE,theta) |
| [P4] Tấm thạch cao duy trì bình nguyên nhiệt nhờ canxi hóa tách 21% H2O kết tinh |
+------------------------------------------------------------------------------------+
- Mệnh đề lý thuyết 1 (P1 - Dịch chuyển trục trung hòa dẻo): Dưới tác động nhiệt 3 mặt (đáy và hai bên bụng), gradient nhiệt độ giữa bản cánh dưới và bản cánh trên (được làm mát cục bộ bởi bản sàn bê tông) tạo ra sự suy giảm không đồng đều của giới hạn chảy $f_{y,\theta}$. Vị trí trục trung hòa dẻo của tiết diện chữ I bị dịch chuyển tịnh tiến về phía bản cánh trên chịu nén, làm biến đổi toàn diện phân bố ứng suất pháp trong phân tố.
- Mệnh đề lý thuyết 2 (P2 - Suy giảm khả năng chịu uốn phi tuyến): Mômen uốn cho phép chịu lửa $M_{fi,t,Rd}$ tại thời điểm $t$ không chỉ là hàm số của suy giảm cường độ vật liệu mà còn chịu chi phối bởi các hệ số điều chỉnh phân bố nhiệt không đều $k_1$ (trên tiết diện) và $k_2$ (dọc chiều dài dầm), thiết lập qua công thức:
$$M_{fi,t,Rd} = \frac{M_{fi,\theta,Rd}}{k_1 \cdot k_2}$$
- Mệnh đề lý thuyết 3 (P3 - Mất ổn định tổng thể phi đàn hồi theo nhiệt độ): Khả năng chịu nén tới hạn $N_{fi,t,Rd}$ của cột thép chữ I trong điều kiện chịu lửa là kết quả đồng thời của sự suy giảm mô đun đàn hồi ($k_{E,\theta} = E_{a,\theta} / E_a$) và giới hạn chảy ($k_{y,\theta} = f_{y,\theta} / f_y$), dẫn đến độ mảnh quy ước trong điều kiện cháy $\bar{\lambda}{\theta} = \bar{\lambda} \sqrt{k{y,\theta} / k_{E,\theta}}$ tăng nhanh phi tuyến, kích hoạt hiện tượng oằn vặn và mất ổn định Euler ở mức tải trọng thấp hơn điều kiện thường.
- Mệnh đề lý thuyết 4 (P4 - Động học phản ứng phân rã nhiệt thạch cao): Cơ chế làm việc của thạch cao chống cháy được mô hình hóa qua hai pha: pha phân hủy nhiệt giải phóng $21%$ khối lượng nước kết tinh ($\text{CaSO}_4 \cdot 2\text{H}_2\text{O} \rightarrow \text{CaSO}_4 \cdot 0.5\text{H}_2\text{O} + 1.5\text{H}_2\text{O}\uparrow$) tạo hiệu ứng thu nhiệt mạnh mẽ duy trì bề mặt lõi thép dưới ngưỡng $140^\circ\text{C}$ trong giai đoạn đầu, tiếp nối bởi pha cách nhiệt xốp của thạch cao khan với hệ số dẫn nhiệt cực thấp $\lambda = 0,20 - 0,25\text{ W/m}^\circ\text{C}$.
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:
- Lý thuyết truyền nhiệt vi phân 3 chiều Fourier:
$$\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}$$
- Lý thuyết dẻo chảy dẻo đa tuyến MISO (Multilinear Isotropic Hardening) phản ánh mối quan hệ ứng suất - biến dạng thực tế của thép ở nhiệt độ cao bao gồm các giai đoạn đàn hồi tuyến tính, chuyển tiếp elip phi tuyến và bình nguyên chảy dẻo.
- Lý thuyết độ tin cậy và tổ hợp tải trọng tai nạn cháy EN 1991-1-2: Thiết lập mức giảm tải danh định thông qua hệ số:
$$\eta_{fi} = \frac{\gamma_{GA} G_k + \psi_{1,1} Q_{k,1}}{\gamma_G G_k + \gamma_{Q,1} Q_{k,1}}$$
với $\gamma_{GA} = 1,0$, $\gamma_G = 1,35$, $\gamma_{Q,1} = 1,5$ và hệ số tổ hợp hoạt tải $\psi_{1,1} = 0,5 - 0,9$.
Điều kiện biên và phạm vi áp dụng (Boundary Conditions): Luận án xác định rõ cấu kiện dầm được phân tích theo sơ đồ uốn thuần túy, liên kết hai đầu dầm là khớp hoặc ngàm lý tưởng; cấu kiện cột được phân tích cho cả sơ đồ khung (chuyển vị ngang tự do) và sơ đồ giằng (chuyển vị ngang bị khống chế) với độ lệch tâm $e = 0\text{ mm}$ (nén đúng tâm) và $e = 200\text{ mm}$ (nén lệch tâm).
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (Positivism), kết hợp phương pháp luận tam giác hóa (Triangulation) đa chiều bao gồm: Giải tích toán học (Mathematical formulation) - Lập trình thuật toán phần tử hữu hạn (MATLAB DT3D) - Mô phỏng số phi tuyến đa trường vật lý (ANSYS Workbench FEA) - Kiểm chứng đối soát với các bộ dữ liệu thực nghiệm chuẩn quốc tế (SP Boras và NIST).
+------------------------------------------------------------------------------------+
| THIẾT KẾ NGHIÊN CỨU TAM GIÁC HÓA |
+------------------------------------------+-----------------------------------------+
| |
+-------------------> 1. THUẬT TOÁN PTHH MATLAB DT3D <-------------------------------+
| * Giải phương trình vi phân truyền nhiệt Fourier 3D |
| * Rời rạc hóa thời gian Backward Euler, ma trận [K] |
| | |
| v |
+-------------------> 2. PHƯƠNG PHÁP ĐƠN GIẢN HÓA SMD <------------------------------+
| * Ứng dụng EN 1993-1-2:2005 xác định Mfi,t,Rd và Nfi,t,Rd |
| * Tích hợp bảng tra hệ số suy giảm ky,theta, kE,theta |
| | |
| v |
+-------------------> 3. MÔ PHỎNG SỐ ANSYS WORKBENCH <-------------------------------+
| * Thermal Analysis (SOLID70/90) gán trường nhiệt độ |
| * Structural Analysis (SOLID185) + Mô hình MISO phi tuyến |
| * Eigenvalue Buckling xác định tải trọng và dạng mất ổn định |
| | |
| v |
+-------------------> 4. KIỂM CHỨNG THỰC NGHIỆM CHUẨN <------------------------------+
| * Thí nghiệm dầm thép không bọc SP Boras (Anderberg) |
| * Thí nghiệm cột thép chèn gạch / bọc thạch cao NIST |
+------------------------------------------------------------------------------------+
Thiết kế nghiên cứu được triển khai theo cấu trúc 3 cấp độ:
- Cấp độ vi mô (Vật liệu): Mô hình hóa đặc tính nhiệt vật lý biến thiên theo nhiệt độ của thép carbon ($\rho_a = 7850\text{ kg/m}^3$, nhiệt dung riêng $C_a(\theta)$, hệ số dẫn nhiệt $\lambda_a(\theta)$, hệ số giãn nở nhiệt $\alpha_a(\theta)$) và các thông số thực nghiệm của vật liệu bọc (vữa chống cháy $\lambda = 0,15 - 0,25\text{ W/mK}$; thạch cao $\lambda = 0,20 - 0,25\text{ W/m}^\circ\text{C}$).
- Cấp độ trung mô (Tiết diện): Phân tích trường nhiệt độ và gradient nhiệt độ trên mạng lưới phần tử 2D của tiết diện chữ I bọc bảo vệ theo chu vi và dạng hộp.
- Cấp độ vĩ mô (Cấu kiện 3D): Khảo sát trạng thái ứng suất, biến dạng dẻo, độ võng dầm, chuyển vị ngang/dọc và hiện tượng mất ổn định tổng thể của cột thép 3D chiều dài $L = 6\text{m} - 8\text{m}$.
Quy trình nghiên cứu rigorous
Quy trình nghiên cứu bao gồm 4 giai đoạn logic khép kín:
[Mô hình cháy ISO 834]
[Quy trình giải tích SMD] [Mô phỏng số ANSYS Workbench]
(Xác định Mfi,t,Rd; Nfi,t,Rd) (Thermal Analysis + Structural + Buckling)
[So sánh, kiểm chứng & đánh giá]
- Giai đoạn 1 - Xác lập điều kiện biên nhiệt: Sử dụng đường cong cháy tiêu chuẩn danh nghĩa ISO 834-1:
$$T = 345 \log_{10}(8t + 1) + 20$$
để mô tả nhiệt độ khí cháy $T_g$. Trao đổi nhiệt tại bề mặt ngoài lớp bọc được thiết lập qua tương tác bức xạ (hệ số phát xạ $\epsilon_s = 0,7 - 0,8$, định luật Stefan-Boltzmann với hằng số $\sigma = 5,67 \times 10^{-8}\text{ W/m}^2\text{K}^4$) và đối lưu nhiệt (hệ số truyền nhiệt đối lưu $h_c = 25\text{ W/m}^2\text{K}$).
- Giai đoạn 2 - Phát triển thuật toán DT3D: Lập trình giải thuật phần tử hữu hạn 3D trên môi trường MATLAB (Phụ lục code PL1 của luận án). Xây dựng hàm hình dạng $N_i$, ma trận truyền nhiệt phần tử $[k_e]$, ma trận truyền nhiệt toàn hệ thống $[K]$ và giải hệ phương trình biến thiên thời gian phi tuyến thông qua sơ đồ sai phân ẩn từng bước (Backward Euler Time Integration Scheme).
- Giai đoạn 3 - Thiết lập quy trình SMD: Xây dựng quy trình chuẩn hóa tính toán khả năng chịu lực dầm chịu uốn theo điều kiện bền và cột chịu nén theo điều kiện ổn định tổng thể dựa trên việc phân chia tiết diện thành $n$ phân tố diện tích $A_i$ với nhiệt độ $\theta_i$ tương ứng từ kết quả DT3D, áp dụng các hệ số suy giảm $k_{y,\theta}, k_{E,\theta}$ từ EN 1993-1-2.
- Giai đoạn 4 - Mô phỏng số ANSYS và Kiểm chứng thực nghiệm: Thực hiện phân tích hai bước (Sequential Thermal-Stress Coupling):
- Mô hình Thermal Analysis: Sử dụng phần tử nhiệt khối 8 nút SOLID70 hoặc 20 nút SOLID90 để tính toán trường nhiệt độ $T(x,y,z,t)$.
- Mô hình Structural Analysis: Chuyển đổi lưới phần tử sang SOLID185, gán trường nhiệt độ từ bước 1 làm tải trọng thể tích, khai báo mô hình vật liệu đàn dẻo MISO có xét đến biến dạng nhiệt ${\epsilon_{th}} = \alpha_a \Delta T$, chạy thuật toán phi tuyến Newton-Raphson.
- Kiểm chứng độ tin cậy: So sánh dữ liệu mô phỏng với các kết quả thực nghiệm dầm thép không bọc và cột thép chèn gạch từ phòng thí nghiệm SP Boras (Thụy Điển).
+------------------------------------------------------------------------------------+
| SƠ ĐỒ MA TRẬN VÀ PHƯƠNG TRÌNH PHẦN TỬ HỮU HẠN TRONG DT3D |
+------------------------------------------------------------------------------------+
| Phương trình cân bằng nhiệt vi phân hệ: |
| [C]{dT/dt} + [K]{T} = {Q} |
| Trong đó: |
| [C] = Tích phân thể tích (rho * Ca * [N]^T * [N]) dV |
| [K] = Tích phân thể tích ([B]^T * [D] * [B]) dV + Tích phân mặt (h * [N]^T * [N])|
| {Q} = Véc-tơ nguồn nhiệt bức xạ và đối lưu bề mặt |
| Sơ đồ tích phân thời gian ẩn: |
| ([C]/dt + [K]_{n+1}) {T}_{n+1} = ([C]/dt) {T}_n + {Q}_{n+1} |
+------------------------------------------------------------------------------------+
Data và phân tích
Toàn bộ dữ liệu tính toán được thực hiện trên các tiết diện thép chữ I định hình phổ biến (I200 đến I500, mác thép CT3, SS400, S275, S355). Ma trận phân tích bao gồm:
- Mẫu dầm: Nhịp $L = 8\text{m}$, liên kết hai đầu ngàm và hai đầu khớp, bọc vữa chống cháy dày $t_{bv} = 15\text{mm}, 25\text{mm}$ (3 mặt chu vi) và tấm thạch cao chống cháy dày $9\text{mm}, 13\text{mm}, 16\text{mm}$ (3 mặt dạng hộp).
- Mẫu cột: Chiều cao $H = 4\text{m} - 6\text{m}$, liên kết hai đầu khớp và một đầu ngàm - một đầu tự do, bọc vữa 4 mặt theo chu vi và bọc thạch cao 4 mặt dạng hộp, chịu tải trọng đúng tâm ($e = 0$) và lệch tâm ($e = 200\text{mm}$).
Phát hiện đột phá và implications
Những phát hiện then chốt
Kết quả nghiên cứu của luận án mang lại 4 phát hiện đột phá có giá trị khoa học và thực tiễn sâu sắc:
+------------------------------------------------------------------------------------+
| BẢNG TỔNG HỢP SO SÁNH HIỆU NĂNG GIỮA CÁC GIẢI PHÁP BỌC |
+----------------------+--------------------+--------------------+-------------------+
| Tiêu chí khảo sát | Không bọc bảo vệ | Bọc thạch cao hộp | Bọc vữa chu vi |
| | | (dày 16mm) | (dày 25mm) |
+----------------------+--------------------+--------------------+-------------------+
| Tốc độ tăng nhiệt | Rất nhanh | Chậm (có bình | Rất chậm, ổn định |
| độ lõi thép | (25-35°C/phút) | nguyên 100-140°C) | (4-6°C/phút) |
+----------------------+--------------------+--------------------+-------------------+
| Nhiệt độ lõi thép | > 750°C (Phá hoại | ~ 320°C - 380°C | ~ 280°C - 340°C |
| tại t = 60 phút | hoàn toàn) | (Thép giữ 85% fy) | (Thép giữ 90% fy) |
+----------------------+--------------------+--------------------+-------------------+
| Giới hạn chịu lửa | 15 - 20 phút | 60 - 90 phút | 120 - 150 phút |
| đạt được (R) | | | |
+----------------------+--------------------+--------------------+-------------------+
| Hiện tượng uốn nhiệt | Biến dạng uốn dẻo | Gradient nhiệt | Phân bố nhiệt |
| & Mất ổn định | cực lớn tức thời | cánh-bụng ~ 120°C | đồng đều quanh chu|
| | | | vi tiết diện |
+----------------------+--------------------+--------------------+-------------------+
- Hiện tượng gradient nhiệt độ phi đối xứng và sự dịch chuyển trục trung hòa dẻo của dầm chữ I: Khi dầm thép chịu tác động cháy 3 mặt (có sàn bê tông phía trên), nhiệt độ bản cánh dưới tăng nhanh hơn bản cánh trên từ $250^\circ\text{C}$ đến $380^\circ\text{C}$ tại thời điểm $t = 60\text{ phút}$. Sự chênh lệch nhiệt độ này làm suy giảm cục bộ cường độ bản cánh dưới ($k_{y,\theta}$ giảm xuống dưới $0,35$), buộc trục trung hòa dẻo dịch chuyển lên phía bản cánh trên một khoảng bằng $15% - 28%$ chiều cao tiết diện dầm. Điều này làm thay đổi cơ chế chịu uốn và đẩy nhanh biến dạng võng phi tuyến ở giữa nhịp dầm.
- Hiệu ứng "Bình nguyên nhiệt độ" (Thermal Plateau) của tấm thạch cao chống cháy: Dữ liệu tính toán từ DT3D và mô phỏng ANSYS đã chứng minh một hiện tượng vật lý then chốt: trong khoảng thời gian từ $t = 10\text{ phút}$ đến $t = 35\text{ phút}$, nhiệt độ bên trong lõi thép bọc thạch cao gần như không tăng và duy trì ổn định ở mức $100^\circ\text{C} - 140^\circ\text{C}$. Đây là hệ quả trực tiếp từ quá trình canxi hóa thu nhiệt phân tách $21%$ nước kết tinh. Chỉ sau khi lượng nước này bay hơi hoàn toàn, nhiệt độ mới tiếp tục tăng với tốc độ khoảng $4,5^\circ\text{C}/\text{phút}$.
- Tác động của hình thức bọc đến hiện tượng mất ổn định tổng thể của cột thép: Đối với cấu kiện cột, hình thức bọc thạch cao dạng hộp (box protection) tạo ra khoang không khí đối lưu bao quanh tiết diện, giúp cân bằng nhiệt độ giữa bản cánh và bản bụng. Trong khi đó, bọc vữa theo chu vi (contour protection) tuy bám dính trực tiếp nhưng có khả năng duy trì giới hạn chịu lực lâu hơn ở các bậc chịu lửa cao ($t = 120 - 150\text{ phút}$). Khi chiều dày lớp vữa đạt $25\text{mm}$, khả năng chịu lực nén tới hạn $N_{fi,t,Rd}$ tại $t = 120\text{ phút}$ vẫn đạt trên $45%$ so với khả năng chịu lực ở nhiệt độ thường, ngăn chặn hoàn toàn hiện tượng oằn vặn vỡ cục bộ.
- Tương quan định lượng và độ an toàn giữa SMD và ANSYS: Phân tích so sánh trên 16 trường hợp dầm và cột cho thấy phương pháp đơn giản hóa SMD luôn đưa ra kết quả khả năng chịu lực thấp hơn (thiên về an toàn hơn) từ $8,2%$ đến $14,6%$ so với mô phỏng ANSYS Workbench. Nguyên nhân là do SMD sử dụng nhiệt độ cực đại trên tiết diện để xác định hệ số suy giảm chung, bỏ qua sự tham gia chịu lực cục bộ của những vùng vật liệu có nhiệt độ thấp hơn. Sai số giữa mô phỏng ANSYS/DT3D với dữ liệu thí nghiệm kiểm chứng quốc tế chỉ dao động trong khoảng $3,8% - 6,5%$, khẳng định độ tin cậy tuyệt đối của các mô hình số.
Implications đa chiều
- Đóng góp cho lý thuyết kết cấu: Cung cấp mô hình phân tích cơ - nhiệt kết hợp cho cấu kiện thép có lớp bọc bảo vệ, làm sáng tỏ cơ chế tương tác giữa gradient nhiệt độ phi tuyến và trạng thái mất ổn định hình học ngoài đàn hồi.
- Đổi mới phương pháp luận: Thiết lập chuỗi công cụ tính toán đồng bộ từ thuật toán vi phân 3 chiều (MATLAB DT3D) đến mô phỏng số phi tuyến (ANSYS), có thể chuyển giao và áp dụng cho các dạng tiết diện thép phức tạp khác như thép ống nhồi bê tông (CFT), thép hình chữ H, hộp hay dầm liên hợp.
- Ứng dụng thực tiễn cho kỹ sư thiết kế: Bộ biểu đồ và bảng tra khả năng chịu lực theo thời gian ($M_{fi,t,Rd}$, $N_{fi,t,Rd}$) được xây dựng trong luận án cho phép các kỹ sư kết cấu tại Việt Nam tra cứu và lựa chọn chính xác chiều dày lớp bọc (thạch cao 9, 13, 16mm; vữa 15, 25, 50mm) tương ứng với bậc chịu lửa yêu cầu của công trình mà không còn phải phụ thuộc vào catalog cảm tính của nhà cung cấp.
- Khuyến nghị chính sách và quy chuẩn: Cung cấp luận cứ khoa học vững chắc phục vụ việc biên soạn, cập nhật các Tiêu chuẩn Quốc gia (TCVN) về thiết kế kết cấu thép chịu lửa, hỗ trợ đắc lực cho công tác thẩm định thiết kế phòng cháy chữa cháy của Cục Cảnh sát PCCC & CNCH (Bộ Công an) và Bộ Xây dựng.
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 một số giới hạn nghiên cứu (Limitations):
- Giả định về tính toàn vẹn liên kết biên: Mô hình tính toán giả thiết liên kết ở hai đầu cấu kiện (ngàm/khớp) là lý tưởng và không thay đổi trong suốt quá trình cháy, chưa xét đến biến dạng co ngót nhiệt của bu-lông cường độ cao và đường hàn tại nút khung thực tế.
- Giả định cơ học về vật liệu bọc: Nghiên cứu giả định lớp thạch cao và vữa chống cháy duy trì liên kết bám dính hoàn hảo với bề mặt thép, chưa mô hình hóa hiện tượng nứt vỡ cơ học, bong tróc (delamination/spalling) dưới tác động của áp suất hơi nước quá nhiệt cục bộ trong đám cháy thực tế.
- Mô hình kịch bản cháy danh nghĩa: Phân tích chủ yếu áp dụng đường cong cháy danh nghĩa ISO 834 (nhiệt độ tăng liên tục theo thời gian), chưa khảo sát giai đoạn suy tàn và làm nguội (cooling phase) của đám cháy thực tế vốn có thể gây ra hiện tượng phá hoại trễ (delayed collapse).
Chương trình nghị sự nghiên cứu tương lai (Future Research Agenda):
- Hướng 1: Mở rộng nghiên cứu ứng xử chịu lửa của các hệ liên kết nút dầm - cột thép (liên kết bu-lông mặt bích, liên kết hàn) có bọc bảo vệ cách nhiệt.
- Hướng 2: Mô hình hóa hiện tượng bong tróc phi tuyến của vật liệu bọc bằng phương pháp phần tử hữu hạn nâng cao kết hợp cơ học phá hủy liên tục (Damage Mechanics).
- Hướng 3: Nghiên cứu ứng xử của khung thép không gian chịu tác động của mô hình cháy tham biến (Parametric Fire) và đám cháy cục bộ (Traveling Fire) có xét chu kỳ tăng - giảm nhiệt độ.
- Hướng 4: Ứng dụng trí tuệ nhân tạo (Machine Learning/Deep Neural Networks) để tối ưu hóa đa mục tiêu: giảm thiểu chiều dày/chi phí lớp bọc bảo vệ nhưng tối đa hóa thời gian chịu lửa của toàn bộ công trình.
Tác động và ảnh hưởng
Luận án tạo ra tác động lan tỏa mạnh mẽ trên nhiều phương diện:
+------------------------------------------------------------------------------------+
| MA TRẬN TÁC ĐỘNG VÀ ẢNH HƯỞNG |
+----------------------+-------------------------------------------------------------+
| Lĩnh vực tác động | Chỉ số & Biểu hiện cụ thể |
+----------------------+-------------------------------------------------------------+
| Học thuật | Đặt nền tảng cho ngành Kỹ thuật An toàn Cháy Kết cấu tại |
| | Việt Nam; công bố trên các tạp chí và hội thảo uy tín |
+----------------------+-------------------------------------------------------------+
| Công nghiệp & R&D | Chuẩn hóa giải pháp thiết kế bọc thạch cao và vữa chống cháy|
| | cho các tổng thầu xây dựng (Coteccons, Hòa Bình, Vinata...) |
+----------------------+-------------------------------------------------------------+
| Kinh tế - Xã hội | Tiết kiệm 15% - 25% chi phí bọc bảo vệ nhờ tối ưu chiều dày;|
| | bảo vệ an toàn sinh mạng hàng triệu người trong nhà cao tầng|
+----------------------+-------------------------------------------------------------+
| Tiêu chuẩn hóa | Đóng góp trực tiếp cơ sở khoa học để hoàn thiện hệ thống |
| | QCVN 06/BXD và biên soạn TCVN Thiết kế Kết cấu Thép Chịu Lửa|
+----------------------+-------------------------------------------------------------+
Về mặt học thuật, công trình đóng góp trực tiếp vào kho tàng tài liệu giảng dạy sau đại học và là tài liệu tham khảo then chốt cho các nhóm nghiên cứu chuyên sâu về Kỹ thuật An toàn Cháy Kết cấu (Structural Fire Engineering) tại các trường đại học khối kỹ thuật xây dựng.
Về mặt công nghiệp, các kết quả nghiên cứu đã được kiểm chứng thực tế tại các đại dự án trọng điểm quốc gia như Tổ hợp Samsung Bắc Ninh, Nhà máy Samsung Thái Nguyên, LG Display Hải Phòng, Trung tâm Hội nghị Quốc gia, Bảo tàng Hà Nội, Tòa nhà Keangnam Hanoi Landmark Tower và Tổ hợp Hóa dầu Long Sơn. Việc định lượng chính xác chiều dày lớp thạch cao và vữa giúp các chủ đầu tư và nhà thầu tiết kiệm hàng chục tỷ đồng chi phí vật liệu mà vẫn đảm bảo tuyệt đối tiêu chuẩn an toàn phòng cháy chữa cháy.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận thuật toán nguồn mở DT3D, phương pháp luận nghiên cứu đa cấp độ và cơ sở dữ liệu mô phỏng số chi tiết để tiếp tục mở rộng các đề tài nghiên cứu chuyên sâu.
- Kỹ sư thiết kế kết cấu công trình: Nhận được quy trình tính toán đơn giản hóa SMD dễ áp dụng, biểu đồ tra cứu khả năng chịu lực ($M_{fi,t,Rd}$, $N_{fi,t,Rd}$) chuẩn hóa, giải phóng khỏi sự phụ thuộc thụ động vào catalog nhà sản xuất.
- Cơ quan Quản lý Nhà nước và Thẩm duyệt PCCC: Có cơ sở dữ liệu khoa học chính xác để đánh giá, thẩm định các giải pháp thiết kế chống cháy cho kết cấu thép theo QCVN 06:2020/BXD.
- Doanh nghiệp Sản xuất Vật liệu Chống cháy: Nắm bắt được cơ chế truyền nhiệt và các yêu cầu nhiệt - cơ học khắt khe để cải tiến công nghệ sản xuất tấm thạch cao cốt sợi, thạch cao bọt và các mác vữa chống cháy nhẹ vermiculite/perlite chất lượng cao.
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 thiết lập Mô hình suy giảm khả năng chịu lực phi đối xứng dưới gradient nhiệt độ 3 chiều cho cấu kiện thép chữ I bọc bảo vệ. Luận án đã mở rộng lý thuyết tính toán của Eurocode 3 (EN 1993-1-2:2005) từ giả định phân bố nhiệt đồng đều sang trường nhiệt phi tuyến thực tế, định lượng hóa chính xác khoảng dịch chuyển của trục trung hòa dẻo ($15% - 28%$) và thiết lập quy trình tính toán $M_{fi,t,Rd}$, $N_{fi,t,Rd}$ tích hợp thuật toán truyền nhiệt DT3D.
2. Sự đổi mới về phương pháp luận so với các nghiên cứu quốc tế trước đây?
So với các nghiên cứu giải tích một chiều cổ điển của Anderberg (1988) hay mô hình 2D trong phần mềm SAFIR (Franssen), luận án đã đổi mới phương pháp luận bằng cách:
- Xây dựng thuật toán phần tử hữu hạn 3D (DT3D trên MATLAB) giải phương trình Fourier không dừng có xét đầy đủ điều kiện biên bức xạ phi tuyến Stefan-Boltzmann và đối lưu;
- Tích hợp chuỗi phân tích ghép đôi Thermal - Structural - Eigenvalue Buckling trên ANSYS Workbench với mô hình vật liệu dẻo đa tuyến MISO;
- Đối soát tam giác giữa thuật toán giải tích - mô phỏng số - thực nghiệm chuẩn quốc tế.
3. Phát hiện nào gây bất ngờ nhất và có bằng chứng định lượng từ số liệu?
Phát hiện bất ngờ nhất là Hiện tượng bình nguyên nhiệt độ (Thermal Plateau) kéo dài 20 - 35 phút của tấm thạch cao chống cháy. Thay vì nhiệt độ thép tăng tuyến tính theo thời gian cháy, phản ứng phân hủy canxi hóa thu nhiệt của $21%$ nước kết tinh trong lõi thạch cao đã khóa chặt nhiệt độ bề mặt thép ở ngưỡng $100^\circ\text{C} - 140^\circ\text{C}$ trong hơn một phần ba thời lượng chịu lửa ban đầu, giúp thép giữ nguyên $100%$ cường độ chịu lực $f_y$ và mô đun đàn hồi $E_a$ cho đến phút thứ 40 của đám cháy tiêu chuẩn.
4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication Protocol) hoàn chỉnh không?
Hoàn toàn có. Luận án cung cấp tường minh toàn bộ giao thức tái lập bao gồm:
- Mã nguồn chương trình MATLAB của thuật toán DT3D trong Phụ lục PL1;
- Tham số hình học, tính chất vật lý của vật liệu thép và lớp bọc (hàm $C_a(\theta), \lambda_a(\theta), \alpha_a(\theta)$);
- Quy trình phân chia lưới phần tử (SOLID70, SOLID185), khai báo mô hình vật liệu phi tuyến MISO và các bước thiết lập tải trọng, điều kiện biên chi tiết trong ANSYS Workbench.
5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào?
Lộ trình 10 năm được định hình qua 3 giai đoạn:
- Giai đoạn 1 (1-3 năm): Mở rộng nghiên cứu khả năng chịu lửa của các nút khung liên kết bu-lông và hàn có bọc bảo vệ cách nhiệt.
- Giai đoạn 2 (4-6 năm): Xây dựng mô hình số song sinh (Digital Twin) mô phỏng ứng xử nhiệt - cơ học thời gian thực của toàn bộ hệ kết cấu công trình chịu kịch bản cháy lan thực tế (Traveling Fire).
- Giai đoạn 3 (7-10 năm): Hoàn thiện và ban hành Tiêu chuẩn Quốc gia (TCVN) toàn diện về Thiết kế Kết cấu Chịu Lửa dựa trên hiệu năng (Performance-based Structural Fire Design) tại Việt Nam.
Kết luận
- Luận án đã giải quyết trọn vẹn và khoa học bài toán đánh giá khả năng chịu lửa của cấu kiện dầm và cột thép tiết diện chữ I được bọc bảo vệ bằng tấm thạch cao và vữa chống cháy, đáp ứng yêu cầu cấp thiết về an toàn công trình tại Việt Nam.
- Xây dựng thành công thuật toán phần tử hữu hạn ba chiều DT3D trong môi trường MATLAB, cho phép mô phỏng chính xác trường nhiệt độ truyền từ đám cháy qua các dạng cấu tạo bọc bảo vệ (dạng hộp và theo chu vi) vào lõi thép với sai số dưới $5%$ so với thực nghiệm quốc tế.
- Đề xuất quy trình tính toán khả năng chịu lực theo phương pháp đơn giản hóa SMD chuẩn mực, đưa ra các công thức xác định mômen chịu uốn cho phép $M_{fi,t,Rd}$ và lực nén cho phép $N_{fi,t,Rd}$ tại các thời điểm chịu lửa tiêu chuẩn ($t = 30, 60, 90, 120, 150\text{ phút}$).
- Thực hiện mô phỏng số phi tuyến nâng cao trên ANSYS Workbench (Thermal Analysis, Structural Analysis, Eigenvalue Buckling), làm sáng tỏ trạng thái ứng suất - biến dạng dẻo, hiện tượng dịch chuyển trục trung hòa dẻo và cơ chế mất ổn định tổng thể của cấu kiện.
- Khẳng định tính chính xác và an toàn của phương pháp SMD (sai lệch an toàn từ $8% - 14%$ so với ANSYS), xác lập cơ sở khoa học để chuyển giao ứng dụng vào thực tiễn thiết kế kết cấu nhà thép tại Việt Nam.
- Đóng góp những luận cứ khoa học then chốt phục vụ công tác chuẩn hóa hệ thống Tiêu chuẩn, Quy chuẩn Quốc gia về an toàn cháy kết cấu, mở ra các hướng nghiên cứu hiện đại về an toàn cháy công trình xây dựng trong tương lai.
