JB-9014炸药超压爆轰产物的状态方程

根据P.K.Tang等提出的对炸药爆轰产物超压状态方程建模时只对JWL状态方程CJ等熵线中 高压指数项做修正的研究思路,首先给定超压状态下内能等熵线的修正项,再根据热力学定律对内能等熵线 求微分而得到沿压力等熵线的修正项。对JB-9014炸药超压爆轰冲击Hugoniot实验数据和声速实验数据同 时进行拟合,得到了3个JWL状态方程在超压爆轰状态下的修正项,并进行了分析与比较。得到超压修正项 的方法简单,3组超压修正项与P.K.Tang的修正项一样,都能很好地拟合超压Hugoniot数据。在实验数据 范围外,对超压状态下的声速-压力实验数据的拟合精度有所差别。     

滑移爆轰作用下飞板运动规律的研究──双 P-M 法

在爆炸载荷作用下飞板运动规律的研究是爆炸焊接中一个十分重要的问题。提出爆炸产物的侧向飞散过程可用普朗特－梅耶（Ｐｒａｎｄｔｌ－Ｍｅｙｅｒ）的绕流理论来描述，而且飞板的存在与弯曲对爆炸产物的稀疏作用也可以用这一理论来描述，即爆炸产物的真实运动状态是这两种运动的综合。进而推导出了描述飞板。     

PBX9502炸药爆轰产物的状态方程

采用vanderWaals等效单组分流体模型和Ross硬球微扰理论软球修正模型,计算了爆轰气相产 物的状态方程;用石墨相、金刚石相、类石墨液相和类金刚石液相4种相态描述凝聚成分,由Gibbs自由能最 小确定了不同状态下的凝聚产物相态。对爆轰产物混合系统采用自由能最小原理,通过化学平衡方程组求解 了炸药爆轰产物系统的平衡组分。使用该理论计算了高含碳炸药PBX9502Chapman-Jouguet(CJ)点的爆轰 参数,计算值与实验值符合很好;同时计算了3条等温线,并与Sesame库比较,发现温度超过1000 ℃时,计 算值与Sesame库的计算结果比较接近。在计算的5802K 等温线上发现了一个拐点,分析发现是由于在此 处游离态的碳发生了相变。     

在两层不同爆速炸药滑移爆轰作用下飞板运动的数值计算研究

本文研究用低爆速炸药复盖高爆速炸药在定滑移爆轰作用下飞板运动的计算方法,并对惯性复盖与炸药复盖的结果作出相互比较。     

双层楔形爆炸反应装甲飞板的运动规律

根据单层爆炸反应装甲(explosive reactive armor, ERA)飞板运动规律,建立了双层楔形ERA飞板的运动模型。利用该模型对双层楔形ERA各飞板的运动规律进行分析,由此计算出飞板之间相互作用的时间。采用X光脉冲闪光测试系统对相应的实验进行记录,结果表明：根据理论模型计算结的相互作用时间与实验测得数据吻合较好,且双层平行ERA较单层ERA的飞板相互作用时间提高了5倍,双层楔形ERA的飞板相互作用时间又较双层平行ERA的提升38%。     

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

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

     

滑移爆轰驱动下飞板运动姿态的连续电阻测试法

飞板运动姿态的测定是爆炸焊接机理研究的基础,针对传统电测方法存在干扰因素多、易产生弯曲波等缺陷,设计了一种适用于野外大当量下爆炸焊接飞板姿态实验的连续电阻测试方法。研制了3种不同结构的梯形支架型连续电阻探针元件,利用有限元程序分析了探针的导通压力和响应时间,在此基础上,对3种探针实施了爆炸焊接实验,实验结果表明:金属丝网型探针元件具有最优的导通效果,各段测试曲线光滑无毛刺。以该探针数据计算获得了待测飞板的运动姿态曲线,并与Richter简化模型下的近似计算公式结果进行了对比,两者基本一致。所述测试方法实现了炸药爆速和飞板变形曲线的连续、可靠和快速测量,为滑移爆轰驱动问题、爆轰产物状态方程等的研究提供了测试方法补充。     

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

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

组合状态方程在数值计算上的应用及分析

利用建立在爆轰产物实验等熵线数据上的组合状态方程以及常用的Jones-Wilkins-Lee (JWL)和Becker-Kistiakowsky-Wilson (BKW)爆轰产物状态方程进行了爆炸数值分析。状态方程数值计算的压力波形与实验波形的对比表明:组合状态方程和JWL状态方程计算的波形的到时、峰值、形状能够与实验吻合;BKW状态方程等熵式计算的数值波形峰值基本能够与实验吻合,但形状吻合不佳。     

工业炸药的原子经济性分析

运用原子经济性概念,引进了工业炸药氧化剂和可燃剂的能量贡献量和能量因子、爆炸产物的单位质量生成热和单位质量比容等指标,分析结果表明:工业炸药应选择能量贡献量、能量因子较大,爆炸产物的单位质量生成热和单位质量比容较大,氢、氧、碳元素含量较高的物质作为原材料,以利于获得性能优越的工业炸药。     

起爆高密度TATB炸药的飞片速度阈值

采用实验与数值模拟相合的方法研究了飞片起爆钝感TATB炸药的性能。设计了飞片可靠起爆钝感TATB炸药的起爆序列,利用全光纤位移干涉测速系统分别测出飞片可靠起爆和未起爆TATB炸药的速度,初步确定飞片起爆钝感TATB炸药的起爆速度阈值。采用DYNA2D程序对飞片起爆TATB炸药过程进行数值模拟,模拟结果与实验结果相符合。

     

爆轰驱动平板的抛射角方程

在运动坐标系中,将材料视作不可压缩流体,考虑流体微元间的相互挤压作用,由动量方程出发,推导出爆轰驱动二维平板过程中物质元的抛射角方程,并讨论了方程的解与矢量三角形几何关系的等价性。     

混合炸药爆速预报的新方法

为了提高混合炸药爆速尤其是含铝混合炸药爆速的计算精度,基于体积加权法思路构造一种新的混合炸药爆速计算公式,对近50种常见的混合炸药配方进行了计算,并与实验值进行对比,结果很好吻合,最大误差低于3%,平均误差低于1%,与BKW和Urizar常用方法相比,计算精度有显著提高。

     

强爆轰驱动飞片的数值模拟研究

在一维流体动力学程序中,使用了JWLT状态方程,对强爆轰驱动飞片的实验模型进行了数值模拟研究,得到了炸药爆轰产物的压力和密度分布,验证了强爆轰的存在,同时得到了二级飞片自由面的速度曲线,与实验测量结果符合较好,为解释实验结果和设计新的强爆轰实验提供了计算依据,同时验证了程序的正确性。     

一般物态方程形式下爆轰产物的一维等熵流动

研究了一般物态方程形式下一维平面和柱、球面对称CJ爆轰波后产物的等熵流动,给出Riemann不变量的表达式、Euler和Lagrange形式的自相似流动方程组,得到一个高精度近似解析解,在实验测量合理的范围内本文理论与实验基本符合。为了具体应用上述理论,给出四种普遍公认的物态方程(律、JWL、HOM和W.C.Davis近来提出的形式)的等熵线和声速曲线,以及较一般情形中产物或惰性介质沿等熵线的比内能表达式。四种等熵线差别不大,表明爆轰数值模拟中不但要应用产物的等熵线,还应注意产物能量状态方程是否适当。本文的理论计算方法,可作为产物一维流动数值模拟正确性的有效校核手段。