孔隙度对PBX炸药冲击起爆影响的实验研究

刘海庆; 段卓平; 白志玲; 温丽晶; 欧卓成; 黄风雷

doi:10.11883/bzycj-2018-0226

孔隙度对PBX炸药冲击起爆影响的实验研究

doi: 10.11883/bzycj-2018-0226

刘海庆^1,,
段卓平^1, ,,
白志玲¹,
温丽晶^{1, 2},
欧卓成¹,
黄风雷¹

1.
北京理工大学爆炸科学与技术国家重点实验室，北京 100081
2.
环境保护部核与辐射安全中心，北京 100082

基金项目: 国家自然科学基金(U1630113)；材料与结构冲击动力学创新小组(11521062)

详细信息

作者简介:
刘海庆（1987－　），男，博士研究生，liuhaiqing0705@163.com

通讯作者:
段卓平（1965－　），男，博士，研究员，博士生导师，duanzp@bit.edu.cn

中图分类号: O381
计量
- 文章访问数: 6148
- HTML全文浏览量: 2125
- PDF下载量: 70
- 被引次数: 0
出版历程
- 收稿日期: 2018-06-25
- 修回日期: 2018-08-15
- 刊出日期: 2019-07-01

Experimental research on effects of porosity on shock initiation of PBX explosive

1.
State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
2.
Nuclear and Radiation Safety Center, Ministry of Environmental Protection, Beijing 100082, China

摘要

摘要: 为了研究孔隙度(装药密度)对PBX炸药冲击起爆爆轰成长的影响，采用炸药冲击起爆锰铜压阻一维拉格朗日实验测试系统，测量了不同孔隙度的PBXC03炸药(HMX的质量分数为87%，TATB的质量分数为7%，黏结剂的质量分数为6%)冲击起爆过程不同拉格朗日位置的压力-时间历史。结果显示：在本文装药范围和加载条件下，孔隙度对PBX炸药冲击起爆爆轰过程的影响不单调，中等密度的炸药冲击起爆和爆轰成长最快，这是热点点火过程与燃烧反应过程共同作用的结果。
- 冲击起爆 /
- PBX炸药 /
- 拉格朗日实验分析测试系统 /
- 孔隙度
Abstract: To investigate the influence of the porosity (charge density) on the shock initiation and detonation of polymer bonded explosives (PBXes), a one-dimensional Lagrangian experimental testing system is adopted to measure the pressure-time histories at different Lagrangian locations of PBXC03 (87% HMX, 7% TATB, 6% binder by weight) with three different porosities (or charge densities), in which, an explosive plane-wave lens is used to generate a high-pressure planar detonation wave loading, and the Manganin piezoresistive pressure gauge measurement technique and the attenuation technique by both air-gap and Al-gap are used. The experimental results show that the detonation grows the fastest in the explosive with a moderate porosity. This work provides more detailed experimental data for the further development of mesoscopic reaction rate models for shock initiation of heterogeneous explosives.
- shock initiation /
- PBX explosives /
- Lagrangian experiment testing system /
- porosity

HTML全文

在近几十年中, 对于爆炸焊接的理论研究以及工业应用得到了迅速发展^[1-2]。在炸药爆轰载荷作用下, 2金属板之间高速碰撞形成金属射流和周期性波纹状的复合界面, 通常被称为界面波。金属复合板强度以及复合率与界面波状态密切相关, 金属爆炸焊接只有在复合界面上形成大小适当的波纹才能保证复合板的结合强度。关于爆炸焊接复合界面波的形成机理^[3-4]有4种经典理论:(1)复板侵彻刻入机理, 复板的周期性刻入基板并在驻点前沿形成驼峰而形成波状界面的特征; (2)亥尔姆霍尔兹失稳模型, 复板再入射流和基板之间存在速度差而形成界面波; (3)应力波模型, 在复板与基板的碰撞点处的可压缩冲击波与在自由边界反射的稀疏波相互作用而形成界面波; (4)卡门涡街模型, 在金属板碰撞点附近材料的卡尔曼涡街流动形成界面波。然而, 这些界面波理论模型还不完善, 迄今为止, 爆炸焊接界面波的形成机理仍没有统一的认识。

