• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊
高级检索 E-mail Alert

破片撞击起爆柱面带壳装药的临界速度修正判据

王昕 蒋建伟 王树有 门建兵

引用本文:
Citation:

破片撞击起爆柱面带壳装药的临界速度修正判据

    作者简介: 王昕(1990-), 女, 博士研究生;
    通讯作者: 蒋建伟, bitjjw@bit.edu.cn
  • 中图分类号: O383

Critical detonation velocity calculation model of cylindrical covered charge impacted by fragment

    Corresponding author: JIANG Jianwei, bitjjw@bit.edu.cn ;
  • CLC number: O383

  • 摘要: 为获得适用于柱面带壳装药的冲击起爆修正判据,以Picatinny工程判据为基础加入修正项进行修正。采用AUTODYN-3D软件对破片撞击柱面带钢壳的B炸药进行数值计算,获得了破片入射角、装药曲率半径对炸药临界起爆速度的影响规律;通过拟合得到修正项表达式,建立了考虑破片入射角、柱壳装药形状函数的炸药起爆临界速度修正判据。判据计算值与实验数据和数值计算值吻合较好,该判据能较好的预测柱形带壳装药的冲击起爆条件。
  • 图 1  破片对柱壳装药作用剖面

    Figure 1.  Tungsten fragment impact cylindrical shell charge at collision point

    图 2  钨球与柱壳装药作用过程的物理及离散化网格模型

    Figure 2.  Physical model and discrete model of tungsten fragment impact cylindrical covered charge

    图 3  钨球以不同入射角撞击柱面带壳装药应力图

    Figure 3.  Stress of tungsten fragment impact cylindrical shell charge with different incidence angles

    图 4  柱壳和平面带壳装药内相同位置观测点压力时程曲线

    Figure 4.  Histories of pressure at same point in cylindrical charge and plate charge

    图 5  f(θ)随入射角度正弦值的变化曲线

    Figure 5.  Relation between f(θ) and sinθ

    图 6  f(h/r)随壳体厚度与曲率半径比值的变化关系

    Figure 6.  Relation between f(h/r) and h/r

    图 7  f(h/r)平均值随壳体厚度与曲率半径比值变化关系

    Figure 7.  Relation between average value of f(θ) and sinθ

    表 1  Comp-B炸药材料参数

    Table 1.  Material parameters of Comp-B explosive

    I/μs-1 b a x G1/(Pa-2·s-1) c d y G2 e g z
    44 0 0.01 4.0 414×10-16 0.222 0.667 2.0 0 0 0 0
    下载: 导出CSV

    表 2  破片、壳体材料模型

    Table 2.  Material model of fragment and casing

    部件 材料 状态方程 强度模型 失效应变
    壳体 Steel 4340 Linear Johnson-Cook Geometric Strain
    破片 Tungsten Shock Johnson-Cook Geometric Strain
    下载: 导出CSV

    表 3  临界起爆速度的数值模拟与理论值对比

    Table 3.  Comparison of critical initiation velocities between simulation and theory

    m/g h/mm θ/(°) vcr/(m·s-1) ε/%
    数值模拟 理论
    3 5 0 2 915 3 074.4 5.2
    4 5 0 2 660 2 545.5 4.5
    下载: 导出CSV

    表 4  各工况下炸药起爆的破片临界速度

    Table 4.  Critical detonation velocity under various conditions

    θ/(°) vcr/(m·s-1)
    r=40 mm r=60 mm r=75 mm r=100 mm r=200 mm r=∞
    0 2 826 2 838 2 828 2 840 2 868 2 915
    15 2 918 2 912 2 915 2 920 2 940 2 955
    30 3 069 3 073 3 084 3 078 3 088 3 092
    45 3 396 3 412 3 407 3 423 3 437 3 541
    55 3 832 3 830 3 835 3 845 3 837 3 843
    下载: 导出CSV

    表 5  不同质量破片临界起爆速度数值模拟与修正判据对比

    Table 5.  Comparison of critical initiation velocities between simulation and rectified criterion value

    工况 m/g h/mm r/mm θ/(°) vcr/(m·s-1) ε/%
    数值模拟 修正判据
    1 4 6 50 0 2 657 2 528 4.59
    2 4 6 75 30 2 928 2 730 6.76
    3 5 6 40 0 2 181 2 205 1.10
    4 5 6 40 15 2 237 2 254 0.76
    5 5 6 40 30 2 414 2 381 1.36
    6 4 5 40 0 2 217 2 182 1.60
    7 4 5 40 30 2 438 2 357 3.44
    下载: 导出CSV
  • [1] 陈海利, 蒋建伟, 门建兵.破片对带铝壳炸药的冲击起爆数值模拟研究[J].高压物理学报, 2006, 20(1):109-112. doi: 10.3969/j.issn.1000-5773.2006.01.021
    CHEN Haili, JIANG Jianwei, MEN Jianbing. Numerical simulation of fragment impacting on charge with aluminum shell[J].Chinese Journal of High Pressure Physics, 2006, 20(1):109-112. doi: 10.3969/j.issn.1000-5773.2006.01.021
    [2] 李小笠, 屈明, 路中华, 等.三种破片对带壳炸药冲击起爆能力的数值分析[J].弹道学报, 2009, 21(4):72-75.
    LI Xiaoli, QU Ming, LU Zhonghua, et al. Numerical analysis of impact initiation ability of three kinds of fragments on shell explosive[J]. Journal of Ballistics, 2009, 21(4):72-75.
    [3] 童宗保, 王金相, 彭楚才, 等.预制破片对屏蔽炸药冲击引爆研究[J].科学技术与工程.2013, 14(7):173-177. DOI:10.3969/j.issn.1671-1815.2014.07.038
    TONG Zongbao, WANG Jinxiang, PENG Chucai, et al. Study on shock initiation of shielded explosive by prefabricated fragment[J]. Science Technology and Engineering, 2013, 14(7):173-177. DOI:10.3969/j.issn.1671-1815.2014.07.038
    [4] HELD M. Initiation phenomenon with shaped charge jets[C]//Proceedings of the 9th International Detonation Symposium. Portland, Oregon, USA, 1989: 1416-1426.
    [5] ROSLUND L A. Initiation of warhead fragments Ⅰ: Normal impacts[R]. Naval Surface Weapons Center, White Oak, 1973.
    [6] LLOYD R M. Conventional warhead systems physics and engineering design[M]. Virginia:American Insititute of Aeronautics and Astronautics, Inc., 1998:502-504.
    [7] 章冠华, 陈大年.凝聚炸药起爆动力学[M].北京:国防工业出版社, 1991.
    [8] 张先锋, 赵有守, 陈惠武.射弹冲击引爆带壳炸药临界条件[J].弹道学报, 2006, 18(4)57-59. doi: 10.3969/j.issn.1004-499X.2006.04.016
    ZHANG Xianfeng, ZHAO Youshou, CHEN Huiwu. The critical condition of shelled explosive initiated by projectile[J]. Journal of Ballistics, 2006, 18(4)57-59. doi: 10.3969/j.issn.1004-499X.2006.04.016
    [9] HELD M. Initiation criteria of high explosive attacked with projectiles or jet densities[J]. Propellants, Explosives, Pyrotechnics, 1996, 21(5):235-237. doi: 10.1002/(ISSN)1521-4087
    [10] 方青, 卫玉章, 张克明.射弹倾斜撞击带盖板炸药引发爆轰的条件[J].爆炸与冲击, 1997, 17(2):153-158.
    FANG Qing, WEI Yuzhang, ZHANG Keming. On the projectile oblique-impact initiation conditions for explosive covered with a plate[J]. Explosion and Shock Waves, 1997, 17(2):154-158.
    [11] 陈卫东, 张忠, 刘家良, 等.破片对屏蔽炸药冲击起爆的数值模拟与分析[J].兵工学报, 2009, 30(9):1187-1191. doi: 10.3321/j.issn:1000-1093.2009.09.007
    CHEN Weidong, ZHANG Zhong, LIU Jialiang, et al. Numerical simulation and analysis of shock initiation of shielded explosive by fragment[J]. Acta Armamentarii, 2009, 30(9):1187-1191. doi: 10.3321/j.issn:1000-1093.2009.09.007
    [12] 杜茂华, 王伟力, 黄勇, 等.舰载超近程反导弹药冲击引爆战斗部的研究[J].工程爆破, 2012, 18(2):14-17.DOI:10.3969/j.issn.1006-7051.2012.02.004.
    DU Maohua, WANG Weili, HUANG Yong, et al. Research on impacting and igniting warhead by super close-in anti-missile ammunition on board[J]. Engineering Blasting, 2012, 18(2):14-17. DOI:10.3969/j.issn.1006-7051.2012.02.004.
    [13] 王昕, 蒋建伟, 王树有, 等.钨球对柱面带壳装药的冲击起爆数值模拟研究[J].兵工学报, 2017, 38(8):1498-1505.DOI:10.3969/j.issn.1000-1093.2017.08.006.
    WANG Xin, JIANG Jianwei, WANG Shuyou, et al. Numerical simulation on the initiation of cylindrical covered charge impacted by tungsten sphere fragment[J]. Acta Armamentarii, 2017, 38(8):1498-1505. DOI:10.3969/j.issn.1000-1093.2017.08.006.
    [14] 梁争峰, 袁宝慧.破片撞击起爆屏蔽B炸药的数值模拟和实验[J].火炸药学报, 2006, 29(1):5-9. doi: 10.3969/j.issn.1007-7812.2006.01.002
    LIANG Zhengfeng, YUAN Baohui. Numerical simulation and experimental study of the initiation of shielded composition b impacted by fragment[J]. Chinese Journal of Explosives and Propellants, 2006, 29(1):5-9. doi: 10.3969/j.issn.1007-7812.2006.01.002
    [15] LEE E L, TARVER C M. Phenomenological model of shock initiation in heterogeneous explosive[J]. Physics of Fluids, 1980, 23(12):2362. doi: 10.1063/1.862940
    [16] IZADIFARD R A, FOROUTAN M. Blastwave parameters assessment at different altitude using numerical simulation[J]. Turkish Journal of Engineering and Enviromental Sciences, 2010, 34(1):25-41.
  • 加载中
