冲击波与破片对波纹杂交夹层板的联合毁伤数值研究

李勇; 肖伟; 程远胜; 刘均; 张攀

doi:10.11883/bzycj-2016-0224

冲击波与破片对波纹杂交夹层板的联合毁伤数值研究

doi: 10.11883/bzycj-2016-0224

李勇¹,
肖伟²,
程远胜¹,
刘均¹,
张攀^1, ,

1.
华中科技大学船舶与海洋工程学院，湖北武汉 430074
2.
中国舰船研究设计中心，湖北武汉 430064

基金项目:

国家自然科学基金项目 51509096

国家自然科学基金项目 51209099

详细信息

作者简介:
李勇(1991-), 男, 硕士

通讯作者:
张攀, panzhang@hust.edu.cn

中图分类号: O383
计量
- 文章访问数: 5624
- HTML全文浏览量: 2009
- PDF下载量: 306
- 被引次数: 0
出版历程
- 收稿日期: 2016-08-11
- 修回日期: 2017-01-19
- 刊出日期: 2018-03-25

Numerical research on response of hybrid corrugated sandwich plates subjected to combined blast and fragment loadings

1.
School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
2.
China Ship Development and Design Center, Wuhan 430064, Hubei, China

摘要

摘要: 通过有限元软件LS-DYNA模拟了波纹杂交夹层板在冲击波与破片联合作用下的响应过程，研究了炸药当量、载荷类型和填充方式对波纹杂交夹层板变形与失效模式的影响，并与实体板、间隔板和波纹夹层板的抗联合毁伤性能进行了对比，讨论了波纹杂交夹层板的能量吸收特性。数值计算结果表明：与冲击波单独作用相比，破片群单独作用和冲击波与破片联合作用对结构造成的毁伤更为严重；当药量较小时，波纹夹层板和波纹杂交夹层板的抗联合毁伤性能优于实体板与间隔板，波纹杂交夹层板的抗联合毁伤性能从全填充、迎爆面填充到背爆面填充逐渐降低；当药量较大时，所有结构均产生破口失效；在能量耗散方面，冲击波单独作用时以波纹芯层吸能为主，破片群单独作用和冲击波与破片联合作用时以上面板吸能为主。
- 抗联合毁伤性能 /
- 冲击波与破片 /
- 波纹杂交夹层板 /
- 失效模式 /
- LS-DYNA
Abstract: The dynamic response of hybrid corrugated sandwich plates subjected to combined effects of blast and fragment loading was analyzed using finite element analysis code LS-DYNA. The effects of charge mass, loading type and filling strategy on deformation/failure pattern of hybrid corrugated sandwich plates were investigated. The comparison of anti combined loadings performance to three equivalent structures (solid plate, double-layered plate and corrugated sandwich plate) was made. Finally, the energy absorption characteristics of hybrid corrugated sandwich panels were discussed. Numerical results show that the damage extent of hybrid corrugated sandwich plates under bare fragment cluster loading or combined blast and fragment loading is more severe than that caused by bare blast loading. When the charge mass is small, the performances of corrugated sandwich plate and hybrid corrugated sandwich plate are superior to equivalent solid plate and double-layered plate. The corrugated sandwich panels with fully filling configuration possess the best damage resistance, followed by that with upper space filling configuration, and that with lower space filling configuration has the worst. All the structures fractured catastrophically when the charge mass is large. For energy absorption, the corrugated core is the main energy absorption part under bare blast loading, while the front face becomes the main energy absorption part under the other two loading conditions.
- anti combined loadings performance /
- blast and fragment /
- hybrid corrugated sandwich plate /
- failure pattern /
- LS-DYNA

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 波纹杂交夹层板有限元模型(1/4模型)

Figure 1. Finite element model of hybrid corrugated sandwich plate (1/4 model)

下载: 全尺寸图片幻灯片

图 2 实体板在联合作用下的毁伤

Figure 2. Failure of solid plate under combined blast and fragment loadings

下载: 全尺寸图片幻灯片

图 3 爆轰波传播过程

Figure 3. Propagation of detonation

下载: 全尺寸图片幻灯片

图 4 冲击波传播过程

Figure 4. Propagation of shock wave

下载: 全尺寸图片幻灯片

图 5 波纹杂交夹层板的响应过程

Figure 5. Response process of hybrid corrugated sandwich plate

下载: 全尺寸图片幻灯片

图 6 波纹杂交夹层板的失效模式

Figure 6. Failure patterns of hybrid plate

下载: 全尺寸图片幻灯片

图 7 联合作用时不同结构的失效模式

Figure 7. Failure patterns of different structures under combined blast and fragment loadings

下载: 全尺寸图片幻灯片

图 8 填充方式对波纹杂交夹层板失效的影响

Figure 8. Effect of filling strategy on the failure patterns of hybrid corrugated sandwich plates

下载: 全尺寸图片幻灯片

表 1 计算工况及数值结果

Table 1. Computational conditions and numerical results

