反应破片对带壳装药的冲击效应

赵宇峰; 段继; 冯顺山

doi:10.11883/bzycj-2024-0063

反应破片对带壳装药的冲击效应

doi: 10.11883/bzycj-2024-0063

赵宇峰^1,,
段继²,
冯顺山^3, ,

1.
中北大学航空宇航学院，山西太原 030051
2.
中北大学机电工程学院，山西太原 030051
3.
北京理工大学爆炸科学与技术全国重点实验室，北京 100081

基金项目: 山西省基础研究计划（202203021222074，20210302124340）

详细信息

作者简介:
赵宇峰（1987－　），男，博士，讲师，zhaoyufeng0613@126.com

通讯作者:
冯顺山（1952－　），男，硕士，教授，ssfeng@bit.edu.cn

中图分类号: O389; TJ55
计量
- 文章访问数: 256
- HTML全文浏览量: 105
- PDF下载量: 128
- 被引次数: 0
出版历程
- 收稿日期: 2024-03-07
- 修回日期: 2024-06-01
- 网络出版日期: 2024-06-06
- 刊出日期: 2024-11-15

A study of the failure of cased charge under impact of reactive fragments

1.
School of Aerospace Engineering, North University of China, Taiyuan 030051, Shanxi, China
2.
School of Mechanical and Electrical Engineering, North University of China, Taiyuan 030051, Shanxi, China
3.
State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China

摘要

摘要: 为了研究反应破片对带壳装药的冲击毁伤效应，通过弹道实验和AUTODYN有限元仿真，结合由等效破片初速和等效格尼速度表征的带壳装药各失效等级判据，获得并对比了惰性破片和反应破片冲击下带壳装药的等效破片初速、等效格尼速度、带壳装药的反应持续时间、鉴证靶破坏情况和炸药层峰值压力，分析了反应破片靶后释能特点对带壳装药失效的影响。结果表明：惰性破片可使带壳装药发生正常爆轰失效；反应破片穿靶后的动能与化学能叠加效应弱，只能使带壳装药发生爆燃失效或爆炸失效，带壳装药的等效格尼速度与格尼速度的比为0.014~0.233，炸药层峰值压力为1.04~3.62 GPa。
- 反应破片 /
- 等效破片初速 /
- 等效格尼速度 /
- 带壳装药 /
- 爆燃
Abstract: Reactive fragments are composed of multifunctional impact reactive structural materials. After reactive fragments penetrate the front target of warhead, the debris cloud generated by the sufficient reaction of reactive material will damage the medium behind the target in the form of kinetic energy-chemical energy coupling damage. Ballistic impact experiments and finite element simulations were conducted to investigate the impact damage effect of reactive fragments on cased charge. Based on the criteria for failure levels of cased charge characterized by equivalent fragments initial velocity and equivalent gurney velocity, the ratio of the equivalent gurney velocity under abnormal detonation conditions to gurney velocity or the ratio of the equivalent fragments initial velocity under abnormal detonation conditions to the fragments initial velocity is used to measure the reaction violence of the cased charge. Equivalent gurney velocity of cased charge under impact of inert fragments and reactive fragments, response duration of cased charge, the damage of the authentication target, and the peak pressure of explosive layer are compared. The influence of energy release characteristics of reactive fragments on the failure of cased charge is also analyzed. The results show that explosive detonate under the impact of inert fragments, while explosive deflagrate or explode under the impact of reactive fragments. The steel verification target only presents significant circular pit during explosive detonation. The explosive detonation process captured by high-speed photography is on the microsecond scale, while the explosive explosion or deflagration process is on the millisecond scale. Under the penetration of six reactive fragments, the corresponding ratio of equivalent gurney velocity to gurney velocity ranges from 0.014 to 0.233, which is far below the ratio of equivalent gurney velocity to gurney velocity under the condition of inert fragments penetrating cased charges. By using AUTODYN, the peak pressure at the observation point on the axis of the cased charge during detonation failure under the penetration of inert fragments ranges from 17.3 to 34.5 GPa, while the peak pressure of cased charge during deflagration failure under the penetration of reactive fragments ranges from 1.04 to 3.62 GPa, which is far below the critical detonation pressure. Based on the ratio of the equivalent gurney velocity to gurney velocity, the peak pressure of explosive and superimposed effect of kinetic energy and chemical energy of reactive fragments, the idea that it is difficult to detonate cased charge under the penetration of reactive fragments is proposed.
- reactive fragments /
- equivalent fragments initial velocity /
- equivalent gurney velocity /
- cased charge /
- deflagration

HTML全文

弹性动力学中将研究对象视为连续的弹性体, 运用数学物理的方法, 解析地研究机械运动, 不仅对数学物理的基础科学研究具有重要意义, 而且对岩土工程、地震工程等工程技术应用具有很大的价值。弹性波的传播及其在全空间中被一些特殊曲面散射的研究已经日趋完善^[1-4]。作为最简单的波动模式, 柱体在含复杂边界介质中的反平面散射研究取得了一定的进展^[5-6]。研究人员陆续研究了无限大弹性半空间、1/4空间中圆柱孔洞、夹杂、衬砌对SH平面波的稳态散射^[7-12]。双相介质半空间中SH波对圆柱孔洞、夹杂和界面裂纹的散射问题也得到了一定的研究^[13-18]。这些研究中主要考虑入射波与孔洞或夹杂在同一介质内的散射问题, 对于平面SH波从介质I穿过界面折射进入介质II, 透射波对介质II中圆孔的散射问题的研究相对匮乏。本文中用Green函数和积分方程, 按复变函数方法对其稳态响应进行分析, 求解双相介质半空间中位移场和应力场的解析表达式, 并采用数值方法计算算例中圆孔边沿动应力集中因子和地表位移幅值的分布。

1. 理论模型