由于界面波的形成涉及多种物理现象, 单纯使用固体力学计算模型对此很难进行分析, 而采用爆炸焊接实验又难以观察到界面波的动态形成过程, 因此, 本文中基于冲击动力学和物质点法^[5]对爆炸焊接界面波的形成过程进行三维数值模拟, 并开展相关的爆炸焊接实验。

1. 物质点法基本理论

采用物质点法求解冲击动力学问题时, 首先定义背景网格, 然后在背景网格内将连续体剖分为具有集中质量的离散质点, 再利用形函数实现质点和背景网格节点之间的2次映射, 通过求解连续方程和运动方程以及材料本构方程计算所有质点的应变和应力值^[5], 在计算过程中背景网格固定不变。与其他数值计算方法相同, 物质点法满足的连续方程、运动方程和能量方程分别为:

$\frac{\mathrm{d} \rho}{\mathrm{d} t}+\rho \nabla \cdot \boldsymbol{v}=0$

(1)

$\rho \boldsymbol{a}=\nabla \cdot \boldsymbol{\sigma}+\rho \boldsymbol{b}$

(2)

$\rho \dot{e}=\boldsymbol{\sigma}: \dot{\boldsymbol{\varepsilon}}$

(3)

式中:ρ为质量密度, a为加速度, b为单位体积力, v为速度, σ为柯西应力张量, e为质量热力学能, 应变率

将连续体离散为N_p个质点后, 为了利用这些质点来描述连续体的变形过程, 假定质点p(p=1, 2, …, N_P)在t时刻位置为x_p^(t), 质量为M_p, 密度为ρ_p^(t), 速度为v_p^(t), 柯西应力张量为σ_p^(t), 将每个质点看作集中质量点, 并且每个离散质点的质量保持不变, 那么连续方程自然满足。为了求解运动方程则需要采用弱形式^[6]:

$\int_{\Omega} \rho {\boldsymbol{w}} \cdot {\boldsymbol{a}} \mathrm{~d} \Omega=-\int_{\Omega} \rho \boldsymbol{\sigma}^{\mathrm{s}}: \nabla {\boldsymbol{w}} \mathrm{~d} \Omega+\int_{\partial \Omega_{c}} \rho \boldsymbol{c}^{\mathrm{s}} \cdot {\boldsymbol{w}} \mathrm{~d} S+\int_{\Omega} \rho \boldsymbol{w} \cdot \boldsymbol{b} \mathrm{d} \Omega$

(4)

式中:σ^s为比应力张量, σ^s=σ/ρ; Ω为当前连续体构形, ∂Ω_c为指定应力边界; w为试函数; 假设边界上应力为c, 且c^s=c/ρ, w在指定位移边界上为零。这样每个质点密度写成δ函数形式^[6]:

$\rho(\boldsymbol{x}, t)=\sum\limits_{p=1}^{N_{p}} M_{p} \delta\left(\boldsymbol{x}-\boldsymbol{x}_{p}^{(t)}\right)$

(5)

背景网格节点包含了多个参量, 坐标x_i^(t)、位移u_i^(t)、速度v_i^(t)和加速度a_i^(t)以及试函数w_i^(t)等, 统一使用变量r_i^(t)(i=1, 2, …, N_n)表示, 则有^[7]:

$\boldsymbol{r}_{p}^{(t)}=\sum\limits_{i=1}^{N_{n}} \boldsymbol{r}_{i}^{(t)} N_{i}\left(\boldsymbol{x}_{p}^{(t)}\right)$

(6)

利用形函数实现从背景网格节点到质点的映射, 更新质点的位移、速度、应变等参量。由于物质点法采用显式积分算法, 将网格节点插值得到的参数代入式(4)可得:

