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

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

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

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

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

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

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

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

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

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

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

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

爆轰产物组成和爆轰参数计算方法的理论研究

为实现爆轰产物组成和爆轰参数的计算,采用拉格朗日乘数法和牛顿迭代的方法预测爆轰产物组成,利用BKW状态方程预测爆轰参数,在0~600 GPa和300~15 000 K压力温度范围内选取金刚石作为碳的生成相;对爆轰产物系统采用最小自由能原理,结合牛顿迭代法求解爆轰产物的化学平衡方程组;对BKW状态方程参数提出修订,取α=0.5,β=0.298,θ=6 620,κ=9.50;采用自编程序实现计算过程。使用此方法和Hugoniot关系计算密度为1.77 g/cm3的PETN爆轰CJ点爆轰参数验证计算精度,结果显示计算与实验结果的误差均小于1%。利用此方法结合Hugoniot关系预测出爆轰CJ点的产物密度为2.43 g/cm3。     

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

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

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

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

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

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

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

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

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

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

内爆炸爆轰产物运动过程二维数值模拟研究(一)——密闭空间

为了研究密闭空间内爆炸爆轰产物运动规律,采用五阶WENO有限差分格式,通过耦合求解爆轰产物质量分数输运方程与可压缩欧拉方程,自主编写了密闭空间内炸药爆炸过程二维数值计算程序。利用所开发的程序数值探讨了密闭空间内爆轰产物的传播过程及炸药形状、产物状态方程对爆轰产物运动的影响规律。研究表明:炸药形状的差异不仅改变了爆炸初期爆轰产物的膨胀形态,同时也改变了爆轰产物与周围空气的后期混合状态。状态方程对产物运动过程同样具有一定的影响,JWL状态方程中高压项可使得爆炸初期产物具有更大的膨胀能力。本文的研究可为进一步探讨考虑燃烧效应后爆轰产物扩张过程及毁伤评估提供一定的参考和指导。     

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

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

用VLW 状态方程计算水的冲击Hugoniot曲线

采用L-J(Lennard-Jones12-6)势函数,结合VLW 状态方程,计算了水的冲击压缩Hugoniot曲线。 将理论计算结果与实验结果进行了比较,当势函数参数/K=120、b0=30.42时,计算值与实验符合较好。与 采用BKW 状态方程计算结果相比,VLW 状态方程能够更准确地描述水的冲击Hugoniot曲线。     

神经网络状态方程在强爆炸冲击波数值模拟中的应用

强爆炸数值模拟的主要挑战在于如何准确地描述爆炸产物状态方程。利用BP神经网络和强爆炸产物状态数据对神经网络产物状态方程进行训练,并将得到的状态方程植入自编的一维球对称数值模拟程序,对强爆炸冲击波参数进行了计算。结果显示,计算得到的冲击波峰值超压、冲击波到时、正压时间与标准值吻合较好,证明将神经网络状态方程应用于强爆炸冲击波数值模拟是可行的。研究结果对确定强爆炸数值模拟方法具有很好的借鉴意义。     

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

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

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

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

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

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

梯恩梯/黑索今(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炸药的结果有一定差别。使用得到的这两种速率函数进行了一维数值模拟计算,计算结果同实验基本一致。