Next Article in Journal
UHMWPE/CaSiO3 Nanocomposite: Mechanical and Tribological Properties
Previous Article in Journal
The Effect of a Polymer-Stabilized Latex Cobinder on the Optical and Strength Properties of Pigment Coating Layers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanism Identification and Kinetics Analysis of Thermal Degradation for Carbon Fiber/Epoxy Resin

College of Airworthiness, Civil Aviation University of China, Tianjin 300300, China
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(4), 569; https://doi.org/10.3390/polym13040569
Submission received: 10 January 2021 / Revised: 7 February 2021 / Accepted: 9 February 2021 / Published: 14 February 2021
(This article belongs to the Section Polymer Physics and Theory)

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.

1. Introduction

Attributed to high resistance to corrosion, high specific strength and stiffness of FRP (Fiber Reinforced Polymer) material, its utilization in infrastructures and manufacturing is more and more popular. For example, concrete slabs with GFRP(Glass Fiber Reinforced Polymer) bars and hallow composite reinforcing system [1] and a novel prefabricated FRP composite jacket with an easy-fit and self-locking mechanical joining system [2] have been developed. In addition, carbon fiber epoxy resin has been widely used in the design and manufacture of aircraft load-bearing and functional structures. However, thermal behavior of FRP is completely different from that of aluminum, because it can soften, decompose and burn to release fumes, especially in fire, leading to threaten the safety of occupants. As a result, both the Federal Aviation Administration (FAA) and the European Aviation Safety Agency (EASA) require us to consider fire protection, flame retardancy, combustion product toxicity, and other thermal issues of composite aircraft structures in airworthiness certification [3]. For this reason, the generic material combustion computational model, ThermaKin [4,5,6,7] and a toxicity assessment model [8,9,10] have been developed by FAA. In addition, Tranchard [11] developed a material thermal response model based on SEMCEF for fire environments and validated it with the T700/M21 composite material used in Airbus 350. Li [12,13] also developed a thermal response model for composites in the fire environment based on ABAQUS.
Thermal response and combustion characteristics of carbon fiber epoxy resin in the fire environment mainly depends on the thermal decomposition of resin matrix. For this reason, Quintiere [14] has conducted reaction kinetic parameter tests, thermophysical property tests and flame retardant performance tests on fiber epoxy resins conforming to the Boeing Material Specification. Tranchard [15,16] also carried out similar tests on T700/M21 to obtain the reaction kinetic parameters of the material and to establish the reaction kinetic model, and to explore decomposition mechanism and the combustion products. Reaction kinetic parameters, including activation energy, reaction order and the pre-exponential factor, are typically obtained by the thermal analysis. The combination of thermal analysis and infrared or mass spectrometry analysis can more accurately analyze the reaction process and mechanism. Considerable progress has been made in modeling of kinetic [17,18,19,20,21,22,23], including multi- or single-step kinetic models that can be utilized in thermal–mechanical model [24], lightning strikes [25], fire behavior [11], curing process design, fiber recycling [26,27], etc.
Moreover, the complexity and even uncertainty of components of composites lead to a large amount of kinetics and mechanism research. Ahamad and Alshehri [28] studied the thermal degradation behavior of UF/epoxy resin blends (UFE) using a thermogravimetric analysis (TGA), Fourier transform infrared (FTIR) spectra and mass spectrometry (MS) to identify decomposition products such as CO2, CO, H2O, HCN, HNCO and NH3, and established a reaction kinetic model using the Coats–Redfern method to clarify the main decomposition mechanisms. Li and Wu [29] studied the pyrolysis performance and flame-retardant mechanism of flame retardant EP with respect to the synthetic amino resin intumescent flame retardant (AIFR), showing that AIFR advances the pyrolysis of EP, reduces the activation energy of pyrolysis, increases the residual carbon and decreases the intensity of the absorption peak in the infrared spectrum. Liu et al. [30] investigated the curing reaction of bisphenol A formaldehyde phenolic epoxy resin (bisANER) with diaminodiphenyl ether (DDE) and the thermal degradation properties of its cured products, and calculated the activation energy of bis-ANER/DDE non-isothermal curing reaction by the Kissinger method and Ozawa–Flynn–Wall method, also determined the thermal degradation mechanism under different atmospheres. Huang et al. [31] used thermogravimetric mass spectrometry to investigate the pyrolysis reaction behavior of the amino-phenolic resin. The kinetic parameters of the pyrolysis reaction of amino-phenolic resin were obtained by the kinetic analysis of the experimental data using the Coats–Redfern integration method, and the decomposition mechanism and pyrolysis products were clarified. Thomas et al. [32] used the standard Flynn–Wall–Ozawa (FWO) method to deduce the Arrhenius parameters from thermograms of biosourced substrates, hence, the energy of activation (Ea) was obtained. Balart et al. [33] studied systematic kinetics of the thermal degradation of recycled ABS in inert atmosphere with using Flynn–Wall–Ozawa (FWO), the Kissinger–Akahira–Sunose (KAS) and the Starink methods, and suggested the degradation process of ABS takes place in a single-step process. Franco-Urquiza et al. [34] performed the thermogravimetric analysis at 5, 10 and 20 °C/min allowed determining the degradation kinetic parameters based on the Friedman and Kissinger models for neat polyester and nanocomposites, which showed that heating rates promoted an increase in the temperature degradation. Krawiec et al. [35] used a testing methodology to determine the toxicometric indicators on belts made of classically used fabric–rubber composite material, and evaluated the degradation kinetics of the polymeric belts by (TGA). Zhang et al. [36] investigated the effect of basalt fiber content in HDPE matrix composites on the thermal decomposition process using dynamic thermogravimetric analysis and utilized improved Coats–Redfern (C-R), Flynn–Wall–Ozawa (F-W-O), Friedman and Kissinger methods to ascertain the specific apparent activation energy (Ea) of each component and composite material.
Therefore, for whether a new material or material modification, it is necessary to carry out thermal analysis, identify the thermal behavior or mechanism of degradation and establish the kinetic model. For carbon fiber epoxy resin used in aerostructure, it usually consists of a fiber, epoxy monomer, curing agent, catalyst and flame retardant, with complex composition and content. Therefore, the thermal decomposition process and mechanism of products from different suppliers are difficult to unify strictly. In this paper, for a new product of carbon fiber reinforced epoxy resin composite, BA3202, developed by AVIC (Aviation Industry Corporation of China) composite corporation for aerostructures, simultaneous thermal analysis (STA), FTIR spectra and mass spectrometry (MS) were conducted to obtain the reaction kinetic parameters, and a degradation kinetic model was established based on the Kissinger/Friedman/Ozawa/Coats–Redfern methods. Moreover, the thermal decomposition products of the material were identified to clarify the reaction mechanism. Compared to the published study in literature, a more comprehensive and thorough insight into thermal behavior of CFRP used in aerostructures is presented, and detailed data can support thermal response modeling, combustion modeling and toxicity assessment of material in fire.

