HNS-Ⅳ炸药的点火增长模型

进行直径?10mm 的圆筒实验,根据冲击Hugoniot关系,拟合了HNS-Ⅳ炸药的未反应时的JWL 状态方程,并通过数值模拟确定了HNS-Ⅳ炸药爆轰产物JWL状态方程参数;测量炸药/窗口界面粒子速度 历程,结合数值模拟获得HNS-Ⅳ炸药的点火增长模型反应速率方程参数。研究表明,HNS-Ⅳ炸药点火增长 模型能够描述反应过程并与实验结果吻合较好。     

爆轰计算JWL状态方程参数的不确定度

针对爆轰计算JWL状态方程参数,开展不确定度量化研究。首先从JWL状态方程性质出发,结 合数学、物理分析确定JWL状态方程中主要参数R1 的概率分布范围,然后利用非嵌入多项式混沌展开方法 对PBX9404炸药驱动飞层外界面速度进行了不确定度量化分析。给出的分析方法与分析流程对爆轰计算及 其他工程数值模拟不确定度量化有一定参考作用。     

PBX-9501强爆轰产物状态方程的数值模拟

利用P.K.Tang的JWL强爆轰产物状态方程改造了DEFEL二维弹塑性流有限元源程序,并用改造后的DEFEL程序模拟了一系列飞片撞击PBX-9501炸药的一维强爆轰过程,给出了炸药驱动二级飞片(铝)的速度。计算得到的产物状态参数与实验测量结果符合较好,验证了增加修正项的新状态方程可较好地描述高能炸药的超压爆轰的高压流体动力学过程。     

PBX-01炸药水中爆轰产物状态方程研究

提出通过水中实验确定炸药的水中爆轰产物JWL状态方程参数的方法;选择PBX-01高能炸药进行水中实验,利用ANSYS/LS-DYNA程序建立炸药的水中实验模型,将实验结果与数值计算结果进行对比,确定PBX-01炸药水中爆轰产物的JWL状态方程参数。研究结果显示,圆筒实验确定的JWL参数在反映炸药水中爆轰产物的膨胀状态时有所不足,水中实验确定的JWL状态方程参数能够更准确地描述PBX-01炸药水中爆轰产物的膨胀过程,因此对水中爆炸的研究需要通过水中爆炸实验建一套状态方程参数。     

RDX基PBX的做功能力及JWL状态方程参数确定

为研究RDX基PBX炸药的做功能力并确定其爆轰产物的JWL状态方程参数,对RDX基PBX炸药和TNT炸药进行?50 mm标准圆筒实验,获得了圆筒膨胀位移和速度的时程曲线,对比得出RDX基PBX炸药的做功能力明显高于TNT炸药;基于能量守恒对实验数据进行非线性拟合,得到2种炸药爆轰产物的JWL状态方程参数。TNT炸药的拟合参数和通过AUTODYN软件计算得到的结果符合较好;将采用上述方法得到的RDX基PBX炸药爆轰产物JWL状态方程参数用于数值模拟,计算结果与实验值吻合较好,符合数值模拟标定JWL状态方程参数的要求。     

一种与爆轰参数封闭的JWL方程参数确定方法

利用数值模拟和理论计算分析了炸药JWL状态方程参数与爆轰参数封闭的重要性;获得了利用圆筒实验测试结果计算炸药爆轰产物绝热等熵指数和爆压的方法;建立了与爆轰参数封闭的JWL状态方程参数的确定方法,并依据公布的圆筒实验数据,应用于两种典型炸药JWL状态方程参数的确定,获得的参数与炸药爆轰参数严格封闭,数值模拟结果与实验结果一致性好,表明炸药JWL参数确定方法合理可靠。     

不同铝粉尺寸含铝炸药加速金属能力的研究

对铝粉直径从几十纳米到几十微米的几种含铝炸药,进行了小尺寸装药条件下炸药加速金属平板实验,用激光速度干涉仪测量了金属平板自由面速度,比较了几种炸药对不同厚度金属平板的加速时间和加速能力。分析了铝粉在炸药爆轰中的反应情况。在小尺寸装药下,铝粉加入炸药和铝粉尺寸大小对炸药爆轰性能有明显影响。纳米铝粉使炸药作功能力有较大提高。含铝炸药爆轰中,铝主要在爆轰后期与爆轰产物反应。     

RDX基含铝炸药爆轰波结构实验研究