$\begin{aligned} \sum\limits_{i=1}^{N_{n}} \boldsymbol{w}_{i}^{(t)} \cdot \sum\limits_{j=1}^{N_{n}} m_{i j}^{(t)} \boldsymbol{a}_{j}^{(t)}=-\sum\limits_{i=1}^{N_{n}} \boldsymbol{w}_{i}^{(t)} & \left.\cdot \sum\limits_{p=1}^{N_{p}} M_{p} \boldsymbol{\sigma}_{p}^{(\mathrm{s}, t)} \cdot \nabla N_{i}(\boldsymbol{x})\right|_{\boldsymbol{x} \phi)}+\sum\limits_{i=1}^{N_{n}} \boldsymbol{w}_{i}^{(t)} \cdot\left(\boldsymbol{c}_{i}^{(t)}+\boldsymbol{b}_{i}^{(t)}\right) \\ m_{i j}^{(t)} & =\sum\limits_{p=1}^{N_{p}} M_{p} N_{i}\left(\boldsymbol{x}_{p}^{(t)}\right) N_{j}\left(\boldsymbol{x}_{p}^{(t)}\right) \\ \boldsymbol{c}_{i}^{(t)} & =\sum\limits_{p=1}^{N_{p}} M_{p} \boldsymbol{c}_{p}^{(\mathrm{s}, t)} N_{i}\left(\boldsymbol{x}_{p}^{(t)}\right) / h \\ \boldsymbol{b}_{i}^{(t)} & =\sum\limits_{p=1}^{N_{p}} M_{p} \boldsymbol{b}\left(\boldsymbol{x}_{p}^{(t)}, t\right) N_{i}\left(\boldsymbol{x}_{p}^{(t)}\right) \end{aligned}$

$\sum\limits_{i=1}^{N_{n}} \boldsymbol{w}_{i}^{(t)} \cdot \sum\limits_{j=1}^{N_{n}} m_{i j}^{(t)} \boldsymbol{a}_{j}^{(t)}=-\sum\limits_{i=1}^{N_{n}} \boldsymbol{w}_{i}^{(t)} \left.\cdot \sum\limits_{p=1}^{N_{p}} M_{p} \boldsymbol{\sigma}_{p}^{(\mathrm{s}, t)} \cdot \nabla N_{i}(\boldsymbol{x})\right|_{\boldsymbol{x} \phi)}+\sum\limits_{i=1}^{N_{n}} \boldsymbol{w}_{i}^{(t)} \cdot\left(\boldsymbol{c}_{i}^{(t)}+\boldsymbol{b}_{i}^{(t)}\right)$

(7)

$m_{i j}^{(t)} =\sum\limits_{p=1}^{N_{p}} M_{p} N_{i}\left(\boldsymbol{x}_{p}^{(t)}\right) N_{j}\left(\boldsymbol{x}_{p}^{(t)}\right)$

(8)

$\boldsymbol{c}_{i}^{(t)} =\sum\limits_{p=1}^{N_{p}} M_{p} \boldsymbol{c}_{p}^{(\mathrm{s}, t)} N_{i}\left(\boldsymbol{x}_{p}^{(t)}\right) / h$

(9)

$\mathit{\boldsymbol{b}}_i^{(t)} = \sum\limits_{p = 1}^{{N_p}} {{M_p}} \mathit{\boldsymbol{b}}\left( {\mathit{\boldsymbol{x}}_p^{(t)},t} \right){N_i}\left( {\mathit{\boldsymbol{x}}_\beta ^{(t)}} \right)$

(10)

将运动方程(7)写成以下形式:

$\sum\limits_{j=1}^{N_{n}} m_{i j}^{(t)} \boldsymbol{a}_{j}^{(t)}=\left(f_{i}^{(t)}\right)^{\mathrm{int}}+\left(\boldsymbol{f}_{i}^{(t)}\right)^{\mathrm{ext}} \quad i=1,2, \cdots, N_{n}$

(11)

$\left(\boldsymbol{f}_{i}^{(t)}\right)^{\mathrm{int}}=-\left.\sum\limits_{p=1}^{N_{p}} M_{p} \boldsymbol{\sigma}_{p}^{(\mathrm{s}, t)} \cdot \nabla N_{i}\right|_{x \varphi t)}$

(12)

