串联战斗部不同介质组合的隔爆能力

刘宏杰; 王伟力; 苗润; 吴世永; 王俊华

doi:10.11858/gywlxb.20180585

串联战斗部不同介质组合的隔爆能力

doi: 10.11858/gywlxb.20180585

刘宏杰^1,,
王伟力^2,*, ,,
苗润¹,
吴世永³,
王俊华⁴

1.
海军航空大学岸防兵学院，山东烟台　264001
2.
海军工程大学兵器工程学院，湖北武汉　430000
3.
海军航空大学基础部，山东烟台　264001
4.
沈阳航空航天大学航空航天工程学部，辽宁沈阳　110136

详细信息

作者简介:
刘宏杰（1993－），男，硕士研究生，主要从事弹药设计与目标毁伤评估研究. E-mail: 1300202650@qq.com

通讯作者:
王伟力（1962－），男，教授，博士生导师，主要从事战斗部毁伤效应研究. E-mail: 858769430@qq.com

中图分类号: O382; TJ410.3
计量
- 文章访问数: 8209
- HTML全文浏览量: 3463
- PDF下载量: 30
出版历程
- 收稿日期: 2018-06-20
- 修回日期: 2018-08-12

Explosive Interruption of Tandem Warhead with Different Multilayer Structures

1.
Coastal Defense College, Naval Aeronautical University, Yantai 264001, China
2.
Academy of Weapons Engineering, Naval Engineering University, Wuhan 430000, China
3.
Basic Department, Naval Aeronautical University, Yantai 264001, China
4.
Faculty of Aerospace Engineering, Shenyang Aerospace University, Shenyang 110136, China

摘要

摘要: 为了有效提升串联切割战斗部隔爆结构衰减爆炸冲击波的性能，解决前级聚能装药结构与后级随进弹的匹配及隔爆问题，在前级切割器和后级随进弹之间加装隔爆结构，使用有限元分析软件ANSYS/LS-DYNA建立模型，进行不同组合结构隔爆性能的数值模拟，比较隔爆能力。模拟结果表明：前级装药起爆后，爆炸冲击波首先向后级随进弹头靠里区域汇聚，而不是向弹头尖端区域汇聚，因此可以适当减薄外层金属隔爆介质头部尖端区域；将外层金属由硬质钢改成铝时，后端壳体应力峰值的变化很小，故确定外层金属介质为铝；铝-聚脲的隔爆能力优于铝-泡沫铝结构，最终确定“软”隔爆介质为聚脲。通过调整铝和聚脲层的厚度，确定了最佳隔爆参数，能够满足实际应用。
- 串联战斗部 /
- 冲击波 /
- 多层介质 /
- 隔爆
Abstract: In this work we improved the explosion shock wave mitigation capability of the interruption structure in the tandem cutting warhead and matched the front shaped charge with the post-stage projectile and the explosive interruption by installing the explosive interruption structure between the front cutter and the rear following projectile, and established the simulation model using ANSYS/LS-DYNA, a finite element analysis software. We also analyzed the mitigation capabilities of the structures in different combinations and compared their explosive interruption capabilities using numerical simulation. The results showed that the explosion shock wave first converged to the inner region of the projectile, rather than to the tip of the warhead, so that the tip area of the outer metal flameproof medium was thinned properly. When the outer metal turned from hard steel to aluminum, little change was observed in the stress peak value of the rear end shell, suggesting that the outer metal medium was aluminum. Judging by the comparative study, the aluminum polyurea flameproof capability was superior to the aluminum foam-aluminum structure. Finally, the "soft" explosive interruption medium was determined to be polyurea, and the best explosion proof parameters were determined by adjusting the thickness of the aluminum and polyurea layer, which serves the needs of practical applications.
- tandem warhead /
- shock wave /
- multilayered media /
- explosive interruption

HTML全文

当前很多学者通过实验和数值仿真研究了聚脲衰减冲击波的特性^[1-3]。曾必强等^[4]通过数值模拟方法研究了不同初始速度的随进弹受前级装药爆炸冲击波影响产生的后抛、侧偏、翻转和力学响应规律，发现随进弹的后抛加速作用主要由前级装药爆炸产生的爆炸冲击波作用到弹体表面而产生，随进弹的侧偏现象则是作用在弹体的冲击波与冲击波后方的爆轰产物共同作用的结果。姜夕博等^[5]、徐森等^[6]研究了冲击波在有机玻璃中的衰减特性。侯海周等^[7]的实验结果表明，药柱冲击波经酚醛隔板传播时，其峰值压力随着隔板厚度的增加呈指数衰减。然而，单层材料的隔爆能力十分有限，需要通过增加隔爆层厚度提升衰减冲击波能力，但是却占用了战斗部空间，增加了战斗部质量。近年来，许多学者对多层组合介质中冲击波的衰减特性开展了大量研究^[8-10]。董永香等^[11]通过分析不同介质层中和层界面的爆炸波形和幅值得出，在脆性材料中放置泡沫材料可以有效减小爆炸波应力峰值，降低自由表面的拉伸应力波幅值，软夹层在多层介质中能够起到明显的削波和增加加载时间的作用，中间夹层改变了爆炸波能量的分配。陈闯等^[12]通过实验与数值模拟相结合的方法，将45钢、铝、有机玻璃进行不同组合，设计了顺序波阻抗梯度、逆序波阻抗梯度以及硬软硬3种多层隔爆结构，结果表明：逆序波阻抗梯度隔爆结构的输出冲击波压力最小，隔爆效果最佳。