图(7)表(5)
计量
  • 文章访问数:  37
  • HTML全文浏览量:  226
  • PDF下载量:  16
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-08-09
  • 录用日期:  2017-10-18
  • 刊出日期:  2019-01-05

破片撞击起爆柱面带壳装药的临界速度修正判据

    作者简介:王昕(1990-), 女, 博士研究生
    通讯作者: 蒋建伟, bitjjw@bit.edu.cn
  • 北京理工大学爆炸科学与技术国家重点实验室, 北京 100081

摘要: 为获得适用于柱面带壳装药的冲击起爆修正判据,以Picatinny工程判据为基础加入修正项进行修正。采用AUTODYN-3D软件对破片撞击柱面带钢壳的B炸药进行数值计算,获得了破片入射角、装药曲率半径对炸药临界起爆速度的影响规律;通过拟合得到修正项表达式,建立了考虑破片入射角、柱壳装药形状函数的炸药起爆临界速度修正判据。判据计算值与实验数据和数值计算值吻合较好,该判据能较好的预测柱形带壳装药的冲击起爆条件。

English Abstract

  • 采用高速破片对来袭弹药进行可靠引爆是有效拦截的重要手段。其问题实质则是破片对带壳装药的冲击起爆问题。学者们针对这一问题已进行了广泛的研究,大多以破片冲击平面带壳装药的简化模型为基础开展分析,包括不同材质、形状破片[1-3]对不同厚度壳体、炸药冲击起爆的影响规律,以及建立的诸多半经验的起爆判据,如著名的Held经验判据[4]、Jacobs-Roslund经验准则[5]、Picatinny冲击引爆解析计算式[6]及非均相炸药起爆判据[7]等。也有学者对经典判据的适用条件进行补充修正,包括张先锋等[8]依据Held[9]起爆裸装药的射弹临界速度计算模型,得到射弹起爆带壳装药判据[9];方青等[10]以Jacobs-Roslund准则为基础,将该准则向斜碰撞上扩展;陈卫东等[11]则推导了不同材质破片冲击起爆屏蔽装药的理论判据等。这些判据能较好预测平面带壳装药的冲击起爆。而实际应用中导弹战斗部多为圆柱形,杜茂华等[12]提到战斗部壳体曲率会影响毁伤元对带壳装药的冲击起爆。已有的数值模拟初步结果也表明,柱面带壳装药与平面带壳装药的冲击起爆特性存在差异[13]。基于平面带壳装药模型的判据已无法准确描述破片对实际弹药的冲击起爆,因此研究适用于柱面带壳装药的临界起爆速度判据非常必要。

    本文中采用AUTODYN-3D数值计算软件,对球形破片以不同入射角起爆柱面带壳装药(简称柱壳装药)问题进行数值计算。在对数值模拟数据分析的基础上,通过建立入射角、柱壳装药形状系数的修正因子,得到基于Picatinny工程判据的柱壳装药临界起爆速度修正判据,并与已有的实验数据和数值模拟结果进行对比,验证该判据的普适性。

    • 现役弹药形状多为圆柱形,将破片对来袭弹药战斗部的相互作用简化为破片对柱壳装药的冲击起爆。图 1所示为破片撞击柱壳装药的物理模型图,其中柱壳装药壳体厚度为h,装药曲率半径为r;破片着速为v,入射角为θ

      图  1  破片对柱壳装药作用剖面

      Figure 1.  Tungsten fragment impact cylindrical shell charge at collision point

      破片对柱壳装药的冲击起爆与破片、装药壳体、炸药以及两者交汇条件等因素相关,主要影响参量如下:

      (1) 破片:直径d,长度l,密度ρp,弹性常数Epγp,屈服极限Yp,形状系数N1

      (2) 壳体:厚度h,密度ρt,弹性常数Etγt,屈服极限Yt,形状系数N2

      (3) 炸药:炸药密度ρe,单位质量炸药释放的化学能Ee,膨胀系数γe

      (4) 弹靶相互作用条件:破片入射角θ,攻角β

      因此,带壳装药的临界起爆速度vcr可写成以下函数表达式:

      $ {v_{{\rm{cr}}}} = f\left( {d,l,{\rho _{\rm{p}}},{E_{\rm{p}}},{\gamma _{\rm{p}}},{Y_{\rm{p}}},{N_1};h,{\rho _{\rm{t}}},{E_{\rm{t}}},{\gamma _{\rm{t}}},{Y_{\rm{t}}},{N_2};{\rho _{\rm{e}}},{E_{\rm{e}}},{\gamma _{\rm{e}}};\theta ,\beta } \right) $

      选取hρeEe为基本量,根据量纲分析π定理,对上式进行无量纲化,得到:

      $ \frac{{{v_{{\rm{cr}}}}}}{{\sqrt {{E_{\rm{e}}}/{\rho _{\rm{e}}}} }} = f\left( {\frac{d}{h},\frac{l}{h},\frac{{{\rho _{\rm{p}}}}}{{{\rho _{\rm{e}}}}},\frac{{{E_{\rm{p}}}}}{{{E_{\rm{e}}}}},{\gamma _{\rm{p}}},\frac{{{Y_{\rm{p}}}}}{{{E_{\rm{e}}}}},{N_1};\frac{{{\rho _{\rm{t}}}}}{{{\rho _{\rm{e}}}}},\frac{{{E_{\rm{t}}}}}{{{E_{\rm{e}}}}},{\gamma _{\rm{t}}},\frac{{{Y_{\rm{t}}}}}{{{E_{\rm{e}}}}},{N_2};{\gamma _e};\theta ,\beta } \right) $

      在破片和柱壳装药材料不变的情况下,上式可简化为:

      $ \frac{{{v_{{\rm{cr}}}}}}{{\sqrt {{E_{\rm{e}}}/{\rho _{\rm{e}}}} }} = f\left( {\frac{d}{h},\frac{l}{h},{N_1},{N_2},\theta ,\beta } \right) $

      因本文中所研究的破片质量较小,忽略破片尺寸变化引起的冲击起爆阈值的差异,并假定破片与柱壳装药交汇时攻角β=0°,则临界起爆速度值只与破片入射角θ和柱壳装药形状系数N2这2个无量纲函数相关。

      无量纲形状系数N2表征柱壳装药的形状对冲击起爆阈值的影响,可用壳体厚度h和曲率半径r组合表示:

      $ {N_2} = \frac{h}{r} $

      上述分析表明,在壳体厚度保持不变的情况下,对柱壳装药的冲击起爆产生影响的破片入射角θ和装药曲率半径r为相互独立的影响参量。下文中基于此开展修正判据的建立。

    • Picatinny工程判据基于屏蔽炸药的敏感系数Kf、屏蔽板厚度h、破片质量m这3个参数建立破片临界起爆速度计算式:

      $ {v_0} = {\left[ {\frac{{{K_{\rm{f}}}\exp \left( {5.37h/{m^{1/3}}} \right)}}{{{m^{2/3}}\left( {1 + 3.3h/{m^{1/3}}} \right)}}} \right]^{1/2}} $

      该式能够较好预测破片垂直撞击平面带壳装药的起爆阈值。本文研究在该工程判据基础上添加入射角和基于曲率半径的柱壳装药形状系数修正项,以期预测破片起爆柱壳装药的临界速度。

      在应用中提出如下假设:

      (1) Picatinny判据适用于小质量规则破片引发的炸药冲击起爆[14],不考虑破片形状变化引起的炸药冲击起爆差异;

      (2) 壳体为钢壳,暂不考虑柱壳装药壳体材质变化引起的炸药冲击起爆差异;

      (3) 破片运动速度与装药轴线垂直,仅考虑撞击位置的变化引起的入射角的改变;

      (4) 忽略破片与柱壳装药作用时,装药绕轴心转动引起的临界起爆速度误差。

      基于上述假设建立判据修正项,记v0为破片以入射角θ=0°起爆平面带壳装药(r=∞)时的临界速度;定义v(N2, θ)为破片以入射角θ≠0°起爆同等壳体厚度、装药曲率半径r≠∞的柱壳装药临界速度。

      定义δ表示破片起爆炸药的临界速度增量,称为修正因子,即:

      $ \delta = \frac{{v\left( {{N_2},\theta } \right) - {v_0}}}{{{v_0}}} $

      变换式(6)得:

      $ v\left( {{N_2},\theta } \right) = \left( {1 + \delta } \right){v_0} $

      式中:v0采用Picatinny工程判据公式计算,v(N2, θ)则由数值计算获得。

      采用入射角θ和柱壳装药形状系数N2来反映破片与柱壳装药交汇条件和柱壳装药结构对临界起爆速度的影响,即δ=f(θ, N2),因θN2相互独立,则δ可表示为:

      $ \delta = f\left( \theta \right)f\left( {{N_2}} \right) $

      $ \delta = f\left( \theta \right)f\left( {h/r} \right) $

      式中:f(θ)为特定装药曲率半径r*条件下,不同入射角θ时的δ值,即:

      $ f\left( \theta \right) = \delta \left( {{r^ * },\theta } \right) $

      保持壳体厚度不变,f(h/r)即为特定碰撞角度θ*条件下,不同曲率半径r时的δ值与特定曲率半径r*δ值之比,即:

      $ f\left( {h/r} \right) = \frac{{\delta \left( {r,{\theta ^ * }} \right)}}{{\delta \left( {{r^ * },{\theta ^ * }} \right)}} $

      以上为基于Picatinny工程判据应用在柱壳装药上的修正判据建立方法。

    • 式(10)~(11)为构造修正因子函数的方法,为获得具体表达式,针对破片与柱壳装药的作用过程进行数值模拟,采用数据拟合方法得到修正因子函数表达式。

      数值模拟采用AUTODYN-3D软件,图 2(a)为破片撞击柱壳装药轴对称数值模型,依据典型预制破片战斗部结构,选取破片为质量3g的钨球,柱壳装药,长度80 mm,壳体厚度6 mm,装填Comp-B炸药。图 2(b)为钨球、壳体和炸药的离散化网格模型。首先开展网格收敛性及优化工作,考虑CPU机时及计算精度,网格数量选取为2 048、2 800和19 200,采用Lagrange方法开展计算。

      图  2  钨球与柱壳装药作用过程的物理及离散化网格模型

      Figure 2.  Physical model and discrete model of tungsten fragment impact cylindrical covered charge

      为描述炸药在冲击作用下的起爆过程,其状态方程采用Lee-Tarver[15]状态方程:

      $ \frac{{{\rm{d}}F}}{{{\rm{d}}t}} = I{\left( {1 - F} \right)^b}{\left( {\mu - a} \right)^x} + {G_1}{\left( {1 - F} \right)^c}{F^d}{p^y} + {G_2}{\left( {1 - F} \right)^e}{F^g}{p^z} $

      式中:F为燃烧质量份数,它在模拟爆轰过程中控制炸药化学能的释放;参数a是临界压缩度参数;参数Ix为点火量冲击强度及持续函数;参数G1d控制点火后早期增长函数;参数G2z为高压反映率相关函数。表 1为Comp-B炸药材料模型参数。

      I/μs-1 b a x G1/(Pa-2·s-1) c d y G2 e g z
      44 0 0.01 4.0 414×10-16 0.222 0.667 2.0 0 0 0 0

      表 1  Comp-B炸药材料参数

      Table 1.  Material parameters of Comp-B explosive

      破片和壳体选用能较好描述材料大应变、高应变率及高温状态的Johnson-Cook强度模型,材料强度模型、状态方程和侵蚀准则列于表 2。参数均取自AUTODYN标准数据库。

      部件 材料 状态方程 强度模型 失效应变
      壳体 Steel 4340 Linear Johnson-Cook Geometric Strain
      破片 Tungsten Shock Johnson-Cook Geometric Strain

      表 2  破片、壳体材料模型

      Table 2.  Material model of fragment and casing

      材料参数的正确与否与计算结果正确性直接相关。采用经典的Picatinny工程判据对已建立的数值模型行验证,表 3为2种质量钨球垂直撞击起爆平面带壳B炸药的临界速度数值模拟与理论值的对比。可看出计算值与理论值误差约5%,即认定数值算法的正确性。

      m/g h/mm θ/(°) vcr/(m·s-1) ε/%
      数值模拟 理论
      3 5 0 2 915 3 074.4 5.2
      4 5 0 2 660 2 545.5 4.5

      表 3  临界起爆速度的数值模拟与理论值对比

      Table 3.  Comparison of critical initiation velocities between simulation and theory

      为得到柱壳装药受破片撞击后的起爆特性,设计30种工况开展计算,即钨球以入射角θ=0°, 15°, 30°, 45°, 55°分别撞击装药曲率半径r=40, 60, 75, 100, 200 mm, ∞时的柱壳装药。图 3所示为典型装药曲率下钨球以v=2 826 m/s,不同θ撞击柱壳装药t=8 μs时刻应力图。可以看出,当v相同时,柱壳装药中的初始传入冲击波及压力随θ的增加而减小,θ=0°时炸药已完全爆轰,而其余入射角条件下无法发生稳定爆轰。而破片撞击起爆不同曲率半径的柱壳装药的过程与平板装药的起爆基本一致,本文中主要从临界起爆速度角度来说明柱壳装药与平面带壳装药的差异。

      图  3  钨球以不同入射角撞击柱面带壳装药应力图

      Figure 3.  Stress of tungsten fragment impact cylindrical shell charge with different incidence angles

      采用升降法获得表 4中炸药起爆的临界破片速度,对比各临界起爆速度可得出如下结论:

      θ/(°) vcr/(m·s-1)
      r=40 mm r=60 mm r=75 mm r=100 mm r=200 mm r=∞
      0 2 826 2 838 2 828 2 840 2 868 2 915
      15 2 918 2 912 2 915 2 920 2 940 2 955
      30 3 069 3 073 3 084 3 078 3 088 3 092
      45 3 396 3 412 3 407 3 423 3 437 3 541
      55 3 832 3 830 3 835 3 845 3 837 3 843

      表 4  各工况下炸药起爆的破片临界速度

      Table 4.  Critical detonation velocity under various conditions

      (1) 入射角θ对柱壳装药冲击起爆影响较大。临界起爆速度随θ的增大而增加,装药曲率r=40 mm条件下,θ=55°时的起爆速度较θ=0°时增加35.6%,较平面带壳装药增加31.5%。

      (2) 装药曲率半径r对柱壳装药的冲击起爆有一定影响。除个别点外,临界起爆速度随r增加基本呈现非线性增大,入射角θ=0°条件下,r=∞时起爆速度较r=40 mm时提高3.2%。

      在平面和柱壳装药(r=∞、r=40 mm)内距壳体炸药边界4、8、12 mm处设置观测点,编号分别为1、2、3。将破片以入射角θ=0°起爆带壳装药时各观测点的压力进行对比,分析平面与柱壳装药冲击起爆阈值的差异。图 4v=3 000 m/s时观测点的压力时程曲线。冲击波到达平面带壳装药的时间早于柱壳装药,但相同位置处柱壳装药内的压力值高于平面装药。

      图  4  柱壳和平面带壳装药内相同位置观测点压力时程曲线

      Figure 4.  Histories of pressure at same point in cylindrical charge and plate charge

      定性分析柱壳装药冲击起爆速度低于平面带壳装药的原因。破片撞击柱壳装药时,会在壳体和炸药界面处发生波系的反射和透射,在炸药中传入透射波。由于柱壳装药曲率半径小于平面带壳装药,传入的波系会因曲率的存在而发生汇聚。相同撞击条件下、炸药内同一位置处的压力值高于平面带壳装药,因此,起爆柱壳装药易于起爆平面带壳装药。

      根据式(6)、(10)~(11)将表 4中临界起爆速度进行归一化处理,以获得柱壳装药临界起爆速度修正项,修正项的获取按文献[16]的方法进行拟合。当r=∞(平面带壳装药)时,令f(h/r)=1,此时比较入射角对临界起爆速度的影响,即δ=f(θ)。图 5所示为f(θ)随入射角正弦sinθ的变化关系,可以看出,其f(θ)随sinθ的增大呈现指数增加。当入射角度大于55°时破片将发生跳飞,临界起爆速度趋于无限大。进一步通过最小二乘法得出f(θ)随sinθ的变化关系:

      $ f\left( \theta \right) = - 0.018 + 0.012{{\rm{e}}^{\frac{{\sin \theta }}{{0.245}}}}\;\;\;\;\;\;0 < \theta < {55^ \circ } $

      图  5  f(θ)随入射角度正弦值的变化曲线

      Figure 5.  Relation between f(θ) and sinθ

      图 6~7分别为不同入射角θ时的f(h/r)和其平均值随壳体厚度与装药曲率半径比值(h/r)的变化关系。由图 7可以看出f(h/r)随h/r的增加而减小,但减小量在2%以内,说明曲率半径对柱壳装药的临界起爆速度有一定影响但影响不大。进一步得出f(h/r)随壳体厚度与装药曲率半径比值的变化关系,采用最小二乘法得到拟合公式为:

      $ f\left( {h/r} \right) = 0.981 + 0.019{{\rm{e}}^{\frac{h}{{ - 0.051\;6r}}}}\;\;\;\;\;\;0 < h/r < 0.15 $

      图  6  f(h/r)随壳体厚度与曲率半径比值的变化关系

      Figure 6.  Relation between f(h/r) and h/r

      图  7  f(h/r)平均值随壳体厚度与曲率半径比值变化关系

      Figure 7.  Relation between average value of f(θ) and sinθ

      综上,得到考虑破片入射角和壳体形状系数的临界起爆速度修正因子的δ的计算公式(9)、(13)~(14),将系数应用于式(5),则得到基于Picatinny工程判据的临界起爆速度修正模型:

      $ \begin{array}{*{20}{c}} {{v_0} = {{\left[ {\frac{{{K_{\rm{f}}}\exp \left( {5.37h/{m^{1/3}}} \right)}}{{{m^{2/3}}\left( {1 + 3.3h/{m^{1/3}}} \right)}}} \right]}^{1/2}}\\ \left[ {1 + \left( { - 0.018 + 0.012{{\rm{e}}^{\frac{{\sin \theta }}{{0.245}}}}} \right)\left( {0.981 + 0.019{{\rm{e}}^{\frac{h}{{ - 0.051\;6r}}}}} \right)} \right]}\\ {0 < \theta < {{55}^ \circ },\;\;\;\;0 \le h/r < 0.15} \end{array} $

    • 由于式(15)为特定破片质量条件下的装药临界起爆速度计算判据,为检验判据是否具有普适性,与文献已公布实验结果及部分计算数据进行对比。

      首先将判据修正值与文献[14]公布的实验数据进行对比,质量为4.65 g破片撞击起爆平面带壳装药的实验值在2 231~2 522 m/s之间,修正判据计算的起爆阈值为2 347 m/s,修正判据计算值与实验数据吻合较好。

      进一步开展不同质量破片撞击不同壳体厚度柱壳装药的数值模拟计算。表 5为各工况下的临界起爆速度计算值及判据计算值。可以看出,计算值与判据修正值误差在7%以内,可认定修正的修正判据普适性良好。该判据同样适用于其他材质破片对柱壳装药的冲击起爆阈值的计算。

      工况 m/g h/mm r/mm θ/(°) vcr/(m·s-1) ε/%
      数值模拟 修正判据
      1 4 6 50 0 2 657 2 528 4.59
      2 4 6 75 30 2 928 2 730 6.76
      3 5 6 40 0 2 181 2 205 1.10
      4 5 6 40 15 2 237 2 254 0.76
      5 5 6 40 30 2 414 2 381 1.36
      6 4 5 40 0 2 217 2 182 1.60
      7 4 5 40 30 2 438 2 357 3.44

      表 5  不同质量破片临界起爆速度数值模拟与修正判据对比

      Table 5.  Comparison of critical initiation velocities between simulation and rectified criterion value

    • 采用数值模拟的方法,获得了破片以不同入射角起爆不同曲率半径柱壳装药的临界速度。通过对数值模拟结果的分析,建立了基于柱壳装药冲击起爆的修正判据,结论如下:

      (1) 入射角对柱壳装药冲击起爆影响较大,临界起爆速度随入射角的增大而非线性增加,装药曲率r=40 mm条件下,θ=55°时的起爆速度较θ=0°时增加35.6%;装药曲率对柱壳装药的冲击起爆有一定影响,临界起爆速度随曲率半径增加基本呈现非线性增大,入射角θ=0°条件下,r=∞时起爆速度较r=40 mm时提高3.2%。

      (2) 通过引入包含入射角θ和柱壳装药形状函数h/r的修正因子,对Picatinny工程判据进行修正,建立了适用于起爆柱面带壳装药临界破片速度的修正判据,判据校验表明修正判据与实验值与数值计算值误差小于7%,能较好的预测柱壳装药的冲击起爆。

参考文献 (16)

目录

    /

    返回文章
    返回