含偏心起爆对EFP战斗部飞行特性的影响*

刘 健 峰,龙 源,纪 冲,赵 长 啸,姜 楠

(1.解 放军理工大 学野战工程 学院,江苏 南 京 210007; 2.北京理工大学爆炸科学与技术国家重点实验室,北京100081; 3.武汉军械士官学校,湖北 武汉430075)

含偏心起爆对EFP战斗部飞行特性的影响*

刘 健 峰1,2,龙 源1,2,纪 冲1,2,赵 长 啸3,姜 楠1

(1.解 放军理工大 学野战工程 学院,江苏 南 京 210007; 2.北京理工大学爆炸科学与技术国家重点实验室,北京100081; 3.武汉军械士官学校,湖北 武汉430075)

为研究起爆不对称性对EFP战斗部飞行特性的影响,对不同偏心量下⌀60 mm 弧锥结合罩 EFP战斗部进行飞行弹道实验。实验结果表明:偏心起爆条件下,当相对偏心量小于3.3% 时,EFP在网靶穿孔接近圆形,弹丸飞行稳定;起爆相对偏心量达到6.7% 时,弹丸飞行过程中摆动幅值增大,降低了对目标的打击精度和毁伤效果。利用 LS-DYNA及 CFX非线性动力学有限元程序对不同起爆偏心量下成型 EFP的空气动力学特性进行数值模拟,描述了偏心起爆影响 EFP成型对称性,改变弹丸在飞行过程中流场的分布特征,从而导致弹丸飞行过程中无规则运动的全过程。

爆炸力学;非对称;空气动力;外弹道;偏心起爆;尾翼弹丸

随着近年来对新型反装甲灵巧弹药的研制,爆炸成 形弹丸(explosively formed projectiles,EFP)战斗部受到越来越广泛的重视和研究。目前正在开发或生产的灵巧弹药多是基于EFP战斗部的末敏弹。这种EFP战斗部的终点弹道能力要求达到1000倍装药口径炸距以上仍可有效打击装甲目标的顶甲或侧甲。在实施远距离攻击时,为保证弹丸经过一个较长的空气弹道阶段后仍能以足够小的攻角准确毁伤目标,要求该类战斗部能够成型具有良好空气动力学特性的EFP。由于EFP是作战实时起爆成型,其部件制造和装配过程中的不对称性以及偏心起爆的位置偏置都会对EFP成型造成影响,使实际成型的弹丸总是或多或少存在不对称性。这些不对称性将直接影响EFP外弹道飞行稳定性及着靶精度等:曹兵等[1]通过实验研究了偏心起爆 对 EFP 成型的影 响,得到了起爆偏 心 量 与 EFP 头 部 横 向 剩 余速 度之间的关系;K.Jach 等[2]通过 数值编 码 对 药 型 罩 施 加 偏 心 载 荷,获 得 尾 翼 偏 置 的 EFP;D.J.Brandeis等[3]通 过 结 构 设 计 研 究 了 各 种 不 对 称 形 状 对 爆 炸 成 形 弹 丸 空 气 动 力 性 能 的 影 响 ;C.Berner等[4]研 究 了尾裙高度、数量及其非对称性对EFP空气动力特性及飞行性能的影响,并且用参数表示了在翻转运动缓慢时非对称性外形结构的的优点;P.Rouge等[5]通过实验和数值计算研 究了 装 药 端面中心点起 爆 以及偏心起爆条件下爆轰产物以及成型弹丸内部压力、速度及变形分布特点。但是,对于由偏心起爆引起EFP外形不对称而对弹丸在飞行弹道(从EFP飞行到命中目标这段距离)运动规律及终点弹道散布水平的影响,尚未见有系统的研究。

本文中基于60 mm 弧锥结合罩EFP装药,通过EFP飞行弹道实验研究不同偏心起爆方式下成型EFP外弹道运动规律特性及终点毁伤效应,利用数值模拟分析不同起爆偏心量下成型EFP的空气动力学特性,研究结果为EFP战斗部起爆参数设置提供有益参考。

1 实验设计

