靶板在爆炸成型弹丸垂直侵彻下的层裂*

李 睿,黄正祥,祖旭东,肖强强,贾 鑫

(南京理工大学机械工程学院,江苏 南京 210094)

层裂是靶板在弹丸冲击载荷作用下的一种重要破坏形式，也是靶后破片形成的组成部分。研究靶板在弹丸冲击载荷作用下的层裂效应对新型装甲设计和靶后破片等领域具有很强的指导意义。因此，准确掌握层裂的形成机理是十分必要的。对层裂效应已展开了大量的研究。Rinehart等[1]系统分析了材料的层裂现象，提出了最大拉应力层裂破坏准则。Rinehart[2]全面介绍了应力波在靶板材料内的传播、相互作用并产生层裂的过程。Ren等[3]提出了无网格数值模拟模型，并对钛铝合金在冲击载荷下的层裂情况进行了数值模拟，发现非弹性脉冲波之间的相互作用在层裂破坏机制中起重要作用。Yu等[4]利用数值模拟和实验相结合的方法研究了金属材料在冲击载荷下的多层层裂情况。杜忠华[5]研究了陶瓷靶板在冲击载荷下层裂的形成机理，并分析了不同的弹丸速度对靶板层裂的影响。陈大年等[6]在NAG模型基础上提出了一种基于空穴聚集的层裂模型。刘飞等[7]采用数值模拟和实验相结合的方法研究了钢板在接触爆炸载荷作用下的层裂效应，认为装药高度、钢板厚度和材质是影响钢板层裂的主要因素。魏波[8]采用SPH方法，开展数值模拟研究了不同飞片速度、不同靶板厚度对钢板和玻璃产生层裂的影响。目前对层裂效应已做了广泛的研究，但是对爆炸成型弹丸(explosively formed projectile,EFP)侵彻靶板产生层裂效应的研究主要集中在数值模拟和实验方面。本文中，基于波动力学和基本假设，从理论上对在EFP冲击载荷作用下靶板发生层裂的过程进行分析，在不考虑EFP强度的基础上建立相应的力学模型，研究靶板厚度和EFP速度对靶板层裂的影响，并通过实验验证理论模型的可靠性。

1 理论模型

1.1 应力波的形成和传播

在EFP垂直侵彻靶板的过程中，由于应力波的相互作用，在靶板背面产生层裂。由于EFP侵彻入孔情况与靶后效应无关，因此本文中暂不考虑，主要研究靶板背面层裂区域的形状和范围。为了便于公式推导和简化计算，作如下几点假设：(1)在侵彻靶板过程中将EFP当作圆柱杆；(2)应力波在靶板内以球面波的形式传播；(3)不考虑应力波的衰减；(4)靶板材料均匀无缺陷。根据适用于高速撞击的A-T模型，EFP侵彻靶板时的侵彻速度为：

(1)

式中：vj为EFP着靶速度，ρj为构成EFP材料的密度，ρt为靶板材料密度，Rt为靶板材料侵彻阻抗。

EFP在侵彻靶板时，在靶板惯性作用下，孔径开始增长，同时产生的应力波在靶板中以球面波的形式传播。孔径增长开始时孔壁的初始压力等于轴向压力，因此可以得到：

(2)

式中：p0为初始孔径增长压力。

径向压力pc与初始孔径增长压力的关系为:

(3)

式中：r0为EFP的初始半径，rc为侵彻孔半径。

由文献[9]可知，在EFP侵彻靶板过程中，随着孔径的增长，应力波以球形波在靶板中传播，不同半径r处的应力波强度为:

(4)

式中：Y为靶板的屈服强度，uc为孔径增长速度。

由文献[10]可知，不同孔径处的孔径增长速度为:

(5)

对式(5)对时间求导，则可得到:

(6)

由Szendrei-Held方程[11]可知侵彻孔径rc随时间t的变化为:

(7)

图1 入射波和反射波相互作用示意图Fig.1 Schematic diagram of interaction between incident and reflected waves

1.2 应力波的反射和层裂的形成

当应力波传播到靶板背面时发生斜反射，则反射波为:

σf=Rσt(r)

(8)

式中：R为反射系数[2]，且：

(9)

