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

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

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

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

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参数确定方法合理可靠。     

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

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

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

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

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

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

梯恩梯/黑索今(35/65)炸药的反应速率函数

本文采用Lagrange分析方法和冲击引爆实验数据研究了TNT/RDX(35/65)炸药的反应速率函数。在冲击压力4.69GPa(持续时间约1s)的作用下,采用F4/203A型锰铜压力计测量了样品内六个不同Lagrange位置上的压力历史。以此作为输入量,采用Lagrange分析方法计算了反应流场中的速度、密度和比内能历史。根据这些结果,再使用HOM、JWL和律爆轰产物状态方程计算了反应进程和反应速率。在此基础上采用最小二乘法拟合了Forest Fire和Cochran点火-成长速率函数中的待定系数。考察了不同产物状态方程对这些系数的影响。用HOM和JWL方程得到的速率函数,基本一致,但用律方程得到的速率比前两者低,这是由于在反应区中有明显熵增的缘故。在较高压力下,本文给出的Forest Fire速率同Mader依据POP关系和反应Hugoniot算出的Comp.B炸药的结果有一定差别。使用得到的这两种速率函数进行了一维数值模拟计算,计算结果同实验基本一致。     

金属约束下定常非理想爆轰的理论研究

为了对金属约束条件下的定常非理想爆轰进行理论研究,对未反应炸药和爆轰产物采用JWL形 式状态方程,对金属采用p(,T)形式状态方程。采用过爆轰前导冲击波的流线偏转角在影响域内沿高度线 性变化的假设,并且由未反应炸药和金属的冲击波极曲线的交点确定爆轰冲击边缘角,则可从未反应炸药斜 冲击波极曲线关系式求出爆轰前导冲击波阵面的形状。采用爆轰前导冲击波阵面之后的流线是直线的假设, 则爆轰流动控制方程由偏微分方程变为沿流线的常微分方程,沿着所有流线求解便给出爆轰化学反应区内 声速线与化学反应结束线的位置。理论分析同时给出约束金属折射冲击波后面的流体的流动状态。理论结 果给出的爆轰化学反应区结构特征和约束金属内的流动状态特征与高精度数值模拟的结果符合良好,说明 本文中给出的理论方法具有良好的合理性和适用性。     

有氧化剂(AP)含铝炸药的爆轰性能

对有氧化剂含铝炸药(RDX/AP/Al/粘合剂=20/43/25/12,下称含铝炸药)爆轰反应的点火增长模型进行研究。用VLW状态方程方法计算了含铝炸药爆轰产物JWL状态方程;用激光速度干涉仪（VISAR）测量炸药/窗口界面粒子速度和炸药驱动金属平板自由表面速度,对试验进行了数值模拟计算,拟合了含铝炸药的反应速率方程。研究结果表明,用VLW状态方程方法和炸药/窗口界面粒子速度确定JWL状态方程和反应速率方程可行,金属平板驱动试验的计算结果与试验结果吻合。     

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

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

RHT-902和Octol炸药爆轰产物JWL状态方程研究

本文对RHT-902和Octol两种炸药做了标准圆筒试验,运用能量守恒原理分析处理试验结果,用解析关系式确定了这两种炸药爆轰产物的JWL状态方程参数。这些参数与国外发表的相应参数基本一致。     

可燃介质中飞行圆球诱导斜爆轰的流场结构

基于带化学反应的二维轴对称Euler方程,利用带有Superbee限制函数的波传播算法,对氢/氧/ 氮预混气中飞行圆球诱导斜爆轰进行了数值模拟。结果表明,在达姆科勒数Da略大于临界情形时,圆球诱导 的驻定爆轰是一个由强过驱斜爆轰、弱过驱斜爆轰、反应激波和惰性激波组合而成的复合结构。波后的圆球 绕流流场内存在2个亚音速区和1个超音速区。在圆球背风表面,还形成了第2道斜激波。     

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

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