火炸药学报    2018, Vol. 41 Issue (5): 489-495   DOI: 10.14077/j.issn.1007-7812.2018.05.012
0

引用本文  

唐秋凡, 屈蓓, 李吉祯, 张林, 张亚俊, 张正中, 樊学忠. 不同TMETN/NG配比对CMDB推进剂塑化特性的影响[J]. 火炸药学报, 2018, 41(5): 489-495. DOI: 10.14077/j.issn.1007-7812.2018.05.012
TANG Qiu-fan, QU Bei, LI Ji-zhen, ZHANG Lin, ZHANG Ya-jun, ZHANG Zheng-zhong, FAN Xue-zhong. Effect of Different TMETN/NG Ratios on the Plasticizing Characteristics of CMDB Propellant[J]. Chinese Journal of Explosives & Propellants, 2018, 41(5): 489-495. DOI: 10.14077/j.issn.1007-7812.2018.05.012

基金项目

国家自然科学基金(No.21401124)

作者简介

唐秋凡(1990-), 男, 硕士, 从事固体推进剂技术研究。E-mail:qiufantang@126.com

文章历史

收稿日期:2018-01-11
修回日期:2018-04-22
不同TMETN/NG配比对CMDB推进剂塑化特性的影响
唐秋凡1, 屈蓓1, 李吉祯1, 张林2, 张亚俊1, 张正中1, 樊学忠1     
1. 西安近代化学研究所, 陕西 西安 710065;
2. 空军驻西北地区军事代表室, 陕西 西安 710043
摘要: 为了提高三羟甲基乙烷三硝酸酯(TMETN)对NC的增塑效果,通过分子动力学模拟不同TMETN/NG配比条件下推进剂体系的回旋半径、自由体积等参数,研究NG部分取代TMETN对CMDB推进剂塑化特性的影响,并通过固化实验以及机械感度测试对模拟结果进行验证。结果表明,在计算范围内,随着NG含量的增加,体系的回旋半径、自由体积均先减小后增加,体系的塑化能力也是先减小后增大,且当NG取代质量分数40%及以上的TMETN时能有效改善体系的增塑效果;推进剂中NG与TMETN质量比大于2:3时,推进剂样品才能正常固化;NG与TMETN质量比控制在3:2~4:1时,既能保证CMDB推进剂固化均匀,又能降低推进剂的机械感度。
关键词: 分子动力学模拟     TMETN     NC     塑化特性     浇铸CMDB推进剂    
Effect of Different TMETN/NG Ratios on the Plasticizing Characteristics of CMDB Propellant
TANG Qiu-fan1, QU Bei1, LI Ji-zhen1, ZHANG Lin2, ZHANG Ya-jun1, ZHANG Zheng-zhong1, FAN Xue-zhong1     
1. Xi'an Modern Chemistry Research Institute, Xi'an 710065, China;
2. Army Representation Office of Air Force in Northwest, Xi'an 710043, China
Abstract: To improve the plasticizing effect of trimethylolethane trinitrate (TMETN) on NC, the cyclotron radius and free volume etc. parameters for propellant system under the conditions of different ratios of TMETN/NG was simulated by molecular dynamics. The effect of partly replacing TMETN with NG on the plasticizing characterizations of CMDB propellant was studied and the simulation results were validated by curing experiment and sensitivity test. Results show that within the calculation range, the cyclotron radius and free volume of the system decrease first and then increase with the increase of NG content. The plasticizing capacity of the system first reduces and then increases and when replacing TMETN of mass fraction as 40% with NG, the plasticizing effect of the system can be effectively improved. The propellant samples can be cured only when the mass ratio of NG and TMETN in the propellant is greater than 2:3. When the mass ratio of NG and TMETN is controlled at 3:2-4:1, it can not only guarantee the uniform curing of CMDB propellant, but also reduce the mechanical sensitivity of the propellant.
Keywords: molecular dynamics simulation     TMETN     NC     plasticizing characteristics     the casting CMDB propellants    
引言