EFP战斗部以自行优化设计的⌀60 mm 变壁厚弧锥结合型紫铜罩战斗部为基础,药型罩口径为56 mm,曲率半径为46 mm,罩顶部壁厚为3.5 mm。装药为圆柱形装药,装药高度为60 mm,战斗部装药结构如图1所示。

实验中采用相同装药结构的EFP战斗部:无偏心量时,采用装药端面中心点起爆方式,起爆点为图2所 示 起 爆 点1;偏 心 起 爆 时,在 距 离 装 药 端 面 中 心 点2、4 mm 处 分 别 设 置 偏 心 起 爆 点2和3,起 爆 点 分别为图2所示起爆点2和3。采用偏心量δ与装药直径d的量纲一比值K表示偏心量相对于装药直径的 大 小,即,则3种 偏 心 起 爆 方 案 的 相 对 偏 心 量 分 别 为0、3.3%和6.7%,图2即 为 3 种不同相对起爆偏心量的设计方案示意图。

图1 ⌀60 mm EFP战斗部实验照片图Fig.1 Experimental picture of⌀60 mm EFP warhead

图2不同偏心量设计方案Fig.2 The design of different initiation eccentricity

图3 所示为偏心起爆条件下 EFP飞行弹道实验现场布置:EFP战斗部水平设置,距离地面1.5 m,沿EFP飞行弹道15 m 处开始,按照5 m 间隔共设置15块网靶,用来捕捉EFP飞行过程中弹形和飞行姿态变化。在离战斗部15 m 处和85 m 处分别设置2组铝箔靶测定 EFP飞行中的速度及速度降。距战斗部中心96 m 处设置1 m×1 m×25 mm 的45钢方形靶,从而获得偏心起爆对EFP的立靶密集度以及侵彻威力的影响。实验之前将战斗部置于特制的木支架上,利用瞄准仪瞄准目标、水平测量仪调整靶板姿态,保证战斗部中心、网靶中心和钢靶中心处在同一水平面上。

图3 EFP战斗部偏心起爆飞行弹道实验布置示意图Fig.3 Setup of EFP flight trajectory experiment

2 偏心起爆影响EFP飞行特性和弹道散布实验

EFP飞行弹道实验中预先在网靶上竖直做一条直线,设置网靶时用铅垂线校正将其作为处理EFP飞行姿态的基准,每发弹丸连续穿过15个网靶。弹丸穿过网靶时,部分网靶破裂,未统计破裂靶纸信息。实验中无偏心时,回收有效靶纸12张;偏心量为3.3%时回收到有效靶纸12张;起爆相对偏心量为6.7%时回收到有效靶纸13张。图4所示为3发 EFP在飞行弹道不同距离上部分网靶穿孔照片。

由图4可以看出网靶较好的记录了EFP飞行过程中姿态的变化情况。实验表明,无偏心量时,弹丸在网靶上的穿孔几乎接近圆形,说明弹丸运行稳定;当相对偏心量为3.3%时,弹丸开始有小幅度波动变化,但是很快弹丸的运动就趋向于稳定,网靶穿孔逐渐接近圆形;当相对偏心量为6.7%时,弹丸在网靶上留下的穿孔接近EFP纵向截面形状,扭曲变形的尾翼在网靶上得到了体现,说明成型 EFP不具有对称性,弹丸在飞行过程中发生剧烈运动,外弹道运行稳定性较差。从 EFP穿孔形状还可以发现EFP弹轴不仅与网靶法线之间的夹角θ(攻角)改变,而且靶纸上弹孔长轴与铅垂线之间夹角γ(摆动角)也不断发生变化。由此可以判断EFP在飞行过程中姿态变化很复杂,在攻角变化的同时还以弹道方向为轴发生摆动。

图4 实验所得EFP在飞行弹道不同距离上网靶穿孔照片Fig.4 Experimental net target perforation photo of EFP at different distances on flight trajectory

EFP在网靶上穿孔的形状由EFP的形状和着靶姿态决定,竖直布置的网靶可以捕捉 EFP的外形和运动特点。通过网靶上的穿孔形状,尺寸及穿孔分布可以测定EFP飞行过程中运动学参数的变化情况[6]。图5~6所示为 EFP飞行攻角和摆动角随飞行 距 离 变化曲线图,图中设 弹 丸 头部穿孔向上 攻 角为正,弹丸头部穿孔向下攻角为负,弹轴偏向铅垂线左侧摆动角为正,弹轴偏向铅垂线右侧摆动角为负。