为了在有限的战斗部空间内尽可能提高隔爆效果并减轻隔爆体质量，本研究在后级随进弹头部壳体表面加装隔爆体，采用ANSYS/LS-DYNA有限元软件进行数值模拟，研究不同结构的隔爆特性，以期使后级随进弹能顺利随进并侵彻后续目标靶。

1. 高压下冲击波在分层材料中的传播

高压下固体材料通常被认为是无黏性的可压缩流体，一般用介质Hugoniot关系曲线处理冲击波在两种不同介质上的反射和透射问题。定义介质的冲击阻抗为介质初始密度（固有）ρ₀和冲击波速度D_S的乘积，可近似用初始密度ρ₀和声速C₀的乘积表示。处理波从介质A传入介质B的最好方法是阻抗匹配法。碰撞界面的右行/左行冲击波满足连续方程和动量守恒。由界面连续条件可知：两介质界面上的粒子速度和压力相同。由一维冲击波理论可知，冲击波在两种界面处的反射波和透射波类型取决于两种介质的阻抗大小。当介质A的阻抗大于介质B时，反射波为稀疏波，反之反射波为冲击波^[13-14]。

当冲击波传播至两层介质分界面处，可以利用介质的Hugoniot关系计算透射波和反射波，如图1所示，曲线1～曲线3分别表示不同介质的Hugoniot关系曲线，介质的阻抗越大，由原点出发与Hugoniot关系曲线上某点连线的斜率越大。图1反映了压力幅值为P_b的冲击波从介质A传入介质B的情形，分界面处的状态既在反射波Hugoniot曲线1′（曲线1′和曲线1呈镜像对称关系）上，也在介质B的右行透射波Hugoniot曲线上，即a点状态。同理，冲击波从介质A传入介质C时，状态从b点变为c点。

图 1 反射波与透射波Hugoniot关系计算曲线

Figure 1. Hugoniot curves of reflection and transmission wave

下载: 全尺寸图片幻灯片

当固体处在高压条件下时，宋博等^[15]提出用p=Au²分析冲击波在多层介质中的传播，其中A可用来表示材料的“相对软硬”程度，A值越大，材料越“硬”，反之，材料越“软”。假设冲击波初始状态为b(p_b, u_b)，则曲线1′可以表示为

$p = {A_1}{(u - 2{u_b})^2}\;\;\;\;u < 2{u_b}$

(1)

则对于c点

$p = {A_1}{(u - 2{u_b})^2} = {A_2}{u^2}\;\;\;\;u < 2{u_b}$

(2)

b点状态为p_b=A₁u_b²，联立（2）式，消去u，求得透射系数T₂为

${T_{\rm{2}}} = \frac{p}{{{p_b}}} = \frac{4}{{{{\left(1 + \dfrac{{{A_1}}}{{{A_2}}}\right)}^2}}}$

(3)

由（3）式可以得出冲击波传播n层介质后的透射系数