随着各种大型舰船/航母、高性能战斗机、武装直升机及多功能战车等高价值平台武器的不断涌现[1-2],对武器系统感度提出了更高要求[3-5]。三羟甲基乙烷三硝酸酯(TMETN)是一种低感度含能增塑剂,其化学结构虽与硝化甘油(NG)相似,但其感度、毒性、挥发性和吸湿性均比NG小的多[6]。国外已将TMETN用于挤压双基系推进剂中,发现TMETN的加入可明显降低此类推进剂的机械感度[7-8]。国内赵凤起[9-10]等也研究了低感的TMETN增塑剂逐步取代双基系推进剂的NG对推进剂机械感度的影响,研究发现,随着NG被TMETN取代,推进剂的撞击感度和摩擦感度均明显降低。

纯TMETN对硝化纤维素(NC)的增塑效果较差,在70℃下很难固化成型,而且由于受实验尺度的限制,对改性双基推进剂塑化的微观结构变化等方面缺乏更深入的了解[11],这使得TMETN很难应用于浇铸改性双基推进剂中。分子动力学模拟方法作为近年来发展较为迅速的微尺度数值计算方法,因其能够准确预测材料的特性且从微观角度揭示材料结构与性能之间的关系,已成为研究火炸药宏观性质本质的一种方法[12-14]

本研究通过分子动力学模型模拟了不同NG含量对TMETN/NC体系塑化特性的影响,研究NG部分取代TMETN对CMDB推进剂塑化特性的影响。通过计算回旋半径、自由体积和径向分布函数等参数从理论上分析在保证浇铸改性双基(CMDB)推进剂塑化工艺过程的条件下,TMETN能取代NG的最大量,并结合实验进行验证,为实现TMETN/NC体系在浇铸改性双基推进剂中的应用及含TMETN不敏感浇铸改性双基推进剂的相关研究提供数据支持。

1 分子动力学模拟

依据NC、TMETN和NG的化学结构式,用美国Accelrys公司开发的Materials Studio 7.0软件建立NC、TMETN和NG分子物理模型。NC、TMETN和NG分子模型见图 1。用Amorphous Cell/Construction构建包含3条NC分子链(每条分子链有20个聚合单元组成)的NC纯物质模型(3条分子链分别标记为NC-A、NC-B和NC-C)以及NG/NC(质量比为1:1)、TMETN/NC(质量比为1:1)和NG/TMETN/NC(质量比分别为1:4:5、2:3:5、3:2:5、4:1:5)共混物模型。

图 1 NC、TMETN和NG分子模型 Figure 1 Molecular models of NC, TMETN and NG

采用Smart minimization对不同混合体系进行5 000步能量最小化优化,优化过程采用convergence level of medium。采用Forcit模块,在系综为NPT系综、298 K和101.325 kPa条件下进行400 ps的分子动力学平衡以获得平衡密度,时间步长为1 fs。NG/NC、TMETN/NC和NG/TMETN/NC的混合模型见图 2(NC-A、NC-B和NC-C分别标记为红色NC链、蓝色NC链和黄色NC链)。

图 2 NG/NC、TMETN/NC和NG/TMETN/NC混合模型 Figure 2 Models of NG/NC, TMETN/NC and NG/TMETN/NC
2 实验 2.1 样品

TMETN和NG,化学纯,西安近代化学研究所;NC(氮含量12%),泸州北方硝化棉公司;RDX,兰州白银银光化学材料厂;配方中其他组分包括弹道稳定剂、催化剂和安定剂均为工业品。

2.2 推进剂样品制备

推进剂样品采用淤浆浇铸工艺制备,配方见表 1。将样品置于70℃烘箱固化72 h,通过观察固化情况来分析不同TMETN/NG配比对CMDB推进剂的塑化效果。

表1 CMDB推进剂配方 Table 1 Formulations of CMDB propellants
2.3 感度测试

采用GJB770B-2005方法602.1对推进剂样品进行摩擦感度测试,测试条件为:表压2.45 MPa,摆角66°,药量20 mg,摩擦感度用爆炸概率P表示。采用GJB770A-1997方法601.2对推进剂样品进行撞击感度测试,测试条件为:落锤质量2 kg,药量30 mg,撞击感度用特性落高H50表示。

3 结果与讨论 3.1 均方位移分析

