活性材料冲击释能行为研究进展

汪德武; 任柯融; 江增荣; 赵宏伟; 陈荣; 郭宝月

doi:10.11883/bzycj-2020-0337

活性材料冲击释能行为研究进展

doi: 10.11883/bzycj-2020-0337

1.
中国人民解放军96901部队，北京 100094
2.
空军工程大学航空机务士官学校，河南信阳 464001
3.
国防科技大学文理学院，湖南长沙 410007

基金项目: 国家自然科学基金（11672328）

详细信息

作者简介:
汪德武（1967－　），男，博士，研究员，Dividwwyy@sohu.com

通讯作者:
陈　荣（1981－　），男，博士，副教授，r_chen@nudt.edu.cn

中图分类号: O389
计量
- 文章访问数: 1334
- HTML全文浏览量: 606
- PDF下载量: 245
- 被引次数: 14
出版历程
- 收稿日期: 2020-09-22
- 修回日期: 2020-11-05
- 网络出版日期: 2021-03-05
- 刊出日期: 2021-03-10

Shock-induced energy release behaviors of reactive materials

1.
96901 Unit of PLA, Beijing 100094, China
2.
Aviation Maintenance NCO Academy, Air Force Engineering University, Xinyang 464001, Henan, China
3.
College of Liberal Arts and Sciences, National University of Defense Technology, Changsha 410007, Hunan, China

摘要

摘要: 活性材料是一种具备释能特性的新型材料，其在冲击导致的高压/高温作用下可以发生化学反应，释放大量的化学能，因此在破片、聚能破甲战斗部等军事领域有广泛的应用潜力。为了实现对活性材料释能过程的设计与控制，推进活性材料武器化应用进程，就必须解答活性材料冲击释能行为中所包含的一系列复杂的力-热-化耦合问题。近40年来，对活性材料的冲击释能行为已开展了大量研究，本文在此基础上系统梳理了活性材料的冲击诱发化学反应机理、动力学以及相关效应的研究现状，重点关注活性材料的冲击释能实验表征技术、冲击诱发化学反应理论模型以及考虑力-热-化耦合的冲击压缩数值模拟方法等3方面的研究进展。总结认为，对活性材料冲击释能行为的研究已经具有一定的积淀，但目前对实验中超快化学反应行为的实时诊断研究还缺乏更加丰富、精细、直观的表征与探索，相关理论与数值模拟研究尚未建立能够完整描述活性材料冲击释能行为的力-热-化理论模型，缺乏能够从宏观尺度描述冲击释能行为的有效方法。因此，超快化学反应实验表征技术、宏观角度的力-热-化机理与模型建立及其数值模拟应用以及具备可调性能的活性材料制备新工艺3方面研究内容将是推进活性材料未来军事化应用的重点关注对象。
- 活性材料 /
- 冲击释能 /
- 冲击诱发化学反应 /
- 力-热-化耦合行为
Abstract: Reactive material is a new type of material with energy-releasing characteristics. It can react chemically and release a large amount of chemical energy under the high pressure and high temperature caused by the impact. Therefore, it has a wide range of potential applications in military fields such as fragments and energy-splitting warheads. In order to realize the design and control of the energy release process of reactive material and promote its weaponized application process, it is necessary to solve a series of complicated mechanical-thermal-chemical coupling problems in the impact-induced energy release behaviors of reactive materials. In the past 40 years, domestic and foreign scholars have carried out a lot of research on the impact-induced energy release behavior of reactive materials. Based on this, this paper systematically combs the research status of the impact-induced chemical reaction mechanism, kinetics and related effects of reactive materials, focusing on the research progress in three aspects: the experimental characterization technology of impact-induced energy release of materials, the theoretical model of impact-induced chemical reaction and the numerical simulation method of shock compression considering the mechanical-thermal-chemical coupling effects. Finally, the summaries are carried out and the future research work, challenges and suggestions are proposed. It is concluded that domestic and foreign scholars have accumulated a certain amount of research on the energy release behaviors of reactive materials, but there is still a lack of richer, finer and intuitive characterization and exploration for the real-time diagnosis of ultra-fast chemical reaction behavior in experiments. However, for the related theoretical and numerical simulation studies, the mechanical thermal chemical theoretical model which can fully describe the impact energy release behavior of active materials has not been established, and there is no effective method to describe the impact energy release behaviors from the macro scale. Therefore, the three aspects of research content, ultra-fast chemical reaction experimental characterization technology, macro-level mechanical-thermal-chemical mechanism and model establishment and its numerical simulation application, and preparation of reactive materials with adjustable properties, will be the focus of attention in promoting the future military application of reactive materials.
- reactive materials /
- shock-induced energy release /
- impact-induced chemical reaction /
- mechanical-thermal-chemical coupling behavior

HTML全文

含铝炸药在军事领域的广泛应用以及其典型的非理想特性，吸引了大量学者对含铝炸药的爆轰驱动特性和铝粉在爆轰产物中的反应过程开展研究。Miller等^[1-3]提出了针对高能非理想炸药和推进剂的反应流模型，该模型将含铝炸药爆轰分为两阶段：理想组分快速爆轰阶段和金属粉末缓慢燃烧阶段，并将该反应流模型应用于二维有限元Lagrangian流体力学代码计算。Frost等^[4-5]、Zhang等^[6]、Ripley等^[7]、Milne等^[8]、Cooper等^[9]、Massoni等^[10]和Kim等^[11]应用多相流体动力学模型（computational fluid dynamics，CFD）数值计算了铝粉在炸药中的反应条件。Kim等^[12]对含铝炸药反应流进行了实验研究，在JWL++状态方程的基础上，研究了RDX/Al/HTPB（50%/35%/15%）的尺寸效应和爆轰特性。

目前，含铝炸药的研究以实验为主要手段，对含铝炸药的爆轰驱动过程一直缺乏理论指导。本文中针对含铝炸药爆轰产物的非理想膨胀过程，提出局部等熵假设，并建立含铝炸药爆轰产物流动的非线性特征线理论模型，为分析含铝炸药的驱动做功能力以及爆轰产物的非等熵流动规律提供一种全新的方法，为深入理解含铝炸药爆轰产物的膨胀过程提供理论方法。

1. 含铝炸药爆轰产物非线性特征线模型的必要假设

1.1 含铝炸药爆轰反应机理假设

含铝炸药的爆轰机理比较复杂，基于目前的研究普遍认为：微米级以及粒度更大的铝粉在爆轰反应区基本不参加反应，铝粉反应主要发生在爆轰产物膨胀驱动过程中^[13]。裴红波等^[14]研究了粒径为9.8 μm的RDX/Al在不同铝粉含铝条件下的反应进程，结果显示，10 μs后铝粉的反应度才开始随时间逐渐增大。陈朗等^[15]对铝粉粒径为几十纳米到几十微米、铝粉含铝为20%的RDX/Al炸药进行了爆轰驱动研究，发现微米级铝粉的反应主要发生在爆轰反应后期。Kim等^[12]设计并研究了混合炸药RDX/Al(65/35)的柱型装药实验（铝粒子直径为微米级），发现大部分铝粒子不参与爆轰反应，反应主要发生于爆轰波阵面后的爆轰产物流动区。

根据文献调研分析以及二次反应理论对含铝炸药爆轰机理的论述，认为对含微米级或微米级以上铝粉的含铝炸药，铝粉在爆轰反应区不发生化学反应或只有极少量铝粉发生反应，因此假设铝粉在炸药爆轰反应区完全不发生化学反应，即铝粉在炸药爆轰反应区表现为惰性。纳米级铝粉活性较高，对于炸药爆轰反应区的影响不可忽略，因此，本文中建立的模型只适用于含有微米级铝粉的含铝炸药。

1.2 铝粉在炸药爆轰产物中的状态假设

（1）含铝炸药爆轰波阵面后的铝粉与爆轰产物以相同的速度运动。

Zhang等^[16]假设在冲击波通过铝粒子的过程中，炸药不发生反应，并采用Murnaghan状态方程描述炸药受到冲击的状态，计算结果发现铝粉与爆轰产物的相对速度与炸药密度有关，对于常用的RDX基和HMX基含铝炸药而言，可近似认为铝粉与爆轰产物一起运动。