图5 攻角随EFP飞行距离变化曲线Fig.5 Attack angle varied with different distance on flight trajectory

图6 摆动角随EFP飞行距离变化曲线Fig.6 Swinging angle varied with different distance on flight trajectory

从图5~6中可以看出,EFP飞行过程中必然伴随着攻角和摆动角的变化,EFP在网靶上穿孔是弹丸在该网靶处攻角和摆动角的合成角度。完整EFP在飞行过程中攻角呈现非线性周期摆动,摆动幅度随弹丸运行距离增加有明显减小的趋势。当炸高大于50 m 后,相对偏心量小于3.3%时弹丸攻角基本保持在10°以内,弹丸飞行稳定,即 EFP 进入攻角小幅波动的稳定飞行阶段;相对偏心量为6.7%时, EFP的攻角和摆动角摆动幅值较大而且摆动收敛速度较小,弹丸在空气中发生非周期性摆动。弹丸运行周期保持在10-2s量级。

EFP偏心起爆的弹道实验中,3种偏心起爆成型的EFP均有效击穿距战斗部中心96 m 处厚25 mm的方形45钢靶。图7所示为以瞄准靶心为原点建立坐标系,3发EFP分别在钢靶上穿孔位置分布。从图7中可以看出0偏心量时,由于重力,风速等自然因素导致弹着点偏离靶 心,EFP 飞 行96 m时 偏 离 靶 心 的 绝 对 距 离 为0.08 m。随着相对偏心量的增大,弹着点偏离靶心的位置逐渐增大,相对偏心量达到6.7%时弹丸偏离靶心的绝对距离为0.44 m,EFP着靶精度明显降低。3发弹丸击穿钢靶后在靶板上留下的穿孔都接近圆形。弹体穿过钢靶时,对靶板造成冲塞式破坏,从靶板坑壁可明显观察到紫铜材质EFP挤凿靶板留下的痕迹,在靶板背面还观察到由于层裂引起靶板背面破片崩落产生的环形区域。

图7 钢靶上弹丸穿孔分布图Fig.7 Perforation distributions of 45 steel target

3 偏心起爆对EFP成型及飞行特性的数值计算

为从理论上更深刻认识偏心起爆对 EFP飞行特性及终点弹道散布的影响,采用 ANSYS/LS-DYNA有限元软件对前述实验中偏心起爆条件下EFP成型过程进行数值模拟,并通过计算流体力学软件CFX对成型EFP的外流场特性进行数值模拟分析。

3.1 偏心起爆对EFP成型特性的影响

利用Lagrange网格对 EFP战斗部计算模型进行划分。为提高分析效率,在文中取1/2物理模型进行数值计算,三维有限元计算模型如图8所示。

药型罩材料为紫铜,其本构方程选用Johnson-Cook模型,状态方程 为 Mie-Grüneisen方程;主装药为8701炸药,炸药本构关系的描述选用 HIGHEXPLOSIVE BURN 形式,状态方程选用JWL 状态方程;计算中添加*CONTACT_SLIDING_ONLY_PENALTY关键字定义炸药和金属罩之间接触算法,药型罩内部采用自动面面接触,从而可以有效地避免网格之间发生穿透、畸形。计算模型参数取自参考文献[7]。

图9所示为不同相对偏心量下 EFP成型特点。0偏心起爆时,当炸药起爆200μs后,EFP外形基本稳定,弹丸的飞行速度维持在1 500 m/s左右。偏心起爆时,爆轰波形相对于装药轴线产生偏斜,作用于药型罩后造成药型罩的不对称压垮。随着偏心量增大,弹丸头部和尾翼分别向相反的方向发生偏转:头部对称轴发生逆时针方向偏转,偏转角度对偏心量变化不敏感;尾翼对称轴发生顺时针方向偏转,偏心量对尾翼对称轴偏转的影响程度较大。相对偏心量控制在1.7%以内时,弹丸外形基本不发生变化。相对偏心量超过3.3%时,弹丸头部和尾翼的相对偏转角增大速率明显加快。

