多破片对柱壳装药冲击起爆速度阈值的数值模拟研究

为了研究实战环境中多个钨球破片对导弹战斗部(柱壳装药)的冲击起爆问题,采用AUTODYN-3D数值模拟软件,基于单破片撞击柱壳装药模型,建立多破片撞击柱壳装药的模型,开展了不同钨球个数、空间碰撞位置间隔(撞击角θ、轴向球心距l)及时间间隔对冲击起爆特性影响的数值模拟,获得了带壳B炸药的起爆速度阈值。结果表明:相同条件下,随着钨球个数的增加、空间碰撞位置间隔的减小,起爆速度阈值逐渐减小,6个钨球同时撞击的起爆速度阈值约为单个作用下的50%;双钨球作用下,柱壳装药相较于平板装药更难以起爆;双钨球间隔撞击柱壳装药时,起爆速度阈值均随着撞击时间间隔的增大而先减小后增大,最小起爆速度阈值约为同时撞击时的95%,且|θ₂|<|θ₁|(θ₁为第1个钨球的撞击角,θ₂为第2个钨球的撞击角)时更容易起爆柱壳装药。     

圆柱形有壳装药侧向飞散速度分布的预估

对圆柱形有壳装药侧向飞散速度的分布,本文提出了一种改进的预估方法,与实验结果有较好相符。这一方法是从装药侧面冲量分布出发,加以适当变换(经验修正),既比用HEMP程序计算来得省钱,也可象用HEMP程序计算那样适合多种起爆方式。为实用目的,文中还给出了侧向飞散最大速度值的估算方法。     

垫层对破片冲击起爆带壳炸药影响的数值模拟

将战斗部简化为带钢壳和垫层的炸药组件,用3种不同厚度的材料作垫层,利用非线性有限元软件AUTODYN对高速破片侵彻、引爆带壳炸药的作用过程进行了数值模拟,得到了高速破片冲击起爆带壳带垫层炸药的波后压力剖面。对不同厚度多种材料作为垫层时带壳炸药的冲击到爆轰距离进行了对比分析,结果表明,通过控制垫层的材料和厚度,可以对破片冲击起爆带壳炸药进行有效抑制。研究结果对战斗部破片起爆的防护设计具有指导意义。

     

不同起爆方式下炸药驱动柱壳膨胀断裂的数值模拟

采用数值模拟方法,研究不同起爆方式下内部炸药对金属柱壳的加载历程差别。计算结果表明,不同起爆方式下柱壳内表所经受的压力载荷、做功历史显著不同,分析了引起差异的原因。提出一种半解耦数值方法,使柱壳网格划分得更细,可以模拟不同起爆方式下柱壳的膨胀破坏过程,通过理论分析说明了这种方法的合理性。

     

杆式弹对厚壁壳体装药冲击起爆机制模拟分析

为研究高速杆式弹冲击厚壁壳体装药的起爆机制,运用冲击物理显式欧拉型动力学SPEED软件,开展了不同弹径和弹长的钨合金杆式弹与厚壁壳体Comp-B装药相互作用过程的数值模拟,采用升降法获得弹体起爆装药临界着速及装药起爆位置变化。研究结果表明:弹体起爆装药临界着速随弹径增大而显著降低,随弹长增大呈先降低后平缓变化的规律;弹体以临界着速起爆装药时,存在2种装药起爆机制,即弹体贯穿壳体后的宏观剪切起爆和未贯穿壳体的低速冲击起爆,且其机制随弹体着速在临界着速以上继续提高会发生转变,最终均会转变为高速冲击起爆机制;装药起爆位置均发生在炸药壳体交界面后一定距离处,相同机制下此距离随弹体着速提高而减小。     

圆柱壳体装药偏心多点起爆下破片速度的分布

针对偏心起爆战斗部破片速度增益的问题,提出爆轰波碰撞形成马赫超压是引起破片速度增加的原因。利用AUTODYN软件,模拟偏心起爆战斗部从壳体径向膨胀、表面产生裂纹到最后形成破片的整个过程,并将模拟得到的破片速度与实验数据对比,两者吻合较好;简化Whitham方法并结合Gurney速度公式得到偏心起爆战斗部定向破片速度和定向区域的计算方法,同时在保证破片初速的前提下,研究偏心多点起爆下起爆点数的选择标准。研究结果表明：偏心多点起爆下定向破片初速增益约34%,定向区域范围约30°,起爆点数的选择与壳体长度和装药口径相关。

     

粘弹塑性热点燃烧模型的冲击起爆理论

