文章快速检索     高级检索
  含能材料  2018, Vol. 26 Issue (6): 483-488.  DOI: 10.11943/j.issn.1006-9941.2018.06.004
0

引用本文  

陈京, 刘萌, 何琦文, 马义丁, 徐抗震, 樊学忠. TMETN/NG混合增塑剂及NC共混体系的性能模拟[J]. 含能材料, 2018, 26(6): 483-488. DOI: 10.11943/j.issn.1006-9941.2018.06.004.
CHEN Jing, LIU Meng, HE Qi-wen, MA Yi-ding, XU Kang-zhen, FAN Xue-zhong. Simulation on the Properties of TMETN/NG Mixed-plasticizers/NC Blends[J]. Chinese Journal of Energetic Materials, 2018, 26(6): 483-488. DOI: 10.11943/j.issn.1006-9941.2018.06.004.

基金项目

国家自然科学基金项目(21703167)

作者简介

陈京(1988-), 男, 助理研究员, 主要从事固体推进剂研究。e-mail: chenjing_mcri@163.com

通信联系人

樊学忠(1962-), 男, 研究员, 主要从事固体推进剂配方及工艺技术研究。e-mail: xuezhongfan@126.com

文章历史

收稿日期:2017-10-21
修回日期:2017-12-27
TMETN/NG混合增塑剂及NC共混体系的性能模拟
陈京 1, 刘萌 1, 何琦文 1, 马义丁 1, 徐抗震 2, 樊学忠 1     
1. 西安近代化学研究所, 陕西 西安 710065;
2. 西北大学化工学院, 陕西 西安 710069
摘要:为了探索高能钝感混合增塑剂与含能粘合剂相互作用对于体系性能的影响, 以硝化棉(NC)为粘合剂, 构建了含有三羟甲基乙烷三硝酸酯(TMETN)、硝化甘油(NG)以及三种不同质量比TMETN/NG混合增塑剂的五种共混模型。采用Materials Studio©软件, 利用分子动力学(MD)方法模拟了各共混体系的相互作用、力学性能、径向分布函数、O—NO2最大引发键长(Lmax)等参数。结果表明, 共混模型中TMETN含量增加时, 共混体系的结合能降低。含有混合增塑剂的共混模型力学性能优于单组分增塑剂共混模型。NG与TMETN质量比为3:1时, 两种硝酸酯增塑剂可以同时与NC建立强的氢键作用, 并为共混体系提供较高的模量与较低的最大引发键长。结果显示,调整混合增塑剂质量比可以增强体系各组分的相互作用, 从而提升共混体系的力学性能及引发键的稳定性。
关键词分子动力学(MD)     混合增塑剂     力学性能     感度     氢键     三羟甲基乙烷三硝酸酯(TMETN)    
Simulation on the Properties of TMETN/NG Mixed-plasticizers/NC Blends
CHEN Jing 1, LIU Meng 1, HE Qi-wen 1, MA Yi-ding 1, XU Kang-zhen 2, FAN Xue-zhong 1     
1. Xi′an Modern Chemistry Research Institute, Xi′an 710065, China;
2. School of Chemical Engineering, Northwest University, Xi′an 710069, China
Abstract: To explore the effect of interactions between insensitive high energy mixed-plasticizers and energetic binders on the system, taking nitrocellulose(NC) as binder, five kinds of blending models which contain 2-methyl-2-propane-1, 3-diyl dinitrate(TMETN), nitroglycerin(NG), and three kinds of TMETN/NG mixed-plasticizers with different mass ratio were constructed. The interactions, mechanical properties, radial distribution functions, and maximum trigger bond(O—NO2) length (Lmax)etc. parameters for each blending system were simulated by molecular dynamics method using Materials studio software. Results show that the binding energy of blending system is decreased when the content of TMETN is increased in the blending models. The mechanical properties of blending model containing mixed-plasticizers are better than those of blending model of single component plasticizer. When the mass ratio of NG to TMETN is 3:1, two kinds of nitrate plasticizers can establish strong hydrogen bond interactions with NC at the same time, and provide higher modulus and lower maximum trigger bond length for blending system, revealing that adjusting the mass ratio of mixed-plasticizers can enhance the interactions between compositions of each system thus the mechanical properties and stability of trigging bond of the blending system are improved.
Key words: molecular dynamics(MD)    mixed-plasticizers    mechanical properties    sensitivity    hydrogen bond    2-methyl-2-propane-1, 3-diyl dinitrate(TMETN)    
1 引言