2. Materials and Methods

Samples of carbon fiber reinforced epoxy resin matrix composite material used for aerostructures were provided by AVIC composite corporation (Beijing, China), in which matrix is mainly DGEBA. The preparation process was prepreg cutting-lay-up-autoclave molding–machining. The size of sample was 4 mm × 4 mm × 1.2 mm with the lamina thickness of prepreg of 0.15 mm, shown as Figure 1. Stacking sequence is [0/45/90/135]s with the fiber volume of 54% and the density of 1.55 g/cm3.
Thermal analysis experiments were performed using the simultaneous thermal analyzer (NETZSCH STA 449 F3, NETZSCH, Selb, Germany), FTIR spectrometer (BRUKER VERTEX 80, BRUKER, Billerica, MA, USA) and mass spectrometry (NETZSCH QMS403C, NETZSCH, Selb, Germany) to investigate the mass changes (TG, Thermogravimetric), derivative of mass changes (DTG, Derivative Thermogravimetric) and gas composition (FTIR/MS), referring to the test standards [37,38,39,40]. The sample used in the simultaneous thermal analyzer was heated from room temperature to 1100 °C at the heating rates of 5, 10, 20 and 40 °C/min, under the atmosphere of air and inert gas, respectively, with the purge flow rate set to 20 mL/min. The primary aim of the TG analyses was to know the thermal degradation characteristics of materials, and to obtain adequate data to perform kinetic analysis in order to obtain the Arrhenius parameters. IR spectra were recorded in the spectral range of 4500–650 cm−1 with a resolution of 4 cm−1 and an average number of scans of 16. The intensities of 17 ions (m/z = 14, 16, 17, 18, 26, 28, 30, 42, 44, 46, 56, 60, 64 and 94) were determined according to the previous work [16,41].

3. Results and Discussion

3.1. Thermal Degradation

Figure 2 shows the TG and DTG curves of carbon fiber reinforced epoxy resin matrix composite heated from room temperature to 1100 °C at 20 °C/min under air and inert atmosphere, respectively. Under inert atmosphere, the TG curve shows three mass loss stages, corresponding to three DTG curve peaks. The first stage of weight loss ranged from room temperature to 286 °C with the mass loss of 0.69%, and the maximum weight loss rate was 0.17%/min with the peak temperature of 244 °C. At this stage, physical escape of low-molecular-weight gases occurs, with the drying effect caused by the temperature rise from room temperature, where water was evaporated rapidly, justified by mass spectrometry in Section 3.3.2. The same drying process is also observed in thermos-gravimetric analysis of epoxy resins by Ahamad [28] and Shao [42].
The second stage ranged from 286 to 515 °C with the mass loss of 23.29%, where the maximum weight loss rate was 4.99%/min at the peak temperature of 410 °C. This range is close to the pyrolysis interval of epoxy resin, where the ether bond is broken to generate BPA (Bisphenol-A), and the phenol comes from the further decomposition of BPA and/or the breakage of the C-C bond connecting the benzene ring after the breakage of unilateral ether bond on the epoxy resin [43], which is an obvious degradation behavior of the resin matrix. The weight loss of 4.54% in the third stage occurred in a wide temperature range from 515 to 1100 °C, which was relatively gentle, and the maximum loss rate was 0.64%/min corresponding to the peak temperature of 537.4 °C. According to the literature [44], the pyrolysis temperature of BPA epoxy resin is higher than 660 °C, so this stage was considered as further pyrolysis of the polymer [28].
Under the air atmosphere, the specimen was almost completely degraded with the mass loss of 90.64% reached at 1100 °C. The existence of four peaks in the DTG curve indicated that the weight loss process was divided into four stages, and the peak temperature was shown in Figure 1. The characteristics of the first weight loss stage in air atmosphere was almost consistent with that in inert gas, which was a typical drying stage. The peak temperature in the second stage was only slightly lower, but the pyrolysis interval was also similar to that under the inert atmosphere, indicating that the degradation of the resin matrix also occurred at this stage. Therefore, it showed that the reaction atmosphere has less influence on the degradation reaction of the resin matrix. However, the percentage of weight loss and weight loss rate in the inert atmosphere were greater.
The temperature range of the third stage was from 497 to 710 °C with the mass loss of 18.38%, the maximum weight loss rate of 2.53%/min and the peak temperature of 684 °C. Thermal oxidation of intermediate transient carbonaceous residues [16], which is a process of intermolecular cross-linking, carbonization, and further removal of small molecules, occurred at this stage [29]. The weight loss in the fourth stage ranged from 710 to 1100 °C with the mass loss of 52.65%; the maximum weight loss rate was 2.55%/min, and the peak temperature was 745.5 °C. This was due to the reaction of the carbon fiber reinforced material with oxygen and the nearly complete oxidative decomposition of the carbon fiber caused by prolonged heating, which was supported by literature [45].

