弱点火条件下炸药燃烧转爆轰的数值模拟

考察颗粒炸药从传导燃烧到对流燃烧再到爆轰的过程.对装填密度为85%的HMX颗粒炸药的燃烧转爆轰过程进行数值模拟,分析传导燃烧、对流燃烧和爆轰的发展过程.点火早期燃烧速度很低,火焰面在8.16 ms之内只前进了不到0.2 mm;形成对流燃烧之后燃烧速度快速增加,只用了0.1 ms就形成了速度为8 165 m·s^-1的稳定爆轰.当炸药颗粒直径或点火压力减小时,形成稳定爆轰所需的时间增加.     

HMX炸药爆轰产物物态方程研究

采用自主研制的考虑产物组分间化学反应的类CHEQ程序研究HMX炸药爆轰产物物态方程,给出初始密度为1.90 g·cm^-3的HMX炸药CJ点的爆速和爆压,计算值与实验值符合较好;计算CJ点的爆轰产物组成,并与BKW和LJD的计算结果进行比较分析;同时给出过CJ点的等熵膨胀线上压强、 γ和产物组分随比体积的变化规律,加深对爆轰产物作功能力的认识.     

以HMX为基的两种压装高密度炸药的燃烧转爆轰实验研究

 为研究以HMX为基的固体高能炸药的燃烧转爆轰性能,采用同轴电探针和压力传感器测试技术对常用的A、B两种压装高密度高能炸药开展燃烧转爆轰实验,研究装药组分和约束条件对压装高密度炸药燃烧转爆轰性能的影响。实验结果表明：这两种压装高密度炸药难以发生燃烧转爆轰;在强约束条件下（45号钢,内径25.4 mm、外径65 mm、长度600 mm）,A压装炸药（HMX质量分数为95%,密度为1.86 g/cm³）基本实现了燃烧转爆轰,爆轰诱导距离约为545 mm;在相同的实验条件下,A压装炸药比B压装炸药（HMX质量分数为87%,密度为1.84 g/cm³）更易于发生燃烧转爆轰,即A压装炸药的安定性相对较差。     

CE/SE方法数值模拟炸药粉尘爆轰

采用CE/SE方法数值模拟悬浮在空气中的RDX炸药粉尘的两相爆轰过程.炸药颗粒在爆轰波阵面后的高温高速气流中加速并升温,释放能量支持爆轰波传播.数值模拟爆轰波管中的粉尘爆轰,得到爆轰波流场中的物理量分布,确定爆轰参数,数值结果与文献符合较好.数值模拟复杂通道中的炸药粉尘爆轰,预测了爆轰波的发展和传播过程以及爆轰波后的流场演化.数值结果表明CE/SE方法能成功模拟气体-固体两相爆轰,为粉尘爆轰的研究提供了新的数值预测手段.     

炸药枪击试验和数值模拟研究

HL-10炸药是一种以RDX为基的含铝炸药,为了研究该炸药在子弹或金属破片撞击作用下的安全性,利用12.7 mm机枪法对钢壳包覆的柱形HL-10装药进行了枪击试验,试验结果表明,炸药没有发生燃烧或爆轰现象,由此可定义该炸药的枪击感度试验反应等级为1级。建立了炸药枪击试验的计算模型和数值计算方法,对子弹撞击和穿透炸药过程进行了三维数值模拟计算,计算结果与试验结果相符,分析了子弹速度对HL-10炸药枪击感度的影响,其结果可为炸药安全性评价分析提供理论根据。为了进一步验证计算模型和方法,对美国的PBX-9404炸药的枪击作用过程进行了计算分析,结果表明,PBX-9404炸药在枪击作用下发生了完全爆轰反应,与Neff等人的试验结果相吻合。     

高密度B炸药的燃烧转爆轰实验

在相同的实验条件下研究了3种不同密度固体压装B炸药（TNT，RDX为40／60）的燃烧转爆轰性能。3种B炸药的密度分别为1．597，1．654g／cm^3与1．681g／cm^3。用时间间隔记录仪、电离式电探针系统记录燃烧波、压缩波或爆轰波到达的时间；用程控电荷放大器、数字示波器及压力传感器记录DDT管不同位置处的压力历程。起爆器、时间间隔记录仪及示波器均由同步机触发。DDT管材料为45号钢，内径φ20mm，外径φ64mm，长度500mm；点火药为小粒黑，质量1．1～1．2g；所使用的电探针为同轴电探针（芯线为0．9mm的铜漆包线，     

DDT管材料对颗粒状RDX床燃烧转爆轰（DDT）影响的实验研究

 采用盖帽探针、离子探针实验研究了点火药点燃时,DDT管材料变化（钢和铝）对颗粒状RDX床的燃烧转爆轰的影响。     

炸药爆轰的连续介质本构模型和数值计算方法

