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

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

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

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

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

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

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

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

炸药爆轰驱动水的初期过程

采用高速转镜分幅相机,利用阴影照相技术研究了PBX-01炸药爆轰产物在水中的初期演变过程, 获得了爆轰产物与水界面在不同时刻的图像。利用LS-DYNA程序对该实验进行了数值模拟,得到了水中冲 击波迹线、炸药爆轰产物与水界面随时间膨胀的位置、速度。数值计算结果与实验结果具有较好的一致性,所 得参数可为进一步设计改进实验提供参考。理论分析表明,水在爆轰产物作用下的最高瞬间温度约342K, 未到达水在实验压力条件下的相变汽化温度。     

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

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

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

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

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

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

爆轰产物状态方程的水下爆炸反演理论研究

由于水下爆炸过程中爆轰产物的信息以水中压缩波形式向外传递,本文旨在探讨如何利用水下爆炸试验数据确定爆轰产物的状态方程。相较于常规圆筒试验,水下爆炸试验具有装置简单成本低、装药尺寸限制少、测定压力范围更广的特点,更适用于大药量或非理想炸药的现场测试。本文从水中冲击波轨迹和波后压力时程曲线出发,发展了由冲击波及其波后流场还原水气界面的逆特征线算法,以及根据水气界面确定爆轰产物状态方程的遗传算法。与水下爆炸正演结果对比,发现逆特征线法可以较准确地还原水气界面的位置和压力参数,有效压力下限可达2 MPa,远低于圆筒试验的测试下限0.1GPa。根据水气界面的反演结果,遗传算法也能稳定地优化JWL方程参数,8种常用炸药的等熵衰减压力误差在爆压至0.01 GPa的区间内都小于3%。结果表明,利用本文的逆特征线算法和遗传算法,理论上可以反演出压力范围较宽、准确性较高的爆轰产物状态方程。     

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

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

柱状含铝炸药水下爆炸近场的特征线法研究

基于之前提出的一种含熵变项的特征线法,通过控制非等熵流中的能量释放来刻画铝粉燃烧的影响,结合简单Chapman-Jouguet模型和JWL-Miller状态方程,计算了柱形含铝炸药水下爆炸的近场参数。对比模拟结果与实验数据,发现这种特征线法可以较好地预测近场冲击波的传播迹线、爆轰产物的膨胀轨迹以及内部压缩波的反射过程。结果表明,这种特征线法可用于含铝炸药水下爆炸的近场计算,进一步可用于评估含铝炸药性能或计算水下能量输出。     

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

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

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

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

PBX-9404的化学反应速率方程及起爆特性

通过对炸药化学反应速率方程的分析,建立了基于Kim弹粘塑性球壳塌缩热点模型原理的三项式化学反应速率方程模型。运用遗传算法确定了反应速率方程相关参数,通过与Forest-Fire反应速率模型数值模拟结果的对比对所建模型的合理性进行了验证。将所建反应速率方程模型嵌入有限元程序对PBX-9404炸药起爆过程进行数值模拟,根据模拟结果分析了炸药遭受冲击加载后炸药孔隙率、颗粒尺寸等变化对炸药冲击起爆过程的影响。     

爆炸冲击波对墙体绕射效应的研究

采用试验与数值模拟相结合的方法,对爆炸冲击波作用于防爆墙的载荷与冲击波绕过 防爆墙的规律进行了研究. 采用压力传感器测压技术获得了防爆墙前后不同距离的压力曲线. 数值模拟基于ALE方法与JWL状态方程,试验结果与数值模拟基本相符,给出了冲击波作用 于防爆墙迎爆面上的载荷并研究了冲击波绕过防爆墙详细过程. 通过综合分析,得到了爆炸 冲击波绕流的内在机理为防爆墙的设计与加固提供依据和参考.     

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

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