3.2. Kinetic Model

Figure 3 shows the TG and DTG curves at different heating rates under inert and air atmosphere. Two inflection points appeared in the TG curves for different heating rates under inert atmosphere, which corresponded to the two peaks of the DTG curve. As the temperature rise rate increases, the weight loss curve shifts to the higher temperature region [46]. However, in air atmosphere, the effect of the temperature rise rate on weight-loss behavior appears more complex due to the oxidative decomposition of carbon fiber. The lower the temperature rise rate is, the longer the time is to reach the same temperature, and so the more sufficient the reaction is. As a result, the two lower heating rates (5 and 10 °C/min) make the temperature at which the weight loss rate reached the peak value lower, so that two more peaks could be seen below 1000 °C. However, under the heating rates of 20 and 40 °C/min, the degradation was not completed within 1100 °C, and no corresponding peaks were observed.
The different methods for solving the kinetic constants of reactions are based on the same basic theoretical formula derived from the Arrhenius formula, as shown in the following equation:
d α d t = A exp ( E A R T ) f ( α ) ,
where α = ( M i M ) / ( M i M e ) is the conversion degree, M i , M , M e are the initial mass, the mass of the specimen at the current moment, and the final mass at the end of the reaction, respectively. f ( α ) is a function of the kinetic mechanism, which is a function of the conversion degree, depending on the specific degradation mechanism. R is the ideal gas constant, with a value of 8.314 J/(mol K).
In this paper, the kinetic parameters of the reaction were solved by the Kissinger [47] differential method, the Friedman [20] differential method and the Ozawa [22] integral method, respectively. These three methods did not require the degradation reaction mechanism of the material to obtain the activation energy. As to calculate the reaction order, it can be assumed that the mechanism function can be expressed as follows: f ( α ) = ( 1 α ) n . The Kissinger method [47] assumes the maximum reaction rate at the peak temperature. The Friedman method [20] assumes that the reaction kinetic constant A, EA and n are independent of the heating rate. While the Ozawa method [22] assumes that the conversion degree at the peak temperature at different heating rates is a constant (iso-conversion method). Additionally, the Coats–Redfern method was used to determine the reaction mechanism function. The data obtained under inert atmosphere are shown in the following Table 1.

3.2.1. Kissinger Method

Differentiation on both sides of (1) yields the following equation:
E A β R T m 2 = A n ( 1 α m ) n 1 exp ( E A R T m ) ,
where T m is the peak temperature from the DTG curve, α m , is the peak conversion corresponding to T m , which is obtained from the TG curve, and β = d T d t is the constant heating rate during the experiment. Kissinger considers n ( 1 α m ) n 1 is independent of β and its value is approximately equal to 1.
The activation energy can be obtained from the above equation by logarithm and derivation:
d ( ln ( β / T m 2 ) ) d ( 1 / T m ) = E A R ,
The thermal weight loss process of carbon fiber reinforced epoxy resin composites under inert atmosphere has two main stages. As shown in Figure 4, the intercept and slope data are obtained.
The relationship between the pre-exponential factor A, and the slope and intercept of the curve is as follows:
ln ( A / ( k ) ) = a ,
When calculating the reaction order n, consider the mechanism function to be expressed as f ( α ) = ( 1 α ) n and n is calculated by the following equation:
n ( 1 α m ) n 1 = 1 + ( n 1 ) 2 R T m E A ,
The results are shown in the Table 2 below:

3.2.2. Friedman Method

Taking the logarithm of (1) gives the following equation
ln ( β d α d T ) = ln ( A ) + n ln ( 1 α ) ( E A R T ) ,
where ln ( β d α d T ) = ln ( d α d t ) .
The relationship between ln ( d α d t ) and 1/T at different conversion degrees was calculated as shown in Figure 5.
Figure 6 shows a line chart of the activation energy with conversion degree, which illustrates that the activation energy is not a constant, but changes as the thermal decomposition reaction proceeds and the conversion degree increases, showing an overall upward trend. Therefore, the thermal decomposition reaction of carbon fiber reinforced epoxy composite under inert gas is a multi-step. Furthermore, the predominant factor was found to be 1.34287 × 1014.

3.2.3. Ozawa Method

The activation energy EA was calculated by the following formula:
E A = R 0.4567 d ( log β ) d ( 1 / T ) ,
The relationship between log β and 1/T was calculated for different conversion degrees as shown in the following Figure 7:
The pre-exponential factor was calculated by the following formula [48,49]:
log A = log β + log E A + 0.434 E A / R T log R 2 log T ,
At the same conversion degree, yielding the values of the corresponding pre-exponential factor for different heating rates, as shown in Figure 8. And the mean values were finally obtained. The activation energies obtained by Ozawa’s method were similar to those obtained by Friedman’s method, which verified the accuracy of the results.

3.2.4. Coats–Redfern Method