式中：α为入射角，β为反射横波的反射角。

EFP侵彻靶板过程中形成的应力波在靶板自由端反射后形成反射拉伸波，在反射波的传播过程中，入射波和反射波在靶内相互作用形成合应力波。应力波和反射拉伸波相互作用的几何关系如图1所示，t0时刻EFP侵彻到A0点，假设此时产生的应力波以入射角α传播到B点，在B点形成反射波并向前传播到某一点C，与t时间后侵彻界面A点产生的应力波在C点相互作用，形成拉伸应力波。按临界应力断裂准则，若拉伸应力波强度达到临界断裂应力σct时，即:

σt≥σct

(10)

则靶板背面出现层裂，在靶板内形成裂纹。

由文献[2]可知，两应力波斜交下的合应力强度为:

(11)

式中：σa为由A点产生的与反射拉伸波σf相互作用的应力波强度，ν为靶板材料泊松比。

由临界断裂准则可知，当该合应力波强度σt达到临界断裂应力σct时，发生层裂。由式(4)、式(10)和式(11)组成方程组，可确定该层裂点距初始入射波阵面的距离r。

1.3 层裂点位置计算

根据图1，在t0时刻EFP侵彻到A点，在入射波从A0点传播到B点再反射传播到某一点C的这段时间里，侵彻界面由A0点侵彻到A点，前进的距离为:

r2=uΔt

(12)

侵彻界面产生的应力波在自由端反射，A0点距反射点B的距离为:

r0=(h-ut)/cosα

(13)

A点距反射点B的距离为:

(14)

由几何关系可知:

(15)

侵彻界面到达A点的时间为:

t=t0+(r0+s)/vt

(16)

式中：vt为RHA材料中纵波的体积波速。

由式(12)～(16)可知，层裂点C距t时刻侵彻界面A点的相对位置为:

(17)

以靶板背面所在平面为横轴、靶板厚度方向为纵轴建立坐标系，坐标原点为侵彻孔径轴线与靶板背面的交点，则层裂点C的坐标可以表示为:

x=rsinθ,y=h-ut-rcosθ

(18)

1.4 理论结果分析

图2 不同EFP速度下靶板背面发生层裂的区域Fig.2 Spallation zones at target backs at different EFP velocities

基于建立的理论模型和Matlab软件，对EFP侵彻有限厚靶板的过程进行计算，EFP材料为高导无氧铜，密度ρj=8.96 g/cm3；轧制均质装甲[12](rolled homogeneous armor，RHA)密度ρt=7.86 g/cm3，屈服强度Y=1.5 GPa，泊松比为ν=0.28。计算过程中EFP半径为21.5 mm，垂直侵彻RHA靶板。

转染48 h后，流式细胞术检测各组SHG-44细胞凋亡情况，结果显示miR-543 mimic组与mimic NC组相比细胞凋亡率显著增加(P＜0.01);miR-543 inhibitor组与inhibitor NC组相比，细胞凋亡率显著减少(P＜0.01)，见表2和图4。由此可见，过表达miR-543可促进SHG-44细胞凋亡，抑制miR-543表达可抑制细胞凋亡。

1.4.1不同速度对靶板层裂的影响

根据理论分析，EFP侵彻靶板过程中产生的应力波的强度与EFP着靶速度密切相关，而靶板能否发生层裂取决于应力波和反射应力波相互作用后的拉应力波强度的大小，因此EFP的着靶速度是影响层裂的主要因素之一。在靶板厚度h=40 mm的情况下，EFP分别以1 550、1 600、1 700、1 750、1 800、1 850、1 900、1 950和2 000 m/s的着靶速度侵彻靶板，靶板背面发生层裂的区域见图2。