固体推进剂未来发展的重要趋势之一是兼具高能及钝感特性[1], 然而含能材料的能量与感度往往是一对矛盾的性质[2]。为开发高能钝感推进剂, 除了引入新型钝感含能材料之外[3], 采用混合增塑剂并调节其高能组分和钝感组分的构成与比例, 可有效平衡推进剂的能量与感度性能。大量研究显示, 混合增塑剂有助于降低复合含能体系的感度并提升力学性能[4-6], 因而在推进剂中得到了广泛的应用[7]

三羟甲基乙烷三硝酸酯(TMETN)自从被法国应用于钝感改性双基推进剂以来[8], 在推进剂领域得到了大量的研究和应用[9-11]。研究显示使用TMETN部分取代硝化甘油(NG)后, 可有效降低推进剂的机械感度[12], 然而完全取代后摩擦感度又有所增加[1]。可见TMETN/NG混合增塑剂体系具有比单组分增塑剂更好的降感作用。

目前, 混合增塑剂的降感原理仍不清楚, 采用实验方法进行研究存在缺乏理论指导、效率低、安全性差等问题, 需结合理论计算手段开展研究。目前对于复合含能体系的模拟主要采用分子动力学(MD)方法[13], 然而尚无研究涉及混合增塑剂体系。为研究TMETN/NG混合增塑剂与硝化棉(NC)粘合剂体系的相容性、力学性能及感度机理, 本工作采用MD模拟的方法, 构建五种共混模型, 用以分析不同比例TMETN/NG混合增塑剂与粘合剂的相互作用, 探索共混体系的感度与力学性能的影响因素, 为研究含有新型混合增塑剂的推进剂感度影响机理提供理论支撑。

2 计算方法 2.1 共混模型构建

首先使用Materials Studio(MS)软件中的Visualizer模块构建各增塑剂与粘合剂的结构模型:依据固体推进剂常用粘合剂硝化棉(NC)的含氮量标准(N质量分数11%~12%), NC分子纤维素母链中的三硝基取代重复单元与二硝基取代重复单元特征比约为1:6(N质量分数为11.7%)。因此依照1:6的特征比进行随机聚合建模, 参考常用分子动力学模型原子总数及NC高分子链长[14-15], 构建2条原子数为610的NC高分子单链模型。而后依次构建NG与TMETN的分子模型, 并将所有组分的模型进行能量优化。

利用Amorphous Cell模块, 在200 Å×200 Å×200 Å的周期箱中构建各组分的纯物质模型及粘合剂-增塑剂共混模型。按照改性双基推进剂中的常用增塑剂-粘合剂质量比, 将共混模型的增塑剂含量控制在57%, NC含量控制在43%。共混模型中由NG或TMETN单一增塑剂与NC组成的双组分模型为两种, 由NG/TMETN混合增塑剂与NC组成的三组分模型为三种, 其增塑剂质量比分别为NG:TMETN=3:1, 1:1, 1:3。各共混模型命名、组分分子数及原子总数见表 1

表 1 共混模型分子数、增塑剂质量比及原子总数 Tab.1 The molecular numbers, plasticizers mass ratios, and total atom numbers of blends
2.2 动力学平衡

将各纯物质模型与共混模型在Forcite模块中进行分子动力学平衡, 选择等温等压系综(NPT), 在298 K和101.325 kPa条件下进行400 ps平衡计算。分子起始速度: Maxwell-Boltzmann随机分布给定, Velocity Verlet算法进行求解, 时间步长为1 fs; 模拟过程温度采用Anderson方法, 范德华力和静电作用分别采用Atom-based和Ewald加和方法, 非键截取半径为0.95 nm并进行截断尾部校正。动力学平衡的后200 ps能量接近平衡, 统计分子运动轨迹及各参数用于共混体系的性能分析。

3 结果与讨论 3.1 组分间结合能分析

各组分分子结构与共混模型平衡结构见图 1。由图 1可以看出, NG、TMETN、NG/TMETN混合增塑剂与NC各组分分子混合均匀, 显示出良好的混溶性。