Coats and Redfern assumed the function of the kinetic mechanism is f ( α ) = ( 1 α ) n , and derived the following expression:
ln ( α T 2 ) = ln ( A R β E A ) [ 1 ( 2 R T E A ) ] ( E A R T ) ,
where α in the logarithmic term of the left equation is the result of integrating 1 / f ( α ) to α and then ignoring higher-order terms under the condition of low conversion rate after polynomial expansion. In order to study the kinetic model of thermal decomposition, assuming that the mechanism function is unknown, the alternative term to α in the logarithmic term is the integral form F ( α ) of 1 / f ( α ) , and considering that 2 R T E A approaches zero, the formula is as follows:
ln ( F ( α ) T 2 ) = ln ( A R β E A ) ( E A R T ) ,
Different mechanism models were established by this method, as shown in Table 3.
The ln ( F ( α ) T 2 ) − 1/T curves for the conversion rates from 0.2 to 0.7 at different heating rates were made to obtain the activation energy values. Pearson’s correlation coefficients were calculated to measure the correlation between the variables, as shown in Table 4. Comparing the activation energy values in Table 4 with those calculated by the three kinetic methods in Table 1. The average EA of F4 is the closest to the previous calculation, which was 211.691 kJ/mol. It is concluded that the pyrolysis reaction of carbon fiber reinforced epoxy resin composite in the inert atmosphere is the closest to that predicted by the F4 model. The expression of the mechanism function is: ( 1 α ) 4 . After averaging the activation energies and pre-exponential factors of the four heating rates, the reaction mechanism model was obtained as follows
d α d t = 8.26 × 10 11 exp ( 211691 R T ) ( 1 α ) 4 ,

3.3. Gas Phase Analysis

3.3.1. FTIR Analysis

The Gram–Schmidt spectrum under inert atmosphere in Figure 2 represented the FTIR absorption intensity of the escaping gas during pyrolysis of the material, showing that there were only two peaks at 425 °C and 557 °C, respectively, which were delayed compared to the corresponding DTG curve. This was due to the time delay between gas detection and gas generation in the FTIR spectrometer. The larger peak at 425 °C indicated that the amount of gas escaping at this stage was higher, and had a higher infrared extinction coefficient.
Figure 9 shows the infrared spectral data of the gas at the characteristic temperatures determined by the DTG curve. At the characteristic temperature of 244 °C, which was the first stage of mass loss, the products were carbon dioxide (CO2: 2300–2400 cm−1) and water (H2O: 3500–4000 cm−1). The same phenomenon is observed by Ahamad [28] in the study of urea-formaldehyde/epoxy resin blends. Thus, this stage was not only a drying stage, but also a physical escape of small molecule gases. Obvious degradation behavior of the resin matrix occurred in the second stage, corresponding to the infrared spectral data at 410 °C. Since the stretching vibration of saturated hydrocarbon C-H is usually close to the frequency absorption of 3000 cm−1, it is speculated that there was methane gas (CH4) according to the three peaks in the range of 2900–3170 cm−1: 3045 cm−1, 3016 cm−1 and 2970 cm−1. In addition to water and carbon dioxide, carbon monoxide (CO: 2000–2200 cm−1 with double peaks at 2048 cm−1 and 2070 cm−1) and phenol (main peaks of C6H5OH: 3650, 1608, 1510, 1260 and 1176 cm−1) were also observed. The above products were considered to be the chain scission and depolymerization reactions of BPA epoxy resins [31], where CO and CO2 are mainly derived from the ether group (R-O-R’), methyl group (-CH3), methylene group (-CH2) and other carbon-containing functional groups. At 537 °C, carbon monoxide almost disappeared and turned into carbon dioxide, while the production of the other gas products decreased compared with that at 410 °C.

3.3.2. MS Analysis

The following gas molecules were detected by mass spectrometry at a heating rate of 20 °C/min under inert atmosphere: CH2 (m/z = 14), CH4 (m/z = 16), H2O (m/z = 18), C2H2 (m/z = 26), CO (m/z = 28), C2H4 (m/z = 28), C3H6 (m/z = 42), CO2 (m/z = 44), C3H8 (m/z = 44) and phenol C6H5OH (m/z = 94).
According to the order of magnitude of ion current intensity, the classification analysis was carried out from high to low. As shown in Figure 10, the left ordinate shows the mass spectra of gas products with a mass-to-charge ratio of 28, which was supposed to be CO and/or C2H4. Additionally, there were three peaks in the curve, which respectively occurred at 125 °C, 437 °C and 920 °C. The peak with the largest ion current intensity appeared between 300 and 500 °C, corresponding to the main decomposition peak in the DTG analysis, which was presumed to be the precipitation of CO gas, also observed by Tranchard [16] and Ahamad [28] and related to the certain amount of oxygen-containing groups in the polymer. The mass spectrum of phenol with the mass-to-charge ratio of 94 was shown on the right ordinate. Its precipitation range was 300–750 °C, which covered the second and third mass loss stages, also corresponded to the peaks of infrared spectrum at 410 °C and 537 °C.
For the detection result of water molecules in Figure 11, a trough appeared before 250 °C, which corresponded to the first peak in the DTG curve, verifying that this was the drying stage. In addition, due to the pyrolysis reaction of the resin, the peak of water molecules appeared at 406 °C, which was ahead of the precipitation of CO, indicating that the decomposition of the epoxy resin would preferentially produce water molecules, but mainly precipitated CO. At this stage, the water produced by the thermal decomposition of epoxy resin mainly come from the breaking of the bonds of -OH functional groups.
For the gas mass-to-charge ratio of 44, CO2, N2O and C3H8 were considered. However, the peak of the curve occurred at 410 °C, ahead of 435 °C, which was the precipitation temperature of the CO peak. Considering the sequence of reactions, CO mainly comes from the reforming of the C-O group in the bisphenol A epoxy resin, and then combined with oxygen ions to form CO2. Therefore, it was speculated that the m/z = 44 mass spectrum curve was mainly represented by the C3H8 gas molecule. Since the reaction was carried out under an inert atmosphere, oxygen ions were supplied only by the material itself, which was also verified in the order of magnitude. It was speculated that the relatively small increasement of the amplitude of this curve afterwards might correspond to the formation of CO2.
The peak temperature of C2H2 gas molecules was close to the temperature of the third peak in the DTG curve, which was supposed to be the product of the second step of decomposition. It can be verified that the decomposition of carbon fiber reinforced epoxy resin matrix composite under inert atmosphere exhibits at least two steps (a major step and a secondary step).