如图 1所示, 连续、均匀、各向同性的弹性半空间, 由介质Ⅰ和Ⅱ组成, 其剪切模量分别为μ₁和μ₂, 质量密度分别为ρ₁和ρ₂, 两种介质的分界为B_V, 这种双相介质半空间的边界平面为B_H。介质Ⅱ中有一个圆形孔洞, 其半径为R, 其圆心O₁距离界面B_V为d, 距离边界B_H为h, 定义平面直角坐标系O₁x₁y₁(图 1), 以出平面方向为z轴正方向, x₁轴、y₁轴和z轴满足右手关系。

图 1 半空间双相介质中的圆孔

Figure 1. Circular cavity in a bi-material half space

下载: 全尺寸图片幻灯片

在介质Ⅰ中有稳态SH平面波入射, 其角频率为ω。入射波与它在边界B_H上的反射和在界面B_V的反射和透射, 造成了介质Ⅰ与介质Ⅱ中质点的反平面运动, 其反平面运动的位移分别为w_Ⅰ和w_Ⅱ, 满足标量波动方程:

$\mu_{1} \nabla w_{\mathrm{I}}=\rho_{1} \frac{\partial^{2} w_{\mathrm{I}}}{\partial t^{2}}, \quad \mu_{2} \nabla w_{\mathrm{II}}=\rho_{2} \frac{\partial^{2} w_{\mathrm{II}}}{\partial t^{2}}$

(1)

分离时间变量和空间变量, 得到对应的Helmholtz方程(介质Ⅰ和介质Ⅱ中位移场的控制方程):

$\nabla w_{\mathrm{I}}+k_{1}^{2} w_{\mathrm{I}}=0, \quad \nabla w_{\mathrm{II}}+k_{2}^{2} w_{\mathrm{II}}=0$

(2)

式中: , 分别为介质Ⅰ和介质Ⅱ中的波数。

位移场w_Ⅰ和w_Ⅱ在半空间边界B_H平面上产生的剪应力都为0, 满足Neumann边界条件; 位移场w_Ⅰ和对应的应力场τ_Ⅰ在界面B_V上与位移场w_Ⅱ和对应的应力场τ_Ⅱ相等, 分别满足Dirichlet边界条件和Neumann边界条件; 应力场τ_Ⅱ在圆形孔洞边界上为0, 满足Neumann边界条件。这4个边界条件与介质Ⅰ和介质Ⅱ中的控制方程(2)组成本文的定解问题。

2. 等价问题

为消除双相介质半空间边界B_H平面的影响, 引入等价问题。如图 2所示, 界面B_V将全空间分为介质Ⅰ和Ⅱ两个半空间。介质Ⅱ中有2个大小相等的圆形孔洞, 关于B_H平面对称, 圆心分别为O₁和O₂, 半径为R。以B_H和B_V平面的交线在坐标平面O₁x₁y₁内的投影O为原点, 建立平面直角坐标系Oxy, 同时以O₂为原点建立平面直角坐标系O₂x₂y₂。按复变函数方法, 引入复变量为虚数单位。对比图 1和图 2可得到复变量z, z₁, z₂的关系:

$\mathrm{z}_{1}=\mathrm{z}-\mathrm{d}+i \mathrm{~h}, \quad \mathrm{z}_{2}=\mathrm{z}-\mathrm{d}-i \mathrm{~h}$

(3)

图 2 等价问题

Figure 2. Equivalent issue

下载: 全尺寸图片幻灯片

在介质Ⅰ中有一对稳态的SH平面波入射, 其角频率均为ω, 初相位相等, 入射角关于平面B_H对称, 其位移w_i1和w_i2为:

