变截面爆炸成型弹丸垂直侵彻装甲钢板靶后破片质量模型

邢柏阳; 刘荣忠; 张东江; 陈亮; 侯云辉; 郭锐

doi:10.11883/bzycj-2018-0187

变截面爆炸成型弹丸垂直侵彻装甲钢板靶后破片质量模型

doi: 10.11883/bzycj-2018-0187

1.
南京理工大学机械工程学院，江苏南京 210094
2.
西安现代控制技术研究所，陕西西安 710065

基金项目: 国家自然科学基金（11372136）；中央高校基本科研业务费专项基金（30916011306）

详细信息

作者简介:
邢柏阳（1992－　），男，博士研究生，xing_boyang@163.com

通讯作者:
郭　锐（1980－　），男，博士，副教授，guorui@njust.edu.cn

中图分类号: O385; TJ012.4; TJ413.+2
计量
- 文章访问数: 5773
- HTML全文浏览量: 2186
- PDF下载量: 84
- 被引次数: 6
出版历程
- 收稿日期: 2018-05-29
- 修回日期: 2018-06-28
- 网络出版日期: 2019-06-25
- 刊出日期: 2019-07-01

A mass model for behind-armor debris generated by normal penetration of a variable cross-section explosively-formed projectile into an armor steel plate

1.
School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu, China
2.
Xi’an Modern Control Technology Research Institute, Xi’an 710065, Shaanxi, China

摘要

摘要: 考虑爆炸成型弹丸(explosively-formed projectile，EFP)变截面的特性，基于流体力学Bernoulli方程和绝热剪切理论，改进了EFP垂直侵彻装甲钢板靶后破片质量模型，结合已有的试验数据和数值仿真方法检验了改进后模型的准确性。在此基础上，分析了靶板厚度和EFP着靶速度对靶板和EFP产生的靶后破片质量的影响规律。结果表明：相比于改进前的模型，改进后的模型能够更准确地解释靶板和EFP产生的靶后破片质量随靶板厚度和EFP着靶速度的变化规律；当EFP着靶速度为1 650 m/s时，随着靶板厚度从30 mm增大到70 mm，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断增强；当靶板厚度为40 mm时，随着EFP着靶速度从1 650 m/s升高到1 860 m/s，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断减弱。
- 爆炸成型弹丸 /
- 变截面特性 /
- 有限厚靶板 /
- 质量模型 /
- 靶后破片
Abstract: The function of an explosively-formed projectile (EFP) is not only penetrating the armor but also destroying the equipment behind the armor by behind-armor debris (BAD). It is necessary to predict the mass of BAD since the mass distribution of BAD is an important measurement to evaluate the level of destruction caused by BAD. A mass model for BAD generated by normal penetration of an EFP into an armor steel plate was improved by considering the variable cross-section characteristic of the EFP, basing on Bernoulli's equation of the flowing fluid and adiabatic shear theory. The accuracy of the model was validated by the experimental data and numerical simulation results. Therefore, the influences of the thickness of the target and the impact velocity of the EFP on the mass of BAD generated by the target and EFP were investigated. The results indicate: (1) compared with the previous model, the improved model can more accurately explain the mass variation of BAD generated by the target and EFP with the thickness of the target and the impact velocity of the EFP; (2) as the impact velocity of the EFP is 1 650 m/s, with the thickness of the target increasing from 30 mm to 70 mm, the effect of the variable cross-section characteristics on the mass of BAD generated by the target and EFP is constantly increasing; (3) as the thickness of the target is 40 mm, with the increase of the initial velocity of the EFP from 1 650 m/s to 1 860 m/s, the effect of the variable cross-section characteristics on the mass of BAD generated by the target and EFP is constantly decreasing.
- explosively-formed projectile /
- variable cross-section characteristic /
- finite thickness target /
- mass model /
- behind-armor debris

HTML全文

爆炸成型弹丸（explosively-formed projectile，EFP）在侵彻装甲后形成的破片会对装甲内部人员和仪器造成毁伤^[1-2]，因此通过靶板厚度和EFP着靶速度预测靶后破片质量、数量以及速度分布特性对EFP靶后毁伤效应评估、工程防护与设计具有重要意义。而准确计算EFP垂直侵彻有限厚装甲钢板时靶板和EFP产生的靶后破片质量，是预测上述分布特性的重要基础。

在侵彻问题中，Alekseevskii^[3]和Tate^[4-7]所建立的A-T模型取得了很好的效果，很多学者对此模型进行了改进和发展^[8-13]，但所涉及侵彻体的主体部分都是细长等截面圆柱杆，没有考虑小长径比和变截面的情况。文献[3-9]在研究长杆弹侵彻问题时忽略了蘑菇头轴向长度对侵彻过程的影响；文献[14-15]在研究EFP侵彻问题时忽略了蘑菇头轴向长度对侵彻过程的影响；文献[15]在研究EFP垂直侵彻RHA（rolled homogenous armour）靶板时，忽略了EFP的强度，直接套用了射流的扩孔模型：这些都是与实际情况不符的。目前对于靶后效应的理论研究均是基于侵彻体的主体部分是细长等截面圆柱杆的情况^[2,16]，对于EFP靶后效应的研究多集中在仿真和试验方面^[17-19]，而在EFP靶后效应的理论研究中均未考虑EFP变截面的特性^[1,14-15]。EFP主要用于末敏弹的战斗部，这类EFP具有变截面的特性，因此不能直接套用主体部分是细长等截面圆柱杆的侵彻理论；由于其长径比小，因此不能忽略蘑菇头轴向长度对侵彻过程的影响；由于其头部形状不规则，因此在侵彻有限厚钢靶板后产生的冲塞体形状不确定；由于其炸高很大（可达1 000倍装药直径），因此不能像射流一样忽略其强度。基于上述原因，使得利用现有的理论计算EFP垂直侵彻有限厚装甲钢板时靶板和EFP产生的靶后破片质量将会带来较大的偏差。