图 1 各组分分子结构及共混模型平衡结构 Fig.1 The molecular structures of compositions and equilibrium structures of blends

结合能是共混体系各组分之间相互作用的重要特征参数, 可以由体系的混合总能量和各组分的单点能来计算。NC粘合剂和硝酸酯增塑剂的结合能Eb可以由式Eb=ENC+Ep-Etotal计算。式中, Etotal为共混体系的总能量, ENC为NC纯物质模型平衡状态的单点能, Ep为增塑剂纯物质模型平衡状态的单点能。本文构建的NC与各增塑剂共混模型的能量数据见表 2

表 2 NC与各增塑剂体系之间的结合能 Tab.2 The binding energy between plasticizers and NC

比较单组分增塑剂与NC的共混体系能量数据, NG与NC的结合能要大于TMETN与NC的结合能, 这是由于同质量的NG增塑剂含有比TMETN更高浓度的—NO2基团, 而—NO2基团可以和NC的—OH基团产生氢键作用[16-17], 并且可与NC中的—NO2基团良好相容。对比三种混合增塑剂体系的结合能数据可以发现, 随着TMETN的加入, NG浓度的降低, 增塑剂与NC的结合能逐渐减小。NC/NG/TMETN-1在三种混合增塑剂共混模型中相互作用最强, 通过表 1中分子数计算可知模型NC/NG/TMETN-1/2/3中的—NO2基团数分别为210/207/201(每个硝酸酯分子含有3个—NO2基团), 因此这个现象同样归因于模型中—NO2基团的浓度差异带来的相互作用差异。

3.2 力学性能模拟

弹性常数是反映粘合剂与增塑剂分子之间相互作用各向异性的表征参数。纯NC与5个共混体系的弹性常数和各向同性力学性能的计算值见表 3

表 3 各共混模型的弹性常数和力学性能 Tab.3 The elastic constants and mechanical properties of blends

表 3可以看出, 纯NC粘合剂具有较高的拉伸模量(E)、剪切模量(K)和体积模量(G), 而与增塑剂混合形成共混体系后, 各模量均有明显降低。这表明增塑剂可增加NC粘合剂的韧性, 与NC塑化后变软的实验现象相符。在三种采用混合增塑剂的共混体系中, 随着TMETN含量的增加, 共混体系的模量呈下降趋势, NC/NG/TMETN-1体系具有相对较高的模量。这是由于—NO2基团浓度的差异, 使得TMETN增塑剂与NC粘合剂之间的相互作用较NG弱, 与表 2的结合能模拟结果显示的规律一致。然而, 相比于单组分增塑剂共混体系, 混合增塑剂共混体系的模量依然有明显的提升, 这是由于混合增塑剂与NC共混之后更好地减少了粘合剂的各向异性[18]。此外, 各共混模型的柯西压(C12C44)均为正值, 表明共混体系均具有一定延展性, 在推进剂配方中具有可加工性。值得注意的是, NC/NG/TMETN-1体系的柯西压显著优于其他共混体系, 可见添加少量的TMETN即可大幅改善NC/NG体系的延展性, 添加过多却会导致延展性下降。

3.3 径向分布函数模拟

径向分布函数(RDF)为某一特定原子附近另外一种原子出现的几率, 其表达式为

$ g(r){\rm{ }} = 4{\rm{ \mathsf{ π} }}{r^2}\rho {\rm{d}}r。$

式中, r为原子间距, ρ为原子数量密度。从径向分布函数的峰值位置可以判断相互作用的类型, 而相互作用力大小可以从峰的高度进行推断[19]。因此可通过径向分布函数对氢键作用进行统计, 从而研究不同增塑剂与NC粘合剂之间的相互作用规律。

NC分子链结构中的主要侧基为—OH与—NO2, 其中—OH氧原子(力场为o2h)可以同时与增塑剂的—NO2氧原子(力场为o12)及NC分子链中的—NO2氧原子(力场为o12)产生氢键作用。由于NC链的分子体积较大且具有一定刚性, 自身之间氢键作用相对较弱, 因此本文主要考察小分子增塑剂与NC之间的氢键作用。