$\left(\boldsymbol{f}_{i}^{(t)}\right)^{\mathrm{ext}}=\boldsymbol{b}_{p}^{(t)}+\boldsymbol{c}_{p}^{(t)}$

(13)

采用行求和方式获得对角质量阵, 这样式(8)和式(10)分别变为:

$m_{i}^{(t)} =\sum\limits_{p=1}^{N_{p}} M_{p} N_{i}\left( {\boldsymbol{x}_{p,t}^{\left( t \right)}} \right)$

(14)

$\boldsymbol{b}_{i}^{(t)} =m_{i}^{(t)} \boldsymbol{b}\left(\boldsymbol{x}_{p}^{(t)}, t\right)$

(15)

从而式(11)变为:

$m_{i}^{(t)} \boldsymbol{a}_{i}^{(t)}=\left(\boldsymbol{f}_{i}^{(t)}\right)^{\mathrm{int}}+\left(\boldsymbol{f}_{i}^{(t)}\right)^{\operatorname{ext}} \quad i=1,2, \cdots, N_{n}$

(16)

基于材料本构模型对质点应力进行更新, 计算应力增量为:

$\Delta \boldsymbol{\sigma}_{p}^{(t+\Delta t)}=\boldsymbol{E}: \Delta \boldsymbol{\varepsilon}_{p}^{(t+\Delta t)}$

(17)

式中:E是切线刚度张量。计算质点当前时刻的应力张量为:

$\boldsymbol{\sigma}_{p}^{(t+\Delta t)}=\boldsymbol{\sigma}_{p}^{(t)}+\Delta \boldsymbol{\sigma}_{p}^{(t+\Delta t)}$

(18)

计算当前时刻质点内能为:

$e_{p}^{(t+\Delta t)}=e_{p}^{(t)}+\boldsymbol{\sigma}_{p}^{(t+\Delta t)}: \Delta \dot{\boldsymbol{\varepsilon}}_{p}^{(t+\Delta t)} /\left[\rho_{p}^{(t+\Delta t)} \Delta t\right]$

(19)

物质点法的理论模型以及数值计算方法都较成熟, 对于爆炸冲击问题通常采用显式积分算法完成。关于物质点法的具体计算步骤和求解过程可以参考文献[8-12]。

2. 爆炸焊接数值计算模型

2.1 爆炸焊接布置形式

实际爆炸焊接工程中, 一般采用飞板与基板平行的方式进行爆炸焊接, 2金属板之间保持一定间距, 为了避免飞板上表面被炸药烧伤, 还要加一薄缓冲层, 如图 1所示, 爆炸焊接装置主要由雷管、炸药层、缓冲层、飞板、支撑以及基板等组成。

图 1 爆炸焊接布置简图

Figure 1. Distribution pattern of explosive welding

下载: 全尺寸图片幻灯片

在爆炸焊接过程中, 雷管起爆后, 炸药会在起爆端不远处达到稳定的爆轰状态, 爆轰波以v_d的速度向炸药未起爆区域传播, 如图 2所示。

图 2 爆炸焊接过程示意图

Figure 2. The process of explosive welding

下载: 全尺寸图片幻灯片

炸药爆炸产物驱动飞板向下运动, 在预置间隙中经历加速过程, 与基板发生碰撞之前获得较高的飞行速度v_p, 称为碰撞速度。随着炸药稳定爆轰的推进, 碰撞点会以稳定速度v_c向未焊合区域行进, 该速度称为碰撞点移动速度, 平行布置时v_c=v_d。在碰撞点附近会产生高压, 同时金属复合界面会受到剧烈的剪切作用, 进而发生强烈的塑性变形, 导致界面材料升温, 甚至发生熔化流动, 并形成向前的微喷射射流使界面完成自清理, 最终实现爆炸焊接。由于应用物质点法对爆炸焊接界面波的模拟研究计算量很大, 不可能对炸药爆轰、基板与飞板碰撞、焊接界面形成的全部过程进行模拟。同时, 基于爆炸焊接的实际情况, 炸药爆轰驱动飞板, 飞板在预置间隙加速后, 爆轰产物对飞板背面施加的压力已经下降为0.1~1GPa量级, 飞板以一定速度与基板碰撞后的压力却高达10~100GPa量级, 远大于炸药爆轰对飞板的压力, 爆轰压力在2金属板接触后对爆炸焊接界面的影响已经很小。因此, 在对爆炸焊接界面的数值模拟过程中可忽略炸药爆炸驱动过程, 单独对高速碰撞过程进行建模, 爆炸焊接的计算模型如图 3所示。