为了获得含铝炸药爆轰反应区附近铝粉的反应情况,对两种RDX/Al炸药和一种RDX/LiF炸药的爆轰波结构进行了测量。实验过程中,利用火炮加载产生一维平面波,通过光子多普勒测速仪测量炸药/LiF窗口的界面粒子速度。结果表明:含铝炸药爆轰波的结构与理想炸药的差异较大,其界面粒子速度曲线没有明显的拐点;反应初期,由于气相产物与添加物之间温度的非平衡性,RDX/Al界面的粒子速度低于RDX/LiF炸药的;随后,由于铝粉反应放能,RDX/Al界面的粒子速度高于RDX/LiF炸药的;微米尺度铝粉在CJ面前几乎不发生反应;2、10 μm等两种粒度铝粉的反应延滞时间小于0.8 μs;在本文中,两种粒度铝粉的反应度为16%~31%。     

CL-20混合炸药的爆轰波结构

基于爆轰数值模拟计算,分析了CL-20混合炸药爆轰反应的特征,设计了炸药与窗口的界面粒子速度测量实验装置;采用激光干涉法,测量了C-1炸药(CL-20/粘合剂/94/6)与窗口的界面粒子速度; 运用先求导、再分段拟合的方法,对界面粒子速度随时间的变化曲线进行了数据处理,确定了炸药爆轰CJ点对应的时间位置;根据CJ点对应的粒子速度,计算获得了炸药的爆轰反应区宽度和CJ爆轰压力。结果显示:密度为1.943 g/cm3的C-1炸药的爆轰反应时间为38 ns,CJ压力为34.2 GPa。     

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

针对含铝炸药爆轰的非理想特性,提出了含铝炸药爆轰产物膨胀的局部等熵假设,建立含铝炸药爆轰驱动的非线性特征线模型,为研究含铝炸药爆轰产物的非等熵流动和膨胀做功提供了一种新的理论分析方法。设计了5、50 μm含铝炸药和含LiF炸药驱动0.5、1 mm厚金属板实验,通过激光位移干涉仪测试金属板运动的速度历程,再通过实验结果计算得到铝粉在爆轰产物中的反应度变化规律,结合含铝炸药爆轰产物的非线性特征线模型,理论计算了含铝炸药驱动金属板的速度历程。对比理论与实验结果,理论方法能够很好地描述铝粉二次反应对炸药做功能力的贡献,同时验证了含铝炸药爆轰驱动的非线性特征线模型的正确性。     

膨胀石墨燃爆剂JWL状态方程参数拟合

用LS-DYNA对膨胀石墨燃爆剂爆炸过程进行数值模拟时,需要提供膨胀石墨燃爆剂的JWL状态方程的参数。把凝聚炸药等熵线物态方程作为目标方程,通过差分进化法拟合得出膨胀石墨燃爆剂JWL状态方程的参数。通过圆筒实验对所求得的参数进行验证,结果表明,理论值与实验值最大误差不超过3.3%,能够满足膨胀石墨燃爆剂爆炸模拟研究的需要。     

爆轰产物状态方程的水下爆炸反演理论研究

由于水下爆炸过程中爆轰产物的信息以水中压缩波形式向外传递，本文旨在探讨如何利用水下爆炸试验数据确定爆轰产物的状态方程。相较于常规圆筒试验，水下爆炸试验具有装置简单成本低、装药尺寸限制少、测定压力范围更广的特点，更适用于大药量或非理想炸药的现场测试。本文从水中冲击波轨迹和波后压力时程曲线出发，发展了由冲击波及其波后流场还原水气界面的逆特征线算法，以及根据水气界面确定爆轰产物状态方程的遗传算法。与水下爆炸正演结果对比，发现逆特征线法可以较准确地还原水气界面的位置和压力参数，有效压力下限可达2 MPa，远低于圆筒试验的测试下限0.1 GPa。根据水气界面的反演结果，遗传算法也能稳定地优化JWL方程参数，8种常用炸药的等熵衰减压力误差在爆压至0.01 GPa的区间内都小于3%。结果表明，利用本文的逆特征线算法和遗传算法，理论上可以反演出压力范围较宽、准确性较高的爆轰产物状态方程。     

低冲击加载下JOB-9003炸药的反应阈值