回旋半径指线性聚合物分子链中每个链节与分子链质心之间距离的统计平均值,是能够直接反映线性分子链构象的特征参数[15-16]。由于NG和TMETN与NC分子间的相互作用会影响NC分子链的构象,同时考虑到温度也是决定NC分子链构象的一个重要因素,因此通过回旋半径来研究NC分子链构象受NG和TMETN分子数量和温度影响的变化,进而得到NG和TMETN与NC分子间相互作用的直观映像[17]。不同TMETN/NG配比的NG/TMETN/NC共混体系模型中,NC分子链(NC-A、NC-B和NC-C)在不同温度下的回旋半径(r)分布趋势图见图 3

图 3 不同TMETN/NG复配体系中NC的回旋半径分布趋势图 Figure 3 Trend charts of cyclotron radius distribution of NC in different TMETN/NG blending systems

图 3可知,随着温度的升高,不同体系模型中NC分子链的回旋半径均增加,可能是因为升高温度使高分子链段运动加剧,导致高分子链更加伸展,回旋半径增加[18]。纯NG体系的回旋半径明显高于纯TMETN体系,说明与TMETN相比,NG的加入更能促进NC链的尺寸增大,回旋半径也随之增加,纯NG对NC的增塑效果要优于纯TMETN。同时,在NG逐渐取代TMETN的过程中,不同体系模型中的NC链(NC-A、NC-B和NC-C)的回旋半径均呈现先减小后增加的趋势,在加入质量分数20%~40%的NG取代TMETN时,回旋半径达到最小值。

NC是由葡萄糖酐环状残基组成的线性聚合物,具有一定的刚性,但其环间的醚链又使链节的内旋转比较容易,所以NC分子链具有一定的内旋自由度,但这种内旋自由度会受到NC分子链中各原子间作用力的抑制[19]。由此可以推测,NG分子与TMETN分子对NC分子链内部的作用力均有一定的弱化作用,而且NG对NC链的弱化作用强于TMETN分子,所以纯NG体系的回旋半径明显高于纯TMETN体系,纯NG对NC的增塑效果要高于纯TMETN。但是在NG/TMETN复配体系中,随着质量分数20%~40%的TMETN被NG取代,NG与TMETN的相互作用会减弱这种弱化作用,导致其内旋自由度减小,NC分子链尺寸减小,回旋半径减小。但随着NG量的增加,由于NG本身对NC分子链内部作用力的强弱化作用,导致其内旋自由度增大,NC分子链尺寸增大,回旋半径又增大。所以,要改善纯TMETN对NC的增塑效果,混合体系NC的回旋半径要大于纯TMETN体系,即NG至少要取代质量分数40%以上的TMETN。

3.2 自由体积分析

一般认为,聚合物材料中的体积可分成两部分,即原子所占有的部分和自由体积部分。自由体积部分包括分子无规则堆砌产生的原子尺寸的缺陷等[20]。自由体积的尺寸、数量以及分布等决定了分子链的活动能力,同时分子链的活动能力也对自由体积有着重要的影响[21-22]。单个NG分子、TMETN分子和NC分子的体积见图 4,分别为0.161、0.204和0.183 nm3。本研究分别在293 K和323 K条件下,对NG/NC、TMETN/NC和NG/TMETN/NC的混合体系进行分子动力学模拟,得到几种模型的体积分布见图 5表 2

图 4 单个NG分子、TMETN分子和NC分子的体积分布 Figure 4 Volume distribution of NG, TMETN and NC molecules
图 5 293 K和323 K下不同共混体系的体积分布 Figure 5 Volume distribution of different blending systems at same temperatures
表2 不同温度下共混体系中分子的占有体积和自由体积 Table 2 Occupancy volume and free volume of molecule in blending system at different temperatures

图 5中白色部分为NC和增塑剂分子的占有体积(Vo),深色部分为分子之间的空隙,定义为自由体积(Vf)。由表 2可知,随温度的升高,所有体系的Vf相应增加;且由于分子间的间隙增大,导致部分重叠体积暴露出来,使探针(半径为0.1 am)所计算的Vo也略有增加[19]。同时由图 5表 2数据可知,在NG逐渐取代TMETN的过程中,自由体积Vf和占有体积Vo均呈现先略有减小后增加的趋势。可见,当体系中NG质量分数超过20%后,起到了部分与升高温度相类似的作用,增加了体系的自由体积,使其在受到外力作用时,容纳NC分子链链段内旋转、跃迁甚至整链滑移的空间增大,形变量相应增加[19]

