PETN、RDX和HMX炸药爆轰参数的数值模拟

用吉布斯自由能最小原理,通过解化学平衡方程组,求解PETN、RDX和HMX炸药爆轰产物系统的平衡组分,计算结果与用BKW和LJD方法计算的结果相近。用自编的程序从碳的石墨相、金刚石相、类石墨液相和类金刚石液相4种相态中确定出炸药爆轰产物中游离碳更可能存在的相态,并用此相态计算碳的吉布斯自由能。以WCA状态方程作为爆轰气相产物的物态方程,对PETN、RDX和HMX炸药爆轰参数作了预言,爆轰CJ点的爆速、爆压和爆温的计算结果与实验值吻合得很好。     

遗传算法在炸药爆轰参数计算中的应用

为解决爆轰参数计算中自变量取值范围要求严格以及收敛性差等问题,在最小自由能法的基础上,引入遗传算法建立了炸药爆轰参数的计算方法,并利用典型单质炸药和混合炸药的实验结果对计算结果进行了验证。结果表明:采用该方法计算得到的单质和混合炸药的爆速、爆压与实验测试结果基本一致,误差在5%以内,满足炸药性能预估的要求;该方法人工干预弱,只需一次性确定少数主要组分物质的量的变化范围,可适应于多配方优化计算。     

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

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

     

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

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

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

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

用 KHT 状态方程计算炸药爆轰参数

调整了部分气态爆轰产物的ＫＨＴ状态方程参数及游离Ｃ的Ｃｏｗａｎ状态方程参数，应用新参数计算了典型单质炸药、混合炸药和燃料空气炸药的爆轰性能。结果表明，应用新参数得出的炸药爆轰性能更加接近实验值；对于含铝炸药，与文献相比，计算结果不仅与实验相符，而且体现了铝为高能添加剂的特点。     

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

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

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

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

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

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

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

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

松弛方法在计算凝聚炸药爆轰问题中的应用

利用松弛近似,将非线性的凝聚炸药爆轰控制方程转化为线性的松弛方程组,并采用五阶WENO格式和五阶线性多步显隐格式对线性松弛方程组进行空间方向和时间方向的离散,由此建立具有高精度和高分辨率性质的计算凝聚炸药爆轰的松弛方法。建立的松弛方法可以避免求解Riemann问题及计算非线性通量的Jacobi矩阵,同时无需分裂处理反应源项。通过对凝聚炸药的平面一维定常爆轰波结构及球面一维聚心、散心爆轰起爆和传播过程的数值模拟,验证了所建立的松弛方法能够很好地计算凝聚炸药爆轰问题。     

挤注炸药的爆轰性能

设计了4种以HMX为主炸药的挤注炸药(ECX)炸药配方,采用25 mm和50 mm 2种装药直径的圆筒试验分别对添加含能增塑剂和惰性增塑剂的两种挤注炸药进行做功能力研究。研究结果表明:含能增塑剂FM-1对ECX有显著的能量贡献;在相同主炸药含量（含HMX为88％）下,含FM-1为9.5％的ECX-05较含惰性增塑剂的ECX-02比动能提高22％～23％,优于国外同类型ECX（EX-08-EL）的水平,其加速金属能力与含HMX95.5％的压装PBX炸药（LX-14）相近。     

JOB-9003和JB-9014炸药平面爆轰驱动飞片的对比研究

采用双灵敏度激光速度干涉仪(VISAR)测量了JOB-9003和JB-9014两种炸药平面一维爆轰驱动飞片的速度历史,比较了二者的做功能力;并用DYNA2D程序进行了数值模拟。在有JO-9159炸药的条件下,实验和计算结果表明,JB-9014炸药的做功能力可以达到JOB-9003炸药的75.0%;无JO-9159炸药条件下的计算结果表明,JB-9014炸药的做功能力约为JOB-9003炸药的71.2%     

可燃液体爆燃特性及其抑制实验

设计了一套可燃液体爆炸特性实验装置,利用该装置和立式爆轰管对RP-5油料、RP-3油料及工业酒精的爆炸特性、1301惰性气体对这3种燃料的抑制进行了研究。结果表明:RP-5油料、RP-3油料及工业酒精爆炸的体积分数范围分别为1.53%～7.73%、0.82%～7.17%及 3.38%～18.25%;酒精云雾爆轰的临界起爆能为2.11 MJ/m2、爆速和爆压分别为1 609 m/s 、1 480 kPa,爆轰波传播的胞格宽度为14.5 mm,长度为16.2 mm。1301惰性气体对RP-5油料、RP-3油料及工业酒精的最小惰化体积分数分别为6.75%、6.8%及 5.56%;二氧化碳和氮气对RP-3油料的最小惰化体积分数分别为45%和49%;1301惰性气体对油料爆炸抑制效果明显好于二氧化碳与氮气。