（2）对于爆轰波阵面后的微米级铝粉，可认为铝粉以气相形式燃烧^[17-18]。

1.3 含铝炸药爆轰产物流动的局部等熵假设

铝粉反应对爆轰产物的热力学影响是不可忽略的，铝粉释放能量增加爆轰产物的内能，进而影响爆轰产物的压力、声速等状态参数，因此，经典理想炸药爆轰产物膨胀的等熵流动理论无法科学描述含铝炸药爆轰产物的流动规律。考虑到爆轰产物中铝粉的氧化反应速率相对较慢，一般可持续反应1～3 ms，若将含铝炸药爆轰产物的膨胀过程沿时间轴划分为许多微小的时间域，在每个微小时间域内，爆轰产物中铝粉反应量较少且爆轰产物体积膨胀较小，那么在每个微小时间域内可近似认为爆轰产物膨胀过程是等熵的，据此本文提出了分析含铝炸药爆轰产物流动的局部等熵假设，其主要假设包括以下内容：

（1）假设铝粉在爆轰反应区内不发生化学反应，即铝粉在炸药爆轰反应区表现为惰性；

（2）考虑到铝粉反应速率相对较慢，假设铝粉反应对爆轰产物的影响具有一定的弛豫效应，也就是说铝粉反应释放的能量需要一段时间后才会对爆轰产物参数（包括当地声速、压强和产物密度）产生影响；

（3）将含铝炸药爆轰产物的膨胀过程沿时间轴分割为有限微小时间域（本文中的计算以0.2 μs为一个微时间域），在每个微时间域内，铝粉反应释放的能量没有立刻对周围爆轰产物状态参数产生影响，因此，可认为任意微时间域内铝粉的反应度不变，爆轰产物可近似为等熵膨胀；

（4）铝粉反应对爆轰产物状态参数的影响体现在相邻下一微时间域的初始时刻，即相邻两个微时间域的交界处，爆轰产物的压强、密度和声速将产生变化，交界处产物粒子速度不变（内能的微小变化不影响粒子速度），在任一微时间域内，爆轰产物遵循等熵流动规律，但微时间域间的熵不同。

以爆轰产物任意一条特征线为例，说明铝粉反应释放能量对爆轰产物流动规律的影响，特征线示意图如图1所示。

图 1 铝粉反应对特征线的影响示意图

Figure 1. Schematic diagram of the influence of aluminum powder reaction on the characteristic lines

下载: 全尺寸图片幻灯片

当铝粉在爆轰产物中发生反应时，由于受到铝粉反应放热的影响，沿特征线的爆轰产物状态参数将不再是常数，爆轰产物的熵也不再是固定值。在微时间域1内（0～t₁），铝粉没有开始反应或反应还没有对爆轰产物参数产生影响，在绝热膨胀条件下，此时间域内爆轰产物的熵S₁恒定不变，特征线为直线，沿特征线的爆轰产物状态参数压强p₁、密度ρ₁、声速c₁等同样恒定不变（经典等熵膨胀理论），也就是说铝粉不反应或不作用时，非等熵模型与等熵模型是等价的；受铝粉反应释放能量的影响，在t₁时刻（两个微时间域的交界时刻），产物状态参数变为p₂、ρ₂、c₂，依据局部等熵假设，在微时间域2内（t₁～t₂），产物状态参数变为p₂、ρ₂、c₂沿特征线维持不变，但受到铝粉反应影响，熵值变为S₂且特征线斜率也发生变化，表明爆轰产物流动规律发生变化。依据局部等熵假设，可分析出其他微时间域的特征线变化。

含铝炸药爆轰产物流动的局部等熵假设是在经典等熵特征线理论的基础上提出的，为了能够科学地理论分析铝粉二次反应对爆轰产物流动的影响规律，将炸药置于无限长刚性圆管中，忽略了侧向稀疏作用的影响，且忽略了炸药对外部环境的热效应（炸药爆轰驱动过程极其迅速，在有效驱动时间内对外部环境的热效应很小）。因此，含铝炸药非线性特性线模型适用于强约束条件下，混合微米级铝粉的含铝炸药爆轰产物的流动规律分析及金属板有效驱动段（金属板从开始运动后的6～8 μs达到最大速度，这段驱动时间称为有效驱动段）的运动规律分析。

2. 含铝炸药爆轰驱动的非线性特征线模型

2.1 含铝炸药爆轰产物状态方程

假设含铝炸药处于无限长刚性圆管中，炸药两侧为真空，因此爆轰产物流动可视为一维流动。上一节必要假设中提到在C-J面后铝粉与爆轰产物以相同的速度运动，忽略了铝粉与理想组分的相互作用，因此，含铝炸药爆轰产物中的理想组分的状态方程可以表示为：

${p_{\rm{e}} } = {w_{\text{b}}}\rho RT$

(1)

式中： ${p_{\text{e}}}$ 为含铝炸药爆轰产物压力， ${w_{\text{b}}}$ 为含铝炸药中炸药组分的初始质量分数， $\rho$ 为含铝炸药爆轰产物密度， $R$ 为气体常数， $T$ 为含铝炸药爆轰产物温度。

铝粉的燃烧状态近似表现为气相燃烧，且铝粉均匀分布于爆轰产物中，因此，铝粉燃烧对压力的贡献可表示为^[12]：

${p_{{\text{Al}}}} = \frac{{{w_{\text{a}}}{\lambda _{{\text{Al}}}}\rho RT}}{{1 - n{A_{\text{n}}}}}$

(2)

式中： ${p_{{\text{Al}}}}$ 为铝粉燃烧对压力的贡献， ${\lambda _{{\text{Al}}}}$ 为已经反应了的铝粉的质量分数， ${w_{\text{a}}}$ 为含铝炸药中铝粉的初始质量分数， $n$ 为单位体积的物质的量， ${A_{\text{n}}}$ 为经验系数且是一个常数。本文计算中 $n{A_{\text{n}}}$ 取值0.75。

应用标准混合准则，含铝炸药爆轰产物的状态方程可表示为：

${p_{{\text{total}}}} = {p_{\rm{e}} } + {p_{{\text{Al}}}} = \left({w_{\text{b}}} + {w_{\text{a}}}\frac{{{\lambda _{{\text{Al}}}}}}{{1 - n{A_{\text{n}}}}}\right)\rho RT{\text{ = }}A\left( {{\lambda _{{\text{Al}}}}} \right)\rho RT$

(3)

式中： $A\left( {{\lambda _{{\text{Al}}}}} \right){\text{ = }}\left({w_{\text{b}}} + {w_{\text{a}}}\dfrac{{{\lambda _{{\text{Al}}}}}}{{1 - n{A_{\text{n}}}}}\right)$ 。

通过热力学变换，可以得到微时间域 $i$ 内含铝炸药爆轰产物的等熵方程：

$\frac{{{p_{{\text{total}}}}}}{{{\rho ^{{\gamma _i}}}}} = {\left[ {\left({w_{\text{b}}} + {w_{\text{a}}}\frac{{{\lambda _{{\text{Al}}i}}}}{{1 - n{A_{\text{n}}}}}\right) \cdot {C_i}R} \right]^{{\gamma _i}}}$

(4)

式中： ${C_i}$ 为常数， ${\lambda _{{\text{Al}}i}}$ 为时间域 $i$ 内的铝粉反应度，下标 $i$ 表示第 $i$ 时间域。

2.2 含铝炸药爆轰产物的非线性特征线模型

忽略爆轰产物流动过程中的热传递，因此，爆轰产物平面绝热运动的流体动力学方程组为：