图 3 爆炸焊接数值计算模型示意图

Figure 3. Calculation model of explosive welding

下载: 全尺寸图片幻灯片

飞板向下以一定速度和倾斜角度β碰撞基板, 选择飞板与基板的倾斜角度应该在爆炸焊接窗口上下限的最大与最小弯折角之间, 焊接窗口上下限是指保证爆炸焊接可以复合的飞板碰撞基板的最大与最小速度, 很多金属的爆炸焊接窗口上下限已经通过实验测定。在已知炸药爆速的前提下, 由飞板碰撞点速度计算公式^[13]:

$v_{\mathrm{p}}=2 v_{\mathrm{d}} \sin (\beta / 2)$

(20)

可以获得本文计算模型飞板的初始倾斜角度, 这样设计的计算模型与轻气炮加速、电磁加速碰撞焊接情况是完全相符的。当确定了飞板与基板材料的爆炸焊接窗口碰撞点速度的上下限以及炸药的爆轰速度后, 由式(20)就能够计算出飞板与基板的碰撞角, 飞板的初始速度、倾斜角度和垂直向下的初始速度也可以由此确定。

2.2 爆炸焊接材料模型

确定爆炸焊接的计算模型后, 假定飞板和基板材料都为普通碳钢, 普通碳钢的爆炸焊接窗口下限为272m/s, 飞板的初始速度设计为300m/s, 碰撞点速度处于焊接窗口中间。由于爆炸焊接整个过程历时较短, 在微秒量级范围内, 所以忽略热传导并假定为绝热过程。2块金属板碰撞过程中, 复合界面材料发生塑性流动。采用物质点法计算与有限元法类似, 为了描述爆炸焊接复合界面处的金属塑性流动变形, 材料模型采用Johnson-Cook塑性模型^[14]:

$\sigma_{y}=\left[A+B\left(\varepsilon_{e}^{p}\right)^{n}\right]\left(1+C \ln \dot{\varepsilon}^{*}\right)$

(21)

式中:σ_y为Von Mises流动应力, ε_e^p为等效塑性应变, 为相对等效塑性应变率, 取1s-¹, B、n为材料应变硬化特性参数, C为材料应变率敏感特性参数; 式中第1个乘积因子表示当 ${{\dot \varepsilon }^*}$ =1.0时流动应力同等效塑性应变之间的关系, 第2个乘积因子表示应变率强化效应。

爆炸焊接过程还涉及到材料在高压条件下的力学行为, 所以对其还要运用相应的状态方程进行描述, 这里材料状态方程选用Mie-Grüneisen状态方程。当材料处于拉伸膨胀状态时, 状态方程形式为:

$p=\rho_{0} c_{\mathrm{g} \mu}^{2} \mu+\left(\gamma_{0}+a \mu\right) E$

(22)

当材料处于压缩状态时, 状态方程形式为:

$p=\frac{\rho_{0} c_{\mathrm{g}}^{2} \mu\left(1+\left(1-\frac{\gamma_{0}}{2}\right) \mu-\frac{a}{2} \mu^{2}\right)}{\left(1-\left(s_{1}-1\right) \mu-s_{2} \frac{\mu^{2}}{\mu+1}-s_{3} \frac{\mu^{3}}{\mu+1}\right)^{2}+\left(\gamma_{0}+a \mu\right) E}$

(23)

式中:c_g、a、s₁、s₂和s₃是与材料冲击压缩特性有关的常数, c_g是冲击波速度(u_s)-质点速度(u_p)曲线的截距, s₁、s₂和s₃是u_s-u_p曲线斜率的系数, a是一阶体积修正系数,γ₀是Grüneisen系数, E为材料比内能, 体积变化率μ=ρ/ρ₀-1。

