高校化学工程学报    2020, Vol. 34 Issue (1): 222-229  DOI: 10.3969/j.issn.1003-9015.2020.01.028
0

引用本文 

张天萍, 甄卫军, 赵玲. 配方对有机硅橡胶非等温硫化动力学的影响[J]. 高校化学工程学报, 2020, 34(1): 222-229.   DOI: 10.3969/j.issn.1003-9015.2020.01.028.
ZHANG Tian-ping, ZHEN Wei-jun, ZHAO Ling. Effects of formulations of silicon rubber on its non-isothermal vulcanization kinetics[J]. Journal of Chemical Engineering of Chinese Universities, 2020, 34(1): 222-229.   DOI: 10.3969/j.issn.1003-9015.2020.01.028.

通讯联系人

赵玲, E-mail:zhaoling@ecust.edu.cn

作者简介

张天萍(1994-), 女, 新疆乌鲁木齐人, 新疆大学硕士生。

文章历史

收稿日期:2019-01-24;
修订日期:2019-03-31。
配方对有机硅橡胶非等温硫化动力学的影响
张天萍 , 甄卫军 , 赵玲     
新疆大学 化学化工学院,石油天然气精细化工教育部重点实验室,新疆 乌鲁木齐 830046
摘要:采用差示扫描量热法系统考察了硫化剂过氧化二异丙苯(DCP)用量、聚有机硅氧烷分子量、填料用量等对有机硅橡胶非等温硫化过程的影响。结果表明:此硫化过程表现出典型自催化反应特征;DCP的加入明显降低了交联起始温度,导致反应前期交联度过高,抑制了后期相对反应程度的增加;聚甲基乙烯基硅氧烷分子量增加也降低交联起始温度,但基本不影响相对反应程度的变化;由于气相白炭黑与硅橡胶基体形成结合胶,其用量对硫化过程影响是非单调的。运用Kissinger法得到含DCP体系的硫化表观活化能为138.77~149.65 kJ·mol-1;Šeatak-Berggren模型较好地拟合了其非等温硫化动力学,反应级数接近1级。
关键词甲基乙烯基硅橡胶    非等温    硫化    配方    反应动力学    自催化模型    
Effects of formulations of silicon rubber on its non-isothermal vulcanization kinetics
ZHANG Tian-ping , ZHEN Wei-jun , ZHAO Ling     
Key Laboratory of Oil & Gas Fine Chemicals of Ministry of Education, College of Chemistry and Chemical Engineering, Xinjiang University, Urumqi 830046, China
Abstract: Non-isothermal vulcanization processes of silicon rubber with different formulations were investigated via a differential scanning calorimeter method. Effects of vulcanizing agent dicumyl peroxide (DCP) dosage, molecular weight of silicon raw rubber and filler content on vulcanization kinetics were systematically analyzed. The results indicate that this vulcanization reaction shows typical autocatalysis features. The initial cross-linking temperature decreased with the addition of DCP, while the increase of relative conversion in the later stage was inhibited by high cross-linking degree in the beginning. Higher molecular weight of poly methyl vinyl silicone resulted in lower initial cross-linking temperature, which had little effects on relative conversion. Due to the formation of bound rubber between silica filler and silicone rubber matrix, the effect of filler content on vulcanization process was non-monotonic. For silicon rubbers with DCP in their formulations, the apparent activation energy of vulcanization reaction obtained by Kissinger method was 138.77~149.65 kJ·mol-1, and all kinetic curves can be well fitted by Šeatak-Berggren model with total reaction order of~1.
Key words: methyl-vinyl silicone rubber    non-isothermal    vulcanization    formulation composition    kinetics    Šeatak-Berggren autocatalysis model    
1 引言

