1
|
Wu B, Weng W, Zeng T, Xia Z, Su Z, Xie F. A Coupled Thermochemical Model for Predicting Fire-Induced Thermal Responses and Decomposition Behavior. Polymers (Basel) 2025; 17:939. [PMID: 40219328 PMCID: PMC11991143 DOI: 10.3390/polym17070939] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2025] [Revised: 03/23/2025] [Accepted: 03/28/2025] [Indexed: 04/14/2025] Open
Abstract
Composite materials are increasingly used in aerospace applications due to their high strength-to-weight ratio, but their fire safety remains a critical concern. This study develops a coupled thermochemical model to predict the thermal response and decomposition behavior of composite materials under high-temperature fire conditions. The framework integrates heat transfer, resin pyrolysis kinetics, and gas generation dynamics, employing the Rule of Mixtures to dynamically update temperature-dependent thermophysical properties (thermal conductivity, specific heat capacity, and density). Decomposition kinetics are governed by an n-th-order Arrhenius equation, explicitly resolving the gas convection effects on energy transport. The governing equations are solved numerically using a hybrid explicit/implicit finite element scheme, ensuring stability under severe thermal gradients. Experimental validation compliant with the 14 CFR Part 25 and ISO 2685 standards demonstrates high predictive accuracy. The model successfully captures key phenomena, including the char layer insulation effects, transient heat flux attenuation, and decomposition-induced property transition. This work establishes a computational foundation for optimizing fire-resistant composites in aerospace applications, addressing critical gaps in the existing models through coupled multiphysics representation.
Collapse
Affiliation(s)
- Bin Wu
- The Second Research Institute of CAAC, Chengdu 610041, China; (B.W.); (Z.X.); (Z.S.); (F.X.)
- Key Laboratory of Aviation Fuel Airworthiness and Green Development of Civil Aviation, Chengdu 610041, China
- School of Safety Science, Institute of Public Safety Research, Tsinghua University, Beijing 100084, China;
| | - Wenguo Weng
- School of Safety Science, Institute of Public Safety Research, Tsinghua University, Beijing 100084, China;
| | - Tai Zeng
- The Second Research Institute of CAAC, Chengdu 610041, China; (B.W.); (Z.X.); (Z.S.); (F.X.)
- Key Laboratory of Aviation Fuel Airworthiness and Green Development of Civil Aviation, Chengdu 610041, China
| | - Zuxi Xia
- The Second Research Institute of CAAC, Chengdu 610041, China; (B.W.); (Z.X.); (Z.S.); (F.X.)
- Key Laboratory of Aviation Fuel Airworthiness and Green Development of Civil Aviation, Chengdu 610041, China
| | - Zhengliang Su
- The Second Research Institute of CAAC, Chengdu 610041, China; (B.W.); (Z.X.); (Z.S.); (F.X.)
- Key Laboratory of Aviation Fuel Airworthiness and Green Development of Civil Aviation, Chengdu 610041, China
| | - Fei Xie
- The Second Research Institute of CAAC, Chengdu 610041, China; (B.W.); (Z.X.); (Z.S.); (F.X.)
- Key Laboratory of Aviation Fuel Airworthiness and Green Development of Civil Aviation, Chengdu 610041, China
| |
Collapse
|
2
|
A Fast Thermal 1D Model to Study Aerospace Material Response Behaviors in Uncontrolled Atmospheric Entries. MATERIALS 2022; 15:ma15041505. [PMID: 35208044 PMCID: PMC8875084 DOI: 10.3390/ma15041505] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/09/2021] [Revised: 01/24/2022] [Accepted: 02/10/2022] [Indexed: 12/04/2022]
Abstract
A preliminary thermal 1D numerical model for studying the demise behavior of stainless steel 316L, silicon carbide (SiC) and carbon fiber reinforced polymer (CFRP) during uncontrolled atmospheric entry is proposed. Test case modeling results are compared to experimental data obtained in the framework of ESA Clean Space initiative: material samples were exposed to different heat flux conditions using the Plasma Wind Tunnel (PWT) facilities at the Institute of Space Systems (IRS) of the University of Stuttgart. This numerical model approximates the heating history of the selected materials by simulating their thermal response and temperature profiles, which have trends similar to the experimental curves that are found. Moreover, when high heat flux conditions are considered, the model simulates the materials’ mass loss due to the ablation process: at the end of the simulation, the difference between the experimental and the modeled results is about 17% for CFRP and 35% for stainless steel. To reduce the model’s uncertainties, the following analysis suggests the need to consider the influence of adequate material thermophysical properties and the physical-chemical processes that affect the samples’ temperature profile and mass loss.
Collapse
|
3
|
Li H, Wang N, Han X, Yuan H, Xie J. Mechanism Identification and Kinetics Analysis of Thermal Degradation for Carbon Fiber/Epoxy Resin. Polymers (Basel) 2021; 13:polym13040569. [PMID: 33672880 PMCID: PMC7917751 DOI: 10.3390/polym13040569] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2021] [Revised: 02/07/2021] [Accepted: 02/09/2021] [Indexed: 11/16/2022] Open
Abstract
For carbon fiber epoxy resin used in aerostructure, thermal degradation mechanism and kinetics play an important role in the evaluation of thermal response and combustion characteristics. However, the thermal decomposition process and mechanism are difficult to unify strictly due to the complexity of the components from different suppliers. In the present study, a product of carbon fiber epoxy resin made by AVIC (Aviation Industry Corporation of China) composite corporation is examined to identify its thermal degradation mechanism and pyrolysis products by measurements, including simultaneous thermal analysis, Fourier transform infrared spectroscopy and mass spectrometry, establish the kinetic model by Kissinger/Friedman/Ozawa/Coats-Redfern methods. The results show thermal degradation occurs in three steps under the inert atmosphere, but in four steps under air atmosphere, respectively. The first two steps in both environments are almost the same, including drying, carbon dioxide escape and decomposition of the epoxy resin. In the third step of inert atmosphere, phenol is formed, methane decreases, carbon monoxide basically disappears and carbon dioxide production increases. However, in air, thermal oxidation of the carbonaceous residues and intermolecular carbonization are observed. Furthermore, thermal degradation reaction mechanism submits to the F4 model. These results provide fundamental and comprehensive support for the application of carbon fiber epoxy resin in aircraft industry.
Collapse
|
4
|
Li H, Wang N, Han X, Fan B, Feng Z, Guo S. Simulation of Thermal Behavior of Glass Fiber/Phenolic Composites Exposed to Heat Flux on One Side. MATERIALS 2020; 13:ma13020421. [PMID: 31963246 PMCID: PMC7014050 DOI: 10.3390/ma13020421] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/02/2019] [Revised: 01/12/2020] [Accepted: 01/14/2020] [Indexed: 11/17/2022]
Abstract
A 3D thermal response model is developed to evaluate the thermal behavior of glass fiber/phenolic composite exposed to heat flux on one side. The model is built upon heat transfer and energy conservation equations in which the heat transfer is in the form of anisotropic heat conduction, absorption by matrix decomposition, and diffusion of gas. Arrhenius equation is used to characterize the pyrolysis reaction of the materials. The diffusion equation for the decomposition gas is included for mass conservation. The temperature, density, decomposition degree, and rate are extracted to analyze the process of material decomposition, which is implemented by using the UMATHT (User subroutine to define a material’s thermal behavior) and USDFLD (User subroutine to redefine field variables) subroutines via ABAQUS code. By comparing the analysis results with experimental data, it is found that the model is valid to simulate the evolution of a glass fiber/phenolic composite exposed to heat flux from one side. The comparison also shows that longer time is taken to complete the pyrolysis reaction with increasing depth for materials from the numerical simulation, and the char region and the pyrolysis reaction region enlarge further with increasing time. Furthermore, the decomposition degree and temperature are correlated with depths, as well as the peak value of decomposition rate and the time to reach the peak value.
Collapse
Affiliation(s)
- Han Li
- School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China;
- College of Airworthiness, Civil Aviation University of China, Tianjin 300300, China; (N.W.); (X.H.); (B.F.); (Z.F.)
- Correspondence:
| | - Nasidan Wang
- College of Airworthiness, Civil Aviation University of China, Tianjin 300300, China; (N.W.); (X.H.); (B.F.); (Z.F.)
| | - Xuefei Han
- College of Airworthiness, Civil Aviation University of China, Tianjin 300300, China; (N.W.); (X.H.); (B.F.); (Z.F.)
| | - Baoxin Fan
- College of Airworthiness, Civil Aviation University of China, Tianjin 300300, China; (N.W.); (X.H.); (B.F.); (Z.F.)
| | - Zhenyu Feng
- College of Airworthiness, Civil Aviation University of China, Tianjin 300300, China; (N.W.); (X.H.); (B.F.); (Z.F.)
| | - Shijun Guo
- School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China;
- Centre of Aeronautics, Cranfield University, Bedfordshire MK43 0AL, UK
| |
Collapse
|
5
|
Toldy A, Szebényi G, Molnár K, Tóth LF, Magyar B, Hliva V, Czigány T, Szolnoki B. The Effect of Multilevel Carbon Reinforcements on the Fire Performance, Conductivity, and Mechanical Properties of Epoxy Composites. Polymers (Basel) 2019; 11:E303. [PMID: 30960287 PMCID: PMC6419153 DOI: 10.3390/polym11020303] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2018] [Revised: 02/01/2019] [Accepted: 02/06/2019] [Indexed: 11/20/2022] Open
Abstract
We studied the effect of a multilevel presence of carbon-based reinforcements-a combination of conventional load-bearing unidirectional carbon fiber (CF) with multiwalled carbon nanotubes (CNT) and conductive CNT-containing nonwoven carbon nanofabric (CNF(CNT))-on the fire performance, thermal conductivity, and mechanical properties of reference and flame-retarded epoxy resin (EP) composites. The inclusion of carbon fibers and flame retardant reduced the peak heat release rate (pHRR) of the epoxy resins. The extent to which the nanoreinforcements reduced the pHRR depended on their influence on thermal conductivity. Specifically, high thermal conductivity is advantageous at the early stages of degradation, but after ignition it may lead to more intensive degradation and a higher pHRR; especially in the reference samples without flame retardant. The lowest pHRR (130 kW/m²) and self-extinguishing V-0 UL-94 rating was achieved in the flame-retarded composite containing all three levels of carbon reinforcement (EP + CNF(CNT) + CNT + CF FR). The plasticizing effect of the liquid flame retardant impaired both the tensile and flexural properties; however, it significantly enhanced the impact resistance of the epoxy resin and its composites.
Collapse
Affiliation(s)
- Andrea Toldy
- Department of Polymer Engineering, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3-9., H-1111 Budapest, Hungary.
| | - Gábor Szebényi
- Department of Polymer Engineering, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3-9., H-1111 Budapest, Hungary.
| | - Kolos Molnár
- Department of Polymer Engineering, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3-9., H-1111 Budapest, Hungary.
- MTA⁻BME Research Group for Composite Science and Technology, Műegyetem rkp. 3., H-1111 Budapest, Hungary.
| | - Levente Ferenc Tóth
- Department of Polymer Engineering, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3-9., H-1111 Budapest, Hungary.
- Soete Laboratory, Department of Electrical Energy, Metals, Mechanical Constructions and Systems, Faculty of Engineering and Architecture, Ghent University, Technologiepark 903., B-9052 Zwijnaarde, Belgium.
| | - Balázs Magyar
- Department of Polymer Engineering, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3-9., H-1111 Budapest, Hungary.
| | - Viktor Hliva
- Department of Polymer Engineering, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3-9., H-1111 Budapest, Hungary.
| | - Tibor Czigány
- Department of Polymer Engineering, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3-9., H-1111 Budapest, Hungary.
- MTA⁻BME Research Group for Composite Science and Technology, Műegyetem rkp. 3., H-1111 Budapest, Hungary.
| | - Beáta Szolnoki
- Department of Organic Chemistry and Technology, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Budafoki út 8, H-1111 Budapest, Hungary.
| |
Collapse
|