假定炸药和爆轰产物处于局部热力学平衡状态, 即它们的压力和温度相同, 利用热力学基本关系建立炸药爆轰过程的连续介质本构模型的一般理论框架. 在此框架下, 炸药爆轰本构模型由一组常微分方程构成, 包括炸药和爆轰产物的状态方程、简单混合法则、化学反应速率方程和能量守恒方程, 易于由成熟的计算方法如梯形法等进行求解. 一组广义Maxwell型非线性固体本构形式的微分方程描述了压力和温度随时间的演化速率与应变率和化学反应速率的关系, 借助简单混合物理论, 其中的系数由炸药和爆轰产物的材料参数确定. 未反应的炸药和爆轰产物采用JWL状态方程, 化学反应率方程采用Lee-Tarver点火-燃烧二项式模型, 模拟PBX-9404炸药的一维冲击波起爆过程和爆轰波传播过程. 计算结果表明了本文给出的本构模型和相应计算方法的有效性. 关键词： 炸药爆轰 本构模型 化学反应率方程 数值模拟     

炸药颗粒压制成型数值模拟

 分析炸药压药过程中细观的力学行为,能够为改进压药工艺和提高炸药元件质量提供理论依据。建立了模压条件下炸药颗粒压制成型的计算模型。模型中炸药颗粒被认为是直径相同的球形颗粒,并按一定规律排列。利用非线性有限元计算方法,对炸药颗粒压制成型过程进行了数值模拟计算,分析了压制过程中炸药颗粒变形、受力和温度变化情况。结果表明：药粒在压缩中存在运动和变形两个阶段。在药粒运动阶段,应力集中主要出现在颗粒与约束面的接触部分;药粒进入了塑性变形后,药粒内部压力迅速升高且压力趋于一致。压缩过程中药粒温度升高,药床接近密实状态时,药床中心处药粒温度最高。     

钢筋混凝土靶板在弹丸冲击及爆炸载荷下响应的数值模拟

 采用HJC混凝土损伤本构模型及LS-DYNA的流固耦合算法,分别对钢筋混凝土靶板在弹丸冲击和爆炸载荷作用下的响应进行了有限元数值模拟,其中模拟参数由实验数据拟合重新获取。将模拟结果与实验结果和经验公式进行对比分析,结果表明：数值模拟再现了弹体贯穿靶板过程中的开坑、隧道及漏斗碎裂区,计算得到的弹体弹道极限及残余速度与实验数据吻合较好;此外,数值模拟也很好地再现了炸药爆炸后冲击波的传播过程以及爆炸载荷作用下混凝土的破坏情况,模拟结果与实验现象具有良好的一致性。     

药柱等静压实验和数值模拟计算

 进行了炸药柱等静压实验,用热电偶测量了药柱压制成型过程中的内部温度,建立了药柱在等静压下的热力耦合模型。对炸药柱保压阶段进行了数值模拟计算,得到等静压条件下药柱内部压力和温度变化,分析了药柱形变,压力和温度分布。结果表明：在等静压下药柱各表面均有明显的向里凹陷的现象,在受力过程中炸药内部温度升高,并存在温差,内部压力分布不均匀。     

高速水射流与凝聚炸药相互作用的数值模拟

 利用LS-DYNA3D软件并采用Lagrange算法，对500 m/s的高速水射流与凝聚炸药的相互作用过程进行了数值模拟；对不同时刻凝聚炸药中的压力分布以及不同径向、轴向处的压力历史进行了分析；确定了所选凝聚炸药的声速。数值模拟的结果与理论分析的结果之间具有很好的一致性。模拟结果表明，在高速水射流与凝聚炸药相互作用的过程中，凝聚炸药药柱截面上的压力随径向距离的增加而降低；药柱中的最大压力出现在水射流头部附近或水射流头部的前方。     

A Possible Interpretation of Burst-Like Characteristics of Explosive Events

Explosive events have been observed to occur consecutively in bursts at intermittent locations along the boundary near the opposite polarity. The aim of the present paper is to explore a possible mechanism to interpret this burst-like characteristic of explosive events. The 2D magnetohydrodynamic (MHD) numerical simulations with resistivity have been carried out to reproduce the intermittent spatial-temporal magnetic reconnection events taking place along the long, compressible current sheet. The observed density enhancements in previously published results have been verified to be associated to magnetic reconnection sites. Late observational evidences, which support present attempts, have also been found, at least in morphological evolution of the consecutive explosive events.     

水中障碍物爆炸毁伤效应的二维数值模拟

 建立了水下障碍物的爆炸毁伤效应的物理模型和计算模型。建立了守恒方程组、混凝土的动态弹塑性本构关系和破坏准则,并确定了其在大应变、高应变率及高压条件下的的屈服准则。基于多物质流体的Euler算法,用面向对象的C++语言自行编制了NM-MMIC通用多物质二维流体弹塑性程序,对不同炸点位置、不同炸药量对轨条砦爆炸破坏过程进行数值模拟,并采用自行开发的后处理软件VISC2D实现了可视化,定量地反应了障碍物迎爆面所承受的压力变化过程,确定了轨条砦毁伤的最佳药量,并对毁伤程度进行了评价,得到了爆炸破障的毁伤规律。研究结果对设计和评价反水中障碍物的弹药具有一定的参考价值。