${T'_n} = \frac{{{2^2}}}{{{{\left(1 + \dfrac{{{A_1}}}{{{A_2}}}\right)}^2}}} \times \dfrac{{{2^2}}}{{{{\left(1 + \dfrac{{{A_2}}}{{{A_3}}}\right)}^2}}} \times \cdots \times \dfrac{{{2^2}}}{{{{\left(1 + \dfrac{{{A_{n - 1}}}}{{{A_{n}}}}\right)}^2}}} = \prod\limits_{i = 1}^{n - 1} {\frac{{{2^2}}}{{{{\left(1 + \dfrac{{{A_i}}}{{{A_{i + 1}}}}\right)}^2}}}} = t_n^2$

(4)

其中

${t_n} = \prod\limits_{i = 1}^{n - 1} {\frac{2}{{1 + \dfrac{{{A_i}}}{{{A_{i + 1}}}}}}} = \dfrac{{{2^{n - 1}}}}{{\displaystyle\prod\limits_{i = 1}^{n - 1} {\left(1 + \dfrac{{{A_i}}}{{{A_{i + 1}}}}\right)} }}$

(5)

由此可知对于相同层数（n）的相同介质来说，如果有某种顺序的排列结构，使得

${S_n} = \prod\limits_{i = 1}^{n - 1} {\left(1 + \frac{{{A_i}}}{{{A_{i + 1}}}}\right)}$

(6)

存在最大值，那么这种结构衰减冲击波的性能最好，即这种结构能最大限度地削弱透射波强度。

2. 多层介质组合衰减冲击波特性分析及数值模拟

第1节的分析仅考虑了不同介质的排列次序对衰减冲击波的影响，而未考虑每层介质材料本身对波的衰减及组合介质不同结构对衰减冲击波的影响。下面通过数值模拟方法首先确定最佳组合介质，以期最大限度地削弱透射波，同时兼顾质量和体积。

2.1 模型建立

为了便于建立计算模型，作以下假设：（1）空气、炸药、药型罩和装药壳体均为连续介质；（2）整个爆炸过程为绝热过程。串联战斗部及隔爆结构如图2所示。为了保证后级随进弹能顺利侵彻，隔爆体A选择金属介质，考虑到计算模型和爆炸载荷具有对称特性，为提高计算效率，建立1/4模型，施加对称边界条件。计算模型由炸药、空气、药型罩、壳体、后级随进弹和隔爆体组成，其中炸药、空气和药型罩定义为Euler网格，壳体、后级随进弹和隔爆体定义为Lagrange网格，隔爆体之间设置面面自动接触。选用多物质ALE（Arbitrary Lagrange-Euler）方法模拟分析前级爆轰场对后级随进弹的影响。起爆方式采用顶端环形起爆。

图 2 串联战斗部及隔爆结构有限元模型

Figure 2. Finite element model of tandem warhead and explosive interruption structure

下载: 全尺寸图片幻灯片

2.2 材料模型及状态方程

炸药选用B炸药，采用MAT_HIGH-EXPLOSIVE-BURN模型和JWL状态方程^[16-17]描述

$p = A\left(1 - \frac{\omega }{{{R_1}\nu }}\right){{\rm{e}}^{ - {R_1}v}} + B\left(1 - \frac{\omega }{{{R_2}\nu }}\right){{\rm{e}}^{ - {R_2}v}} + \frac{{\omega E}}{\nu }$

(7)

式中：p为爆轰压力；E为炸药比内能；v为相对比容；A、B、R₁、R₂、ω为参数，取值列于表1。表1中：ρ为密度，D为B炸药爆速，p_CJ为CJ面压力，V₀为初始相对比容。

表 1 B炸药材料参数

Table 1. Material parameters of composition B

ρ/(g·cm^-3)	D/(m·s^-1)	p_CJ/GPa	A/GPa	B/GPa	R₁	R₂	ω	E₀/GPa	V₀
1.713	7 500	28.6	524.2	7.678	4.2	1.1	0.34	8.499	1.0

下载: 导出CSV

| 显示表格

药型罩材料选用紫铜，采用MAT_Steinberg材料模型和Grüneisen状态方程^[18]描述

$p = \dfrac{{{\rho _0}{C^2}\mu \left[ {1 + \left(1 - \dfrac{{{\gamma _0}}}{2}\right)\mu - \dfrac{a}{2}{\mu ^2}} \right]}}{{{{\left[ {1 - ({S_1} - 1)\mu - {S_2}\dfrac{{{\mu ^2}}}{{\mu + 1}} - {S_3}\dfrac{{{\mu ^3}}}{{{{(\mu + 1)}^2}}}} \right]}^2}}} + ({\gamma _0} + a\mu )E$

(8)

式中：E为单位体积内能；ρ₀为材料初始密度；C为u_s与u_p的交点，其中u_s为冲击波速度，u_p为粒子速度；S₁、S₂和S₃是u_s-u_p曲线的斜率；γ₀为Grüneisen系数；a为对γ₀的修正；μ=(1/V)-1，V为当前空气的相对体积。

壳体选用30CrMnSiNi2A钢，用MAT_Johnson-Cook材料模型和Grüneisen状态方程描述；铝选用MAT_Johnson-Cook材料模型和Grüneisen状态方程描述，参数见表2，其中A_JC、B_JC、n、m为Johnson-Cook材料模型参数，T_m为材料熔化温度，T₀为常压下的熔化温度，a为对Grüneisen系数γ₀的修正。空气采用MAT_Null材料模型和Linear-Polynomial状态方程描述。

表 2 30CrMnSiNi2A钢和金属铝计算参数

Table 2. Material performance parameters of 30CrMnSiNi2A steel and aluminum

Material	A_JC/MPa	B_JC/MPa	n	C	m	T_m/K	T₀/K	S₁	γ₀	a
Steel	1 280	420	0.30	0.030	1.00	1 793	294	1.490	2.17	0.46
Al	265	426	0.34	0.015	1.00	775	294	1.345	2.13	0.10

下载: 导出CSV

| 显示表格

聚脲选用MAT_PIECEWISE_LINEAR_PLASTICITY材料模型^[19]描述。该模型常用来模拟弹塑性材料，屈服准则由断裂应变定义，应力-应变曲线由屈服强度、剪切模量定义或者通过8组塑性应变和应力的数组定义，应变率的影响采用Cowper-Symbols模型，弹性体材料各参数为^[20-21]: ρ=1 020 kg/m³，杨氏模量E=230 MPa，泊松比ν=0.4，初始屈服强度σ₀=14 MPa。

3. 数值模拟结果与分析

用B1、B2、B3和B4分别表示隔爆体A/B、B/C、C/D、D/后级弹壳体的界面，计算中两个界面处均采用面面自动接触条件。表3列出了5种工况。

表 3 不同工况及简便书写方式

Table 3. Different conditions and indications

Condition	Material
Condition	A	B	C	D
Ⅰ	Steel	Polyurea	Aluminum foam	Aluminum
Ⅱ	Aluminum	Polyurea	Aluminum foam	Aluminum
Ⅲ	Aluminum	Aluminum foam	Aluminum foam	Aluminum
Ⅳ	Aluminum	Aluminum foam	Aluminum foam	Aluminum foam
Ⅴ	Aluminum	Polyurea	Polyurea	Polyurea

下载: 导出CSV

| 显示表格

3.1 爆轰波在介质中的衰减

工程实践表明：具有软夹层的硬软硬多层介质具有较明显的抗冲击和抗爆炸能力，其中硬层抗击近距离的爆炸波作用，发挥其刚度与强度效应，而软层由多孔泡沫材料组成，波阻抗低，强度和刚度低，吸波能力强。如果外加爆炸波的持续时间足够长，软夹层不再发挥作用，软夹层下面的硬层将受到较大的冲击作用。

这里假定几种材料的Hugoniot曲线分别为p=A₁u(硬质钢)，p=A₂u(铝)，p=A₃u(泡沫铝)，p=A₄u (聚脲)，并且有A₁>A₂>A₃>A₄。

3.2 硬软硬组合介质的隔爆能力分析

（1）工况I

对“硬质钢-聚脲-泡沫铝-铝”结构，有

${S_{n1}} = \left(1 + \frac{{{A_1}}}{{{A_4}}}\right)\left(1 + \frac{{{A_3}}}{{{A_4}}}\right)\left(1 + \frac{{{A_2}}}{{{A_4}}}\right)$

(9)

装药起爆后，冲击波首先撞击复合隔爆结构中的硬质钢，该碰撞过程为非弹性碰撞。在该冲击速度下钢板发生弹塑性变形，历时很短，通常钢板发生局部的弹塑性变形，而其他区域受碰撞挤压作用影响较小。当冲击波传播至B1界面时，发生反射和透射，聚脲作为一种弹性涂覆体，其冲击阻抗小于硬质钢，界面B1处反射稀疏波，然而聚脲层的厚度仅为1 mm，衰减冲击波的能力有限。冲击波从介质B传入介质C时，在界面B2处反射冲击波，接下来讨论中间泡沫铝层的塑性变形和卸载冲击波强度的变化规律。泡沫铝在冲击载荷作用下产生递增硬化的变形即破坏波（压实波），其非线性应力-应变曲线可简化成3个阶段：弹性段、屈服塑性段、压实段，其力学性能除了与泡沫的泡孔结构有关外，还与基体材料性质、应变率和密度等有关，而且对静水压有很强的依赖性。应力波在泡沫铝中经历以下两个步骤：一是泡沫铝压实过程中的冲击波传播；二是泡沫铝压实后，应力波在包括外层硬质钢、聚脲涂覆层、压实泡沫铝以及后端铝介质的多次反射和透射。冲击波从介质C传入介质D时，在界面B3处反射冲击波，经过前面3层介质的衰减作用后，铝介质的冲击压力较低，铝介质作为一种软金属材料，在冲击载荷作用下易发生变形和破坏，因而能最大限度地吸收能量，最大程度地减小后级装药的压力峰值。

为了检验工况I衰减冲击波的能力，查看随进弹壳体最大应力集中的区域以及装药压力峰值的区域和时刻。t=264 μs时，装药内部靠近轴线区域压力（见图3(a)）达到最大值，为143.6 MPa；t=189 μs时，外侧壳体应力集中区域的应力峰值如图3(b)所示，可见该区域未发生塑性变形；应力峰值附近某单元应力-时间变化关系如图3(c)所示。

图 3 工况I中装药压力及壳体应力分析

Figure 3. Loading pressure and shell stress of working condition I

下载: 全尺寸图片幻灯片

由图4可知：t=30 μs时，外层钢介质头部尖端靠里区域首先发生弹性变形；随着冲击波在介质中继续传播，t=38 μs时，钢介质弹性变形区域向后级随进弹方向移动；t=45 μs时，应力集中区域为钢介质最薄处；t=51 μs时，钢介质最薄处应力峰值继续增加；t=60 μs时，应力峰值达到最大值，为1 130 MPa。前级装药采用中空结构，装药起爆后爆轰波向后运动并汇聚的过程中，首先作用到钢介质尖端靠里区域，考虑到节约质量成本，钢介质头部尖端区域可适当减薄。

图 4 典型时刻外层钢介质的应力云图

Figure 4. Stress nephogram of steel at typical times

下载: 全尺寸图片幻灯片

（2）工况II

与工况I相比，为了控制质量，外侧介质A由密度较大的硬质钢变为密度较小的铝。透射系数影响因子为

${S_{n{\rm{2}}}} = \left(1 + \frac{{{A_2}}}{{{A_4}}}\right)\left(1 + \frac{{{A_3}}}{{{A_4}}}\right)\left(1 + \frac{{{A_2}}}{{{A_4}}}\right)$

(10)

比较S_n₁和S_n₂的大小，得到S_n₁>S_n₂，可知工况I所示结构中衰减冲击波压力的性能优于工况II所示结构。

装药起爆后，冲击波首先撞击复合隔爆结构中的金属铝，该碰撞过程为非弹性碰撞。金属铝的屈服强度小，在高速冲击载荷下铝板发生破坏，部分区域发生弹塑性变形。冲击波在聚脲层、泡沫铝和后端铝介质传播衰减分析方法与工况I类似。

由图5可知：t=120 μs时，装药尖端压力达到最大值，为421.2 MPa；相同时刻，外侧壳体应力出现最大值1 062.0 MPa，见图5(b)；应力峰值附近某单元的应力-时间曲线如图5(c)所示，后级随进弹壳体未进入塑性变形段。硬软硬结构存在不利的一面：随着作用时间的增加，对于多层介质而言，在波动效应后期结构振动效应增强。

图 5 工况II装药压力及壳体应力分析

Figure 5. Loading pressure and shell stress of working condition II

下载: 全尺寸图片幻灯片

（3）工况III

工况III所示隔爆结构为铝-泡沫铝-铝硬软硬结构，与工况II相比，将介质B换为泡沫铝。含有低阻抗软夹层的多层介质能够减小应力幅值，增加波形脉宽（即增加作用时间），从而减小由于应力波峰超过材料强度而引起局部破坏的可能性。

由图6可知：t=117 μs时，装药尖端压力达到最大值363.5 MPa，与中间有聚脲层时对应的压力峰值时刻基本一致；t=174 μs时，外侧壳体应力峰值为1 062.0 MPa，峰值区域如图6(b)所示；应力峰值附近单元的应力-时间曲线如图6(c)所示，后级随进弹壳体未进入塑性变形段。

图 6 工况III装药压力及壳体应力

Figure 6. Loading pressure and shell stress of working condition III

下载: 全尺寸图片幻灯片

3.3 硬软两层组合介质的隔爆能力分析

（1）工况IV

隔爆介质为两层，外层为金属铝，里层为泡沫铝，透射系数影响因子为

${S_{n{\rm{4}}}} = 1 + \frac{{{A_2}}}{{{A_3}}}$

(11)

对比工况IV（见图7）和工况III：将厚端铝介质改为泡沫铝后，增加了泡沫铝的厚度，在爆炸冲击波作用下，泡沫铝内衬层能够吸收更多的能量，随进弹装药压力峰值比硬软硬结构中装药应力峰值出现得晚些；冲击波继续传播过程中，泡沫铝衰减冲击波的能力不再增加，最终泡沫铝被压实，密度增加，弹性模量也增加，紧挨软夹层的金属层受到较大的冲击作用，继续衰减冲击波。

图 7 工况IV装药压力及壳体应力

Figure 7. Loading pressure and shell stress of working condition IV

下载: 全尺寸图片幻灯片

（2）工况V

隔爆介质分层排列顺序为金属铝和弹性涂覆体聚脲，其透射系数影响因子为

${S_{n{\rm{5}}}} = 1 + \frac{{{A_2}}}{{{A_4}}}$

(12)

计算可得S_n₅>S_n₄。

考虑将隔爆结构简化为金属铝-聚脲结构（即工况V），对比工况IV（即将泡沫铝改为弹性介质聚脲），则弹性体厚度明显增加。工况V中装药压力及壳体应力见图8：t=156 μs时，装药尖端压力达到最大值272.0 MPa，与工况IV中出现压力峰值时刻一致；t=192 μs时，外侧壳体应力峰值为1 050.0 MPa；应力峰值附近单元的应力-时间曲线显示，后级随进弹壳体未进入塑性变形段。聚脲作为一种喷涂材料，受制作工艺限制，喷涂有效厚度有限，该工况的聚脲层厚度为4 cm，远大于聚脲层的最大有效厚度。

图 8 工况V装药压力及壳体应力

Figure 8. Loading pressure and shell stress of working condition V

下载: 全尺寸图片幻灯片

由图9可知：采用工况II和工况III所示隔爆结构时，后级随进弹出现加速度峰值的时间最早，且峰值最大；工况I所示隔爆结构中后级弹的加速度峰值最小，工况IV和工况I出现加速度峰值的时间基本一致；工况V所示结构延缓了冲击波作用到后级随进弹的时间，出现加速度峰值的时间最晚。由图10可知，采用上述5种隔爆结构时，后级随进弹速度增加到一定值时不再增加。

图 9 后级随进弹加速度-时间变化曲线

Figure 9. Acceleration vs. time for the following stage’s incoming projectile

下载: 全尺寸图片幻灯片

图 10 后级随进弹速度-时间变化曲线

Figure 10. Velocity vs. time for the following stage’s incoming projectile

下载: 全尺寸图片幻灯片

3.4 总结分析

采用工况I的隔爆结构时，后端装药的压力峰值出现得最晚，且压力峰值最小，仅为143.6 MPa；受质量成本的限制，工况II中介质A由密度较大的硬质钢换成密度较小的金属铝，金属铝的屈服强度较低，爆炸冲击载荷下容易发生破坏，部分区域发生弹塑性变形，装药压力峰值区域为装药尖端区域，远早于A介质为硬质钢的情形，且峰值压力为421.2 MPa，远大于工况I中对应的峰值压力，但后级壳体应力峰值与工况I的应力峰值相近，均未进入塑性变形区；工况III与工况II相比，介质B由弹性涂覆体聚脲改为泡沫铝，即隔爆结构转变为“铝-泡沫铝-铝”硬软硬结构，由于介质B的厚度仅为1 mm，对冲击波反射和透射的影响不大，因此硬软硬结构出现装药压力峰值的时刻与工况II基本一致，装药峰值压力小于工况II的情形；工况IV和工况V都降低了隔爆层数，分别为铝-泡沫铝结构和铝-聚脲结构，对比透射系数影响因子，前者小于后者，因此得出铝-聚脲结构隔爆效果优于铝-泡沫铝结构，但工况V中聚脲层的厚度过大，无法实现实际应用。接下来将改进隔爆结构，使之既能充分提升隔爆性能，节约体积和质量成本，又能实现实际应用。

4. 改进结构的数值模拟结果与讨论

由第3节分析可知，装药起爆后，冲击波首先作用到最外层隔爆体头部靠里区域。为了控制总质量，适当减薄金属层尖端区域的厚度，改进后隔爆结构如图11所示，保持后级弹尖端与前级聚能装药结构后端的距离不变，即图中d。表4列出了4种改进结构。

图 11 串联战斗部及改进隔爆结构有限元模型

Figure 11. Finite element model of tandem warhead and improved explosive interruption structure

下载: 全尺寸图片幻灯片

表 4 改进结构

Table 4. Improved structures

Condition	Material
Condition	A	B	C
Ⅰ′	Steel	Aluminum	Polyurea
Ⅱ′	Steel	Aluminum	Aluminum foam
Ⅲ′	Aluminum	Aluminum	Polyurea
Ⅳ′	Aluminum	Aluminum	Polyurea（2 cm）

下载: 导出CSV

| 显示表格

4.1 逆序波阻抗结构的隔爆能力分析

（1）工况I′

工况I′下透射系数影响因子为

${S'_{n1}} = \left(1 + \frac{{{A_1}}}{{{A_2}}}\right)\left(1 + \frac{{{A_2}}}{{{A_4}}}\right)$

(13)

由图12可知：t=163 μs时，装药尖端压力达到最大值114.8 MPa；t=243 μs时，外侧壳体应力峰值为1 188.0 MPa；应力峰值附近单元的应力-时间曲线显示，后级随进弹壳体尖端区域进入塑性变形区。

图 12 工况I′装药压力及壳体应力

Figure 12. Loading pressure and shell stress of working condition I′

下载: 全尺寸图片幻灯片

（2）工况II′

工况II′下，透射系数影响因子为

${S'_{n2}} = \left(1 + \frac{{{A_1}}}{{{A_2}}}\right)\left(1 + \frac{{{A_2}}}{{{A_3}}}\right)$

(14)

比较可得，工况I′对应的透射系数大于工况II′，说明工况I′所示隔爆结构的衰减冲击波能力大于工况II′。由图13可知：t=246 μs时，装药尖端压力达到最大值125.7 MPa；t=297 μs时，外侧壳体应力峰值为1 209.0 MPa；应力峰值附近单元的应力-时间曲线显示，后级随进弹壳体尖端区域进入塑性变形区。

图 13 工况II′装药压力及壳体应力分析

Figure 13. Loading pressure and shell stress of working condition II′

下载: 全尺寸图片幻灯片

4.2 改进结构硬软两层组合介质的隔爆能力分析

（1）工况III′

将最外层密度较大的硬质钢改为密度较小的金属铝，采用两层隔爆层，即金属铝-聚脲结构。由图14可知：t=159 μs时，装药内部尖端区域压力达到最大值170.8 MPa；t=195 μs时，外侧壳体应力峰值为1 005.0 MPa；应力峰值附近单元的应力-时间曲线显示，后级随进弹壳体未进入塑性变形区。

图 14 工况III′装药压力及壳体应力分析

Figure 14. Loading pressure and shell stress of working condition III′

下载: 全尺寸图片幻灯片

（2）工况IV′

考虑到聚脲喷涂的厚度有限，将3 cm聚脲层减薄至2 cm，研究其隔爆性能。由图15可知：t=144 μs时，装药内部尖端区域压力达到最大值296.3 MPa；t=177 μs时，外侧壳体应力峰值为1 035.0 MPa；应力峰值附近单元的应力-时间曲线显示，后级随进弹壳体未进入塑性变形区。

图 15 工况IV′装药压力及壳体应力分析

Figure 15. Loading pressure and shell stress of working condition IV′

下载: 全尺寸图片幻灯片

由图16和图17可知：采用工况IV′所示隔爆结构时，后级随进弹出现加速度峰值的时间最早，且峰值最大；工况I′所示隔爆结构中，后级随进弹的加速度峰值最小；工况IV′和工况III′出现加速度峰值的时间基本一致；工况II′所示结构延缓了冲击波作用到后级随进弹的时间，出现加速度峰值的时间最晚。由图17可知，采用以上所述4种隔爆结构时，后级随进弹速度增加到一定值时将不再增加。

图 16 后级随进弹加速度-时间变化曲线

Figure 16. Acceleration vs. time for the following stage’s incoming projectile

下载: 全尺寸图片幻灯片

图 17 后级随进弹加速度-时间曲线

Figure 17. Velocity vs. time for the following stage’s incoming projectile

下载: 全尺寸图片幻灯片

4.3 改进结构的隔爆能力总结

总结改进后的工况：采用工况I′对应的隔爆结构时，后端装药压力峰值最小，仅为114.8 MPa；受质量成本的限制，工况II′中将介质C换成泡沫铝，延缓了冲击波作用到达后级随进弹的时间。但是在工况I′和工况II′对应的隔爆结构中，由于外层硬质钢挤压变形，使得后级随进战斗部壳体头部发生变形，进入塑性变形区；工况IV′与工况III′相比，介质C对应弹性涂覆体聚脲的厚度减薄，装药压力较工况III′所示结构有所提高，但增幅不大，后端壳体应力峰值的变化不大。

5. 结　论

采用数值模拟方法研究不同介质组合衰减冲击波的能力，通过对比后级弹壳体应力峰值及装药压力峰值，考察不同介质组合的隔爆能力，得出以下结论：（1）前级中空结构装药起爆后，爆炸冲击波首先向后级随进弹头靠里区域汇聚，而不是向弹头尖端区域汇聚；（2）在“硬软硬”隔爆结构中，将外层介质A由密度较大的硬质钢改为金属铝后，后级弹壳体应力峰值的变化很小，故确定隔爆结构中外层金属介质为铝，从而大大节约质量成本；（3）“硬软”两层组合介质隔爆结构中，铝-聚脲的隔爆能力优于铝-泡沫铝结构，最终确定“软”隔爆介质为聚脲，但聚脲的喷涂厚度有限，需要进一步减薄聚脲层的厚度；（4）改进结构中，减薄了外层金属介质层尖端的厚度，若保持隔爆层厚度不变，通过增加金属铝的厚度降低聚脲层的厚度，最终确定聚脲层的厚度为2 cm，从而既节约质量和体积成本，又具有优良的隔爆性能，能够充分削弱作用在后级随进弹上的冲击波强度，使后级随进弹的加速度和速度较小。

图 1 反射波与透射波Hugoniot关系计算曲线

Figure 1. Hugoniot curves of reflection and transmission wave

下载: 全尺寸图片幻灯片

图 2 串联战斗部及隔爆结构有限元模型

Figure 2. Finite element model of tandem warhead and explosive interruption structure

下载: 全尺寸图片幻灯片

图 3 工况I中装药压力及壳体应力分析

Figure 3. Loading pressure and shell stress of working condition I

下载: 全尺寸图片幻灯片

图 4 典型时刻外层钢介质的应力云图

Figure 4. Stress nephogram of steel at typical times

下载: 全尺寸图片幻灯片

图 5 工况II装药压力及壳体应力分析

Figure 5. Loading pressure and shell stress of working condition II

下载: 全尺寸图片幻灯片

图 6 工况III装药压力及壳体应力

Figure 6. Loading pressure and shell stress of working condition III

下载: 全尺寸图片幻灯片

图 7 工况IV装药压力及壳体应力

Figure 7. Loading pressure and shell stress of working condition IV

下载: 全尺寸图片幻灯片

图 8 工况V装药压力及壳体应力

Figure 8. Loading pressure and shell stress of working condition V

下载: 全尺寸图片幻灯片

图 9 后级随进弹加速度-时间变化曲线

Figure 9. Acceleration vs. time for the following stage’s incoming projectile

下载: 全尺寸图片幻灯片

图 10 后级随进弹速度-时间变化曲线

Figure 10. Velocity vs. time for the following stage’s incoming projectile

下载: 全尺寸图片幻灯片

图 11 串联战斗部及改进隔爆结构有限元模型

Figure 11. Finite element model of tandem warhead and improved explosive interruption structure

下载: 全尺寸图片幻灯片

图 12 工况I′装药压力及壳体应力

Figure 12. Loading pressure and shell stress of working condition I′

下载: 全尺寸图片幻灯片

图 13 工况II′装药压力及壳体应力分析

Figure 13. Loading pressure and shell stress of working condition II′

下载: 全尺寸图片幻灯片

图 14 工况III′装药压力及壳体应力分析

Figure 14. Loading pressure and shell stress of working condition III′

下载: 全尺寸图片幻灯片

图 15 工况IV′装药压力及壳体应力分析

Figure 15. Loading pressure and shell stress of working condition IV′

下载: 全尺寸图片幻灯片

图 16 后级随进弹加速度-时间变化曲线

Figure 16. Acceleration vs. time for the following stage’s incoming projectile

下载: 全尺寸图片幻灯片

图 17 后级随进弹加速度-时间曲线

Figure 17. Velocity vs. time for the following stage’s incoming projectile

下载: 全尺寸图片幻灯片

表 1 B炸药材料参数

Table 1. Material parameters of composition B

ρ/(g·cm^-3)	D/(m·s^-1)	p_CJ/GPa	A/GPa	B/GPa	R₁	R₂	ω	E₀/GPa	V₀
1.713	7 500	28.6	524.2	7.678	4.2	1.1	0.34	8.499	1.0

下载: 导出CSV

表 2 30CrMnSiNi2A钢和金属铝计算参数

Table 2. Material performance parameters of 30CrMnSiNi2A steel and aluminum

Material	A_JC/MPa	B_JC/MPa	n	C	m	T_m/K	T₀/K	S₁	γ₀	a
Steel	1 280	420	0.30	0.030	1.00	1 793	294	1.490	2.17	0.46
Al	265	426	0.34	0.015	1.00	775	294	1.345	2.13	0.10

下载: 导出CSV

表 3 不同工况及简便书写方式

Table 3. Different conditions and indications

Condition	Material
Condition	A	B	C	D
Ⅰ	Steel	Polyurea	Aluminum foam	Aluminum
Ⅱ	Aluminum	Polyurea	Aluminum foam	Aluminum
Ⅲ	Aluminum	Aluminum foam	Aluminum foam	Aluminum
Ⅳ	Aluminum	Aluminum foam	Aluminum foam	Aluminum foam
Ⅴ	Aluminum	Polyurea	Polyurea	Polyurea

下载: 导出CSV

表 4 改进结构

Table 4. Improved structures

Condition	Material
Condition	A	B	C
Ⅰ′	Steel	Aluminum	Polyurea
Ⅱ′	Steel	Aluminum	Aluminum foam
Ⅲ′	Aluminum	Aluminum	Polyurea
Ⅳ′	Aluminum	Aluminum	Polyurea（2 cm）

下载: 导出CSV

参考文献(21)

[1]	LEBLANC J, SHILLINGS C, GAUCH E, et al. Near field underwater explosion response of polyurea coated composite plates [J]. Experimental Mechanics, 2016, 56(4): 569–581. doi: 10.1007/s11340-015-0071-8
[2]	HARIS A, LEE H P, TAN V B C. An experimental study on shock wave mitigation capability of polyurea and shear thickening fluid based suspension pads [J]. Defence Technology, 2018, 14(1): 12–18. doi: 10.1016/j.dt.2017.08.004
[3]	DAI L, WU C, AN F, et al. Experimental investigation of polyurea-coated steel plates at underwater explosive loading [J]. Advances in Materials Science and Engineering, 2018: 1264276.
[4]	曾必强, 姜春兰, 严翰新, 等. 串联攻坚战斗部前级爆轰场对随进弹随进影响分析 [J]. 兵工学报, 2010(Suppl 1): 162–166 ZENG B Q, JIANG C L, YAN H X, et al. Analysis for effects of precursory detonation field on projectile following course in tandem warhead [J]. Acta Armamentarii, 2010(Suppl 1): 162–166
[5]	姜夕博, 饶国宁, 徐森, 等. 冲击波在有机玻璃中衰减特性的数值模拟与实验研究 [J]. 南京理工大学学报(自然科学版), 2012, 36(6): 1059–1064 doi: 10.3969/j.issn.1005-9830.2012.06.028 JIANG X B, RAO G N, XU S, et al. Numerical simulation and experimental research on shock wave attenuation properties in PMMA [J]. Journal of Nanjing University of Science and Technology, 2012, 36(6): 1059–1064 doi: 10.3969/j.issn.1005-9830.2012.06.028
[6]	徐森, 刘大斌, 彭金华, 等. 药柱冲击波在有机玻璃中的衰减特性研究 [J]. 高压物理学报, 2010, 24(6): 431–437 doi: 10.11858/gywlxb.2010.06.005 XU S, LIU D B, PENG J H, et al. Study on the shock wave attenuation of the booster charge in the PMMA gap [J]. Chinese Journal of High Pressure Physics, 2010, 24(6): 431–437 doi: 10.11858/gywlxb.2010.06.005
[7]	侯海周, 彭金华, 胡毅亭. 爆炸冲击波在酚醛层压材料中衰减特性的实验研究 [J]. 火工品, 2016(2): 13–16 doi: 10.3969/j.issn.1003-1480.2016.02.005 HOU H Z, PENG J H, HU Y T. Experimental study of shock wave attenuation properties in phenolic cotton fabric material [J]. Initiators & Pyrotechnics, 2016(2): 13–16 doi: 10.3969/j.issn.1003-1480.2016.02.005
[8]	TEDESCO J W, LANDIS D W. Wave propagation through layered systems [J]. Computers & Structures, 1989, 32(3/4): 625–638.
[9]	PETEL O E, JETTÉ F X, GOROSHIN S, et al. Blast wave attenuation through a composite of varying layer distribution [J]. Shock Waves, 2011, 21(3): 215–224. doi: 10.1007/s00193-010-0295-6
[10]	石少卿, 刘仁辉, 汪敏. 钢板-泡沫铝-钢板新型复合结构降低爆炸冲击波性能研究 [J]. 振动与冲击, 2008, 27(4): 143–146 doi: 10.3969/j.issn.1000-3835.2008.04.037 SHI S Q, LIU R H, WANG M. Shock wave reduction behavior of a new compound structure composed of a foam aluminum layer between two steel plates [J]. Journal of Vibration and Shock, 2008, 27(4): 143–146 doi: 10.3969/j.issn.1000-3835.2008.04.037
[11]	董永香, 冯顺山, 李学林. 爆炸波在硬-软-硬三明治介质中传播特性的数值分析 [J]. 弹道学报, 2007, 19(1): 59–63 doi: 10.3969/j.issn.1004-499X.2007.01.017 DONG Y X, FENG S S, LI X L. Numerical analysis of propagation characteristics of explosive wave in the hard-soft-hard sandwich media [J]. Journal of Ballistics, 2007, 19(1): 59–63 doi: 10.3969/j.issn.1004-499X.2007.01.017
[12]	陈闯, 王晓鸣, 李文彬, 等. 多层介质阻抗匹配对隔爆效果的影响 [J]. 振动与冲击, 2014, 33(17): 105–110 CHEN C, WANG X M, LI W B, et al. Influence of multilayered media impedance matching on explosion interruption effect [J]. Journal of Vibration and Shock, 2014, 33(17): 105–110
[13]	北京工业学院八系. 爆炸及其作用 [M]. 北京: 国防工业出版社, 1994: 105-112.
[14]	孙承纬. 爆炸物理学 [M]. 北京: 科学出版社, 2011: 205-207.
[15]	宋博, 胡时胜, 王礼立. 分层材料的不同排列次序对透射冲击波强度的影响 [J]. 兵工学报, 2000, 21(3): 272–274 doi: 10.3321/j.issn:1000-1093.2000.03.021 SONG B, HU S S, WANG L L. Influence on the transmitted intensity of shock wave through different tactic orders of layered materials [J]. Acta Armamentarii, 2000, 21(3): 272–274 doi: 10.3321/j.issn:1000-1093.2000.03.021
[16]	汪文革, 杨世军, 韩永要, 等. 基于ANSYS/LS-DYNA的聚能射流侵彻装甲钢的有限元分析 [J]. 兵工自动化, 2008, 27(3): 39–41 doi: 10.3969/j.issn.1006-1576.2008.03.015 WANG W G, YANG S J, HAN Y Y, et al. Finite element analysis of shaped charge jet penetrating into target based on ANSYS/LS-DYNA [J]. Ordnance Industry Automation, 2008, 27(3): 39–41 doi: 10.3969/j.issn.1006-1576.2008.03.015
[17]	TOGAMI T C, BAKER W E, FORRESTAL M J. A split Hopkinson bar technique to evaluate the performance of accelerometers [J]. Journal of Applied Mechanics, 1996, 63(2): 353–356. doi: 10.1115/1.2788872
[18]	FRANZEN R R, SCHNEIDEWIND P N. Observations concerning the penetration mechanics of tubular hypervelocity penetrators [J]. International Journal of Impact Engineering, 1991, 11(3): 289–303. doi: 10.1016/0734-743X(91)90040-M
[19]	WASMUND T L. New model to evaluate weapon effects and platform vulnerability: AJEM [J]. Wstiac Newsletter, 2001, 2: 1–3.
[20]	LEE W H, PAINTER J W. Material void-opening computation using particle method [J]. International Journal of Impact Engineering, 1999, 22(1): 1–22. doi: 10.1016/S0734-743X(98)00041-4
[21]	甘云丹, 宋力, 杨黎明. 弹性体涂覆钢板抗冲击性能的数值模拟 [J]. 兵工学报, 2009(Suppl 2): 15–18 GAN Y D, SONG L, YANG L M. Numerical simulation for anti-blast performances of steel plate coated with elastomer [J]. Acta Armamentarii, 2009(Suppl 2): 15–18