本文中普通钢状态方程模型参数^[15-16]为:ρ=7 850kg/m³, a=0.43, G=124MPa, A=400MPa, B=220MPa, n=0.08, C=0.025, c_g=4 573m/s, s₁=1.33, s₂=0, s₃=0, γ₀=1.07。

3. 爆炸焊接界面波形成数值模拟

假设爆炸焊接计算模型的尺寸板X方向板长100mm, Y方向板厚20mm, Z方向板宽16mm, 在Z方向的板两侧面施加固定约束条件。利用物质点法对本文中设计的爆炸焊接计算模型进行三维数值模拟之前, 需要划分背景网格和质点单元等, 三维背景网格单元为六面体立方体单元, 网格单元大小为2mm, 在X、Y、Z3个方向划分的背景网格单元数量分别为100、100和50, 每个背景网格单元内的质点数量划分为8个, 飞板和基板的材料模型选择Johnson-Cook模型, 材料在高压下的力学行为选用Mie-Grüneisen状态方程进行描述。爆炸焊接计算模型的三维前处理结果见图 4。

图 4 爆炸焊接数值计算模型前处理

Figure 4. Preprocess of explosive welding

下载: 全尺寸图片幻灯片

爆炸焊接数值模拟是利用自主研发的物质点法程序SPS2.0完成的, 该程序是应用C++编制实现的, 并提供了基本的后处理功能, 图 5给出了爆炸焊接数值模拟的三维变形图, 图 6给出了爆炸焊接数值模拟的中间剖面变形图。爆炸焊接中间剖面的等效应力σ_e后处理云图, 见图 7。爆炸焊接中间剖面的压力p场后处理云图, 见图 8。

图 5 爆炸焊接复合界面波的形成

Figure 5. Interfacial wave formation of explosive welding

下载: 全尺寸图片幻灯片

图 6 爆炸焊接复合界面波x－y中间剖面

Figure 6. Interfacial wave formation of explosive welding at x－ysection

下载: 全尺寸图片幻灯片

图 7 爆炸焊接x－y中间剖面等效应力

Figure 7. Effective stress of explosive welding at x－ysection

下载: 全尺寸图片幻灯片

图 8 爆炸焊接x－y中间剖面压力

Figure 8. Pressure of explosive welding at x－ysection

下载: 全尺寸图片幻灯片

为了验证物质点法计算的准确性, 本文中设计了普通钢板爆炸焊接实验:飞板和基板采用平行布置, 基板和飞板的长度均为1 000mm, 宽度均为400mm, 厚度均为20mm; 炸药采用铵油, 爆速为2 400m/s; 2金属板之间的支撑间距为10~12mm, 保证飞板碰撞速度v_p在爆炸焊接窗口内, 约300m/s。爆炸焊接实验后, 取复合板中间部分做为试样进行研磨, 利用电子金相显微镜获得复合界面金相显微组织照片, 如图 9所示。

图 9 爆炸焊接复合界面波状态金相照片

Figure 9. Metallograph of interfacial wave formation

下载: 全尺寸图片幻灯片

将爆炸焊接复合界面的金相组织以及界面波纹形状(图 9)与数值模拟结果(图 6)进行比较可知, 两者的界面波形状基本一致, 由此说明采用物质点法对爆炸焊接界面波形成问题的模拟结果具有一定的准确性和可靠性。

分析复合界面的金相照片, 可以发现界面处的金属材料有熔化现象, 这与数值模拟结果(图 7)是相符的。在飞板与基板碰撞点位置附近计算的等效应力等于零, 说明飞板与基板在冲击碰撞位置的材料几乎处于流体状态。为什么在碰撞点附近区域出现熔化流体现象, 其中重要原因是该处的压力最高。数值计算获得的压力大小在12GPa以上, 远高于材料强度。通过数值模拟得到的爆炸焊接压力分布图(图 8)可以直观地说明这一点。

4. 结论