$\left\{ \begin{array}{l} {w_{i1}}(z,t) = {{\rm{W}}_0}\exp \left[ {\frac{{{\rm{i}}{k_1}}}{2}\left( {z{{\rm{e}}^{ - {\rm{i}}\alpha }} + \bar z{{\rm{e}}^{{\rm{i}}\alpha }}} \right) - {\rm{i}}\omega t} \right]\\ {w_{{\rm{i}}2}}(z,t) = {W_0}\exp \left[ {\frac{{{\rm{i}}{k_1}}}{2}\left( {z{{\rm{e}}^{{\rm{i}}\alpha }} + \bar z{{\rm{e}}^{ - {\rm{i}}\alpha }}} \right) - {\rm{i}}\omega t} \right] \end{array} \right.$

(4)

式中:W₀为入射波的位移振幅, α为入射位移波w_i1传播方向与x轴的夹角, 为z的共轭复变量。

这样, 本文中定解条件变成等价的边值问题, 由B_V界面上位移和应力连续的Dirichlet和Neumann边界条件, 及介质Ⅱ中圆形孔洞边界B₁和B₂柱面上剪应力自由的Neumann边界条件组成。

介质Ⅰ中的入射位移波w_i1和w_i2在界面B_V将发生反射和折射, 折射波将透射进入介质Ⅱ, 在圆形孔洞的边界B₁和B₂柱面发生散射, 散射波也将在界面B_V发生反射和折射, 为了分析界面B_V对散射波的反射和折射效应, 引入Green函数方法。

3. Green函数

按照Green函数方法, 考虑点源荷载在界面B_V上分别对介质Ⅰ和介质Ⅱ的稳态位移响应, 其分离时间变量后, 即为Green函数。

在介质Ⅰ中, 点源荷载δ(z-iy₀)exp(-iωt)产生的Green函数G_Ⅰ(z, iy₀)是半空间中的基本解, 表达式为

$G_{\mathrm{I}}\left(z, \mathrm{i} y_{0}\right)=\frac{\mathrm{i}}{2 \mu_{1}} H_{0}^{(1)}\left(k_{1}\left|z-\mathrm{i} y_{0}\right|\right)$

(5)

式中:iy₀是坐标面Oxy上点源荷载作用点对应的复数值, 位于界面B_V上, 与x轴的距离为为0阶第一种Hankel函数, Hankel函数是第三类Bessel函数。

如图 3所示, 在介质Ⅱ中, 点源荷载δ(z-iy₀)exp(-iωt)产生的Green函数G_Ⅱ(z, iy₀), 按照叠加原理, 有:

$G_{\text {II }}\left(z, \mathrm{i} y_{0}\right)=G\left(z, \mathrm{i} y_{0}\right)+w_{\mathrm{Gl}}(z)+w_{\mathrm{Gl}}^{\prime}(z)+ \\ \ \ \ \ \ \ w_{\mathrm{G} 2}(z)+w_{\mathrm{G} 2}^{\prime}(z)$

(6)

$G\left(z, \mathrm{i} y_{0}\right)=\frac{\mathrm{i}}{2 \mu_{2}} H_{0}^{(1)}\left(k_{2}\left|z-\mathrm{i} y_{0}\right|\right)$

(7)

图 3 介质Ⅱ中的Green函数

Figure 3. Green function in medium Ⅱ

下载: 全尺寸图片幻灯片

式中:G(z, iy₀)是介质Ⅱ中没有孔洞时, 半空间中的基本解; w_G1(z)和w_G2(z)分别是点源荷载δ(z-iy₀)对圆形孔洞边界B₁和B₂柱面产生的散射波的位移场, G2(z)是w_G1(z)和w_G2(z)关于界面B_V对称的镜像散射波, 按照波函数展开法, 它们满足:

$\left\{\begin{array}{l} w_{\mathrm{G} 1}(z)=\sum_{n=-\infty}^{+\infty} A_{n} H_{n}^{(1)}\left(k_{2}\left|z_{1}\right|\right)\left(\frac{z_{1}}{\left|z_{1}\right|}\right)^{n} \\ w_{\mathrm{G} 2}(z)=\sum_{n=-\infty}^{+\infty} B_{n} H_{n}^{(1)}\left(k_{2}\left|z_{2}\right|\right)\left(\frac{z_{2}}{\left|z_{2}\right|}\right)^{n} \\ w_{\mathrm{G} 1}^{\prime}(z)=\sum_{n=-\infty}^{+\infty}(-1)^{n} A_{n} H_{n}^{(1)}\left(k_{2}\left|z_{1}^{\prime}\right|\right)\left(\frac{z_{1}^{\prime}}{\left|z_{1}^{\prime}\right|}\right)^{-n} \\ w_{\mathrm{G} 2}^{\prime}(z)=\sum_{n=-\infty}^{+\infty}(-1)^{n} B_{n} H_{n}^{(1)}\left(k_{2}\left|z_{2}^{\prime}\right|\right)\left(\frac{z_{2}^{\prime}}{\left|z_{2}^{\prime}\right|}\right)^{-n} \end{array}\right.$

(8)

式中: 是n阶第一种Hankel函数, A_n和B_n分别是圆形孔洞边界B₁和B₂柱面散射波的Fourier-Hankel波函数级数展开式的待定系数, z′₁=z₁+2d, z′₂=z₂+2d, 分别为以圆心O₁和O₂关于界面B_V的镜像点O′₁和O′₂建立的直角坐标系O′₁x′₁y′₁和O′₂x′₂y′₂的复数形式。

介质Ⅱ对应半空间中的基本解G(z, iy₀), 位移波w_G1和w′_G1的矢量和, 位移波w_G2和w′_G2的矢量和, 它们的应力场都满足界面B_V上的剪应力自由条件。位移场G_Ⅱ(z, iy₀)的应力场必须满足圆形孔洞边界B₁和B₂柱面上的剪应力自由条件做Fourier展开, 得到等价的剪应力自由条件, 用于确定Fourier-Hankel波函数系数A_n和B_n, 即定解条件。根据三角函数系的正交性, 得到和定解条件等价的两组方程。截断Fourier-Hankel波函数级数展开式, 取对应的方程, 组成定解方程组, 按Gauss消去法解得A_n和B_n的具体数值, 这样, 就确定了介质Ⅱ中Green函数G_Ⅱ(z, iy₀)位移场的具体表达式。

4. 位移波场

按照分区契合的方法, 研究SH波在双相介质中的传播与散射。界面B_V将双相介质分为介质Ⅰ和介质Ⅱ, 分别研究介质Ⅰ和介质Ⅱ, 定义介质Ⅰ的边界为平面, 介质Ⅱ的边界为平面, 在和平面上分别加以沿y轴分布的方向相反、大小相等的一组力f(y), 表征介质Ⅰ和介质Ⅱ的相互作用力, 如图 4所示, 其中点表示垂直纸面向外, 叉表示垂直纸面向里。

图 4 位移波场的契合

Figure 4. Conjunction of displacement wave field

下载: 全尺寸图片幻灯片

介质Ⅰ中的位移场w_Ⅰ(z)由3部分组成:入射波的位移场、反射波的位移场、介质Ⅱ对介质Ⅰ的作用力f(y)产生的位移场, 其表达式为:

$w_{\mathrm{I}}(z)=w_{i 1}(z)+w_{i 2}(z)+w_{r 1}(z)+$

$w_{r 2}(z)+\int_{-\infty}^{+\infty} G_{\mathrm{I}}\left(z, \mathrm{i} y_{0}\right) f\left(y_{0}\right) d y_{0}$

(9)

式中:w_i1(z)和w_i2(z)是入射波的位移场, 根据式(4), 略去时间谐和项exp(-iωt), 得到入射波位移场的表达式:

$\left\{\begin{array}{l} w_{\mathrm{i} 1}(z)=W_{0} \exp \left[\frac{\mathrm{i} k_{1}}{2}\left(z \mathrm{e}^{-\mathrm{i} \alpha}+\bar{z} \mathrm{e}^{\mathrm{i} \alpha}\right)\right. \\ w_{\mathrm{i} 2}(z)=W_{0} \exp \left[\frac{\mathrm{i} k_{1}}{2}\left(z \mathrm{e}^{\mathrm{i} \alpha}+\bar{z} \mathrm{e}^{-\mathrm{i} \alpha}\right)\right. \end{array}\right.$

(10)

w_r1(z)和w_r2(z)是反射波的位移场, 表达式为:

$\left\{\begin{array}{l} w_{\mathrm{rl}}(z)=W_{1} \exp \left[\frac{\mathrm{i} k_{1}}{2}\left(z \mathrm{e}^{-\mathrm{i}(\pi-\alpha)}+\bar{z} \mathrm{e}^{\mathrm{i}(\pi-\alpha)}\right)\right. \\ w_{\mathrm{r} 2}(z)=W_{1} \exp \left[\frac{\mathrm{i} k_{1}}{2}\left(z \mathrm{e}^{\mathrm{i}(\pi-\alpha)}+\bar{z} \mathrm{e}^{-\mathrm{i}(\pi-\alpha)}\right)\right. \end{array}\right]$

(11)

式中:W₁是反射波的振幅。讨论介质Ⅰ边界平面上的分布力f(y)的影响, 当y=y₀时, 即取坐标面Oxy上复坐标iy₀点, 加以集中力f(y₀)δ(z-iy₀), G_Ⅰ(z, iy₀)是介质Ⅰ中产生的Green函数, G_Ⅰ与f(y₀)乘积的积分项是介质Ⅱ对介质Ⅰ的作用力在介质Ⅰ中产生的位移场。

介质Ⅱ中的位移场w_Ⅱ(z)同样由3部分组成:折射波的位移场、散射波的位移场, 以及介质Ⅰ对介质Ⅱ的作用力-f(y)产生的位移场, 表达式为:

$w_{\mathrm{II}}(z)=w_{\mathrm{t} 1}(z)+w_{\mathrm{t} 2}(z)+w_{\mathrm{s} 1}(z)+w_{\mathrm{s} 1}^{\prime}(z)+$

$w_{\mathrm{s} 1}(z)+w_{\mathrm{s} 2}^{\prime}(z)-\int_{-\infty}^{+\infty} G_{\mathrm{II}}\left(z, \mathrm{i} y_{0}\right) f\left(y_{0}\right) d y_{0}$

(12)

w_t1(z)和w_t2(z)是折射波的位移场, 表达式为:

$\left\{\begin{array}{l} w_{\mathrm{t} 1}(z)=W_{2} \exp \left[\frac{\mathrm{i} k_{2}}{2}\left(z \mathrm{e}^{-\mathrm{i} \beta}+\bar{z} \mathrm{e}^{\mathrm{i} \beta}\right)\right. \\ w_{\mathrm{t} 2}(z)=W_{2} \exp \left[\frac{\mathrm{i} k_{2}}{2}\left(z \mathrm{e}^{\mathrm{i} \beta}+\bar{z} \mathrm{e}^{-\mathrm{i} \beta}\right)\right. \end{array}\right]$

(13)

式中:W₂是折射波的振幅, β是折射波的角度(见图 4)。按照界面B_V上位移和应力的连续性条件, 角度α和β满足Snell定律:

k₁sinα=k₂sinβ

入射波振幅W₀、反射波振幅W₁和折射波振幅W₂, 满足

$\left\{\begin{array}{l} W_{0}+W_{1}=W_{2} \\ W_{0} k_{1} \mu_{1} \cos \alpha-W_{1} k_{1} \mu_{1} \cos \alpha=W_{2} k_{2} \mu_{2} \cos \beta \end{array}\right.$

(15)

w_s1(z)和w_s2(z)是介质Ⅱ中圆形孔洞边界B₁和B₂柱面对折射波的散射所产生的位移场, w′_s1(z)和w′_s2(z)是散射波w_s1(z)和w_s2(z)分别对界面B_V反射产生的镜像散射波的位移场, 按照波函数展开法, 得到:

$\left\{\begin{array}{l} w_{\mathrm{s} 1}(z)=\sum\limits_{n=-\infty}^{+\infty} C_{n} H_{n}^{(1)}\left(k_{2}\left|z_{1}\right|\right)\left(\frac{z_{1}}{\left|z_{1}\right|}\right)^{n} \\ w_{\mathrm{s} 2}(z)=\sum\limits_{n=-\infty}^{+\infty} D_{n} H_{n}^{(1)}\left(k_{2}\left|z_{2}\right|\right)\left(\frac{z_{2}}{\left|z_{2}\right|}\right)^{n} \\ w_{\mathrm{s} 1}^{\prime}(z)=\sum\limits_{n=-\infty}^{+\infty}(-1)^{n} C_{n} H_{n}^{(1)}\left(k_{2}\left|z_{1}^{\prime}\right|\right)\left(\frac{z_{1}^{\prime}}{\left|z_{1}^{\prime}\right|}\right)^{-n} \\ w_{\mathrm{s} 2}^{\prime}(z)=\sum\limits_{n=-\infty}^{+\infty}(-1)^{n} D_{n} H_{n}^{(1)}\left(k_{2}\left|z_{2}^{\prime}\right|\right)\left(\frac{z_{2}^{\prime}}{\left|z_{2}^{\prime}\right|}\right)^{-n} \end{array}\right.$

(16)

式中:C_n和D_n是待定的Fourier-Hankel波函数系数。讨论介质Ⅱ边界平面上分布力-f(y)的影响, 当y=y₀时, 取坐标面Oxy上复坐标iy₀点, 加以集中力-f(y₀)δ(z-iy₀), G_Ⅱ(z, iy₀)是介质Ⅱ中产生的Green函数, G_Ⅱ与-f(y₀)乘积的积分项是介质Ⅰ对介质Ⅱ的作用力在介质Ⅱ中产生的位移场。

散射波w_s1和w_s2与镜像散射波w′_s1和w′_s2的矢量和满足平面上应力自由的条件, 其与折射波w_t1和w_t1的矢量和应满足圆形孔洞边界B₁和B₂柱面上的应力自由条件。做Fourier级数展开, 得到确定系数C_n和D_n的定解方程组; 相应地截断定解方程组和散射波的Fourier-Hankel波函数级数, 得到C_n和D_n的具体形式, 这样就求得了介质Ⅱ中由折射波在圆形孔洞边界上所形成的散射波。

5. 积分方程

由(9)和(12)式, 按照界面B_V上位移和应力的连续条件, 即契合条件, 得到:

$w_{\mathrm{I}}(\mathrm{i} y)=w_{\mathrm{II}}(\mathrm{i} y), \quad \tau_{\mathrm{I}}(\mathrm{i} y)=\tau_{\mathrm{II}}(\mathrm{i} y)$

(17)

式中:w_Ⅰ(iy)和w_Ⅱ(iy)分别是介质Ⅰ和介质Ⅱ中的位移沿界面B_V的分布, τ_Ⅰ(iy)和τ_Ⅱ(iy)分别是介质Ⅰ和介质Ⅱ中应力分量沿界面B_V的分布。作用在和平面上的一组大小相等、方向相反的力f(y), 本身代表介质Ⅰ和介质Ⅱ之间的相互作用力, 满足应力的契合条件。将第4节中构造的位移波的位移场和应力场代入方程(17), 入射波和反射波在介质Ⅰ中的位移场和应力场与折射波在介质Ⅱ中的位移场和应力场满足界面B_V上的契合条件, 而介质Ⅰ中的Green函数、介质Ⅱ中的Green函数、散射波与其镜像散射波的合位移波在界面B_V上剪应力均为0, 也满足界面B_V上的应力契合条件。因此, 界面B_V上的契合条件(17)等价于介质Ⅰ和介质Ⅱ中的位移连续性条件:

$\int_{-\infty}^{+\infty}\left[G_{\mathrm{I}}\left(\mathrm{i} y, \mathrm{i} y_{0}\right)+G_{\mathrm{II}}\left(\mathrm{i} y, \mathrm{i} y_{0}\right)\right] f\left(y_{0}\right) d y_{0}=w_{\mathrm{s} 1}(\mathrm{i} y)+w_{\mathrm{sl}}^{\prime}(\mathrm{i} y)+w_{\mathrm{s} 2}(\mathrm{i} y)+w_{\mathrm{s} 2}^{\prime}(\mathrm{i} y)$

(18)

这是积分限为无穷的第一种Fredholm积分方程。在y轴上对界面B_V取N+1个节点y_k, 同样地, 对积分区间取N+1个节点y_l, 按数值方法对(18)式进行离散求解:

$\sum_{l=1}^{N+1} Q_{l}\left[G_{\mathrm{I}}\left(\mathrm{i} y_{k}, \mathrm{i} y_{l}\right)+G_{\mathrm{II}}\left(\mathrm{i} y_{k}, \mathrm{i} y_{l}\right)\right] f\left(y_{l}\right) \approx w_{\mathrm{s} 1}\left(\mathrm{i} y_{k}\right)+w_{\mathrm{s} 1}^{\prime}\left(\mathrm{i} y_{k}\right)+w_{\mathrm{s} 2}\left(\mathrm{i} y_{k}\right)+w_{\mathrm{s} 2}^{\prime}\left(\mathrm{i} y_{k}\right)$

(19)

式中:Q_l是求积系数, 伴随求积节点y_l, 取决于数值积分所采用的方法。按Gauss消去法, 计算节点y_l上f(y_l)的数值, 之后用Lagrange插值拟合界面B_V上的分布力f(y)

6. 算例分析

按照第4节中的讨论, 结合第5部分求得的界面B_V上的分布力f(y), 分别按照(9)式和(12)式, 得到介质Ⅰ和介质Ⅱ中的位移场w_Ⅰ(z)和w_Ⅱ(z), 求得应力场τ_Ⅰ(z)和τ_Ⅱ(z), 它们足以描述介质Ⅰ和介质Ⅱ中全部质点的运动状态。

为直观的讨论圆孔散射所造成的动应力集中效应, 引入动应力集中因子(z)定量表征介质Ⅱ中的圆孔对透射SH波的散射:

$\tau_{\theta \xi}^{*}(z)=\frac{\left|\tau_{\theta \xi \text { І }}(z)\right|}{W_{2} k_{2} \mu_{2}}$

(20)

$\tau_{\theta \xi 一{\mathrm{II}}}(z)=\frac{1}{\left|z_{1}\right|}\left[\frac{\partial w_{\mathrm{II}}(z)}{\partial z_{1}} z_{1}+\frac{\partial w_{\mathrm{II}}(z)_{-}}{\partial \bar{z}_{1}} z_{1}\right]$

(21)

式中:τ_θξⅡ是介质Ⅱ中应力张量场τ_Ⅱ的周向分量。是介质Ⅱ中z点在O₁x₁y₁坐标系中的周向应力分量的时间幅值与透射位移波w_t1在同一点产生的最大动应力的比值, 当|z₁|=R时, τ_θξⅡ(z)是B₁柱面上的应力, 是B₁柱面上的动应力集中因子。

为直观的讨论圆孔散射对水平边界B_H位移的影响, 引入标准化的位移幅值W^*(z):

$W^{*}(z)=\left\{\begin{array}{l} \left|w_{\mathrm{I}}(z)\right| / W_{0}, z \in \text { Medium } \mathrm{I} \\ \left|w_{\mathrm{II}}(z)\right| / W_{0}, z \in \text { Medium } \end{array}\right.$

(22)

水平边界B_H平面的存在极大地影响了圆孔的散射, 包括2部分:与圆孔到B_H平面的距离h无关的镜像透射位移波w_t2(z)对圆孔的散射, 以及随着h的增大而衰减的镜像散射波w_s2(z)对圆孔的散射。水平边界对圆孔散射的影响在弹性动力学反平面半空间问题中得到了广泛的讨论, 相较而言, 在双相介质半空间中, 垂直界面B_V对介质Ⅱ中圆孔散射的作用更值得注意。

算例以圆孔半径R为参考长度, 给定h=2R, d=2R, 取入射角α=0, 数值计算双相介质的不同波数、剪切模量对圆孔散射的具体影响, 定义波数比k^*=k₁/k₂, 剪切模量比μ^*=μ₁/μ₂。图 5和图 6分别给出由圆孔散射所造成的圆孔边沿的动应力集中因子的分布和受其影响的水平边界B_H平面上标准化的位移幅值的分布。

图 5 动应力集中因子在圆孔边沿的分布

Figure 5. Distribution of dynamic stress concentrations factors around the cylindrical cavity

下载: 全尺寸图片幻灯片

图 6 标准化位移幅值在水平边界上的分布

Figure 6. Distribution of normalized displacement amplitude along the horizontal bound

下载: 全尺寸图片幻灯片

7. 结论

双相介质半空间中圆孔的反平面稳态散射受到很多因素的影响, 诸如, 圆孔和边界的距离、和界面的距离, 双相介质的波数比、剪切模量比, 平面SH波的入射角度、频率, 等等。不同的圆孔位置对散射的影响可以参照1/4空间圆孔问题的研究, 当圆孔距离界面足够远, 界面另一侧介质的影响可以忽略, 界面附近的圆孔, 对穿过界面的透射平面SH波的散射受到界面另一侧介质的强烈影响, 集中体现在介质的波数和剪切模量。作为SH波在双相介质中传播的问题, 应充分考虑界面两侧介质对平面SH波透射的影响, 给定入射角度和频率。

在相对低频的透射波作用下, 圆孔散射造成更大的动应力集中, 在相对高频的透射波作用下, 圆孔散射造成的动应力分布受到界面另一侧介质的影响更大。在一定范围内, 界面另一侧介质的剪切模量很大程度地决定着圆孔散射所造成的动应力的数值与分布, 这种效应随着界面另一侧介质的质量密度的增大而减小。在相对高频的透射波作用下, 圆孔散射造成局部更大的地表位移, 也更受界面另一侧介质的影响, 随着界面另一侧介质波数的增大, 圆孔散射所造成的局部地表位移的变化也越大, 这种效应随着界面另一侧介质的剪切模量的增大也略微地增大。

图 1 轴对称等壁厚带壳装药结构一端轴线起爆条件下破片初速随弹轴的分布

Figure 1. Axial distribution of fragments velocity of axial symmetry and equal wall thickness cased charge under one axis detonating condition

下载: 全尺寸图片幻灯片

图 2 实验布局示意图及现场布置

Figure 2. Schematic of terminal ballistic experimental setup and experiment environment

下载: 全尺寸图片幻灯片

图 3 带壳装药

Figure 3. Experimental target

下载: 全尺寸图片幻灯片

图 4 仿真模型示意图

Figure 4. Schematics of simulation models

下载: 全尺寸图片幻灯片

图 5 破片冲击带壳装药后的钢鉴证靶损伤

Figure 5. Damage of steel identification target after fragments impacting cased charge

下载: 全尺寸图片幻灯片

图 6 破片冲击带壳装药后试件失效的典型高速摄影照片

Figure 6. Typical high-speed photography of specimen failure after fragments impacting cased charge

下载: 全尺寸图片幻灯片

图 7 破片冲击带壳装药仿真结果

Figure 7. Simulation results of fragments impacting cased charge

下载: 全尺寸图片幻灯片

图 8 惰性破片冲击带壳装药时壳体观测点的峰值压力时程曲线

Figure 8. Peak pressure histories at the observation points of the casing when inert fragments impacting cased charge

下载: 全尺寸图片幻灯片

图 9 反应破片作用下炸药各观测点的峰值压力

Figure 9. Peak pressure at observation points of explosives under the impact of reactive fragments

下载: 全尺寸图片幻灯片

图 10 破片动能时程曲线

Figure 10. Kinetic energy histories of the fragments

下载: 全尺寸图片幻灯片

表 1 反应破片参数

Table 1. Parameters of reactive fragments

破片序号	破片类型	破片尺寸/mm	材料组成	质量/g
1	惰性破片	$\varnothing$ 12×12	45钢	10.44
2	反应破片	$\varnothing$ 12×18	DU/PTFE	16.46
3	反应破片	$\varnothing$ 12×18	Al/PTFE	4.32
4	反应破片	$\varnothing$ 12×18	DU/PTFE	16.27
5	反应破片	$\varnothing$ 12×18	DU/PTFE	7.32
6	反应破片	$\varnothing$ 12×18	DU/PTFE	7.40
7	反应破片	$\varnothing$ 12×18	DU/PTFE	12.53

下载: 导出CSV

表 2 Al/PTFE材料的强度模型参数和状态方程参数^[19]

Table 2. Strength model parameters and equation of state parameters of Al/PTFE^[19]

强度模型参数						状态方程参数
A₀/MPa	B₀/MPa	n	C₀	m	T_m/K	C₁/(m·s⁻¹)	S₁	γ₀
8.044	250.6	1.8	0.4	1.426	500	1450	2.2584	0.9

下载: 导出CSV

表 3 反应破片和惰性破片冲击带壳装药的实验结果对比

Table 3. Comparison of experimental results between reactive fragments and inert fragments impacting cased charge

破片序号	着靶速度/(m·s⁻¹)	v_xSi/(m·s⁻¹)	$\sqrt{{{2E_{{\mathrm{S}}_i}}}}$ /(m·s⁻¹)	k_Ei	鉴证靶破坏情况	炸药失效等级
1	1210	1198	2598.7	0.964	圆形凹坑	正常爆轰
2	934	290	629.1	0.233	无凹坑	爆炸
3	1298	17	36.9	0.014	无凹坑	爆燃
4	869	186	403.5	0.160	无凹坑	爆燃
5	888	159	344.9	0.128	无凹坑	爆燃
6	875	26	56.4	0.021	无凹坑	爆燃
7	927	132	286.3	0.106	无凹坑	爆燃

下载: 导出CSV

表 4 惰性破片和反应破片冲击带壳装药的实验与仿真结果对比

Table 4. Comparison of experimental and simulation results of inert fragments and reactive fragments impacting cased charge

破片序号	v_xSi			$\sqrt{{{2E_{{\mathrm{S}}_i}}}}$			k_Ei
破片序号	实验/(m·s⁻¹)	仿真/(m·s⁻¹)	误差/%	实验/(m·s⁻¹)	仿真/(m·s⁻¹)	误差/%	实验	仿真	误差/%
1	1198.0	1211.3	1.11	2598.7	2627.5	1.11	0.9640	0.9740	1.04
3	17.0	18.2	7.06	36.9	39.5	7.04	0.0137	0.0146	6.57

下载: 导出CSV

参考文献(24)

[1]	FENG S S, WANG C L, HUANG G Y. Experimental study on the reaction zone distribution of impact-induced reactive materials [J]. Propellants, Explosives, Pyrotechnics, 2017, 42(8): 896–905. DOI: 10.1002/prep.201600274.
[2]	冯顺山, 李伟, 周彤, 等. 内置冲击反应材料弹丸壳体侵彻损伤效应研究 [J]. 兵工学报, 2016, 37(S2): 61–68. FENG S S, LI W, ZHOU T, et al. Research on damage effect of a warhead filled with reactive materials [J]. Acta Armamentarii, 2016, 37(S2): 61–68.
[3]	顾阳晨, 王金相, 陈兴旺, 等. 高速动能破片和包覆活性材料对屏蔽装药的串联毁伤效应 [J]. 含能材料, 2021, 29(7): 607–616. DOI: 10.11943/CJEM2020328. GU Y C, WANG J X, CHEN X W, et al. Tandem damage effect of high-speed kinetic fragments and coated active materials on shielded charges [J]. Chinese Journal of Energetic Materials, 2021, 29(7): 607–616. DOI: 10.11943/CJEM2020328.
[4]	赵宏伟, 余庆波, 邓斌, 等. 活性破片终点毁伤威力试验研究 [J]. 北京理工大学学报, 2020, 40(4): 375–381. DOI: 10.15918/j.tbit1001-0645.2019.092. ZHAO H W, YU Q B, DENG B, et al. Experimental study on terminal demolition lethality of reactive fragments [J]. Transactions of Beijing Institute of Technology, 2020, 40(4): 375–381. DOI: 10.15918/j.tbit1001-0645.2019.092.
[5]	罗普光. 钨锆含能破片冲击反应特性及其对屏蔽炸药的毁伤研究 [D]. 北京: 北京理工大学, 2017: 72–78. LUO P G. Study on impact-initiated characters of W/Zr energetic fragments and damage effect to shielded explosives [D]. Beijing: Beijing Institute of Technology, 2017: 72–78.
[6]	周杰, 何勇, 何源, 等. 含能毁伤元冲击引爆模拟战斗部试验研究 [J]. 含能材料, 2016, 24(11): 1048–1056. DOI: 10.11943/J.ISSN.1006-9941.2016.11.003. ZHOU J, HE Y, HE Y, et al. Experimental study on shock initiation of simulative warhead by energetic kill element [J]. Chinese Journal of Energetic Materials, 2016, 24(11): 1048–1056. DOI: 10.11943/J.ISSN.1006-9941.2016.11.003.
[7]	梁君夫. 活性破片作用屏蔽装药引爆增强效应研究 [D]. 北京: 北京理工大学, 2016: 63–68. LIANG J F. Research on enhanced initiation behavior of reactive material fragment impacting covered explosive [D]. Beijing: Beijing Institute of Technology, 2016: 63–68.
[8]	李旭锋, 李向东, 顾文彬, 等. 含能破片引爆带壳炸药过程的数值模拟 [J]. 爆炸与冲击, 2014, 34(2): 202–208. DOI: 10.11883/1001-1455(2014)02-0202-07. LI X F, LI X D, GU W B, et al. Numerical simulation on detonating shelled explosives by energetic fragments [J]. Explosion and Shock Waves, 2014, 34(2): 202–208. DOI: 10.11883/1001-1455(2014)02-0202-07.
[9]	王海福, 郑元枫, 余庆波, 等. 活性破片引爆屏蔽装药机理研究 [J]. 北京理工大学学报, 2012, 32(8): 786–789, 823. DOI: 10.3969/j.issn.1001-0645.2012.08.004. WANG H F, ZHENG Y F, YU Q B, et al. Study on initiation mechanism of reactive fragment to covered explosive [J]. Transactions of Beijing Institute of Technology, 2012, 32(8): 786–789, 823. DOI: 10.3969/j.issn.1001-0645.2012.08.004.
[10]	何源, 何勇, 潘绪超, 等. 含能破片冲击引爆屏蔽炸药研究 [J]. 南京理工大学学报, 2011, 35(2): 187–193. DOI: 10.3969/j.issn.1005-9830.2011.02.008. HE Y, HE Y, PAN X C, et al. Initiation of shielded high explosive impacted by energetic fragment [J]. Journal of Nanjing University of Science and Technology, 2011, 35(2): 187–193. DOI: 10.3969/j.issn.1005-9830.2011.02.008.
[11]	叶小军, 李向东. 含能破片撞击引燃屏蔽炸药的实验研究 [J]. 弹箭与制导学报, 2009, 29(6): 131–134. DOI: 10.3969/j.issn.1673-9728.2009.06.036. YE X J, LI X D. Experimental study on reactive fragments ignited charge covered with a metal plateafter the impact [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2009, 29(6): 131–134. DOI: 10.3969/j.issn.1673-9728.2009.06.036.
[12]	GURNEY R W. The initial velocities of fragments from bombs, shell, and grenades [R]. Aberdeen: Ballistic Research Laboratory, 1943.
[13]	LI W, HUANG G Y, FENG S S. Effect of eccentric edge initiation on the fragment velocity distribution of a cylindrical casing filled with charge [J]. International Journal of Impact Engineering, 2015, 80: 107–115. DOI: 10.1016/j.ijimpeng.2015.01.007.
[14]	冯顺山, 蒋建伟, 何顺录, 等. 偏轴心起爆破片初速径向分布规律研究 [J]. 兵工学报, 1993, 14(S1): 12–16. FENG S S, JIANG J W, HE S L, et al. On the pattern of radial distribution pattern of initial velocities of fragments under asymmetrical initiation [J]. Acta Armamentarii, 1993, 14(S1): 12–16.
[15]	冯顺山, 崔秉贵. 战斗部破片初速轴向分布规律的实验研究 [J]. 兵工学报, 1987, 8(4): 60–63. FENG S S, CUI B G. An experimental investigation for the axial distribution of initial velocity of shells [J]. Acta Armamentarii, 1987, 8(4): 60–63.
[16]	冯顺山, 赵宇峰, 边江楠, 等. 动能侵彻体冲击带壳炸药装药的爆燃失效 [J]. 含能材料, 2019, 27(6): 521–527. DOI: 10.11943/CJEM2018215. FENG S S, ZHAO Y F, BIAN J N, et al. Deflagration failure of explosive cased charge under impact of kinetic energy penetrators [J]. Chinese Journal of Energetic Materials, 2019, 27(6): 521–527. DOI: 10.11943/CJEM2018215.
[17]	Naval Sea Systems Command. Hazard assessment tests for non-nuclear munitions: MIL-STD-2105E [S]. Washington: Naval Sea Systems Command, 2022.
[18]	孙业斌. 爆炸作用与装药设计 [M]. 北京: 国防工业出版社, 1987: 76–78. SUN Y B. Explosion effect and charge design [M]. Beijing: National Defense Industry Press, 1987: 76–78.
[19]	RAFTENBERG M N, MOCK W, KIRBY G C. Modeling the impact deformation of rods of a pressed PTFE/Al composite mixture [J]. International Journal of Impact Engineering, 2008, 35(12): 1735–1744. DOI: 10.1016/j.ijimpeng.2008.07.041.
[20]	陈刚, 陈忠富, 徐伟芳, 等. 45钢的J-C损伤失效参量研究 [J]. 爆炸与冲击, 2007, 27(2): 131–135. DOI: 10.11883/1001-1455(2007)02-0131-05. CHEN G, CHEN Z F, XU W F, et al. Investigation on the J-C ductile fracture parameters of 45 steel [J]. Explosion and Shock Waves, 2007, 27(2): 131–135. DOI: 10.11883/1001-1455(2007)02-0131-05.
[21]	陈刚, 陈忠富, 陶俊林, 等. 45钢动态塑性本构参量与验证 [J]. 爆炸与冲击, 2005, 25(5): 451–456. DOI: 10.11883/1001-1455(2005)05-0451-06. CHEN G, CHEN Z F, TAO J L, et al. Investigation and validation on plastic constitutive parameters of 45 steel [J]. Explosion and Shock Waves, 2005, 25(5): 451–456. DOI: 10.11883/1001-1455(2005)05-0451-06.
[22]	吴曼林, 刘玉存. 冲击波感度试验(SSGT)的数值模拟 [J]. 火工品, 2004(2): 16–19. DOI: 10.3969/j.issn.1003-1480.2004.02.005. WU M L, LIU Y C. Numerical modeling of shock sensitivity experiments (SSGT) [J]. Initiators & Pyrotechnics, 2004(2): 16–19. DOI: 10.3969/j.issn.1003-1480.2004.02.005.
[23]	章冠人, 陈大年. 凝聚炸药起爆动力学 [M]. 北京: 国防工业出版社, 1991: 109. ZHANG G R, CHEN D N. Initiation dynamics of condensed explosives [M]. Beijing: National Defense Industry Press, 1991: 109.
[24]	王成龙, 黄广炎, 冯顺山. 贫铀材料冲击破碎和能量释放后效研究 [C]//战斗部与毁伤技术学术交流会. 北京, 2017: 785–791. WANG C L, HUANG G Y, FENG S S. Research on impact crushing and energy release aftereffects of depleted uranium materials [C]//Symposium on Warheads and Damage Technology. Beijing, 2017: 785–791.

施引文献

资源附件(0)

访问统计

图(10) / 表(4)

计量

文章访问数: 256
HTML全文浏览量: 105
PDF下载量: 128
被引次数: 0

1. 理论模型
2. 等价问题
3. Green函数
4. 位移波场
5. 积分方程
6. 算例分析
7. 结论

反应破片对带壳装药的冲击效应

doi: 10.11883/bzycj-2024-0063

作者简介: 赵宇峰（1987－ ），男，博士，讲师，zhaoyufeng0613@126.com

通讯作者: 冯顺山（1952－ ），男，硕士，教授，ssfeng@bit.edu.cn

计量

出版历程

A study of the failure of cased charge under impact of reactive fragments

1. 理论模型

2. 等价问题

3. Green函数

4. 位移波场

5. 积分方程

6. 算例分析

7. 结论

计量

出版历程

目录

1. 理论模型

2. 等价问题

3. Green函数

4. 位移波场

5. 积分方程

6. 算例分析

7. 结论

作者简介:
赵宇峰（1987－　），男，博士，讲师，zhaoyufeng0613@126.com

通讯作者:
冯顺山（1952－　），男，硕士，教授，ssfeng@bit.edu.cn