炸药起爆后,由于偏心量存在,弹丸获得一个垂直于弹轴的横向速度,运动方向沿偏心起爆点指向装药中心。随着偏心量增加,弹丸的横向速度效应越来越明显。图10所示为不同相对偏心量下弹丸横向速度的变化情况。当相对偏心量小于3.3%时,弹丸横向速度随偏心量变化程度较小,当相对偏心量为3.3%时,弹丸获得横向速度仅为6.1 m/s;相 对 偏 心 量超 过 5% 时 弹丸 的 横 向 速 度 增加 速 率 明显 加快,当相对偏心量达到6.7%时,弹丸的横向速度达到17.3 m/s,经过百米 量级中间弹道飞行后弹 丸 的射偏量增大,命中目标的准确性概率降低,这与实验中3种偏心起爆条件下 EFP命中96 m 处45#钢靶的规律特点吻合较好。

图8 1/2有限元计算模型Fig.8 The 1/2 structure of the finite-element model

图9 不同相对起爆偏心量影响EFP成型Fig.9 Influence of different relative initiation eccentricity on EFP forming

图10 不同相对偏心量对弹丸横向速度效应的影响Fig.10 Influence of different relative initiation eccentricity on EFP lateral velocity

3.2 偏心起爆对EFP空气动力特性影响

为进一步研究偏心起爆条件下成型EFP的空气动力学特性,将在LS-DYNA 中成型的 EFP通过计算流体力学软件CFX对弹丸外流场特性进行数值模拟。数 值 计 算 基 本 方 程 为 三 维 Navier-Stokes方程,湍流模型采用SST模型,采用“双时间步”的二阶隐式格式求解非定常过程。气动力计算条件为标准状态,即大气密度为ρ=1.25 kg/m3,大气压力为101 325 Pa,大 气 温 度 288.15 K,此 时 声 速 为340.29 m/s。以 EFP 稳 定 成 型 的 外 形[8]及 气 动 参数作为气动力学参数分析的初始条件。表1为气动力初始计算状态及基本参数,表中S为最大迎风面积,L为EFP长度。

表1 不同工况气动力初始计算状态及基本参数Table 1 The initial calculation state and basic parameters of different aerodynamic conditions

图11所示为利用计算流体力学软件CFX得到不同起爆偏心量下EFP飞行弹道流场分布情况。由数值计算结果可以看出EFP在空气中超音速飞行,将会在弹丸头部产生弓形脱体正激波。此时弹丸头部区域压力相当高,造成弹体阻力增大,速度衰减加快。弧锥结合形药型罩经过爆炸成型后形成前折式尾翼对流经弹体表面气流起到进一步的阻碍作用,在尾翼前端形成了气体的滞止回流区。随着偏心量的增大,EFP头部和尾翼相对扭曲程度增大,弹丸飞行过程中所受气动阻力增加,EFP速度降迅速增大,弹丸的终点毁伤能力下降。

图11 偏心起爆获得EFP的飞行弹道流场分布Fig.11 The distribution of EFP flow field on flight trajectory

从弹丸尾翼流场分布情况可以看出,随着偏心量增大,弹丸尾翼不对称性增加:尾翼附近,气体发生不对称的分离流动,周围流场分布特点由对称性逐渐发生畸变,弹体压心向弹体前部偏移;弹丸上下表面的压强差增大,产生法向力,在弹体上形成一垂直弹轴的力矩,使弹丸在一定攻角范围内做无规则的俯仰运动更加剧烈,增加了弹体本身的不稳定性,不利于弹体飞行。弹体尾部由于激波和回流区相互干扰以及尾翼结构的不对称性使弹尾流场分布不规则。当偏心量逐渐增大时,弹体尾翼流场分布逐渐由圆形向倒三角形变化,流场分布失去对称性,弹丸运动不稳定性增加。

尾翼弹丸在空气中运动的稳定性可以用弹丸抵抗外界干扰以保持自身飞行稳定的稳定储备量来衡量。为了定量的认识偏心起爆对EFP飞行稳定性的影响,引入稳定储备量B来表征不同相对偏心起爆量条件下,EFP在中间弹道飞行阶段的稳定性。稳定储备量B是指弹丸阻力中心与质心位置的相对距离,即

