季戊四醇丙烯醛树脂是由一丙烯醛季戊四醇和二丙烯醛季戊四醇组成的一种透明黏稠状液体树脂。该树脂的固化产物具有粘结强度高, 固化收缩率低, 耐酸碱和化学稳定性好等优点, 其最重要的应用领域是固体推进剂和改性炸药。以季戊四醇丙烯醛树脂为粘结剂的固体推进剂燃速与燃温高, 强度大, 在点火时遇到高温和冲击不龟裂, 可长时间承受-65 ℃的低温和135 ℃的高温, 该树脂还能够改善端羟基聚丁二烯的力学性能[1-2]。
尽管国内外对季戊四醇丙烯醛树脂的合成及应用进行了许多研究[3-5], 然而, 针对该树脂体系的固化反应动力学及流变特性的研究却鲜有报道。但其固化反应动力学和流变特性对于了解树脂本身的特性以及固化工艺优化都有重要价值。目前, 研究树脂固化动力学的方法很多, 包括经典的傅里叶变换红外光谱(FTIR)、热重分析法、动态差示扫描量热法(DSC)、萃取法以及动态扭振法等, 其中动态DSC法是最常用的固化动力学研究方法, 采用动态DSC法可研究固化反应过程和固化反应机理, 并可计算出固化反应的动力学参数。近年, 一些学者开展了利用动态DSC结果确定树脂体系的固化工艺研究, 取得了较好的结果[6-7]。
本研究利用动态DSC法对季戊四醇丙烯醛树脂的非等温固化过程进行了研究, 考察了升温速率、固化时间对固化度的影响; 采用Kissinger, Crane方程[8-10]建立了该树脂体系的固化动力学模型, 并对该模型与实验结果的一致性进行了验证。同时研究了季戊四醇丙烯醛树脂固化过程的流变特性, 获得了粘度模型, 以期对季戊四醇丙烯醛树脂体系固化工艺的确定提供依据。
2 实验部分 2.1 实验试剂季戊四醇丙烯醛树脂, 黎明化工研究院, 平均摩尔质量500 g·mol-1, 密度1.2 g·mL-1, 粘度20 mPa·s (50 ℃); 硫酸二乙酯(DES), 分析纯, 成都联合化工试剂研究所。
2.2 测试仪器及实验方法差示扫描量热仪: PE Diamond, 功率补偿型。DV-Ⅲ+型旋转粘度计: Brookfield公司, 粘度测量范围为1 ~100000 mPa·s, 测试标准为GB/T 7193-2008。恒温水浴: Thermo公司, 控温精度0.01 K。实验过程中所有的试剂称量都使用Mettler-Toledo梅特勒公司生产的New Classic MS天平, 其精确度为0.0001 g。
动态DSC曲线的测定:样品质量10 mg左右, 氮气气氛, 气流速度为30 mL·min-1; 测定季戊四醇丙烯醛树脂体系在升温速率为10, 15, 20, 30 ℃·min-1下的固化行为, 参比样为α-Al2O3, 铝坩埚。
恒温粘度的测量:将季戊四醇丙烯醛树脂升温至实验测定温度, 然后加入催化剂(质量为季戊四醇丙烯醛树脂的1/50), 在恒温水浴中将体系混合均匀。将样品池加热到所设定的温度, 并开动转子使其在空气中归零, 然后将混合均匀的样品迅速加入到样品池中, 测定测试温度为68, 74, 80, 86 ℃下不同时间树脂体系的粘度。
3 实验结果与分析 3.1 季戊四醇丙烯醛树脂体系热分析采用动态DSC法测定了季戊四醇丙烯醛树脂体系在固化过程中的热效应, 建了固化动力学模型。季戊四醇丙烯醛树脂体系在4种升温速率下的固化放热DSC曲线见图 1。从图 1可以得出不同升温速率下体系的特征固化温度, 见表 1。图 1和表 1表明, 随着升温速率(β)的增大, 反应起始温度(Ti)、峰顶温度(Tp)和反应终止温度(Te)都会升高, 反应热(Q)降低。
对表 1中的反应热与升温速率进行线性拟合, 外推至升温速率为零, 可以得到体系完全固化时的反应热为: 149.01 J·g-1。由于热固性树脂的固化度只与反应热相关, 因此固化度可定义为[7]:
$ \alpha = {H_{\rm{t}}}/{H_{\rm{u}}} $ | (1) |
式中, Ht为反应热, J·g-1; Hu为树脂完全固化的总反应热, J·g-1。由式(1)可以推断出体系在不同升温速率下的固化度, 结果如表 2。
由表 2可以看出, 升温速率越快, 其固化程度越低。这是因为在固化反应过程中, 已经形成的网络结构会限制未反应官能团的自由运动, 减小其有效碰撞的几率, 升温速率越快, 固化反应越快, 已反应的部分对未反应的官能团的限制就越明显, 固化程度就越低。
3.2 固化动力学模型参数计算及动力学方程的验证Kissinger微分法是处理动力学DSC曲线数据, 获取各表观动力学参数的常用方法[8-9]。假设季戊四醇丙烯醛树脂固化符合n阶固化动力学模型, 则Kissinger方程可以描述为:
$ \ln \left( {\beta /T_{\rm{p}}^2} \right) = \ln \left( {AR/{E_{\rm{a}}}} \right)-{E_{\rm{a}}}/\left( {R{T_{\rm{p}}}} \right) $ | (2) |
式中, β为升温速率常数, ℃·min-1; R为理想气体常数, 8.314 J·mol-1·K-1; Tp为反应放热峰温度, K; A为指前因子, K; Ea为反应活化能, kJ·mol-1。
以ln(β/Tp2)对1/Tp作图, 并经线性拟合, 如图 2所示。由斜率可求得反应活化能: Ea=14.32 kJ·mol-1, 由截距可得指前因子: A=3.84×103 K。
固化反应级数n采用Crane经验方程[10]进行估算:
$ \ln \beta =-{E_{\rm{a}}}/\left( {nR{T_{\rm{p}}}} \right) + C $ | (3) |
式中, n是反应级数。以-lnβ对1/Tp作图, 线性回归求得固化反应级数n=0.95。
根据Kamal模型, 固化反应动力学方程可表示为[11]:
$ {\rm{d}}\alpha /{\rm{d}}t = k\left( T \right)f\left( \alpha \right) $ | (4) |
式中, α是固化度, 由式(1)计算; k(T)是固化速度常数, 通常用Arrhenius公式表示[12]:
$ k\left( T \right) = A\exp \left( {-{E_{\rm{a}}}/RT} \right) $ | (5) |
n阶固化模型表达式为:
$ {\rm{d}}\alpha /{\rm{d}}t = A\exp \left( {-{E_{\rm{a}}}/RT} \right){\left( {1-\alpha } \right)^{\rm{n}}} $ | (6) |
将固化反应动力学参数代入式(6)中得到季戊四醇丙烯醛树脂体系固化反应动力学方程:
$ {\rm{d}}\alpha /{\rm{d}}t = 3.84 \times {10^3}\exp \left( {-14.32 \times {{10}^3}/RT} \right){\left( {1-\alpha } \right)^{0.95}} $ | (7) |
为了验证式(7)的有效性, 测试了季戊四醇丙烯醛树脂体系在359 K下的恒温DSC曲线, 获得了不同反应时间的固化度。图 3给出了359 K下的实验结果与反应动力学方程模拟曲线的对比。由图 3可以看出, 式(7)的理论模拟曲线与实验结果比较吻合, 相关系数为0.9724, 说明季戊四醇丙烯醛树脂体系固化反应动力学方程是适用的。
季戊四醇丙烯醛树脂性能与其固化度密切相关, 通常固化度越高, 热性能越好, 而固化反应进行的程度受固化工艺影响很大, 所以, 固化反应的工艺过程对最终的固化产物的性能起决定性作用。图 4给出了由式(7)计算得到的不同固化温度下固化度与反应时间的关系。由图 4可见, 固化反应温度越低, 要到达某一固化度所用时间越长, 但最终产品所能达到的固化度越高。
粘度是衡量固化体系主要工艺参数之一。研究了341, 347, 353, 359 K温度下, 季戊四醇丙烯醛树脂体系粘度随时间的变化规律, 结果如图 5所示。由图 5可见, 反应温度越高, 体系的初始粘度越低, 但粘度增加越快, 凝胶点出现的时间越短。这是因为当温度升高时, 大分子链段活动能力增加, 流动性增大, 即粘度减小, 但温度继续升高, 聚合交联反应加快, 则树脂流动性随温度升高而降低, 粘度增大。
双阿累尼乌斯化学流变模型(Dual Arrhenius)是在研究热固性树脂体系粘度变化规律的基础上提出的半经验公式, 研究证明它适用于总体为一级反应和多数非一级动力学反应的热固性树脂体系[13-14]。本研究采用Dual Arrhenius模型研究了季戊四醇丙烯醛树脂在整个工艺温度范围内的粘度。该流变模型如下:
$ \eta = {\eta _0}\left( T \right) + \Delta {\eta _{\rm{C}}} $ | (8) |
式中, η0 (T)表示恒温温度T时体系的初始粘度, mPa·s; ηc表示化学变化对粘度的贡献, mPa·s。
$ {\eta _0}\left( T \right) = {\eta _\infty }\exp \left( {{k_1}/RT} \right) $ | (9) |
$ \Delta {\eta _C} = \exp \left( {{k_2}t} \right) $ | (10) |
式中, η∞为温度无限高时体系的理想粘度, mPa·s; k1, k2为树脂的化学流变参数, J·mol-1; R为气体常数, 8.314 J·mol-1·K-1; T为温度, K; t为反应时间, min。
图 6给出了季戊四醇丙烯醛树脂体系初始粘度随温度变化的曲线。由图 6可知, 两者线性关系较好, 相关系数为0.9986。并采用方程(9)对温度与体系初始粘度数据进行拟合, 得出此时η∞和k1分别为: 5.30×10-6 mPa·s和4.71×104 J·mol-1。结合式(9)和式(10)通过最小二乘法拟合得到了季戊四醇丙烯醛树脂体系在不同温度时的k2, 结果如表 3。
令k2 (T) =a+bT, 对表 3的数据中温度和k2值进行线性拟合, 线性相关系数为0.9886, 表明两者线性关系良好, 得到的线性拟和方程为:
$ {k_2} =-1.152 + 3.44 \times {10^{-3}}T $ | (11) |
将式(9)和式(10)带入式(8)可得季戊四醇丙烯醛树脂反应体系粘度计算式为:
$ \eta = {\eta _\infty }\exp \left( {{k_1}/RT} \right) + \exp \left( {{k_2}t} \right) $ | (12) |
将η∞, k1和k2带入(12)中可得体系在不同温度条件下的粘度计算式为:
$ \begin{array}{l} \eta = 5.30 \times {10^{-6}}\exp \left( {4.71 \times {{10}^4}/RT} \right) + \\ \;\;\;\;\;\;\exp \left( {\left( {-1.152 + 3.44 \times {{10}^{-3}}T} \right)t} \right) \end{array} $ | (13) |
为验证粘度方程(13), 分别将方程(13)的计算粘度值和实验结果进行了对比。图 7给出了在350 K时季戊四醇丙烯醛树脂体系粘度随时间变化的模拟曲线以及实验结果, 结果显示两者吻合较好, 相关系数为0.9782, 说明方程(13)能较好地预测季戊四醇丙烯醛树脂体系在不同温度下粘度随时间的变化。
(1) 采用动态DSC法研究了季戊四醇丙烯醛树脂体系的固化过程, 利用Kissinger方程和Crane方程计算了固化反应动力学参数, 反应活化能Ea=14.32 kJ·mol-1, 指前因子A=3.84×103 K, 反应级数n=0.95。
(2) 建立了固化反应动力学模型: dα/dt=3.84×103exp(-14.32×103/RT)(1-α)0.95, 利用该模型对固化反应温度、固化度与时间的关系进行了理论分析。
(3) 通过双阿累尼乌斯化学流变模型, 获得了季戊四醇丙烯醛树脂的粘度随时间及温度变化方程: η=5.30×10-6exp(4.71×104/RT)+exp((-1.152+3.44×10-3T)t)。该方程可反应固化过程中温度对体系粘度变化的影响规律, 所建立的模型与实验结果相关系数为0.9782, 两者一致性良好。
[1] |
于剑昆. 123树脂的合成、改性及在推进剂和炸药中的应用[J].
化学推进剂与高分子材料, 2011, 9(1): 31-39. YU Jian-kun. The synthesis, modifying and application in explosives and propellants of 123-resin[J]. Propellants and Macromolecule Materials, 2011, 9(1): 31-39. |
[2] |
汤业朋, 董海山, 陈淼鸿. RM-123型挤注炸药[J].
爆炸与冲击, 1981(2): 110-112. TANG Ye-peng, DONG Hai-shan, CHEN Miao-hong. RM-123 extrusion explosives[J]. Explosion and Shock Waves, 1981(2): 110-112. |
[3] |
Fischer R F, Smith C W. Cyclic acrolein acetals[J].
Journal of Organic Chemistry, 1960, 25(3): 319-324. DOI:10.1021/jo01073a002 |
[4] |
Kutsenko N I, Pertsov L D, Kalinkin S F, et al. Synthesis of diallyldiene pentaerythritol[J].
Plast Massy, 1967(2): 47-48. |
[5] |
古忠云, 袁昭慧, 李尚斌. 丁羟胶性能及其改性研究[J].
火炸药学报, 2002(3): 56-57. GU Zhong-yun, YUAN Zhao-hui, LI Shang-bin. Study on the properties and improvement of HTPB[J]. Chinese Journal of Explosives & Propellants, 2002(3): 56-57. |
[6] |
张泽义, 李笑江, 王晗, 等. PEG/ IPDI与PEG/TDI固化反应动力学研究[J].
含能材料, 2007, 15(4): 320-323. ZHANG Zhe-yi, LI Xiao-jiang, WANG Han, et al. Curing reaction kinetics of polyethylene glycol with isophorone diisocyanate and tolulene diisocyanates[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2007, 15(4): 320-323. |
[7] |
刘晶如, 罗运军. 非等温DSC研究Al /HTPB / TDI体系的固化反应动力学[J].
含能材料, 2009, 17(1): 83-86. LIU Jing-ru, LUO Yun-jun. Curing kinetics of HTPB/TDI/Al system by non-isothermal DSC[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2009, 17(1): 83-86. |
[8] |
Kissinger H E. Reaction kinetics in different thermo analysis[J].
Analytical Chemistry, 1957, 29: 1702-1706. DOI:10.1021/ac60131a045 |
[9] |
Ozawa T. A new method of analyzing thermo-gravimetric data[J].
Bulletin of the Chemical Society of Japan, 1965, 3(11): 1881-1886. |
[10] |
Crane L W. Analys is of curing kinetics in polymer composites[J].
Polymer Science, Polymer Letter Edition, 1973, 11: 533-543. DOI:10.1002/pol.1973.130110808 |
[11] |
段轶锋, 王小群, 刘羽中, 等. 几种聚醚胺改性蒙脱土对环氧树脂固化过程的影响[J].
化学学报, 2012, 70(10): 1179-1186. DUAN Yi-feng, WANG Xiao-qun, LIU Yu-zhong. Effect of montmorillonite modified by poly(oxypropylene) amines on curing process of epoxy nanocomposites[J]. Acta Chim Sinica, 2012, 70(10): 1179-1186. DOI:10.3969/j.issn.0253-2409.2012.10.005 |
[12] |
刘宏, 单国荣, 潘鹏举. 聚酯树脂粉末涂料的固化行为[J].
化工学报, 2012, 63(4): 1315-1320. LIU Hong, SHAN Guo-rong, PAN Ju-peng. Curing behavior of powder coating of polyester resin[J]. Chinese Journal of Chemical Engineering, 2012, 63(4): 1315-1320. |
[13] |
Roller M B. Rheology of curing thermosets: A review[J].
Polymer Engineering Science, 1986, 26: 432-440. DOI:10.1002/(ISSN)1548-2634 |
[14] |
Halley P J, Mackay M E. Chemorheology of thermosets-an overview[J].
Polymer Engineering Science, 1996, 36: 593-608. DOI:10.1002/(ISSN)1548-2634 |
Thenon-isothermal curing kinetics of acrolein-pentaerythritol resin system were studied by dynamic differential scanning calorimetry(DSC). And the curing kinetic model for acrolein-pentaerythritol resin system was established as well as the rheological model based on the viscosity experiments at 341, 347, 353 K and 359 K.