陈明华,黄伟佳,张洋洋,孙也尊
(1.中国人民解放军32181部队, 西安 710032; 2.中国人民解放军92337部队,辽宁 大连 116023;3.中国人民解放军32382部队,北京 100072)
热分析技术已经广泛应用于含能材料的性能研究中,汤崭等通过DSC的实验对单质炸药HMX的热性能进行研究,利用得到的实验数据,计算了HMX得热分解参数和机理函数,为HMX的安全性和应用提供理论依据;朱一举等利用TG的方法研究了湿热条件下的RDX/AP-NEPE推进剂的热安全性及力学性能的影响,表明了初期老化对表观活化能影响较大。目前,端羟基聚丁二烯复合固体推进剂广泛应用于多种弹药发动机中,复合固体推进剂的热性能研究是其生产、运输和储存等方面安全性能及相容性和寿命预测评估的基础。因此研究该型推进剂的热分解性能具有重要的意义。本文对该型推进剂进行了TG和DSC试验,对其热性能进行了测试,并计算推进剂的活化能利用DSC试验测得的数据求取热爆炸临界温度和自加速分解温度等热安全性参数;利用TG热分析得到该型推进剂热分解的反应速率d/d和反应深度等参数。采用Malek法推断推进剂的最概然机理函数,较为全面地评价了该型推进剂的热安定性和热安全性,为该型推进剂的安定性与寿命评估提供了重要依据。
HTPB/RDX/Al/AP四组元推进剂的主要组成成分为:端羟基聚丁二烯(HTPB)、铝(Al)、黑索金(RDX)和高氯酸铵(AP)其他成分包括固化剂异佛尔酮二异氰酸酯(IPDI)和增塑剂癸二酸二辛脂(DOS)等。
试验仪器:DSC 8000差示扫描量热仪,美国PE公司;Pyris-1型热重分析仪,美国PE公司。
试验条件:高纯氮气(99.999%);气压为0.3 MPa;气体流速20 mL·min,普通铝池卷边;升温速率分别为2.5 ℃·min、5 ℃·min、10 ℃·min和15 ℃·min,从100 ℃升温至 450 ℃。
..反应动力学参数计算
本文选用Flynn-Wall-Ozawa法和Kissinger法求解活化能,其方程式为:
1) 积分法
(5)
2) 微分法
(6)
式中,()为动力学机理函数,为反应转化率;为峰温(K);为指前因子;=8.314 J·mol·K,为摩尔气体常数;为升温速率(℃·min);为反应活化能 (10J·mol)。
运用以上的2种方法求解2个放热峰的反应活化能,结果如表1所示。根据平均活化能计算得到指前因子=5.645×10。
表1 推进剂峰温活化能(DSC)Table 1 activation energy of propellant peak temperature (DSC)
..热分解特征参数
图1为该型推进剂在设定的升温速率下得到的DSC曲线。
图1 推进剂DSC曲线
从图1可以看出推进剂在加热过程中出现了2个放热峰,分别在200~237 ℃和337~385 ℃范围内,说明推进剂在这两个温度范围内发生了剧烈的分解。在240~248 ℃温度范围内出现一个吸热峰。大约在400 ℃后推进剂不再出现吸(放)热峰,说明推进剂反应结束。由此可以看出:升温速率越快,放热峰峰值越大,说明放热越迅速;放热峰面积也随之增大,说明放热反应越充分。表2给出了各个升温速率下推进剂的外推起始温度、分解终止温度、峰顶温度和分解焓变Δ。
表2 各个升温速率下推进剂的DSC试验结果Table 2 DSC test results of propellants at various heating rates
根据多项式回归的方法得到更精确的外推起始温度。
(1)
表2中的峰温和温升速率代入公式当升温速率=0时,求得=193.80 ℃。
根据Zhang-Hu-Xie-Li法,可求热爆炸临界温度。
(2)
(3)
指升温速率=0时的放热峰值温度;、、为拟合系数。将表2中的峰温和温升速率代入式(3),求得,将其代入式(2),求得。
再根据公式:
(4)
求得自加速分解温度。
推进剂的第一个峰温在204.30~220.50 ℃之间,并且计算得到推进剂的外推起始温度=193.8 ℃,热爆炸临界温度=206.59 ℃,自加速分解温= 204.86 ℃,说明了推进剂的热安全性能较为良好。RDX约在200~250 ℃温度范围内的熔融状态下开始发生热分解,AP的热分解温度在330~400 ℃范围内。说明推进剂各个组分的混合体系并没有改变各成分的主要热分解行为,呈现出分阶段分解的规律。而推进剂热分解的外推起始温度、热爆炸临界温度和自加速分解温度与RDX接近,说明了RDX的加入,决定了整个混合体系第一阶段的主要热分解行为。从而可以确定推进剂加入RDX在增加了推进剂的燃烧性能和比冲的同时,也一定程度上降低了推进剂的热安全性。
..反应动力学参数计算
本文使用等转化率的动力学方法计算推进剂TG实验的活化能。
1) 积分法(Flynn-Wall-Ozawa法)
(7)
2) 微分法(Starink法)
(8)
表3 等转化率法计算得到活化能ETable 3 The activation energy E calculated by the equal conversion method
..热分解特征参数
从图2和图3可以看出,推进剂的热分解失重分为3个阶段,在150~220 ℃的范围内失重约质量的20%,这个阶段热分解反应速率最快,持续时间较短;在220~375 ℃的范围失重约质量的50%,在这个阶段推进剂持续较为匀速的热分解;在375~515 ℃的范围内失重约为20%,在这个阶段之后热分解结束;最终总的失重约为质量的90%,剩余约10%的残渣。对比DSC的结果可知,第1阶段主要为为RDX的热分解失重,同时,部分Al也将会与RDX发生氧化还原反应;并且RDX和AP最为强氧化剂和HTPB橡胶将发生“后固化”和强氧化作用,主要产物含碳的气相氧化物;第2个阶段主要为AP的热分解;第3个阶段为转变晶型后AP的高温分解和HTPB橡胶的热分解;剩余残渣主要为Al。表4给出了各个升温速率下推进剂的分解起始温度、最大损失速率、及其对应的温度和质量损失Δ%。
图2 推进剂TG曲线
图3 升温速率为15 ℃·min-1的推进剂TG-DTG曲线Fig.3 TG-DTG diagram of propellant with heating rate of 15 ℃·min-1
表4 TG各个升温速率下推进剂的TG试验结果Table 4 TG test results of propellants at Various heating rates at TG
可以看出,随着升温速率的增大,最大损失速率也依次增大,反应越来越剧烈,其对应的温度也逐渐升高,对比DSC的第一个放热峰峰温,可以看出两者相差较小,说明DSC与TG试验结果相符。
Malek法是定义函数()和()来确定机理函数()和()的一种方法。通过数学推导后得到()和()与()和()的关系式,利用TG和DTG的实验数据,从运用等转化率法计算活化能开始逐步得到动力学机理函数,其优点在于不需要额外的假设条件。
例如,在家属经商办企业情况清查工作中,各单位领导绝大多数都能够如实上报,但在基层干部中却出现了一些漏报瞒报情况。在专项检查中,基层班站长违规违纪问题也逐渐显现,部分班站长吃请收受礼品、截留员工薪酬工资、私自处理废旧生产物资,工作方式简单粗暴,虽有各种原因,但却反映出基层干部队伍建设存在一定短板。这就要求我们加强基层干部的依法合规培训和监督检查力度,将班站长教育管理、作风建设同履职能力培训、纪委举报案件调查相结合,多措并举、从严管理、压实责任,提升基层干部队伍素质的整体水平。
1) 定义函数()的表达式:
(9)
(10)
其中
(11)
当=0.5时:
(12)
2) 定义函数()的表达式:
()=()()
(13)
(14)
根据Luke近似式可知:
(15)
(16)
将升温速率为15 ℃·min的TG试验中各个反应深度的温度、反应速率d/d以及平均活化能代入上面的公式中,得到表5和表6。
表5 y(α)-α的数据Table 5 The data of y(α)-α
表6 Z(α)-α的数据Table 6 The data of Z(α)-α
对表5和表6的数据进行多项式拟合得到()-和()-的拟合曲线(如图4、图5)及关系多项式。其中()-曲线的置信水平为0.904,()-的置信水平为0.899,具有较高的置信水平。
()=2.199-12.836+32.088-22.202
(17)
()=0.179-1.041+2.593-1.792
(18)
图4 y(α)-α拟合曲线
图5 Z(α)-α拟合曲线
常见的动力学机理函数共有41种。对比分析可知,推进剂的最概然机理函数为47号函数,自催化的级反应,()=(1-)(1+)。与TG和DTG图对比分析可知,在0.1左右,dd最大,出现一个陡峰,反应最为迅速。观察反应速率变化发现0.3~0.65,出现了2个较大的峰值,在0.75~0.85又再出现一个反应速率峰值,与机理函数拟合曲线结果一致。说明该型推进剂的机理函数为自催化的级反应,=2。
由于自催化的级反应的积分函数不存在解析式,故根据反应速率方程,联立式(11)可得:
(19)
由于经过换算()和()的方程式相差很小,直接以()作为()()的乘积,最后求解得到:
(20)
1) 通过TG和DSC测试,计算得=206.59℃,自加速分解温度=204.86 ℃,推进剂活化能约在(187.30~205.51)×10J·mol,说明其热安全性和热安定性较为良好。
2) RDX的加入一定程度上降低了推进剂的热安全性和热安定性。
3) 由Malek法发现该型推进剂热分解是一个持续的分阶段的自催化级反应,=2,计算得到推进剂的动力学方程分别为: