Processability Map of a Recyclable Thermoplastic for Structural Applications. Kinetics of Induction Period
Sihem Zaidi, Daniel Sánchez-Rodríguez, Jordi Farjas, Raquel Verdejo, Miguel Ángel López-Manchado, Daniel Trias, Josep Costa

TL;DR
This study creates a processability map for an acrylic thermoplastic to guide its use in industrial manufacturing by analyzing its polymerization and decomposition.
Contribution
A novel methodology is introduced to quantify the induction period of polymerization using DSC data.
Findings
A time–temperature–transformation (TTT) map was constructed to define the processing window.
The polymerization process involves an induction period followed by rapid reaction.
The method supports efficient integration of acrylic resins into composite manufacturing.
Abstract
In this study, we conducted a comprehensive kinetic analysis of the polymerization and thermal decomposition of an acrylic thermoplastic to establish a time–temperature–transformation (TTT) processability map. The process exhibits an induction period followed by rapid polymerization, indicating a complex mechanism. We introduce a novel methodology to quantify this induction period using differential scanning calorimetry (DSC) under both isothermal and dynamic conditions. The resulting kinetic parameters were used to construct the TTT map, which delineates a broad and practical processing window. This map enables precise control of the polymerization process, minimizes the risk of thermal degradation, and supports the efficient integration of acrylic resins into industrial composite manufacturing.
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
1
2
3
4
5
6
7
8- —Ministerio de Ciencia, Innovaci?n y Universidades10.13039/100014440
- —Ministerio de Ciencia, Innovaci?n y Universidades10.13039/100014440
- —NextGenerationEU10.13039/100031478
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsThermal and Kinetic Analysis · Epoxy Resin Curing Processes · Fiber-reinforced polymer composites
Introduction
1
Wind energy is essential to achieving Europe’s net-zero carbon target by 2050.? However, the nonrecyclable composite blades of wind turbines are generating growing landfill waste, projected to reach 43 million tons globally by 2050about 25% in Europe. This underscores the urgent need for recyclable blade materials.?
In this work, we investigate the polymerization and thermal degradation kinetics of Akelite, a fully recyclable thermoplastic resin patented by the Institute of Polymer Science and Technology (ICTP-CSIC).? Akelite is a low-viscosity liquid acrylic resin, making it particularly suitable for the fabrication of fiber-reinforced polymer composites. Its low viscosity enables fiber impregnation at room temperature, allowing the use of conventional manufacturing techniques originally developed for thermosetting epoxy-based composites. Moreover, the resin exhibits excellent storage stability, remaining useable for over six months when kept below 25 °C. In addition to its favorable processing characteristics, Akelite demonstrates outstanding mechanical performance, positioning it as a promising eco-friendly alternative for structural applications such as wind turbine blades.? To support its industrial implementation, this study focuses on characterizing both the polymerization and degradation kinetics of the resin.
To characterize the polymerization reaction and thermal degradation kinetics of Akelite, we employ isoconversional methods. These model-free approaches do not rely on predefined reaction mechanisms, making them more versatile and easier to apply than traditional model-based methods. Isoconversional techniques have been successfully used to study both polymerization and degradation processes in a wide range of polymeric systems. ?−? ? ? ? ? ? ? ? ? Furthermore, they provide accurate kinetic predictions, enabling the construction of time–temperature–transformation (TTT) diagrams, which are valuable tools for optimizing processing conditions. ?−? ? ? ?
The main monomer in Akelite is methyl methacrylate (MMA), while benzoyl peroxide is used as an initiator to promote free-radical polymerization. This compound initiates a chain reaction through thermal decomposition, which involves the homolytic cleavage of the O–O bond, generating two free radicals. These radicals react with the monomers, initiating a propagation process in which the polymer chains grow by successive monomer additions. In the presence of a peroxide initiator, the bulk free-radical polymerization of MMA exhibits a characteristic sharp increase in reaction rate following an initial induction period, a phenomenon known as the Trommsdorff or “gel” effect. ?−? ?
Initially, the low viscosity of the reaction medium facilitates rapid termination of macroradicals via recombination, which suppresses the polymerization rate. As the reaction progresses and the degree of polymerization reaches approximately 35–40%, ?,? the viscosity of the system increases significantly. This rise in viscosity restricts the mobility of macroradicals, thereby reducing termination events. Meanwhile, monomer molecules remain sufficiently mobile to sustain propagation. As a result, the reaction rate increases abruptly due to the accumulation of active radicals. The kinetics of the system is therefore mainly governed by the induction period, which corresponds to the time required for the suppression of macroradical termination.
Several authors have modeled the polymerization kinetics of MMA in the presence of initiators. ?,?,?−? ? ? ? However, accurately capturing the Trommsdorff effect remains a significant challenge due to the complexity of the underlying mechanisms and the need to determine numerous physical parameters. Moreover, all existing models incorporate free parameters that must be adjusted to accurately capture the reaction kinetics. Consequently, a considerable degree of empiricism underlies these models. Furthermore, the reaction kinetics is highly sensitive to the specific formulation of the resin, including factors such as initiator concentration, other additional monomers, and the presence of predissolved poly(methyl methacrylate) (PMMA).
Model-free methods generally offer a simplified approach to reaction kinetics analysis. However, as will be demonstrated, the presence of the Trommsdorff effect significantly limits the effectiveness of isoconversional techniques in accurately characterizing the polymerization kinetics of MMA. Despite this limitation, we propose a straightforward methodology to quantify the induction period. Based on this characterization, we have developed a processing map that facilitates the identification of optimal manufacturing conditions while avoiding thermal degradation. These findings provide valuable guidance for wind turbine blade manufacturers, supporting the evaluation and adoption of recyclable materials as viable alternatives to conventional thermosetting systems.
Materials and Methods
2
Materials
2.1
The material under study is a combination of several acrylic monomers with methyl methacrylate (MMA) being the main monomer, and benzoyl peroxide is used as the initiator to promote bulk free-radical polymerization. The resin exhibits a maximum glass transition temperature of approximately 120 °C. During synthesis, the viscosity of the resin is adjusted to around 140 cP at room temperature, making it suitable for composite manufacturing processes via infusion. This resin has been patented by ICTP-CSIC under the name Akelite.?
To investigate the polymerization behavior, the resin is mixed with 3 wt % of initiator. To ensure homogeneity, the mixture is placed in an ultrasonic bath for 15 min. It is then transferred into hermetically sealed aluminum crucibles (40 μL capacity, TA Instruments) to prevent MMA evaporation during thermal analysis. The typical sample mass ranges from 0.5 to 1 mg, which minimizes the risk of self-heating and ensures proper thermal equilibration during differential scanning calorimetry (DSC) experiments.?
Differential Scanning Calorimetry
2.2
The polymerization reaction was monitored using a TA Instruments Q2000 differential scanning calorimeter (DSC). Experiments were conducted under both dynamic and isothermal conditions. Dynamic scans were conducted from 0 to 200 °C at heating rates ranging from 1.25 to 7.5 °C/min. Isothermal experiments were performed at temperatures between 90 and 120 °C, with durations varying from 1 h at 120 °C to 2.5 h at 90 °C. All measurements were conducted under an inert atmosphere, maintained by a continuous flow of high-purity nitrogen at 50 mL/min. An empty aluminum crucible, similar to that of the sample was used as a reference.
In isothermal experiments, the glass transition temperature (Tg) was used as an indicator of the degree of transformation. Following the isothermal step, each sample was subjected to a second DSC scan from 0 to 200 °C at a heating rate of 20 °C/min. The Tg was determined according to standard methodology, using the midpoint between the extrapolated onset and end point of the transition.
Thermogravimetric Analysis
2.3
Thermal degradation of the cured resin was analyzed using a Mettler Toledo thermobalance (model TGA/DSC1) via thermogravimetric analysis (TGA). A series of dynamic experiments were carried out from 50 to 700 °C at heating rates ranging from 1.25 to 20 °C/min, under an inert atmosphere maintained by a constant flow of high-purity nitrogen at 60 mL/min. Approximately 6 mg of sample was used in each test, placed in 150 μL Al_2_O_3_ crucibles.
Kinetic Analysis
2.4
The kinetic parameters of polymerization and degradation, including the apparent activation energy and the pre-exponential factor, were determined using Friedman’s isoconversional method. ?,? Isoconversional methods are based on the assumption that, for a given degree of conversion (α), the reaction rate depends solely on temperature. ?,? The Friedman method is expressed as
where α is the degree of conversion (0 ≤ α ≤ 1), T is the absolute temperature, R is the universal gas constant, and E α is the activation energy corresponding to a specific conversion degree.
For most thermally activated processes, the energy barrier remains nearly constant throughout the transformation.? Therefore, variations in E α should not be interpreted as changes in a single energy barrier, but rather as the result of multiple overlapping processes, each characterized by its own activation energy, or as a deviation from Arrhenius-type behavior. ?,?,? In such complex systems, E α is best interpreted as an apparent activation energy rather than a true intrinsic value.
Friedman method is based on the integration of the eq
Thus, the reaction kinetics can be described by two parameters that depend on the degree of transformation: the apparent activation energy E α and the product Af(α). The first step in determining these parameters is to discretize the degree of conversion using a small, constant increment Δα, such that α_ j _ = jΔα. The objective of the kinetic analysis is then to determine the set of values E _ j _ and Af(α)|_ j . This can be achieved by performing experiments at different heating rates β i dT/dt (nonisothermal conditions) or at different fixed temperatures (isothermal conditions). ?,? For each heating rate β i , the transformation rate dα/dt| j,i _ and the corresponding temperature T _ j,i _ at a given degree of conversion α_ j _, are determined. According to eq, a plot of versus 1/T _ j,i _ should yield a straight line.
From the slope and intercept of the linear fit, the values of E _ j _ and Af(α)|_ j _, can be determined, respectively. Once these parameters are known, the evolution of the curing process under any given temperature program can be predicted by integrating eq
where u is the degree of transformation.
Results and Discussion
3
This section presents the kinetic characterization of Akelite polymerization and thermal degradation. First, the polymerization kinetics is analyzed under both dynamic and isothermal conditions, with particular emphasis on characterizing the induction period associated with the Trommsdorff effect. Subsequently, the thermal degradation kinetics is examined. Finally, a time–temperature-transformation (TTT) processability map is constructed by integrating both processes.
Polymerization Kinetics
3.1
To determine the evolution of the degree of polymerization as a function of temperature, it is assumed that the heat flow released by the sample is directly proportional to the transformation rate. Therefore, the transformation rate dα/dt can be obtained directly from the DSC signal, normalized by the total enthalpy of the reaction?
where ΔH is the enthalpy obtained from the integration of the DSC curve and dH/dt is the heat flow measured by the DSC. Figurea shows the evolution of the rate of transformation and degree of polymerization obtained from DSC measurements made when the sample is heated at a constant rate, specifically five different heating rates have been analyzed: 1.25, 2.5, 3.75, 5, and 7.5 °C/min.
(a) Solid lines represent the rate of transformation (dα/dt) and the degree of polymerization (α) as functions of temperature, obtained by integrating DSC curves. Symbols indicate isoconversional predictions based on eq . (b) Activation energy of the polymerization reaction determined using Friedman isoconversional analysis (left axis), along with the coefficient of determination (R 2) of the linear fit (right axis).
The result of Friedman’s isoconversional analysis, based on eq, is presented in Figureb. From a degree of polymerization of approximately 30% (α = 0.3), the calculated activation energy drops to near zero, and the coefficient of determination shows that the data no longer follows a linear trend. This indicates that beyond this conversion level, the reaction is no longer governed by thermally activated kinetics and cannot be described by the isoconversional approach.
As a self-consistency check of the isoconversional model, eq was used to simulate the polymerization behavior under the same dynamic conditions as those used in the experiments. The results, shown as symbols in Figurea, demonstrate good agreement with the experimental data, despite the limitations of the model beyond α = 0.3. This indicates that the isoconversional approach can still provide reliable predictions for nonisothermal polymerization processes within the explored heating rates range (1.25 to 7.5 °C/min).
The evolution of the polymerization rate, shown in Figurea, reveals two clearly distinguishable stages. The initial stage is characterized by a smooth and gradual increase in the process rate, followed by a sudden and pronounced acceleration. This behavior is a direct consequence of the Trommsdorff effect.
At the beginning of the reaction, the low viscosity of the resin facilitates efficient termination of macroradicals, maintaining a relatively low polymerization rate. However, after an induction period, the viscosity of the system increases significantly, which impedes the mobility of macroradicals and suppresses termination reactions. As a result, the concentration of active radicals rises sharply, leading to an abrupt increase in the polymerization rate. Beyond this point, the reaction is no longer governed by thermal activation, but rather by the accumulation of reactive species due to the inhibited termination.
A noteworthy observation is that the induction period corresponds to a specific degree of conversion. ?,? As shown in Figurea, the reaction kinetics exhibit a sudden acceleration at α = 0.3 across all experiments. This transition point marks the onset of the autoaccelerated regime.
Figureb further highlights this behavior by revealing three distinct kinetic regions. Below α = 0.3, the activation energy remains approximately constant at around 90 kJ/mol, and the high correlation coefficient indicates a strong linear relationship, consistent with thermally activated kinetics. Beyond this point, the activation energy drops sharply, and the correlation coefficient decreases significantly, suggesting a breakdown of the linear relationship between and 1/T _ j,i _. This deviation confirms that the reaction is no longer governed by thermal activation, a feature that can be attributed to the autoaccelerated reaction.
To further investigate the mechanism change, we performed a master plot analysis using the normalized conversion function f(α)g(α)/f(0.5)g(0.5), which was obtained from the experimental data through the relationship dα/dt|α/dα/dt|0.5.? This master curve is characteristic of the transformation mechanism. The results is shown in Figure, where a clear change in its trend is observed at α ≈ 0.3. This deviation indicates a transition in the polymerization mechanism, which is in agreement with the Friedman analysis shown in Figure. The consistency between both independent kinetic analyses reinforces the conclusion that the polymerization process undergoes a mechanistic change at this conversion level, marking the onset of the autoaccelerated stage.
Master plot derived from DSC data at different heating rates (Figure a).
At higher degrees of conversion, particularly above α = 0.9, the correlation coefficient increases once again, approaching unity. This trend suggests a renewed thermal activation of the process, with the activation energy rising from 20 to 40 kJ/mol. Such behavior aligns with previous studies, ?,? which report that at advanced stages of conversion, the reaction rate declines because the reaction becomes diffusion-controlled. ?,?,?
Therefore, to characterize the kinetics, it is essential to determine how the induction period depends on temperature. Given that this induction period corresponds to a fixed degree of transformation during the initial stage, it can be inferred that the induction period itself is thermally activated.
To characterize the induction period, a series of isothermal experiments were conducted, as shown in Figure. According to eq, plotting the logarithm of the reaction onset time against the reciprocal of the absolute temperature yields a straight line. From the linear fit, both the activation energy and the pre-exponential factor can be determined.
Continuous lines represent DSC isothermal measurements, while dashed lines correspond to isothermal isoconversional predictions derived from eq .
The resulting fit, presented in Figure, provides an activation energy of E _ i _ = 118 kJ/mol, a pre-exponential factor of A _ i _ = 4.22 × 10^–14^ s, and a coefficient of determination R ^2^ = 0.995. Based on the linear fit and eq, it is possible to predict the time required to initiate the chain reaction at any temperature. For instance, at 30 °C, the estimated induction time is approximately five months, which aligns with the known stability of Akelite at room temperature.
Plot of the logarithm of the reaction onset time versus the reciprocal of temperature.
Equation can also be used to estimate the temperature at which the chain reaction initiates under nonisothermal conditions. Specifically, the contribution to the induction period during an infinitesimal time interval dt is given by
Considering that the heating rate is βdT/dt, the temperature T _ i _ at which the autoaccelerated stage starts is the solution of the equation
which can be reduced to
where, p(x) is the integral of temperature?
and is the integration variable.
This integral does not have an exact solution, but there is a very precise sixth degree Padé series that allows its analytical evaluation.?
Using eqs and (?), we calculated the temperature T _ i _ corresponding to various heating rates and compared it with the temperature at which an abrupt change is observed in the DSC experiments shown in Figurea. The results are presented in Figure, which demonstrates that the kinetic parameters derived from the isothermal experiments successfully reproduce the behavior observed under constant heating rate conditions.
Logarithm of the heating rate as a function of the temperature at which the autoaccelerated reaction is triggered, as predicted by eqs and () (solid line) and as observed experimentally (symbols).
It is important to emphasize that the line shown in Figure is not a fitted curve, but rather a prediction based on eqs and (?), using the values of E _ i _ and A _ i _ obtained from the isothermal data. The agreement between predicted and observed behavior confirms that the kinetic analysis of the induction period is reliable within the explored temperature range.
It is worth noting that eq bears a formal resemblance to the relationship underlying the Flynn-Wall-Ozawa (FWO) method, ?,? which is widely used in thermal analysis for determining activation energies. The FWO method is based on the general rate equation for thermally activated processes under constant heating rate conditions.
Integration of this eq leads to
where .
When g(α) = 1, this expression becomes formally identical to eq. However, the physical meanings are fundamentally different: the FWO method describes the kinetics of thermally activated reactions through the conversion degree α, whereas eq relates to the induction period preceding the reaction. Despite this conceptual difference, the mathematical similarity proves useful for understanding the linear behavior observed in Figure.
The FWO method achieves linearization through Doyle’s approximation of the temperature integral?
Substituting this approximation and taking the natural logarithm yields
This relationship predicts that a plot of ln β versus 1/T should exhibit a linear trend with slope −1.0518 E/R. By analogy, the same approximation justifies the linear behavior observed in Figure.
Kinetics of Thermal Degradation
3.2
TGA is particularly suitable for characterizing degradation kinetics because it monitors the mass evolution. In this case, it is usually assumed that the degree of degradation is linearly dependent on mass.
where m _ i _ and m _ f _ are the initial and final sample masses, respectively.
Figure a presents the degradation profiles at various heating rates, alongside the results of Friedman’s isoconversional analysis (Figureb). In this case, a good agreement is observed over a broad conversion range (ε > 0,06). Specifically, the comparison between the experimental curves and the predictions obtained from eq, shown in Figurea, confirms that the isoconversional model accurately reproduces the measured behavior. This demonstrates that the isoconversional approach provides a reliable framework for predicting the evolution of degradation under varying thermal conditions.
(a) Solid lines represent the degradation rate (dε/dt) and the degree of degradation (ε) as functions of temperature, obtained from TGA measurements. Symbols indicate the isoconversional predictions derived from eq . (b) Activation energy of the degradation reaction determined via Friedman isoconversional analysis (left axis), along with the coefficient of determination (R 2) of the linear fit (right axis).
Processing Map
3.3
From the kinetic analysis we can construct a transformation time–temperature diagram (TTT) that takes into account the kinetics of polymerization and degradation. For this purpose and with the help of eq, we have determined the time required for the autoaccelerated reaction to occur. We have also included the time required for reaching a certain degree of degradation. The result is presented in Figure. The separation between polymerization and degradation, especially at elevated temperatures, is remarkable. This wide processing window significantly facilitates the use of Akelite resin for the manufacture of laminates. In particular, polymerization at elevated temperatures is in principle feasible for relatively short times.
Isothermal processing map of Akelite resin. Red lines indicate the time required to reach a certain degree of thermal degradation at various temperatures. The blue line represents the incubation time at each temperature. Dots correspond to the isothermal polymerization cycles used to prepare samples A, B, and C.
To evaluate the reliability of the processing map, three samples were subjected to isothermal treatment for a fixed duration of 30 min at different temperatures: 85 °C (sample A), 105 °C (sample B), and 140 °C (sample C). We selected this time interval because, according to the TTT diagram prediction, sample A will not reach the incubation time, whereas samples B and C will. In addition, for sample C, even though the temperature is considerably higher, the prediction indicates that no significant degradation would occur. Samples B and C were selected to demonstrate the wide processing window available for Akelite. These isothermal curing conditions are represented as blue dots in Figure.
Figure shows the DSC signals recorded during these polymersization cycles. Based on eq, the estimated incubation times for samples A, B, and C are 9050, 1100, and 45 s, respectively. For a 30 min curing cycle, sample A does not reach the autoaccelerated regime, as evidenced by the nearly flat DSC signal, indicating an almost negligible transformation rate throughout the cycle. In contrast, samples B and C exhibit sharp exothermic peaks, confirming that the curing process is completed in both cases. Notably, after the curing is complete, the DSC signal remains flat, suggesting that no degradation occurs under these conditions.
Isothermal DSC measurements corresponding to the curing cycles labeled (A) (85 °C), (B) (105 °C), and (C) (140 °C) as shown in Figure .
Visual inspection revealed that sample A had undergone almost no polymerization, whereas samples B and C were clearly solid, indicating a significant degree of polymerization. To confirm these observations and to quantify the extent of polymerization, a second measurement was performed using DSC. This involved heating the samples from room temperature to 200 °C at a rate of 20 K/min to determine their glass transition temperatures (Tg). As expected, no glass transition was detected for sample A. This is consistent with the reported Tg of MMA at −126 °C, ?,? which lies below the experimental temperature range, further indicating that no polymerization had occurred. In the same experiment, however, an exothermic signal was observed between 100 and 200 °C, corresponding to the polymerization process, as expected for a nonpolymerized sample. In contrast, samples B and C exhibited glass transition temperatures of 124 and 118 °C, respectively. Given that the maximum Tg for fully cured Akelite is approximately 120 °C, these results confirm that both samples B and C are effectively fully cured.
Conclusions
4
In this work, the polymerization and thermal degradation kinetics of Akelite resin were thoroughly analyzed, with particular attention to the role of the initiator in introducing an induction period that governs the overall reaction kinetics. Isothermal experiments enabled the characterization of the temperature dependence of this induction period, providing key insights into the incubation behavior associated with the Trommsdorff effect. The proposed methodology also allows predicting the temperature above which the polymerization reaction becomes self-accelerating under constant heating conditions. The proposed methodology for analyzing induction time offers a valuable framework for understanding the curing kinetics of resins influenced by this phenomenon.
In addition, thermal degradation kinetics were investigated through isoconversional analysis, yielding a detailed understanding of the resin’s thermal stability.
Based on these kinetic data, a Time–Temperature–Transformation (TTT) diagram was constructed to visualize the relationship between induction onset and degradation thresholds. This diagram serves as a practical guide for selecting suitable curing cycles, and the resulting processing map reveals a broad window within which composite panels can be manufactured safely and efficiently, minimizing the risk of thermal degradation and enhancing process reliability.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1European Commission Directorate-General for Communication, What is the European Green Deal?; Publications Office, 2019.
- 2Jensen P. D.Purnell P.Velenturf A. P. M.Highlighting the need to embed circular economy in low carbon infrastructure decommissioning: The case of offshore wind Sustain. Prod. Consum.20202426628010.1016/j.spc.2020.07.012 · doi ↗
- 3Verdejo Márquez, R. ; López Manchado, M. A. ; Santiago Bethencourt, A. Fibre-reinforced polymers based on thermoplastic matrices. 4,230,695 A 1, 2021.
- 4Carnicero R.Cano L.Cruz I.Lopez-Manchado M. A.Verdejo R.Analysis of a small recyclable wind turbine blade IOP Conf Ser. Mater. Sci. Eng.2023129301200810.1088/1757-899X/1293/1/012008 · doi ↗
- 5Vyazovkin S.Achilias D.Fernandez-Francos X.Galukhin A.Sbirrazzuoli N.ICTAC Kinetics Committee recommendations for analysis of thermal polymerization kinetics Thermochim. Acta 202271417924310.1016/j.tca.2022.179243 · doi ↗
- 6Vyazovkin S.Burnham A. K.Favergeon L.Koga N.Moukhina E.Pérez-Maqueda L. A.Sbirrazzuoli N.ICTAC Kinetics Committee recommendations for analysis of multi-step kinetics Thermochim. Acta 202068917859710.1016/j.tca.2020.178597 · doi ↗
- 7Vyazovkin S.Sbirrazzuoli N.Mechanism and Kinetics of Epoxy–Amine Cure Studied by Differential Scanning Calorimetry Macromolecules 1996291867187310.1021/ma 951162 w · doi ↗
- 8Vafayan M.Hosain Beheshty M.Ghoreishy M. H. R.Abedini H.Advanced integral isoconversional analysis for evaluating and predicting the kinetic parameters of the curing reaction of epoxy prepreg Thermochim. Acta 2013557374310.1016/j.tca.2013.01.035 · doi ↗
