Effect of Near-Field Overpressure Enhancement of Reactive Material on Low Collateral Damage Ammunition
-
摘要: 为了实现低附带毁伤弹药的近场爆炸威力增强效应,提出在分装式低附带毁伤弹药的重金属颗粒嵌层中加入活性元成分,以期增强近场超压与比冲量。开展了不同含量活性元的静爆实验,利用自由场压力测试系统测得爆炸后近场及中远场的冲击波压力曲线。结果表明:在重金属颗粒嵌层中加入一定含量的活性元后,冲击波超压峰值和比冲量在37.5倍装药直径处分别提高31.6%和21.3%。根据实验结果,利用数值模拟确定了Miller反应速率模型参数,讨论了活性元后燃反应能量释放规律以及活性元组分反应度随时间的变化关系,在充分燃烧的理想情况下,活性元二次燃烧持续时间可达300 ms,且活性元含量区间极有可能存在最优配比。研究结果可为分装式低附带毁伤武器的近场冲击波区域增强效应及其工程化设计提供参考。Abstract: To realize the effect of the near-field shock wave enhancement of low collateral damage ammunition, this study proposes to add reactive material into the embedding layer of heavy metal particle in the sub-package low collateral damage ammunition to enhance the near-field overpressure and specific impulse. Static explosion experiments with different contents of reactive material were carried out, and the pressure curves of the near field and the mid-far field shock wave after the explosion were measured by the free field pressure test system. The results show that: after adding a certain content of reactive material into the heavy metal particle intercalation, the peak value of the overpressure and the specific impulse of shock wave at 37.5 times charge diameter are increased by 31.6% and 21.3%, respectively. According to the experimental results, the parameters of the Miller reaction rate model were determined by numerical simulation, and the law of energy release after-combustion reaction of the active element and the time-dependent change of the reactivity of the reactive material components were discussed. The duration of the secondary combustion can reach 300 ms under the ideal situation of sufficient combustion, and an optimal proportion of the range of reactive material content is likely exist. This research provides considerable development direction and application prospects for the regional enhancement effect of near-field shock wave and its engineering design of sub-packaged low collateral damage weapons.
-
Mn+1AXn相材料是三元层状化合物(简称MAX相材料),其中:M代表Ti、V、Zr等过渡金属元素;A代表A组元素;X代表C或者N;n=1, 2, 3, ···[1]。MAX相这一概念最早由Barsoum[2]提出,这类材料普遍具有陶瓷材料和金属材料的双重特性,可在高压、高温、强腐蚀等极端条件下稳定存在[3-4],并表现出较好的稳定性和抗氧化性,具有极其重要的研究价值和广阔的发展前景[5],因此探究高压等极端状态下的晶体性质变化具有重要意义。
近年来,有关三元层状Mn+1AXn相材料的研究很多,主要集中在211相、312相和413相[6]。随着研究的深入,Mn+1A2Xn双“A”层221相、322相等结构被陆续得到,第一种双“A”层MAX相化合物Mo2Ga2C由Hu等[7]于2015年成功制备,2016年Thore等[8]通过第一性原理计算预测出V2Ga2C的存在,V2Ga2C、Ti3Au2C2等双“A”型MAX相的理论预测和实验制备极大地丰富了MAX族化合物[9]。V2Ga2C 是典型的由理论预测得到的新型双“A”层MAX相材料, Thore等[10]根据声子谱没有虚频判定V2Ga2C具有稳定结构,研究发现这类双“A”层的MAX相材料普遍具有较强的金属性,如更好的机械延展性、易于加工等[11]。目前,常压下V2Ga2C的研究日趋丰富,受限于实验条件的复杂性,高压下V2Ga2C的结构、电子、弹性等性能研究较为困难,为此基于密度泛函理论的第一性原理计算能够很好地解决这一问题。
本研究通过第一性原理对V2Ga2C六方结构的能带结构、态密度等电子结构和弹性性能等力学性质进行计算,根据玻恩稳定准则等相关理论,预测高压状态下V2Ga2C结构的力学稳定性,并对高压下V2Ga2C的晶体结构、电子结构和弹性性质等进行分析,为新型双“A”型Mn+1A2Xn相的相关研究提供理论参考。
1. 计算方法与参数
采用第一性原理计算方法,运用基于密度泛函理论的Materials Studio软件中的CASTEP量子力学程序[12-13],选用倒易点阵空间表征的Cepeley-Alder超软赝势[14]。利用总能量的平面波赝势替代离子势,并通过广义梯度近似(Generalized gradient approximation, GGA)中的PBE(Perdew, Burke and Ernzerhof)[15-16]方法对电子间的相互作用和相关势进行校正。为确保总能量和原子间的作用力最小化,采用Broyden-Fletcher-Goldfarb-Shanno(BFGS)算法,布里渊区K点网格数为17 × 17 × 3,平面波截断能选取550 eV。进行原胞的几何优化(Geometry)时,能量收敛标准为5×10−6 eV/atom,最大作用力为0.01 eV/Å,应力偏差小于0.02 GPa,自洽场收敛精度为5×10−7 eV/atom。
2. 结果与讨论
2.1 晶体结构
V2Ga2C为六方晶系,空间群是P63/mmc,每个晶胞有10个原子,晶体结构见图1。V2Ga2C晶体与常见的V2GaC的结构和性质类似,晶胞键角
α =β =90°,γ =120°,不同的是V2Ga2C有双Ga层结构及不同的晶胞键长a(a=b)、c,V2Ga2C晶胞中各原子坐标为V(1/3, 2/3, 0.0645)、Ga(1/3, 2/3, 0.6814)、C(0, 0, 0)。经优化计算,得到V2Ga2C的晶胞参数为a=b=2.950 Å,c=17.807 Å,与Thore等[8]计算得到的数据(a=b=3.064 Å,c=18.153 Å)基本一致,即本研究构建的模型是准确可行的。为了研究高压对V2Ga2C晶胞结构的影响,在0~70 GPa压强范围内以10 GPa为间隔进行结构优化,得到V2Ga2C晶胞的相对晶格参数变化情况,见图2。从图2可以看出,随着压强增大,晶格常数a、c和体积V均有不同程度的减小,同时在压强范围内V2Ga2C晶胞表现出较好的可压缩性,其中相对键长比a/a0和c/c0从1逐渐减小到0.9019和0.9331,相对晶格参数c/a从0 GPa的6.0362上升到70 GPa的6.2455,c轴较a轴随压强增大收缩得较慢,且键长的减小导致了晶胞体积V的缩小,上述晶胞参数的变化均体现了V2Ga2C的各向异性。此外,根据计算得到的V2Ga2C晶胞在不同压力下的晶格参数及相对晶格常数a/a0、c/c0、c/a和相对晶胞体积V/V0的变化趋势平缓,可判定在0~70 GPa压力范围内V2Ga2C很难发生相变,即本研究利用图1的V2Ga2C结构探究压力对其电子性质、弹性性质的影响是合理准确的。
2.2 高压下的力学稳定性
力学稳定性是晶体材料稳定存在的重要因素。为研究压强对V2Ga2C晶胞力学稳定性的影响,从0 GPa开始,以每10 GPa为一个间隔进行结构优化,通过不同压强下的弹性常数预测V2Ga2C晶胞的力学稳定性。通过各个压强状态下V2Ga2C晶胞的结构优化,得到0~80 GPa不同压强状态下的晶体结构,各压强状态下的弹性常数见表1。
表 1 不同压强下V2Ga2C的弹性常数Table 1. Pressure dependences of elastic constants for V2Ga2CPressure/GPa C11/GPa C33/GPa C44/GPa C12/GPa C13/GPa 0 275.07 309.04 88.92 65.69 48.31 10 325.89 437.85 103.70 79.38 101.23 20 392.18 492.57 112.70 113.70 122.66 30 475.27 582.90 99.80 176.94 174.81 40 462.17 636.80 103.00 156.69 180.10 50 524.87 696.25 51.41 228.45 221.66 60 558.67 743.34 29.53 220.05 219.34 70 640.36 854.90 4.53 280.99 275.36 80 618.66 892.67 −84.84 305.13 305.24 弹性常数是晶体对作用力反应最直观的数据体现。根据V2Ga2C的晶胞结构,V2Ga2C晶体的弹性常数具有对称性,即C11= C22,C13=C31=C32=C23,C12=C21,C44=C55。由表1数据可知,弹性常数C11、C33随压强的增大逐渐增大,C44先增大后逐渐变小直至减小到负数,C12、C13和C66也有不同程度的增大。V2Ga2C晶胞为六方晶系,因此可以通过玻恩稳定准则[17]、正交系统的力学稳定性公式[18]以及弹性常数的变化规律,预测V2Ga2C六方三元层状化合物的力学稳定性。
玻恩稳定准则可写为
C12>0,C33>0,C44>0,C11−C12>0,(C11+C12)C33−2C213>0 (1) 验证V2Ga2C晶胞正交系统力学稳定性的公式为
Cij>0(i=j,0⩽i⩽6),C11+C22−2C12>0,C11+C33−2C13>0,C22+C33−2C23>0,C11+C22+C33+2C12+2C13+2C23>0 (2) 将表1的弹性常数代入式(1)、式(2),可知V2Ga2C晶胞在0~70 GPa符合式(1),在80 GPa时不符合式(1)。因此六方V2Ga2C晶胞的的弹性常数在0~70 GPa压强范围内处于力学稳定状态,80 GPa下V2Ga2C晶胞结构不稳定。
2.3 弹性性质
为了研究压强对V2Ga2C晶体弹性性质的影响,在不同的压强下对晶胞结构进行优化,在此基础上计算不同压强状态下的弹性常数(见表1)。弹性常数C11、C22和C33分别表示晶胞受压沿a、b和c轴的线性压缩阻力,C11、C22较小而C33最大,说明V2Ga2C在a、b轴上容易压缩,在c轴上难压缩;弹性常数C44、C55和C66与材料抗剪切变形能力有关,C44还与硬度有关,随压强增大而减小的C44表明V2Ga2C材料抵抗形变的能力一般。
根据Voigt-Reuss-Hill近似理论[19],V2Ga2C的体积弹性模量B的最大值BV、最小值BR和平均值BH,以及剪切弹性模量G的最大值GV、最小值GR和平均值GH可以通过式(3)~式(8)得到
BV=2(C11+C12)+C33+4C139 (3) BR=(C11+C12)C33−2C213C11+C12+2C33−4C13 (4) BH=BV+BR2 (5) GV=C11+C12+2C33−4C13+12C55+12C663 (6) GR=52[(C11+C12)C33−2C213]C55C663BVC55C66+[(C11+C12)C33−2C213](C55+C66) (7) GH=GV+GR2 (8) 根据Pugh准则[17]可以鉴别晶体的韧脆性,BH/GH < 1.74为脆性材料,相反为韧性材料。根据表1弹性常数和式(3)~式(8),可以计算体积模量BH和剪切模量GH,得到压强与BH/GH的关系曲线,见图3。从图3可以明显看出,压强小于20.15 GPa时,V2Ga2C表现为脆性材料,压强为20.15~70.00 GPa时表现为韧性材料。此外,通过式(9)可以预测维氏硬度(HV)的变化(其中K = GH/BH,HV的单位为GPa),维氏硬度随压强的变化见图4。从图4中曲线的变化趋势可以看出,维氏硬度随着压强的增大逐渐变小,原因是V2Ga2C晶胞的键长和a、b轴随压强的增大急剧压缩,故维氏硬度随之减小。
HV=0.92K1.137G0.708H (9) 然而遗憾的是,目前公开发表的有关V2Ga2C在高压状态下的力学性能研究报道较少,难以与本计算得到的理论预测进行对比分析。
2.4 电子性质
为探究压强与V2Ga2C电子性质的关系,在0~70 GPa的压强范围内通过GGA-PBE密度泛函理论计算,得到V2Ga2C的能带结构图、电子总态密度图,其中0 eV处的虚线表示费米能。
下面以0、35和70 GPa的能带结构为例进行分析,如图5所示。从图5可以明显看出,0 GPa下V2Ga2C无带隙,35 GPa下仍无带隙,直到接近力学稳定临界状态的70 GPa下仍未产生带隙,总体上能带曲线仅有很小幅度的变化。由此可知,在力学稳定范围内,V2Ga2C均无带隙,且压强的增加对能带结构的影响很小,即V2Ga2C材料为导体材料且压强对其影响较小或几乎没有影响。
电子态密度也是V2Ga2C电子性质的重要组成部分,选取0、35和70 GPa状态下的电子态密度分析压强与电子态密度的关系,见图6。由图6可知,随着压强增大,V2Ga2C的总态密度在费米能级附近变动较小,对电子性质影响较小。
3. 结 论
基于密度泛函理论的第一性原理,研究了压强对V2Ga2C晶体的力学稳定性及压强对V2Ga2C结构、弹性和电子性质的影响。根据玻恩稳定准则预测了V2Ga2C结构稳定存在的压强区间为0~70 GPa,并通过正交系统的力学稳定公式验证结果可靠。同时,研究了0~70 GPa压强下V2Ga2C的晶体结构、弹性性质与电子结构,压强使V2Ga2C压缩,体积、相对晶胞参数a/a0、c/c0等均有不同程度的减小,都体现了V2Ga2C具有各向异性;随着压强增大,通过弹性常数可知V2Ga2C在a、b轴上较c轴易压缩,且在20.15 GPa时从韧性转变为脆性,其硬度也随之变小;从V2Ga2C在各个压强状态下的态密度和能带结构可知,在力学稳定的条件下压强对V2Ga2C材料的电子性质影响不大。然而目前有关V2Ga2C材料的研究较少,希望本研究结果可以为双“A”型MAX相材料的实验制备和理论研究提供参考。
-
表 1 活性重金属颗粒嵌层配方
Table 1. Embedded formulation of reactive heavy metal particle
Serial number Mass fraction of heavy metal particle embedded/% WC diameter/µm WC RM Additive LCD-1 90 0 10 150–250 LCD-2 82 10 8 150–250 表 2 JWL状态方程参数
Table 2. Parameters of JWL equation of state
Explosives and RM A/GPa B/GPa R1 R2 ω Q/(kJ·g−1) E/(kJ·g−1) Composition B/0% Al 524.63 7.678 4.2 1.1 0.34 0 4.95 Composition B/10% Al 524.63 7.678 4.2 1.1 0.34 15.67 4.95 表 3 冲击波超压峰值数值模拟结果与实验对比
Table 3. Comparison of numerical simulation results of the peak overpressure of shock wave with experiments
R/m Δp/kPa δ/% Δpe Δps 1.5 169.2 163.6 3.3 2.0 86.0 89.5 4.4 3.0 40.0 42.8 3.3 -
[1] SIRAK M. Air Force focuses on new low-collateral-damage warhead for small bomb [N]. Defense Daily, 2006, 230(10): 1. [2] 刘意. 高密度惰性金属炸药爆轰与粒子流形成过程研究 [D]. 北京: 北京理工大学, 2015: 7–37.LIU Y. Detonation of dense inert metal explosive and the formation of particles flow [D]. Beijing: Beijing Institute of Technology, 2015: 7–37. [3] 申超. 重金属粉末嵌层CFRP壳体内爆下低附带毁伤特性表征 [D]. 北京: 北京理工大学, 2015: 9–22.SHEN C. The low collateral damage characterization of CFRP shell structure with heavy mental powder embedded as a layer under implosion [D]. Beijing: Beijing Institute of Technology, 2015: 9–22. [4] FROST D L, ORNTHANALAI C, ZAREI Z, et al. Particle momentum effects from the detonation of heterogeneous explosives [J]. Journal of Applied Physics, 2007, 101(11): 113529. doi: 10.1063/1.2743912 [5] FROST D L, GRÉGOIRE Y, PETEL O, et al. Particle jet formation during explosive dispersal of solid particles [J]. Physics of Fluids, 2012, 24(9): 091109. doi: 10.1063/1.4751876 [6] XUE K, YU Q Q, BAI C H. Dual fragmentation modes of the explosively dispersed granular materials [J]. The European Physical Journal E, 2014, 37(9): 88. doi: 10.1140/epje/i2014-14088-y [7] 黄德雨, 王坚茹, 陈智刚, 等. 炸药配比对陶瓷低附带毁伤战斗部能量输出的影响 [J]. 弹箭与制导学报, 2012, 32(1): 108–110, 130. doi: 10.3969/j.issn.1673-9728.2012.01.032HUANG D Y, WANG J R, CHEN Z G, et al. Influence of explosive radio on energy output of ceramic low collateral damage warhead [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2012, 32(1): 108–110, 130. doi: 10.3969/j.issn.1673-9728.2012.01.032 [8] 霍奕宇, 王坚茹. 某新型破片材料在低附带毁伤弹药中的应用研究 [J]. 科学技术与工程, 2015, 15(35): 175–178. doi: 10.3969/j.issn.1671-1815.2015.35.031HUO Y Y, WANG J R. Application of the research on the low collateral damage ammunition about new materials fragment [J]. Science Technology and Engineering, 2015, 15(35): 175–178. doi: 10.3969/j.issn.1671-1815.2015.35.031 [9] 霍奕宇, 王坚茹, 陈智刚, 等. 碳纤维壳体壁厚对陶瓷球初速及性能的影响 [J]. 爆破器材, 2016, 45(1): 30–33, 38. doi: 10.3969/j.issn.1001-8352.2016.01.007HUO Y Y, WANG J R, CHEN Z G, et al. Influence of thickness of carbon fiber shell on initial velocity and capability of ceramic ball [J]. Explosive Materials, 2016, 45(1): 30–33, 38. doi: 10.3969/j.issn.1001-8352.2016.01.007 [10] 李俊承, 樊壮卿, 梁斌, 等. 一种低附带弹药金属颗粒定向加载技术 [J]. 爆炸与冲击, 2018, 38(4): 869–875. doi: 10.11883/bzycj-2016-0376LI J C, FAN Z Q, LIANG B, et al. Experimental study on directed loading metal particles of low collateral damage ammunition [J]. Explosion and Shock Waves, 2018, 38(4): 869–875. doi: 10.11883/bzycj-2016-0376 [11] 刘俊, 姚文进, 郑宇, 等. 低附带毁伤弹药的炸药/钨粉质量比对钨粉抛撒特性的影响 [J]. 含能材料, 2015, 23(3): 258–264. doi: 10.11943/j.issn.1006-9941.2015.03.011LIU J, YAO W J, ZHENG Y, et al. Effect of explosive/tungsten powder mass ratio for LCD ammunition on dispersal characteristics of tungsten powder [J]. Chinese Journal of Energetic Materials, 2015, 23(3): 258–264. doi: 10.11943/j.issn.1006-9941.2015.03.011 [12] 左腾, 皮爱国. 亚毫米级金属微粒的爆炸驱动: 模型与数值计算 [J]. 兵工学报, 2016, 37(Suppl 2): 165–170.ZUO T, PI A G. Detonation driving of submillimeter-sized metal particles: model and numerical calculation [J]. Acta Armamentarii, 2016, 37(Suppl 2): 165–170. [13] 杨世全, 孙传杰, 钱立新, 等. 非金属壳体低附带战斗部实验与破片飞散分析 [J]. 高压物理学报, 2018, 32(4): 045103. doi: 10.11858/gywlxb.20170573YANG S Q, SUN C J, QIAN L X, et al. Experimentation and fragment flight analysis of low-collateral-damage warhead with nonmetal shell [J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 045103. doi: 10.11858/gywlxb.20170573 [14] 冯吉奎, 皮爱国, 刘源, 等. 爆炸驱动亚毫米级金属颗粒群的飞散特性 [J]. 高压物理学报, 2019, 33(6): 065104. doi: 10.11858/gywlxb.20190741FENG J K, PI A G, LIU Y, et al. Scattering characteristics of sub-millimeter metal particle group driven by explosion [J]. Chinese Journal of High Pressure Physics, 2019, 33(6): 065104. doi: 10.11858/gywlxb.20190741 [15] 仲倩, 王伯良, 邓金榜. 化学活性材料对温压炸药冲击波参数影响规律研究 [C]//中国力学大会2011暨钱学森诞辰100周年纪念大会论文集. 哈尔滨: 中国力学学会, 2011.ZHONG Q, WANG B L, DENG J B. Study on effect rule of chemically active material on blast wave of thermobaric explosive[C]//Proceedings of the Chinese Society of Theoretical and Applied Mechanics. Harbin: Chinese Society of Mechanics, 2011. [16] 王晓峰, 冯晓军. 温压炸药设计原则探讨 [J]. 含能材料, 2016, 24(5): 418–420. doi: 10.11943/j.issn.1006-9941.2016.05.00XWANG X F, FENG X J. Discussion on the design principle for thermobaric explosives [J]. Chinese Journal of Energetic Materials, 2016, 24(5): 418–420. doi: 10.11943/j.issn.1006-9941.2016.05.00X [17] JIANG C, LU G, MAO L, et al. Effects of aluminum content on the energy output characteristics of CL-20-based aluminized explosives in a closed vessel [J]. Shock Waves, 2021, 31(2): 141–151. doi: 10.1007/s00193-021-01001-1 [18] 隋树元, 王树山. 终点效应学 [M]. 北京: 国防工业出版社, 2000.SUI S Y, WANG S S. Terminal effects [M]. Beijing: National Defense Industry Press, 2000. [19] 田少康, 李席, 刘波, 等. 一种RDX基温压炸药的JWL-Miller状态方程研究 [J]. 含能材料, 2017, 25(3): 226–231. doi: 10.11943/j.issn.1006-9941.2017.03.009TIAN S K, LI X, LIU B, et al. Study on JWL-miller equation of state of RDX-based thermobaric explosive [J]. Chinese Journal of Energetic Materials, 2017, 25(3): 226–231. doi: 10.11943/j.issn.1006-9941.2017.03.009 [20] 张宝銔, 张庆明, 黄风雷. 爆轰物理学 [M]. 北京: 兵器工业出版社, 2001.ZHANG B P, ZHANG Q M, HUANG F L. Detonation physics [M]. Beijing: Weapon Industry Press, 2001. [21] MILLER P J, GUIRGUIS R H. Experimental study and model calculations of metal combustion in Al/Ap underwater explosives [J]. MRS Online Proceedings Library, 1992, 296(1): 299–304. doi: 10.1557/PROC-296-299 [22] 张奇, 白春华, 刘庆明, 等. 一次引爆燃料空气炸药及其爆炸效应研究 [J]. 实验力学, 2000, 15(4): 448–453. doi: 10.3969/j.issn.1001-4888.2000.04.014ZHANG Q, BAI C H, LIU Q M, et al. Investigation on single igniting fuel-air explosive and its explosion effect [J]. Journal of Experimental Mechanics, 2000, 15(4): 448–453. doi: 10.3969/j.issn.1001-4888.2000.04.014 [23] 黄菊, 王伯良, 仲倩, 等. 温压炸药能量输出结构的初步研究 [J]. 爆炸与冲击, 2012, 32(2): 164–168. doi: 10.11883/1001-1455(2012)02-0164-05HUANG J, WANG B L, ZHONG Q, et al. A preliminary investigation on energy output structure of a thermobaric explosive [J]. Explosion and Shock Waves , 2012, 32(2): 164–168. doi: 10.11883/1001-1455(2012)02-0164-05 -