薄片炸药与固体靶冲量耦合的计算模型

薄片炸药可用于模拟强脉冲X射线辐照下空间结构的响应,其与固体靶的冲量耦合作用是进行模拟载荷设计的基础。本文中基于有限元模拟技术,研究不同直径薄片炸药的边界稀疏作用和安装支座的反射作用对薄片炸药与固体靶冲量耦合的影响,结果表明:边界稀疏作用对薄片炸药的比冲量起到减弱的作用,由边界稀疏作用引起的薄片炸药比冲量的变化量和薄片炸药直径的倒数近似成线性关系;薄片炸药安装支座的反射作用对薄片炸药的比冲量起到增强的作用,由支座反射作用引起的薄片炸药比冲量的变化量和薄片炸药直径的平方近似成反比;随着薄片炸药直径的增加,薄片炸药比冲量趋于定值,边界稀疏作用和支座反射作用可以忽略;存在一个临界直径,当薄片炸药的直径为临界直径时,薄片炸药的边界稀疏作用和支座反射作用达到平衡,其比冲量大小和薄片炸药的直径为无穷大时的比冲量值相等。     

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

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

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

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

一类爆轰合成用乳化炸药的爆轰反应特征

结合热分解特性、爆轰产物状态及数值计算探讨所制备的含Fe、Mn(Zn)元素乳化炸药的爆轰反 应特征。对该系列乳化炸药及其爆轰产物分别进行DSC和DTG实验、XRD检测和TEM 表征,通过比较不 同类型炸药及其爆轰产物的热分解特征、爆轰产物的成分和形貌、晶型畸变度等研究该类炸药爆轰反应特征; 并通过数值计算研究和佐证了不同温度和压强状态下,爆轰固体产物可能存在状态。结果表明:硝酸铵有利 于炸药爆轰时形成相对比较均匀和平稳的爆轰反应结构,而这种结构有利于可重复性合成比较单一和均匀 的爆轰产物。高密度炸药的爆轰产物比低密度炸药的相对纯净。低密度炸药不完全爆轰,导致了爆轰产物杂 质较多;数值计算结果表明,爆轰固体产物分布相图可以辅助分析爆轰固体产物可能存在状态及Mn(Zn)的 爆轰合成机制。     

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

采用电探针及压力传感器测试技术对密度为1.597 g/cm3的固体B炸药（TNT/RDX=40/60）的燃烧转爆轰性能进行研究。实验结果表明,在较强的约束条件下（45号钢管,内径20 mm,外径64 mm,长500 mm）,B炸药形成了DDT现象,诱导爆轰距离为295~310 mm。     

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

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

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。     

钝感炸药的超压爆轰与冲击起爆过程数值模拟

采用Hybrid反应率结合修正的JWL方程,研究了LX-17、超细TATB等钝感炸药的冲击起爆（SDT）过程,并计算了爆轰波的对碰现象。结果表明,该方法计算钝感炸药的冲击起爆过程与实验数据符合较好;计算爆轰波对碰区的峰值压力提高了10％。     

JOB-9003和JB-9014炸药平面爆轰驱动飞片的对比研究

采用双灵敏度激光速度干涉仪(VISAR)测量了JOB-9003和JB-9014两种炸药平面一维爆轰驱动飞片的速度历史,比较了二者的做功能力;并用DYNA2D程序进行了数值模拟。在有JO-9159炸药的条件下,实验和计算结果表明,JB-9014炸药的做功能力可以达到JOB-9003炸药的75.0%;无JO-9159炸药条件下的计算结果表明,JB-9014炸药的做功能力约为JOB-9003炸药的71.2%     

松弛方法在计算凝聚炸药爆轰问题中的应用

利用松弛近似,将非线性的凝聚炸药爆轰控制方程转化为线性的松弛方程组,并采用五阶WENO格式和五阶线性多步显隐格式对线性松弛方程组进行空间方向和时间方向的离散,由此建立具有高精度和高分辨率性质的计算凝聚炸药爆轰的松弛方法。建立的松弛方法可以避免求解Riemann问题及计算非线性通量的Jacobi矩阵,同时无需分裂处理反应源项。通过对凝聚炸药的平面一维定常爆轰波结构及球面一维聚心、散心爆轰起爆和传播过程的数值模拟,验证了所建立的松弛方法能够很好地计算凝聚炸药爆轰问题。     

混合炸药爆速预报的新方法

为了提高混合炸药爆速尤其是含铝混合炸药爆速的计算精度,基于体积加权法思路构造一种新的混合炸药爆速计算公式,对近50种常见的混合炸药配方进行了计算,并与实验值进行对比,结果很好吻合,最大误差低于3%,平均误差低于1%,与BKW和Urizar常用方法相比,计算精度有显著提高。

     

爆炸复合中起爆区不定常段的确定