首先对NC/NG和NC/TMETN的平衡共混模型进行径向分布函数分析。图 2显示了NC的羟基氧原子(o2h)和两种增塑剂的硝基氧原子(o12)之间的径向分布函数曲线。由图 2可知, 两种体系均约在2.5 Å处存在尖锐的峰, 表明NG与TMETN都在与NC混合的过程中形成了较强的氢键作用, 位于约4.9 Å的峰显示增塑剂与NC之间也存在相对较弱的范德华力。

图 2 NC羟基氧原子(o2h)与不同单组分增塑剂硝基氧原子(o12)之间的径向分布函数 Fig.2 The radial distribution function between O atom from NC hydroxyl-group and O atom from single plasticizer nitro-group

为了研究混合增塑剂与NC之间相互作用的影响, 通过径向分布函数曲线分别分析共混模型中两种增塑剂与NC之间的相互作用, 结果见图 3。由图 3可知, 在共混体系NC/NG/TMETN-1(NG:TMETN=3:1)中, 两种增塑剂的硝基氧原子(o12)与NC中的羟基氧原子(o2h)的径向分布函数曲线在~2.5处均出现尖锐的峰, 这表明增塑剂各组分与粘合剂之间均存在较强的氢键作用。而在共混模型NC/NG/TMETN-2和NC/NG/TMETN-3中, 仅有一种增塑剂能够与NC形成明显的氢键, 另外一种增塑剂与NC之间仍以范德华力为主(图 3)。在NG与TMETN质量比为1:1时, TMETN的浓度相比于NC/NG/TMETN-1有所上升, 然而无法与NC产生明显的氢键作用。当NG与TMETN质量比达到1:3时, NG与NC的氢键作用变的不明显。

图 3 各共混模型中NC羟基氧原子(o2h)与不同增塑剂硝基氧原子(o12)之间的径向分布函数 Fig.3 The radial distribution function between O atom from NC hydroxyl-group and O atom from mixed-plasticizer nitro-group in blends

对比以上三种共混模型可以发现, 混合增塑剂体系各组分与NC的氢键作用强弱并非简单随—NO2基团浓度增加而增强, 而与增塑剂相互比例有直接关系。这一规律表明, 在特定比例下, 混合增塑剂的—NO2浓度可以与粘合剂的—OH基团达到匹配浓度, 从而产生较强的相互作用。此结果与共混体系力学性能模拟结果保持一致。

3.4 最大引发键长模拟

含能材料的引发键是分子中最弱的化学键, 在外界热、机械等刺激下将首先断裂, 进而引发含能材料的热分解、燃烧或爆炸。最大引发键长(Lmax)已经被用于各类含能材料感度性能的判据, 并尤其适用于多组分体系的感度预估。因此本文统计了各共混体系中硝酸酯增塑剂引发键(O—NO2键)在达到动态平衡后(200~400 ps)的键长数据, 用以表征不同共混体系的感度差异。

各共混体系中硝酸酯增塑剂的最大引发大键长见表 4。由表 4可见, NC/NG体系的敏感组分NG最大键长1.576 Å, 明显大于NC/TMETN模型的1.563 Å。可见使用低感度的TMETN替代NG与NC混合后, 可使共混体系的感度下降, 这与现有的研究结论一致[12]。随着TMETN在增塑剂中质量比的增加(25%~75%), NG的最大引发键长逐渐下降, 这可能是由于体系中NG浓度的降低, 影响了NG与体系中其他分子的相互作用, 同时减少了模型引发键的统计基数。基于同样的原因, TMETN的最大引发键长随着其浓度的增加而逐渐增加。

表 4 各共混模型中NG与TMETN引发键的最大键长(Lmax) Tab.4 The maximum trigger bond lengths of NG and TMETN in blends

在三种含有混合增塑剂的共混体系中, NG与TMETN均有最大引发键长的统计结果, 无法进行各体系间的直接对比。因此选取共混体系中所有增塑剂最大引发键长的最大值, 进行比较并分析共混体系的感度大小, 可以看出使用TMETN部分替代NG与NC混合后, 体系的最大引发键长明显降低, 而随着TMETN含量增加至全部取代, 共混体系的感度有所增加, 这与文献[1]的研究结论一致。当NG:TMETN=3:1时, 共混体系具有最低的最大引发键长, 这意味着NC/NG/TMETN-1可能是本文设计的5种共混体系中感度最低的一个。