性能优异的有机硅橡胶及其复合材料已广泛应用于工业、农业、交通运输、军事、建筑及日用品等领域[1]。有机硅橡胶一般是分子量为40~80万的高分子量聚有机硅氧烷(即生胶)加入补强填料和其它各种添加剂,采用有机过氧化物等为硫化剂在高温下生成,最常用的补强填料是白炭黑,它可使硅橡胶的强度增加数十倍。不同配方的有机硅橡胶加工性能不同,生胶分子量、硫化剂、填料的种类和用量等都对制备过程中发生的硫化反应产生影响[2-3]。聚合物固化反应及其动力学研究已有大量报道,它们大多集中在环氧树脂、不饱和聚酯和酚醛树脂、聚氨酯等体系[4-10],通常采用差示扫描量热(differential scanning calorimeter,DSC)的实验测量方法,并结合经典理论Kissinger方程[11]和Crane方程[12],基于n级模型及自催化模型建立反应动力学。对有机硅橡胶硫化过程,国内外学者通常更关注新型硅氧烷和新型硫化剂的制备,以及不同硅橡胶配方及相关硫化工艺对产品性能的影响[13-14],而对硫化动力学的研究大多针对特定组成含量的配方体系。HONG等[15]采用非等温DSC法对线性结构的聚二甲硅氧烷进行了硫化动力学分析,用不同方法计算了硅氢加成反应的表观活化能,发现自催化反应模型和修正的Chang模型的模拟曲线与实验值相吻合。VASILAKOS等[16]采用非等温DSC法对蒙脱土/乙烯基端聚硅氧烷纳米复合材料的硫化机理及相关动力学进行了分析,发现由于不同类型填料表面官能团的差异造成发生的物理或化学交联不同,从而影响反应进程。SONG等[17]采用扭矩法研究了不同类型的二氧化硅填料对甲基乙烯基硅橡胶硫化动力学的影响,发现有较高比表面积且含量较高的二氧化硅更有利于结合胶的形成,结合胶含量越高对硫化反应的抑制越明显,这是由于其在硫化过程中形成三维网状结构,阻隔了橡胶基体之间的联通性。总体上,这些有机硅橡胶硫化动力学的研究大多侧重于配方中单一组分变化的考察,而配方中主要成分生胶、硫化剂和填料的用量配比不同、生胶分子量高低均会影响硫化过程动力学,本文采用DSC法系统考察这些因素对硅橡胶非等温硫化动力学的影响,研究结果可以为不同配方硅橡胶加工参数的优化确定提供依据。

2 实验部分 2.1 实验原料

聚甲基乙烯基硅氧烷(110-2),生胶分子量分别为500 000、600 000和670 000 g·mol-1;羟基硅油,羟基含量6.0%~12.0%,黏度(295 K):≤ 20 mm2·s-1;气相白炭黑(HS-200),比表面积185~225 m2·g-1。以上原料均由合盛硅业股份有限公司提供,其中气相白炭黑粉体在使用前在75 ℃恒温鼓风干燥箱中预处理3 h去除水分,以便更好地分散;硫化剂过氧化二异丙苯(dicumyl peroxide,DCP),纯度99%,由上海阿拉丁生化科技股份有限公司提供。

2.2 实验过程

采用微量双螺杆挤出机(HAAKE MinLab2,德国哈克公司)制备了不同配方的混炼胶样品,螺杆转速为30 r·min-1,混炼5 min[18-19]。配方以甲基乙烯基硅橡胶100质量份(parts per hundreds of rubber, phr)和润滑剂羟基硅油4 phr为基准,硫化剂DCP分别选取0、1、2、3、4phr,填料气相白炭黑分别选取0、10、20、30 phr,不同配方组成见表 1

表 1 硅橡胶配方组成 Table 1 Formulations of silicone rubbers

随后,采用DSC(TA DHR-2,美国TA公司)分析了常压氮气氛围中不同配方混炼胶样品的非等温硫化过程,温度从350升至520 K,升温速率分别为5、10、15、20 K·min-1。不同硫化时间的相对反应程度α可由t时间内的放热量ΔH(t)与反应的总放热量ΔHtotal的比值确定。