因此，本文将在考虑了EFP蘑菇头轴向长度、EFP强度的基础上，对冲塞体形状进行假设，并运用理论的方法分别分析考虑和不考虑EFP变截面的特性时，靶板和EFP产生的靶后破片质量。本文将借助文献[20]中的的试验数据验证理论分析中的关键参量，借助文献[18-19]中的的试验数据验证仿真方法的可信性，并利用此仿真方法对某典型EFP垂直侵彻装甲钢板时靶板和EFP产生的靶后破片质量进行计算，利用所得的仿真结果分析和判断EFP变截面的特性对靶板和EFP产生的靶后破片质量的影响。

1. 质量模型

1.1 基本假设

冲塞体是当EFP侵彻有限厚靶板到一定深度时，在靶板背面形成的一个内部布满裂纹的金属块^[16]，部分消蚀的EFP内部也布满裂纹。冲塞体和部分消蚀的EFP内部的裂纹会随着侵彻的进行，逐渐延伸扩展，当侵彻完成，冲塞体和部分消蚀的EFP将在靶板后形成若干碎块，即为靶板和EFP产生的靶后破片。为了简化侵彻过程，本文中做如下假设：（1）EFP开坑阶段所需的时间短、消耗的能量少，因此不计EFP开坑对侵彻过程的影响。（2）初始扩孔的轴向压力与径向压力相等，径向压力与扩孔面积成反比。（3）针对本文所涉及的EFP垂直侵彻装甲钢板的情况，假设冲塞体为圆台且母线与底面夹角45°。理由是：Backman等^[21]指出脆性靶板产生靶后破片时靶后出孔较大呈喇叭状；Yarin等^[16]在研究钨合金长杆弹侵彻装甲板时在其理论中指出靶后出孔较大呈喇叭状；Arnold等^[22]和Dalzell等^[17]分别运用仿真的方法发现射流和EFP垂直侵彻装甲板时靶后出孔较大呈喇叭状；王昕等^[18]、李睿等^[15]在运用试验的方法研究EFP垂直侵彻装甲板时发现侵彻出孔呈锥形，外侧直径较大；张先锋等^[14]在研究EFP对有限厚靶板侵彻过程及后效中认为冲塞崩落块为圆台且母线与底面夹角45°。

1.2 变截面EFP对装甲钢板的侵彻模型

邢柏阳等^[23]、Xing等^[24]在之前的研究中发现，密实结构EFP具有旋转对称的特性，可以用任意通过其旋转轴的截面表示其形状。本文中根据某典型EFP的激光高速照相测量结果（如图1上半部分所示），绘制了EFP垂直侵彻装甲钢板的过程示意图，如图2所示，H₀为靶板厚度，H为EFP未消蚀部分和消蚀部分交界面距离靶板底面的距离，R表示侵彻孔（蘑菇头）半径，u为EFP的侵彻速度，v表示EFP未消蚀部分的速度，r为EFP未消蚀部分和消蚀部分交界面处的EFP截面半径，l为EFP未消蚀部分和消蚀部分交界面距离EFP尾部的距离，P为EFP的侵彻深度。

图 1 某典型爆炸成型弹丸的形状

Figure 1. The shape of a certain typical explosively-formed projectile

下载: 全尺寸图片幻灯片

图 2 爆炸成型弹丸（EFP）垂直侵彻靶板的过程

Figure 2. The process of an explosively-formed projectile (EFP) normally penetrating into a target

下载: 全尺寸图片幻灯片

由于EFP具有变截面的特性且没有函数表达式可以描述其横截面半径与其轴向位置的对应关系，因此只能根据EFP的激光高速照相测量结果获取横截面半径（r_i）与其轴向位置（l_i）关系的离散点数据，记为r（l_i），l为EFP未消蚀部分的长度，由于离散点的数量足够大，因此 ${l_{i + 1}} - {l_i} \ll l$ ，可以认为dl=l_i+1−l_i，结合牛顿第二定律，可知EFP未消蚀部分的减速度满足以下关系，

$\quad\quad\quad\quad\quad\frac{{{\rm{d}}v}}{{{\rm{d}}t}} = - \frac{{{Y_{\rm{p}}}r_i^2}}{{\int_0^l {{\rho _{\rm{p}}}r_i^2{\rm{d}}l} }}\quad\quad\quad\quad\quad\quad\quad(1)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad$

式中：i为自某典型EFP尾部起的离散点序号，由于EFP具有变截面的特性，使得没有函数表达式可以描述其横截面半径与其轴向位置的对应关系，因此通过对激光高速照相图像识别而获得离散点；l_i为第i个离散点处与某典型EFP尾部的距离；r_i为某典型EFP的l_i处外轮廓横截面半径；dl为某典型EFP上任意两个离散间的距离，i足够大， ${\rm{d}}l = {l_{i + 1}} - {l_i} \ll l$ ；t为侵彻时间，从t=0时刻开始侵彻；ρ_p为EFP的密度；Y_p为EFP的强度因子。

不考虑EFP变截面的特性时，即认为EFP轴上所有位置处的横截面半径均相同，其值为前述所有离散点横截面半径（r_i）的平均值：

$r' = \frac{1}{k}\sum\limits_{i = 1}^k {{r_i}}$

(2)

式中：k为离散点的数量。此时EFP未消蚀部分的减速度关系即为A-T模型^[3-7]中的减速度关系：

$\frac{{{\rm{d}}v}}{{{\rm{d}}t}} = - \frac{{{Y_{\rm{p}}}}}{{{\rho _{\rm{p}}}l}}$

(3)

考虑和不考虑EFP变截面的特性时，EFP的消蚀速度、侵彻速度均满足A-T模型^[3-7]：

$\frac{{{\rm{d}}l}}{{{\rm{d}}t}} = - (v - u)$

(4)

