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

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

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

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

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状态方程和反应速率方程可行,金属平板驱动试验的计算结果与试验结果吻合。     

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

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

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

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

含铝炸药水下滑移爆轰实验研究

为了研究含铝炸药水中爆炸近场范围内沿装药径向的冲击波特性及爆轰产物的膨胀规律,设计了水下滑移爆轰实验装置。采用高速扫描相机和阴影照相技术记录了一种RDX基含铝炸药的水下滑移爆轰过程,获得了清晰的水下滑移爆轰过程的扫描图像。通过对图像进行分析,得到了水中爆炸近场范围内沿装药径向冲击波的传播迹线及爆轰产物气泡的膨胀迹线,进而分析出冲击波传播速度、阵面压力及爆轰产物气泡的膨胀位移等参数的变化规律。此外,根据气泡的膨胀迹线,标定了该含铝炸药爆轰产物的JWL状态方程参数,并将其与Φ50mm圆筒实验确定的参数值在p-V图上进行了对比。结果显示,两者偏差较小,证实了该方法的可行性。     

含铝炸药爆轰产物JWL状态方程研究

对ＲＯＴ ９０１含铝炸药进行了爆轰产物ＪＷＬ状态方程研究。对ＲＯＴ ９０１炸药进行了２５ｍｍ和５０ｍｍ两种装药直径的圆筒试验,测试了圆筒膨胀过程的Ｒ ｔ关系。由冲击Ｈｕｇｏｎｉｏｔ关系和参考的ＰＯＰ关系拟合了ＲＯＴ ９０１炸药的反应速率。对两种装药直径的圆筒试验进行了综合分析和二维数值模拟计算,综合评估了ＲＯＴ ９０１的作功能力和圆筒试验的特性。通过对两种装药直径圆筒试验圆筒膨胀过程的解析分析计算和二维数值模拟,确定得到了ＲＯＴ ９０１含铝炸药爆轰产物ＪＷＬ状态方程参数。     

基于圆筒实验的RDX/Al炸药反应进程

对RDX炸药和2种铝粉质量分数分别为15%、30%的RDX基含铝炸药进行?50mm圆筒实验,研究铝粉含量对炸药做功能力的影响,根据格尼公式分析铝粉与爆轰产物的反应进程。结果表明:在圆筒实验记录的时间范围内,铝粉质量分数为15%的含铝炸药做功能力最强,RDX炸药次之,铝粉质量分数为30%炸药做功能力最弱;34μs时刻,铝粉质量分数为15%的炸药,铝粉的反应度为0.49,而铝粉质量分数为30%炸药铝粉的反应度仅为0.21,含铝炸药中铝粉的反应时间在50~200μs之间。     

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

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

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

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

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

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

小隔板试验冲击起爆和熄灭流场的实验与Lagrange分析

根据小隔板试验的拉氏传感器测试结果,用拉氏分析方法求得了粗颗粒包复TNT炸药中的冲击波成长过程和熄灭过程的流场分布,用实验结果确定了点火增长模型中的待定参数。利用拉氏分析方法时,对未反应炸药和已反应产物的状态方程都采用JWL状态方程,使确定的参数可以直接用于有限元程序LSDYAN计算。通过拉氏分析得到了各迹线上的反应度,结果表明,冲击波在熄灭过程中仍能引起炸药的局部化学反应。     

圆筒试验JWL状态方程参数的贝叶斯标定

在研究数值模拟的输入参数引入的不确定性时,通常需要人为给定每个输入参数的概率分布,且输入参数概率分布的选择可能会对分析结果产生直接影响。利用贝叶斯方法标定了圆筒试验JWL状态方程参数,得到了标定参数的估计值和后验分布,并研究了不同统计模型假设对标定参数的估计值和后验分布的影响。贝叶斯后验分布融合了基准试验的试验数据的信息,因此将其作为不确定度量化分析时输入参数的初始概率分布,可以尽量减少分布选择引入的认知不确定性。