(1) 爆炸焊接界面波的形成机理是由于飞板与基板的高速碰撞, 在碰撞点处压力最高, 材料发生熔化形成流体涡旋变形并且迅速冷却而复合为一体, 这个过程随着碰撞点的移动不断地向前推进, 使得2金属板复合形成周期性波状界面。

(2) 物质点法的三维数值模拟结果与爆炸焊接实验结果比较符合, 由此说明物质点法用于模拟爆炸焊接界面波的形成过程具有实际工程应用价值。

(3) 此外, 基于本文数值模拟结果也证明了物质点法用于爆炸焊接和冲击动力学问题计算具有突出的优点:一是避免了有限元法在复合界面处难以克服的网格畸变问题; 二是与其他无网格法相比, 物质点法计算效率更高, 本文数值模拟质点数量10万多, 台式机计算时间约1h, 而利用光滑粒子流体动力学法则需要几天时间。因此, 物质点法在爆炸焊接、高速冲击、塑性大变形等力学问题计算中具有广泛的应用前景。

图 1 炸药冲击起爆锰铜压阻一维拉格朗日实验测试系统^[11-13]

Figure 1. One-dimensional Lagrangian system with manganin piezoresistive pressure gauges for testing explosive shock initiation^[11-13]

下载: 全尺寸图片幻灯片

图 2 PBXC03炸药样品

Figure 2. PBXC03 samples

下载: 全尺寸图片幻灯片

图 3 封装的锰铜压阻传感器

Figure 3. Packaged manganin piezoresistive pressure gauges

下载: 全尺寸图片幻灯片

图 4 示波器记录的典型实验信号

Figure 4. Typical voltage signals recorded by an oscilloscope

下载: 全尺寸图片幻灯片

图 5 不同孔隙度的PBXC03炸药冲击起爆过程的实验测试结果

Figure 5. Experimental pressure-time curves of the shock initiation processes of the PBXC03 with different porosities.

下载: 全尺寸图片幻灯片

图 6 同一加载压力下不同实验中PBXC03炸药0 mm位置的压力-时间曲线对比

Figure 6. Pressure-time curves at 0 mm of the PBXC03 with different porosities in different experiments under the same loading pressure

下载: 全尺寸图片幻灯片

图 7 不同孔隙度的PBXC03炸药前导冲击波阵面压力成长历史和前导冲击波迹线。

Figure 7. Pressure growth histories on the shock wave front and precursory shock wave trajectories in the PBXC03 with different porosities

下载: 全尺寸图片幻灯片

参考文献(15)

