《Polymer Engineering & Science》：Extension of Isoconversional Analysis to Diffusion-Controlled Epoxy Cure Kinetics: Comparison With Model-Fitting Approaches
编辑推荐：
摘要翻译：环氧树脂的固化动力学通常采用模型拟合或无模型（等转化）方法描述。在模型拟合方法中，Kamal–Sourour模型被广泛接受，但需要调整以描述高转化率下向扩散控制区域的转变。等转化方法虽然更灵活，但传统上被认为无法捕捉扩散控制下的固化演化过程。在本研究
摘要翻译：环氧树脂的固化动力学通常采用模型拟合或无模型（等转化）方法描述。在模型拟合方法中，Kamal–Sourour模型被广泛接受，但需要调整以描述高转化率下向扩散控制区域的转变。等转化方法虽然更灵活，但传统上被认为无法捕捉扩散控制下的固化演化过程。在本研究中，研究人员将等转化分析扩展至扩散控制区域，并平行修改了Kamal–Sourour模型，引入扩散项以进行一致性比较。研究人员首先证明等转化方法能够再现Kamal–Sourour模型衍生的所有现象学特征。随后将两种方法应用于商用RTM6环氧基树脂的实验固化数据，结果表明等转化分析在捕捉复杂实验现象方面具有更优能力。此外，所提出的等转化分析扩展即使在扩散主导固化速率时也能提供准确预测。研究结果表明，一旦扩展至考虑扩散控制，等转化方法在再现实验固化现象方面优于模型拟合方法。可靠的动力学模型对于优化工业制造过程中的固化周期至关重要；所提出的方法为开发环氧基复合材料结构的此类模型提供了准确且灵活的框架。
论文解读文章
**研究背景与问题**
环氧树脂因其优异的力学、热学和化学性能，在现代工程中广泛应用，尤其是在航空航天领域。准确表征固化动力学对于模拟和优化制造工艺至关重要，能够预测固化反应随温度的演化，支持固化周期模拟和加工图谱开发。然而，环氧树脂的固化过程涉及复杂反应路径，包括伯胺-环氧加成、仲胺-环氧加成以及醚化等副反应，且随着固化进行，体系由化学动力学控制转变为扩散控制。传统模型拟合方法中，Kamal–Sourour（K–S）模型被广泛采用，但该模型无法描述高转化率下扩散控制的影响，需要引入经验修正。等转化方法（如Friedman微分法）虽更具灵活性，但长期以来被认为无法捕捉扩散控制阶段的固化动力学。因此，本研究旨在扩展等转化分析方法至扩散控制区域，并与改进的K–S模型进行系统比较，以验证其适用性和优越性。
**研究内容与结论**
研究人员首先通过数值模拟评估Friedman等转化方法对理想K–S模型（含扩散项）的再现能力。模拟采用RTM6树脂的文献动力学参数，在无扩散条件下，等转化分析能精确拟合K–S模型，活化能随转化率从约73 kJ/mol降至58 kJ/mol，对应两个竞争反应机制。引入扩散后，低升温速率（≤1.25°C/min）下高转化率阶段（α>0.8）出现明显偏离，R2下降，表明传统等转化方法失效。为此，研究人员提出两步法：首先利用高升温速率（β≥2.5°C/min）实验数据通过等转化分析确定纯化学动力学，然后基于阻力串联模型（1/k = 1/k_Iso + 1/k_D）拟合扩散速率常数k_D（采用Macedo–Litovitz模型，结合Arrhenius项和自由体积项）。数值验证表明，拟合的扩散活化能（137.4 kJ/mol）与预设值（138 kJ/mol）高度吻合，模型能准确预测等温和非等温条件下的固化行为。
随后，研究人员对商用RTM6树脂进行实验验证。差示扫描量热（DSC）实验在0.6–20°C/min升温速率和等温条件（145–170°C）下进行，并测量不同转化率下的玻璃化转变温度以确定DiBenedetto参数。等转化分析显示，在α<0.8区域拟合极佳，但α>0.8时活化能上升且R2下降，实验数据在低升温速率0.6°C/min下出现明显扩散影响。仅采用β≥2.5°C/min数据进行等转化拟合，所得模型能准确预测所有升温速率下的固化曲线（包括未参与拟合的0.6和1.25°C/min），仅在α>0.9时低升温速率下偏离归因于扩散。与之对比，K–S模型拟合（CMA-ES算法优化）无法捕捉实验数据末尾的肩峰，表明其不能描述醚化等副反应。研究人员进一步指出，实验活化能在高转化区仍上升，且所有峰在末尾出现肩峰，提示存在第三种反应机制（醚化），而K–S模型无法涵盖。通过等温实验（145–170°C）拟合扩散项，获得扩散活化能143 kJ/mol，与文献值一致。最终模型预测的等温和非等温固化曲线与实验数据高度吻合，TTT图预测的34个等转化时间点的平均相对误差仅2.4%，标准差3.6%。
**关键技术方法**
本研究采用的主要技术方法包括：Friedman微分等转化分析（基于方程dα/dt = A exp(-Eα/RT) f(α)），用于提取活化能和指数前因子；差分扫描量热法（DSC）测量固化放热和玻璃化转变温度（Tg）；DiBenedetto方程拟合Tg-α关系；Macedo–Litovitz扩散模型（k_D = A_D exp(-E_D/RT) exp(-b/f)）描述扩散速率常数；阻力串联模型（1/k = 1/k_Iso + 1/k_D）整合化学与扩散贡献；协方差矩阵自适应进化策略（CMA-ES）算法优化模型参数。样本为Hexcel公司的RTM6商用环氧树脂，实验在TA Instruments Q2000 DSC上进行，氮气保护，升温速率范围0.6–20°C/min，等温实验在140–180°C进行，Tg测量按照ASTM E1356标准。
**研究结果**
**1. 数值模拟：无扩散情况**
通过模拟K–S模型（不含扩散）的六种升温速率数据，Friedman分析获得R2≈1，等转化预测与模拟数据完全一致，活化能从73 kJ/mol降至58 kJ/mol，证明等转化方法能准确描述多步反应动力学。
**2. 数值模拟：含扩散情况**
引入扩散项后，低升温速率（0.6和1.25°C/min）下α>0.8时等转化分析R2下降，预测等温固化时高转化率段速率偏离，表明传统方法失效。通过排除低升温速率数据，可分离纯化学动力学，验证了两步法的可行性。
**3. 扩散修正方法验证**
采用阻力串联模型，以高升温速率数据确定k_Iso，以等温模拟数据拟合k_D，获得扩散活化能137.4 kJ/mol（预设138 kJ/mol），R2=0.9998，模型预测与模拟数据在等温和非等温条件下均高度一致。
**4. 实验数据分析：化学动力学**
实验DSC数据表明，低升温速率（0.6°C/min）末尾出现明显扩散效应，而高升温速率（≥2.5°C/min）数据无扩散影响。等转化分析（β≥2.5°C/min）成功拟合化学动力学，并能预测未参与拟合的升温速率曲线，仅α>0.9时0.6°C/min数据出现偏差。K–S模型拟合虽能描述整体趋势，但无法捕捉末尾肩峰，R2低于等转化方法，表明其不能反映醚化等复杂机制。
**5. 实验数据分析：扩散拟合与验证**
利用等温实验（145–170°C）拟合扩散项，获得AD=1.96×10^20 s^-1，ED=143 kJ/mol，R2=0.9831。模型预测与等温（140°C、180°C）和非等温实验均吻合，TTT图预测误差极小，验证了模型的普适性。
**讨论与结论**
讨论部分强调，传统等转化方法在扩散控制下失效，而所提出的两步法通过分离化学和扩散贡献，有效解决了这一局限。实验结果表明，等转化框架能捕捉醚化等副反应，而K–S模型无法实现，体现了等转化方法的灵活性。最终结论指出，扩散修正的等转化分析为环氧固化建模提供了稳健且广泛适用的框架，无需对反应机制进行先验假设，即可实现可靠预测，对于工业固化周期优化具有重要意义。该研究发表在《Polymer Engineering》期刊上。