本文作者建议一种粘弹塑性燃烧模型作为冲击起爆的物理图象,并由此导出了著名的冲击起爆判据p~2=常数。     

线起爆膨胀柱壳实验加载及诊断技术

基于电爆炸丝引爆炸药、继而驱动尼龙对金属柱壳进行碰撞加载的方式,在金属柱壳中部、半柱高范围内实现了一维柱面膨胀加载。同时,基于沿轴向的加载（或径向速度）一致性和沿环向的加载（或径向速度）轴对称性,提出了一维柱面加载的有效性判据。相比于滑移爆轰加载,一维柱面加载方式具有应力状态相对简单、易简化为二维轴对称问题分析的优点。在柱壳断裂诊断方面,建立了分布式表面速度诊断方法来监测柱壳圆周范围内的初始断裂。其原理为：均匀承载壳体断裂引起的局部承载失效将导致均匀速度曲线簇出现分叉（或演化趋势变化）。与高速分幅照相诊断方法相比,分布式表面速度诊断方法可准确获取柱壳圆周范围内的初始断裂信息（含断裂时刻和断裂位置）。利用建立的线起爆膨胀柱壳实验加载和诊断技术,获得了304钢和45钢柱壳的一维柱面动态拉伸初始断裂性能数据（含断裂应变、平均应变率）,其中,45钢柱壳的断裂应变（或延展性）低于304钢柱壳的。

     

一次起爆FAE中心装药的实验研究

为了探讨不同组分的中心装药对不同FAE燃料可能的一次起爆情况,采用物理方法对中心装药的配方进行设计,并对不同的FAE燃料进行了一次起爆实验研究,结果表明含有诱发剂的中心装药可以实现FAE燃料的一次起爆。为了进一步弄清中心装药中诱发剂含量对FAE燃料一次起爆过程的影响,利用不同诱发剂含量的中心装药对自行研制的新型复合燃料进行了一次起爆实验研究,CCD像机现场拍摄记录及压力测试结果表明,中心装药中的诱发剂含量会直接影响FAE燃料抛撒、延迟点火时间、爆炸波超压大小以及爆炸作用效果,诱发剂含量存在一个最优值,中心装药中加入适量诱发剂会增强FAE的一次起爆效果。     

偏心起爆对战斗部装药能量分配增益的影响

为研究不同方式的偏心起爆对炸药装药能量分配及增益的影响,建立了偏心起爆战斗部的计算模型,通过局部装填比这一变量,给出了偏心起爆战斗部破片的初速计算公式。采用数值模拟与试验验证结合的方法,对六分位条件下不同偏心起爆方式的破片速度增益和能量增益进行了对比,得出以中心起爆为基准,分别以邻位双线、连位三线、间位双线、偏心单线方式起爆,定向方位内破片的速度增益依次增大;邻位双线起爆时,目标方向破片速度增益达25.47%,定向区域破片动能占总能量的24.57%,能量增益超过40%。     

球形装药动态爆炸冲击波超压场计算模型

为获得球形装药动态爆炸冲击波超压场计算模型,对静态爆炸冲击波超压Baker计算公式加入修正因子进行修正,并建立了构造包含装药运动速度、对比距离和方位角的修正因子函数的方法。为获得修正因子的函数表达式,采用高精度显式欧拉流体动力学软件SPEED针对具有典型运动速度的球形装药空中爆炸过程进行了数值模拟,得到了沿装药不同对比距离和方位角处的动态爆炸冲击波超压峰值。在对数值模拟结果处理的基础上,经过数据拟合获得了动态爆炸冲击波超压场计算模型。校验结果表明,该模型能较准确描述动态爆炸冲击波超压分布,具有普适性。     

金属圆柱壳受大质量低速冲击的屈曲变形

对钢质和铜质金属圆柱壳的轴向冲击动力响应进行了实验研究,记录了两种不同材料圆柱壳在大质量低速冲击下的冲击力时程曲线,得到其屈曲模态。采用高速摄像及模拟技术给出了钢质圆柱壳渐进屈曲的全过程,为理解钢质圆柱壳的屈曲机理提供了直观的结果。黄铜质圆柱壳在大质量低速冲击下, 出现整个壳面滿布屈曲波纹的塑性动力屈曲现象,说明高速冲击不是产生塑性动力屈曲的充要条件。像铜这样具有高密度的韧性材料,在大质量低速冲击下,会在轴向产生持续的压缩塑性流作用而出现塑性动力屈曲现象。     

圆柱形钨破片撞击铝靶时的跳飞临界角