$\alpha {\rm{ = }}\frac{{\Delta H\left( t \right)}}{{\Delta {H_{{\rm{total}}}}}}$ (1)
3 结果与讨论 3.1 硅橡胶硫化过程特征

有机硅橡胶在不同升温速率下进行硫化的典型DSC曲线见图 1,可以得到反应起始温度To,反应峰值温度TP和反应结束温度Tf 3个特征温度,将这些特征温度对不同升温速率β线性拟合,外推至β = 0时的ToTPTf分别为有机硅橡胶体系的交联起始温度、硫化温度和反应后处理温度[20-21]

图 1 常压N2氛围中有机硅橡胶硫化过程DSC曲线 Fig.1 DSC curves of silicone rubber vulcanization processes in N2 atmosphere

表 2为不同配方硅橡胶的特征温度汇总表,可以看出硫化剂DCP对交联初始温度的影响最为显著,不加DCP时交联初始温度很高为409 K,一旦加了DCP其下降幅度达到十多度,并随着DCP用量的增加进一步降低;生胶分子量越大,交联初始温度也有所下降;白炭黑含量对交联初始温度影响非单调。硫化温度和后处理温度受各因素影响的变化幅度相对较小。

表 2 不同升温速率下有机硅橡胶的ToTPTf Table 2 To, TP and Tf of silicone rubbers at different heating rates

图 2是升温速率为10 K·min-1时不同配方硅橡胶硫化相对反应程度随时间变化的关系图,可以发现在硅橡胶硫化过程中相对反应程度的上升均表现出了S型反应特征,即在硫化起始阶段变化缓慢,随后迅速增大,而接近相对反应程度100%时又变小,这是典型的自催化反应特征,其原因可归结为有机硅橡胶发生过氧化物交联反应[22]。在反应的初期,DCP受热分解产生自由基,反应处于诱导期,交联反应缓慢;自由基上的未配对P电子活性很大,硅氧键键长较长,键对侧基转动的位阻小,很容易受到自由基的进攻,其开始催化交联反应,引起连锁反应加快了反应速率;而在反应末期,由于组分的减少,反应速率又减慢。

图 2 10 K·min-1升温速率下硅橡胶硫化过程相对反应程度的变化 Fig.2 Profiles of relative conversion of silicon rubber vulcanization under heating rate of 10 K·min-1

此外,从图 2(a)可知,是否加硫化剂DCP对相对反应程度变化影响显著,随着硫化剂用量增加,相对反应程度上升慢,这是由于硫化剂用量的增大,交联在较低的初始温度即可发生,从而导致反应前期交联度过高,反应活性组分减少且传质阻力增大,反而抑制了后期交联反应的发生。与此同时随着DCP用量的增加硫化时间增长,该现象可以归结为聚有机硅氧烷的交联反应是以硫化剂DCP的分解为控制步骤[22]。因此只要能达到所需的交联度,应尽量减少硫化剂的用量。由图 2(b)可知,聚有机硅氧烷生胶分子量的大小基本不影响硫化反应过程中相对反应程度的变化,这可能是由于交联初始温度随生胶分子量增加有所下降,而传质阻力随生胶分子量增加又有所增加,两者相抵消导致影响程度小,文献也曾报道同一种链结构的高分子其交联效率不随分子量的大小而变化[22]。从图 2(c)可知,气相白炭黑填料用量对反应程度变化影响是非单调的,即相比较无填料体系,填料含量为10和20份时前期相对反应程度上升减慢,而当填料含量为30份时,前期相对反应程度变化又加快,这是因为硅生胶与填料混合后有相当一部分被SiO2表面吸附形成结合胶,影响了有机硅化合物的交联反应;当填料超过一定量时,填料及结合胶将生胶阻隔开限制了反应的进一步进行[17, 23]

