Research advances in circumferential multiple linear explosively-formed projectile technology
-
摘要: 传统的破片式防空反导战斗部爆炸后产生的破片杀伤元数量虽多，却不能有效击毁来袭的不敏感弹药，存在威力不足问题，因而限制了其发展。周向多线性爆炸成型弹丸(multiple linear explosively-formed projectile, MLEFP)战斗部爆炸后在周向产生多个高速、大质量、大长径比的对折型线性爆炸成型弹丸(linear explosively-formed projectile, LEFP)，具备击穿、击爆厚壁壳体不敏感弹药的能力，因此在中近程防空反导作战中具备广阔的应用前景。从线性毁伤元的发展和对折型LEFP的成型技术出发，重点分析了炸药装药、药型罩等关键部件影响线性毁伤元成型的研究成果，对比了3种毁伤元初速工程计算模型的理论依据、优缺点等，概括了近年来对折型LEFP侵彻试验结果，最后总结了周向MLEFP战斗部及其毁伤元未来的发展方向。Abstract: The traditional fragmental warhead used in air defense and antimissile can not effectively destroy the incoming insensitive ammunition due to the weak fragmentations, which limits its development. After explosion of the circumferential multiple linear explosively-formed projectile (MLEFP) warhead, a number of folded linear explosively formed projectiles (LEFP) are produced in the circumferential direction with high speed, large mass and large length to diameter ratio. These projectiles can penetrate thick cases and ignite insensitive ammunitions, so the MLEFP warhead has a great application prospect in the medium-short range air defense and anti-missile system. Based on the development of linear projectile and the forming method of new LEFPs, the influences of charge and liner are focused on. Theories, merits and demerits of three initial velocity models are compared. The results of penetration tests using folded LEFPs in recent years are summarized. Finally, the future development direction of circumferential MLEFP warhead and linear projectile is analyzed.
-
延性金属在强动载荷下的动态破坏是冲击动力学中重要的基础科学问题之一，随着国防建设、汽车工业以及航空航天等领域的发展，深入地了解延性金属在高速冲击下的动态损伤有着重要的意义[1]。层裂是延性金属动态损伤的一种典型表现形式，层裂实验中得到的自由面速度曲线包含了延性金属在动态破坏过程中的丰富信息[2-3]。在对延性金属层裂行为的研究中，平板撞击（planar plate impact, PPI）实验是目前主要采用的一种研究方法，实验过程中内部波系传播如图1（a）所示：飞片在t0时刻以一定速度撞击样片，二者相撞时会分别传入冲击波，当冲击波到达样片自由面和飞片后表面后会形成反射稀疏波，两束稀疏波会在样片内部相遇并产生拉伸区域，由此引起的内部损伤逐渐累积，直至材料失效断裂。Curran等[3]基于对大量实验结果的分析，指出延性金属的层裂是由于内部的微损伤形核、生长和聚集后，形成层裂飞片的复杂动力学过程。
图 1 平板撞击实验原理及自由面速度曲线示意图Figure 1. The principle of the flat plate impact experiment and a schematic diagram of a typical free-surface velocity curveAntoun等[4]对层裂的研究历史和现状作了权威性总结，指出了层裂研究还有许多问题尚待解决。Thomas等[5]采用对称碰撞的方式研究了铁在低于冲击相变阈值下的层裂行为，并对铁在低速冲击下的Hugoniot弹性极限与层裂强度进行了分析；陈永涛等[6]采用等厚对称碰撞实验，研究了纯铁在冲击压力超过相变阈值时的层裂行为，发现样品中出现了二次层裂现象；翟少栋等[7]通过轻气炮加载与激光加载两种方式研究了铝的层裂特性，结果表明铝的层裂强度具有明显的应变率效应；Koller等[8]通过改变冲击波波形，对无氧铜在不同冲击波形式下的层裂特性进行了研究，结果表明冲击波波形与拉伸脉冲的持续时间会显著影响无氧铜样片中的层裂损伤演化。实验研究在宏观尺度上提供了大量宝贵的实验数据，但由于其自身的特点，难以直接测量延性金属内部物理量随时间的变化数据，只能通过如自由面速度曲线等间接数据分析材料内部损伤演化过程。
随着计算机的发展与应用，数值模拟方法已经成为实验研究的有力补充。目前在冲击力学的数值模拟研究中主要有两种方法，网格法与无网格法，其中无网格法以光滑粒子流体动力学法（smootfied particle hydrodynamics SPH）[9]为代表得到了广泛应用。张凤国等[10]对延性金属在三角波加载下多次层裂破坏问题进行了数值模拟研究，对比分析了应力瞬时断裂准则和Tuler-Butcher准则的差异，结果表明影响层裂片厚度的主要因素是材料自身的物理性质和加载条件；种涛等[11]基于一维流体动力学数值模拟方法，分析了铁层裂损伤与冲击相变之间的相互影响；Glam等[12]采用数值模拟方法研究了温度效应对铝的层裂强度与Hugonoit弹性极限的影响，结果表明铝的初始温度与层裂强度正相关。Libersky等[13]率先在SPH方法中引入了材料强度效应；徐金宏等[14]利用一种改进SPH方法对靶板厚度、弹丸速度和形状对碎片云的影响进行了研究；Zhou等[15]在对铝球超高速撞击靶板的研究中，利用SPH方法成功地进行了数值模拟分析，结果与实验数据有较好的吻合性。贺年丰等[16]采用SPH方法对爆轰条件下锡的层裂行为进行了数值模拟，结果表明SPH方法得到的结果与实验数据有很好的一致性；可以看出，在宏观尺度上，数值模拟方法与实验有很好的吻合性，并且数值模拟的加载范围与条件更加灵活，计算快速、结果直观，可以获得材料宏观尺度上的动态破坏信息，但数值模拟结果的准确定取决于采用的受状态方程、强度模型和失效模型等因素，并且由于单元或粒子尺度影响，从微观尺度研究层裂机理的能力有限。
在微观尺度，分子动力学方法（molecular dynamics, MD）已经逐渐成为一种在微观尺度研究材料性能的主流方法。席涛等[17]通过分子动力学方法研究分析了超高应变率下铜的层裂特性，获得了铜在109 ～1010 s−1时的层裂强度；Rawat等[18]研究了应变率对单晶铁中的损伤演化的影响，结果表明随着应变率的增加，空洞更容易成核，但在低应变率条件下空洞的生长速度更快；Yang等[19]研究了应变率、温度、晶粒尺寸等因素对铝动态损伤发展的影响，揭示了不同影响因素下损伤发展的过程；Chen等[20]通过分子动力学方法研究了晶界特性对动态损伤的影响，结果表明损伤分布与晶界取向和空洞形核应力相关。可以看出：分子动力学模拟可以揭示动态损伤的微观机理，给出材料内部应力与损伤演化信息，有助于从根本上了解金属材料的动态损伤机制，但受到计算能力的限制，无法在全尺寸上复现层裂现象[21-23]。
典型自由面速度曲线如图1（b）所示，从中可以分析得到延性金属材料动态损伤的宏观响应的一系列重要信息，如加载应力幅值、拉伸应变率、层裂强度和层裂片厚度等，因此自由面速度曲线是研究层裂演化过程的最重要参考数据之一[24]。实际上延性金属的层裂过程十分复杂，是一种内部损伤随时空演化的动力学过程，材料从受到拉伸应力影响开始直至完全破坏要经历一个损伤累积的过程，期间会伴随能量的耗散和性能的退化，这都会影响材料内部波的传播过程，从而影响自由面速度曲线的形式。因此，自由面速度曲线是微观上的损伤演化过程和宏观响应场之间相互作用的一种宏观表象，由于宏观尺度上的实验和数值模拟很难得到反应内部损伤演化的信息，这对充分认识宏观响应与微观损伤发展之间的关系造成了阻碍，也就导致对自由面速度曲线中的一些特征理解存在争议。裴晓阳等[2]基于理论分析和微观损伤演化动力学理论，对延性金属层裂自由面速度曲线的某些特征所对应的物理内涵进行了阐述，定性地指出了宏观上自由面速度曲线的响应特征与内部损伤演化的对应关系，但是并没有从微观尺度上的给出材料层裂过程中的损伤演化过程，缺少在微观尺度的机理分析与验证。
目前的研究大多只局限于单一尺度下层裂现象模拟或微观损伤演化分析，少有结合宏观与微观尺度下对层裂自由面速度曲线特征的多尺度视角下的研究，缺乏对自由面速度曲线的微观视角解析。本文以延性金属钽为研究对象，在宏观尺度上对钽在高速冲击下的层裂行为进行数值模拟，通过与实验数据的比较，验证模型的合理性；在微观制度上通过分子动力学方法分析钽在层裂过程中的损伤演化过程。进而将宏观尺度与微观尺度的模拟结果有机结合起来，从多尺度视角对延性金属钽的层裂行进行研究。揭示了自由面速度曲线上的典型特征所蕴含的损伤演化规律和及物理内涵，研究结果可以为全面了解层裂这一复杂的动力学行为提供参考。
1. 计算模型与方法
为了便于与实验数据进行对比验证，宏观尺度全尺寸模型参照平板撞击实验建立[25]，选取样片自由面中心位置数据绘制自由面速度曲线，如图2（a）所示；微观尺度分子动力学模型选取样片中心拉伸区域进行模拟，MD模型如图2（b）所示。
1.1 宏观尺度数值模拟
1.1.1 状态方程与材料本构
Mie-Grüneisen状态方程能够较为准确地描述金属在冲击载荷下的动力学行为，其表达式如下[16]：
p=ρ0c20μ(1+μ)[1−(S1−1)μ]+γρ0e (1) μ=ρρ0−1 (2) 式中：ρ为冲击后材料密度，ρ0为初始材料密度，S1为材料常数，γ为Grüneisen系数，e为质量内能。Mie-Grüneisen状态方程参数如表1所示[26]。
表 1 Mie-Grüneisen状态方程参数Table 1. Parameters for Mie-Grüneisen equation of state材料 ρ0 /(kg·m−3) c0/(m·s−1) S1 γ Ta 16690 3340 1.20 1.67 考虑了三种常用的描述延性金属在强动载荷下的动态破坏行为的模型，分别是Johnson-Cook（JC）模型[27]、Steinberg-Cochran-Guinan（SCG）模型[28]以及Zerilli-Armstrong（ZA）模型[29]。
（1） JC模型
JC模型假设材料各向同性，综合考虑了应变、应变率硬化和温升软化等因素，并且将应变、应变率硬化效应和温升软化效应解耦，JC模型的表达式为：
σ=(A+Bεn)(1+Cln˙ε˙ε0)(1−T∗m) (3) T∗=T−TrTm−Tr (4) 式中：A、B、C、n、m为材料常数，Tr为参考温度，Tm为材料融化温度，
˙ε0 为参考应变率。JC模型参数如表2所示[30]。表 2 Johnson-Cook模型参数Table 2. Parameters for the Johnson-Cook model材料 A/MPa B/MPa n C m Ta 142 164 0.31 0.057 0.88 （2） ZA模型
基于热激活和位错动力学理论，Zerilli等[29]提出了一种描述体心立方金属（BCC）的强度模型，考虑了温度、应变率与温度的影响，ZA模型的表达式为：
σ=C0+k1λ−1/2+C2exp(−C3T+C4Tln˙ε)+C5εn (5) 式中：C0、k1、C2、C3、C4、C5为材料常数，T为绝对温度，ε为塑性应变，λ为晶粒尺寸。ZA模型参数如表3所示[29]。
表 3 Zerilli-Armstrong模型参数Table 3. Parameters for the Zerilli-Armstrong model材料 C0/MPa k1/（MPa·m−3/2） C2/MPa C3/K−1 C4/K−1 C5/MPa n Ta 1125 10 178 5.35×10−3 0.327×10−3 310 0.44 （3） SCG模型
Steinberg等[28]提出的SCG模型考虑了静水压力对剪切模量和屈服强度的强化效应以及温度软化效应，同时忽略了应变率效应，并假设金属在融化状态下可以忽略材料强度。屈服强度Y和剪切模量G分别为：
G=G0[1+(G′pG0)pη1/3+(G′TG0)(T−300)] (6) Y=Y0(1+βεp)nGG0 (7) 式中：η为压缩比，εp为等效塑性应变，T为材料温度，p为压力，G0和Y0分别为初始状态下的剪切模量和屈服强度（T=300 K，p=0，ε=0），
G′p 、G′T 、β和n为材料常数。Y满足：
Y0(1+βεp)n≤Ymax (8) 式中：Ymax为最大屈服强度。SCG模型参数如表4所示[28]。
表 4 Steinberg-Cochran-Guinan模型参数Table 4. Parameters for the Steinberg-Cochran-Guinan model材料 G0/GPa Y0/GPa Ymax/GPa β n G′p G′T/（MPa ·K−1） Tm0/K Ta 69 0.77 1.10 10 0.1 1.005 −8.97 4340 1.1.2 层裂模型
Grady[31]基于断裂力学理论提出了层裂的能量平衡破碎模型，将层裂强度ps与层裂碎片平均尺寸和应变率联系起来。对于延性金属，层裂强度随冲击加载过程而改变，表达式为[32]：
ps=(2ρc20Yεc)1/2 (9) 式中：ρ为材料密度，Y为屈服强度，c0为材料体积声速，εc为材料临界失效应变，对金属材料一般可取为0.15[31]。
1.2 自由面速度曲线分析
层裂强度是表征材料在强动载荷下性能的一个重要指标，在材料动态损伤研究中具有重要参考意义。Novikov[33]基于声学近似的方法给出了计算层裂强度的公式：
σspall=12ρ0cbΔus (10) 式中：ρ0为材料初始密度，cb为材料初始体声速，Δus为自由面速度曲线最大值与第一次速度拉回时的差值。
层裂区域平均拉伸应变率
˙εs 的计算公式为[34]：˙εs=ΔusΔts12cb (11) 式中：Δus与Δts分别为自由面上速度最大值和第一个极小值之间的速度差和时间差。
自由面速度回跳速率，即极小值到峰值速度之间的斜率，计算公式为：
˙εr=ΔurΔtr12cb (12) 式中：Δur与Δtr分别为自由面上速度第一个极小值与其后峰值之间的速度差和时间差。
层裂片厚度h为：
h=Δt⋅cl2 (13) 式中：Δt为第一个层裂振荡周期持续时间，cl为弹性纵波声速。
1.3 微观尺度MD模型
钽原子间的相互作用采用Ravelo等[35]提出的嵌入原子势函数，MD模拟基于Lammps软件[36]进行。MD模型尺寸为16.53 nm×16.53 nm×16.53 nm，共250000个原子，在模型x、y、z方向施加周期性边界条件。加载之前采用NPT系综对模型进行30 ps的弛豫，确保体系在施加荷载之前处于平衡状态；加载过程采用NVE系综，在x、y、z方向施加相等的应变率，为了保持与宏观尺度模拟的一致性，施加应变率数值与宏观尺度模型中的拉伸应变率相同。
2. 宏观尺度模型验证
分别采用Lagrange网格法与SPH无网格法结合三种强度模型对钽的层裂进行了模拟，并对二者得到的结果进行了比较分析。表5中列出了用于验证的模型尺寸参数、撞击速度、强度模型等信息。模型为二维轴对称，Lagrange网格尺寸为0.05 mm，SPH粒子尺寸为0.1 mm，样片直径为50 mm，撞击速度为306 m/s。
表 5 用于验证的模型编号及参数设置Table 5. Model number and parameter settings for validation模型编号 飞片厚度/mm 样片厚度/mm 强度模型 方法 V-01 3 4.95 JC Lagrange V-02 3 4.95 JC SPH V-03 3 4.95 ZA Lagrange V-04 3 4.95 ZA SPH V-05 3 4.95 SCG Lagrange V-06 3 4.95 SCG SPH 不同模型得到的自由面速度曲线如图3所示。从模拟结果与实验数据对比可以发现，Lagrange网格法与SPH无网格法在数值模拟中各有优势，在0～2 μs时，样片自由面速度曲线速度由0 m/s增加到最大值，Lagrange方法得到的模拟结果与实验数据有较好的吻合性，并且Lagrange方法得到的结果可以观察到明显的Hugoniot弹性极限信号，而由SPH方法得到的结果Hugoniot弹性极限信号不是特别明显；在自由面速度开始下降后，Lagrange方法得到的模拟结果中自由面速度回拉速率比实验数据更快，SPH方法得到的结果与实验数据较为吻合。结果表明采用SPH无网格法在钽层裂模拟中是可行的，结合两种方法可以得到完整的自由面速度曲线数据。此外，结果表明SCG模型在钽的层裂模拟中更为适合。
为了分析不同加载应变率下的层裂特性，通过改变加载条件改变了层裂区域的拉伸应变率，表6给出了模型加载条件参数，拉伸应变率范围为2.13×104～5.40×104 s−1，压力范围为6.19～12.25 GPa。不同拉伸应变率下的钽层裂宏观尺度数值模拟结果如图4所示，从图中可以看出，自由面最大速度与飞片撞击速度基本相同；在相同的加载速度下，通过改变飞片厚度自由面速度平台宽度也随之改变，二者之间呈现正比关系，飞片厚度越大，自由面速度曲线平台宽度越宽。
表 6 不同应变率下平面撞击层裂模型参数与结果Table 6. Parameters of planar plate impact simulations and results under various strain rates模型编号 飞片厚度/mm 样片厚度/mm 加载速度/(m·s−1) p/GPa ˙εs/s−1 σspall/GPa ˙εr/s−1 S-01 2 4.95 306 8.84 5.40×104 4.92 3.57×104 S-02 2 4.95 250 7.05 4.69×104 4.70 3.25×104 S-03 3 4.95 410 12.25 3.92×104 4.14 2.32×104 S-04 3 4.95 306 8.84 3.28×104 3.97 1.74×104 S-05 3 4.95 210 6.19 2.68×104 3.71 1.69×104 S-06 4 4.95 306 8.84 2.31×104 3.34 1.34×104 
图 4 宏观尺度下不同拉伸应变率自由面速度曲线Figure 4. Free surface velocity under various tensile strain rates at macro-scale3. 结果分析
3.1 层裂强度
层裂强度是金属材料动态损伤行为的重要参数之一，通过对模拟结果的分析同时参照文献[3]给出的铜层裂强度与拉伸应变率的关系
σspall=0.15˙ε0.2 ，我们得到了钽的层裂强度与拉伸应变率之间的关系：σspall=0.042˙ε0.44 (14) 图5给出了钽的层裂强度与拉伸应变率的关系，从图中可以发现，层裂强度不是一个仅受到加载速度影响的物理量，在加载速度相同的条件下，S-06、S-04、S-01的层裂强度随着拉伸应变率的升高而增加，对于S-05、S-03、S-02尽管加载速度相差较大，但是同样表现出层裂强度与拉伸应变率之间的正比关系。这表明层裂强度呈现出一种率相关特性，即随拉伸应变率的增长而增加。
图 5 层裂强度与拉伸应变率关系Figure 5. Relationship between spall strength and tensile strain rate图6给出了对数坐标下
σspall 与˙εs 之间的关系，从图中可以看出在对数坐标下层裂强度表现出了明显的率相关特性。我们对比了更高拉伸应变率下的实验数据[37]，更高的拉伸应变率下层裂强度同样是随应变率提高而增加，同时这也验证了数值模拟结果的准确性。
图 6 对数坐标下层裂强度与拉伸应变率的关系Figure 6. Relationship between spall strength and tensile strain rate in logarithmic coordinates前文分析中的层裂强度是根据式（11）进行计算得到的，也是目前研究中广泛采用的声学近似计算方法，但并不是唯一的计算方法。Stepanov[33]指出，对于延性金属，平板撞击过程中冲击应以弹性纵波速度cl传播。而在它前面的入射稀疏塑性波应以体积声速cb传播，当重点考虑弹性波的影响时，有效声速ce为：
ce=2clcb/(cl+cb) (15) 据此，可得到层裂强度的一种修正形式：
σ(1)spall=ρ0cbΔu(11+cb/cl) (16) 式中：
ρ0 为材料密度，Δu 为自由面上速度最大值和第一个极小值之间的速度差。除此之外，Kanel[38]指出在层裂强度计算中还应考虑层裂片厚度影响，当计及层裂片厚度影响时：σ(2)spall=12ρ0cbΔu+δ2 (17) δ2=(hspcb−hspcl)⋅|˙u1˙u2||˙u1|+˙u2⋅12ρ0cb (18) 式中：
˙u1 为样品自由面速度-时间曲线上Pullback信号出现前的速度变化率，˙u2 为层裂反弹信号的速度变化率，hsp为层裂片厚度。由3种公式计算得到的层裂强度如表7所示，分析发现式（17）与式（11）得到的结果最大相差8%，式（18）与式（17）得到的结果相差不大。表明通过自由面速度曲线来计算层裂强度与选用的计算方法相关，而且不同计算模型之间存在一定差异，在分析时应该得到关注。表 7 不同计算公式得到的层裂强度Table 7. Spall strengths obtained by different formulas模型编号 σspall/GPa σ(1)spall/GPa σ(2)spall/GPa S-01 4.92 5.25 5.36 S-02 4.71 5.02 5.32 S-03 4.13 4.40 4.48 S-04 3.97 4.23 4.32 S-05 3.72 3.96 4.07 S-06 3.35 3.57 3.58 3.2 自由面速度曲线特征多尺度分析
从宏观尺度下层裂特性的分析可以看出，自由面速度曲线反映了层裂过程中材料内部损伤演化与宏观响应场的相互作用，全面和准确地认识自由面速度曲线对研究层裂现象具有重要的意义。基于目前的研究手段，对层裂过程中损伤演化的原位观测依然难以实现，无法直观的准确的了解层裂过程中材料的宏观动力学响应与微结构演化之间的关系，进而导致对自由面速度曲线的理解也存在一些争议，例如：Pullback信号的意义、损伤演化的形式、回跳速率的意义等。针对以上几个方面的问题，通过分子动力学方法给出层裂过程中损伤演化的情况，定性的从微观尺度对自由面速度曲线上的一些典型特征进行解读。
3.2.1 Pullback信号
Pullback信号是用来判断层裂是否发生的重要依据，但此时材料是否完全分离尚存争议。Zurek等[39]认为Pullback信号出现时材料已经完全分离，但Kanel等[40]等观察到在层裂面损伤很小时，依然有Pullback信号出现；Bonora等[2]认为Pullback信号的出现与微孔洞成核相关，但未给出验证。下文以S-04加载情况为例，对其自由面速度曲线中Pullback信号意义从微观尺度进行了分析。
图7～8给出了自由面速度下降后的时程曲线与层裂区MD模型内部应力及损伤演化情况，在分析中考虑了自由面速度曲线中Pullback信号在选取时的误差，将Pullback信号出现时间扩展为0.05 μs时间范围进行分析。从图7中可以看出，在Pullback出现的时间范围内，层裂区域内的应力达到峰值，在经过Pullback信号后，层裂区域内的应力迅速下降，出现应力松弛；对比此时层裂区域内的损伤发展可以发现，Pullback信号出现的时间范围对应的是损伤发展的初始阶段，也就是微孔洞成核阶段。在图8中可以直观的观察到内部损伤演化的情况，可以看出，在Pullback信号时间内微孔洞开始逐渐产生，此时内部孔洞体积占比非常小，这都表明钽在Pullback信号出现的时间范围内层裂区域实际上并未完全分离，而且Pullback信号对应的实际上是内部微孔洞形核的过程，从微观尺度证明了Bonora等[2]的推测。此外，在通过自由面速度曲线计算层裂强度时，Pullback对应的速度降低幅值（图1b中Δu）是一个重要的依据，通过前面的分析表明在钽的层裂中自由面速度曲线Pullback信号对应的是微孔洞形核的过程，这表明Δu实际上反映了了微孔洞形核的条件，那么通过Δu得到的层裂强度更准确地说应该是抵抗损伤起源的强度或微孔洞成核的一种表征。
图 7 Pullback信号及其对应的MD模型应力与空洞演化情况Figure 7. Pullback signal and the stress and void volume evolution in the MD model
图 8 Pullback信号持续时间Δt内空洞演化情况Figure 8. Void evolution during Pullback signal duration Δt层裂区域的损伤演化在宏观尺度的试验下难以直接观测，借助微观尺度的MD模型，可以从微观尺度对钽层裂过程中的损伤演化进程进行分析。在图9～11中我们给出了损伤演化过程中内部孔洞所占体积分数、微孔洞数量与内部孔洞情况。从图9中可以看出层裂区域内微损伤的演化可以分为3个阶段：孔洞形核（S1）、长大（S2）与聚集（S3）。在内部压力达到峰值附近时，微孔洞开始形核，此时孔洞数量快速增加，但是孔洞体积占比非常小，在S1阶段空洞数量演化情况如图10所示；随后进入孔洞长大阶段，此时应力开始快速下降，而应力松弛意味着出现能量耗散，耗散的能量为孔洞长大提供动力，孔洞的快速增长导致孔洞占比迅速升高，孔洞数量在这个阶段开始降低，这表明在S2阶段中孔洞体积的增加是由于孔洞长大所致，而不是新孔洞出现造成的；在应力下降到稳定状态后，孔洞开始进入聚集阶段，此时孔洞体积占比增长速度降低，并且孔洞数量趋于稳定，这表明在损伤发展后期是层裂区域内大孔洞再聚集最终导致层裂的过程，S2阶段的空洞长大与S3阶段的空洞聚集过程如图11所示。
图 9 空洞体积演化与应力及空洞数量的关系Figure 9. Relationship of the evolution of void volume with stress and void numbers
图 11 阶段2与3的空洞长大与聚集情况Figure 11. Void growing and coalescence during stage 2 and stage 33.2.2 回跳速率
在以往的分析中，往往只是关注Pullback信号及其前段自由面速度曲线，对Pullback信号出现后的自由面速度回跳曲线少有涉及。以自由面速度曲线中Pullback信号为起点，分析了自由面速度回跳曲线的变化情况，不同应变率条件下的自由面速度回跳曲线如图12所示。从图中可以看出，自由面速度回跳曲线随加载速度和飞片厚度没有明显的变化趋势，但是随着拉伸应变率的增加，曲线逐渐变得陡峭，斜率逐渐增加。
图13给出了层裂强度与回跳速率之间的关系。可以看出，层裂强度与回跳速率之间近似呈线性增长的关系，而层裂强度具有率相关特性，这表明自由面速度回跳曲线也与应变率相关。
自由面速度回跳与拉伸应变率的关系如图14所示。从图14中可以看出，自由面速度回跳速率随着拉伸应变率的升高而增长，二者近似呈线性变化趋势。因此，可以推断钽在层裂过程中样片内部的损伤演化速率也是随拉伸应变率的增长而增加。
图 14 不同拉伸应变率下自由面速度回跳速率Figure 14. Free surface bounce rate under various tensile strain rates从图9中发现Pullback信号后自由面速度回跳时间段内层裂区域内损伤快速增长的阶段，微孔洞体积占比迅速增加，这表明自由面速度回跳与损伤演化之间存在对应关系。Kanel等[40]基于特征线方法分析了回跳速率和损伤演化速率的关系，指出二者之间是正比关系。图15给出了层裂区域内不同应变率下应力演化的情况，不同应变率下损伤演化的情况如图16所示。
从前文的分析中，我们发现回跳速率与拉伸应变率二者近似呈线性变化趋势，拉伸应变率越高，回跳速率越大。不同应变率下应力演化情况如图15所示，从宏观上波系分析可知自由面速度回跳速率与应力松弛速度相关，从层裂区域中的应力变化情况可以发现，应力下降速度随应变率的升高而加快，这表明回跳速率与应力松弛速度是呈正相关的。应力松弛速度越快，代表能量耗散的速度越快，而能量耗散的速度与损伤演化的速度相关。图16给出了不同应变率下损伤演化的情况，可以明显看出随着应变率的增加，微孔洞在长大阶段的速度越来越快。综上所述，自由面速度回跳速率是层裂区域损伤演化速率的宏观表征。
4. 结 论
本文以延性金属钽为研究对象，通过改变飞片厚度与加载速度改变加载应变率，对不同加载条件下的自由面速度曲线特性进行了分析。重要的是，结合宏观和微观尺度视角揭示了自由面速度曲线与微损伤演化之间的联系，特别是对自由面速度曲线上的一些典型特征的物理内涵从微观尺度进行了解读，结果表明：
（1）对钽层裂的宏观尺度模拟中，Lagrange方法在模拟钽的弹塑性转变过程有较好的表现，能够观测到明显的Hugoniot弹性极限信号，SPH方法得到的结果与实验数据更为符合，能够较为准确地描述层裂全过程中的自由面速度曲线变化情况，证明SPH方法研究钽的层裂行为是可行的；
（2）在宏观尺度上，钽的层裂强度随加载应变率的升高而增加，二者近似呈指数关系；分析表明，层裂强度与回跳速率之间存在正比关系，表明自由面速度回跳速率也是层裂行为的一种反映，应该得到关注；通过观察发现了回跳速率与拉伸应变率之间近似为线性关系，随拉伸应变率的升高而增长。
（3）微观尺度分析表明，自由面速度曲线是层裂区域内的微孔洞形核、长大、聚集直至完全断裂的微损伤演化过程的宏观反映。Pullback信号的临界点是微孔洞的形核在宏观上的表征，自由面速度曲线在Pullback信号前的下降幅值代表了微孔洞形核的条件，而由此计算得到的层裂强度反映了微孔洞形核强度，Pullback信号后速度回跳速率反映了层裂区域内损伤演化的速率。
-
图 6 不同成型装药形成不同毁伤元分类
Figure 6. Classification of types of projectiles by different shaped charges
图 12 单方向曲率线性圆缺药型罩示意图
Figure 12. Schematic diagram of a linear circular charge liner with unidirectional curvature
方法 LEFP成型速度/（m·s−1） 误差/% Gurney理论 1 303 9 数值模拟 1 490 1.6 静爆试验 1 433 / 
下载: 导出CSV
周向MLEFP战斗部 成型形态 侵彻后钢靶剖面 无外衬 
有外衬 
下载: 导出CSV
表 3 3种LEFP初速模型对比
Table 3. Comparison among the three models for calculating the initial velocity of an LEFP
初速模型 基础理论 最大误差/% 优点 缺点 公式（2） Gurney方法，
动量守恒12.6 形式简单，
使用方便简单合计不同组分误差大，
需增加修正系数公式（3） C-J爆轰理论，
动量守恒4 形式简洁，
精度高适用范围窄 公式（5） 有效装药理论，
Gurney方法6 精度较好，
适用范围广需计算获取大量初速数据，
使用不便
下载: 导出CSV
-
[1] ROSLUND L A. Initiation of warhead fragments I: Normal impacts: NOLTR 73-124 [R]. White Oak, MD: Naval Surface Weapons Center, 1973. [2] HELD M. Initiation criteria of high explosives at different projectile or jet densities [J]. Propellants, Explosives, Pyrotechnics, 1996, 21(5): 235–237. DOI: 10.1002/prep.19960210505. [3] 童宗保, 王金相, 彭楚才, 等. 预制破片对屏蔽炸药冲击引爆研究 [J]. 科学技术与工程, 2014, 14(7): 173–177. DOI: 10.3969/j.issn.1671-1815.2014.07.038.TONG Z B, WANG J X, PENG C C, et al. Study on the initiation of shielded explosive impacted by prefabricated fragment [J]. Science Technology and Engineering, 2014, 14(7): 173–177. DOI: 10.3969/j.issn.1671-1815.2014.07.038. [4] 魏继锋, 魏锦, 王树山. 离散杆战斗部技术研究 [J]. 战术导弹技术, 2015(2): 101–105. DOI: 10.16358/j.issn.1009-1300.2015.02.18.WEI J F, WEI J, WANG S S. Research on technology of discrete rod warhead [J]. Tactical Missile Technology, 2015(2): 101–105. DOI: 10.16358/j.issn.1009-1300.2015.02.18. [5] CEDER R, LEUS V, DAVID I. Exploring the v2d initiation criterion within the Lee-Tarver model [C]// Proceedings of the 30th International Symposium on Ballistics. Long Beach, 2017: 1887−1896. DOI: 10.12783/ballistics2017/16966. [6] 苗润源, 梁争峰. 防空反导用多爆炸成型弹丸技术现状与发展 [J]. 飞航导弹, 2017(7): 89–93. DOI: 10.16338/j.issn.1009-1319.2017.07.19.MIAO R Y, LIANG Z F. Current situation and development of multiple explosive shaped projectile technology for air defense and antimissile [J]. Aerodynamic Missile Journal, 2017(7): 89–93. DOI: 10.16338/j.issn.1009-1319.2017.07.19. [7] 奥尔连科. 爆炸物理学[M]. 孙承纬, 译. 北京: 科学出版社, 2011 [8] AHMED M, MALIK A Q. A review of works on shaped charges [J]. Engineering, Technology and Applied Science Research, 2017, 7(5): 2098–2103. DOI: 10.48084/etasr.1532. [9] WEIMANN K. Research and development in the area of explosively formed projectiles charge technology [J]. Propellants, Explosives, Pyrotechnics, 1993, 18(5): 294–298. DOI: 10.1002/prep.19930180511. [10] 松全才, 金韶华, 李文. 俄罗斯切割索的现状 [J]. 爆炸与冲击, 1997, 17(4): 382–388.SONG Q C, JIN S H, LI W. Manufacture of cutting explosives in Russia [J]. Explosion and Shock Waves, 1997, 17(4): 382–388. [11] HAYES G A. Linear shaped-charge (LSC) collapse model [J]. Journal of Materials Science, 1984, 19(9): 3049–3058. DOI: 10.1007/BF01026984. [12] SEOKBIN L, PAUL W. An investigation of the characteristics of linear shaped charges used in demolition [C]// Proceedings of the 29th Annual Conference on Explosives and Blasting Technique, 2003: 297−306. [13] CHASE J B, KUKLO R M, SHAW L L, et al. High resolution diagnostics of a linear shaped charge jet [C]// 18th International Symposium and Exhibition on Ballistics. San Antonio: Lawrence Livermore National Lab, 1999, 1: 442–448. [14] 曾新吾, 薛鸿陆. 线型聚能装药的理论研究 [J]. 爆炸与冲击, 1988, 8(2): 97–105.ZENG X W, XUE H L. A theoretical study on linear shaped charge [J]. Explosion and Shock Waves, 1988, 8(2): 97–105. [15] 李彬峰, 潘国斌. 聚能罩壁厚对切割深度的影响 [J]. 爆破器材, 1999, 28(6): 16–19. DOI: 10.3969/j.issn.1001-8352.1999.06.005.LI B F, PAN G B. A study of the influence of wall thickness of shaped charge liner on the depth of cutting [J]. Explosive Materials, 1999, 28(6): 16–19. DOI: 10.3969/j.issn.1001-8352.1999.06.005. [16] 祝逢春, 邓振礼, 胡瑜. 线性聚能装药切割航空炸弹可靠性评估 [J]. 火工品, 2000, 11(2): 24–26. DOI: 10.3969/j.issn.1003-1480.2000.02.007.ZHU F C, DENG Z L, HU Y. Reliability assessment on cutting aerial bombs using linear shaped charge [J]. Initiators and Pyrotechnics, 2000, 11(2): 24–26. DOI: 10.3969/j.issn.1003-1480.2000.02.007. [17] 苟瑞君, 赵国志, 杜忠华, 等. 线性成型装药的威力研究 [J]. 爆破器材, 2005, 34(5): 25–28. DOI: 10.3969/j.issn.1001-8352.2005.05.009.GOU R J, ZHAO G Z, DU Z H, et al. Study on the power of linear shaped charge [J]. Explosive Materials, 2005, 34(5): 25–28. DOI: 10.3969/j.issn.1001-8352.2005.05.009. [18] 苟瑞君, 赵国志. 线性成型装药起爆方式的比较 [J]. 火工品, 2006(1): 42–45. DOI: 10.3969/j.issn.1003-1480.2006.01.013.GOU R J, ZHAO G Z. Comparison of linear shaped charge initiation manners [J]. Initiators and Pyrotechnics, 2006(1): 42–45. DOI: 10.3969/j.issn.1003-1480.2006.01.013. [19] 陶钢, 朱鹤荣, 石连捷, 等. 爆炸成型弹丸药型罩结构分析 [J]. 弹道学报, 1995, 7(3): 84–86,93.TAO G, ZHU H R, SHI L J, et al. On the analysis of the shape for the EFP’s liner [J]. Journal of Ballistics, 1995, 7(3): 84–86,93. [20] 苟瑞君. 线性爆炸成型侵彻体形成机理研究[D]. 南京: 南京理工大学, 2006GOU R J. Forming mechanism of linear explosively formed penetrator [D]. Nanjing: Nanjing University of Science and Technology, 2006. [21] 杜忠华, 段卫毅. 起爆方式对LEFP成型及侵彻影响的数值模拟研究 [J]. 南京理工大学学报(自然科学版), 2009, 33(6): 756–759. DOI: 10.3969/j.issn.1005-9830.2009.06.009.DU Z H, DUAN W Y. Numerical simulation on formation and penetration effect of LEFP by blast ways [J]. Journal of Nanjing University of Science and Technology (Natural Science), 2009, 33(6): 756–759. DOI: 10.3969/j.issn.1005-9830.2009.06.009. [22] 段卫毅. 线性爆炸成型侵彻体成型机理与侵彻研究[D]. 南京: 南京理工大学, 2009DUAN W Y. Study on the formation mechanism and penetration of linear explosively formed penetrators [D]. Nanjing: Nanjing University of Science and Technology, 2009. [23] 高接东. 基于变壁厚药型罩的LEFP成型机理研究[D]. 南京: 南京理工大学, 2013GAO J D. Study on LEFP forming mechanism based on variable wall thickness liner [D]. Nanjing: Nanjing University of Science and Technology, 2013. [24] JOO J, CHOI J, KOO M, et al. Numerical analysis of penetration reduction of a long rod penetrator impacted by a linear explosively formed penetrator [C]// Proceedings of the 30th International Symposium on Ballistics. Long Beach: Destech, 2017: 1887−1896. [25] JOO J, CHOI J, KOO M, et al. Numerical analysis on penetration reduction of a WHA penetrator by an impact of linear explosively formed penetrator (LEFP) [J]. Journal of the Korea Institute of Military Science and Technology, 2017, 20(3): 384–392. [26] LI Y S, HUANG Z H, SHI A S, et al. Deformation and fracture failure of a high-speed long rod intercepted by linear explosively formed penetrators sequence [J]. Materials, 2020, 13(22): 5086. DOI: 10.3390/ma13225086. [27] 周涛, 王康康, 杜忠华, 等. 起爆方式对线性成型装药爆炸威力的影响 [J]. 火炸药学报, 2014, 37(2): 37–42. DOI: 10.3969/j.issn.1007-7812.2014.02.008.ZHOU T, WANG K K, DU Z H, et al. Effect of initiation way on the blasting power of linear shaped charge [J]. Chinese Journal of Explosives and Propellants, 2014, 37(2): 37–42. DOI: 10.3969/j.issn.1007-7812.2014.02.008. [28] 刘杰, 杜忠华, 王锋, 等. 基于爆轰波碰撞形成LEFP的研究 [J]. 工程力学, 2014, 31(8): 235–242. DOI: 10.6052/j.issn.1000-4750.2013.03.0198.LIU J, DU Z H, WANG F, et al. Linear explosively-formed penetrators based on detonation wave collision [J]. Engineering Mechanics, 2014, 31(8): 235–242. DOI: 10.6052/j.issn.1000-4750.2013.03.0198. [29] 王锋. 基于爆轰波对撞原理的新型EFP研究[D]. 南京: 南京理工大学, 2012WANG F. Research on new EFP based on detonation wave collision principle [D]. Nanjing: Nanjing University of Science and Technology, 2012 [30] CAO Y W, SUN X Y, YUAN B H, et al. Study on design and performance of linear EFP warhead [C]// Proceedings of the 28th International Symposium on Ballistics. Atlanta Georgia: DEStech Publishing, 2014: 173−178. [31] 张明丛. 周向MLEFP成型机理研究[D]. 南京: 南京理工大学, 2015ZHANG M C. Study on forming mechanism of circumferential MLEFP [D]. Nanjing: Nanjing University of Science and Technology, 2015. [32] 张明丛, 杜忠华, 周涛, 等. 药型罩结构参数对周向多线性爆炸成型弹丸成型及侵彻能力的影响 [J]. 火炸药学报, 2016, 39(1): 60–65,69. DOI: 10.14077/j.issn.1007-7812.2016.01.011.ZHANG M C, DU Z H, ZHOU T, et al. Influence of liner structural parameters on formation and penetration of the circumferential multiple linear explosive formation penetration [J]. Chinese Journal of Explosives and Propellants, 2016, 39(1): 60–65,69. DOI: 10.14077/j.issn.1007-7812.2016.01.011. [33] 王康康. 基于线性装药的大威力EFP研究[D]. 南京: 南京理工大学, 2015WANG K K. Research on EFP based on linear charge [D]. Nanjing: Nanjing University of Science and Technology, 2015. [34] 李鹏, 袁宝慧, 李刚, 等. 一种杆式周向多爆炸成型弹丸战斗部仿真及实验研究 [J]. 南京理工大学学报, 2017, 41(6): 681–685. DOI: 10.14177/j.cnki.32-1397n.2017.41.06.003.LI P, YUAN B H, LI G, et al. Simulation and experimental study on warhead of rod-shaped circumferential multiple explosively formed penetrator [J]. Journal of Nanjing University of Science and Technology, 2017, 41(6): 681–685. DOI: 10.14177/j.cnki.32-1397n.2017.41.06.003. [35] 李鹏, 袁宝慧. 爆炸成型杆式弹丸战斗部技术研究[C]// 第15届全国战斗部与毁伤技术学术交流会论文集. 北京, 2017, 11: 305-307LI P, YUAN B H. Research on the warhead numerical simulation of a rod type circumferential multiple explosively formed projectile [C]// Proceedings of the 15th National Symposium on Warhead and Damage Technology. Beijing, 2017, 11: 305−307. [36] 威廉·普·沃尔特斯, 乔纳斯·埃·朱卡斯. 成型装药原理及其应用[M]. 王树魁, 贝静芬, 译. 北京: 兵器工业出版社, 1992.WALTERS, ZUKAS. Fundamentals of shaped charges [M]. WANG S K, BEI J F, trans. Beijing: Weapons Industry Press, 1992. [37] 王树山. 终点效应学[M]. 北京: 科学出版社, 2019.WANG S S. Terminal effects [M]. Beijing: Science Press, 2019. [38] 张守中. 爆炸与冲击动力学[M]. 北京: 兵器工业出版社, 1993. [39] 苗润源, 李有仓, 梁争峰. 药型罩曲率半径对周向线性多爆炸成型弹丸成型性能及毁伤效应的影响 [J]. 科学技术与工程, 2017, 17(31): 235–238. DOI: 10.3969/j.issn.1671-1815.2017.31.038.MIAO R Y, LI Y C, LIANG Z F. Liner radius’ influence on circumferential linear multiple explosively formed projectiles' forming performance and damage effects [J]. Science Technology and Engineering, 2017, 17(31): 235–238. DOI: 10.3969/j.issn.1671-1815.2017.31.038. [40] 苗润源, 梁争峰, 程淑杰. 速度梯度对周向MLEFP的成型性能影响分析[C]// 第十六届战斗部与毁伤技术学术交流会论文集. 昆明, 2019: 276−278.MIAO R Y, LIANG Z F, CHENG S J. Analysis of the influence of velocity gradient on the forming performance of circumferential MLEFP [C]// Proceedings of the 16th Academic Conference on Warhead and Damage Technology. Kunming, 2019: 276−278. [41] 陈曦, 刘杰, 杜忠华, 等. 两端起爆下曲率药型罩线性装药的成型特性 [J]. 火炸药学报, 2018, 41(5): 441–446. DOI: 10.14077/j.issn.1007-7812.2018.05.003.CHEN X, LIU J, DU Z H, et al. Forming characteristics of linear charge with curved liner under two-end initiation [J]. Chinese Journal of Explosives and Propellants, 2018, 41(5): 441–446. DOI: 10.14077/j.issn.1007-7812.2018.05.003. [42] 李鹏, 袁宝慧, 李刚, 等. 一种杆式周向多爆炸成型弹丸战斗部数值模拟研究 [J]. 弹箭与制导学报, 2017, 37(1): 69–72. DOI: 10.15892/j.cnki.djzdxb.2017.01.016.LI P, YUAN B H, LI G, et al. Research on the warhead numerical simulation of a rod type circumferential multiple explosively formed projectile [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2017, 37(1): 69–72. DOI: 10.15892/j.cnki.djzdxb.2017.01.016. [43] 李鹏, 李刚, 袁宝慧, 等. 一种杆式多爆炸成型侵彻体战斗部 [J]. 爆炸与冲击, 2018, 38(4): 883–890. DOI: 10.11883/bzycj-2016-0356.LI P, LI G, YUAN B H, et al. A rod-shaped explosively formed penetrator warhead [J]. Explosion and Shock Waves, 2018, 38(4): 883–890. DOI: 10.11883/bzycj-2016-0356. [44] 李鹏, 袁宝慧, 丁刚, 等. LEFP战斗部毁伤元成型及初速分析研究 [J]. 科学技术与工程, 2017, 17(10): 188–191. DOI: 10.3969/j.issn.1671-1815.2017.10.032.LI P, YUAN B H, DING G, et al. The study of linear explosively formed penetrators warhead kill element on forming and initial velocity [J]. Science Technology and Engineering, 2017, 17(10): 188–191. DOI: 10.3969/j.issn.1671-1815.2017.10.032. [45] 史云鹏, 袁宝慧, 梁争峰, 等. 线形EFP药型罩设计 [J]. 火炸药学报, 2007, 30(3): 37–40,44. DOI: 10.3969/j.issn.1007-7812.2007.03.010.SHI Y P, YUAN B H, LIANG Z F, et al. The design of linear EFP liner [J]. Chinese Journal of Explosives & Propellants, 2007, 30(3): 37–40,44. DOI: 10.3969/j.issn.1007-7812.2007.03.010. [46] 李鹏, 李刚, 袁宝慧, 等. 旋转对爆炸成型杆式侵彻体毁伤威力的影响 [J]. 爆炸与冲击, 2018, 38(3): 616–621. DOI: 10.11883/bzycj-2016-0263.LI P, LI G, YUAN B H, et al. Influence of rotation on damage power of an explosively-formed rod-like penetrator [J]. Explosion and Shock Waves, 2018, 38(3): 616–621. DOI: 10.11883/bzycj-2016-0263. [47] 李鹏, 袁宝慧, 李刚, 等. 新型多线形爆炸成型弹丸战斗部的弹丸成型及毁伤研究 [J]. 火炸药学报, 2017, 40(1): 65–68;80. DOI: 10.14077/j.issn.1007-7812.2017.01.013.LI P, YUAN B H, LI G, et al. Study on the formation of penetrators and damage of a new type of multiple linear explosively formed penetrator warhead [J]. Chinese Journal of Explosives and Propellants, 2017, 40(1): 65–68;80. DOI: 10.14077/j.issn.1007-7812.2017.01.013. 期刊类型引用(3)
1. Zheng-qing Zhou，Ze-chen Du，Xiao Wang，Hui-ling Jiang，Qiang Zhou，Yu-long Zhang，Yu-zhe Liu，Pei-ze Zhang. Micro defects formation and dynamic response analysis of steel plate of quasi-cracking area subjected to explosive load. Defence Technology. 2024(01): 580-593 . 
必应学术2. 王路生，罗龙，刘浩，杨鑫，丁军，宋鹍，路世青，黄霞. 冲击速度对单晶镍层裂行为的影响规律及作用机制. 物理学报. 2024(16): 123-134 . 百度学术3. 杜欣，袁福平，熊启林，张波，阚前华，张旭. CoCrFeMnNi高熵合金冲击波响应与层裂强度的分子动力学研究. 力学学报. 2022(08): 2152-2160 . 百度学术其他类型引用(2)
-
下载:
必应学术
计量
- 文章访问数: 676
- HTML全文浏览量: 499
- PDF下载量: 150
- 被引次数: 5