3.3.3. Effect of the Heating Rate on Gas Phase Products

The increase of the heating rate leads to the phenomenon of thermal hysteresis, which is specifically manifested as the peak temperatures of CO, H2O, CO2 and CH4 ion flow all move toward the high temperature region with the increase of the heating rate. The production of H2O, CO2 and CH4 all increase with the increase of heating rate as shown in Figure 12, which indicates that faster temperature rise promotes the reaction. The production of CO had a relatively obvious weak value at the low heating rate, but had little change at the heating rates of 10 °C, 20 °C and 40 °C/min.

4. Conclusions

Thermal degradation of a new product of carbon fiber reinforced epoxy resin composite, namely BA3202, occurs in three steps under inert atmosphere, but in four steps under air atmosphere, respectively. The first two steps were almost the same. The first step is the physical changes, mainly water drying and carbon dioxide escape. The second step is the decomposition of the epoxy resin matrix, and the products are mainly methane, water, carbon monoxide, carbon dioxide and phenol. The third step in the inert atmosphere is the further pyrolysis of the polymer, in which phenol is formed, methane decreases, carbon monoxide basically disappears and carbon dioxide production increases. However, in air, the thermal oxidation of the carbonaceous residues and the intermolecular carbonization are observed. Furthermore, the last step in air is the oxidative decomposition of the carbon fibers. In addition, it is concluded that faster temperature rise can promote the reaction by the experimental comparison in different heating rates. These results provide a comprehensive and thorough understanding into the thermal behavior of CFRP in the thermal environment or fire, and can be used to support material selection and airworthiness certification of aviation structure.
Moreover, the activation energy value obtained by mechanism F4 is closed to that of the Kissinger/Friedman/Ozawa’s method. Therefore, it can be inferred that the thermal degradation reaction mechanism of CFRP for this study obeys the F4 model, which indicates the order of reaction is 4. As a result, the kinetic equation is obtained to support the application of carbon fiber reinforced epoxy matrix composite, such as modeling and assessment of thermal protection, fire resistance and lightning strikes.

Author Contributions