Kissinger法已经广泛地应用于确定聚合物交联反应过程中的表观活化能(Ea)[20, 24]

$\ln \left( {\frac{\beta }{{T_{\rm{P}}^2}}} \right) = \ln \left( {\frac{{AR}}{{{E_{\rm{a}}}}}} \right) - \frac{{{E_{\rm{a}}}}}{{R{T_{\rm{P}}}}}$ (2)

其中TP为峰值温度(K),Ea为表观硫化反应活化能(J·mol-1),A为指前因子,R为气体常数,β为升温速率(K·min-1)。由ln(β/TP2)对1/Tp作线性回归方程即可得出EaA的值。表 3为计算得到的不同配方硅橡胶的EaA。不加DCP的表观活化能低是由于交联反应在相对高的温度下发生。一旦加了DCP,随着有机硅橡胶分子量的增大和填料含量的增加,Ea增加幅度并不大。但其随着DCP含量变化较明显,这是由于过氧化物产生自由基先加到乙烯基上引起交联,乙烯基用完后,剩余的自由基需克服更大的能垒在甲基上进一步交联[23]

表 3 不同配方有机硅橡胶硫化反应的EaA Table 3 Ea and A values of vulcanization of silicon rubbers with different formulations
3.2 机理函数及硫化动力学参数的确定

有机硅橡胶硫化反应表现出了明显的自催化反应特征,通常自催化反应动力学模型(3)中反应级数mn可以通过Málek法确定[25],该方法引入了两个特殊形状因子函数y(α)和z(α)确定Šeatak-Berggren自催化模型参数[26]

$\frac{{{\rm{d}}\alpha }}{{{\rm{d}}t}} = A{{\rm{e}}^{ - x}}{\alpha ^m}{\left( {1 - \alpha } \right)^n}$ (3)
$y\left( \alpha \right) = \left( {\frac{{{\rm{d}}\alpha }}{{{\rm{d}}t}}} \right){e^x}$ (4)
$z\left( \alpha \right) = {\rm{\pi}} \left( x \right)\left( {\frac{{{\rm{d}}\alpha }}{{{\rm{d}}t}}} \right)\frac{T}{\beta }$ (5)

式中,x = Ea/RT,表观活化能Ea已由Kissinger法求得;π(x)为温度的积分形式,是温度的函数,可用SENUM和YANG提出的方程近似表达[27],即

$\pi \left( x \right) = \frac{{{x^3} + 18{x^2} + 88x + 96}}{{{x^4} + 20{x^3} + 120{x^2} + 240x + 120}}$ (6)

αMy(α)对应于最大值的交联度,$\alpha _{\rm{P}}^\infty $z(α)对应于最大值的交联度。不同配方硅橡胶硫化反应的y(α)~ αz(α)~α关系式中的详细参数如表 4所示。

表 4 不同配方有机硅橡胶y(α)和z(α)的曲线峰值αMαP Table 4 Peak curve values of αM and αP of silicone rubbers with different formulations

通过y(α)和z(α)的形状及其特征值αM$\alpha _{\rm{P}}^\infty $可以推断动力学机理函数,由Málek提出的判据,如果0 < αM < $\alpha _{\rm{P}}^\infty $,并且$\alpha _{\rm{P}}^\infty \ne 0.632$,则可确定有机硅橡胶体系非等温硫化反应动力学机理函数符合自催化模型,可将动力学方程(3)变形为

$\ln \left[ {\left( {\frac{{{\rm{d}}\alpha }}{{{\rm{d}}t}}} \right)\exp \left( {\frac{{{E_{\rm{a}}}}}{{RT}}} \right)} \right] = \ln A + n\ln \left[ {{\alpha ^{{\alpha _{\rm{M}}}/\left( {1 - {\alpha _M}} \right)}}\left( {1 - \alpha } \right)} \right]$ (7)