工况编号	W/g	n	填充方式	中心点最大形变/mm		破损形式/塑性应变			变形能/kJ
工况编号	W/g	n	填充方式	f	b	f	b	c	f	b	c	foam
HP-1	166	121	F	-	-	破口	破口	破口	16.644	9.460	4.088	2.424
HP-2	166	121	U	-	-	破口	破口	破口	16.744	8.176	3.604	1.556
HP-3	166	121	L	-	-	破口	破口	破口	15.920	7.416	4.068	0.812
HP-4	0	121	F	-	-	破口	破口	破口	15.824	8.972	3.848	2.240
HP-5	166	0	F	-	19.15	破口	破口	损伤	6.404	7.240	2.056	1.916
HP-6	111	121	F	-	23.12	破口	破口	损伤	14.304	8.044	3.180	2.104
HP-7	111	121	U	-	23.47	破口	破口	损伤	14.600	7.328	2.516	1.336
HP-8	111	121	L	-	-	破口	破口	破口	13.756	6.412	3.180	0.624
HP-9	0	121	F	-	22.98	破口	破口	损伤	13.600	7.276	3.012	1.924
HP-10	111	0	F	-	15.40	破口	破口	损伤	5.140	6.308	1.488	1.776
HP-11	55	121	F	-	15.97	破口	破口	0.420	9.616	5.376	1.184	1.556
HP-12	55	121	U	-	16.61	破口	破口	损伤	9.472	4.652	0.896	1.024
HP-13	55	121	L	-	17.21	破口	破口	损伤	9.100	4.384	1.348	0.612
HP-14	0	121	F	-	15.60	破口	破口	0.395	9.544	5.284	1.084	1.492
HP-15	55	0	F	22.57	7.60	0.420	损伤	0.289	2.568	3.904	0.548	1.476
EP-1	166	121	-	-	-	破口	破口	破口	16.052	7.344	3.356	-
EP-2	111	121	-	-	23.16	破口	破口	损伤	13.160	6.200	2.280	-
EP-3	55	121	-	-	15.82	破口	破口	0.273	9.096	3.560	0.868	-
GP-1	166	121	-	-	-	破口	破口	-	20.244	-	5.852	-
GP-2	111	121	-	-	32.28	破口	破口	-	16.828	-	5.204	-
GP-3	55	121	-	-	22.07	破口	0.394	-	11.892	-	1.840	-
SP-1	166	121	-	-	-	破口	-	-	28.584	-	-	-
SP-2	111	121	-	-	-	破口	-	-	23.744	-	-	-
SP-3	55	121	-	23.41	-	损伤	-	-	13.876	-	-	-

下载: 导出CSV

参考文献(19)