已有研究表明, 强氢键作用可能使得炸药引发键键长变短, 感度降低[20]。3.3节径向分布函数的模拟结果显示, 共混体系中NC/NG/TMETN-1的增塑剂与NC具有最强的氢键作用, 而最大引发键长分析可知该体系也具有最低的感度。可见, 在NC/NG/TMETN共混体系中, 强的氢键可以对增塑剂的引发键产生稳定作用, 从而降低共混体系模型的感度。

4 结论

(1) 构建了含有TMETN、NG以及不同质量比TMETN/NG混合增塑剂的5种共混模型, 均具有良好的混溶性。随着TMETN的逐步取代NG, 共混体系的硝基浓度下降, 结合能呈降低趋势。径向分布函数结果表明当NG:TMETN质量比为3:1时, 增塑剂的—NO2基团与NC的—OH基团可以良好匹配, 使两种增塑剂同时与NC产生较强的氢键作用。

(2) 采用混合增塑剂的共混体系的力学性能显著优于单增塑剂体系。当NG:TMETN质量比为3:1时共混体系的综合力学性能最佳。同时该模型组分间强的氢键作用有效的稳定了增塑剂中O—NO2引发键, 降低了共混体系的感度, 对复合含能体系组分间相互作用对于模型感度的影响机理做出了解释。

(3) 硝酸酯增塑剂与NC粘合剂共混后产生的氢键作用可能对体系的力学性能及感度性能产生直接影响。采用混合增塑剂并调节增塑剂配比, 可以增强组分间氢键作用, 优化共混体系的性能。