Methodology, H.L.; formal analysis, N.W.; investigation, X.H.; resources, H.L.; data curation, X.H.; writing—original draft preparation, H.Y.; writing—review and editing, H.L.; supervision, J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Al-Rubaye, M.; Manalo, A.; Alajarmeh, O.; Ferdous, W.; Lokuge, W.; Benmokrane, B.; Edoo, A. Flexural behaviour of concrete slabs reinforced with GFRP bars and hollow composite reinforcing systems. Compos. Struct. 2020, 236, 111836. [Google Scholar] [CrossRef] [Green Version]
  2. Mohammed, A.A.; Manalo, A.C.; Ferdous, W.; Zhuge, Y.; Vijay, P.V.; Pettigrew, J. Experimental and numerical evaluations on the behaviour of structures repaired using prefabricated FRP composites jacket. Eng. Struct. 2020, 210, 110358. [Google Scholar] [CrossRef]
  3. Composite Aircraft Structure; FAA AC 20-107B; Federal Aviation Administration: Washington, DC, USA, 2010.
  4. Stoliarov, S.I.; Crowley, S.; Walters, R.N.; Lyon, R.E. Prediction of the burning rates of charring polymers. Combust. Flame 2010, 157, 2024–2034. [Google Scholar] [CrossRef]
  5. Linteris, G.T.; Lyon, R.E.; Stoliarov, S.I. Prediction of the gasification rate of thermoplastic polymers in fire-like environments. Fire Saf. J. 2013, 60, 14–24. [Google Scholar] [CrossRef]
  6. Stoliarov, S.I.; Leventon, I.T.; Lyon, R.E. Two-dimensional model of burning for pyrolyzable solids. Fire Mater. 2014, 38, 391–408. [Google Scholar] [CrossRef]
  7. Mckinnon, M.B.; Ding, Y.; Stoliarov, S.I.; Crowley, S.; Lyon, R.E. Pyrolysis model for a carbon fiber/epoxy structural aerospace composite. J. Fire Sci. 2016, 35, 36–61. [Google Scholar] [CrossRef]
  8. Mouritz, A.P. Review of Smoke Toxicity of Fiber-Polymer Composites Used in Aircraft. J. Aircr. 2009, 46, 737–745. [Google Scholar] [CrossRef]
  9. Speitel, L.C. Toxicity Assessment of Combustion Gases and Development of a Survival Model; DOT/FAA/AR-95/5; Federal Aviation Administration: Washington, DC, USA, 1995. [Google Scholar]
  10. Speitel, L.C. Fractional effective dose model for post-crash aircraft survivability. Toxicology 1996, 115, 167–177. [Google Scholar] [CrossRef]
  11. Tranchard, P.; Samyn, F.; Duquesne, S.; Estèbe, B.; Bourbigot, S. Modelling Behaviour of a Carbon Epoxy Composite Exposed to Fire: Part II—Comparison with Experimental Results. Materials 2017, 10, 470. [Google Scholar] [CrossRef] [Green Version]
  12. 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, 421. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Li, H.; Fan, B.; Wang, N.; Han, X.; Feng, Z.; Guo, S. Thermal response study of carbon epoxy laminates exposed to fire. Polym. Compos. 2020, 41, 4757–4770. [Google Scholar] [CrossRef]
  14. Quintiere, J.G.; Walters, R.N.; Crowley, S. Flammability Properties of Aircraft Carbon-Fiber Structural Composite; DOT/FAA/AR-07/57; Federal Aviation Administration: Washington, DC, USA, 2007. [Google Scholar]
  15. Tranchard, P.; Samyn, F.; Duquesne, S.; Estèbe, B.; Bourbigot, S. Modelling Behaviour of a Carbon Epoxy Composite Exposed to Fire: Part I—Characterisation of Thermophysical Properties. Materials 2017, 10, 494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Tranchard, P.; Duquesne, S.; Samyn, F.; Estèbe, B.; Bourbigot, S. Kinetic analysis of the thermal decomposition of a carbon fibre-reinforced epoxy resin laminate. J. Anal. Appl. Pyrolysis 2017, 126, 14–21. [Google Scholar] [CrossRef]
  17. Brown, M.E.; Maciejewski, M.; Vyazovkin, S.; Nomen, R.; Sempere, J.; Burnham, A.; Opfermann, J.; Strey, R.; Anderson, H.L.; Kemmler, A.; et al. Computational aspects of kinetic analysis: Part a: The ictac kinetics project-data, methods and results. Thermochim. Acta 2000, 355, 125–143. [Google Scholar] [CrossRef]
  18. Coats, A.W.; Redfern, J.P. Kinetic parameters from thermogravimetric data. Nature 1964, 201, 68–69. [Google Scholar] [CrossRef]
  19. Simon, P. The isoconversional method for determination of energy of activation at constant heating rates. J. Anal 1983, 27, 95–102. [Google Scholar]
  20. Friedman, H.L. Kinetics of thermal degradation of char-forming plastics from thermogravimetry: Application to a phenolic plastic. J. Polym. Sci. Part C Polym. Symp. 1964, 6, 183–195. [Google Scholar] [CrossRef]
  21. Opfermann, J.; Kaisersberger, E. An advantageous variant of the ozawa-flynn-wall analysis. Thermochim. Acta 1992, 203, 167–175. [Google Scholar] [CrossRef]
  22. Ozawa, T. A new method of analyzing thermogravimetric data. Bull. Chem. Soc. Jpn. 1965, 38, 1881–1886. [Google Scholar] [CrossRef] [Green Version]
  23. Šesták, J.; Berggren, G. Study of the kinetics of the mechanism of solid-state reactions at increasing temperatures. Thermochim. Acta 1971, 3, 1–12. [Google Scholar] [CrossRef]
  24. Burns, L.A.; Feih, S.; Mouritz, P.A. Fire-under-load Testing of Carbon Epoxy Composites. In Proceedings of the 47th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 5–8 January 2009. [Google Scholar]
  25. Chippendale, R.D. Modelling of the Thermal Chemical Damage Caused to Carbon Fibre Composites. Ph.D. Thesis, University of Southampton, Southampton, UK, 2013. [Google Scholar]
  26. Deng, J.; Xu, L.; Liu, J.; Peng, J.; Han, Z.; Shen, Z.; Guo, S. Efficient method of recycling carbon fiber from the waste of carbon fiber reinforced polymer composites. Polym. Degrad. Stab. 2020, 182, 109419. [Google Scholar] [CrossRef]
  27. Wang, B.; Wang, X.; Xu, N.; Shen, Y.; Lu, F.; Liu, Y.; Huang, Y.; Hu, Z. Recycling of carbon fibers from unsaturated polyester composites via a hydrolysis-oxidation synergistic catalytic strategy. Compos. Sci. Technol. 2021, 203, 108589. [Google Scholar] [CrossRef]
  28. Ahamad, T.; Alshehri, S.M. Thermal degradation and evolved gas analysis: A polymeric blend of urea formaldehyde (UF) and epoxy (DGEBA) resin. Arab. J. Chem. 2013, 7, 1140–1147. [Google Scholar] [CrossRef] [Green Version]
  29. Li, G.; Wu, F. Flame-Retardant Mechanism of Epoxy Resin Treated with Amino Resinous Intumescent Flame Retardant. Plast. Sci. Technol. 2010, 15–18. [Google Scholar]
  30. Liu, Y.; Du, Z.; Zhang, C.; Li, H. Curing Reaction and Thermal Degradation of Bisphenol A Type Novolac Epoxy Resin with 4, 4′-Diaminodiphenyl Ether. Chin. J. Appl. Chem. 2007, 3, 256–260. [Google Scholar]
  31. Huang, N.; Liu, L.; Wang, X. Pyrolysis and Kinetics of Phenolic Resin by TG-MS Analysis. Aerosp. Mater. Technol. 2012, 2, 109–112. [Google Scholar]
  32. Thomas, A.; Moinuddin, K.; Tretsiakova-Mcnally, S.; Joseph, P. A Kinetic Analysis of the Thermal Degradation Behaviours of Some Bio-Based Substrates. Polymers 2020, 12, 1830. [Google Scholar] [CrossRef]
  33. Balart, R.; Garcia-Sanoguera, D.; Quiles-Carrillo, L.; Montanes, N.; Torres-Giner, S. Kinetic Analysis of the Thermal Degradation of Recycled Acrylonitrile-Butadiene-Styrene by non-Isothermal Thermogravimetry. Polymers 2019, 11, 281. [Google Scholar] [CrossRef] [Green Version]
  34. Franco-Urquiza, E.A.; May-Crespo, J.F.; Escalante Velázquez, C.A.; Mora, R.P.; González García, P. Thermal Degradation Kinetics of ZnO/polyester Nanocomposites. Polymers 2020, 12, 1753. [Google Scholar] [CrossRef] [PubMed]
  35. Krawiec, P.; Wargula, L.; Małoziec, D.; Kaczmarzyk, P.; Dziechciarz, A.; Czarnecka-Komorowska, D. The Toxicological Testing and Thermal Decomposition of Drive and Transport Belts Made of Thermoplastic Multilayer Polymer Materials. Polymers 2020, 12, 2232. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, X.; Huang, R. Thermal Decomposition Kinetics of Basalt Fiber-Reinforced Wood Polymer Composites. Polymers 2020, 12, 2283. [Google Scholar] [CrossRef]
  37. GB/T 19267.12-2008. Physical and Chemical Examination of Trace Evidence in Forensic Sciences—Part 12: Thermoanalysis; China National Standardization Administration: Beijing, China, 14 August 2008. [Google Scholar]
  38. ASTM E2550-17. Standard Test Method for Thermal Stability by Thermogravimetry; ASTM International: West Conshohocken, PA, USA, 2017. [Google Scholar]
  39. GB/T 19267.1-2008. Physical and Chemical Examination of Trace Evidence in Forensic Sciences—Part 1: Infrared Absorption Spectrometry; China National Standardization Administration: Beijing, China, 14 August 2008. [Google Scholar]
  40. GB/T 19267.7-2008. Physical and Chemical Examination of Trace Evidence in Forensic Sciences—Part 7: Gas Chromatography/Mass Spectrometry; China National Standardization Administration: Beijing, China, 14 August 2008. [Google Scholar]
  41. Marker, T.R.; Speitel, L.C. Development of a Laboratory-Scale Test for Evaluating the Decomposition Products Generated Inside an Intact Fuselage during a Simulated Postcrash Fuel Fire; DOT/FAA/AR-TN07/15; Federal Aviation Administration: Washington, DC, USA, 2008. [Google Scholar]
  42. Shao, J.; Yan, R.; Chen, H.; Wang, B.; Lee, D.; Liang, D. Pyrolysis Characteristics and Kinetics of Sewage Sludge by Thermogravimetry Fourier Transform Infrared Analysis. Energy Fuels 2008, 2, 38–45. [Google Scholar] [CrossRef]
  43. Ding, H.; Zhang, J.; Cai, P. Thermal Degradation of DGEBA/EDA Epoxy Resin. Polym. Mater. Sci. Eng. 2011, 27, 83–85. [Google Scholar]
  44. Yang, X.; Chen, R. Study on curing kinetics of bisphenol a epoxy resin by pyrolysis gas chromatography (I)—Determination of epoxy value. China Adhes. 1987, 3, 24–25. [Google Scholar]
  45. Feih, S.; Mouritz, A.P. Tensile properties of carbon fibres and carbon fibre–polymer composites in fire. Compos. Part A 2012, 43, 765–772. [Google Scholar] [CrossRef] [Green Version]
  46. Xu, Y.; Yang, Y.; Zhang, Y.; Wang, Z. Pyrolysis characteristics of unidirectional carbon fiber/epoxy prepreg. Acta Mater. Compos. Sin. 2018, 35, 2442–2448. [Google Scholar]
  47. Kissinger, H.E. Reaction Kinetics in Differential Thermal Analysis. Anal. Chem. 1957, 29, 1702–1706. [Google Scholar] [CrossRef]
  48. Bai, Y. Material and Structural Performance of Fiber-Reinforced Polymer Composites at Elevated and High Temperatures. Ph.D. Thesis, School of Architecture, Civil and Environmental Engineering, Lausanne, Switzerland, 2009. [Google Scholar]
  49. Lee, J.Y.; Shim, M.J.; Kim, S.W. Thermal decomposition kinetics of an epoxy resin with rubber-modified curing agent. J. Appl. Polym. Sci. 2001, 81, 479–485. [Google Scholar] [CrossRef]
