RDX炸药爆轰产物物态方程

采用基于统计物理的爆轰产物物态方程模型,研究了不同初始密度下RDX炸药爆轰产物性质,发现炸药初始密度小于1.20g/cm3时,炸药爆轰后只有气相产物生成,随着炸药初始密度进一步增大,爆轰后有固相碳析出。炸药初始密度小于1.60g/cm3时,各气相产物按均匀混合处理,计算的爆速、爆压和实验值符合较好。随着炸药初始密度进一步增大,计算的爆速、爆压与实验值出现较大的偏差,经分析发现此处可能出现了超临界流体相分离现象,在物态方程模型中考虑了这种效应后,得到较好的结果。     

爆轰产物驱动飞片运动数值模拟

化爆加载装置是研究材料动态响应特性的重要冲击加载装置之一。如果炸药与飞片直接贴附，将难以获得较低的飞片速度，而且飞片中载荷的上升时间较短。为研究材料在冲击压缩低压状态和准等熵压缩条件下的动态响应特性，对利用爆轰产物经过空气隙驱动飞片运动的加载装置进行了二维数值模拟，获得了一些规律性认识。     

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

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

基于统计物理的爆轰产物物态方程研究

用van der Waals等效单组分流体模型和Ross硬球微扰理论软球修正模型，计算爆轰气相产物的状态方程；用石墨相、金刚石相、类石墨液相和类金刚石液相４种相态描述凝聚成分，由Gibbs自由能最小确定不同状态下的凝聚产物相态.对爆轰产物混合系统采用Gibbs自由能最小原理，通过化学平衡方程组求解炸药爆轰产物系统的平衡组分，计算结果与Becker-Kistiakowsky-Wilson (BKW)和Lennard-Jones-Devonshire结果相近.使用该理论对炸药的爆轰参数做了预言，与BKW，Jo 关键词： 爆轰产物 物态方程 化学平衡方程组     

新型爆轰产物物态方程

 从位力（Virial）理论和相似假设出发，建立了爆轰产物物态方程，命名为VLW物态方程。本物态方程既适用于高压下炸药爆轰性能的计算，也适用于较高压强与中等压强状态下火炮发射药与火箭推进剂燃烧性能的计算。     

混合物冲击温度的数值模拟

 介绍温度的数值模拟方法，采用三维动力学有限元程序模拟冲击压缩状态，分析了数值模拟的准确性。模拟了铜和铝两种金属单质组成混合物的冲击压缩性能，在分析了热传导影响后，计算了混合物的冲击温度和冲击绝热线。     

凝聚炸药中超压爆轰的实验研究

 采用飞片碰撞技术，在TNT/RDX（40/60）炸药中获得了2.5倍于正常爆轰的最大超压值，得到了超压爆轰下爆轰产物物态方程p=Aρ^{k+A₁(p-p_J)}（p－爆压，单位GPa，ρ－密度，单位kg/m³，A=ρ_J/ρ^k_J，p_J=27.06 GPa，ρ_J=2.3×10³ kg/m³，k=2.77，A₁=2.7×10^-3 GPa^-1，下表J代表正常爆轰状态）。该方程还可以较好地描述超压爆轰产物的二次冲击状态。     

H32.685N16.8430O3.0和H8.154N4.5770O3.0炸药爆轰参数的数值模拟

用Gibbs自由能最小原理,通过解化学平衡方程组,求解HNO型炸药H32·685N16·8430O3·0(70/30Hydrazine/HydrazineNitrate)和H8·154N4·5770O3·0(21/79Hydrazine/HydrazineNitrate)爆轰产物系统的平衡组分,计算结果与用BKW和LJD方法计算的结果相近。以Ree修正的WCA状态方程并考虑Ross软球修正的硬球微扰理论作为爆轰气相产物的物态方程,用自编的程序对H32·685N16·8430O3·0和H8·154N4·5770O3·0炸药爆轰参数作了预言,爆轰CJ点的爆速、爆压和爆温的计算结果与实验值吻合得很好。     

