3. Results and Discussion
3.1 Non-Isothermal Curing Behavior of T300 Prepreg
To evaluate the curing behavior of the T300 carbon fiber/epoxy resin prepreg, non-isothermal DSC analyses were conducted under four heating rates: 5, 10, 15, and 20°C/min. The DSC thermograms exhibited a single prominent exothermic peak for each heating rate (Fig. 1a), indicating that the curing process follows a one-step dominant reaction. At low temperatures, minor fluctuations in the baseline were observed, corresponding to the glass transition of the uncured resin, where the matrix transitions from a rigid to a viscoelastic state.
After the first dynamic scan, a second scan was performed to confirm reaction completion. As shown in Fig. 1b, no exothermic peak was observed, indicating that the matrix was fully cured during the first cycle. However, a weak endothermic peak appeared near 100°C, which may be attributed to the melting or evaporation of low molecular weight additives or residual volatiles.
The average glass transition temperature of the fully cured matrix (Tg,∞) was determined to be 108.1°C in accordance with ASTM D3418. Curing enthalpies, calculated from the integrated area under the exothermic peaks, varied with heating rate: 115.4 J/g at 5°C/min, 78.85 J/g at 10°C/min, 164.22 J/g at 15°C/min, and 105.0 J/g at 20°C/min. These variations highlight the competition between kinetic and thermal diffusion effects during the curing process.
3.1.1 Kinetic Model Identification and Reaction Mechanism
In order to clarify the mechanism underlying the curing behavior, phenomenological kinetic models were employed. The curing process of thermosetting epoxy resins can typically be described by two types of models: n-level reaction model and autocatalytic reaction15. In the former, the maximum reaction rate occurs at the onset of the process, while in the latter, a delay is observed due to catalytic accumulation.
Analysis of the DSC curves as shown in Fig.
1 suggested that the T300 system exhibits an autocatalytic curing behavior, as the maximum rate is not observed at the beginning. Accordingly, the reaction model is described by:
where α is the degree of cure, A is the pre-exponential factor, Ea is the activation energy, and m and n are empirical reaction orders.
3.1.2 Activation Energy from Kissinger and Ozawa Approaches
To quantitatively determine the apparent activation energy, the Kissinger and Ozawa methods were applied to the non-isothermal DSC data. The Kissinger and Ozawa methods are classic kinetic approaches in thermal analysis, both capable of calculating the activation energy of curing reactions. The Kissinger method assumes:
The Kissinger equation is derived through a series of integrations and calculations based on this Eq. 1
6:
In the equation, Tp represents the peak temperature of the DSC curve, β represents the heating rate of the reaction, A represents the pre-exponential factor, Ea represents the activation energy of the reaction, and R = 8.314. Plotting ln(β/Tp2) against 1/Tp and performing linear regression, the intercept of the fitted line can be used to calculate the pre-exponential factor A, while the slope can be used to calculate the activation energy Ea of the curing reaction.
The Kissinger fitting results are depicted by the blue line in Fig. 2. The slope and intercept of the fitting curve are 9.33×103 and 11.82, respectively, indicating a well-established linear relationship. Based on the fitted data, the activation energy and pre-exponential factor are calculated to be 77.57 kJ/mol and 1.27×109 min− 1, respectively.
The Ozawa method can similarly be employed to determine the activation energy of curing kinetics, with its expression as follows
17–19:
During the temperature rise of the material, the transformation from temperature to time is obtained:
In the equation,
β represents the heating rate. By performing variable integration of this equation and then substituting it into the conversion rate function integral, the Ozawa equation can be derived based on the approximate solution of the Boltzmann factor integral and the Doyle approximation:
Based on the above equation, plotting ln(β) against 1/Tp and performing linear regression yields the fitting results shown by the red line in Fig. 2, with a slope of 1.02×104. The activation energy of the reaction is calculated to be 80.61 kJ/mol. As the Ozawa equation involves the kinetic model f(α), it is not possible to calculate the pre-exponential factor A using the slope of the fitted curve in the absence of a clear kinetic model.
These close results, obtained from both the Kissinger and Ozawa approaches, confirm the reliability of the model-fitting methods for capturing the global curing kinetics of this system.
3.1.3 Conversion-Dependent Activation Energy from Friedman Analysis
Both the Kissinger and Ozawa methods can only yield a single value for the activation energy. However, during the curing reaction, chemical reactions proceed continuously, and the activation energy of the reaction varies at different degrees of cure. Therefore, the Friedman method can be further employed to calculate the activation energy at different degrees of cure in the reaction process. The expression for this method is as follows
20,21:
Figure 3 depicts the curves of the degree of cure as a function of temperature at different heating rates. It can be observed that with higher heating rates, the temperature required to reach the same degree of cure is higher.
A
The data from Fig.
3 corresponding to various heating rates were analyzed by extracting
ln(dα/dt) values at fixed conversion levels ranging from 0.1 to 0.9, and plotting them against the reciprocal of the corresponding temperatures (1/
Tα). Linear regression was then applied to these plots to determine the apparent activation energy as a function of the degree of cure. The results, presented in Fig.
4, reveal a continuous decline in activation energy throughout the curing process. Initially, the activation energy is 75.03 kJ/mol at α = 0.1. As the curing process proceeds, the resin matrix gradually transfers from a rigid state to a viscoelastic gel, leading to the enhanced molecular mobility and improved collisions between epoxy groups and curing agents. As a result, the activation energy steadily decreases with increasing conversion, reaching the minimum value of 35.66 kJ/mol at α = 0.9.
3.1.4 Model Validation via Málek and Sestak–Berggren Approach
Although the Kissinger, Ozawa, and Friedman methods provide key kinetic parameters such as the activation energy (
Ea) and pre-exponential factor (
A), they are insufficient for identifying the specific reaction model
f(α) governing the curing process of T300 carbon fiber prepreg. To determine the type of curing reaction in the resin system, given the activation energy (
Ea), the Málek method can be employed to determine which kinetic equation is suitable for describing the reaction process in the resin system. The Málek method primarily relies on two functions,
y(α) and
z(α) 22,23:
A
where
x represents
Ea/
RT, where
Ea is the activation energy,
R is the gas constant,
T represents the absolute temperature of the reaction,
β represents the heating rate of the reaction, and
π(x) denotes the expression for temperature integral. Senum and Yang proposed that
π(x) can be approximated using the following fourth-order function
24,25:
The average activation energy obtained from the Friedman method was substituted into the functions y(α) and z(α) for calculation and normalization. The relationship among y(α), z(α) and α is shown in Fig. 5.
The values of αp∞, αM, and αp at various heating rates β, as calculated using the Málek method, are summarized in Table 1. The respective degrees of cure corresponding to the maximum values of the functions y(α) and z(α) are denoted as αM and αp∞. Based on Málek 's research 26, when αp∞ ≠ 0.632 and αM∈(0, αp), the curing reaction conforms to the self-catalytic Sestak-Berggren model. According to these criteria, it is observed that the curing kinetics of the epoxy resin system studied in this paper can be described by the self-catalytic Sestak-Berggren model.
Table 1
The values of αp∞, αM, and αp obtained from Malek model in terms of β.
Heating rates β (°C/min) | αp∞ | αM | αp |
|---|
5 | 0.380 | 0.367 | 0.372 |
10 | 0.377 | 0.374 | 0.381 |
15 | 0.353 | 0.346 | 0.350 |
20 | 0.360 | 0.348 | 0.353 |
The values of
αp∞,
αM, and
αp calculated using the Málek method at different heating rates
β can be utilized to establish dynamic kinetic analysis on the fundamental reaction rate equation.
Sestak and Berggren described the reaction mechanism function.
According to Arrhenius law:
Substituting Eqs. (
12) and (
13) into the Eq. (
11) yields the dynamic kinetic calculation formula:
After simplification, the equation can be derived as follows:
In the equation, m and n represent two reaction orders. The reaction order n can be calculated by plotting ln[(dα/dt)ex] against ln[αp(1-α)] and fitting the slope, while the relationship between reaction orders m and n is given by m = pn, where p = αM/(1-αM) 27. The fitting calculation results of A, reaction orders m and n under different heating rates are summarized in Table 2.
Table 2
The relevant kinetic parameters of the calculated curing reaction of the prepreg
Heating rate (°C/min) | A | m | n |
|---|
5 | 0.7036 | 0.8533 | 1.361 |
10 | 1.043 | 0.8104 | 1.292 |
15 | 1.592 | 0.8053 | 1.433 |
20 | 1.908 | 0.7936 | 1.441 |
Figure 6 presents a comparison between the experimental DSC data and the kinetic simulation results. The fitted curves based on the Sestak–Berggren model exhibit good agreement with the experimental data across different temperatures, indicating that this model accurately captures the curing behavior of the epoxy resin system studied. The consistency between the calculated and measured reaction rates further validates the applicability of the Sestak–Berggren model to describe the curing kinetics.
3.2 Isothermal Curing and Diffusion-Controlled Mechanism
3.2.1 Isothermal DSC Analysis of Prepreg
To further examine the cure behavior under practical manufacturing conditions, isothermal DSC experiments were conducted at 120°C, 130°C, 140°C, and 150°C-all above the Tg,∞.
The results of the first round of isothermal DSC scans are shown in Fig. 7. The initial scans reveal that the prepreg system requires an induction period at the beginning of isothermal curing, followed by a gradual increase in the rate of exothermic curing until the reaction is complete. It is evident from the graph that the curing reaction is exothermic, as evidenced by the appearance of an exothermic peak in all DSC curves.
Figure 8 illustrates the variation of curing degree over time. It can be observed that the curing rate increases with the increase of isothermal curing temperatures. At 120°C, 130°C, 140°C, and 150°C, the time required to reach 90% curing degree is 33.4 min, 22.3 min, 18.5 min, and 14.4 min, respectively.
3.2.2 Application of the Diffusion-Controlled Kamal Model
Kamal et al. investigated the curing kinetics of epoxy resins and proposed the Kamal model as follows
28:
The variables m and n represent the reaction orders, while k1 and k2 denote the rate constants for uncatalyzed and catalyzed reactions, respectively.
To further examine the influence of diffusion control on the curing system during the later stages of the reaction, researchers introduced a diffusion factor based on the free volume theory
29:
In the equation,
C represents the restriction constant, and
αc represents the critical degree of cure. By incorporating this diffusion factor into the Kamal model, the diffusion-controlled Kamal equation can be expressed as follows:
The relationship between the degree of curing and time at 120°C, as shown in Fig. 9, was obtained by fitting the experimental data using the diffusion-controlled Kamal equation (Eq. (18)). The model was successfully fitted to isothermal experimental data, as shown in Fig. 9. Table 2 presents the fitted data for the relationship between cure rate and degree of cure. The simulation accurately reproduced the transition from chemically controlled to diffusion-limited regimes, demonstrating that the modified Kamal model can effectively describe the complete curing behavior of the prepreg under isothermal conditions.
Table 2
Fitted data for the relationship between cure rate and cure degree
Temperature (°C) | K1 | K2 | m | n | C | |
|---|
120 | -0.00105 | 0.218 | 0.4325 | 0.7578 | 1.393 | 0.3815 |
130 | 0.02884 | 0.6018 | 0.6495 | 0.7393 | 2.824 | 0.1263 |
140 | -0.04209 | 0.4739 | 0.2813 | 0.5526 | 4.636 | 0.431 |
150 | -0.08672 | 0.8402 | 0.291 | 0.7765 | 4.377 | 0.4327 |