$\frac{{\rm{1}}}{{\rm{2}}}{\rho _{\rm{p}}}{(v - u)^2} + {Y_{\rm{p}}} = \frac{{\rm{1}}}{{\rm{2}}}{\rho _{\rm{t}}}{u^2} + {\sigma _{\rm{t}}}$

(5)

式中： ${\rho _{\rm{t}}}$ 、 ${\sigma _{\rm{t}}}$ 分别为靶板的密度和强度因子。

由于本文所研究的某典型EFP的材料为紫铜，因此 ${Y_{\rm{p}}} \text{＜} {\sigma _{\rm{t}}}$ ，由式（5）可以得到：

$u = \frac{{v - \sqrt {\displaystyle\frac{{{\rho _{\rm{t}}}}}{{{\rho _{\rm{p}}}}}{v^2} + \left( {1 - \displaystyle\frac{{{\rho _{\rm{t}}}}}{{{\rho _{\rm{p}}}}}} \right)v_{{\rm{cr}}}^2} }}{{1 - \displaystyle\frac{{{\rho _{\rm{t}}}}}{{{\rho _{\rm{p}}}}}}}$

(6)

${v_{{\rm{cr}}}} = \sqrt {\frac{{2({\sigma _{\rm{t}}} - {Y_{\rm{p}}})}}{{{\rho _{\rm{p}}}}}}$

(7)

式中：v_cr为临界侵彻速度。当v＞v_cr时，EFP可以侵彻靶板；当v≤v_cr时，EFP停止侵彻靶板。

分析图2中的受力关系，得到EFP对靶板的轴向作用力满足：

${F_{\rm{t}}} = {\rm{\pi }}{R^2}\left( {\frac{1}{2}{\rho _{\rm{t}}}{u^2} + {\sigma _{\rm{t}}}} \right)$

(8)

由假设（2）可以得到初始扩孔的径向压力p₀、侵彻过程中扩孔的径向压力p与扩孔半径的关系：

${p_{\rm{0}}} = \frac{1}{2}{\rho _{\rm{p}}}{({v_{\rm{0}}} - {u_{\rm{0}}})^2} + {Y_{\rm{p}}}$

(9)

$p = \frac{{{r^2}}}{{{R^2}}}{p_{\rm{0}}}$

(10)

式中：u₀为t=0时刻的侵彻速度，v₀为EFP着靶速度，r为EFP消蚀部分和未消蚀部分交界面处的横截面半径。可以发现，r是影响EFP-靶板交界面处压力的关键因素。另外，r的取值视情况而定，当考虑EFP变截面的特性时取r_i，当不考虑EFP变截面的特性时取 $r'$ ，后文中r的取值规则与此相同。

侵彻过程中扩孔的径向压力可以使用含有靶板强度因子的流体力学Bernoulli方程得到:

$p = \frac{1}{2}{\rho _{\rm{t}}}{\left( {\frac{{{\rm{d}}R}}{{{\rm{d}}t}}} \right)^2} + {\sigma _{\rm{t}}}$

(11)

由式（9）～（11）可以得到扩孔速度：

$\frac{{{\rm{d}}R}}{{{\rm{d}}t}} = \sqrt {\frac{{{r^2}}}{{{R^2}}}\frac{{{\rho _{\rm{p}}}{{({v_{\rm{0}}} - {u_{\rm{0}}})}^2} + 2{Y_{\rm{p}}}}}{{{\rho _{\rm{t}}}}} - \frac{{2{\sigma _{\rm{t}}}}}{{{\rho _{\rm{t}}}}}}$

(12)

由式（12）可以发现，扩孔速度与R负相关，与r和v₀正相关。当考虑EFP变截面的特性时，随着侵彻进行：结合图1可以发现r先快速增大再缓慢变化，然后又快速增大，再缓慢变化；结合式（12）可知R先快速增大（bI阶段），再缓慢增大（bII阶段，R≈R_bII（常数）），然后又快速增大（bIII阶段），再缓慢增大（bIV阶段）。当不考虑EFP变截面的特性时，随着侵彻的进行，r恒为定值，由于扩孔速度恒为正值，但r/R的值逐渐减小，因此R先快速增大（dI阶段），然后一直保持缓慢增大（dII阶段，R≈R_dII（常数）），理论分析并计算可得R_bII≈R_dII。

1.3 靶板和EFP产生的靶后破片质量

冲塞体形成时刻 t_c的各参数关系如图3所示，网格表示冲塞体内部的碎块，H_c为冲塞体高度，R_c为冲塞体形成时刻的侵彻孔（蘑菇头）半径，u_c为冲塞体形成时刻的EFP的侵彻速度，v_c为冲塞体形成时刻的EFP未消蚀部分的速度，r_c为冲塞体形成时刻的EFP未消蚀部分和消蚀部分交界面处的EFP截面半径，l_c为冲塞体形成时刻的EFP未消蚀部分和消蚀部分交界面距离EFP尾部的距离，P_c为冲塞体形成时刻的EFP的侵彻深度。

图 3 冲塞体形成时刻

Figure 3. Forming time of a plug

下载: 全尺寸图片幻灯片

根据假设（3），冲塞体为圆台且母线与底面夹角45°，因此冲塞体形成时受到的轴向剪切力如式（13）左边所示，此时其值与EFP对靶板轴向的作用力相等，如式（13）右边所示，即：

$\sqrt 2 {\text{π}}{H_{\rm{c}}}({H_{\rm{c}}} + 2{R_{\rm{c}}}){\tau _{\max }} = {\text{π}}R_{\rm{c}} ^2\left( {\frac{1}{2}{\rho _{\rm{t}}}u_{\rm{c}}^2 + {\sigma _{\rm{t}}}} \right)$

(13)

式中：τ_max为靶板的最大剪切强度。τ_max的计算方法如下：

