基于非平衡多相模型的含铝炸药爆速研究

通过对爆轰反应区中铝粉压力、温度弛豫时间的计算,发现爆轰产物压力是平衡的,温度是不平衡 的。对传统CJ模型进行了改进,考虑爆轰产物的多相性和产物温度间的非平衡性,提出了一种新的计算含 铝炸药爆速的模型。采用该模型对几种含铝炸药的爆速进行了计算,并与已有的实验数据进行了对比,计算 结果优于传统CJ模型,与实验的误差在2%以内。     

CL-20混合炸药的爆轰波结构

基于爆轰数值模拟计算,分析了CL-20混合炸药爆轰反应的特征,设计了炸药与窗口的界面粒子速度测量实验装置;采用激光干涉法,测量了C-1炸药(CL-20/粘合剂/94/6)与窗口的界面粒子速度; 运用先求导、再分段拟合的方法,对界面粒子速度随时间的变化曲线进行了数据处理,确定了炸药爆轰CJ点对应的时间位置;根据CJ点对应的粒子速度,计算获得了炸药的爆轰反应区宽度和CJ爆轰压力。结果显示:密度为1.943 g/cm3的C-1炸药的爆轰反应时间为38 ns,CJ压力为34.2 GPa。     

客车车身骨架动应力研究的现状与未来

综述了多体系统动态模型的发展及动态应力历程求解的进展,给出了国际上的最新研究成果,并以CJ6121GCHK型客车车骨架为例,给出了计算其动应力历程的步骤。     

基于LADM模型的爆轰参数计算

炸药爆轰以高低两种速度稳定传播。经典的ZND模型没有对应于低速爆轰的数学解,无法计算低速爆轰参数。而LADM爆轰模型(Least Action Detonation Model)有四组数学解析解。其中一组解对应于高速爆轰,可取代有争议的CJ假设;另一组解对应于低速爆轰,可计算低速爆轰速度。本文用这两组解计算了四种常用炸药的高/低速爆轰参数,不需要CJ假设等任何假设。结果表明:高速爆轰速度的实测值与计算值的相对误差范围是0.069%~1.38%;炸药的绝热指数γ1为1.239~1.312,接近多原子气体的理论值1.286;爆轰产物的绝热指数γ2为1.425~1.462,略大于双原子气体的理论值1.400。     

非理想爆轰的推飞片效率

本文采用一定形式的状态方程,模拟爆轰波后稀疏波区内介质的放能,初步解析地讨论了非理想爆轰波的推飞片效率问题,并与相应CJ爆轰波的情况进行了对比。由本文的结果可见:(1)对CJ点之后释放能量的非理想爆轰波,波后压力分布与瞬时释放能量的多少和CJ点之后释放能量的多少和快慢有关。(2)非理想爆轰的推飞片效率,除炸药与飞片质量比之外,不仅依赖于CJ点之前瞬时释放出来的能量Q,还依赖于CJ之后释放的能量q和标志释放快慢的参数n。(3)非理想爆轰波、包括弱爆轰的推飞片效率,对所讨论的参数与CJ爆轰相似,均存在一个最佳的选择。     

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

为实现爆轰产物组成和爆轰参数的计算,采用拉格朗日乘数法和牛顿迭代的方法预测爆轰产物组成,利用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。     

用三项式点火增长反应速率模型计算JB-9014的反应区

在一维流体动力学编码SSS程序中,利用三项式点火增长反应速率模型对钝感炸药JB-9014的反应区进行了计算。计算中,未反应炸药采用固态HOM状态方程,产物采用气态JWL状态方程,计算得到了钝感炸药JB-9014化学反应区的峰值压力、CJ压力、反应区宽度和反应时间等参数,计算结果与实验结果符合较好。分析表明,三项式点火增长反应速率模型可用于研究钝感炸药的反应区结构。     

基于PDV的JOB-9003炸药爆轰反应区测量

炸药的反应区数据对爆轰过程的精密建模具有重要意义,为了得到JOB-9003炸药的反应区信息,采用光子多普勒测速仪(PDV)对JOB-9003炸药的爆轰反应区进行了实验研究。实验中利用火炮发射高速蓝宝石飞片冲击起爆被测炸药,在炸药后表面安装镀膜氟化锂(Li F)窗口测量炸药一维稳态爆轰时的界面粒子速度,测试过程的时间分辨率小于1 ns,测速相对不确定度小于2%。通过读取界面粒子速度时程曲线的拐点来确定CJ点,根据阻抗匹配公式计算炸药的CJ压力。研究结果表明,JOB-9003炸药界面粒子速度时程曲线上存在较为明显的拐点,JOB-9003炸药的化学反应时间为(11±2)ns,对应的化学反应区宽度为(0.075±0.014)mm,JOB-9003炸药的CJ爆压为(35.6±0.9)GPa,冯诺依曼(Von Neumann)峰处的压力为(47.9±1.2)GPa。     