对于尾翼弹丸,要保证其能够在空气中良好的稳定飞行,其稳定 储备量B必须满足:B>15%[9]。

图12所示为不同偏心量条件下成型EFP的稳定储备量,从图中可以看出0偏心量的弹丸具有相当高的的稳定储备量,为30.29%;当相对偏心量小于5%时,弹丸的稳定储备量保持在20%以上,弹丸飞行过程中所受的稳定力矩大于翻转力矩,即弹丸在飞行过程中抵抗外界干扰能力较强,能够稳定保持自身飞行稳定。当相对偏心量达到6.7%时,弹丸稳定储量仅为13.71%,弹丸飞行过程中的质心和压心距离较小,飞行过程中抵抗外界干扰能力迅速降低,弹丸飞行过程中易受外界条件干扰而发生大幅度摆动甚至翻转。EFP在飞行过程中的攻角和摆动角摆动幅值较大而且摆动收敛速度较小,弹丸攻角长久偏离平衡位置,速度降增大。

图12 不同相对偏心量的弹丸的稳定储备量Fig.12 The stabilization storage of projectile with different relative initiation eccentricity

综上所述,一方面偏心起爆使EFP获得垂直于弹轴的横向速度,影响弹道射偏量的大小,降低了对目标的打击精度;另一方面偏心起爆通过影响EFP对称成型,改变弹丸在飞行过程中流场对称性,使弹丸在飞行过程中的不稳定性增加,弹丸速度降增大,减弱了EFP对终点目标的毁伤效果。

4 结 论

本文中利用数值模拟分析不同起爆偏心量下成型EFP的空气动力学特性,得出如下结论:

(1)偏心起爆时,爆轰波对药型罩的不对称压垮导致 EFP头部和尾翼对称轴线相对偏斜,弹丸获得垂直于弹轴的横向速度。随着偏心量的增大,EFP横向速度效应明显,弹道射偏量增大,降低了对目标的打击精度。

(2)EFP飞行弹道实验和数值计算结果表明:相对偏心量小于3.3%时,EFP 在网靶穿孔接近圆形,说明弹丸在此偏心量内产生的不对称变形对 EFP飞行稳定基本不发生影响;相对偏心量位于3.3%~ 6.7%时,弹丸不对称变形 程度增大,弹体周围流场 分布逐渐失去对称性,当相对偏心量达到 6.7% 时, EFP在一路网靶上留下的穿孔接近弹丸纵向截面形状,说明弹丸在飞行过程中发生较大幅值的摆动,弹丸速度衰减增大,降低了 EFP对目标的终点毁伤能力。因此,EFP起爆的相对偏心量应当控制在3.3%以内。

(3)利用偏心起爆条件下成型 EFP的空气动力学数值计算结果和飞行弹道试验结果有效分析了偏心起爆对EFP飞行弹道稳定性的影响,为评价EFP飞行弹道的稳定性提供了一种简单有效的方法,研究结果为EFP战斗部起爆参数设置以及对目标的精确打击提供有益参考。

[1]曹 兵 ,高 森 烈.偏 心 起 爆 对 EFP 成 形 形 态 影 响 的 实 验 研 究[J].弹 道 学 报,1997,9(1):27-30. Cao Bing,Gao Sen-lie.The effect of eccentric initiation on explosive formed projectile shape[J].Journal of Ballistics,1997,9(1):27-30.

[2]Jach K,Mroczkowski M,Sarzynski A,et al.3D free particle computer modeling of explosive formation of projectiles[C]∥Proceedings of the 16th International Symposium on Ballistics.San Francisco,USA:International Ballistics Committee,1996:557-563.

[3]Brandeis D J.Effect of shape and asymmetry on the aerodynamic performance of explosively formed projectiles[C]∥Proceedings of the 13th International Symposium on Ballistics.Stockholm,Sweden:International Ballistics Committee,1992:137-144.