对2岩石炸药的定常爆轰形成过程进行了实验研究;对起爆端稀疏波效应进行了计算。实验结果表明,距起爆端约10倍装药厚度以内为不定常爆轰区。距离起爆端小于5倍装药厚度时,飞板速度及碰撞角增幅较大,5～10倍装药厚度,两者增幅趋缓,10倍装药厚度以后,爆速、飞板速度及碰撞角逐步进入定常状态。飞板碰撞角的工程计算与实验结果的一致性较好。     

爆轰波对碰区产物驱动金属圆管的研究

利用闪光X光照相、光学高速分幅照相实验,测得了金属圆管在两个爆轰波对碰后产物驱动外壳的膨胀、变形过程,该过程表明爆轰波对碰作用使得对碰处金属圆管相对贯穿断裂时间明显提前,对碰处金属圆管径向膨胀速度较非对碰处明显增加。利用拉格朗日二维流体动力学程序TTD2C,数值模拟了对碰驱动过程,计算结果与实验结果符合很好,且用Taylor断裂判据得到的金属圆管相对贯穿断裂时间也明显提前。     

一种以TATB/HMX为基炸药的到爆轰距离

利用光纤探针/光电转换器/示波器技术和楔形炸药块方法,采取电炮驱动薄飞片冲击加载产生高压短脉冲激励的装置,测量了一种以TATB/HMX为基炸药的到爆轰距离,研究了到爆轰距离与初始冲击波压力幅值和脉宽的关系,给出了该炸药的POP曲线和表达式。实验结果有助于了解压力脉宽对钝感炸药的冲击起爆和爆轰成长的影响。对于短脉冲冲击加载,入射压力脉宽对炸药的到爆轰距离影响明显,相同实验条件下,压力脉宽越长,炸药的到爆轰距离越短;相同压力脉宽下,加载压力越高,炸药的到爆轰距离越短。     

PBX-01炸药水中爆轰产物状态方程研究

提出通过水中实验确定炸药的水中爆轰产物JWL状态方程参数的方法;选择PBX-01高能炸药进行水中实验,利用ANSYS/LS-DYNA程序建立炸药的水中实验模型,将实验结果与数值计算结果进行对比,确定PBX-01炸药水中爆轰产物的JWL状态方程参数。研究结果显示,圆筒实验确定的JWL参数在反映炸药水中爆轰产物的膨胀状态时有所不足,水中实验确定的JWL状态方程参数能够更准确地描述PBX-01炸药水中爆轰产物的膨胀过程,因此对水中爆炸的研究需要通过水中爆炸实验建一套状态方程参数。

     

高-低爆速圆板炸药串联爆轰引起平面爆轰波的变凸现象

为分析高-低爆速圆板炸药于空中串联爆轰时平面爆轰波在传播过程中演变成凸面波的现象,分别对平面爆轰波阵面后爆轰产物状态参数的分布、有效药量与柱面装药几何尺寸的关系、爆轰产物的状态方程及强爆轰关系式进行了讨论。并以100 mm50 mm的TNT与RC炸药串联爆轰为例,描述了平面爆轰波演变为凸面爆轰波的过程,预估了爆轰波的平面范围和波形差。预估结果与实验结果基本符合。     

一点起爆半球形装药装置中散心爆轰波传播的数值模拟研究

运用自行研制的二维流体动力学程序TDY2D,对一点起爆半球壳装药装置中散心爆轰波的传播过程进行数值模拟计算.对冲击波到达炸药和飞层各界面的波形的计算值与实验值进行了对比,两者符合较好,从而验证了该程序的准确性和有效性;对爆轰波的传播过程进行了分析研究,获得了一些散心爆轰波传播的规律性认识,为进一步研究散心爆轰波特性提供了数值分析基础.     

起爆高密度TATB炸药的飞片速度阈值

采用实验与数值模拟相合的方法研究了飞片起爆钝感TATB炸药的性能。设计了飞片可靠起爆钝感TATB炸药的起爆序列,利用全光纤位移干涉测速系统分别测出飞片可靠起爆和未起爆TATB炸药的速度,初步确定飞片起爆钝感TATB炸药的起爆速度阈值。采用DYNA2D程序对飞片起爆TATB炸药过程进行数值模拟,模拟结果与实验结果相符合。

     

爆速对爆炸焊接铝/不锈钢复合管界面及结合性能的影响

通过在粉状乳化炸药中添加不同比例的密度调节剂,配制了爆速范围为1 450~2 550 m/s的低爆速炸药;采用该爆速炸药进行了铝/不锈钢复合管爆炸焊接实验,结合最小碰撞速度理论,对实验结果及其界面微观结构和结合强度进行了测试和分析,确定该复合管爆炸焊接的合适爆速约为1 950~2 150 m/s,其结合质量能够满足后续加工要求;同时发现界面由介于直线与波形之间的波状形态组成,且呈现不太规则的扁平波状结合,经分析,炸药爆速、复合管的爆炸焊接环境和爆炸产物飞散条件对界面结合波形及熔化层厚度有很大影响。