3.3 径向分布函数分析

径向分布函数(RDF)是用来描述所研究的粒子周围其他粒子的分布情况的物理量,它能有效的反映模型中不同分子或者原子之间的相互作用强弱[23-24]。通过分析293K条件下H—N以及H—O之间的径向分布函数,研究随着NG含量的增加,TMETN与NC之间及TMETN与NG之间的相互作用,结果见图 6图 7。TMETN分子中的H、O和N原子用H(1)、O(1)和N(1)表示,NG和NC的H、O和N原子分别用H(2)、O(2)和N(2)表示。

图 6 TMETN与NC之间的径向分布函数 Figure 6 Radial distribution function between TMETN and NC
图 7 TMETN与NG之间的径向分布函数 Figure 7 Radial distribution function between TMETN and NG

图 6图 7可知,随着NG含量的增加,TMETN与NC之间的径向分布函数峰值(1H-2O中0.22 nm附近和0.26 nm附近的前两个峰以及1H-2N中0.27 nm附近和0.33 nm附近的前两个峰)逐渐减小,说明TMETN与NC之间的相互作用随着NG含量的增加逐渐减弱;同时,TMETN与NG之间的径向分布函数峰值(1H-2O中0.21 nm附近和0.26 nm附近的前两个峰以及1H-2N中0.27 nm附近和0.34 nm附近的前两个峰)也逐渐减小,说明TMETN与NG之间的相互作用随着NG含量的增加也在逐渐减弱。可能是少量NG的加入减弱了TMETN与NC、TMETN与NG之间的相互作用,使回旋半径、自由体积均呈现微弱的减小趋势。但随着NG含量继续增加,NG对NC的强增塑作用使回旋半径、自由体积又呈现增加的趋势。

3.4 NG对TMETN推进剂样品固化效果的影响

分别对L4~L6号CMDB推进剂样品固化结果进行观察,截面图如图 8所示。并测试固化后样品的感度,结果如表 3所示。

图 8 L4~L6推进剂截面图 Figure 8 Section of propellants L4-L6
表3 不同CMDB推进剂样品感度测试结果 Table 3 The mechanical sensitivities of different propellant samples

L1~L3推进剂样品未能固化;从表 3图 8可以看出,L4推进剂样品虽然已固化,但在切面图上可以看出很明显的分层现象(上层基本没有白色的斑点),说明固化不均匀;而L5和L6推进剂样品切面图无明显分层,说明固化较为均匀。可见,推进剂中NG与TMETN质量比大于2:3时,推进剂样品能正常固化,但NG与TMETN质量比为3:2时固化效果不理想。当NG与TMETN质量比大于3:2时,推进剂样品能均匀固化,这一结论与模拟结果一致。对比感度测试结果可知,随着TMETN含量的增加,推进剂摩擦感度和撞击感度均得到了改善。综合来看,NG与TMETN质量比控制在3:2~4:1,既能保证推进剂工艺性能,同时又能降低推进剂的机械感度。

4 结论

(1) 通过分子动力学模拟的回旋半径、自由体积等参数分析发现,NG取代TMETN的过程中,随着NG含量的增加,各模拟参数均呈现先减小后增加的趋势。要改善纯TMETN对NC的增塑效果,NG/TMETN体系的回旋半径、自由体积均要大于纯TMETN体系,即NG至少要取代质量分数40%以上的TMETN。

(2) 加入少量的NG会减弱TMETN与NC、TMETN与NG之间的相互作用,这两方面的作用减小了NC的回旋半径以及整个体系的自由体积,使容纳NC分子链链段运动的空间减小,降低了NC分子链的链段移动性。在宏观上使推进剂(NG与TMETN质量比小于2:3)表现出未能固化的现象。

(3) 推进剂中NG与TMETN质量比大于2:3时,推进剂样品能正常固化;NG与TMETN质量比大于3:2时,推进剂样品能均匀固化。这与分子模拟结果一致。

(4) NG与TMETN质量比控制在3:2~4:1,既能保证CMDB推进剂工艺性能,同时又能改善推进剂的感度性能。


