A Study on Curing Kinetics of Nano-Phase Modified Epoxy Resin

In this paper, DSC curves at different heating rates were measured by DSC, the characteristic curing temperature was determined, and the optimum curing conditions were obtained. The KAS method, Friedman method and FWO method were used to analyze the DSC curves respectively. The kinetic parameters and the reaction mechanism function of the curing system were obtained, and the results of different analytical methods were compared and analyzed. Result from fitting and verification of the curing kinetic model for the curing system of the nano-phase modified epoxy resin further demonstrates that the nano-particle could play a catalytic role in the curing reaction of the epoxy resin and could reduce the apparent activation energy of the system, thus it is considered as a breakthrough in the field of resin research.

the physical changes and chemical reactions, by means of a certain mathematical approach to obtain the corresponding reaction of the kinetic parameters and reaction mechanism [15][16][17][18][19][20][21] . There are two methods of thermal analysis kinetics: isothermal and non-isothermal methods. The non-isothermal method has many advantages over the conventional isothermal method 22 . So it gradually becomes the main method of thermal analysis kinetics. Non-isothermal method is divided into the single scanning rate non-isothermal method and the multiple scanning rate non-isothermal method 23,24 .
Multiple scanning rate non-isothermal method refers to using different heating rates measured under a number of dynamic curves to carry out the method of dynamic analysis 25,26 . According to the different mathematical methods, it can be divided into two categories of differential and integral methods. At present, in the method of research dynamics, the differential method is mainly represented by Kissinger-Akahira-Sunose (KAS) method and Friedman method, and the integral method is mainly represented by Flynn-Wall-Ozawa (FWO) method 20,23,[27][28][29][30] .
The KAS method is based on the assumption that the maximum rate of the curing reaction occurs at the peak temperature of the curing reaction exothermic peak, so the KAS method is also called as the maximum rate method 31 . Friedman method and FWO method are used in a same conversion of several kinetic curves of the data. They are is also called as the Isoconversional method, which can obtain more reliable values of the activation energy E in the premise of not involving the kinetic model function eliminating the effects of model function in calculating Arrhenins parameters. So it is also known as a model free method. As the results show that the multiple scanning rate method and the equal conversion method can describe the complex solid -state reaction kinetics more effectively.  The epoxy resin kinetic model parameters are typically obtained by differential scanning calorimetry (DSC) 32 . The advantages of this method include use of less sample, convenient operation and strong ability of data analysis, particularly suitable for thermosetting resin curing reaction kinetics of the study.
The 21st century features robust development of macromolecular materials, among which the epoxy resin is the pilot in the development and application in this area. Epoxy resin is widely used in electrics, semiconductor electronics, special coating aeronautics, astronautics, etc., owing to its excellent properties. The curing kinetics analysis on epoxy resin is helpful to better understand the characteristics of its system, the curing conditions and the impact imposed on the curing activities by the dimensions of curing molds. It can also help to establish a quantitative relation between the curing conditions, epoxy resin, as well as the chemical structure of curing agents and the properties of curing objects to reduce the experiments in order to optimize the curing process.
As an amorphous viscous liquid, epoxy resin is plastic when heated, without clear melting point, becoming soft when heated, turning sticky when gradually melted, insoluble in water, and will not harden by itself. Therefore, it almost has no value to use independently. It must first react with curing agents to create a polymer with a three-dimensional structure which does not melt or dissolve, and then it can be usable. When a certain  amount of curing agents is added, it gradually hardens and turns into chemical materials with distinct properties. Therefore, curing agents must be added to react under certain conditions to turn into a cross-linked network structure in order to create formula resin. In this way, the epoxy material can display excellent properties and be truly usable. The common amine curing agents are: ethylenediamine, diethylene polyamine, and polyethylene polyamine, etc.