式中$p=\alpha_{\mathrm{M}} /\left(1-\alpha_{\mathrm{M}}\right)=m / n$,以$\ln \left[(\mathrm{d} \alpha / \mathrm{d} t) \exp \left(E_{\mathrm{a}} / R T\right)\right]$$\ln \left[ {{\alpha ^p}\left( {1 - \alpha } \right)} \right]$作线性回归,可得到斜率,计算出nmA的值列于表 5。由表可知,加了DCP后该反应总级数基本接近为1,与文献报道一致[28]

表 5 不同配方有机硅橡胶自催化反应动力学参数 Table 5 Kinetic parameters of autocatalytic reaction of silicone rubbers with different formulation

图 3~5为不同配方的硅橡胶硫化动力学实验数据与自催化模型计算数据的对比。加了DCP后,自催化模型可以较好地拟合实验数据,反应初期出现了一些偏差,是由于随着升温速率的提高,放热的累积效应导致热惯性大[29]。而未添加DCP时,见图 3(a),后期拟合效果较差,这是由于未添加硫化剂时热惯性的延迟导致。

图 3 不同DCP量有机硅橡胶硫化动力学实验数据与自催化模型拟合结果对比 Fig.3 Comparison of vulcanization kinetic experimental data with model predicated values for silicon rubber with different DCP dosages
图 4 不同分子量有机硅橡胶硫化动力学实验数据与自催化模型拟合结果对比 Fig.4 Comparison of vulcanization kinetic experimental data with model predicated values for silicon raw rubber with different molecular weights
图 5 不同气相白炭黑添加量的有机硅橡胶硫化动力学实验数据与自催化模型拟合结果对比 Fig.5 Comparison of vulcanization kinetic experimental data with model predicated values for silicon rubber with different silica filler contents
4 结论

以甲基乙烯基硅橡胶100 phr为基准,在生胶分子量500 000~670 000 g·mol-1、硫化剂DCP用量0~4 phr和填料气相白炭黑用量0~30 phr范围内,采用DSC法研究了不同配方组成有机硅橡胶的非等温硫化过程,结果表明:

(1) 不同配方有机硅橡胶硫化过程相对反应程度的变化趋势均呈现出典型的自催化反应特征。DCP用量对硫化过程影响最为显著,不加DCP的交联起始温度高,随着DCP用量的增加,较低的交联起始温度导致反应前期交联度过高反而抑制了后期相对反应程度的上升,硫化时间明显延长;生胶分子量大小对硫化过程相对反应程度变化的影响较小;气相白炭黑含量对硫化过程则表现出了非单调的影响。

(2) 采用Kissinger方程法计算得到了表观活化能Ea。未添加DCP时高温硫化的Ea为113.26 kJ·mol-1;添加了DCP和填料的体系,随着生胶分子量的增加、DCP含量和填料含量的增加,Ea有所增加,在138.77~149.65 kJ·mol-1

(3) 自催化Šeatak-Berggren模型可以较好地拟合不同配方有机硅橡胶的非等温硫化动力学。通过Málek最概然机理函数法确定了动力学参数,添加了DCP配方体系的反应级数均接近1级。

符号说明:

A  —指前因子    T  —温度,K

dα/dt  —反应速率,min-1    To  —反应起始温度,K (β = 0时,交联起始温度)

Ea  —表观活化能,J·mol-1    Tp  —反应峰值温度,K (β = 0时,硫化温度)

ΔHt  —时间内的放热量,J·g-1    Tf  —反应结束温度,K (β = 0时,反应后处理温度)

ΔHtotal  —非等温硫化过程总放热量,J·g-1    t  —硫化反应时间,min

m, n  —反应级数    α  —相对反应程度

R  —气体常数,8.3145 J·mol-1·K-1    β  —升温速率,K·min-1

