A model for rigid sharp-nosed projectile perforating metallic targets considering free-surface and cracking effects
-
摘要: 将靶体视为不可压缩材料,假定空腔膨胀产生塑性-弹性响应分区,构造了靶背自由表面效应的衰减函数。将衰减函数乘以可压缩幂次硬化材料的阻力方程,得到了弹体贯穿金属靶板的阻力函数。基于弹性衰减-塑性衰减-开裂3个阶段,建立了同时考虑靶体可压缩性、靶背自由表面和开裂影响的弹体贯穿有限厚金属靶板的分析模型,推导得出了弹体瞬时速度的解析方程,并采用数值方法计算得到了弹体的过载、瞬时速度和残余速度。通过与6组实验数据和已有模型的对比得到,当靶板厚度和弹体冲击速度在一定范围内时,需要考虑自由表面效应的影响。Abstract: Treating the target as the incompressible material, by assuming that the cavity expansion produce plastic-elastic response region, the decay function for the back free-surface effect of target is constructed. The forcing function of metallic targets for perforation is obtained by multiplying the forcing function of compressible Strain-Harding targets with the decay function. Based on the three stage of elastic-decaying, plastic-decaying and cracking, the analytical model considering the compressibility, the back free-surface effect of target and cracking is established, and the analytical equation of instantaneous velocity of projectile is deduced. The deceleration, instantaneous and residual velocity of projectile is calculated by numerical methods. Through comparison with six sets of experimental data and other existing models, with the target thickness and impact velocity in a certain range, the free-effect should be considered.
-
多功能含能结构材料(multifunctional energetic structural material, MESM)是综合利用化学能和动能以提高战斗部毁伤效能的新型功能材料[1]。不同于传统惰性破片单一的动能杀伤作用,当MESM破片以一定速度撞击靶板时,会激发剧烈的爆炸/燃烧反应,破片贯穿靶板后,通过动能和化学能的联合作用,实现对靶后目标产生更大的杀伤破坏[2]。
MESM对特定目标的毁伤是现今较活跃的领域之一,研究主要集中在其释能特性和侵彻毁伤效应上。Zhang等[3-4]和Xiong等[5-6]系统研究了Al/Ni基MESM的冲击释能特性及反应机理,并推导了温度控制冲击诱发化学反应的热化学模型[4, 6],结果表明MESM的释能反应特性与破片冲击引发的材料温升有关,反应规律符合Avrami-Erofeev方程。Wang等[7]和Luo等[8]等分别研究了高聚物、W/ZrMESM的冲击释能特性。Xu等[2,9]研究了PTFE/Al/W破片对装甲铝板的侵彻毁伤效应,通过弹道试验结果拟合了PTFE/Al/W破片对装甲铝板的极限穿透速度计算公式,并研究了破片对双层铝板的结构破坏作用。
W/ZrNiAlCu亚稳态合金复合材料是新一代MESM,具有较强的侵彻能力及释能特性,但其对轧制均质装甲(RHA)靶板的侵彻释能特性研究仍属空白。基于此,本文通过弹道枪侵彻实验和高速摄影,测量破片对RHA靶板的撞击速度和靶板贯穿形成的冲塞体速度,记录破片的侵彻释能过程,研究W/ZrNiAlCu亚稳态合金复合材料破片的侵彻释能特性,为该亚稳态合金复合材料的进一步研究和应用提供理论和实验依据。
1. 侵彻实验
1.1 实验方案
实验试样为W/ZrNiAlCu亚稳态合金复合材料破片,材料通过液态真空压力浸渗工艺制备:将熔体ZrNiAlCu亚稳态合金通过高纯氩气吹入由W颗粒粉末烧结而成的基体中,保压0.5 h保证熔体完全浸渗W基体孔隙,随后放入饱和食盐水中淬火,制备出W/ZrNiAlCu亚稳态合金复合材料。通过机加工,得到边长8 mm的立方体破片,破片密度为12.9 g/cm3。
侵彻实验的实验布置如图1所示。破片及弹托由14.5 mm弹道枪发射,弹托刻有凹槽,保证出枪口后弹托、破片分离;10.5 mm厚RHA靶板(GY4)中点位于射击线上,距弹道枪口5 m,一组测速靶置于RHA靶板前方,距靶板1 m,用以测量破片初速,另一组测速靶置于RHA靶板后方,距靶板1 m,用以测量靶板剪切形成的冲塞体速度,每组测速靶由相距0.5 m的两块断通靶纸组成,断通靶纸与计时器相连接;高速摄影机置于靶板侧方,用以记录W/ZrNiAlCu亚稳态合金复合材料破片对RHA靶板的侵彻释能过程。
1.2 实验结果
图2为4个典型撞击速度vi下W/ZrNiAlCu亚稳态合金复合材料破片对RHA靶板的侵彻释能过程,可以看出,破片撞击速度对破片侵彻过程中的释能现象有较大影响。图2(a)中,由于撞击速度较低,破片未贯穿靶板。破片撞击靶板引发了材料的释能反应,靶板前方可看到2束对称的火光并逐渐向外扩散。图2(b)~2(d)中,随着撞击速度增加,靶板前方光束变长,火光范围增大、亮度变高。图2(c)和图2(d)中均可观察到4道光束。当立方体破片撞击靶板时,破片的棱角处应力集中导致材料破碎并向外飞散,因此侧向观察可以看到2~4束对称的火光,火光由材料碎片燃烧产生,并伴随破片飞散而逐渐扩散。当破片撞击速度增大时,撞击能增加,碎片数量增多、初始动能变大、材料温度更高、释能反应更充分,因此靶板前方火光范围增大、亮度变高。
当撞击速度vi=1 067或1 259 m/s时(图2(b)和2(c)),靶板后方呈现一道圆柱型光束,光束顶端可以看到靶板剪切破坏形成的冲塞体。在撞击并贯穿靶板过程中破片破碎,部分碎片穿过靶板并继续飞行,破片碎片由于初始速度、质量、形状的差异,形成碎片云并呈圆柱状分布。当vi=1 508 m/s时(图2(d)),靶板后方圆柱型光束的末端出现椭球状光亮,并在圆柱型光柱消散后继续燃烧(t=3.75 ms)。当破片撞击速度较大时,部分破片碎片温度较高,贯穿靶板并与空气混合后发生燃烧,形成椭球状碎片云,该部分破片温度高,反应充分,燃烧时间较长,在靶后方形成了持续一定时间的椭球状光亮。在一定撞击速度下,W/ZrNiAlCu亚稳态合金复合材料破片既可有效贯穿10.5 mm厚的RHA靶板,又可在靶后发生明显的燃烧反应,是一种具有较强侵彻能力的多功能含能结构材料。需要注意的是,由于侵彻过程中破片的破碎不可忽略,不能使用传统惰性破片侵彻公式解释W/ZrNiAlCu亚稳态合金复合材料破片的侵彻规律。
表1为不同撞击速度vi下W/ZrNiAlCu亚稳态合金复合材料破片对RHA靶板的侵彻结果,m为破片质量,vp为测得的冲塞体速度,p为根据连续介质力学公式计算得到的破片理论冲击压力[10]。可以看到,vi=973 m/s时:破片贯穿靶板;而vi=995 m/s时,破片未贯穿靶板。考虑到破片质量差别及弹道极限速度vbl附近数据量较小,可以认为W/ZrNiAlCu亚稳态合金复合材料破片对10.5 mm厚RHA靶板的弹道极限速度vbl处于810~1 067 m/s范围内。
表 1 破片侵彻实验结果Table 1. Experimental results of fragments penetration编号 m/g vi/(m·s−1) vp/(m·s−1) p/GPa 是否穿透 1# 6.65 726 − 16.46 否 2# 6.54 810 − 18.63 否 3# 6.39 995 − 23.62 否 4# 6.68 973 34 23.01 是 5# 6.80 1 067 137 25.63 是 6# 6.72 1 144 176 27.83 是 7# 6.37 1 159 131 28.26 是 8# 6.65 1 240 195 30.64 是 9# 6.51 1 295 317 32.28 是 10# 6.56 1 359 397 34.21 是 11# 6.68 1 508 497 38.85 是 2. 理论分析
2.1 弹道极限
破片侵彻靶板过程中加速度历程较为复杂,因此采用质量守恒和能量守恒为基础的整体法求解弹道极限速度,做以下三点假设:
(1)破片首先与自由冲塞体发生塑性碰撞,此过程不考虑靶板的绝热剪切效应[11];
(2)破片撞击靶板后部分破碎,靶前飞散碎片不参与侵彻过程;
(3)冲塞体速度vp约为靶后碎片飞出靶板瞬间统计平均速度vr的1.35倍[11]。
破片贯穿靶板前后的能量守恒方程为:
12mv2i=12mf(mm+mpvi)2+12mbv2r+12mpv2p+Efn+Wp (1) 式中:mf为靶板前方飞溅的破片碎片质量之和,mb为破片贯穿靶板后残余质量,mp为靶板冲塞体质量,Efn为破片与靶板塑性碰撞损失能量,Wp为冲塞绝热剪切耗能。弹靶组合确定时,mf仅与撞击速度有关,采用二次多项式的形式拟合:
mf=m[a(viv∗bl)2+b(viv∗bl)+c] (2) 式中:a、b、c为待定系数,
v∗bl 为不考虑质量损失时的弹道极限速度,mf满足质量守恒m=mf+mb。式(1)中,Efn、Wp的表达式分别为[11]:Efn=mpm+mp⋅12mv2i (3) Wp=4dh2⋅Y√3 (4) 式中:d为破片边长;h为靶板厚度;Y=1.2 GPa,为靶板失效应力。由弹道极限速度的定义,当vi= vbl时,vp=vr=0,令a=b=c=0,将式(3)和(4)代入式(1)得vbl*=937.8 m/s;将表1实验数据代入式(1)拟合得到参数a= −2.02, b= 6.47, c=−4.41,解式(1)得vbl=987.1 m/s,在实验所确定的弹道极限速度范围内。
图3为W/ZrNiAlCu亚稳态合金复合材料破片侵彻10.5 mm厚RHA靶板的vi-vp曲线,其中蓝色实线为考虑破片质量损失所得曲线,红色实线为未考虑破片质量损失所得曲线,星型为实验数据,可以看到,vp与vi正相关,蓝色实线始终位于红色实线下方。破片侵彻能力与其初始动能有关,故破片撞击速度vi越大,冲塞体速度vp越大,侵彻能力增强。破片在侵彻靶板过程中的质量损失将导致动能降低,进而降低破片的侵彻能力,因此蓝色实线位于红色实线下方,理论弹道极限速度vbl>
v∗bl ,破片质量损失对侵彻能力有较大影响。图上星状点为实验数据,可以看到,理论计算结果与实验数据吻合较好,表明所做假设合理,式(1)可以较好地解释包含质量损失的破片对中厚靶板的侵彻规律。2.2 释能特性
由Grünesien物态方程和Rankin-Hugoniot能量方程得到固体p-V形式的物态方程[6]:
p(V)=Vγ(V)pc(V)−Ec(V)Vγ(V)−12(V0−V) (5) 式中:p为材料压力,V为比容,
γ 为Grünesien系数,pc、Ec分别为材料冷压和冷能。通过Meyer势,金属材料的冷压和冷能可以表示为材料常数Q-q的形式[6]:Ec(δ)=3Qρ0K{1qexp[q(1−δ−1/3)]−δ−1/3−1q+1} (6) pc(δ)=Qδ2/3{exp[q(1−δ−1/3)]−δ2/3} (7) 式中:
δ=ρ/ρ0K=V0K/V ,为材料的相对于0 K时密度的压缩度,ρ0K 、V0K 为0 K时的密度和比容。Dugdale -MacDonald关系式的Q-q形式为[6]:
γ(V)=16⋅q2δ−1/3exp[q(1−δ−1/3)]−6δqexp[q(1−δ−1/3)]−2δ (8) 采用胡金彪等给出的解析法可以求得材料常数Q、q[6]:
s[1+(s2−18γ20Ks−1)αVT0]=112⋅q2+6q−18q−2 (9) {C0[1+(2s−γ20K4−1)αVT0]}2=Q(q−2)3ρ0K (10) 式中:C0、s为Hugoniot参数,其中C0理论上等于材料的零压体积声速;
γ0K 为0 K时的Grünesien系数;αV 为材料的体积膨胀系数;T0为起始温度,即实验室室温,取298 K。Boslough根据固态材料绝热线和等熵线的关系,假定材料等容比热为常数,得到冲击波温度表达式[12]:
TH=T0exp(−∫VHV0VγdV)+VγpH−pSCV (11) 式中:TH为击波温度,pH为击波压力,V0、VH分别为初态、击波压缩状态下的比容,CV为材料的等容比热,pS为与冲击绝热线同一起始状态出发的等熵线压力,其解析式为[12]:
pS(η)=ρ0C20exp(γ0η)∫η0(1+(s−γ0))x(1−sx)3exp(−γ0x)dx (12) 式中:
η=1−V/V0 为压缩度,x为压缩度的积分变量,ρ0 为起始状态下材料密度,γ0 为初始Grünesien系数。将式(5)、式(11)联立,即可得到击波压力与击波温度的关系。假设含能结构材料的化学反应过程仅受击波温度控制,研究表明,含能结构材料在高升温速率下的固态反应动力学满足n维Avrami-Erofeev方程,经过整理,Avrami-Erofeev方程可以表示为击波温度对反应效率一次微分的形式[6]:
dTdy=RuT2Ea[12y−nln(1−y)+n−1n(1−y)[−ln(1−y)]] (13) 式中:y为反应效率,Ea为表观活化能,Ru为理想气体常数,n为与边界条件和反应机理有关的系数。
表2为计算所需W/ZrNiAlCu亚稳态合金复合材料参数,图4为破片冲击压力与反应效率关系的理论曲线,可以看到,在40 GPa冲击压力范围内,材料并未达到完全反应,当p=40 GPa,材料反应效率y=0.541;材料反应效率与冲击压力正相关,因此,随着冲击速度增加,破片冲击释能所产生的火光范围、火光亮度增大,与实验现象相符合。
表 2 W/ZrNiAlCu亚稳态合金复合材料参数[13]Table 2. Parameters of W/ZrNiAlCu metastable reactive alloy composite材料 C0/(m∙s−1) s γ0K αV/(10−5·K−1) CV/(J·g−1·K−1) Ea/(kJ·mol−1) n W/ZrNiAlCu 3 417 1.732 1.629 1.559 0.182 459.25 0.347 3. 结 论
(1)W/ZrNiAlCu亚稳态合金复合材料破片可以贯穿10.5 mm厚的RHA靶板,破片撞击、贯穿靶板过程中导致材料破碎并激发燃烧反应,靶前、靶后均可观察到大范围火光;随着撞击速度增加,火光范围增大,亮度提高,燃烧反应时间变长。
(2)推导了W/ZrNiAlCu亚稳态合金复合材料破片侵彻规律方程,得到了破片侵彻10.5 mm厚RHA靶板的理论弹道极限速度vbl=987.1 m/s,理论计算结果与实验数据吻合较好;破片侵彻过程中的质量损失对其侵彻能力有较大影响,在相同撞击速度下,W/ZrNiAlCu亚稳态合金复合材料破片侵彻形成冲塞体速度小于不考虑质量损失的理想情况下冲塞体速度。
(3)由材料的Grunesien物态方程和Avrami-Erofeev方程得到材料反应效率y与冲击压力p的理论关系,在15~40 GPa范围内,材料反应效率随着冲击压力的增加而增加,当p=40 GPa时,y=0.541,因此实验过程中材料反应并不完全,出现随着撞击速度增加,火光范围增大、亮度提高的现象。
-
表 1 刚性弹贯穿金属靶板实验的弹体参数
Table 1. Experimental projectile parameters in perforation of rigid projectiles into metallic targets
表 2 刚性弹贯穿金属靶板实验的靶体参数
Table 2. Experimental target parameters in perforation of rigid projectiles into metallic targets
No. 文献 弹头形状 靶体材料 H/mm E/GPa ρ/(g·cm-3) Y/MPa γ n 1 [13] 锥头弹 6061-T6铝 25.4 68.9 2.71 276 0.33 0.051 2 [15] 锥头弹 5083-H131铝 12.7 70.3 2.66 276 1/3 0.084 3 [15] 锥头弹 5083-H131铝 50.8 70.3 2.66 276 1/3 0.084 4 [15] 锥头弹 5083-H131铝 76.2 70.3 2.66 276 1/3 0.084 5 [16] 卵头弹 6061-T651铝 26.3 69 2.71 262 0.33 0.085 6 [17] 卵头弹 5083-H116铝 20 71 2.66 240 0.33 0.108 -
[1] 陈小伟.穿甲/侵彻问题的若干工程研究进展[J].力学进展, 2009, 39(3):316-351. doi: 10.3321/j.issn:1000-0992.2009.03.006Chen Xiaowei. Advances in the penetration/perforation of rigid projectiles[J]. Advances in Mechanics, 2009, 39(3):316-351. doi: 10.3321/j.issn:1000-0992.2009.03.006 [2] 肖云凯, 吴昊, 方秦, 等.穿甲弹冲击多层金属靶板终点弹道参数的研究进展[J].防护工程, 2014, 36(3):52-68. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=fhgc201403010Xiao Yunkai, Wu Hao, Fang Qin, et al. Review on the terminal ballistic parameters of the armor-piercing penetrator striking the multi-layered metallic shields[J]. Projective engineering, 2014, 36(3):52-68. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=fhgc201403010 [3] 蒋志刚, 曾首义, 周建平.尖头弹丸撞击下金属靶板弹道极限的两种工程模型[J].应用力学学报, 2005, 22(1):21-25. doi: 10.3969/j.issn.1000-4939.2005.01.005Jiang Zhigang, Zeng Shouyi, Zhou Jianping. The two engineering models of ballistic limits for sharp-nosed projectile striking metallic plates[J]. Chinese Journal of Applied Mechanics, 2005, 22(1):21-25. doi: 10.3969/j.issn.1000-4939.2005.01.005 [4] Forrestal M J, Warren T L. Perforation equations for conical and ogive nose rigid projectiles into aluminum target plates[J]. International Journal of Impact Engineering, 2009, 36(2):220-225. doi: 10.1016/j.ijimpeng.2008.04.005 [5] Chen Xiaowei, Li Qingming. Perforation of a thick plate by rigid projectiles[J]. International Journal of Impact Engineering, 2003, 28(7):743-759. doi: 10.1016/S0734-743X(02)00152-5 [6] Wen Heming. Predicting the penetration and perforation of targets by projectiles at normal incidence[J]. Mechanics of Structures and Machines, 2002, 30(4):543-577. doi: 10.1081/SME-120015076 [7] 吴乔国.不同材料靶板的抗弹性能研究[D].合肥: 中国科学技术大学, 2012. http://cdmd.cnki.com.cn/Article/CDMD-10358-1012503653.htm [8] 孙炜海.锥头弹丸正撞击下金属靶板破坏模式的理论和数值模拟研究[D].合肥: 中国科学技术大学, 2009. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y1498383 [9] 蒋志刚, 曾首义, 周建平.中等厚度金属靶板的三阶段贯穿模型[J].兵工学报, 2007, 28(9):1046-1052. doi: 10.3321/j.issn:1000-1093.2007.09.005Jiang Zhigang, Zeng Shouyi, Zhou Jianping. A three-stage model for the perforation of moderately thick metallic plates[J]. Acta Armamentarii, 2007, 28(9):1046-1052. doi: 10.3321/j.issn:1000-1093.2007.09.005 [10] Luk V K, Forrestal M J, Amos D E. Dynamic spherical cavity expansion of strain-hardening materials[J]. Journal of Applied Mechanics, 1991, 58(1):1-6. doi: 10.1115/1.2897150 [11] Warren T L, Poormon K L. Penetration of 6061-T6511 aluminum targets by ogive-nosed VAR 4340 steel projectiles at oblique angles: Experiments and simulations[J]. International Journal of Impact Engineering, 2001, 25(10):993-1022. doi: 10.1016/S0734-743X(01)00024-0 [12] 孔祥振, 方秦, 吴昊.考虑靶体自由表面和开裂区影响的可变形弹体侵彻脆性材料的终点弹道分析[J].兵工学报, 2014, 35(6):814-821. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bgxb201406010Kong Xiangzhen, Fang Qin, Wu Hao. Ballistics study of the deformable projectile penetrating brittle material targets with the free-surface and crack region effects[J]. Acta Armamentarii, 2014, 35(6):814-821. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bgxb201406010 [13] Rosenberg Z, Forrestal M J. Perforation of aluminum plates with conical-nosed rods-additional data and discussion[J]. Journal of Applied Mechanics, 1988, 55(1):236-238. doi: 10.1115/1.3173639 [14] Wen Heming. Deformation and t earing of clamped work-hardening beams subjected impulsive loading[J]. International Journal of Impact Engineering, 1996, 18(4):425- 433. doi: 10.1016/0734-743X(95)00053-D [15] Forrestal M J, Luk V K, Brar N S. Perforation of aluminum armor plates with conical-nose projectiles[J]. Mechanics of Materials, 1990, 10(1):97-105. doi: 10.1016-0167-6636(90)90020-G/ [16] Piekutowski A J, Forrestal M J, Poormon K L, et al. Perforation of aluminum plates with ogive-nose steel rods at normal and oblique impacts[J]. International Journal of Impact Engineering, 1996, 18(7):877-887. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=416180fc47231590075ba96a78c39250 [17] Børvik T, Forrestal M J, Warren T L. Perforation of 5083-H116 aluminum armor plates with ogive-nose rods and 7.62 mm APM2 bullets[J]. Experimental Mechanics, 2010, 50(7):969-978. doi: 10.1007/s11340-009-9262-5 [18] Forrestal M J, Brar N S, Luk V K. Penetration of strain-hardening targets with rigid spherical-nose rods[J]. Journal of Applied Mechanics, 1991, 58(1):7-10. doi: 10.1115/1.2897183 [19] 周辉.弹塑性材料中的空穴膨胀理论及其在侵彻力学中的应用[D].合肥: 中国科学技术大学, 2004. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y615875 [20] Tucker M T, Horstemeyer M F, Whittington W R, et al. The effect of varying strain rates and stress states on the plasticity, damage, and fracture of aluminum alloys[J]. Mechanics of Materials, 2010, 42(10):895-907. doi: 10.1016/j.mechmat.2010.07.003 [21] 肖云凯, 方秦, 吴昊, 等.刚性尖头弹扩孔贯穿金属靶板理论模型的讨论[J].振动与冲击, 2016, 35(1):195-215. http://d.old.wanfangdata.com.cn/Periodical/zdycj201601031Xiao Yunkai, Fang Qin, Wu Hao, et al. Discussions on the ductile hole expansion perforation models for the rigid sharp-nosed projectile perforating metallic plate[J].Journal of Vibration and Shock, 2016, 35(1):195-215. http://d.old.wanfangdata.com.cn/Periodical/zdycj201601031 [22] Forrestal M J, Frew D J, Hickerson J P, et al. Penetration of concrete targets with deceleration-time measurements[J]. International Journal of Impact Engineering, 2003, 28(5):479-497. doi: 10.1016/S0734-743X(02)00108-2 [23] 陈小伟, 金建明.动能深侵彻弹的力学设计(Ⅱ):弹靶的相关力学分析与实例[J].爆炸与冲击, 2006, 26(1):71-78. doi: 10.3321/j.issn:1001-1455.2006.01.012Chen Xiaowei, Jin Jianming. Mechanics of structural design of EPW (Ⅱ): Analyses on the design of EPW projectiles, concrete targets and examples[J]. Explosion and Shock Waves, 2006, 26(1):71-78. doi: 10.3321/j.issn:1001-1455.2006.01.012 [24] Wu Hao, Fang Qin, Peng Yong, et al. Hard projectile perforation on the monolithic and segmented RC panels with a rear steel liner[J]. International Journal of Impact Engineering, 2014, 76:232-250. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=68d38b07a3a413e4666ef584536d4659 [25] 黄徐利, 陈小伟, 梁冠军.尖头弹穿甲韧性金属靶模型分析[J].爆炸与冲击, 2011, 31(5):490-496. doi: 10.11883/1001-1455(2011)05-0490-07Huang Xuli, Chen Xiaowei, Liang Guanjun. Analysis on perforation of ductile metallic plates by a rigid sharp-nosed projectile[J].Explosion and Shock Waves, 2011, 31(5):490-496. doi: 10.11883/1001-1455(2011)05-0490-07 -