[1]	侯海量, 张成亮, 李茂, 等.冲击波和高速破片联合作用下夹芯复合舱壁结构的毁伤特性[J].爆炸与冲击, 2015, 35(1):116-123. doi: 10.11883/1001-1455(2015)01-0116-08 HOU Hailiang, ZHANG Chengliang, LI Mao, et al. Damage characteristics of sandwich bulkhead under the impact of shock and high-velocity fragments[J]. Explosion and Shock Waves, 2015, 35(1):116-123. doi: 10.11883/1001-1455(2015)01-0116-08
[2]	何翔, 庞伟宾, 曲建波, 等.防护门在空气冲击波和破片作用下的破坏[J].爆炸与冲击, 2004, 24(5):475-479. http://www.bzycj.cn/CN/abstract/abstract9989.shtml HE Xiang, PANG Weibin, QU Jianbo, et al. Protective door damaged by air shock wave and fragment arisen from explosion in prototype tunnel[J]. Explosion and Shock Waves, 2004, 24(5):475-479. http://www.bzycj.cn/CN/abstract/abstract9989.shtml
[3]	LEPPÄNEN J. Experiments and numerical analyses of blast and fragment impacts on concrete[J]. International Journal of Impact Engineering, 2005, 31(7):843-860. doi: 10.1016/j.ijimpeng.2004.04.012
[4]	金乾坤.破片和冲击波毁伤圆柱靶的数值仿真[J].兵工学报, 2006, 27(2):215-218. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bgxb200602006 JIN Qiankun. Simulation of cylindrical shell damage by fragments and shock waves[J]. Acta Armamentarii, 2006, 27(2):215-218. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bgxb200602006
[5]	RAKVÅG K, UNDERWOOD N, SCHLEYER G, et al. Transient pressure loading of clamped metallic plates with pre-formed holes[J]. International Journal of Impact Engineering, 2013, 53(3):44-55. https://www.sciencedirect.com/science/article/pii/S0734743X12001662
[6]	张成亮, 朱锡, 侯海量, 等.爆炸冲击波与高速破片对夹层结构的联合毁伤效应试验研究[J].振动与冲击, 2014, 33(15):184-188. http://www.cnki.com.cn/Article/CJFDTotal-ZDCJ201415033.htm ZHANG Chengliang, ZHU Xi, HOU Hailiang, et al. Tests for combined damage effect of blast waves and high-velocity fragments on composite sandwich plates[J]. Journal of Vibration and Shock, 2014, 33(15):184-188. http://www.cnki.com.cn/Article/CJFDTotal-ZDCJ201415033.htm
[7]	KONG X S, WU W G, LI J, et al. Experimental and numerical investigation on a multi-layer protective structure under the synergistic effect of blast and fragment loadings[J]. International Journal of Impact Engineering, 2014, 65(2):146-162. https://www.researchgate.net/publication/259515939_Experimental_and_Numerical_Investigation_on_a_Multi-layer_Protective_Structure_under_the_Synergistic_Effect_of_Blast_and_Fragment_Loadings
[8]	蒋建伟, 侯俊亮, 门建兵, 等.爆炸冲击波作用下预制孔靶板塑性变形规律的研究[J].高压物理学报, 2014, 28(6):723-728. doi: 10.11858/gywlxb.2014.06.013 JIANG Jianwei, HOU Junliang, MEN Jianbing, et al. Study on deformation of perforated plates under blast loading[J]. Chinese Journal of High Pressure Physics, 2014, 28(6):723-728. doi: 10.11858/gywlxb.2014.06.013
[9]	杜志鹏, 李晓彬, 夏利娟, 等.反舰导弹攻击舰船舷侧防护结构过程数值仿真[J].哈尔滨工程大学学报, 2006, 27(4):484-487. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=hebg200604001&dbname=CJFD&dbcode=CJFQ DU Zhipeng, LI Xiaobin, XIA Lijuan, et al. Numerical simulation of anti-ship missile attack warship broadside pro-cess[J]. Journal of Harbin Engineering University, 2006, 27(4):484-487. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=hebg200604001&dbname=CJFD&dbcode=CJFQ
[10]	NYSTRÖM U, GYLLTOFT K. Numerical studies of the combined effects of blast and fragment loading[J]. International Journal of Impact Engineering, 2009, 36(8):995-1005. doi: 10.1016/j.ijimpeng.2009.02.008
[11]	HATCH-AGUILAR T, NAJJAR F, SZYMANSKI E. Computational hydrocode study of target damage due to frag-ment-blast impact [C]//26th International Symposium of Ballistics. Livermore, California, USA: Lawrence Livermore National Laboratory, 2011: 1-17. http://digital.library.unt.edu/ark:/67531/metadc840112/
[12]	ZHANG C Z, CHENG Y S, ZHANG P, et al. Numerical investigation of the response of I-core sandwich panels subjected to combined blast and fragment loading[J]. Engineering Structures, 2017, (151):459-471. https://www.researchgate.net/publication/319645640_Numerical_investigation_of_the_response_of_I-core_sandwich_panels_subjected_to_combined_blast_and_fragment_loading
[13]	段新峰, 程远胜, 张攀, 等.冲击波和破片联合作用下I型夹层板毁伤仿真[J].中国舰船研究, 2015, 10(6):45-59. doi: 10.3969/j.issn.1673-3185.2015.06.008 DUAN Xinfeng, CHENG Yuansheng, ZHANG Pan, et al. Numerical analysis of the damage on I-core sandwich panels subjected to combined blast and fragment loading[J]. Chinese Journal of Ship Research, 2015, 10(6):45-59. doi: 10.3969/j.issn.1673-3185.2015.06.008
[14]	李茂, 朱锡, 侯海量, 等.冲击波和高速破片对固支方板的联合作用数值模拟[J].中国舰船研究, 2015, 10(6):60-67. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgjcyj201506009 LI Mao, ZHU Xi, HOU Hailiang, et al. Numerical simulation of steel plates subjected to the impact of both impact waves and fragments[J]. Chinese Journal of Ship Research, 2015, 10(6):60-67. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgjcyj201506009
[15]	ZHANG P, CHENG Y S, LIU J, et al. Experimental study on the dynamic response of foam-filled corrugated core sandwich panels subjected to air blast loading[J]. Composites Part B: Engineering, 2016(105):67-81. https://www.deepdyve.com/lp/elsevier/experimental-study-on-the-dynamic-response-of-foam-filled-corrugated-A8Ed0bJ5WD
[16]	CHENG Y S, ZHOU T Y, WANG H, et al. Numerical investigation on the dynamic response of foam-filled corrugated core sandwich panels subjected to air blast loading[J]. Journal of Sandwich Structures & Materials, 2017. DOI: 10.1177/1099636217700350.
[17]	CHENG D S, HUNG C W, PI S J. Numerical simulation of near-field explosion[J]. Journal of Applied Science and Engineering, 2013, 16(1):61-67. https://www.researchgate.net/publication/288763805_Numerical_simulation_of_near-field_explosion
[18]	ZHANG P, CHENG Y S, LIU J, et al. Experimental and numerical investigations on laser-welded corrugated-core sandwich panels subjected to air blast loading[J]. Marine Structures, 2015, 40:225-246. doi: 10.1016/j.marstruc.2014.11.007
[19]	JING L, XI C Q, WANG Z H, et al. Energy absorption and failure mechanism of metallic cylindrical sandwich shells under impact loading[J]. Materials & Design, 2013, 52(24):470-480. https://www.sciencedirect.com/science/article/pii/S0261306913005220