简化CH4/O2/N2基元反应模型在爆轰模拟中的应用

基于改进的全耦合TVD格式和简化的基元化学反应模型,建立了二维可燃预混合气爆轰波传播过程的模型并编制了相应的CFD程序,对CH4／O2／N2预混气爆轰波的形成和传播进行了数值模拟。使用不同的计算初始条件来考察模型的计算模拟能力,计算结果表明,该模型能够精确捕捉爆轰波及其波后的反射过程,同时,也能够预测不同计算初始条件下爆轰波的传播轨迹。计算结果与相同条件下CJ爆轰理论模型的计算结果进行了比较,两者符合较好,这表明简化的CH4／O2／N2机理能够较好地应用于可燃气爆轰过程的数值计算中。     

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

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

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

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

数字相关方法在测量CJ6120LC6HK型客车前后轮轴心距离的应用研究

本文从计算机图像处理技术中的数字相关方法出发，在测量领域中与实际工程背景相结合，提出一种新的量程大，精度高，非接触式的方法，用来测量CJ6120LC6HK型客车前后轮轴心距离，本文在计算机图像处理技术中对数字相关方法在实际测量应用中所涉及到的相关函数，测量精度和测量误差进行了分析和研究，在用1m标尺精确定位的基础上，对CJ6120LC6HK型客车前后轮轴心距离进行了实际测量实验，得到了比较理想的结果，为大距离精密实际测量提供了一个良好的方法。     

根据爆轰产物的微观描述预估凝聚炸药的爆轰压力

40年代以来,从凝聚炸药的化学成分和密度直接计算它的爆轰性质有几种不同的方法,它们的共同点都是采用稳定流动的守恒定律和在爆轰的不连续面上的CJ条件,不同点在于对爆轰产物这样的非常稠密的气体和固体的混合物用什么样的状态方程来描述,本文综述几种状态方程的特点和它们所预估的爆轰压力值之间的差异。      

快速热点火熔奥梯铝炸药燃烧转爆轰实验研究

以熔铸型含铝混合炸药熔奥梯铝为对象,研究铸装含铝混合炸药快速热点火后的燃烧转爆轰特性。建立了快速热点火燃烧转爆轰实验平台,由实验装置(加热装置、约束钢管、炸药)、压力测试系统、光纤测速系统组成;加热装置加热15 mm厚45钢钢板,峰值温度大于1 100℃,温升速率为85～95℃/s。开展了快速热点火带壳熔奥梯铝炸药燃烧转爆轰实验,由加热装置加热约束钢管内熔奥梯铝炸药,炸药化学反应阵面压力和传播速度分别由压电性高压压力传感器和光纤探针测定;实测阵面压力约1 GPa,传播速度最大约2 600 m/s。由光纤数据获得炸药化学反应阵面传播轨迹,通过特征线方法获得冲击形成点,半定量给出冲击形成距离大于850 mm;并比较了管体破片质量实测值与炸药完全爆轰时破片平均质量计算值,实测值远小于计算值。综合实测化学反应阵面传播速度和压力、冲击形成距离分析、破片质量比较,可确定熔奥梯铝炸药没有发生完全爆轰,其化学反应状态为爆燃。另外,采用Adams和Pack模型、CJ燃烧模型,都能够半定量的预估冲击形成距离和燃烧波后压力,为实验设计提供依据,但CJ燃烧模型的计算结果更接近于实测值。     

含铝炸药能量释放的简化模型

为了在水下爆炸效应中反映出非理想爆轰特性的影响,建立了含铝炸药非理想爆轰能量释放的简化模型。该模型以CJ爆轰理论和二次反应理论为基础,把含铝炸药化学反应划分为快速反应和慢速反应,以释放的化学能和慢反应速率常数作为非理想特征参数,并应用于一维数值模拟。计算结果与基本实验结果一致,冲击波峰值的计算误差不大于10％,衰减时间常数的误差小于5％,冲击波能与实验值也具有良好的一致性。简化模型合理地描述了含铝炸药非理想爆轰的主要过程及非理想特性,可应用于含铝炸药的设计和爆炸效应的分析。     

惰性介质对钝感炸药爆轰的约束效应