发展了一种研究炸药反应阈值的实验方法和分析技术:采用火药炮发射飞片的加载技术产生低冲击加载压力,应用电磁粒子速度计测量JOB-9003炸药后界面与PMMA之间界面粒子速度。通过分析界面粒子速度曲线,得到了低冲击加载下炸药与PMMA之间的界面粒子速度历史,获得了入射压力与未反应和反应后的界面粒子速度之间的u_p-p关系。JOB-9003炸药在低冲击加载下的化学反应阈值和点火阈值分别为1.42、2.62 GPa。     

GPU-accelerated smoothed particle hydrodynamics modeling of jet formation and penetration capability of shaped charges

The prediction of the penetration of three-dimensional (3D) shaped charge into steel plates is a challenging task. In this paper, the smoothed particle hydrodynamics (SPH) method is applied to simulate the jet formation generated by the shaped charge detonation and its damage to steel plates. The Jones–Wilkins–Lee (JWL) equation of state (EOS), Tillotson EOS, and elastic–perfectly plastic constitutive model were incorporated into SPH for the modeling of explosive detonation and dynamic behavior of metal material. The compute unified device architecture (CUDA) parallel programming interface has been employed in SPH to improve the computational efficiency of SPH. Firstly, the constitutive models and EOSs are validated by 3D TNT slab detonation and aluminum–aluminum (Al–Al) high velocity impact. Then the jet formation of the shaped charge detonation and its penetration into the steel plates are investigated using the graphics processing unit (GPU)-accelerated SPH methodology. The numerical results of these test cases are compared against the published experimental data or analytical result, which shows that the GPU-accelerated SPH methodology is capable of tackling the 3D shaped charge detonation and penetration involving millions of particles with high computational efficiency.     

基于通用炸药状态方程分析飞板运动规律的特征线法

从小扰动波(马赫波)的物理概念出发,导出了不依赖流体状态方程表达形式的平面二维超声速定 常流的特征线方程;重新定义了以流体密度为单自变量的Prantl-Meyer函数,形成了求解平面二维超声速定 常流的封闭方程组。还利用这种通用物态方程的特征线差分解法,针对滑移爆轰驱动飞板运动问题构建了爆 轰产物流场内部和飞板边界特征线差分法格式。对TNT炸药和乳化炸药采用JWL状态方程和多方方程进 行了对比计算。结果表明,炸药爆轰对飞板的驱动能力与状态方程表示的炸药的做功能力是一致的。     

一种基于Hugoniot关系的爆轰产物等熵状态方程

对136组不同炸药的爆轰产物压力-粒子速度实验数据进行分段拟合,得到一个过C-J点的爆轰产物Hugoniot经验关系;对该经验关系进行Riemann积分,得到一个描述爆轰产物压力相对比容关系的爆轰产物等熵状态方程,该方程的参数仅为炸药的初始比容和C-J状态量,与传统经验等熵状态方程相比,不需要进行实验标定,因此可节约标定方程的实验成本和计算成本。为验证方程的合理性,采用该方程在压力相对比容平面上给出了Comp-B、HMX、PETN、ANFO、TNT以及LX-14炸药的爆轰产物等熵膨胀曲线,发现与采用JWL状态方程给出的相应炸药爆轰产物等熵膨胀曲线符合较好。     

钝感炸药的超压爆轰与冲击起爆过程数值模拟

采用Hybrid反应率结合修正的JWL方程,研究了LX-17、超细TATB等钝感炸药的冲击起爆（SDT）过程,并计算了爆轰波的对碰现象。结果表明,该方法计算钝感炸药的冲击起爆过程与实验数据符合较好;计算爆轰波对碰区的峰值压力提高了10％。     

ANALYTICAL MODEL OF CERAMIC/METAL ARMOR IMPACTED BY DEFORMABLE PROJECTILE

A new analytical model was established to describe the complex behavior of ceramic/metal armor under impact of deformable projectile by assuming some hypotheses. Three aspects were taken into account: the mushrooming deformation of the projectile, the fragment of ceramic tile and the formation and change of ceramic conoid and the deformation of the metal backup plate. Solving the set of equations, all the variables were obtained for the different impact velocities: the extent and particle velocity in rigid zone; the extent, cross-section area and particle velocity in plastic zone; the velocity and depth of penetration of projectile to the target; the reduction in volume and compressive strength of the fractured ceramic conoid; the displacement and movement velocity of the effective zone of backup plate. Agreement observed among analytical result, numerical simulation and experimental result confirms the validity of the model, suggesting the model developed can be a useful tool for ceramic/metal armor design.