[4]Berner C,Fleck V.Pleat and asymmetry effects on the aerodynamics of explosively formed penetrators[C]∥Proceedings of the 18th International Symposium on Ballistics.San Antonio,USA:International Ballistics Committee, 1999:237-245.

[5]Rouge P,Weimann K.Consequences on EFP formation of an intentionally created default[C]∥Proceedings of the 10th International Symposium on Ballistics,International Ballistics Committee.San Diego,USA,1987:277-288.

[6]翁 佩 英,任 国 民,于 骐 .弹 药 靶 场 试 验[M].北 京 :兵 器 工 业 出 版 社 ,1995:32-42.

[7]周 翔 .爆 炸 成 形 弹 丸 战 斗 部 的 相 关 技 术 研 究[D].南 京:解 放 军 理 工 大 学,2006.

[8]刘 健 峰,龙 源,纪 冲,等.不 同 药 型 罩 壁 厚 组 合 影 响 同 轴 EFP 成 型 规 律 研 究[J].工 程 爆 破,2012,18(4):9-13. Liu Jian-feng,Long Yuan,Ji Chong,et al.Influence of different wall-thickness of arc-cone liner in collinear EFP forming performance[J].Engineering Blasting,2012,18(4):9-13.

[9]魏 惠 之.弹 九 设 计 理 论[M].北 京:国 防 工 业 出 版 社,1985:80-82.

[10]Bender D,Chhouk B.Explosively formed penetrators(EFP)with canted fins[C]∥Proceedings of the 19th International Symposium on Ballistics.Interlaken,Switzerland:International Ballistics Committee,2001:755-762.

[11]赵 长啸,龙源 ,余道 强,等.切 割 式 多 爆 炸 成 形 弹 丸 成 形 及 对 钢 靶 的 穿 甲 效 应 [J].爆 炸 与 冲 击,2013,33(2):186-193. Zhao Chang-xiao,Long Yuan,Yu Dao-qiang,et al.Formation of incised multiple explosively-formed projectiles and their armor-piercing effect against steel targe[J].Explosion and Shock Waves,2003,33(2):186-193.

[12]纪 冲,龙 源,余道强 ,等 .切割 式双模战斗 部毁伤元成 型及侵彻钢 靶特性研究[J].高压物 理学 报,2012,26(10):508-516. Ji Chong,Long Yuan,Yu Dao-qiang,et al.Experimental and numerical study on the formation and penetration properties of dual-mode warhead[J].Chinese Journal of High Pressure Physics,2012,26(10):508-516.

Effect of eccentric initiation on the flight characteristics and ballistic dispersion of EFP

Liu Jian-feng1,2,Long Yuan1,2,Ji Chong1,2,Zhao Chang-xiao3,Jiang Nan1
(1.College of Filed Engineering,PLA University of Science and Technology, Nanjing 210094,Jiangsu,China; 2.State key Laboratory of Explosion Science and Technology,Beijing Institute of Technology, Beijing 100081,China; 3.Wuhan Ordnance NCO Academy of PLA,Wuhan 430075,Hubei,China)

Flight trajectory experiments were conducted with⌀60 mm diameter EFP warhead in order to investigate the effect of eccentric initiation on the flight characteristics of EFP.The results indicate that the projectile flight stability with perforations nearly circular in the net target when eccentricity is less than 2 mm.The accuracy of strike is lower and terminal effects are less harmful because the projectile overturns during the flight,with the eccentric value of 4 mm.The forming of EFP and its flight characteristics with different eccentricity was numerically analyzed by using LS-DYNA and CFX.EFP is asymmetrical under the conditions of the eccentric initiation.It changes the symmetry of the flow field which leads to instability during flight.The research results provide a reference to detonation parameter of EFP warhead.

mechanics of explosion;asymmetrical;aerodynamic;exterior trajectory;eccentric initiation;projectile with stabilizing fins

O383;TJ410.33国标学科代码:13035

:A

10.11883/1001-1455-(2015)03-0335-08

(责任编辑 王易难)

2013-10-29;

2014-05-20

北京理工大学爆炸科学与技术国家重点实验室开放基金项目(KFJJ10-2 M)

刘健 峰(1988— ),男,博士研 究生,18260098162@163.com。