应力波以球面波的形式传入靶板，随着入射角的增大，到达自由表面的路程不断增加，反射拉应力波的强度随入射角的增大而减小。因此层裂的厚度随入射角的增大而增大。综合效果，使靶板背面出现弯月形的层裂区域，如图2所示。在计算条件下，当EFP的着靶速度不同时，弯月形层裂区域的形状不同。随着EFP着靶速度的增大，弯月形层裂区域的厚度不断减小，层裂区域的长度不断增大。EFP速度从1 550 m/s增大到1 900 m/s时，靶板背面层裂区厚度从21.8 mm减小到3.9 mm，弯月形层裂区长度从86.4 mm增大到133.4 mm。这是因为增大EFP着靶速度，而入射波到自由表面的衰减量不变，因此反射波强度增大，于是层裂时刻提前，靶板背面的层裂厚度减小。当速度大于1 900 m/s时，层裂区域的面积急剧衰减，并且层裂区域形状也不再是弯月形。当速度达到2 000 m/s时，层裂长度为12 mm，厚度为2 mm，此时对靶后效应影响甚微。

图3 不同厚度靶板背面发生层裂的区域Fig.3 Spallation zones at target backs with different thicknesses

1.4.2不同厚度对靶板层裂的影响

根据理论分析，在给定EFP参数的基础上，靶板厚度直接影响反射应力波的强度，因此靶板厚度是靶板发生层裂的重要影响因素。在EFP着靶速度为1 800 m/s情况下，分别侵彻30、35、40、45、50、55和60 mm厚的靶板，靶板背面发生层裂的区域如图3所示。

2 实验验证

2.1 实验布置

采用∅140 mm的无壳体成型装药。针对EFP在7倍炸高下穿透40 mm厚RHA靶板进行实验。聚能装药如图4所示，装药直径为140 mm，装药长度为102 mm，炸药为JH-2。药型罩材料为紫铜，变壁厚，采用8#电雷管起爆。EFP药型罩如图5所示，经数值计算，EFP半径为21.5 mm，长度为83.5 mm，速度为1 800 m/s。为了完整地回收靶后崩落的碎片，采用注满水的水箱来回收靶后破片；成型装药以垂直方式布设，以保证EFP垂直入射RHA靶板；靶板距水面1 000 mm，实验现场布置见图6。在上述实验情况下，共进行2发验证实验。

图4 成型装药Fig.4 Photograph of shaped charge

图5 EFP药型罩Fig.5 Diagram of EFP liner

图6 侵彻实验现场布置Fig.6 Layout of penetration experiment

2.2 实验结果分析

由于EFP侵彻靶板过程中在靶板背面形成弯月形的层裂区域，因此层裂的厚度不均匀，理论结果和实验结果如图7所示。分别测量图7中实验所得的层裂参数，层裂厚度取层裂区域最外侧处裂纹消失处的厚度，理论计算和实验结果见表1。

利用∅140 mm成型装药形成EFP对40 mm厚装甲钢靶板进行侵彻，实验后靶板发生层裂。由表中结果可知，理论计算结果与实验结果吻合较好，层裂的厚度和长度与实验1的误差分别为4.76%和2.89%；与实验2的误差分别为3.30%和1.81%。验证了EFP侵彻有限厚靶板形成层裂理论分析模型是正确的。

图7 实验与理论计算结果Fig.7 Experimental and theoretical results

层裂厚度/mm理论计算实验1实验 2层裂长度/mm理论计算实验1实验 28.88.49.1114.0110.8116.1

3 结 论

(1)利用应力波传播特性，理论分析了EFP侵彻有限厚靶板形成层裂的过程，得到了计算层裂点位置的表达式。

(2)编程计算了在靶板厚度h=40 mm不变的情况下，不同的EFP着靶速度vj对层裂的影响，计算结果表明，速度在1 500～1 900 m/s范围内，随着EFP着靶速度的增大，弯月形层裂区厚度不断减小，长度不断增大。EFP着靶速度大于1 900 m/s时，层裂区域形状不再呈弯月形，层裂区域急剧缩小。

(3)利用Matlab软件编程计算了在vj=1 800 m/s不变的情况下，不同靶板厚度对层裂的影响，计算结果表明，靶板厚度在35～60 mm范围内，随着靶板厚度的增大，弯月形层裂区厚度不断增大，长度不断减小。当靶板厚度小于35 mm时，层裂区域形状不再呈弯月形，层裂区域急剧缩小。

(4)根据EFP侵彻实验，得到在vj=1 800 m/s且h=40 mm情况下，靶板发生层裂时的参数。同时利用EFP垂直侵彻有限厚装甲钢靶板形成层裂理论对在相同条件下靶板层裂情况进行计算，两者结果吻合较好。