参考文献
[1]
刘栓虎, 程根旺, 骆广梁. 高能推进剂钝感含能材料研究现状[J]. 化学推进剂与高分子材料, 2010, 8(2): 5-10.
LIU Shuan-hu, CHENG Gen-wang, LUO Guang-liang. Research status of insensitive energetic materials used in high-energy propellant[J]. Chemical Propellants & Polymeric Materials, 2010, 8(2): 5-10.
[2]
Zeman S, Jungová M. Sensitivity and performance of energetic materials[J]. Propellants, Explosives, Pyrotechnics, 2016, 41(3): 426-451. DOI:10.1002/prep.201500351
[3]
付小龙, 樊学忠. 钝感推进剂配方研究及发展趋势[J]. 火炸药学报, 2014, 37(5): 1-8.
FU Xiao-long, Fan Xue-zhong. Research and development trends of insensitive solid propellant formulation[J]. Chinese Journal of Explosives & Propellants, 2014, 37(5): 1-8.
[4]
刘长波, 朱天兵, 马英华, 等. BDNPF/A增塑剂的性能及其应用[J]. 化学推进剂与高分子材料, 2010, 8(1): 23-27.
LIU Chang-bo, ZHU Tian-bing, MA Ying-hua, et al. Properties and its application of BDNPF/A plasticizer[J]. Chemical Propellants & Polymeric Materials, 2010, 8(1): 23-27.
[5]
徐司雨, 仪建华, 贾永杰, 等. 动态气氛下混合酯发射药的热分解性能[J]. 火炸药学报, 2007, 30(4): 70-72.
XU Si-yu, YI Jian-hua, JIA Yong-jie, et al. Thermal decomposition behavior of mixed nitric ester gun propellant under dynamic atmosphere[J]. Chinese Journal of Explosives & Propellants, 2007, 30(4): 70-72.
[6]
何利明, 何伟, 罗运军, 等. 几种钝感含能增塑剂之间的相容性[J]. 固体火箭技术, 2015, 38(4): 523-527.
HE Li-ming, HE Wei, LUO Yun-jun, et al. Compatibility of some insensitive energetic plasticisers with each other[J]. Chinese Journal of Explosives & Propellants, 2015, 38(4): 523-527.
[7]
Gore G M, Tipare K R, Divekar C N, et al. Studies on effect of incorporation of BDNPF/A on burning rates of RDX/AP/AI filled CMDB propellants[J]. Journal of Energetic Materials, 2002, 20(3): 255-278. DOI:10.1080/07370650208244824
[8]
赵凤起, 李上文, 宋洪昌, 等. 国外新型钝感双基推进剂的研究[J]. 飞航导弹, 1999(9): 29-32.
ZHAO Feng-qi, LI Shang-wen, SONG Hong-chang, et al. Research on new insensitive double base propellant abroad[J]. Winged Missiles Journal, 1999(9): 29-32.
[9]
Nair U R, Asthana S N, Rao A S, et al. Advances in high energy materials[J]. Defence Science Journal, 2010, 60(2): 137-151. DOI:10.14429/dsj
[10]
张超, 秦能, 李宏岩, 等. 改性双基推进剂冲击波感度研究[J]. 火工品, 2015(2): 14-17.
ZHANG Chao, QIN Neng, LI Hong-yan, et al. Study on the shock wave sensitivity of modified double base solid propellants[J]. Initiators & Pyrotechnics, 2015(2): 14-17.
[11]
赵凤起, 徐司雨, 仪建华, 等. 含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(Hanneng Cailiao), 2006, 14(6): 406-410.
[12]
赵凤起, 杨栋, 李上文, 等. 以NC和TMETN为基的微烟推进剂机械感度研究[J]. 火炸药学报, 1999, 22(4): 8-4.
ZHAO Feng-qi, YANG Dong, LI Shang-wen, et al. Impact and friction sensitivity of minimum smoke propellant based on NC and TMETN[J]. Chinese Journal of Explosives & Propellants, 1999, 22(4): 8-4.
[13]
朱伟, 刘冬梅, 肖继军, 等. 多组分高能复合体系的感度判据、热膨胀和力学性能的MD研究[J]. 含能材料, 2014, 12(5): 582-587.
ZHU Wei, LIU Dong-mei, XIAO Ji-jun, et al. Molecular dynamics study on sensitivity criterion, thermal expansion and mechanical properties of multi-component high energy systems[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2014, 12(5): 582-587.
[14]
孟玲玲, 齐晓飞, 王江宁, 等. 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, 2015, 38(3): 86-89.
[15]
齐晓飞, 严启龙, 刘萌, 等. DIANP对NC溶塑作用的实验与模拟[J]. 含能材料, 2016, 24(3): 269-273.
QI Xiao-fei, YAN Qi-long, LIU Meng, et al. Experiment and simulation on plastication of NC by DIANP[J]. Chinese Journal of Energetic Materials(Hanneng Cailiao), 2016, 24(3): 269-273. DOI:10.11943/j.issn.1006-9941.2016.03.010
[16]
齐晓飞, 张晓宏, 李吉祯, 等. 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.
[17]
王艳群, 王宝山, 李伟. 硝化甘油在α-Al2O3(0001)和γ-Al2O3(110)表面吸附的理论计算[J]. 高等学校化学学报, 2017, 38(8): 1383-1389.
WANG Yan-qun, WANG Bao-shan, LI Wei. Theoretical investigation on the adsorption of nitroglycerin on the α-Al2O3(0001) and γ-Al2O3(110) surfaces[J]. Chemical Journal of Chinese Universities, 2017, 38(8): 1383-1389. DOI:10.7503/cjcu20170228
[18]
FU X L, FAN X Z, JU X H, et al. Molecular dynamic simulations on the interaction between an HTPE polymer and energetic plasticizers in a solid propellant[J]. Rsc Advances, 2015, 65(5): 52844-52851.
[19]
Xenides D, Randole B R, RODE B M. Hydrogen bonding in liquid water: An ab initio QM/MM MD simulation study[J]. Journal of Molecular Liquids, 2006, 123(2): 61-7.
[20]
李满, 王艳红, 黄红英, 等. 氢键和分子-离子相互作用对奥克托今(HMX)引发键影响的理论研究[J]. 原子与分子物理学报, 2016, 33(1): 46-52.
LI Man, WANG Yan-hong, HUANG Hong-ying, et al. A theoretical study on effect of hydrogen bond and molecular-ion interaction on Octogen(HMX) trigger bond[J]. Journal of Atomic and Molecular Physics, 2016, 33(1): 46-52.
图文摘要

Taking nitrocellulose (NC) as binder, five kinds of blending models which contain 2-methyl-2-[(nitrooxy)methyl]propane-1, 3-diyl dinitrate(TMETN), nitroglycerin(NG), and three kinds of TMETN/NG mixed-plasticizers with different mass ratio were constructed. The interactions, mechanical properties, radial distribution functions, and maximum trigger bond (O—NO2) length (Lmax)etc. parameters for each blending system were simulated by molecular dynamics method using materials studio software.