HMX炸药燃烧转爆轰数值模拟

以两相流模型为基础,气相产物状态方程采用基于统计物理的类CHEQ的计算结果,建立了HMX炸药的燃烧转爆轰数学模型。采用CE/SE方法模拟了颗粒度为125μm的HMX炸药的燃烧转爆轰过程,得到了爆轰参数及流场变化规律。模拟了装填密度对HMX炸药燃烧转爆轰的影响,并与实验进行了对比。数值模拟结果表明,在相同的点火条件下,爆轰成长距离在一定范围内随装填密度呈"U"形变化。     

爆轰产物物态方程（II）—爆轰的ZND理论不成立吗？（下）

用Ｊｏｎｅｓ公式检验爆轰的ＺＮＤ理论是另一种检验方法。这种检验方法不需要知道爆轰产物物态方程的形式,只需知道其存在即可,百密一疏,美国科学家用这种方法检验得出的结论是ＺＮＤ理论不成立。更稹密地运用这种方法得出的检验结论是爆轰的ＺＮＤ理论成立。     

水下爆炸冲击波数值仿真研究

由于在水下爆炸冲击波的数值仿真研究中,水的状态方程、人工黏性系数和网格尺寸对数值计算结果影响很大,采用常规TNT炸药的水下爆炸为例,以冲击波的峰值压力和比冲量为衡量指标,研究了这3个主要影响因素对数值仿真结果的影响。首先,通过采用常用的5种水的状态方程进行系列仿真,给出了各种状态方程的适用范围;其次,讨论了人工黏性系数对计算结果的影响,并给出了一次与二次人工黏性系数的建议取值范围;最后,通过对不同炸药当量及不同网格尺寸开展系列运算,从而得到不同炸药当量在满足工程计算精度要求下所对应的建议网格尺寸,并得到了不同炸药当量所对应的建议网格尺寸的表达式。     

爆炸作用驱动液体抛撒初始阶段数值仿真

 针对灭火战斗部引战配合设计,研究爆炸作用驱动液体（水基灭火剂）抛撒初始阶段的特性。依据典型灭火战斗部的结构建立数值仿真模型,采用LS-DYNA软件进行仿真计算,分别获得了加载系数（驱动装药与液体的质量比）对壳体破裂时间和壳体破裂初期液体/空气界面运动速度的影响规律,以及液体/空气界面后续的运动特性。进行了同条件下爆炸驱动水基灭火剂抛撒实验,得到了不同加载系数条件下壳体破裂初期液体/空气界面的运动速度。实验结果与数值仿真结果吻合良好,表明所采用的数值仿真方法具有正确性和有效性,对相关问题的研究具有重要参考价值。     

强爆炸早期火球现象的一维数值研究

采用一维球对称辐射流体力学方程组,构建直接求解辐射输运方程的一种数值方法,得到强爆炸早期火球发展过程中火球阵面、弹壳冲击波的形成发展过程.计算结果与Brode结果进行对比,二者一致,证明数值方法的合理性.数值模拟爆炸源区不同当量、不同密度以及不同爆炸高度等条件下强爆炸早期火球发展过程中火球阵面、弹壳冲击波的形成发展过程,对计算结果进行比较分析.     

面层混凝土对跑道内爆毁伤影响的仿真研究

 采用考虑混凝土拉伸破坏的TCK失效模型,运用AUTODYN软件对跑道内爆毁伤进行了数值模拟。研究了炸点置于面层混凝土底部时,不同装药量和不同面层混凝土厚度对爆破漏斗坑参数的影响规律。结果表明：数值模拟结果较好地反映了混凝土因拉伸损伤引起的破裂和抛掷漏斗现象。基于数值模拟得到的漏斗坑参数值,建立了漏斗坑半径、空腔半径与面层厚度的无量纲表达式,计算得到的漏斗坑参数值得到了实验验证。     