$\left\{ \begin{array}{l} \dfrac{{\partial \rho }}{{\partial t}} + u\dfrac{{\partial \rho }}{{\partial x}} + \rho \dfrac{{\partial u}}{{\partial x}} = 0 \\ \dfrac{{\partial u}}{{\partial t}} + u\dfrac{{\partial u}}{{\partial x}} + \dfrac{1}{\rho }\dfrac{{\partial p}}{{\partial x}} = 0 \\ \dfrac{{\partial S}}{{\partial t}} + u\dfrac{{\partial S}}{{\partial x}} = 0\end{array} \right.$

(5)

对流场状态参数全微分，得到：

${\text{d}}p = {\left(\frac{{\partial p}}{{\partial \rho }}\right)_S}{\text{d}}\rho + {\left(\frac{{\partial p}}{{\partial S}}\right)_\rho }{\text{d}}S$

(6)

经变换并联立方程组(5)第1式，得到：

$\frac{{{\text{d}}p}}{{{\text{d}}t}} - {\left(\frac{{\partial p}}{{\partial E}}\right)_\rho }\left(\frac{{T{\text{d}}S}}{{{\text{d}}t}}\right) + \rho {c^2}\frac{{\partial u}}{{\partial x}} = 0$

(7)

引入系数 $\lambda$ ，联立方程组(5)第2式和式(7)，得到：

$\lambda \left[\frac{{\partial u}}{{\partial t}} + \left(u + \frac{{\rho {c^2}}}{\lambda }\right)\frac{{\partial u}}{{\partial x}}\right] + \left[\frac{{\partial p}}{{\partial t}} + \left(u + \frac{\lambda }{\rho }\right)\frac{{\partial p}}{{\partial x}}\right] - {\left(\frac{{\partial p}}{{\partial E}}\right)_\rho }\left(\frac{{T{\text{d}}S}}{{{\text{d}}t}}\right) = 0$

(8)

令 $\lambda = \pm \rho c$ ，得到：

$\pm \rho c\left[\frac{{\partial u}}{{\partial t}} + (u \pm c)\frac{{\partial u}}{{\partial x}}\right] + \left[\frac{{\partial p}}{{\partial t}} + (u \pm c)\frac{{\partial p}}{{\partial x}}\right] - {\left(\frac{{\partial p}}{{\partial E}}\right)_\rho }\left(\frac{{T{\text{d}}S}}{{{\text{d}}t}}\right) = 0$

(9)

当 $\dfrac{{{\text{d}}x}}{{{\text{d}}t}} = u \pm c$ 时，方程(9)可变为：

$\left\{ \begin{array}{l} \dfrac{{{\text{d}}x}}{{{\text{d}}t}} = u \pm c \\ \dfrac{{{\text{d}}p}}{{{\text{d}}t}} \pm \rho c\left(\dfrac{{{\text{d}}u}}{{{\text{d}}t}}\right) - {\left(\dfrac{{\partial p}}{{\partial E}}\right)_\rho }\left(\dfrac{{T{\text{d}}S}}{{{\text{d}}t}}\right) = 0 \end{array} \right.$

(10)

根据含铝炸药爆轰产物的局部等熵假设，在每个小时间域内，爆轰产物的流动规律可近似认为是等熵流动。因此，对于微时间域 $i$ ，方程组(10)可表示为：

$\left\{ \begin{array}{l} \dfrac{{{\text{d}}x}}{{{\text{d}}t}} = {u_i} \pm {c_i} \\ {\text{d}}{u_i} \pm \dfrac{1}{{{\rho _i}{c_i}}}{\text{d}}{p_i}{\text{ = }}0 \end{array} \right.$

(11)

对方程组(11)第2式求积分，得到：

${u_i} \pm \int {\frac{{{\text{d}}p}}{{\rho c}}} = {u_i} \pm \int {\frac{{c{\text{d}}\rho }}{\rho }}$

(12)

根据爆轰产物的等熵方程(4)，声速 $c_i^2 = {\gamma _i}{\left[ {\left({w_{\text{b}}} + {w_{\text{a}}}\dfrac{{{\lambda _{{\text{Al}}i}}}}{{1 - n{A_{\text{n}}}}}\right){C_i}R} \right]^{{\gamma _i}}}\rho _i^{{\gamma _i} - 1}$ ，于是可以得到：

$\int {\frac{{{c_i}{\text{d}}{\rho _i}}}{{{\rho _i}}}} = \int {\frac{2}{{{\gamma _i} - 1}}} {\text{d}}{c_i} = \frac{2}{{{\gamma _i} - 1}}{c_i}$

(13)

由此可以得到爆轰产物的特征线方程：

$\left\{ \begin{array}{l} \dfrac{{{\rm{d}}x}}{{{\rm{d}}t}} = {u_i} \pm {c_i}\\ {u_i} \pm \dfrac{2}{{{\gamma _i} - 1}}{c_i} = {\rm{constant}} \end{array} \right.$

(14)

式中： ${u_i}$ 和 ${c_i}$ 分别表示微时间域 $i$ 内爆轰产物的粒子速度和当地声速，由于在微时间域 $i$ 内产物流动是局部等熵的，因此，沿特征线 ${u_i}$ 和 ${c_i}$ 的值保持不变。

2.3 含铝炸药爆轰驱动金属板的非线性特征线计算

假设长为L的含铝炸药和质量为M的金属平板置于无限长刚性圆管中，炸药两侧为真空环境。采用一侧端面引爆炸药，炸药对金属平板驱动的特征线图如图2所示。

图 2 含铝炸药爆轰产物驱动金属板的特征线示意图

Figure 2. Schematic diagram of the characteristic lines of the aluminized explosive

下载: 全尺寸图片幻灯片

设含铝炸药引爆时刻为0时刻，金属板后爆轰产物的流动方程可表示为：

$\left\{ \begin{array}{l} x = ({u_i} + {c_i})t + {F_i} \\ x = ({u_i} - {c_i})t + F_i^{'} \end{array} \right.$

(15)

式中： ${F_i}$ 为与铝粉反应度相关的特征线参数， $F_i^{'}$ 为与铝粉反应度和金属板运动相关的特征线参数。

根据经典力学牛顿第二定律可知：

$M\frac{{{\text{d}}v}}{{{\text{d}}t}} = {A_{\text{r}}}{p_{\text{m}}}$

(16)

式中： $M$ 为金属板的质量， $v$ 为金属板的速度， ${A_{\text{r}}}$ 为金属板的横截面面积， ${p_{\text{m}}}$ 为金属板内表面处爆轰产物的压强。基于局部等熵假设，在微时间域 $i$ 内，可以得到如下关系：

$\frac{{p{}_{\text{m}}}}{{{p_i}}} = {\left( {\frac{{{c_{\text{m}}}}}{{{c_i}}}} \right)^\gamma }$

(17)

式中： ${p_{\text{m}}}$ 和 ${c_{\text{m}}}$ 分别为紧挨金属板内表面处爆轰产物的压强和当地声速， ${p_i}$ 为微时间域 $i$ 内右传波到达金属板内表面前爆轰产物的初始压强。将式(17)代入式(16)，得到：

$\frac{{{\text{d}}v}}{{{\text{d}}t}} = \frac{{{A_{\text{r}}}{p_i}}}{{Mc_i^\gamma }}c_{\text{m}}^\gamma$

(18)

爆轰产物中的每一道右传波都以各自的 $u + c$ 速度传播，并且在微时间域内沿特征线传播速度保持不变。当右传波追赶上运动的金属板时，右传波将发生反射，在此瞬间，爆轰产物的粒子速度立即由 ${u_i}$ 降低为金属板壁面的速度 ${u_{\text{m}}}$ ，声速也立即由 ${c_i}$ 变为金属板后产物的声速 ${c_{\text{m}}}$ ，从而得到：

${u_i} + {c_i} = u{}_{\text{m}} + {c_{\text{m}}}$

(19)

对方程组(15)中的第1式对时间 $t$ 求导得到：

$\frac{{{\text{d}}x}}{{{\text{d}}t}} = \left( {{u_{\text{m}}} + {c_{\text{m}}}} \right) + \left( {\frac{{{\text{d}}{u_{\text{m}}}}}{{{\text{d}}t}} + \frac{{{\text{d}}{c_{\text{m}}}}}{{{\text{d}}t}}} \right)t$

(20)

由于金属板内表面处的产物速度 ${u_{\text{m}}}$ 与金属板运动的速度 $v$ 相等，并且 $\dfrac{{{\rm{d}}x}}{{{\rm{d}}t}} = v = {u_{\rm{m}}}$ ，因此将式(18)代入式(20)，得到：

$\frac{{{\text{d}}{c_{\text{m}}}}}{{{\text{d}}t}} + \frac{{{c_{\text{m}}}}}{t} + \psi c_{\text{m}}^\gamma = 0$

(21)

式中： $\psi = \dfrac{{{A_{\text{r}}}{p_i}}}{{Mc_i^\gamma }}$ 。考虑到炸药处于无限长刚性圆管中，可认为爆轰产物在短时间内始终处于高压状态，因此有 $\gamma \approx 3$ 。应用第2节中提到的参数变换法求解微分方程(21)，得到：

$\frac{{1 + 2\psi c_{\text{m}}^2t}}{{c_{\text{m}}^2{t^2}}} = {\vartheta _1}{\text{ = constant}}$

(22)

在微时间域1，即 $0\sim{t}_{1}$ 时间段内，金属板后的爆轰产物的初始状态参数为：

$\left\{ \begin{array}{l} {p_1} = {p_{\text{H}}} = \dfrac{1}{4}{\rho _0}D_{{\text{Al}}}^2 \\ {c_1} = {c_{\text{H}}} = \dfrac{3}{4}{D_{{\text{Al}}}} \end{array} \right.$

(23)

其边界条件为：

$\left\{ \begin{array}{l} t = \dfrac{L}{{{D_{{\text{Al}}}}}} \\ {u_{\text{m}}} = 0 \\ {c_{\text{m}}} = {D_{{\text{Al}}}} \end{array} \right.$

(24)

式中： ${\rho _0}$ 为含铝炸药的初始密度， ${p_{\text{H}}}$ 和 ${c_{\text{H}}}$ 为爆轰波阵面上产物的压力和声速， $L$ 为炸药的长度， ${D_{{\text{Al}}}}$ 为含铝炸药的爆速。将边界条件(24)代入式(22)，得到：

${\vartheta _1} = \frac{{1 + \eta }}{{{L^2}}}{\text{ = constant}}$

(25)

式中： $\eta = \dfrac{{32}}{{27}}\dfrac{m}{M}$ ，其中 $m$ 为炸药质量， $m = {A_{\text{r}}}L{\rho _0}$ 。

将式(25)代入式(22)，得到：

${c_{\text{m}}} = \frac{L}{t}{\xi _1}$

(26)

式中： ${\xi _1} = {\left[ {1 + \eta \left( {1 - \dfrac{L}{{{D_{{\text{Al}}}}t}}} \right)} \right]^{ - 0.5}}$ 。

将式(26)代入式(18)，并积分得到：

$v = {D_{{\text{Al}}}}\left( {1 + \frac{{2\left( {{\xi _1} - 1} \right)}}{{\eta {\xi _1}}} - \frac{{L{\xi _1}}}{{{D_{{\text{Al}}}}t}}} \right)$

(27)

式（27）表示在微时间域1内含铝炸药对质量为M的金属板的驱动规律。

基于局部等熵假设，结合以上方法，可以得到任意微时间域 $i$ 内，爆轰产物对金属板的驱动速度如下：

$\left\{ \begin{array}{l} v - v\left( {{t_{i - 1}}} \right) = \left( {\dfrac{{L{\xi _i}\left( {{t_{i - 1}}} \right)}}{{{t_{i - 1}}}} + \dfrac{{2{D_{{\text{Al}}}}}}{\eta }\dfrac{1}{{{\xi _i}\left( {{t_{i - 1}}} \right)}} - \dfrac{{L{\xi _i}}}{t} - \dfrac{{2{D_{{\text{Al}}}}}}{{\eta {\xi _i}}}} \right) \\ {\xi _i} = {\left( {{L^2}{\vartheta _i} - \dfrac{{\eta L}}{{{D_{{\text{Al}}}}t}}} \right)^{ - 0.5}} \\ {\vartheta _i} = \dfrac{{1 + \dfrac{{\eta c_{\text{m}}^2({t_{i - 1}}){t_{i -1}}}}{{L{D_{{\text{Al}}}}}}}}{{c_{\text{m}}^2({t_{i - 1}})t_{i - 1}^2}} = {\text{constan}}{{\text{t}}} \end{array} \right.$

(28)

式中： ${t_{i{\text{ - }}1}}$ 表示微时间域 $i$ 的起始时刻， $v\left( {{t_{i - 1}}} \right)$ 表示微时间域 $i$ 的起始时刻金属板的速度。

3. 含铝炸药爆轰产物非线性特征线模型的实验验证

3.1 含铝炸药驱动金属板的实验结果

选择RDX基含铝炸药作为驱动炸药，炸药中铝粉的平均直径为5、50 μm两种，分别驱动厚0.5、1 mm的紫铜板。炸药参数如表1所示。

表 1 实验所用的炸药参数和铜板尺寸

Table 1. Parameters of the explosive and size of the copper plate

编号	炸药	铝粉（或LiF）质量分数/%	炸药密度/（g·cm⁻³）	铝粉直径/μm	金属板尺寸/mm×mm
1	RDX/Al/黏结剂	20	1.82	5	$\varnothing$ 50×1
2	RDX/Al/黏结剂	20	1.82	50	$\varnothing$ 50×1
3	RDX/LiF/黏结剂	20	1.80	−	$\varnothing$ 50×1
4	RDX/Al/黏结剂	20	1.82	5	$\varnothing$ 50×0.5
5	RDX/Al/黏结剂	20	1.82	50	$\varnothing$ 50×0.5
6	RDX/LiF/黏结剂	20	1.80	−	$\varnothing$ 50×0.5
注：炸药尺寸为 $\varnothing$ 50 mm×50 mm，黏结剂含量占炸药总量的5%。

下载: 导出CSV

| 显示表格

实验布置示意图如图3所示。

图 3 炸药驱动金属板实验布置示意图

Figure 3. Schematic diagram of the test arrangement of the explosive driven metal plate

下载: 全尺寸图片幻灯片

含铝炸药和含LiF炸药驱动金属板实验结果如图4～5所示。

图 4 含铝炸药和含LiF炸药驱动0.5 mm厚铜板的测试结果

Figure 4. The velocity of a 0.5 mm thick copper plate driven by the aluminized explosive and LiF explosive

下载: 全尺寸图片幻灯片

图 5 含铝炸药和含LiF炸药驱动1 mm厚铜板的测试结果

Figure 5. The velocity of a 1 mm thick copper plate driven by the aluminized explosive and LiF explosive

下载: 全尺寸图片幻灯片

3.2 实验条件下铝粉的反应度变化规律

在同等实验条件下，以含铝炸药驱动金属板的动能减去含LiF炸药驱动金属板的动能，就得到了铝粉反应释放能量对金属板所做的有用功：

$\Delta E\left( t \right) = \frac{1}{2}M\left[ {v_{{\text{Al}}}^2\left( t \right) - v_{{\text{LiF}}}^2\left( t \right)} \right]$

(29)

式中： ${v_{{\text{Al}}}}$ 和 ${v_{{\text{LiF}}}}$ 分别表示含铝炸药和含LiF炸药驱动金属板的运动速度。

设平板驱动实验的效率为 $\varpi$ ，铝粉的氧化反应热为 ${Q_{{\text{Al}}}}$ ，铝粉的反应度为 $\lambda \left( t \right)$ ，据此可得到：

$\Delta E\left( t \right){\text{ = }}\varpi {Q_{{\text{Al}}}}m{w_{\text{a}}}\lambda \left( t \right)$

(30)

式中： $m$ 表示含铝炸药的质量， ${w_{\text{a}}}$ 表示含铝炸药中铝粉的质量分数。

将式(29)代入式(30)，得到：

$\varpi {Q_{{\text{Al}}}}m{w_{\text{a}}}\lambda \left(t \right){\text{ = }}\frac{1}{2}M\left[ {v_{{\text{Al}}}^2\left( t \right) - v_{{\text{LiF}}}^2\left( t \right)} \right]$

(31)

通过计算化学生成焓得到 ${Q_{{\text{Al}}}}$ =20.126 kJ/g，炸药驱动金属板的实验效率 $\varpi$ 取为0.18^[19]，据此可计算出铝粉的反应度变化，计算得到的反应度变化规律如图6～7所示：

图 6 含铝炸药驱动1 mm厚金属的铝粉反应度

Figure 6. Reaction degree of the detonation products of aluminium powder behind the 1 mm thick copper plate

下载: 全尺寸图片幻灯片

图 7 50 μm含铝炸药驱动0.5 mm金属板的铝粉反应度

Figure 7. Reaction degree of the detonation products of aluminium powder behind the 0.5 mm thick copper plate

下载: 全尺寸图片幻灯片

3.3 非线性特征线理论计算结果与实验结果对比

基于含铝炸药爆轰产物的非等熵流动模型，应用非线性特征线计算方法计算5、50 μm含铝炸药对0.5、1 mm金属板的驱动过程，含铝炸药参数见表2。

表 2 含铝炸药和含LiF炸药的参数

Table 2. Parameters of the aluminized explosives and the LIF explosive

炸药组成	铝粉质量分数/%	铝粉直径/μm	炸药密度/（g·cm⁻³）	炸药爆速/（km·s⁻¹）
RDX/Al	20	5	1.8	8.223
RDX/Al	20	50	1.8	8.223
RDX/LiF	−	−	1.8	8.223

下载: 导出CSV

| 显示表格

含铝炸药和含LiF炸药驱动0.5、1 mm厚金属板的实验和计算结果对比如图8～9所示。

图 8 50 μm含铝炸药和含LiF炸药驱动0.5 mm金属的实验和计算结果

Figure 8. The test and calculated velocities of a metal plate driven by the aluminized explosive and LiF explosive with a particle diameter of 50 μm

下载: 全尺寸图片幻灯片

图 9 不同粒度含铝炸药和含LiF炸药驱动1 mm金属板的实验和计算结果

Figure 9. The test and calculated velocities of a 1 mm metal plate driven by the aluminized explosives and LiF explosives with different particle diameters

下载: 全尺寸图片幻灯片

从图8中的对比结果可以看出，非线性特征线模型的计算结果与实验结果保持了较好的一致性。在6.3～8 μs的驱动时间段内，计算结果稍大于实验测试结果，8 μs后计算结果与实验结果一致性较高。需要指出的是，此结果仅针对爆轰波到达金属后的6 μs的时间内，由于强约束壳体膨胀破裂的时间约为6 μs，在此时间段内爆轰产物受稀疏波的影响较小，因此理论结果与实验结果较一致，而当强约束壳体破裂时，产物开始泄露稀疏波将对驱动将产生很大影响，计算结果将高于实验测试结果。

图9为5 μm含铝炸药、50 μm含铝炸药和含LiF炸药对1 mm厚金属板的实验和计算结果。计算结果较好地体现了铝粉反应对炸药做功能力的贡献以及铝粉粒度不同对炸药做功能力的影响，但由于受到1 mm金属板内应力波的影响，驱动加速段实验测试结果波动较大。但总体而言，非线性特征线模型能够较好地体现铝粉后期反应对做功能力的贡献，且前期驱动结果与强约束条件下的驱动金属板实验结果保持了较好的一致性。

4. 结　论

主要对微米级铝粉含铝炸药爆轰产物的非等熵膨胀规律进行了理论分析，目标是建立一种理论分析模型，描述包含铝粉氧化反应的爆轰产物流动规律，理论研究铝粉反应对炸药做功能力的贡献，为理论分析含铝炸药爆轰产物的流动规律提供一种方法。

针对微米级铝粉含铝炸药爆轰产物的非理想特性，在经典理想炸药爆轰产物等熵流动理论的基础上，通过科学分析和计算，得到如下结论：

（1）根据微米级铝粉含铝炸药的二次反应特性，提出含铝炸药爆轰产物膨胀过程的局部等熵假设，应用特征线方法能够科学构建含铝炸药爆轰产物的准解析模型；

（2）基于理想炸药爆轰产物的近似状态方程和铝粉反应对压力的贡献方程，能够得到含铝炸药爆轰产物状态方程；

（3）基于局部等熵假设和含铝炸药爆轰产物状态方程，应用非线性特征线方法计算了含铝炸药爆轰产物的流动规律，并推导出含铝炸药爆轰驱动的非线性特征线方程组，分析了含铝炸药爆轰驱动金属板和金属板后产物的非等熵膨胀过程，非线性特征线模型能够定性分析微米级铝粉在爆轰产物膨胀区的二次反应特性；

（4）设计了5、50 μm含铝炸药和含LiF炸药在强约束条件下驱动0.5、1 mm厚金属板实验，通过激光位移干涉仪测试金属板运动的速度历程，通过实验结果间接计算了铝粉在爆轰产物中的反应度，结合含铝炸药爆轰产物的非等熵流动模型，理论计算了不同炸药驱动金属板的速度历程，对比实验结果，理论方法能够很好的描述铝粉二次反应对炸药做功能力的贡献，且与实验结果能够保持较好的一致性，验证了非等熵流动模型对于微米级铝粉含铝炸药的正确性。

图 1 惰性/活性破片对航空煤油打击（上）与航空电子设备打击（下）效果对比^[14]

Figure 1. Comparison of the effect of inert/active chips against aviation kerosene (upper) and avionics (lower)^[14]

下载: 全尺寸图片幻灯片

图 2 平板撞击实验^[28]

Figure 2. Plate impact experiment^[28]

下载: 全尺寸图片幻灯片

图 3 Al/Ni材料的输入压力-冲击波速度曲线^[31]

Figure 3. Input stress-shock wave velocity curves for Al/Ni materials^[31]

下载: 全尺寸图片幻灯片

图 4 Ni/Al实验测量状态方程与模拟计算的状态方程对比^[35]

Figure 4. Ni/Al experimental measurement equation of state is compared with the simulation calculation equation of state^[35]

下载: 全尺寸图片幻灯片

图 5 直接弹道实验

Figure 5. Direct trajectory experiment

下载: 全尺寸图片幻灯片

图 6 Ti/Si粉末微观结构和XRD分析结果^[51]

Figure 6. Microstructure and XRD results of Ti/Si mixture^[51]

下载: 全尺寸图片幻灯片

图 7 直接弹道实验回收试样分析结果^[59]

Figure 7. Analysis results of recovered samples in direct trajectory experiment^[59]

下载: 全尺寸图片幻灯片

图 8 Graham^[61]提出的冲击诱发化学反应构图

Figure 8. Shock-induced chemical reaction model proposed by Graham^[61]

下载: 全尺寸图片幻灯片

图 9 反应动力学参数计算结果^[66]

Figure 9. Calculation results of kinetic parameter of reaction^[66]

下载: 全尺寸图片幻灯片

图 10 离散元模型的冲击反应算例^[63]

Figure 10. An example of impact response by the discrete element method model^[63]

下载: 全尺寸图片幻灯片

图 11 两种Ni/Al粉末的冲击波前沿压力面分布^[78]

Figure 11. Pressure surface profiles of two types of Ni/Al particles^[78]

下载: 全尺寸图片幻灯片

图 12 活性材料的多尺度模型冲击释能数值模拟计算流程^[84-85]

Figure 12. Flow of numerical simulation of RMs impact release energy with the multi-scale model^[84-85]

下载: 全尺寸图片幻灯片

图 13 中尺度计算模型模拟的Al/Ni粉末混合物在冲击压缩状态下的温度云图^[86]

Figure 13. Temperature cloud maps of Al/Ni powder mixture under impact compression simulated by the mesoscale calculation model^[86]

下载: 全尺寸图片幻灯片

图 14 分子动力学模拟计算结果：在Al/Ni颗粒界面接触处出现了金属相变与金属间化学反应^[87]

Figure 14. Molecular dynamics simulation results: metal phase transitions and intermetallic chemical reactionsoccurred at the interface of Al/Ni particles^[87]

下载: 全尺寸图片幻灯片

参考文献(87)

[1]	WILLIS M J, JASON T D. Effect of aluminum particle size on the impact initiation of pressed rods [J]. Shock Compression of Condensed Matter, 2007, CP955: 971–974.
[2]	HAMBLING D. Behind the reactive materials revolution [N/OL]. [2008–05–07] [2020–09–22] https://www.wired.com/2008/05/reactive-revolu-3/.
[3]	张先锋, 赵晓宁. 多功能含能结构材料研究进展 [J]. 含能材料, 2009, 17(6): 731–739. DOI: 10.3969/j.issn.1006-9941.2009.06.021. ZHANG X F, ZHAO X N. Review on multifunctional energetic structural materials [J]. Chinese Journal of Energetic Materials, 2009, 17(6): 731–739. DOI: 10.3969/j.issn.1006-9941.2009.06.021.
[4]	XIONG W, ZHANG X F, ZHENG L, et al. The shock-induced chemical reaction behaviour of Al/Ni composites by cold rolling and powder compaction [J]. Journal of Materials Science, 2019, 54(1): 6651–6667. DOI: 10.1007/s10853-019-03357-3.
[5]	CHENG J L, HNG H H, LEE Y W, et al. Kinetic study of thermal- and impact-initiated reactions in Al-Fe₂O₃ nanothermite [J]. Combustion and Flame, 2010, 157(12): 2241–2249. DOI: 10.1016/j.combustflame.2010.07.012.
[6]	PERRY W L, TAPPAN B C, REARDON B L, et al. Energy release characteristics of the nanoscale aluminum-tungsten oxide hydrate metastable intermolecular composite [J]. Journal of Applied Physics, 2007, 101(6): 064313. DOI: 10.1063/1.2435797.
[7]	WANG L, LIU J X, LI S K, et al. Insensitive high-energy energetic structural material of tungsten-polytetrafluoroethylene-aluminum composites [J]. AIP Advances, 2015, 5(11): 117142. DOI: 10.1063/1.4936557.
[8]	FENG B, FANG X, LI Y C, et al. Influence of processing techniques on mechanical properties and impact initiation of an Al-PTFE reactive material [J]. Central European Journal of Energetic Materials, 2016, 13(4): 989–1004. DOI: 10.22211/cejem/61496.
[9]	HERBOLD E B, JORDAN J L, THADHANI N N. Effects of processing and powder size on microstructure and reactivity in arrested reactive milled Al+Ni [J]. Acta Materialia, 2011, 59(17): 6717–6728. DOI: 10.1016/j.actamat.2011.07.029.
[10]	XIONG W, ZHANG X F, WU Y, et al. Influence of additives on microstructures, mechanical properties and shock-induced reaction characteristics of Al/Ni composites [J]. Journal of Alloys and Compounds, 2015, 648: 540–549. DOI: 10.1016/j.jallcom.2015.07.004.
[11]	WEI H Y, YOO C S. Dynamic responses of reactive metallic structures under thermal and mechanical ignitions [J]. Journal of Materials Research, 2012, 27(21): 2705–2717. DOI: 10.1557/jmr.2012.302.
[12]	刘晓俊, 任会兰, 宁建国. Zr-W多功能含能结构材料的制备及动态压缩特性 [J]. 复合材料学报, 2016, 33(10): 2297–2303. DOI: 10.13801/j.cnki.fhclxb.20160315.001. LIU X J, REN H L, NING J G. Preparation and dynamic compression properties of Zr-W multifunctional energetic structural material [J]. Journal of Composite Materials, 2016, 33(10): 2297–2303. DOI: 10.13801/j.cnki.fhclxb.20160315.001.
[13]	REN K R, CHEN R. LIN Y L, et al Probing the impact energy release behavior of Al/Ni-based reactive metals with experimental and numerical methods [J]. Metals, 2019, 9(5): 499. DOI: 10.3390/met9050499.
[14]	Committee on Advanced Energetic Materials and Manufacturing Technologies, National Research Council. Advanced energetic materials [M]. Washington D C: The National Academies Press, 2004. DOI: 10.17226/10918.
[15]	丁亮亮. PELE弹活性内芯配方与弹体结构设计及毁伤机理研究 [D]. 长沙: 国防科技大学, 2019.
[16]	WADDELL J T, BOOTES T H, BUDY G D, et al. Reactive shaped charge, reactive liner, and method for target penetration using a reactive shaped charge: US8037829B1 [P]. 2011-08-18.
[17]	BATSANOV S S, DORONIN G S, KLOCHKOV S V, et al. Synthesis reactions behind shock fronts [J]. Combustion, Explosion, and Shock Waves, 1986, 22(6): 765–768. DOI: 10.1007/BF00751890.
[18]	BAUER F. Ferroelectric properties and shock response of a poled PVF₂ polymer and of VF₂/C₂F₃H copolymers [C] // Shock Waves in Condensed Matter. Boston, MA: Springer, 1986: 483–496. DOI: 10.1007/978-1-4613-2207-8_69.
[19]	GRAHAM R A, ANDERSON M U, BAUER F, et al. Piezoelectric polarization of the ferroelectric polymer PVDF from 10 MPato 10 GPa: studies of loading-path dependence [C]// SCHMIDT S C, DICK R D, FORBES J W, et al. Shock Compression of Condensed Matter. Amsterdam: Elsevier, 1992: 883–886. DOI: 10.1016/B978-0-444-89732-9.50203-X.
[20]	NELLIS W J, MITCHELL A C, MCCANDLESS P C, et al. Electronic energy gap of molecular hydrogen from electrical conductivity measurements at high shock pressures [J]. Physical Review Letters, 1992, 68(19): 2937–2940. DOI: 10.1103/PhysRevLett.68.2937.
[21]	BOSLOUGH M B. Shock-induced chemical reactions in nickel-aluminum powder mixtures: radiation pyrometer measurements [J]. Chemical Physics Letters, 1989, 160(5–6): 618–622. DOI: 10.1016/0009-2614(89)80074-0.
[22]	BARKER L M, HOLLENBACH R E. Laser interferometer for measuring high velocities of any reflecting surface [J]. Journal of Applied Physics, 1972, 43(11): 4669–4675. DOI: 10.1063/1.1660986.
[23]	BARKER L M, HOLLENBACH R E. Interferometer technique for measuring the dynamic mechanical properties of materials [J]. Review of Scientific Instruments, 1965, 36(11): 1617–1620. DOI: 10.1063/1.1719405.
[24]	MOORE D S, SCHMIDT S C. Experimental molecular spectroscopy in shock-compressed materials [C] // Proceedings the 5th APS Topical Conference on ‘Shock Waves in Condensed Matter’. Monterey, CA, USA: APS, 1987.
[25]	GUPTA Y M. Time-resolved optical spectroscopy under shock loading: electronic and chemical changes in liquid carbon disulfide [J]. High Pressure Research, 1992, 10(5–6): 717–732. DOI: 10.1080/08957959208225324.
[26]	DLOTT D D. Ultrafast spectroscopy of shock waves in molecular materials [J]. Annual Review of Physical Chemistry, 1999, 50(1): 251–278. DOI: 10.1146/annurev.physchem.50.1.251.
[27]	TROTT W M, KNUDSON M D, CHHABILDAS L C, et al. Measurements of spatially resolved velocity variations in shock compressed heterogeneous materials using a line-imaging velocity interferometer [J]. AIP Conference Proceedings, 2000, 505(1): 993–998. DOI: 10.1063/1.1303635.
[28]	EAKINS D E, THADHANI N N. The shock-densification behavior of three distinct Ni+Al powder mixtures [J]. Applied Physics Letters, 2008, 92(11): 111903. DOI: 10.1063/1.2896653.
[29]	BENNETT L S, SORRELL F Y, SIMONSEN I K, et al. Ultrafast chemical reactions between nickel and aluminum powders during shock loading [J]. Applied Physics Letters, 1992, 61(5): 520–521. DOI: 10.1063/1.107874.
[30]	IYER K R, BENNETT L S, SORRELL F Y, et al. Solid state chemical reactions at the shock front [J]. AIP Conference Proceedings, 1994, 309(1): 1337–1340. DOI: 10.1063/1.46457.
[31]	EAKINS D, THADHANI N N. Shock-induced reaction in a flake nickel+spherical aluminum powder mixture [J]. Journal of Applied Physics, 2006, 100(11): 113521. DOI: 10.1063/1.2396797.
[32]	GRAHAM R A, ANDERSON M U, HORIE Y, et al. Pressure measurements in chemically reacting powder mixtures with the Bauer piezoelectric polymer gauge [J]. Shock Waves, 1993, 3(2): 79–82. DOI: 10.1007/BF02115887.
[33]	XU X, THADHANI N N. Investigation of shock-induced reaction behavior of as-blended and ball-milled Ni+Ti powder mixtures using time-resolved stress measurements [J]. Journal of Applied Physics, 2004, 96(4): 2000–2009. DOI: 10.1063/1.1773380.
[34]	MCQUEEN R G, MARSH S P, TAYLOR J W, et al. The equation of state of solids from shock wave studies [M]// KINSLOW R. High Velocity Impact Phenomena. New York: Academic Press, 1970.
[35]	KELLY S C, THADHANI N N. Shock compression response of highly reactive Ni+Al multilayered thin foils [J]. Journal of Applied Physics, 2016, 119(9): 095903. DOI: 10.1063/1.4942931.
[36]	AMES R G. Vented chamber calorimetry for impact-initiated energetic materials [C] // Proceedings of the 43rd AIAA Aerospace Sciences Meeting and Exhibit. Reno, USA: AIAA, 2005: 1013. DOI: 10.2514/6.2005-279.
[37]	AMES R G. Energy release characteristics of impact-initiated energetic materials [C] // MRS Online Proceeding Library. UK: Cambridge University Press, 2005: 0896-H03-08. DOI: 10.1557/PROC-0896-H03-08.
[38]	王海福, 刘宗伟, 俞为民, 等. 活性破片能量输出特性实验研究 [J]. 北京理工大学学报, 2009, 29(8): 663–666. DOI: 1001-0645(2009)08-0663-04. WANG H F, LIU Z W, YU W M, et al. Experimental investigation of energy release characteristics of reactive fragments [J]. Transactions of Beijing Institute of Technology, 2009, 29(8): 663–666. DOI: 1001-0645(2009)08-0663-04.
[39]	ZHANG X F, SHI A S, QIAO L, et al. Experimental study on impact-initiated characters of multifunctional energetic structural materials [J]. Journal of Applied Physics, 2013, 113(8): 083508. DOI: 10.1063/1.4793281.
[40]	LUO P G, WANG Z C, JIANG C L, et al. Experimental study on impact-initiated characters of W/Zr energetic fragments [J]. Materials and Design, 2015, 84: 72–78. DOI: 10.1016/j.matdes.2015.06.107.
[41]	XIONG W, ZHANG X F, TAN M T, et al. The energy release characteristics of shock-induced chemical reaction of Al/Ni composites [J]. The Journal of Physical Chemistry C, 2016, 120(43): 24551–24559. DOI: 10.1021/acs.jpcc.6b06530.
[42]	ZHANG S, LIU J X, YANG M, et al. Effects of multi-component co-addition on reaction characteristics and impact damage properties of reactive material [J]. Materials and Design, 2018, 153: 1–8. DOI: 10.1016/j.matdes.2018.04.077.
[43]	MASON B A, GROVEN L J, SON S F. The role of microstructure refinement on the impact ignition and combustion behavior of mechanically activated Ni/Al reactive composites [J]. Journal of Applied Physics, 2013, 114(11): 113501. DOI: 10.1063/1.4821236.
[44]	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.
[45]	EAKINS D E, THADHANI N N. Investigation of shock-induced reactions in a Ni+Al powder mixture [J]. AIP Conference Proceedings, 2006, 845(1): 1153–1156. DOI: 10.1063/1.2263528.
[46]	THADHAN N N. Shock-induced and shock-assisted solid-state chemical reactions in powder mixtures [J]. Journal of Applied Physics, 1994, 76(4): 2129–2138. DOI: 10.1063/1.357624.
[47]	NESTERENKO V F, MEYERS M A, CHEN H C, et al. Controlled high-rate localized shear in porous reactive media [J]. Applied Physics Letters, 1994, 65(24): 3069–3071. DOI: 10.1063/1.112509.
[48]	VANDERSALL K S, THADHANI N N. Investigation of“shock-induced”and“shock-assisted”chemical reactions in Mo+2Si powder mixtures [J]. Metallurgical and Materials Transactions A, 2003, 34(1): 15–23. DOI: 10.1007/s11661-003-0204-z.
[49]	YU L H, MEYERS M A. Shock synthesis and synthesis-assisted shock consolidation of suicides [J]. Journal of Materials Science, 1991, 26(3): 601–611. DOI: 10.1007/BF00588294.
[50]	DUNBAR E, THADHANI N N, GRAHAM R A. High-pressure shock activation and mixing of nickel-aluminium powder mixtures [J]. Journal of Materials Science, 1993, 28(11): 2903–2914. DOI: 10.1007/BF00354693.
[51]	THADHANI N N, GRAHAM R A, ROYAL T, et al. Shock-induced chemical reactions in titanium-silicon powder mixtures of different morphologies: time-resolved pressure measurements and materials analysis [J]. Journal of Applied Physics, 1997, 82(3): 1113–1128. DOI: 10.1063/1.365878.
[52]	ROYAL T E, NAMJOSHI S, THADHANI N N. Mechanistic processes influencing shock chemistry in powder mixtures of the Ti-Si, Ti-Al, and Ti-B systems [J]. Metallurgical and Materials Transactions A, 1996, 27(7): 1761–1771. DOI: 10.1007/BF02651925.
[53]	HORIE Y, GRAHAM R A, SIMONSEN I K. Synthesis of nickel aluminides under high-pressure shock loading [J]. Materials Letters, 1985, 3(9–10): 354–359. DOI: 10.1016/0167-577X(85)90075-8.
[54]	AKASHI T, SAWAOKA A, SAITO S, et al. Structural changes of boron nitride caused by multiple shock-compressions [J]. Japanese Journal of Applied Physics, 1976, 15(5): 891. DOI: 10.1143/JJAP.15.891.
[55]	BOSLOUGH M B. A thermochemical model for shock-induced reactions (heat detonations) in solids [J]. The Journal of Chemical Physics, 1990, 92(3): 1839–1848. DOI: 10.1063/1.458066.
[56]	MEYERS M A. Dynamic behaviour of materials [M]. New York: John Wiley & Sons, 1994. DOI: 10.1002/9780470172278.
[57]	BOLME C A, MCGRANE S D, MOORE D S, et al. Single shot measurements of laser driven shock waves using ultrafast dynamic ellipsometry [J]. Journal of Applied Physics, 2007, 102(3): 033513. DOI: 10.1063/1.2767376.
[58]	VANDERSALL K S, THADHANI N N. Time-resolved measurements of the shock-compression response of Mo+2Si elemental powder mixtures [J]. Journal of Applied Physics, 2003, 94(3): 1575–1583. DOI: 10.1063/1.1586968.
[59]	HUANG C M, CHEN J, BAI S X, et al. Enhancement of energy release performance of Al-Ni composites by adding CuO [J]. Journal of Alloys and Compounds, 2020, 835: 155271. DOI: 10.1016/j.jallcom.2020.155271.
[60]	EAKINS D E, THADHANI N N. Shock compression of reactive powder mixtures [J]. International Materials Reviews, 2009, 54(4): 181–213. DOI: 10.1179/174328009X461050.
[61]	GRAHAM R A. Issues in shock-induced solid state chemistry [R]. Albuquerque: Sandia National Lab, 1989.
[62]	YANO K, HORIE Y. A numerical study of shock-induced particle velocity dispersion in solid mixtures [J]. Journal of Applied Physics, 1998, 84(3): 1292–1298. DOI: 10.1063/1.368197.
[63]	TAMURA S, HORIE Y. Discrete meso-dynamic simulation of thermal explosion in shear bands [J]. Journal of Applied Physics, 1998, 84(7): 3574–3580. DOI: 10.1063/1.368532.
[64]	AUSTIN R A. Numerical simulation of the shock compression of microscale reactive particle systems [D]. Atlanta: Georgia Institute of Technology, 2005.
[65]	BOLKHOVITINOV L G, BATSANOV S S. Theory of solid-state detonation [J]. Combustion, Explosion, and Shock Waves, 2007, 43(2): 219–221. DOI: 10.1007/s10573-007-0030-5.
[66]	ZHANG X F, QIAO L, SHI A S, et al. A cold energy mixture theory for the equation of state in solid and porous metal mixtures [J]. Journal of Applied Physics, 2011, 110(1): 013506. DOI: 10.1063/1.3603018.
[67]	张先锋, 赵晓宁, 乔良. 反应金属冲击反应过程的理论分析 [J]. 爆炸与冲击, 2010, 30(2): 145–151. DOI: 10.11883/1001-1455(2010)02-0145-07. ZHANG X F, ZHAO X N, QIAO L. Theory analysis on shock-induced chemical reaction of reactive metal [J]. Explosion and Shock Waves, 2010, 30(2): 145–151. DOI: 10.11883/1001-1455(2010)02-0145-07.
[68]	HORIE Y, KIPP M E. Modeling of shock-induced chemical reactions in powder mixtures [J]. Journal of Applied Physics, 1988, 63(12): 5718–5727. DOI: 10.1063/1.340309.
[69]	DO I P H, BENSON D J. Micromechanical modeling of shock-induced chemical reactions in heterogeneous multi-material powder mixtures [J]. International Journal of Plasticity, 2001, 17(4): 641–668. DOI: 10.1016/S0749-6419(00)00065-6.
[70]	BENSON D J. The calculation of the shock velocity-particle velocity relationship for a copper powder by direct numerical simulation [J]. Wave Motion, 1995, 21(1): 85–99. DOI: 10.1016/0165-2125(94)00044-6.
[71]	NIEH T G, LUO P, NELLIS W, et al. Dynamic compaction of aluminum nanocrystals [J]. Acta Materialia, 1996, 44(9): 3781–3788. DOI: 10.1016/1359-6454(96)83816-X.
[72]	AUSTIN R A, MCDOWELL D L, BENSON D J. Numerical simulation of shock wave propagation in spatially-resolved particle systems [J]. Modelling and Simulation in Materials Science and Engineering, 2006, 14(4): 537–561. DOI: 10.1088/ 0965-0393/14/4/001.
[73]	BENSON D J. An analysis by direct numerical simulation of the effects of particle morphology on the shock compaction of copper powder [J]. Modelling and Simulation in Materials Science and Engineering, 1994, 2(3A): 535–550. DOI: 10.1088/0965-0393/2/3A/008.
[74]	BAER M R, TROTT W M. Mesoscale studies of shock loaded tin sphere lattices [J]. AIP Conference Proceedings, 2004, 706(1): 517–520. DOI: 10.1063/1.1780290.
[75]	BAER M R. Modeling heterogeneous energetic materials at the mesoscale [J]. Thermochimica Acta, 2002, 384(1–2): 351–367. DOI: 10.1016/S0040-6031(01)00794-8.
[76]	BENSON D J, CONLEY P. Eulerian finite-element simulations of experimentally acquired HMX microstructures [J]. Modelling and Simulation in Materials Science and Engineering, 1999, 7(3): 333–354. DOI: 10.1088/0965-0393/7/3/304.
[77]	EAKINS D E, THADHANI N N. Discrete particle simulation of shock wave propagation in a Ni+Al powder mixture [J]. Journal of Applied Physics, 2007, 101(4): 043508. DOI: 10.1063/1.2431682.
[78]	EAKINS D E, THADHANI N N. Mesoscale simulation of the configuration-dependent shock-compression response of Ni+Al powder mixtures [J]. Acta Materialia, 2008, 56(7): 1496–1510. DOI: 10.1016/j.actamat.2007.12.009.
[79]	CHAWLA N, SIDHU R S, GANESH V V. Three-dimensional visualization and microstructure-based modeling of deformation in particle-reinforced composites [J]. Acta Materialia, 2006, 54(6): 1541–1548. DOI: 10.1016/j.actamat.2005.11.027.
[80]	SIDHU R S, CHAWLA N. Three-dimensional (3D) visualization and microstructure-based modeling of deformation in a Sn-rich solder [J]. Scripta Materialia, 2006, 54(9): 1627–1631. DOI: 10.1016/j.scriptamat.2006.01.013.
[81]	AYYAR A, CHAWLA N. Microstructure-based modeling of the influence of particle spatial distribution and fracture on crack growth in particle-reinforced composites [J]. Acta Materialia, 2007, 55(18): 6064–6073. DOI: 10.1016/j.actamat.2007.06.044.
[82]	ZHAI J, TOMAR V, ZHOU M. Micromechanical simulation of dynamic fracture using the cohesive finite element method [J]. Journal of Engineering Materials and Technology, 2004, 126(2): 179–191. DOI: 10.1115/1.1647127.
[83]	YANG Y, GOULD R D, HORIE Y, et al. Shock-induced chemical reactions in a Ni/Al powder mixture [J]. Applied Physics Letters, 1997, 70(25): 3365–3367. DOI: 10.1063/1.119172.
[84]	QIAO L, ZHANG X F, HE Y, et al. Multiscale modelling on the shock-induced chemical reactions of multifunctional energetic structural materials [J]. Journal of Applied Physics, 2013, 113(17): 173513. DOI: 10.1063/1.4803712.
[85]	乔良. 多功能含能结构材料冲击反应与细观特性关联机制研究[D]. 南京: 南京理工大学, 2013.
[86]	REDING D J, HANAGUD S. Multiscale chemical reactions in reactive powder metal mixtures during shock compression [J]. Journal of Applied Physics, 2010, 105(2): 024912. DOI: 10.1063/1.2976313.
[87]	XIONG Y N, XIAO S F, DENG H Q, et al. Investigation of the shock-induced chemical reaction (SICR) in Ni+Al nanoparticle mixtures [J]. Physical Chemistry Chemical Physics, 2017, 27(19): 17607–17617. DOI: 10.1039/C7CP03176A.

施引文献

期刊类型引用(10)

1.	张甲浩，郭萌萌，周晟，余庆波. 活性横向增强弹靶后横向效应实验与数值模拟. 兵工学报. 2024(07): 2270-2281 . 百度学术
2.	任思远，武强，张品亮，宋光明，陈川，龚自正，李正宇. 活性弹丸超高速撞击蜂窝夹芯板双层结构的损伤特性. 爆炸与冲击. 2024(07): 83-96 . 本站查看
3.	王存洪，曹玉武，陈进，孔霖，孙兴昀. 金属型含能材料力学行为研究进展. 爆炸与冲击. 2023(07): 3-24 . 本站查看
4.	侯先苇，张先锋，熊玮，谈梦婷，刘闯，戴兰宏. 活性无序合金冲击的释能特性及在毁伤元中应用研究进展. 爆炸与冲击. 2023(09): 4-42 . 本站查看
5.	刘振扬，卢冠成，胡蝶，陈鹏万，郑元枫，王海福. PTFE/Al活性破片碰撞三层间隔靶冲击响应行为. 北京理工大学学报. 2023(10): 1068-1077 . 百度学术
6.	丁建，朱顺官. 氟聚物基活性材料释能及毁伤特性研究进展. 含能材料. 2023(08): 844-856 . 百度学术
7.	丁建，朱顺官，陈顺禄，程建辉. 后效毁伤增强聚能战斗部设计与试验. 兵器装备工程学报. 2023(11): 29-37 . 百度学术
8.	任耶平，刘睿，陈鹏万，郭岩松，胡启文，葛超，王海福. Al/PTFE活性材料冲击载荷作用下响应特性研究. 爆炸与冲击. 2022(06): 50-61 . 本站查看
9.	熊玮，张先锋，陈海华，刘闯，谈梦婷. Al/Ni粉末复合材料冲击反应细观机制. 兵工学报. 2022(08): 1823-1834 . 百度学术
10.	李向荣，侯聪，张金忠，张凌云，罗鑫. 碰撞速度对活性弹丸作用多层靶影响研究. 兵器装备工程学报. 2021(09): 66-72 . 百度学术