参考文献
[1]
范召东, 孙全吉, 黄艳华, 等. 有机硅材料在航空工业的应用[J]. 有机硅材料, 2015, 29(6): 491-498.
FAN Z D, SUN Q J, HUANG Y H, et al. Application of silicone materials in aviation industry[J]. Silicone Material, 2015, 29(6): 491-498.
[2]
ZHAO S, XIE S, LUI X, et al. Covalent hybrid of graphene and silicon dioxide and reinforcing effect in rubber composites[J]. Journal of Polymer Research, 2018, 25(225): 1-10.
[3]
UTARA S, JANTACHUM P, SUKKANEEWAT B. Effect of surface modification of silicon carbide nanoparticles on the properties of nanocomposites based on epoxidized natural rubber/natural rubber blends[J]. Journal of Applied Polymer Science, 2017, 134(45289): 1-14.
[4]
LYU J X, HU D D, LIU T, et al. Non-isothermal kinetics of epoxy resin curing reaction under compressed CO2[J]. Journal of Thermal Analysis and Calorimetry, 2018, 131(2): 1499-1507.
[5]
HU D D, LYU J X, LIU T, et al. Solvation effect of CO2 on accelerating the curing reaction process of epoxy resin[J]. Chemical Engineering and Processing-Process Intensification, 2018, 127: 159-167.
[6]
YU J, YE H, LIU Y. Synthesis of spirocyclosiloxanes for transparent copolymer thermosets[J]. Journal of Applied Polymer Science, 2018, 135(46370): 1-6.
[7]
杨泽, 胡冬冬, 刘涛, 等. 高压气体氛围中的聚氨酯非等温固化动力学[J]. 化工学报, 2018, 69(11): 4728-4736.
YANG Z, HU D D, LIU T, et al. Non-isothermal curing kinetics of polyurethane under high-pressure gas atmosphere[J]. CIESC Journal, 2018, 69(11): 4728-4736.
[8]
董红星, 赵庆晨, 张进, 等. CMC/E44/DDM体系固化过程及动力学[J]. 化工学报, 2014, 65(11): 4392-4399.
DONG H X, ZHAO Q C, ZHANG J, et al. Curing process and kinetics of CMC/E44/DDM systems[J]. CIESC Journal, 2014, 65(11): 4392-4399.
[9]
严芳芳, 方波, 卢拥军, 等. 有机锆交联羟丙基瓜胶流变动力学[J]. 高校化学工程学报, 2015, 29(4): 802-808.
YAN F F, FANG B, LU Y J, et al. Rheokinetics of organic zirconium crosslinked hydroxypropyl guar gels[J]. Journal of Chemical Engineering of Chinese Universities, 2015, 29(4): 802-808.
[10]
HONG C, WANG X, PAN Z, et al. Curing thermodynamics and kinetics of unsaturated polyester resin with different chain length of saturated aliphatic binary carboxylic acid[J]. Journal of Thermal Analysis and Calorimetry, 2015, 122(1): 427-436.
[11]
KISSINGER H E. Reaction kinetics in differential thermal analysis[J]. Analytical Chemistry, 1957, 29(11): 1702-1706.
[12]
CRANE L W, DYNES P J, KAELBLE D H. Analysis of curing kinetics in polymer composites[J]. Journal of Polymer Science Part Polymer Letters Edition, 1973, 11(8): 533-540.
[13]
NOURI N M, SAADATBAKHSH M. Fabrication method of large-scale and mechanically durable superhydrophobic silicon rubber/aerogel coating on fibrous substrates[J]. Journal of Coatings Technology and Research, 2017, 14(2): 477-488.
[14]
SONG Y H, ZENG L B, ZHENG Q. Understanding the reinforcement and dissipation of natural rubber compounds filled with hybrid filler composed of carbon black and silica[J]. Chinese Journal of Polymer Science, 2017, 35(11): 1436-1446.
[15]
HONG I K, LEE S. Cure kinetics and modeling the reaction of silicone rubber[J]. Journal of Industrial and Engineering Chemistry, 2013, 19(1): 42-47.
[16]
VASILAKOS S P, TARANTILI P A. In situ monitoring by DSC and modeling of curing of vinyl polysiloxanes in layered silicate nanocomposites[J]. Journal of Thermal Analysis and Calorimetry, 2017, 127(3): 2049-2058.
[17]
SONG L X, LI Z H, CHEN L, et al. The effect of bound rubber on vulcanization kinetics in silica filled silicone rubber[J]. RSC Advances, 2016, 6(103): 101470-101476.
[18]
BAI J W, LIAO X, HUANG E B, et al. Control of the cell structure of microcellular silicone rubber/nanographite foam for enhanced mechanical performance[J]. Materials and Design, 2017, 133: 288-298.
[19]
YANG Q, YU H T, SONG L X, et al. Solid-state microcellular high temperature vulcanized (HTV) silicone rubber foam with carbon dioxide[J]. Journal of Applied Polymer Science, 2017, 134(44807): 1-9.
[20]
房佳丰, 郑铭焕, 覃柳琪, 等. 加成型液体硅橡胶硫化动力学的非等温DSC研究[J]. 有机硅材料, 2018, 32(1): 1-7.
FANG J F, ZHENG M H, QIN L Q, et al. Study on non-isothermal DSC of curing kinetics of additional liquid silicone rubber[J]. Silicone Material, 2018, 32(1): 1-7.
[21]
牛芊, 周彩荣, 詹自力. 甲硫氨酸的非等温热分解动力学的研究[J]. 高校化学工程学报, 2018, 32(5): 1112-1118.
NIU Q, ZHOU C R, ZHAN Z L. Study on non-isothermal decomposition kinetics of methionine[J]. Journal of Chemical Engineering of Chinese Universities, 2018, 32(5): 1112-1118.
[22]
来国桥. 有机硅化学与工艺[M]. 北京: 化学工业出版社, 2011.
LAI G Q. Organic silicon chemistry and technology[M]. Beijing: Chemical Industry Press, 2011.
[23]
李光亮. 有机硅高分子化学[M]. 北京: 科学出版社, 1998.
LI G L. Organic silicon polymer chemistry[M]. Beijing: Science Press, 1998.
[24]
LOPEZ L M, COSGROVE A B, HERNANDEZORTIZ J P, et al. Modeling the vulcanization reaction of silicone rubber[J]. Polymer Engineering and Science, 2007, 47(5): 675-683.
[25]
JIŘÍ M. A computer program for kinetic analysis of non-isothermal thermoanalytical data[J]. Thermochimica Acta, 1989, 138(2): 337-346.
[26]
JAROSLAV Š, BERGGREN G. Study of the kinetics of the mechanism of solid-state reactions at increasing temperatures[J]. Thermochimica Acta, 1971, 3(1): 1-12. DOI:10.1016/0040-6031(71)85051-7
[27]
SENUM G I, YANG R T. Rational approximations of the integral of the Arrhenius function[J]. Journal of Thermal Analysis, 1977, 11(3): 445-447. DOI:10.1007/BF01903696
[28]
刘慧娟, 展喜兵, 祝文亲, 等. 液体乙烯基硅树脂的制备及固化反应动力学[J]. 化工学报, 2014, 65(4): 1509-1516.
LIU H J, ZHAN X B, ZHU W Q, et al. Preparation and curing reaction kinetics of liquid vinyl silicone resin[J]. CIESC Journal, 2014, 65(4): 1509-1516.
[29]
宋秋生, 徐蕤, 杨洋, 等. BPFER/DSX体系固化反应动力学与热性能[J]. 固体火箭技术, 2012, 35(2): 262-266.
SONG Q S, XU R, YANG Y, et al. Curing kinetics and thermal properties of BPFER/DSX system[J]. Journal of Solid Rocket Technology, 2012, 35(2): 262-266.