参考文献
[1]
Sikder A K, Sikder N. A review of advanced high performance, insensitive and thermally stable energetic materials emerging for military and space applications[J]. Journal of Hazardous Materials, 2004, A112: 1–15.
[2]
Badgujar D M, Talawar M B, Asthana S N, et al. Advances in science and technology of modern energetic materials:an overview[J]. Journal of Hazardous Materials, 2008, 151: 289–305. DOI:10.1016/j.jhazmat.2007.10.039
[3]
付小龙, 樊学忠. 钝感推进剂配方研究及发展趋势[J]. 火炸药学报, 2014, 37(5): 1–8.
FU Xiao-long, FAN Xue-zhong. Research and development trends of insensitive solid propellant formalation[J]. Chinese Journal of Explosives & Propellants (Huozhayao Xuebao), 2014, 37(5): 1–8. DOI:10.3969/j.issn.1007-7812.2014.05.001
[4]
庞爱民, 郑剑. 高能固体推进剂技术未来发展展望[J]. 固体火箭技术, 2004, 27(4): 289–293.
PANG Ai-min, ZHENG Jian. Prospect of the research and development of high energy solid propellant technology[J]. Journal of Solid Rocket Technology, 2004, 27(4): 289–293. DOI:10.3969/j.issn.1006-2793.2004.04.012
[5]
庞维强, 张教强, 国际英, 等. 21世纪国外固体推进剂的研究与发展趋势[J]. 化学推进剂与高分子材料, 2005, 3(3): 16–21.
PANG Wei-qiang, ZHANG Jiao-qiang, GUO Ji-ying, et al. Research and development trends of solid propellants abroad in the 21st century[J]. Chemical Propellants & Polymeric Materials, 2005, 3(3): 16–21. DOI:10.3969/j.issn.1672-2191.2005.03.004
[6]
赵凤起, 徐司雨, 仪建华, 等. 含TMETN的钝感推进剂燃烧特性数值模拟[J]. 含能材料, 2006, 14(6): 406–410.
ZHAO Feng-qi, XU Si-yu, YI Jian-hua, et al. Numerical simulation for combustion characteristics of insensitive propellant containing trimethylolethane trinitrate (TMETN)[J]. Chinese Journal of Energetic Materials, 2006, 14(6): 406–410. DOI:10.3969/j.issn.1006-9941.2006.06.002
[7]
Oyymi Y, Kimura E. Insensitive munitious and combustion characteristics of GAP/AN composition propellant[J]. Propellants, Explosives, Pyrotechnics, 1996, 121(5): 271–276.
[8]
石小兵, 庞维强, 蔚红建. 钝感推进剂研究进展及发展趋势[J]. 化学推进剂与高分子材料, 2007, 5(2): 24–32.
SHI Xiao-bing, PANG Wei-qiang, WEI Hong-jian. Research progress and development trends of insensitive propellant[J]. Chemical Propellants & Polymeric Materials, 2007, 5(2): 24–32. DOI:10.3969/j.issn.1672-2191.2007.02.005
[9]
赵凤起, 杨栋, 李上文, 等. 以NC和TMENT为基的微烟推进剂机械感度研究[J]. 火炸药学报, 1999, 22(4): 5–8.
ZHAO Feng-qi, YANG Dong, LI Shang-wen. Impact and friction sensitivity of minimum smoke propellant based on NC and TMETN[J]. Chinese Journal of Explosives & Propellants (Huozhayao Xuebao), 1999, 22(4): 5–8. DOI:10.3969/j.issn.1007-7812.1999.04.002
[10]
赵凤起, 陈沛, 李上文, 等. 钾盐消焰剂与TMETN和燃烧催化剂相互作用的实验研究[J]. 含能材料, 2001, 9(3): 100–103.
ZHAO Feng-qi, CHEN Pei, LI Shang-wen, et al. Interaction of potassium salt flame suppressors with TMETN and burning catalysts during decomposition[J]. Energetic Materials, 2001, 9(3): 100–103. DOI:10.3969/j.issn.1006-9941.2001.03.002
[11]
孟玲玲, 齐晓飞, 王江宁, 等. DNTF对NC塑化特性的分子动力学模拟及实验研究[J]. 火炸药学报, 2015, 38(3): 86–89.
MENG Ling-ling, QI Xiao-fei, WANG Jiang-ning, et al. Molecular dynamics simulation and experimental study of DNTF on plasticizing properties of NC[J]. Chinese Journal of Explosives & Propellants (Huozhayao Xuebao), 2015, 38(3): 86–89.
[12]
Pan R, Liu X, Zhang A, et a1. Molecular simulation on structure property relationship of polyimides with methylene spacing groups in biphenyl side chain[J]. Computation Materials Science, 2007, 39(4): 887–895. DOI:10.1016/j.commatsci.2006.10.019
[13]
Tang Q F, Fan X Z, Li J Z, et al. Experimental and theoretical studies on stability of new stabilizers for N-Methyl-P-nitroaniline derivative in CMDB propellants[J]. Journal of Hazardous Materials, 2017, 327: 187–196. DOI:10.1016/j.jhazmat.2016.12.049
[14]
Zhao Yu, Zhang Xiao-hong, Zhang Wei, et al. Simulation and experimental on the solvation interaction between the GAP matrix and insensitive energetic plasticizers in solid propellants[J]. The Journal of Physical Chemistry A, 2016, 120: 765–770. DOI:10.1021/acs.jpca.5b10540
[15]
Chen P K, Yao L, Liu Y Y, et al. Experimental and theoretical study of dilute polyacrylamide solutions:effect of salt concentration[J]. J Mol Model, 2012, 18: 3153–3160. DOI:10.1007/s00894-011-1332-9
[16]
Liu Y Y, Chen P K, Luo J L, et al. Molecular simulation of dilute polyacrylamide solutions[J]. Acta Phys Chim Sin, 2010, 26: 2907–2914.
[17]
齐晓飞, 张晓宏, 李吉祯, 等. NC/NG共混体系的分子动力学模拟研究[J]. 兵工学报, 2013, 34(1): 93–99.
QI Xiao-fei, ZHANG Xiao-hong, LI Ji-zhen, et al. Molecular dynamics simulation of NC/NG blends[J]. Acta Armamentarii, 2013, 34(1): 93–99.
[18]
Rubinstein M, Colby R. Polymer Physics[M]. Oxford: Oxford University Press, 2002.
[19]
齐晓飞, 张晓宏, 郭昕, 等. NENA对NC溶塑作用的实验与模拟[J]. 固体火箭技术, 2013, 36(4): 516–520.
QI Xiao-fei, ZHANG Xiao-hong, GUO Xin, et al. Experiments and simulation on plastication of NENA on NC[J]. Journal of Solid Rocket Technology, 2013, 36(4): 516–520.
[20]
周成飞. 聚氨酯自由体积及其对性能的影响[J]. 聚氨酯工业, 2014, 29(3): 1–5.
ZHOU Cheng-fei. The free volume of polyurethane and effect on its properties[J]. Polyurethane Industry, 2014, 29(3): 1–5. DOI:10.3969/j.issn.1005-1902.2014.03.001
[21]
卓启疆. 聚合物自由体积[M]. 成都: 成都科技大学出版社, 1990.
[22]
倪凯, 廖霞, 赵雪燕, 等. 等规聚丙烯在户外老化过程中化学重结晶作用和自由体积变化的研究[J]. 塑料工业, 2015, 43(7): 69–72.
NI Kai, LIAO Xia, ZHAO Xue-yan, et al. Investigation on chemi-recrystallization and free volume changes of isotactic polypropylene during outdoor weathering[J]. China Plastics Industry, 2015, 43(7): 69–72. DOI:10.3969/j.issn.1005-5770.2015.07.017
[23]
Fu Xiao-long, Fan Xue-zhong, Ju Xue-hai, et al. Molecular dynamic simulations on the interaction between an HTPE polymer and energetic plasticizers in a solid propellant[J]. RSC Advances, 2015, 5: 52844–52851. DOI:10.1039/C5RA05312A
[24]
焦东明, 杨月诚, 强洪夫, 等. 铝粉氧化对端羟基聚丁二烯界面吸附影响的分子模拟[J]. 火炸药学报, 2009, 32(6): 79–83.
JIAO Dong-ming, YANG Yue-cheng, QIANG Hong-fu, et al. Molecular simulation of effect of aluminum powder oxidation on interface adsorption for HTPB[J]. Chinese Journal of Explosives & Propellants (Huozhayao Xuebao), 2009, 32(6): 79–83. DOI:10.3969/j.issn.1007-7812.2009.06.019