采用理论计算、数值模拟与实验相结合的方法,研究了直径5.7 mm、长6.7 mm的圆柱形破片以800~1 200 m/s的速度撞击2~10 mm厚铝靶时的跳飞特性。建立了破片斜侵彻有限厚靶板的跳飞临界角理论模型,计算得到破片跳飞临界角与破片入射速度、靶板厚度的关系,并与模拟值、实验值对比,三者吻合较好。结果表明：破片撞靶速度相同时,随着靶板厚度的增加,破片的跳飞临界角减小。靶板厚度相同的情况下,在所计算的速度范围内,入射速度越大,破片跳飞临界角越大。速度在800~1 200 m/s时,破片撞击2 mm厚靶板的跳飞临界角为81°~81.25°;撞击4 mm厚靶板的跳飞临界角为72.5°~76.25°。

     

Motion of the detonation products of a point-ignited charge in a cylindrical shell

A numerical solution is obtained for the two-dimensional nonsteady problem of the motion of detonation products from a cylindrical high-explosive charge enclosed in a shell with the initiation of detonation at a central point in the end of the charge. The detonation products propagate in vacuum. The strength of the shell is not considered. A three-term equation of state is used for the detonation products.     

Critical velocity of sandwich cylindrical shell under moving internal pressure

Critical velocity of an infinite long sandwich shell under moving internal pressure is studied using the sandwich shell theory and elastodynamics theory. Propagation of axisymmetric free harmonic waves in the sandwich shell is studied using the sandwich shell theory by considering compressibility and transverse shear deformation of the core, and transverse shear deformation of face sheets. Based on the elastodynamics theory, displacement components expanded by Legendre polynomials, and position-dependent elastic constants and densities are introduced into the equations of motion. Critical velocity is the minimum phase velocity on the desperation relation curve obtained by using the two methods. Numerical examples and the finite element （FE） simulations are presented. The results show that the two critical velocities agree well with each other, and two desperation relation curves agree well with each other when the wave number k is relatively small. However, two limit phase velocities approach to the shear wave velocities of the face sheet and the core respectively when k limits to infinite. The two methods are efficient in the investigation of wave propagation in a sandwich cylindrical shell when k is relatively small. The critical velocity predicted in the FE simulations agrees with theoretical prediction.     

圆柱装药在偏心定向起爆时水中近场压力特性

对UHL-5装药在侧向9点偏心定向起爆和几何中心点起爆时的水下爆炸近场峰值压力分布特性进行实验研究和数值模拟,得到2种起爆条件下不同距离与方位角测点处冲击波峰值压力和侧向9点起爆时的定向增益区域.结果表明,对于5 kg圆柱形UHL-5装药,采用侧向9点偏心定向起爆方式时,在爆距750 mm、方位角90°范围内的流场区域...     

内爆和柱壳条件下无氧铜的层裂特性

利用任意反射面位移干涉系统(DISAR)激光测速技术,成功地获得了滑移内爆加载和柱壳结构条件下无氧铜的内表面（自由面）速度剖面,并对其层裂特性进行了初步分析。结果表明:(1)在固定炸药和改变无氧铜圆管壁厚条件下,层裂片厚度随着圆管壁厚h的减小而增加;以圆管壁厚h为参照进行归一化,则相对层裂片厚度(/h)随相对装药厚度(he/h)的增大而增加,这种规律与以往对20钢的研究结果一致,但圆管发生层裂的临界条件,却显示出明显的材料相关性。(2)初步来看,无氧铜的层裂强度对结构的依赖性不明显,而与加载脉冲的幅值和宽度相关。(3)受无氧铜粘性和Taylor波衰减的影响,无氧铜的层裂强度随管壁厚度的增加而略有降低;同时,材料分散性也对此有一定影响。     

圆柱形爆轰波的二维数值模拟

基于带化学反应的二维Euler方程,对圆柱形爆轰波的直接起爆和传播过程进行了二维数值模拟研究,拟分析起爆条件和初始压强对圆柱形爆轰波形成和传播的影响。研究发现,圆柱形爆轰波起爆成功向外传播的过程中,新的三波结构的生成标志着爆轰波进入稳定传播阶段。在起爆能量足够的情况下,起爆半径(曲率)的大小决定着三波结构初始形成时的数目和传播半径,起爆压强对其基本不产生影响;起爆半径大(曲率小)时,三波结构初始形成时的传播半径大、数目多,圆柱形爆轰波进入稳定传播阶段的传播距离长;数值模拟中,初始压强的提高,有助于圆柱形爆轰在较短的传播距离内进入稳定传播阶段。