Experimental Study
Curing Principle. Ammonium hexafluoroantimonate is a kind of Cationic curing agent by heated, which could trigger low temperature and curing speed, and dissolve in epoxy resin. Curing products have better adhesion and corrosion resistance, anti-aging properties.
The curing mechanism is shown in Fig. 1. Firstly, epoxygroups in the epoxy group react with ammonium hexafluoroantimonate to capture its active hydrogen, to generate oxygenium ions, with hydroxyl groups of oxygen ions as the initiator, and epoxy groups continue to react, to open the way for chain ring polymerization of curing epoxy resin. As long as the epoxy repair agent is in contact with the curing agent, the curing reaction can be automatically carried out. Because the curing  agent can be well dissolved in the epoxy resin, it could ensure that micro-cracks derive from self-healing, and achieve a high polymerization rate and polymerization degree. Figure 2 shows the thermogravimetric curve of the curing agent. It can be seen that the weight loss starting temperature is 254.4 °C, which indicates that the ammonium hexafluoroantimonate curing agent has good thermal stability and meets its requirement of stability in the self-repairing system. Curing agent in the matrix can exist with stability for a long term.   The curing reaction exothermic peak of the epoxy resin and the nano rubber modified epoxy resin system are related to the curing temperature and time in different heating rates, as shown in Figs 3-6, in which the arrow direction is the increasing direction of the heating rate. It can be seen from Figs 3-6 that the curing reaction characteristic temperature of each curing system is closely related to the heating rate. From the curve of the   exothermic peak with time, it can be seen that with the increase of the heating rate, the exothermic peak becomes steeper gradually, and the curing time is obviously shortened. From the change curve of the exothermic peak to the temperature (Figs 3a and 5c), it can be seen that the curing initiation temperature and the peak temperature of the system move toward the high temperature with the increase of the heating rate.   This shows that the curing reaction is not only a thermodynamic process, but also a dynamic process. At a lower heating rate, the reaction system has sufficient time to react, so the reaction starts at the lower temperature. With the increase of the reaction rate, the ratio dH/dt of the system becomes larger, which means the increase of    Table 3. The apparent activation energy and pre-exponential factors of the curing system.       the thermal effect per unit time and the thermal inertia. The greater the temperature difference is, the higher the temperature of the curing reaction exothermic is. After the addition of nano-particles, the exothermic peak can be seen to be sharpened (Figs 5c and 6d), and the addition of the nano-particles to the pure epoxy resin increases from the original bimodal to the single peak. Based on the peak value of heat peak, it can be seen that the amount of heat released during the reaction is obviously reduced and the reaction time is shorter. It shows that the curing reaction becomes gentler and the reaction rate is accelerated. The nano-particles can promote the curing of epoxy resin.

Analysis of Curing System of Curing System.
Kissinger-Akahira-Sunose (KAS) equation is as follows: In the formula, β is the heating rate, T P is the curing peak temperature, E is the apparent activation energy, A is the pre-index factor, R is the universal gas constant. We can use , you can get a straight line from the slope of the straight line, and obtain the intercept from the apparent activation energy E and pre-exponential factors A. According to the DSC curve of the graph and the peak temperature T p at different heating rates for each curing system the linear regression of 1  The E and A of each curing system are obtained from the slope and intercept of the straight line in the figure.
The results are shown in Table 1. The different curing systems E and A can be seen by KAS method. The results show that the E and A of the curing system are obviously decreased after the addition of the nanoparticles, which indicates that the reaction system is more rapid after adding the nanoparticles.

Analysis of Curing System by Friedman Method. Friedman equation is as follows:
In the formula, the physical meaning of β, E, A, R, α is the conversion rate, and f(α) is the differential form of the reaction mechanism function. Friedman equation E and A can be obtained by the following method. The specific steps are as follows: Since the same conversion rate α is selected at different heating rates β, tthe integral form of   Table 2.
It can be seen from Figs 11 and 12 and Table 2 that the activation energy at the beginning of the reaction is slightly higher for the pure epoxy resin, and the apparent activation energy at the low conversion rate stage is greatly reduced because the curing agent needs a process of dissociation that could produce a curing effect, during which to the energy needs to be dissipated. When the conversion rate is more than 20%, the E and A of the system increase with the increase of the conversion rate, and the apparent activation energy deviation also gradually rises. When the conversion rate is more than 90%, the apparent activation energy can be steeply increased, indicating that the epoxy resin has a gel effect, the viscosity of the system increases sharply, and the reaction functional groups close to the gel point are bound in the crosslinking network, the energy required to react with each other between the reaction functional groups increases, so that the latter stage of the curing reaction requires a great amount of energy until the end of the curing. The apparent activation energy of epoxy resin after nano-phase modification is stable, and it is lower than that of pure epoxy resin in different conversion stages. It shows that the curing reaction is gentle, and the presence of nanoparticles for the curing reaction plays a catalytic role. The change trend of lgA is similar to apparent activation energy E.
The E and A data obtained by Friedman analysis in Table 2 are averaged to obtain the apparent activation energy E and the pre-exponential factor A for each solidified system. The results are shown in Table 3. It can be seen that the addition of organic nanoparticles, the apparent activation energy of the system and the pre-exponential factors have a significant reduction, indicating that the nano-particles on the curing play a catalytic role.