[1]	AN C, LI H, YE B, et al. Preparation and characterization of ultrafine HMX/TATB explosive co-crystals [J]. Central European Journal of Energetic Materials, 2017, 14(4): 876–887. DOI: 10.22211/cejem/77125.
[2]	TALAWAR M B, AGARWAL A P, ANNIYAPPAN M, et al. Method for preparation of fine TATB (2-5 microm) and its evaluation in plastic bonded explosive (PBX) formulations [J]. Journal of Hazardous Materials, 2006, 137(3): 1848. DOI: 10.1016/j.jhazmat.2006.05.031.
[3]	WANG Z, GUO X, WU F, et al. Preparation of HMX/TATB composite particles using a mechanochemical approach [J]. Propellants, Explosives, Pyrotechnics, 2016, 41(2): 327–333. DOI: 10.1002/prep.201500136.
[4]	GREBENKIN K F. Comparative analysis of physical mechanisms of detonation initiation in HMX and in a low-sensitive explosive (TATB) [J]. Combustion, Explosion, and Shock Waves, 2009, 45(1): 78–87. DOI: 10.1007/s10573-009-0011-y.
[5]	MASSONI J, SAUREL R, BAUDIN G, et al. A mechanistic model for shock initiation of solid explosives [J]. Physics of Fluids, 1999, 11(3): 710–736. DOI: 10.1063/1.869941.
[6]	BOURNE N K, MILNE A M. Shock to detonation transition in a plastic bonded explosive [J]. Journal of Applied Physics, 2004, 95(5): 2379–2385. DOI: 10.1063/1.1644632.
[7]	COOPER M, TROTT W. On the development of a modified wedge test for shock-to-detonation transition in explosives using ORVIS [C] // 17th Biennial International Conference of the APS Topical Group on Shock Compression of Condensed Matter. Chicago, Illinois, 2011.
[8]	URTIEW P A, TARVER C M. Shock initiation of energetic materials at different initial temperatures: review [J]. Combustion, Explosion, and Shock Waves, 2005, 41(6): 766–776. DOI: 10.1007/s10573-005-0085-0.
[9]	VANDERSALL K S, TARVER C M, GARCIA F, et al. Shock initiation experiments on PBX9501 explosive at 150 ℃ for ignition and growth modeling [C] // American Physical Society Meeting on Shock Compression of Condensed Matter, Baltimore, MD, United States, 2005.
[10]	白志玲, 段卓平, 景莉, 等. 飞片冲击起爆高能钝感高聚物粘结炸药的实验研究 [J]. 兵工学报, 2016, 37(8): 1464–1468. DOI: 10.3969/j.issn.1000-1093.2016.08.018. BAI Zhiling, DUAN Zhuoping, JING Li, et al. Experimental research on initiation of insensitive high energy plastic bonded explosive by flyer impact [J]. Acta Armamentarii, 2016, 37(8): 1464–1468. DOI: 10.3969/j.issn.1000-1093.2016.08.018.
[11]	温丽晶, 段卓平, 张震宇, 等. 不同加载压力下炸药冲击起爆过程实验和数值模拟研究 [J]. 兵工学报, 2013, 34(3): 283–288. DOI: 10.3969/j.issn.1000-1093.2013.03.005. WEN Lijing, DUAN Zhuoping, ZHANG Zhenyu, et al. Experimental and numerical study on the shock initiation of pbxc03 explosive under the different loading pressure [J]. Acta Armamentarii, 2013, 34(3): 283–288. DOI: 10.3969/j.issn.1000-1093.2013.03.005.
[12]	温丽晶, 段卓平, 张震宇, 等. HMX基和TATB基PBX炸药爆轰成长差别的实验研究 [J]. 爆炸与冲击, 2013, 33(S): 135–139. WEN Lijing, DUAN Zhuoping, ZHANG Zhenyu, et al. Experimental research on differences of detonation growth process between HMX-based and TATB-based plastic bonded explosives [J]. Explosion and Shock Waves, 2013, 33(S): 135–139.
[13]	WEN L J, DUAN Z P, ZHANG L S, et al. Effects of HMX particle size on the shock initiation of PBXC03 explosive [J]. International Journal of Nonlinear Sciences and Numerical Simulation, 2012, 13(2): 189–194. DOI: 10.1515/ijnsns.2011.129.
[14]	GIBSON L L, DATTELBAUM D, BARTRAM B, et al. Shock initiation sensitivity and Hugoniot-based equation of state of composition-B Obtained using in-situ electromagnetic gauging [J]. Journal of Physics Conference Series, 2014, 500(19): 192004. DOI: 10.1088/1742-6596/500/19/192004.
[15]	DUAN Z P, LIU Y, PI A G, et al. Foil-like manganin gauges for dynamic high pressure measurements [J]. Measurement Science and Technology, 2011, 22(7): 075206. DOI: 10.1088/0957-0233/22/7/075206.

施引文献

资源附件(0)

访问统计

图(7)

计量

文章访问数: 6148
HTML全文浏览量: 2125
PDF下载量: 70
被引次数: 0

1. 物质点法基本理论
2. 爆炸焊接数值计算模型
2.1 爆炸焊接布置形式
2.2 爆炸焊接材料模型
3. 爆炸焊接界面波形成数值模拟
4. 结论

孔隙度对PBX炸药冲击起爆影响的实验研究

doi: 10.11883/bzycj-2018-0226

作者简介: 刘海庆（1987－ ），男，博士研究生，liuhaiqing0705@163.com

通讯作者: 段卓平（1965－ ），男，博士，研究员，博士生导师，duanzp@bit.edu.cn

计量

出版历程