首先用改进冲击波极曲线理论分析惰性介质对钝感炸药爆轰的约束作用类型。改进冲击波极曲线基于爆轰ZND模型建立在前导冲击波上,并且未反应炸药采用JWL状态方程,惰性介质采用p（ρ,T）形式状态方程。理论考察声速小于炸药CJ爆速且压缩性不同的6种典型惰性介质约束情况。然后用带三项式Lee-Tarver化学反应率的二维Lagrange流体力学方法数值模拟考察约束相互作用。数值考察约束介质的影响因素有:压缩性、厚度、典型双层介质组合约束。从数值结果看出,由介质压缩性的不同给出的约束作用方式共7种:其中6种出现在介质声速小于炸药CJ爆速条件下,可运用冲击波极曲线理论;另外一种出现在介质声速大于炸药CJ爆速条件下,不能使用冲击波极曲线理论。同时,介质厚度、双层介质组合方式也能够影响爆轰前导冲击波阵面形状以及爆轰化学反应流动状态。     

气相爆轰的胞格结构

气体爆轰现象是Berthelot等人在1881年发现的。Chapman(1899年)与Jouguet(1917年)首先提出了爆轰波的一维流体动力学理论《称CJ理论),他们把爆轰波阵面简化为强间断,假定流动是一维定常的,利用流体力学的Euler方程组解出了波阵面上的爆轰参数。     

爆轰波与激波对撞的实验研究

对乙炔氧气混合气体中爆轰波与激波的正面对撞现象的实验研究是以高速摄影获取两波对撞的x-t纹影图,以烟迹板记录对撞中的爆轰胞格图案,并基于激波理论和经典CJ爆轰理论求解了两波对撞的稳态解并探寻其规律. 研究发现透射波系包括一道激波和爆轰波,以及紧随爆轰波后的稀疏波区,这一结果对应于一维理论分析中的CJ解. 透射波系基本不受初始压强影响;初始温度也只成比例地改变流场整体速度,温度越高,速度越快;对波系起实质影响作用的是入射激波强度,激波越强,则整个透射流场呈现偏向激波的趋势;理论分析还指出,稀疏波区的出现不可避免,当激波强度趋于声波稀疏波区趋于消失,激波越强则疏波区趋于扩大. 两波对撞存在一个有限的转变阶段,透射爆轰首先减缓,接着迅速迸发为过驱爆轰,然后再逐渐平衡为CJ爆轰. 对于强不稳定的燃气,对撞后爆轰波在空间上的发展极不均衡,一些区域发生火焰面与诱导激波的严重脱离,随后的火焰面失稳发展为诱导激波区内的爆轰波,实验观察到了这种爆轰在烟迹板上留下的极为精细的迹线.      

5.12汶川大地震诱发大型崩滑灾害动力特征初探

对乙炔氧气混合气体中爆轰波与激波的正面对撞现象的实验研究是以高速摄影获 取两波对撞的x-t纹影图,以烟迹板记录对撞中的爆轰胞格图案,并基于激波 理论和经典CJ爆轰理论求解了两波对撞的稳态解并探寻其规律. 研究发现透射波系包括一道激波和爆轰波, 以及紧随爆轰波后的稀疏波区,这一结果对应于一维理论分析中的CJ解. 透射波系基本不受 初始压强影响;初始温度也只成比例地改变流场整体速度,温度越高,速度越快;对波系起 实质影响作用的是入射激波强度,激波越强,则整个透射流场呈现偏向激波的趋势;理论分 析还指出,稀疏波区的出现不可避免,当激波强度趋于声波稀疏波区趋于消失,激波越强则 疏波区趋于扩大. 两波对撞存在一个有限的转变阶段,透射爆轰首先减缓,接着迅速迸发为 过驱爆轰,然后再逐渐平衡为CJ爆轰. 对于强不稳定的燃气,对撞后爆轰波在空间上的发展 极不均衡,一些区域发生火焰面与诱导激波的严重脱离,随后的火焰面失稳发展为诱导激波 区内的爆轰波,实验观察到了这种爆轰在烟迹板上留下的极为精细的迹线.     

The criterion of the existence or inexistence of transverse shock wave at wedge supported oblique detonation wave

A simplified theoretic method and numerical simulations were carried out to investigate the characterization of propagation of transverse shock wave at wedge supported oblique detonation wave.After solution validation,a criterion which is associated with the ratio Φ (u 2 /u CJ) of existence or inexistence of the transverse shock wave at the region of the primary triple was deduced systematically by 38 cases.It is observed that for abrupt oblique shock wave (OSW)/oblique detonation wave (ODW) transition,a transverse shock wave is generated at the region of the primary triple when Φ < 1,however,such a transverse shock wave does not take place for the smooth OSW/ODW transition when Φ > 1.The parameter Φ can be expressed as the Mach number behind the ODW front for stable CJ detonation.When 0.9 < Φ < 1.0,the reflected shock wave can pass across the contact discontinuity and interact with transverse waves which are originating from the ODW front.When 0.8 < Φ < 0.9,the reflected shock wave can not pass across the contact discontinuity and only reflects at the contact discontinuity.The condition (0.8 < Φ < 0.9) agrees well with the ratio (D ave /D CJ) in the critical detonation.