预制破片与轻质壳体阻抗匹配对破片初速及完整性的影响

周涛 沈飞 王辉

周涛, 沈飞, 王辉. 预制破片与轻质壳体阻抗匹配对破片初速及完整性的影响[J]. 高压物理学报, 2018, 32(4): 045104. doi: 10.11858/gywlxb.20170602
引用本文: 周涛, 沈飞, 王辉. 预制破片与轻质壳体阻抗匹配对破片初速及完整性的影响[J]. 高压物理学报, 2018, 32(4): 045104. doi: 10.11858/gywlxb.20170602
ZHOU Tao, SHEN Fei, WANG Hui. Influence of Impedance Matching between Prefabricated Fragments and Light Shell on Initial Velocity and Completeness of Fragments[J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 045104. doi: 10.11858/gywlxb.20170602
Citation: ZHOU Tao, SHEN Fei, WANG Hui. Influence of Impedance Matching between Prefabricated Fragments and Light Shell on Initial Velocity and Completeness of Fragments[J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 045104. doi: 10.11858/gywlxb.20170602

预制破片与轻质壳体阻抗匹配对破片初速及完整性的影响

doi: 10.11858/gywlxb.20170602
基金项目: 

国防基础预研项目 00402020202

详细信息
    作者简介:

    周涛(1979-), 男, 博士, 副研究员, 主要从事爆炸力学及战斗部设计研究.E-mail:zt2756@sohu.com

  • 中图分类号: O389;TJ760.2

Influence of Impedance Matching between Prefabricated Fragments and Light Shell on Initial Velocity and Completeness of Fragments

  • 摘要: 依据杀伤战斗部装药对破片爆轰加载过程的特征,设计了与其较为相似的滑移爆轰单元结构实验模型,采用闪光X射线照相方法获得了预制破片和轻质壳体在两种典型排布顺序下的破片初速及破损情况,并结合应力波传播理论对实验结果进行了分析。结果表明:破片外置时,初始应力波由低阻抗金属材料向高阻抗金属材料传播,破片受到壳体传入的冲击波及空气传入的拉伸波作用,初速较高,轻微破损;破片内置时,初始应力波由高阻抗金属材料向低阻抗金属材料传播,虽然破片受到爆轰产物传入的冲击波及壳体反射的拉伸波作用,但初速相对偏低,易发生破损,甚至有明显层裂现象。

     

  • 全预制破片类杀伤战斗部的轻质壳体常采用高强度铝合金加工而成,主要用于提高战斗部结构强度、控制战斗部外形等,同时也便于装药过程的实施,而预制破片可布置于壳体内侧或外侧[1]。由于壳体材料与破片材料(如低碳钢、钨合金等)的阻抗存在较大差异,在爆炸加载过程中,不同的阻抗顺序将造成应力波的传播过程出现较大差别[2-3],可能导致破片的初速及完整性发生改变,进而影响战斗部的威力。因此,合理设计预制破片与轻质壳体的阻抗匹配顺序对于该类杀伤战斗部的优化设计具有重要意义。

    采用数值模拟方法研究该问题时,轻质壳体在加载初期由于变形过大而被失效删除,常难以对比不同阻抗顺序造成的差异[4]。因此,最直接的方式是对两种工况下的破片速度进行实际测量,并回收破片进行完整性分析。由于常用的断通靶等电测方式难以直接获得破片初速,一些科研人员基于圆筒结构实验模型,采用更高精度的转镜式高速扫描法及激光干涉测速法获取了爆炸加载过程中不同金属破片的初速,并分析了金属材料与炸药的匹配性[5-6];但该方法不适用于干扰因素较多的预制破片工况,且无法观测破片加速过程中的形态。而脉冲X射线测量方法虽然可观测破片加速过程中的形态,但常由于破片飞散方向杂乱难以准确获得破片的初速。正是由于上述观测及实验方法的限制,导致目前很少见到预制破片与轻质壳体的阻抗匹配方面的报道。

    本研究拟基于预制破片与轻质壳体的典型阻抗排布顺序,依据炸药装药爆轰加载释能特点,设计与杀伤战斗部爆轰加载过程较为相似的滑移爆轰单元结构实验模型(即采用滑移爆轰方式对“单条”破片进行加载),从而获得不同阻抗排布顺序下的破片初速及形态,并结合冲击波在不同介质间的传播规律,分析破片及壳体中的冲击波压力及粒子速度变化过程,以期能为相关杀伤战斗部的设计提供参考。

    图 1显示了杀伤战斗部全预制破片和轻质壳体的两种典型放置工况,即破片置于轻质壳体外侧和内侧。在爆轰加载过程中,破片外置时,爆轰产物首先冲击轻质壳体,使壳体中形成较强的冲击波,到达预制破片界面时,其冲击波主要由低阻抗介质向高阻抗介质传播;而破片内置时,冲击波则由高阻抗介质向低阻抗介质传播。因此,虽然这两种工况的装填比基本相同,但冲击波传播过程存在明显差异,进而可能影响破片的初速及变形。

    图  1  预制破片与轻质壳体的典型放置工况
    Figure  1.  Typical arrangements of prefabricated fragments and light shell

    然而,炸药爆轰加载过程并不是持续的冲击波驱动作用,由于爆轰产物膨胀过程分为高、中、低压3个阶段,且对战斗部破片及壳体的有效加载阶段主要位于高压及中压阶段,因此,一些研究认为炸药装药对破片及壳体的爆轰加载主要依靠两种作用力的共同作用,即冲击波的驱动力和气体爆轰产物的膨胀力[7]。冲击波驱动力可以使破片及壳体在加载初期速度迅速提高,但持续时间较短;而气体爆轰产物膨胀推动力虽然峰值相对偏低,但持续时间较长。以B炸药滑移爆轰驱动∅50 mm圆筒(即标准圆筒试验)为例。图 2显示了圆筒壁的加速历程。为了将两种作用力的作用效果区分开,采用us表示冲击波驱动力单独加载下圆筒壁的速度,ug表示气体产物膨胀力单独加载下圆筒壁的速度,ut(ut=us+ug)为两种作用力共同加载下圆筒壁的速度。从图 2中可以看出,整个加载过程可以明显地划分为两个阶段,即冲击碰撞加速阶段和气体产物膨胀力驱动加速阶段[7-8]。由于气体产物的膨胀力相对偏低,可假定在破片外置和内置两种工况下该阶段的加速效果差异较小,其差异主要在冲击碰撞加速阶段,因此可通过相应的单元实验对该阶段的加载规律进行分析。

    图  2  圆筒试验中圆筒壁膨胀速度-时间历程
    Figure  2.  Time histories of expansionvelocity of cylinder wall

    由于常用的断通靶测量方法在爆炸近场所受干扰因素较多,因此拟采用脉冲X射线系统对破片的分布及加速过程进行观测。从光学成像原理方面分析,脉冲X射线测量过程实际是将三维的实验模型投影到二维平面,因此,为了避免底片中的图像过于杂乱,只能对局部单位模块进行研究。为了尽可能地模拟实际战斗部的爆轰加载方式,采用滑移爆轰方式对“单条”破片进行加载,其中,起爆点为圆柱形装药一端面的中心处(如图 3所示),在柱形装药的中段,其径向冲量接近常量,则沿传爆方向的破片具有相同的加速度,且能够规则排列。此外,考虑到稀疏波的影响,破片距装药首、尾位置均空余一段距离,以保证破片所受爆炸冲量一致。这种加载方式对爆轰产物没有有效的约束,对破片的加载主要依靠冲击波驱动力作用。

    图  3  滑移爆轰驱动破片示意
    Figure  3.  Sketch of sliding detonation-driven fragments

    在滑移爆轰加载过程中,相邻破片运动的初始时间间隔固定,破片的飞散轨迹如图 4所示。根据X射线照片的判读结果可计算出各个破片的速度,其计算公式为

    图  4  破片飞散轨迹示意
    Figure  4.  Flying trajectory of fragments
    vi=ζ(SiSi1)ζL/D=SiSi1L/D
    (1)

    式中:vi为第i个破片的速度,ζ为实物尺寸与图像尺寸的比例,L为底片中相邻破片质心间的距离,Si为底片中第i个破片距装药下表面的距离,D为炸药的爆速。由于(1)式中不显含ζ,所以在获取破片的飞行速度时,不需要考虑底片中ζ的具体值。但有时在实验底片中无法显现装药的下表面,故难以判读破片运动的真实距离,只能获得其相对距离,同时考虑到不同破片的速度有时会出现波动,因此,实际分析时可根据(2)式计算破片的初速,即

    v=ni=2(SiSi1)(n1)L/D=DSnS1(n1)L=Dk
    (2)

    式中:n表示破片的总数,Sn-S1和(n-1)L分别表示底片中第n个破片与第1个破片之间的纵向和横向的距离。因此,(2)式也可以表示为:破片速度等于底片中破片曲线的斜率k与炸药爆速的乘积。

    实验时,破片及轻质壳体的材料、尺寸均根据实际战斗部进行选择。破片采用93钨合金方形破片,尺寸为6.5 mm×6.5 mm×3.8 mm,质量约为2.81 g;轻质壳体材料选用LY12铝合金,厚度为1.8 mm,宽度略大于破片。将破片排布并粘接在条状铝合金板上,然后再与炸药侧面粘接,以模拟破片内置及外置两种工况。实验所用主装药均为B炸药,密度为1.65 g/cm3,爆速约为7.8 km/s。由于装药的长度大于底片盒的宽度,故实验时将装药竖立放置,破片向侧向运动(实验布局如图 5所示);同时,由于采用的数字成像板不能二次曝光,图像中难以同时出现装药的基准线,因此起爆前先拍一发静止相片以判断装药的轴线是否与竖直悬挂的底片平行,然后再将成像板重新装入底片盒进行实验。

    图  5  实验布局
    Figure  5.  Experiment layout

    图 6为爆轰加载实验所获得的X射线摄影底片。图 6(a)显示了破片外置时钨破片和铝壳加速过程中的状态,其中并没有独立显示出铝壳体的迹线,且钨破片的下边界不太整齐,表明铝壳体与钨破片粘接在一起,两者速度基本一致;钨破片上方有部分细小的碎渣,其图像灰度与钨破片近似,表明爆轰加载过程中钨破片发生了轻微的破损。图 6(b)显示了破片内置时钨破片与铝壳加速过程中的状态,可以看出,钨破片下表面较为完好,但上表面处既存在一些细小的碎渣,同时还有较完整的条状碎裂体,其图像灰度与钨破片近似,可能是由于钨破片直接与炸药接触,其压力较高致使破片发生了较明显的层裂;钨破片图像上方出现了一条灰度相对较浅的线条,这是铝壳体的迹线(X射线摄影中,图像灰度与物质的密度相关),表明加载过程中铝壳体的速度明显高于钨破片。

    图  6  实验所获底片
    Figure  6.  Experimental film

    底片图像的横向距离(设坐标为x)可用于计算时间,纵向距离(设坐标为y)可用于计算破片的位移。由于每张底片的放大比不同,所以这些数据仅表示在底片中所处的位置,其数值并没有实际的物理意义,仅用于计算曲线的斜率。采用最小二乘法获得曲线的斜率k后,再按照(2)式计算破片的初速v

    图 7为底片的判读结果。对数据进行拟合计算时,考虑到边侧少量破片在加载时受到稀疏波等因素的影响,其速度明显偏低,故将其删去后再拟合,所获结果列于表 1

    图  7  底片判读结果
    Figure  7.  Results from experimental film
    表  1  不同工况下的破片及壳体速度
    Table  1.  Velocity of fragments and shell under different conditions
    Working condition Shell Fragments
    k v/(m·s-1) k v/(m·s-1)
    Tungsten fragments situated outside aluminum shell 0.126 983 0.126 983
    Tungsten fragments situated inside aluminum shell 0.143 1 115 0.110 858
    下载: 导出CSV 
    | 显示表格

    表 1中可以看出,破片内置时,钨破片的速度明显偏低,铝壳体的速度相对较高,且破片与壳体的动能之和也低于破片外置时的状态。这可能是因为此时钨破片直接与炸药接触,其上表面发生了相对严重的层裂(如图 6(b)所示),其碎裂部分的质量较大,且速度相对偏高,因此从整体上看,破片外置时,钨破片的完整性及速度均具有明显优势。轻质壳体和破片的完整性及初速均涉及其受力过程,故对应力波在这两种介质中的传播过程进行分析。

    为了简化分析过程,不考虑应力波在同一介质中的衰减,主要分析在不同介质间传播时的应力变化。分析过程中,首先需要描述爆轰产物作用过程的p(u)曲线以及铝合金壳体、钨合金破片的p(u)曲线,即波阵面压力与波阵面后粒子速度的关系曲线。若炸药以滑移爆轰方式驱动壳体及破片,则爆轰产物的p(u)曲线可根据(3)式计算[9]

    u=2Dγγ21[1(ppCJ)γ12γ]
    (3)

    式中:pCJ=ρ0D2/(γ+1),ρ0D分别为炸药密度和爆速;B炸药的多方指数γ=2.706[9]。铝合金和钨合金的p(u)曲线根据(4)式计算

    P=ρu(c0+λu)
    (4)

    式中:ρc0分别为金属材料的密度及声速,λ为本构方程参数。其具体数值列于表 2

    表  2  金属材料参数[9]
    Table  2.  Parameters of metals[9]
    Material ρ/(g·cm-3) c0/(km·s-1) λ
    Tungsten alloy 17.5 3.832 1.497
    Aluminum alloy 2.785 5.238 1.338
    下载: 导出CSV 
    | 显示表格

    图 8为应力波在不同介质中传播过程示意图,其中:3条实线分别表示B炸药爆轰产物的等熵线以及铝合金壳体、钨合金破片的p(u)曲线,爆轰产物的等熵线与金属材料p(u)曲线的交点坐标即为金属材料中的入射冲击波压力和波阵面后的粒子速度;金属介质间应力波的传播路径主要沿冲击材料在该状态点的反演线(如图 8中虚线所示),其计算表达式为[10]

    图  8  不同介质间应力波传播分析(①~⑥为界面编号)
    Figure  8.  Analysis of stress wave transmission betweendifferent materials (①-⑥:Interface No.)
    p=ρ(2u0u)[c0+λ(2u0u)]
    (5)

    式中:u0表示碰撞前的粒子速度。反演曲线与被冲击材料p(u)曲线的交点即表示碰撞后的状态。此外,由于空气的密度及声速均远低于金属材料,因此图 8中采用横坐标轴近似其p(u)曲线。

    图 8表 3可知, 破片置于轻质壳体外侧时:钨破片一端受到从壳体传入的冲击波作用,另一端受到空气反射的拉伸波作用,使其速度容易快速提升;而轻质壳体一端受到爆轰产物传入的冲击波作用,另一端受到钨破片反射的冲击波作用,两个作用方向相反,限制其速度过快提升,不会过多地占用炸药所转换的动能,同时没有拉伸波作用,轻质壳体也不容易快速破裂,在实际战斗部结构中可以增加对爆轰产物的约束,有利于对后期爆轰产物膨胀驱动能量的利用。

    表  3  材料界面处的状态参数
    Table  3.  State parameters of material interfaces
    Interface No. u/(km·s-1) p/GPa
    0.931 16.82
    0.385 29.70
    0.770
    0.310 23.30
    0.502 8.26
    1.004
    下载: 导出CSV 
    | 显示表格

    破片置于轻质壳体内侧时:虽然钨破片受到爆轰产物的冲击波及铝合金壳体反射的拉伸波作用,但从整体上看,其界面的平均粒子速度低于破片外置时的状态;而铝合金壳体由于外侧没有高阻抗介质阻挡,所受到的拉伸波作用可使其速度快速提升,但也易造成过早破裂。该工况在一些具有特定要求的战斗部中具有某些方面的优势,设计时可考虑在破片与装药之间增添一些性能较好的缓冲体,以提升破片的完整性;在轻质壳体与破片之间,也可适当增设缓冲材料,以减小轻质壳体与破片之间的速度差,不仅可以避免壳体过多地占用能量,也可减小壳体所受到的冲击力,延缓其破裂时间,从而增加对爆轰产物的约束时间,提升爆轰产物驱动能量的利用率。当然,本研究只针对冲击碰撞阶段的加速过程,未进行两个加速阶段的系统研究,也未涉及轻质壳体的动态破裂机理问题,还需要进一步深入研究。

    (1) 预制破片置于轻质壳体外侧时,初始应力波由低阻抗金属材料向高阻抗金属材料传播,使破片受到壳体传入的冲击波及空气传入的拉伸波作用,其速度易迅速提升,且破损轻微;而轻质壳体没有受到拉伸波作用,不易快速破裂,对后期爆轰产物膨胀驱动能量的利用较为有利。

    (2) 预制破片置于轻质壳体内侧时,初始应力波主要由高阻抗金属材料向低阻抗金属材料传播,虽然破片受到爆轰产物传入的冲击波及壳体反射的拉伸波作用,但其初速低于破片外置时的状态,且易发生破损,甚至有明显层裂现象;而铝合金壳体由于受到空气传入的拉伸波作用,其速度可快速提升,但也易造成其过早破裂,可能减少对爆轰产物的约束时间,降低炸药驱动能量的利用率。

  • 图  预制破片与轻质壳体的典型放置工况

    Figure  1.  Typical arrangements of prefabricated fragments and light shell

    图  圆筒试验中圆筒壁膨胀速度-时间历程

    Figure  2.  Time histories of expansionvelocity of cylinder wall

    图  滑移爆轰驱动破片示意

    Figure  3.  Sketch of sliding detonation-driven fragments

    图  破片飞散轨迹示意

    Figure  4.  Flying trajectory of fragments

    图  实验布局

    Figure  5.  Experiment layout

    图  实验所获底片

    Figure  6.  Experimental film

    图  底片判读结果

    Figure  7.  Results from experimental film

    图  不同介质间应力波传播分析(①~⑥为界面编号)

    Figure  8.  Analysis of stress wave transmission betweendifferent materials (①-⑥:Interface No.)

    表  1  不同工况下的破片及壳体速度

    Table  1.   Velocity of fragments and shell under different conditions

    Working condition Shell Fragments
    k v/(m·s-1) k v/(m·s-1)
    Tungsten fragments situated outside aluminum shell 0.126 983 0.126 983
    Tungsten fragments situated inside aluminum shell 0.143 1 115 0.110 858
    下载: 导出CSV

    表  2  金属材料参数[9]

    Table  2.   Parameters of metals[9]

    Material ρ/(g·cm-3) c0/(km·s-1) λ
    Tungsten alloy 17.5 3.832 1.497
    Aluminum alloy 2.785 5.238 1.338
    下载: 导出CSV

    表  3  材料界面处的状态参数

    Table  3.   State parameters of material interfaces

    Interface No. u/(km·s-1) p/GPa
    0.931 16.82
    0.385 29.70
    0.770
    0.310 23.30
    0.502 8.26
    1.004
    下载: 导出CSV
  • [1] 卢芳云, 蒋邦海, 李翔宇, 等.武器战斗部投射与毁伤[M].北京:科学出版社, 2013.
    [2] DRUMMOND W E.Explosive induced shock wave.Part Ⅰ.plane shock waves[J].Journal of Applied Physics, 2004, 28(12):1437-1441. https://www.researchgate.net/publication/224537576_Explosive_Induced_Shock_Waves_Part_I_Plane_Shock_Waves
    [3] RAM O, SADOT O.A simple constitutive model for prodicting the pressure histories developed behind rigid porous media impinged by shock wave[J].Journal of Fluid Mechanics, 2013, 718:507-523. doi: 10.1017/jfm.2012.627
    [4] 李卫平, 孙红, 张海丰.基于ALE方法的预制破片战斗部数值研究[J].弹箭与制导学报, 2012, 32(6):93-95. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-AGLU200708002004.htm

    LI W P, SUN H, ZHANG H F.Numerical simulation of prefabricated-fragment warheads exploding based on ALE method[J].Journal of Projectiles, Rockets, Missiles and Guidance, 2012, 32(6):93-95. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-AGLU200708002004.htm
    [5] MARK S, SCOTT I J.Dynamics of high sound-speed metal confiners driven by non-ideal high-explosive detonation[J].Combustion and Flames, 2015, 162:1857-1867. doi: 10.1016/j.combustflame.2014.12.007
    [6] MOSTERT F J, SNYMAN I M, OLIVIER M. An experimental technique to characterize the dynamic response of materials, or material combinations, to explosive blast[C]//26th International Symposium on Ballistics. Miami: International Ballistics Society, 2011: 1334-1337. http://adsabs.harvard.edu/abs/2011JAM....78e1008M
    [7] LINDSAY C M, BUTLER G C, RUMCHIK C G, et al.Increasing the utility of the copper cylinder expansion test[J].Propellants, Explosives, Pyrotechnics, 2010, 35(5):433-439. doi: 10.1002/prep.201000072
    [8] 沈飞, 王辉, 袁建飞, 等.铝含量对RDX基含铝炸药驱动能力的影响[J].火炸药学报, 2013, 36(3):50-53. http://www.cnki.com.cn/Article/CJFDTotal-BGXB201502014.htm

    SHEN F, WANG H, YUAN J F, et al.Influence of Al content on the driving ability of RDX-based aluminized explosives[J].Chinese Journal of Explosives and Propellants, 2013, 36(3):50-53. http://www.cnki.com.cn/Article/CJFDTotal-BGXB201502014.htm
    [9] 奥尔连科Л П. 爆炸物理学[M]. 孙承纬, 译. 北京: 科学出版社, 2011.
    [10] 经福谦, 陈俊祥.动高压原理与技术[M].北京:科学出版社, 2006.
  • 加载中
图(8) / 表(3)
计量
  • 文章访问数:  6650
  • HTML全文浏览量:  2856
  • PDF下载量:  286
出版历程
  • 收稿日期:  2017-07-06
  • 修回日期:  2017-07-20

目录

/

返回文章
返回