${\tau _{\max }} = {\tau _0}\gamma _i^n\left[ {1 + {{\left( {\frac{{{{\dot \gamma }_{\rm{m}}}}}{{q'}}} \right)}^{1/q}}} \right]\exp \left\{ { - \frac{{\alpha \beta }}{{{\rho _{\rm{t}}}{C_V}}}\left[ {1 + {{\left( {\frac{{{{\dot \gamma }_{\rm{m}}}}}{{q'}}} \right)}^{1/q}}} \right]\frac{{{\tau _0}}}{{n + 1}}\gamma _i^{n + 1}} \right\}$

(14)

${\dot \gamma _{\rm{m}}} = \frac{{1 - n}}{n}\frac{u}{R}$

(15)

${\gamma _i} = {\left\{ {\frac{{n{\rho _{\rm{t}}}{C_V}}}{{{\tau _0}\alpha \beta \left[ {1 + {{\left( {\displaystyle\frac{{{{\dot \gamma }_{\rm{m}}}}}{{p'}}} \right)}^{1/q}}} \right]}}} \right\}^{1/(n + 1)}}$

(16)

式中：τ₀为γ=1时材料的剪切强度；n为材料的应变硬化指数；q和 $q'$ 为材料的应变率敏感指数；α为应力的温度软化系数；C_V为材料的定容比热容；β为塑形功转化为热量的比例系数（Taylor-Quinney系数^[25]）； ${\gamma _i}$ 为剪应变； ${\dot \gamma _{\rm{m}}}$ 为平均剪应变率。

式（12）～（13）表明：r可以影响冲塞体形成时刻的侵彻孔（蘑菇头）半径（R_c），并进一步影响冲塞体形成时刻的冲塞体高度（H_c）。

考虑EFP蘑菇头轴向方向的长度对侵彻过程的影响，有如下关系：

${P_{\rm{c}}} + {H_{\rm{c}}} - {H_0} - {R_{\rm{c}}} = 0$

(17)

${P_{\rm{c}}} = \int_0^{{t_{\rm{c}}}} {u{\rm{d}}t}$

(18)

式（17）～（18）表明r将会影响冲塞体形成时刻的侵彻深度（P_c），并最终决定冲塞体的形成时刻（t_c）。因此冲塞体质量，即靶板产生的靶后破片质量，满足以下关系：

${m_{\rm{t}}} = {\rho _{\rm{t}}}\left\{ {\frac{1}{3}{\text{π}}{H_{\rm{c}}}[R_{\rm{c}}^2 + {{({R_{\rm{c}}} + {H_{\rm{c}}})}^2} + {R_{\rm{c}}}({R_{\rm{c}}} + {H_{\rm{c}}})] - \frac{2}{3}{\text{π}}R_{\rm{c}}^3} \right\}$

(19)

由式（13）和（19）可以得到：

${m_{\rm{t}}} = \frac{{\text{π}}}{3}{\rho _{\rm{t}}}R_{\rm{c}}^3f({u_{\rm{c}}})$

(20)

在本文讨论的靶板厚度（30 mm≤H₀≤70 mm）和EFP着靶速度（1 650 m/s≤v₀≤1 860 m/s）范围内，f（u_c）＞0恒成立，因此m_t与R_c总是呈正相关，并且当R_c变化幅度较大时，R_c是影响m_t的主要因素；当R_c变化幅度较小时，u_c是影响m_t的主要因素。

冲塞体形成时EFP总的剩余质量即为EFP产生的靶后破片质量，即：

${m_{\rm{p}}} = {\rho _{\rm{p}}}{\text{π}}\int_0^{{l_{\rm{c}}} + {R_{\rm{c}}}} {{r^2}} {\rm{d}}l$

(21)

基于上述理论分析可以发现，r通过影响扩孔速度，进而影响冲塞体形成时刻的各个参量（R_c、H_c、P_c、t_c等），最终影响靶板和EFP产生的靶后破片质量。本文将基于某典型EFP，对不同靶板厚度和EFP着靶速度条件下，考虑和不考虑EFP变截面的特性时，所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，进行分析和比较，具体地验证改进后的模型可以减小与实际情况的偏差。

2. 模型验证

2.1 侵彻孔半径模型验证

Yarin等^[16]在通过试验研究钨合金长杆弹侵彻装甲板时发现，在统计靶板的质量损失时结果波动很大，接近30%，因此本文中将利用试验中可以稳定测量的参数验证力学模型的准确性。由于Held等^[20]在侵彻孔半径模型方面所做的工作非常经典，并且射流侵彻水介质和EFP侵彻装甲钢都属于侵彻体强度小于靶板强度的模式，因此将文献[20]中的初始条件代入式（12）中，得到的侵彻孔半径（R）随侵彻时间（t）的变化曲线如图4中的红线所示，黑点为Held等^[20]的理论结果。

图 4 侵彻孔半径随时间变化关系

Figure 4. The time history of the crater radius

下载: 全尺寸图片幻灯片

由图4可以发现，利用本文理论模型所得的侵彻孔半径（R）随侵彻时间（t）的变化关系与Held等^[20]的理论结果十分吻合，表明本文理论模型具有较高的准确性。

2.2 仿真模型验证

Dalzell等^[17]采用AUTODYN-3D^TM对EFP侵彻靶板形成破片云的过程进行了数值模拟，指出SPH（smoothed particle hydrodynamics）算法较Lagrange算法和Euler算法在研究靶后破片方面更具优势。在试验中所统计的靶板质量损失波动很大^[16]，并且，经常发现靶板和EFP产生的靶后破片会粘结在一起。因此，本文首先利用试验中可以稳定测量的结果^[18-19]验证本文仿真方法的可信性，进而利用得到验证的仿真方法对某典型EFP进行仿真计算，验证理论模型所得的、但在试验统计中波动较大的结果。

本文使用AUTODYN-3D^TM的SPH算法，仿真模型中EFP和靶板的材料种类、状态方程以及本构方程取自文献[19,26]。文献[18-19]的试验中，EFP战斗部由药型罩、主装药、传爆药、壳体组成，其装药直径（D）为56 mm，药型罩为壁厚3 mm的球缺形紫铜罩，装药为长径比0.86的压装JH-2，壳体厚度为0.045D。EFP垂直侵彻厚度为20 mm的45钢板，EFP的形状、尺寸如图5下半部分的X光照片所示，试验中测得EFP着靶速度为2 120 m/s，45钢板的长、宽均为200 mm。在仿真中构建与文献[18-19]的试验中完全一致的初始条件，仿真中EFP的形状、尺寸如图5上半部分所示，并对靶板边界施加固定约束，EFP侵彻45钢板前的仿真模型如图6所示。

图 5 依据文献[18-19]构建的EFP模型

Figure 5. The EFP model developed according to references [18-19]

下载: 全尺寸图片幻灯片

图 6 依据文献[18-19]构建的仿真模型

Figure 6. The simulation model developed by according to references [18-19]

下载: 全尺寸图片幻灯片

仿真所得侵彻80 µs和95 µs时的靶后破片云形状如图7（a）、（b）上半部分所示，文献[18]中给出的相应时刻的X光照片如图7（a）、（b）下半部分所示，对比发现仿真所得的靶后破片云图和文献[18]中试验所得的X光照片非常吻合。另外，仿真所得侵彻95 µs时的靶后破片云径向膨胀速度v_J=531 m/s,靶后破片云头部膨胀速度v_T=1 690 m/s与文献[18]中的试验结果v_J=525 m/s，v_T=1 678 m/s的误差不超过2%，表明本文所用的仿真方法是可信的。

图 7 本文仿真所得的靶后破片云与文献[18]的对比

Figure 7. Behind-armor debris clouds obtained by this paper compared with that in reference [18]

下载: 全尺寸图片幻灯片

3. 对比分析

图1中某典型EFP战斗部装药直径（D）为125 mm，装药长度为0.84D，药型罩为壁厚0.064D的球缺形紫铜罩，主装药、传爆药均为8701，周侧壳体厚度为0.048D，底部壳体厚度为0.12D，由西安现代控制技术研究所设计并由辽沈工业集团有限公司生产。某典型EFP的最大半径为28.5 mm，最大长度为98 mm，其形状及尺寸均在图1中标明，其激光高速照相测量结果如图1上半部分所示，其仿真中的模型如图1下半部分所示，并对靶板边界施加固定约束。使用上述仿真方法计算某典型EFP以着靶速度1 650 m/s垂直侵彻厚度为70 mm的装甲钢板，仿真得到图8所示的靶板出孔形状，其中黑色虚线表示靶板出孔侧面，黄色虚线表示冲塞体形成前EFP的侵彻路径，可以发现，靶板出孔与靶板底面夹角接近45°，与文献[14-18,21-22]比较吻合，因此假设（2）是具有可信性的。

图 8 某典型EFP垂直侵彻靶板的出孔形状(单位：mm)

Figure 8. The exit hole shape of a target normally penetrated by a certain typical EFP (unit in mm)

下载: 全尺寸图片幻灯片

质量模型中的参数全部取自文献[7,13,19,25-26]，Y_p和R_t的计算方法由文献[7]给出，主要参数为：τ₀=0.907 GPa，n=0.25， $q'$ =8.0×10⁷ s⁻¹，q=7.33，α=6.8×10⁻⁴，β=0.9，C_V=452 J/（kg·K），靶板的杨氏模量E_t=210 GPa，Tate给出的计算Y_p时的系数λ=0.7，ρ_p=8.93 g/cm³，ρ_t=7.85 g/cm³。

3.1 靶板厚度对靶后破片质量的影响

由式（19）和（21）可以发现，冲塞体的形状参数直接决定了靶板产生的靶后破片质量，间接影响了EFP产生的靶后破片质量。由式（18）可以发现，冲塞体的形状参数与靶板厚度密切相关，因此靶板厚度是影响靶后破片质量的主要因素之一。在某典型EFP着靶速度为1 650 m/s的情况下，垂直侵彻厚度为30、40、50、60、70 mm的装甲钢板，靶板和EFP产生的靶后破片质量如图9所示。

图 9 不同靶板厚度条件下的靶后破片质量

Figure 9. Mass of behind-armor debris for different thicknesses of targets

下载: 全尺寸图片幻灯片

通过编写Matlab^TM程序进行图像识别，读取某典型EFP的激光高速照相测量结果（如图7上半部分所示），辨识其边界点，从而获取横截面半径与其轴向位置关系的离散点数据，对式（20）和（21）进行数值求解，可以得到考虑EFP变截面的特性时靶板和EFP产生的靶后破片质量，分别记为m_t，bjm、m_p，bjm。利用所得的离散点数据以及式（2），求出r'=16.23 mm，对式（20）和（21）进行数值求解，可以得到不考虑EFP变截面的特性时靶板和EFP产生的靶后破片质量，分别记为m_t，djm、m_p，djm。使用本文仿真方法所得的靶板和EFP产生的靶后破片质量，分别记为m_t，sim、m_p，sim。

由图9可以发现，相比于不考虑EFP变截面特性的结果，考虑EFP变截面特性的结果与仿真结果更吻合，能够更加准确地解释靶板和EFP产生的靶后破片质量随靶板厚度的变化规律。

考虑EFP变截面的特性时，随着靶板厚度的增大，侵彻时间不断增长，在本节研究的靶板厚度范围内R_c逐渐由bII阶段过渡到bIII阶段，最后进入bIV阶段。由于bIII阶段R_c快速增大，由理论分析可知，此时R_c是影响m_t的主要因素，结合式（20）可知靶板产生的靶后破片质量先基本不变，然后快速增加，最后缓慢增加。随着靶板厚度的增大，EFP的消蚀长度不断增大，因此EFP产生的靶后破片质量不断减小。

不考虑EFP变截面的特性时，由理论分析并计算可得R_bII≈R_dII，并且在本节研究的靶板厚度范围内R_c始终处于dII阶段，因此对于厚度较小的靶板，考虑与不考虑EFP变截面所得的靶板产生的靶后破片质量接近；但随着靶板厚度的增大，考虑EFP变截面的R_c开始快速超过不考虑EFP变截面的R_c，因此考虑与不考虑EFP变截面所得的靶板产生的靶后破片质量差不断增大。由于此时R_c的变化极其缓慢，由理论分析可知，u_c是影响m_t的主要因素，随着靶板厚度的增大，u_c不断减小，但由于u_c对m_t的影响比较小，因此不考虑EFP变截面的特性时，靶板产生的靶后破片质量缓慢减小。由于不考虑EFP变截面时EFP的头部质量大于考虑EFP变截面时EFP的头部质量，这使得在侵彻开始阶段不考虑EFP变截面时EFP损失的质量较大，因此不考虑EFP变截面所得的EFP产生的靶后破片质量总是小于考虑EFP变截面所得的结果。

计算以上5种工况条件下，考虑和不考虑EFP变截面的特性时，所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，如表1所示。举例，考虑EFP变截面的特性时所得的靶板产生的靶后破片质量与仿真结果的偏差的具体计算方法为：ε_t，bjm=|（m_t，bjm−m_t，sim）/m_t，sim|×100%，另外3个偏差的计算方法与此类似。

表 1 不同靶板厚度条件下靶后破片质量的偏差

Table 1. Mass deviations of behind-armor debris for different thicknesses of targets

H₀/mm	ε_t，bjm/%	ε_t，djm/%	ε_p，bjm/%	ε_p，djm/%
30	4.1	13.5	3.5	32.0
40	1.8	12.3	2.9	40.3
50	0.8	25.9	4.1	45.4
60	7.4	55.4	7.4	54.8
70	6.0	70.1	10.4	64.3

下载: 导出CSV

| 显示表格

由表1可以发现，考虑EFP变截面的特性时所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，均明显小于相同靶板厚度条件下不考虑EFP变截面的特性时所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，并且，随着靶板厚度的增大，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断增强。表明，相比于改进前的模型，改进后的模型能够更准确地解释靶板和EFP产生的靶后破片质量随靶板厚度的变化规律。

3.2 EFP着靶速度对靶后破片质量的影响

由式（12）可以发现，EFP着靶速度对扩孔速度有显著的影响，而扩孔速度与冲塞体形成的临界条件密切相关，因此EFP着靶速度是影响靶后破片质量的主要因素之一。在装甲钢板厚度为40 mm的情况下，某典型EFP分别以1 650、1 680、1 740、1 800、1 860 m/s的着靶速度进行垂直侵彻，靶板和EFP产生的靶后破片质量如图10所示。

图 10 不同EFP着靶速度条件下的靶后破片质量

Figure 10. Mass of behind-armor debris for different impact velocities of EFPs

下载: 全尺寸图片幻灯片

由图10可以发现，相比于不考虑EFP变截面特性的结果，考虑EFP变截面特性的结果与仿真结果更吻合，能够更准确地解释靶板和EFP产生的靶后破片质量随EFP着靶速度的变化规律。

考虑EFP变截面的特性时，由于靶板厚度固定为40 mm，当EFP着靶速度为1 650 m/s时，侵彻时间最长，但是由3.1节中的分析可知，此时的R_c处于bII阶段，即R_c缓慢增大，此时影响R_c的主要因素是EFP的着靶速度，次要因素是侵彻时间。因此在本节研究的EFP着靶速度范围内，随着EFP着靶速度的增高，R_c不断增大，同时u_c不断增大，结合式（20）可知靶板产生的靶后破片质量不断增大。由于侵彻时间不断减短，因此EFP的消蚀长度不断减小，因此EFP的剩余长度不断增大，使得EFP产生的靶后破片质量不断增加。

同理，不考虑EFP变截面的特性时，此时的R_c处于dII阶段，影响R_c的主要因素是EFP的着靶速度。因此在本节研究的EFP着靶速度范围内，随着EFP着靶速度的增高，靶板产生的靶后破片质量不断增大。由于侵彻时间不断减短，使得EFP产生的靶后破片质量不断增加。

计算以上5种工况条件下，考虑和不考虑EFP变截面的特性时，所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，如表2所示。

表 2 不同EFP着靶速度条件下靶后破片质量的偏差

Table 2. Mass deviations of behind-armor debris for different impact velocities of EFPs

v₀/（m·s⁻¹）	ε_t，bjm/%	ε_t，djm/%	ε_p，bjm/%	ε_p，djm/%
1 650	1.8	12.3	2.9	40.3
1 680	3.9	11.2	2.4	38.6
1 740	4.9	9.6	1.2	37.1
1 800	2.7	5.0	1.2	36.1
1 860	2.2	2.3	0.4	34.1

下载: 导出CSV

| 显示表格

由表2可以发现：考虑EFP变截面的特性时，所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差，均明显小于相同EFP着靶速度条件下不考虑EFP变截面的特性时所得的靶板和EFP产生的靶后破片质量与仿真结果的偏差；并且，随着EFP着靶速度的提高，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断减弱，当EFP着靶速度达到1 860 m/s时，考虑和不考虑EFP变截面的特性所得的靶板产生的靶后破片质量十分接近。表明，相比于改进前的模型，改进后的模型能够更准确地解释靶板和EFP产生的靶后破片质量随EFP着靶速度的变化规律。

4. 结　论

本文中建立了变截面爆炸成型弹丸（explosively-formed projectile，EFP）垂直侵彻装甲钢板靶后破片质量模型，并得到以下结论：（1）相比于改进前的模型，改进后的模型能够更准确地解释靶板和EFP产生的靶后破片质量随靶板厚度和EFP着靶速度的变化规律，可以减小与实际情况的偏差；（2）当EFP着靶速度为1 650 m/s时，随着靶板厚度从30 mm增大到70 mm，靶板产生的靶后破片质量先保持基本不变，然后快速增大，再缓慢增大，EFP产生的靶后破片质量不断减小，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断增强；（3）当靶板厚度为40 mm时，随着EFP着靶速度从1 650 m/s提高到1 860 m/s，靶板和EFP产生的靶后破片质量均不断增大，EFP变截面的特性对靶板和EFP产生靶后破片质量的影响不断减弱。

图 1 某典型爆炸成型弹丸的形状

Figure 1. The shape of a certain typical explosively-formed projectile

下载: 全尺寸图片幻灯片

图 2 爆炸成型弹丸（EFP）垂直侵彻靶板的过程

Figure 2. The process of an explosively-formed projectile (EFP) normally penetrating into a target

下载: 全尺寸图片幻灯片

图 3 冲塞体形成时刻

Figure 3. Forming time of a plug

下载: 全尺寸图片幻灯片

图 4 侵彻孔半径随时间变化关系

Figure 4. The time history of the crater radius

下载: 全尺寸图片幻灯片

图 5 依据文献[18-19]构建的EFP模型

Figure 5. The EFP model developed according to references [18-19]

下载: 全尺寸图片幻灯片

图 6 依据文献[18-19]构建的仿真模型

Figure 6. The simulation model developed by according to references [18-19]

下载: 全尺寸图片幻灯片

图 7 本文仿真所得的靶后破片云与文献[18]的对比

Figure 7. Behind-armor debris clouds obtained by this paper compared with that in reference [18]

下载: 全尺寸图片幻灯片

图 8 某典型EFP垂直侵彻靶板的出孔形状(单位：mm)

Figure 8. The exit hole shape of a target normally penetrated by a certain typical EFP (unit in mm)

下载: 全尺寸图片幻灯片

图 9 不同靶板厚度条件下的靶后破片质量

Figure 9. Mass of behind-armor debris for different thicknesses of targets

下载: 全尺寸图片幻灯片

图 10 不同EFP着靶速度条件下的靶后破片质量

Figure 10. Mass of behind-armor debris for different impact velocities of EFPs

下载: 全尺寸图片幻灯片

表 1 不同靶板厚度条件下靶后破片质量的偏差

Table 1. Mass deviations of behind-armor debris for different thicknesses of targets

H₀/mm ε_t，bjm/% ε_t，djm/% ε_p，bjm/% ε_p，djm/%

30 4.1 13.5 3.5 32.0
40 1.8 12.3 2.9 40.3
50 0.8 25.9 4.1 45.4
60 7.4 55.4 7.4 54.8
70 6.0 70.1 10.4 64.3

下载: 导出CSV

表 2 不同EFP着靶速度条件下靶后破片质量的偏差

Table 2. Mass deviations of behind-armor debris for different impact velocities of EFPs

v₀/（m·s⁻¹） ε_t，bjm/% ε_t，djm/% ε_p，bjm/% ε_p，djm/%

1 650 1.8 12.3 2.9 40.3
1 680 3.9 11.2 2.4 38.6
1 740 4.9 9.6 1.2 37.1
1 800 2.7 5.0 1.2 36.1
1 860 2.2 2.3 0.4 34.1

下载: 导出CSV

参考文献(26)

[1]	KIM H S, ARNOLD W, HARTMANN T, et al. A model for behind armor debris from EFP impact [C] // BAKER E, TEMPLETON D. 26th International Symposium on Ballistics. USA: DEStech Publications Inc., 2011.
[2]	YOSSIFON G, YARIN A L. Behind-the-armor debris analysis [J]. International Journal of Impact Engineering, 2002, 27(8): 807–835. DOI: 10.1016/S0734-743X(02)00009-X.
[3]	ALEKSEEVSKⅡ V P. Penetration of a rod into a target at high velocity [J]. Combustion, Explosion and Shock waves, 1966, 2(2): 63–66. DOI: 10.1007/BF00749237.
[4]	TATE A. A theory for the deceleration of long rods after impact [J]. Journal of the Mechanics and Physics of Solids, 1967, 15(6): 387–399. DOI: 10.1016/0022-5096(67)90010-5.
[5]	TATE A. Further results in the theory of long rod penetration [J]. Journal of the Mechanics and Physics of Solids, 1969, 17(3): 141–150. DOI: 10.1016/0022-5096(69)90028-3.
[6]	TATE A. Long rod penetration models: Part I: a flow field model for high speed long rod penetration [J]. International Journal of Mechanical Sciences, 1986, 28(8): 535–548. DOI: 10.1016/0020-7403(86)90051-2.
[7]	TATE A. Long rod penetration models: Part Ⅱ: extensions to the hydrodynamic theory of penetration [J]. International Journal of Mechanical Sciences, 1986, 28(9): 599–612. DOI: 10.1016/0020-7403(86)90075-5.
[8]	孙庚辰, 吴锦云, 赵国志, 等. 长杆弹垂直侵彻半无限厚靶板的简化模型 [J]. 兵工学报, 1981(4): 1–8. SUN Gengchen, WU Jinyun, ZHAO Guozhi, et al. A simplified model of the penetration of the long-rod penetrator against the plates with semi-infinite thickness at normal angle [J]. Acta Armamentarii, 1981(4): 1–8.
[9]	ROSENBERG Z, MARMOR E, MAYSELESS M. On the hydrodynamic theory of long-rod penetration [J]. International Journal of Impact Engineering, 1990, 10(1/2/34): 483–486. DOI: 10.1016/0734-743X(90)90081-6.
[10]	WALKER J D, ANDERSON Jr C E. A time-dependent model for long-rod penetration [J]. International Journal of Impact Engineering, 1995, 16(1): 19–48. DOI: 10.1016/0734-743X(94)00032-R.
[11]	ZHANG Liansheng, HUANG Fenglei. Model for long-rod penetration into semi-infinite targets [J]. Journal of Beijing Institute of Technology (English Edition), 2004, 13(3): 285–289. DOI: 10.15918/j.jbit1004-0579.2004.03.011.
[12]	LAN Bin, WEN Heming. Alekseevskii-Tate revisited: an extension to the modified hydrodynamic theory of long rod penetration [J]. Science China Technological Sciences, 2010, 53(5): 1364–1373. DOI: 10.1007/s11431-010-0011-x.
[13]	何雨. 长杆弹撞击下金属靶板侵彻与穿透的进一步研究 [D]. 合肥: 中国科技大学, 2013: 54−57.
[14]	张先锋, 陈惠武, 赵有守. EFP对有限厚靶板侵彻过程及后效研究 [J]. 爆炸与冲击, 2006, 2(4): 323–327. DOI: 10.11883/1001-1455(2006)04-0323-05. ZHANG Xianfeng, CHEN Huiwu, ZHAO Youshou. Investigation of process and after effect of EFP penetration into target of finite thickness [J]. Explosion and Shock Waves, 2006, 2(4): 323–327. DOI: 10.11883/1001-1455(2006)04-0323-05.
[15]	李睿, 黄正祥, 祖旭东, 等. 靶板在爆炸成型弹丸垂直侵彻下的层裂 [J]. 爆炸与冲击, 2018, 38(5): 1039–1044. DOI: 10.11883/1001-1455(2018)05-1039-06. LI Rui, HUANG Zhengxiang, ZU Xudong, et al. Theoretical analysis on the spallation of EFP vertical penetration target [J]. Explosion and Shock Waves, 2018, 38(5): 1039–1044. DOI: 10.11883/1001-1455(2018)05-1039-06.
[16]	YARIN A L, ROISMAN I V, WEBER K, et al. Model for ballistic fragmentation and behind-armor debris [J]. International Journal of Impact Engineering, 2000, 24(2): 171–201. DOI: 10.1016/S0734-743X(99)00048-2.
[17]	DALZELL M W, HAZELL P J, MEULMAN J H. Modelling behind-armour debris formed by the perforation of an EFP through a steel target [C] // CARLEONE J, ORPHAL D. 20th International Symposium on Ballistics. USA: DEStech Publications Inc., 2002: 23−27.
[18]	王昕, 蒋建伟, 王树有, 等. 爆炸成型弹丸侵彻钢靶的后效破片云实验研究 [J]. 兵工学报, 2018, 39(7): 1284–1290. DOI: 10.3969/j.issn.1000-1093.2018.07.005. WANG Xin, JIANG Jianwei, WANG Shuyou, et al. Experimental research on fragments after explosively- formed projectile penetrating into steel target [J]. Acta Armamentarii, 2018, 39(7): 1284–1290. DOI: 10.3969/j.issn.1000-1093.2018.07.005.
[19]	WANG Yangyang, JIANG Jianwei, MENG Jiayu, et al. Effect of add-on explosive reactive armor on EFP penetration [C] // WOODLEY C, CULLIS I. 29th International Symposium on Ballistics. USA: DEStech Publications Inc., 2016: 2395−2406.
[20]	HELD M, HUANG N S, JIANG D, et al. Determination of the crater radius as a function of time of a shaped charge jet that penetrates water [J]. Propellants, Explosives, Pyrotechnics, 1996, 21(2): 64–69. DOI: 10.1002/prep.19960210203.
[21]	BACKMAN M E, GOLDSMITH W. The mechanics of penetration of projectiles into targets [J]. International Journal of Engineering Science, 1978, 16(1): 1–99. DOI: 10.1016/0020-7225(78)90002-2.
[22]	ARNOLD W, ROTTENKOLBER E. Physics of behind armor debris threat reduction [J]. International Journal of Impact Engineering, 2006, 33: 53–61. DOI: 10.1016/j.ijimpeng.2006.09.021.
[23]	邢柏阳, 郭锐, 刘荣忠, 等. 内嵌结构对末敏弹EFP成型影响研究 [J]. 弹箭与制导学报, 2016, 36(5): 37–40. DOI: 10.15892/j.cnki.djzdxb.2016.05.010. XING Boyang, GUO Rui, LIU Rongzhong, et al. Study on influence of the embedded structure on the terminal sensitive projectile EFP forming [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2016, 36(5): 37–40. DOI: 10.15892/j.cnki.djzdxb.2016.05.010.
[24]	XING Boyang, LIU Rongzhong, GUO Rui, et al. Influence of the embedded structure on the EFP formation of compact terminal sensitive projectile [J]. Defence Technology, 2017, 13(4): 310–315. DOI: 10.1016/j.dt.2017.05.006.
[25]	BAI Y L, JOHNSON W. Plugging: physical understanding and energy absorption [J]. Metals Technology, 1982, 9(1): 182–190. DOI: 10.1179/030716982803285945.
[26]	赵方宣, 沈兆欣, 刘宁, 等. 靶板材料对聚能射流跳弹角影响的数值模拟与试验 [J]. 含能材料, 2016, 24(1): 33–37. DOI: 10.11943/j.issn.1006-9941.2016.01.005. ZHAO Fangxuan, SHEN Zhaoxin, LIU Ning, et al. Numerical simulation and experimental research on the effect of target material on the ricochet angle of shaped charge jet [J]. Chinese Journal of Energetic Materials, 2016, 24(1): 33–37. DOI: 10.11943/j.issn.1006-9941.2016.01.005.

施引文献

期刊类型引用(2)

1.	黄炫宁，李伟兵，郭腾飞，李文彬，王晓鸣. EFP垂直侵彻靶板后效破片云形状描述. 含能材料. 2020(11): 1068-1075 . 百度学术
2.	邢柏阳，侯云辉，李泰华，张东江，刘荣忠，郭锐. 爆炸成型弹丸垂直侵彻装甲钢靶后破片动能分析. 兵工学报. 2019(10): 2014-2021 . 百度学术