点起爆炸药驱动平板飞片运动的数值模拟研究

 为研究点起爆发散爆轰波驱动飞片的运动特性,自行研制了二维流体动力学程序TDY2D,采用由实验所标定炸药的状态方程参数,对3种长径比相差较大的点起爆炸药驱动平板飞片运动的模型进行了数值模拟计算,并对结果进行了分析。计算所得的飞片运动历史和击靶波形均与实验符合较好,具有较高的计算精度,从而验证了该程序的准确性和有效性。     

爆轰产物JWL状态方程应用研究

 通过对凝聚态炸药爆轰产物JWL状态方程的研究,提出了一种确定JWL状态方程参数的γ拟合法,即利用γ律状态方程来拟合JWL状态方程参数。这种新方法不需进行圆筒试验,与现有方法比较,具有经济、安全、方便、准确的特点。利用这种方法拟合确定了4种常用炸药TNT、C-4、PETN、HMX的JWL状态方程参数,通过与圆筒试验法得出的JWL状态方程p-V曲线进行对比,证明γ拟合法具有较高的精度,完全可以满足计算爆炸力学的应用。     

不同方式地下爆炸地震耦合效应的数值模拟

通过一种Euler与Lagrange相结合的有限差分法,开展地下爆炸地震耦合效应的数值模拟.比较花岗岩介质中填实和空腔两种爆炸方式在地震耦合效应方面的差异,进一步分析初始空腔大小在填实以及空腔爆炸地震效应中的作用以及空腔爆炸解耦因子随频率的变化特性.计算结果表明,花岗岩介质中完全空腔解耦因子能达到50~60倍.     

应用Monte Carlo方法模拟爆轰产物状态方程的研究进展

 应用Monte Carlo方法模拟爆轰产物状态方程,区别于建立在球形分布分子势基础上的传统爆轰产物状态方程计算方法,从原子水平上直接模拟实际炸药的爆轰产物状态方程,克服了传统方法中混合产物状态方程不精确的缺点。综述了该方法的历史背景及发展现状。对模拟中宏观量的描述、统计平均的计算、势函数的选取及边界条件进行了详尽描述。     

水下爆炸气泡脉动周期的简便计算方法

将爆轰产物状态方程与气泡运动方程相结合,提出了一种水下爆炸气泡脉动周期的数值计算方法,能够快速计算不同炸药水下爆炸的气泡脉动周期。结果表明:采用JWL状态方程与采用γ律状态方程的气泡周期计算方法相比,前者计算TNT炸药水下爆炸气泡脉动周期的误差更小,计算误差小于3%;通过对比RS211炸药水下爆炸实验结果,进一步验证了该数值计算方法对含铝炸药水下爆炸气泡脉动周期的计算同样具有良好的适用性。     

Numerical Simulation of the Response of an Aluminum Target Impacted by an Intense Laser Driven Flyer

1　Ｉｎｔｒｏｄｕｃｔｉｏｎ　　Ｔｈｅｉｍｐａｃｔｔｅｃｈｎｉｑｕｅｏｆｆｌｙｅｒｄｒｉｖｅｎｂｙｃｈｅｍｉｃａｌｅｘｐｌｏｓｉｖｅａｎｄｇａｓｇｕｎｈａｓｂｅｅｎｕｓｅｄｔｏｓｔｕｄｙｈｉｇｈ ｐｒｅｓｓｕｒｅｅｑｕａｔｉｏｎ ｏｆ ｓｔａｔｅ (ＥＯＳ)ｆｏｒｍａｎｙｙｅａｒｓ[1] .Ｂｅｃａｕｓｅｔｈｅｆｌｙｅｒｉｓａｔａｌｏｗｅｎｅｒｇｙｄｅｎｓｉｔｙａｎｄｅｎｔｒｏｐｙｅｎｈａｎｃｅｍｅｎｔｓｔａｔｅｉｎｔｈｉｓｃａｓｅ ,ｉｔｃａｎｅａｓｉｌｙｒｅｃｏｖｅｒｉｔｓｉｎｉｔｉａｌｓｔａｔｅｂｅｆｏ…