Figure 1. Specimen in the crucible and placement in the equipment.
Figure 1. Specimen in the crucible and placement in the equipment.
Polymers 13 00569 g001
Figure 2. TG and DTG curves for the carbon/epoxy composite material carried out under inert and air atmosphere at 20 °C/min.
Figure 2. TG and DTG curves for the carbon/epoxy composite material carried out under inert and air atmosphere at 20 °C/min.
Polymers 13 00569 g002
Figure 3. (a) TG curves and (c) DTG curves under inert atmosphere; (b) TG curves and (d) DTG curves under air atmosphere for T300/Epoxy composites at 5, 10, 20 and 40 °C/min.
Figure 3. (a) TG curves and (c) DTG curves under inert atmosphere; (b) TG curves and (d) DTG curves under air atmosphere for T300/Epoxy composites at 5, 10, 20 and 40 °C/min.
Polymers 13 00569 g003aPolymers 13 00569 g003b
Figure 4. Plots of ln ( β / T m 2 ) 1 / T m by Kissinger method under inert and air atmosphere.
Figure 4. Plots of ln ( β / T m 2 ) 1 / T m by Kissinger method under inert and air atmosphere.
Polymers 13 00569 g004
Figure 5. Plots of ln ( d α d t ) -1/T by the Friedman method under inert atmosphere.
Figure 5. Plots of ln ( d α d t ) -1/T by the Friedman method under inert atmosphere.
Polymers 13 00569 g005
Figure 6. Activation energy as a function of fractional mass loss by the Friedman method under inert atmosphere.
Figure 6. Activation energy as a function of fractional mass loss by the Friedman method under inert atmosphere.
Polymers 13 00569 g006
Figure 7. Plots of log β -1/T by the Ozawa method under inert atmosphere.
Figure 7. Plots of log β -1/T by the Ozawa method under inert atmosphere.
Polymers 13 00569 g007
Figure 8. Activation energy and pre-exponential factor as function of fractional mass loss by Ozawa method under inert atmosphere.
Figure 8. Activation energy and pre-exponential factor as function of fractional mass loss by Ozawa method under inert atmosphere.
Polymers 13 00569 g008
Figure 9. Spectra collected at 244, 410 and 537 °C of the degradation products of composites released during TGA (20 °C/min–inert atmosphere).
Figure 9. Spectra collected at 244, 410 and 537 °C of the degradation products of composites released during TGA (20 °C/min–inert atmosphere).
Polymers 13 00569 g009
Figure 10. Intensity curves of m/z = 28 and 94 degradation products of composite released during TGA (20 °C/min-inert gas).
Figure 10. Intensity curves of m/z = 28 and 94 degradation products of composite released during TGA (20 °C/min-inert gas).
Polymers 13 00569 g010
Figure 11. Intensity curves of m/z = 16, 18, 26 and 44 degradation products of the composite released during TGA (20 °C/min-inert atmosphere).
Figure 11. Intensity curves of m/z = 16, 18, 26 and 44 degradation products of the composite released during TGA (20 °C/min-inert atmosphere).
Polymers 13 00569 g011
Figure 12. (a) CO, (b) H2O, (c) CO2 and (d) CH4 release against temperature with the heating rate during pyrolysis.
Figure 12. (a) CO, (b) H2O, (c) CO2 and (d) CH4 release against temperature with the heating rate during pyrolysis.
Polymers 13 00569 g012
Table 1. Kinetic parameters by Kissinger, Friedman and Ozawa methods.
Table 1. Kinetic parameters by Kissinger, Friedman and Ozawa methods.
Method α E A / ( kJ / mol ) A /min‒1
Kissinger 264.6713.088 × 1024
216.5083.14 × 1013
Friedman0.2211.9571.343 × 1014
0.3215.000
0.4217.153
0.5222.009
0.6221.851
0.7285.503
Ozawa0.2181.6701.775 × 1014
0.3194.2799.12 × 1014
0.4198.6115.1 × 1014
0.5202.8172.89 × 1014
0.6208.2231.56 × 1014
0.7225.2637.54 × 1013
Table 2. Kinetic parameters by Kissinger method.
Table 2. Kinetic parameters by Kissinger method.
Reactions Step 1Step 2
AtmosphereArAirArAir
Activation energyEA/(kJ·mol−1)264.671301.884216.508313.307
Pre-exponential factorA/min−13.088 × 10241.96 × 10283.14 × 10131.75 × 1021
Ordern0.940.930.960.945
Table 3. Thirteen kinds of commonly used pyrolysis kinetics models.
Table 3. Thirteen kinds of commonly used pyrolysis kinetics models.
Mechanism F ( α ) Reaction Mechanism
F1 ln ( 1 α ) One-dimensional random nuclear reaction
F2 1 1 α Two-dimensional random nucleation reaction
F3 1 ( 1 α ) 2 Three-dimensional random nuclear reaction
Fn [ ( 1 α ) 1 n 1 ] / ( n 1 ) nth order
A2 [ ln ( 1 α ) ] 1 / 2 Avrami random nucleation n = 2
A3 [ ln ( 1 α ) ] 1 / 3 Avrami random nucleation n = 3
A4 [ ln ( 1 α ) ] 1 / 4 Avrami random nucleation n = 4
R2 1 ( 1 α ) 1 / 2 Phase boundary reaction, cylindrical symmetry
R3 1 ( 1 α ) 1 / 3 Phase boundary reaction, spherical symmetry
D1 α 2 One-dimensional diffusion
D2 ( 1 α ) ln ( 1 α ) + α Two-dimensional diffusion, cylindrical symmetry, (Valensi equation)
D3 [ 1 ( 1 α ) 1 / 3 ] 2 Three-dimensional diffusion, spherical symmetry Jander equation
D4 ( 1 2 α / 3 ) ( 1 α ) 2 / 3 Three-dimensional diffusion, spherical symmetry, Ginstling–Brounshtein equation
Table 4. Correlation coefficients corresponding to various mechanism function for the plot ln(F/T2)-1/T at different heating rates.
Table 4. Correlation coefficients corresponding to various mechanism function for the plot ln(F/T2)-1/T at different heating rates.
Activation Energy (kJ/mol)β = 5 °C/minβ = 10 °C/minβ = 20 °C/minβ = 40 °C/minR
F185.900990.205993.732197.1478−0.98982
F246.078748.159250.015952.5525−0.98913
F3103.327107.678111.615116.918−0.99107
F4199.137207.792215.374224.461−0.99922
A237.365439.423341.074542.6674−0.98643
A321.186922.495823.5218924.5073−0.98121
A413.097614.032014.745615.4273−0.97261
R272.419176.184979.219281.9901−0.98259
R376.729180.668483.860486.8359−0.98531
D1132.311139.064144.484149.150−0.97677
D2147.134154.499160.467165.816−0.98217
D3164.628172.696179.304185.485−0.98728
D4152.932160.530166.711172.335−0.98407
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, H.; Wang, N.; Han, X.; Yuan, H.; Xie, J. Mechanism Identification and Kinetics Analysis of Thermal Degradation for Carbon Fiber/Epoxy Resin. Polymers 2021, 13, 569. https://doi.org/10.3390/polym13040569

AMA Style

Li H, Wang N, Han X, Yuan H, Xie J. Mechanism Identification and Kinetics Analysis of Thermal Degradation for Carbon Fiber/Epoxy Resin. Polymers. 2021; 13(4):569. https://doi.org/10.3390/polym13040569

Chicago/Turabian Style

Li, Han, Nasidan Wang, Xuefei Han, Haoran Yuan, and Jiang Xie. 2021. "Mechanism Identification and Kinetics Analysis of Thermal Degradation for Carbon Fiber/Epoxy Resin" Polymers 13, no. 4: 569. https://doi.org/10.3390/polym13040569

APA Style

Li, H., Wang, N., Han, X., Yuan, H., & Xie, J. (2021). Mechanism Identification and Kinetics Analysis of Thermal Degradation for Carbon Fiber/Epoxy Resin. Polymers, 13(4), 569. https://doi.org/10.3390/polym13040569

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop