高能炸药爆速曲率关系的一种计算方法

研究平面爆轰波爆轰化学能和状态方程之间的关系,得到平面爆轰波关系式和CJ条件.基于柱坐标系下的一维流体力学方程组,采用坐标变换方法得到含曲率的爆轰波关系式及广义声速条件.在反应区内采用压力平衡和温度平衡的混合规则,在给定的复杂状态方程和复杂反应率下求解,得到PBX9404和PBX9502炸药的爆速和曲率的关系.已知炸药...     

炸药爆轰的连续介质本构模型和数值计算方法

假定炸药和爆轰产物处于局部热力学平衡状态, 即它们的压力和温度相同, 利用热力学基本关系建立炸药爆轰过程的连续介质本构模型的一般理论框架. 在此框架下, 炸药爆轰本构模型由一组常微分方程构成, 包括炸药和爆轰产物的状态方程、简单混合法则、化学反应速率方程和能量守恒方程, 易于由成熟的计算方法如梯形法等进行求解. 一组广义Maxwell型非线性固体本构形式的微分方程描述了压力和温度随时间的演化速率与应变率和化学反应速率的关系, 借助简单混合物理论, 其中的系数由炸药和爆轰产物的材料参数确定. 未反应的炸药和爆轰产物采用JWL状态方程, 化学反应率方程采用Lee-Tarver点火-燃烧二项式模型, 模拟PBX-9404炸药的一维冲击波起爆过程和爆轰波传播过程. 计算结果表明了本文给出的本构模型和相应计算方法的有效性. 关键词： 炸药爆轰 本构模型 化学反应率方程 数值模拟     

炸药枪击试验和数值模拟研究

HL-10炸药是一种以RDX为基的含铝炸药,为了研究该炸药在子弹或金属破片撞击作用下的安全性,利用12.7 mm机枪法对钢壳包覆的柱形HL-10装药进行了枪击试验,试验结果表明,炸药没有发生燃烧或爆轰现象,由此可定义该炸药的枪击感度试验反应等级为1级。建立了炸药枪击试验的计算模型和数值计算方法,对子弹撞击和穿透炸药过程进行了三维数值模拟计算,计算结果与试验结果相符,分析了子弹速度对HL-10炸药枪击感度的影响,其结果可为炸药安全性评价分析提供理论根据。为了进一步验证计算模型和方法,对美国的PBX-9404炸药的枪击作用过程进行了计算分析,结果表明,PBX-9404炸药在枪击作用下发生了完全爆轰反应,与Neff等人的试验结果相吻合。     

CTVD格式数值计算非均质炸药爆轰问题

将高分辨率激波捕捉格式CTVD格式拓展应用到非均质炸药爆轰的数值模拟问题.增加了化学反应率控制方程,引入Lee-Tarver点火成长模型,未反应的固体炸药和化学反应气体产物都使用JWL形式状态方程.数值模拟了非均质固体炸药PBX-9404和TATB的冲击起爆问题.获得了较高的爆轰波分辨率和光滑解区的数值精度,对具有复杂物态方程形式的固体炸药爆轰问题,CTVD格式具有简单实用、高效和高分辨的特点.     

猛炸药RHT-901和钝感炸药IHE-2的爆轰波直角绕射

 采用高速转镜分幅相机和电探针技术研究了猛炸药RHT-901和钝感炸药IHE-2的爆轰波直角绕射图像和不同位置上的爆轰波传播时间。从研究得出，两种炸药都在拐角顶点附近绕射，爆轰波传播时间增长，爆速变小。但是两种炸药绕射爆轰波的状态不一样，钝感炸药IHE-2中爆轰波绕过直角时，在拐角顶点附近约10 mm范围内炸药未完全反应，猛炸药RHT-901中爆轰波绕过直角时未出现类似现象。两者相比，钝感炸药中绕射爆轰波速度变化大，波阵面曲率半径小，而猛炸药的绕射爆轰波速度变化小，波阵面曲率半径大。这说明炸药的爆轰波绕射与炸药的冲击感度、反应区宽度有关。     

利用电磁法研究HMX与TATB混合钝感炸药的冲击起爆特性

为研究含有少量奥克托金(HMX)且以三氨基三硝基苯(TATB)为基的高能钝感炸药PBX-3的冲击起爆反应增长规律,采用火炮驱动蓝宝石飞片的方法和铝基组合式电磁粒子速度计技术进行了一维平面冲击实验。通过实验测量撞击表面及内部不同深度处的冲击波后粒子速度,得到PBX-3炸药的Hugoniot关系。根据冲击波示踪器所测数据绘制了炸药到爆轰的时间-距离(x-t)图,获得了反映炸药冲击起爆性能的Pop关系。将入射压力为12.964 GPa时达到爆轰的6条速度曲线修整成相同零点,通过读取6条曲线的分离点即反应区末端的C-J点,计算出化学反应区时间和宽度。     

20 GPa斜波压缩下PBX-14炸药的动力学响应

利用磁驱动加载实验技术和激光干涉测速技术,开展了未反应固体TATB基PBX-14炸药的斜波压缩实验,获得了20 GPa峰值压力下PBX-14炸药的后表面速度波剖面实验数据。基于阻抗匹配修正的迭代Lagrange数据处理方法处理实验数据,获得了0~20 GPa压力范围内PBX-14炸药的压力-相对比容关系、高压声速-粒子速度关系等动力学特性参数。结合等熵状态方程和由实验获得的动力学参数,对PBX-14炸药的斜波压缩实验过程开展了一维流体动力学数值模拟,计算结果与实验结果吻合良好,验证了本实验方法、数据处理方法及选取的物理模型的正确性。     

凝聚炸药爆轰波在高声速材料界面上的折射现象分析

凝聚炸药爆轰在边界高声速材料约束下传播时,爆轰波会在约束材料界面上产生复杂的折射现象.本文针对凝聚炸药爆轰波在高声速材料界面上的折射现象展开理论和数值模拟分析.首先通过建立在爆轰ZND模型上的改进爆轰波极曲线理论给出爆轰波折射类型,然后发展一种求解爆轰反应流动方程的基于特征理论的二阶单元中心型Lagrange计算方法来数值模拟典型的爆轰波折射过程.从改进爆轰波极曲线理论和二阶Lagrange方法数值模拟给出的结果看出,凝聚炸药爆轰波在高声速材料界面上的折射类型有四种:反射冲击波的正规折射、带束缚前驱波的非正规折射、带双Mach反射的非正规折射、带λ波结构的非正规折射.     

JO-9159炸药的初始密度对爆轰波反应区宽度影响的实验研究

 利用电磁法对两种具有不同初始密度的JO-9159炸药的爆轰反应区宽度进行了测量。结果表明,对于JO-9159这种爆轰反应区非常窄的炸药,既使爆轰达到了定常状态,化学反应仍然与炸药的初始物理状态有关。     

弹粘塑性双球壳塌缩热点反应模型

 基于Kim的弹粘塑性单球壳塌缩模型,考虑PBX炸药中的粘结剂效应,假设炸药和粘结剂均为弹粘塑性材料,建立了弹粘塑性双球壳塌缩热点反应模型,给出了炸药球壳在冲击压力作用下的速度、应变、温度和化学反应速率的时空分布,以及新的热点反应速率理论表达式。把新的热点反应项与Kim的低压下慢反应项和张震宇提出的高压反应速率方程相结合,得到了新的冲击起爆三项式细观反应速率模型。把该模型加入DYNA2D中,模拟了PBX-9501炸药的一维冲击起爆过程,结果表明：该模型除了可以解释炸药颗粒度和孔隙度的影响外,还可以较好地描述粘结剂强度和含量对PBX炸药冲击起爆感度的影响。     

The SURF model and the curvature effect for PBX 9502

SURF is a high explosive burn model based on the ignition & growth concept of hot-spot reaction. For the TATB based explosive PBX 9502, the model has been calibrated to shock-to-detonation transition experiments. To apply the SURF model for propagating detonation waves, the rate has to be extended to a higher pressure regime than is sampled by shock initiation experiments. The experimentally measured curvature effect – detonation speed as a function of front curvature or D _n(κ) – provides the appropriate data for calibrating the propagation regime. The calibration to the curvature effect is based on the ODEs for the reaction zone profile of a detonation wave in conjunction with a shooting algorithm to determine the rate model parameters, for a given κ, needed to obtain a specified detonation speed. A complication for calibrating PBX 9502 rate models arises from the kink in the experimentally measured D _n(κ) curve. This results from the fast and slow reactions that TATB exhibits. To account for this, we use an extension of the SURF model that utilises a sequence of two reactions. The first, with a fast rate, is due to molecular decomposition and is described by the original SURF formulation. The second, with a slow rate, is due to carbon clustering and is used to contribute additional energy from the formation of carbon bonds. The wave profile equations are generalised to the SURF-plus model. Model parameters are then determined for the propagation regime to fit the curvature effect data. The extended model is applicable to both the shock initiation regime and the propagating detonation wave regime.     

Experimental measurement of energy release from an initiating layer in an insensitive explosive

When subjected to a shock of insufficient strength to trigger prompt reaction, heterogenous condensed phase explosives can form regions where significant amounts of the explosive remain unreacted for times much greater than the reaction time of the detonating explosive. This phenomena is observed for the explosive PBX 9502 (95 wt% TATB) both for planar and oblique input shocks. In this work, we build on previous results by performing cylinder expansion (CYLEX) tests where the explosive charge is comprised of a faster core of PBX 9501 (95 wt% HMX) inside a slower annulus of PBX 9502. The detonation in the faster PBX 9501 drives an oblique shock into the adjacent PBX 9502, and an annular transverse initiating layer (IL) results. In the test geometry, the IL travels steadily down the length of the test after a short run distance. At radial positions beyond the IL, an annular region of detonating PBX 9502 is observed. Using standard CYLEX test diagnostics, we infer the total energy release of this experiment. By making the assumptions that (1) the combined energy release is comprised of contributions from detonating PBX 9501, detonating PBX 9502, and the IL in the PBX 9502 and (2) mass-specific energy release for the detonating explosives is approximately the same as typically observed for each explosive, the IL energy release and reaction efficiency can be computed. Results are compared to prior results for a similar geometry, and indicate that while shock deadened PBX 9502 does not detonate promptly, it does eventually release a significant portion of its chemical potential energy over longer timescales on the order of 10 µs.     

Sensitivity of run-to-detonation distance in practical explosives

The dynamics of detonation in a granular explosive following a piston impact is examined computationally for a two-phase model. Different choices are considered for equations of state and reaction rate. Of special interest is the behaviour of the run-to-detonation distance as a function of the initial porosity of the explosive, for which new experimental information has recently become available. It is found that this response can vary both qualitatively and quantitatively depending upon the constitutive input to the model. Computations based upon up-to-date equation-of-state and reaction-rate information for the explosive PBX-9501 show that the response of the run-to-detonation distance as a function of the initial porosity is in the shape of an inverted U, which is in qualitative agreement with the latest experiments. Mechanisms responsible for this behaviour are identified.     

Experimental and modeling analysis of detonation in circular arcs of the conventional high explosive PBX 9501

We examine the diffraction dynamics of a two-dimensional (2D) detonation in a circular arc of the conventional HMX-based, high performance, solid explosive PBX 9501, for which the detonation reaction zone length scale is estimated to be of the order of 100–150 µm. In this configuration, a steady propagating detonation will develop, sweeping around the arc with constant angular speed. We report on results from three PBX 9501 arc experiments, exploring the variation in linear speed on the inner and outer arc surfaces for the steady wave along with the structure of the curved detonation front, as a function of varying inner surface radius and arc thickness. Comparisons of the properties of the motion of the steady wave for each arc configuration are then made with a spatially-distributed PBX 9501 reactive burn model, calibrated to detonation performance properties in a 2D planar slab geometry. We show that geometry-induced curvature of the detonation near the inner arc surface has a significant effect on the detonation motion even for conventional high explosives. We also examine the detonation driving zone structure for each arc case, and thus the subsonic regions of the flow that determine the influence of the arc geometry on the detonation propagation. In addition, streamline paths and reaction progress isolines are calculated. We conclude that a common approximation for modeling conventional high explosive detonation, wherein the shock-normal detonation speed is assumed equal to the Chapman–Jouguet speed, can lead to significant errors in describing the speed at which the detonation propagates.