Analysis of Curing System by FWO Method. Flynn-Wall-Ozawa (FWO) equation is as follows:
As we know, G(α) is the integral form of reaction mechanism function in the formula. There are formulas that uses FWO method to find E, A. Due to the same conversion rate α is chosen at different heating rates β, G(α) is a constant value in the formula, and in this situation, lgβ and 1 T have a linear relationship. As the Friedman method is adopted, the corresponding slope and intercept can be used to find the E and A values corresponding A. For each α can be a set of E and A values, for all numerical averages, you can determine the reasonable E and A value.
According to the FWO equation, the DSC curve of the graph is analyzed, and the linear regression curves of Figs 13 and 14 are obtained by linear fitting with lgβ and 1 T , linear regression and computer. For each curing system, using the slope and intercept of the regression curve at each conversion rate, the corresponding conversion rates E and A were obtained. The results are shown in Figs 15 and 16 and Table 4.
As it can be seen from Table 4, for each curing system, with the increase in conversion rate, the system E and A as a whole shows an increasing trend, and the apparent activation energy deviation gradually increases. Trends are consistent with Friedman's analysis.
The Eand A values obtained by the FWO method in Table 4 were averaged to obtain the apparent activation energ E and the pre-exponential factor A of the solidified system as a whole. The results are shown in Table 5.
Therefore, FWO method is used to analyze the system of E. The addition of organic nano particles, the apparent activation energy and pre-exponential factor decreases significantly, which indicated that nanoparticles to the same catalytic have an effect on curing up.
Determination of curing reaction series. The reaction series of the curing reaction can be obtained by the Crane equation.
Which the physical meaning of β, T p , E, R is the same as above, the reaction series is n. The above formula can be integrated, According to the formula, lnβ to 1 T p linear fit, the results shown in Fig. 9. As can beseen from Fig. 17, the curing system lnβ has a good linear relationship with 1 T p , according to the Crane equation, the apparent activation energy of the cured system was analyzed by KAS method, Friedman method and FWO method, and the reaction order of curing system was obtained by different methods, as shown in Table 6. It can be seen that the apparent activation energy obtained by the KAS method, the Friedman method and the FWO method are not all integers, indicating that the curing reaction is a complex reaction and the results of the three methods are also correct.
SCiEntifiC RepoRTS | (2018) 8:3045 | DOI:10.1038/s41598-018-21208-0 Model validation. The commonly used kinetic equation for nonhomogeneous systems under non-isothermal conditions is shown as follows: A exp E RT f( ) (6) In this formula, d dT α for the reaction rate, A for the pre-exponential factor, E for the apparent activation energy, f(α) for the reaction model coefficient, R for the ideal gas constant. β = dt/dT for the heating rate. The above equation can be transformed into a linear equation: The equation can be seen as = + y a a x 0 1 , in which = x 1/T can be analyzed using multiple linear programming. The results are shown in Table 7.
Based on the function expression of different reaction models f(α), the computer being used to analyze the data fits with the commonly used dynamic mechanism function to obtain the most suitable curing kinetic model.
Five different reaction models in the table were used to fit the previous data using the multivariate least squares regression method. The correlation coefficient (r) is defined as follows: s k sk 2 k sk 2 s In the formula: Y jk is the measured value, Y jk is the regression value, S is the number of measurements under different conditions, N S is the number of data measured at one time. The experimental data were fitted and the results are shown in Table 8, the results are shown in Table 9.
Comparing the results of 5 commonly used dynamic mechanism function fitting, we can see that the curing reaction of epoxy resin modified by nano phase belongs to the PT kinetic mechanism function.
Comparing on fitting curves and test curves, the result is shown in Figs 18-22, the solid line is the fitting curve, and the dot matrix curve is the actual measurement data.
From the chart it can clearly seen that the curing system of DSC curve to simulate the PT kinetic mechanism functions coincide with test results of DSC. In the meantime, reaction mechanism corresponds to the self catalytic reaction shows that nano modified epoxy resin and six hexafluoroantimonate ammonium curing reaction proceeds in accordance with the autocatalytic reaction and it is consistent with the previous reaction mechanism.

Conclusion
The curing curves of the two epoxy resin systems show that the peak of heating becomes steeper with the heating rate increasing, and the curing time is shorter. The curing characteristic temperature of the system moves to the high temperature.
After the addition of nanoparticles, the heat release was reduced and the time was shortened, indicating that the nanoparticles played a catalytic role in the reaction.
Three kinds of analytical methods were used to study the curing kinetics of the two curing systems, and the kinetic apparent activation energy E and the pre -exponential factor A were obtained.
The results show that the apparent activation energy and pre -exponential factors of the system are greatly reduced after the addition of the nanoparticles, which indicates that the nanoparticles can promote the curing reaction of the system and reduce the apparent activation energy of the system.
The kinetic model of nano-phase modified epoxy resin was fitted by multiple linear analysis method and different reaction mechanism. The results showed that the curing reaction mechanism was autocatalytic reaction.