柱状装药预制破片缩比战斗部爆炸冲击波和破片的作用时序

夏冰寒 王金相 周楠 陈兴旺 卢孚嘉

引用本文:
Citation:

柱状装药预制破片缩比战斗部爆炸冲击波和破片的作用时序

    作者简介: 夏冰寒(1995-),男,硕士研究生,主要从事毁伤与防护理论与技术研究.E-mail: 117121011409@njust.edu.cn;
    通讯作者: 王金相, wjx@njust.edu.cn
  • 中图分类号: E932.4; O382.1

Blast Wave and Time Sequence of Prefabricated Fragments for Scaled Warhead with Cylindrical Charge

    Corresponding author: WANG Jinxiang, wjx@njust.edu.cn ;
  • CLC number: E932.4; O382.1

  • 摘要: 针对柱状装药的周向预制破片战斗部,结合无量纲分析方法和爆炸驱动理论,确定了影响破片和冲击波相遇位置的关键参数,给出了由缩比战斗部推广预测原型战斗部爆炸产生的破片冲击波作用时序的方法。采用ANSYS/LS-DYNA有限元软件进行数值模拟,对比验证了理论分析和数值试验结果,分析了战斗部缩比比例对冲击波和破片作用时序的影响。结果表明:缩比模型与原型战斗部爆炸产生的破片和冲击波的相遇位置之比和相遇时间之比主要取决于两模型的质量比,在不考虑破片速度衰减时,两模型中载荷相遇位置之比和相遇时间之比等于其质量比的0.33次方。受破片速度衰减影响,该方法仅适用于质量缩比不小于0.2的模型。
  • 图 1  战斗部模型示意图

    Figure 1.  Schematic diagram of the warhead model

    图 2  有限元模型

    Figure 2.  Finite element analysis model

    图 3  模型2中破片和冲击波的传播

    Figure 3.  Fragmentation and shock wave trajectory in Model 2

    图 4  模型2中破片和冲击波的传播距离与时间的关系

    Figure 4.  Propagation of blast wave and fragments as a function of time in air in Model 2

    表 1  破片和冲击波相遇距离问题中相关物理量及其单位和量纲

    Table 1.  Parameters and their units and dimensions related to the location of the two encounters

    ObjectParametersSymbolUnitDimension
    FragmentMassmfkgM
    ExplosiveMassmekgM
    Densityρekg∙m–3ML–3
    Chemical energy released per unit mass of explosiveEem2∙s–2L2T–2
    Expansion indexγe1SI
    AirInitial pressurepakg∙m–1∙s–2ML–1T–2
    Initial densityρakg∙m–3ML–3
    Adiabatic indexγa1SI
    下载: 导出CSV

    表 2  缩比战斗部尺寸

    Table 2.  Scaled warhead size

    Model r/cm h/cm d/cm me/g mf/gMass shrinkage
    ratio
    Dimension shrinkage
    ratio
    11.966 5.0230.126 100 62.920.10.464
    22.476 6.3280.159 200125.830.20.585
    32.835 7.2440.181 300188.750.30.669
    43.120 7.9720.200 400251.670.40.737
    53.93110.0460.252 800503.340.80.928
    64.23710.8220.2711 000629.151.01.000
    下载: 导出CSV

    表 3  TNT炸药材料参数及JWL状态方程参数

    Table 3.  Parameters of TNT material and JWL equation of state

    ρ/(kg∙m–3)D/(m∙s–1)pCJ/GPaE/(GJ∙m–3)A1/GPaB1/GPaR1R2ωV
    1 6406 93019.46.23093.094.4850.790.301
    下载: 导出CSV

    表 4  空气材料参数及状态方程参数

    Table 4.  Equation of state parameters of air

    ρ/(kg∙m–3)E/MPaC0/MPaC1C2C3C4C5C6
    1.250.25–0.100000.40
    下载: 导出CSV

    表 5  破片材料参数及状态方程参数

    Table 5.  Equation of state parameters of fragments

    ρ/(kg∙m–3)νσ/MPa${\dot\varepsilon } $
    7.830.31 0750.9
    下载: 导出CSV

    表 6  理论结果与仿真结果对比

    Table 6.  Comparison of theoretical and simulation results

    ModelMass reduction
    ratio
    ${\left( {\dfrac{m}{M}} \right)^{0.33}}$Meeting timeMeeting distance
    Value/µsReductionDeviation/%Value/cmReductionDeviation/%
    10.10.4641400.40313.1270.397 14.4
    20.20.5851880.531 9.2360.529 9.5
    30.30.6692290.642 4.0450.651 2.7
    40.40.7372700.763 3.5520.764 3.6
    50.80.9293340.941 1.3640.945 1.7
    61.01.000355 1.000 68 1.000
    下载: 导出CSV
  • [1] HU W, CHEN Z. Model-based simulation of the synergistic effects of blast and fragmentation on a concrete wall using the MPM [J]. International Journal of Impact Engineering, 2006, 32(12): 2066–2096. doi: 10.1016/j.ijimpeng.2005.05.004
    [2] LEPPÄNEN J. Concrete subjected to projectile and fragment impacts: modelling of crack softening and strain rate dependency in tension [J]. International Journal of Impact Engineering, 2006, 32(11): 1828–1841. doi: 10.1016/j.ijimpeng.2005.06.005
    [3] NYSTRÖM U, GYLLTOFT K. Numerical studies of the combined effects of blast and fragment loading [J]. International Journal of Impact Engineering, 2009, 36(8): 995–1005. doi: 10.1016/j.ijimpeng.2009.02.008
    [4] LEPPÄNEN J. Experiments and numerical analyses of blast and fragment impacts on concrete [J]. International Journal of Impact Engineering, 2005, 31(7): 843–860. doi: 10.1016/j.ijimpeng.2004.04.012
    [5] 张成亮, 朱锡, 侯海量, 等. 爆炸冲击波与高速破片对夹层结构的联合毁伤效应试验研究 [J]. 振动与冲击, 2014, 33(15): 184–188.
    ZHANG C L, ZHU X, HOU H L, et al. Tests for combined damage effect of blast waves and high-velocity fragments on composite sandwich plates [J]. Journal of Vibration and Shock, 2014, 33(15): 184–188.
    [6] 李茂, 朱锡, 侯海量, 等. 冲击波和高速破片对固支方板的联合作用数值模拟 [J]. 中国舰船研究, 2015, 10(6): 60–67. doi: 10.3969/j.issn.1673-3185.2015.06.009
    LI M, ZHU X, HOU H L, et al. Numerical simulation of steel plates subjected to the impact of both impact waves and fragments [J]. Chinese Journal of Ship Research, 2015, 10(6): 60–67. doi: 10.3969/j.issn.1673-3185.2015.06.009
    [7] 侯海量, 张成亮, 李茂, 等. 冲击波和高速破片联合作用下夹芯复合舱壁结构的毁伤特性 [J]. 爆炸与冲击, 2015, 35(1): 116–123. doi: 10.11883/1001-1455(2015)01-0116-08
    HOU H L, ZHANG C L, LI M, et al. Damage characteristics of sandwich bulkhead under the impact of shock and high-velocity fragments [J]. Explosion and Shock Waves, 2015, 35(1): 116–123. doi: 10.11883/1001-1455(2015)01-0116-08
    [8] LLOYD R. Conventional warhead systems physics and engineering design [M]. Reston: American Institute of Aeronautics and Astronautics, 1998.
    [9] 梁为民, 张晓忠, 梁仕发, 等. 结构内爆炸破片与冲击波运动规律试验研究 [J]. 兵工学报, 2009(Suppl 2): 223–227.
    LIANG W M, ZHANG X Z, LIANG S F, et al. Experimental research on motion law of fragment and shock wave under the condition of internal explosion [J]. Acta Armamentarii, 2009(Suppl 2): 223–227.
    [10] 安振涛, 王超, 甄建伟, 等. 常规弹药爆炸破片和冲击波作用规律理论研究 [J]. 爆破, 2012, 29(1): 15–18. doi: 10.3963/j.issn.1001-487X.2012.01.004
    AN Z T, WANG C, ZHEN J W, et al. Theoretical research on action law of fragment and shock wave of traditional ammunition explosion [J]. Blasting, 2012, 29(1): 15–18. doi: 10.3963/j.issn.1001-487X.2012.01.004
    [11] 郑红伟, 陈长海, 侯海量, 等. 破片尺寸对空爆冲击波及破片传播过程的影响仿真分析 [J]. 中国舰船研究, 2017, 12(6): 73–80. doi: 10.3969/j.issn.1673-3185.2017.06.011
    ZHENG H W, CHEN C H, HOU H L, et al. Simulation analysis of effects of single fragment size on air-blast wave and fragment propagation [J]. Chinese Journal of Ship Research, 2017, 12(6): 73–80. doi: 10.3969/j.issn.1673-3185.2017.06.011
    [12] 史志鑫, 尹建平, 王志军. 预制破片的形状对破片飞散性能影响的数值模拟研究 [J]. 兵器装备工程学报, 2017(12): 31–35. doi: 10.11809/scbgxb2017.12.008
    SHI Z X, YIN J P, WANG Z J. Numerical simulation of the influence of prefabricated fragments shape on fragment scattering performance [J]. Journal of Ordnance Equipment Engineering, 2017(12): 31–35. doi: 10.11809/scbgxb2017.12.008
    [13] 曾首义, 蒋志刚, 陈斌, 等. 冲击波与破片共同作用探讨 [C]//中国土木工程学会防护工程分会理事会暨学术会议, 2006: 263–267.
    ZENG S Y, JIANG Z G, CHEN B, et al. Discussion on the interaction between shock wave and fragmentation [C]//China Civil Engineering Society Protection Engineering Branch Council and Academic Conference, 2006: 263–267.
  • [1] 陈兴周兰伟李向东胡诤哲张瑞 . 破片式战斗部破片与冲击波相遇位置研究. 高压物理学报, 2018, 32(6): 065101-1-065101-9. doi: 10.11858/gywlxb.20180591
    [2] 陈海利蒋建伟门建兵 . 破片对带铝壳炸药的冲击起爆数值模拟研究. 高压物理学报, 2006, 20(1): 109-112 . doi: 10.11858/gywlxb.2006.01.021
    [3] 耿荻马天宝宁建国 . 定向战斗部爆炸驱动规律研究. 高压物理学报, 2013, 27(5): 685-690. doi: 10.11858/gywlxb.2013.05.005
    [4] 胡八一董庆东韩长生张海平 . TC4钛合金及40Cr钢破片中绝热剪切带的TEM分析. 高压物理学报, 1996, 10(1): 37-43 . doi: 10.11858/gywlxb.1996.01.006
    [5] 陈长海侯海量朱锡胡年明李典 . 破片式战斗部空中爆炸下冲击波与破片的耦合作用. 高压物理学报, 2018, 32(1): 015104-1-015104-9. doi: 10.11858/gywlxb.20170849
    [6] 章冠人 . 冲击波作用下材料的极化率. 高压物理学报, 1990, 4(3): 161-166 . doi: 10.11858/gywlxb.1990.03.001
    [7] 章冠人 . 冲击波作用下铁电体的击穿. 高压物理学报, 1995, 9(1): 1-4 . doi: 10.11858/gywlxb.1995.01.001
    [8] 王长利周刚马坤陈春林赵南冯娜 . 聚能装药水下爆炸冲击波载荷规律. 高压物理学报, 2017, 31(4): 453-461. doi: 10.11858/gywlxb.2017.04.014
    [9] 李海涛朱锡黄晓明牟金磊 . 水下爆炸冲击波作用下空化区域形成的特性研究. 高压物理学报, 2008, 22(2): 181-186. doi: 10.11858/gywlxb.2008.02.019
    [10] 杨世源金孝刚董玉斌 . 冲击波作用下多组分粉体的细化与均化研究. 高压物理学报, 2001, 15(1): 48-53 . doi: 10.11858/gywlxb.2001.01.007
    [11] 刘宏杰王伟力苗润吴世永王俊华 . 串联战斗部不同介质组合的隔爆能力. 高压物理学报, 2019, 33(1): 015104-1-015104-13. doi: 10.11858/gywlxb.20180585
    [12] 王新征张松林邹广平 . 内部短药柱爆炸作用下钢筒破裂特征的数值分析. 高压物理学报, 2010, 24(1): 61-66 . doi: 10.11858/gywlxb.2010.01.011
    [13] 张博李伟兵李文彬王晓鸣高旭东郑宇 . 偏心起爆战斗部随机破片数值仿真. 高压物理学报, 2012, 26(4): 442-448. doi: 10.11858/gywlxb.2012.04.013
    [14] 张绍兴李翔宇丁亮亮张震宇卢芳云 . 聚焦式战斗部破片轴向飞散控制技术. 高压物理学报, 2018, 32(1): 015103-1-015103-8. doi: 10.11858/gywlxb.20170512
    [15] 杨世全孙传杰钱立新卫剑锋 . 非金属壳体低附带战斗部实验与破片飞散分析. 高压物理学报, 2018, 32(4): 045103-1-045103-5. doi: 10.11858/gywlxb.20170573
    [16] 金孝刚王弘 . 冲击波处理45钢的微结构研究. 高压物理学报, 1993, 7(4): 254-259 . doi: 10.11858/gywlxb.1993.04.003
    [17] 王同权张若棋汤文辉肖亚斌 . 强流质子束辐照形成的冲击波. 高压物理学报, 2003, 17(2): 117-121 . doi: 10.11858/gywlxb.2003.02.007
    [18] 章冠人 . 冲击波和分形集团的凝聚. 高压物理学报, 1997, 11(4): 241-244 . doi: 10.11858/gywlxb.1997.04.001
    [19] 肖强强黄正祥顾晓辉 . 冲击波影响下的聚能射流侵彻扩孔方程. 高压物理学报, 2011, 25(4): 333-338 . doi: 10.11858/gywlxb.2011.04.008
    [20] 徐森刘大斌彭金华王建灵郭伟金朋刚贾宪政 . 药柱冲击波在有机玻璃中的衰减特性研究. 高压物理学报, 2010, 24(6): 431-437 . doi: 10.11858/gywlxb.2010.06.005
  • 加载中
图(4)表(6)
计量
  • 文章访问数:  229
  • 阅读全文浏览量:  218
  • PDF下载量:  0
出版历程
  • 收稿日期:  2019-05-17
  • 录用日期:  2019-05-28
  • 网络出版日期:  2020-01-21
  • 刊出日期:  2020-02-01

柱状装药预制破片缩比战斗部爆炸冲击波和破片的作用时序

    作者简介:夏冰寒(1995-),男,硕士研究生,主要从事毁伤与防护理论与技术研究.E-mail: 117121011409@njust.edu.cn
    通讯作者: 王金相, wjx@njust.edu.cn
  • 1. 南京理工大学瞬态物理国家重点实验室,江苏 南京 210094
  • 2. 南京森林警察学院,江苏 南京 210023
  • 3. 北京北方车辆集团有限公司,北京 100072

摘要: 针对柱状装药的周向预制破片战斗部,结合无量纲分析方法和爆炸驱动理论,确定了影响破片和冲击波相遇位置的关键参数,给出了由缩比战斗部推广预测原型战斗部爆炸产生的破片冲击波作用时序的方法。采用ANSYS/LS-DYNA有限元软件进行数值模拟,对比验证了理论分析和数值试验结果,分析了战斗部缩比比例对冲击波和破片作用时序的影响。结果表明:缩比模型与原型战斗部爆炸产生的破片和冲击波的相遇位置之比和相遇时间之比主要取决于两模型的质量比,在不考虑破片速度衰减时,两模型中载荷相遇位置之比和相遇时间之比等于其质量比的0.33次方。受破片速度衰减影响,该方法仅适用于质量缩比不小于0.2的模型。

English Abstract

  • 带壳装药战斗部爆炸时,产生空气冲击波的同时,还伴随着大量的高速弹片,共同对目标产生破坏作用。在早期研究该类战斗部的毁伤效果时,通常在远距离处只考虑破片的作用,而在近距离处将其解耦成近场冲击波的作用和高速破片群的作用[1-2]。近年来,国内外学者开始对破片和冲击波的联合毁伤效应展开了研究。Nyström等[3]研究了爆炸产生的冲击波和破片对混凝土的毁伤效果,分别对破片、冲击波的单独作用和联合作用展开分析,得到了联合作用毁伤效果大于两者单独作用毁伤效果的结论。Leppänen[4]利用试验和数值模拟方法研究了破片和冲击波对混凝土的联合毁伤作用,分析了破片密度、炸药量和起爆方式对毁伤的影响,得出的结论为破片和冲击波联合作用下的毁伤大于两种毁伤元单独作用时的毁伤。国内学者张成亮等[5]、李茂等[6]、侯海量等[7]开展了预制破片战斗部对一些组合防护结构的毁伤特性研究,分析了高速破片和爆炸冲击波对不同结构破坏模式的影响,结果表明,高速破片和冲击波的联合作用会加剧目标结构的破坏。

    由于冲击波与破片的运动规律大不相同,因此在不同的距离下战斗部爆炸产生的破片和冲击波作用在目标上的先后顺序也不同。确定冲击波与破片同时到达的距离,是研究两种毁伤效应耦合的重要基础。对于破片冲击波的作用时序问题,国内外学者均有研究,Nyström等[3]模拟了炸弹爆炸后破片和冲击波的相遇位置,从而对破片和冲击波的协同作用进行研究。Lloyd[8]分析了导弹战斗部爆炸后不同时间破片相对于冲击波的位置。梁为民等[9]在模拟爆腔内完成了模拟弹对目标靶板的爆炸破坏效应试验,研究了战斗部在结构内爆炸条件下破片和冲击波的运动演化过程,分析了不同比例距离和装药系数下,破片与冲击波的运动规律。安振涛等[10]对破片和冲击波的运动时序问题进行了理论分析,得出了冲击波和破片同时到达目标的距离,该距离大小与破片形状和单个破片质量关系不大。在对冲击波和破片的作用时序问题研究中,通常采用理论计算、数值仿真和试验等方式开展[9-13]。考虑到对于某些较大型战斗部,在仿真计算时会出现网格过多、计算时间过长的问题,实验研究时对实验场地要求很高,且操作困难,耗资较大,因此有必要进行缩比模型的相似律研究。目前相关的研究报道较少,且主要针对冲击波和目标结构的缩比,对于缩比战斗部爆炸产生的破片和冲击波的运动规律,尤其是两者作用时序的规律研究还不够充分。

    针对柱状装药预制破片缩比战斗部爆炸产生的破片和冲击波的传播过程与作用时序展开探究,考虑到实验难度较大,各项测量难以实施,故采用量纲分析方法,结合爆炸驱动理论与现有经验公式,对影响破片和冲击波作用时序的因素及其缩比后的相似准则进行分析,确定影响破片和冲击波相遇位置的关键参数,并结合数值模拟结果验证并分析战斗部缩比比例对破片和冲击波作用时序的影响。

    • 建立如图1所示的柱状装药预制破片战斗部模型,模型由炸药和破片2部分组成。炸药采用柱状TNT炸药,轴向密集排布540枚扇形破片,破片材料均为45钢,战斗部整体模型如图1(a)所示,其截面如图1(b)所示,其中rhd分别为装药半径、装药高度和破片厚度。由于炸药爆炸所产生的能量中破片的变形能约占总能量的1%[10],因此破片的变形、破坏和质量损失忽略不计,只考虑破片的动能、爆炸产物的动能和内能。决定破片和冲击波相遇位置的控制参数主要来自3个方面:破片的质量mf,炸药的质量me、装药密度ρe、单位质量炸药释放的化学能Ee、爆炸产物的膨胀指数γe,空气的初始压力pa、初始密度ρa、绝热指数γa。以上物理量的单位与量纲如表1所示。

      图  1  战斗部模型示意图

      Figure 1.  Schematic diagram of the warhead model

      ObjectParametersSymbolUnitDimension
      FragmentMassmfkgM
      ExplosiveMassmekgM
      Densityρekg∙m–3ML–3
      Chemical energy released per unit mass of explosiveEem2∙s–2L2T–2
      Expansion indexγe1SI
      AirInitial pressurepakg∙m–1∙s–2ML–1T–2
      Initial densityρakg∙m–3ML–3
      Adiabatic indexγa1SI

      表 1  破片和冲击波相遇距离问题中相关物理量及其单位和量纲

      Table 1.  Parameters and their units and dimensions related to the location of the two encounters

      若以R表示破片和冲击波的相遇距离,那么存在函数关系式

      $ R = f\left( {{m_{\rm{f}}};m_{\rm{e}},{\rho _{\rm{e}}},{E_{\rm{e}}},{\gamma _{\rm{e}}};{p_{\rm{a}}},{\rho _{\rm{a}}},{\gamma _{\rm{a}}}} \right) $

      上述物理量以L、M、T为基本量纲,可取meρeEe为基本量,则式(1)可化为下面的无量纲关系

      $\frac{R}{{{\rho _{\rm{e}}}^{ - 1/3}{m_{\rm{e}}}^{ - 1/3}}} = f\left( {\frac{{{m_{\rm{f}}}}}{m};1,1,1,{\gamma _{\rm{e}}};\frac{{{p_{\rm{a}}}}}{{{\rho _{\rm{e}}}{E_{\rm{e}}}}},\frac{{{\rho _{\rm{a}}}}}{{{\rho _{\rm{e}}}}},{\gamma _{\rm{a}}}} \right)$

      若采用同种炸药在空气中做缩比模型实验,那么有6个有关的控制参数与原模型相同,即

      $\left( {{\rho _{\rm{e}}},{E_{\rm{e}}},{\gamma _{\rm{e}}};{p_{\rm{a}}},{\rho _{\rm{a}}},{\gamma _{\rm{a}}}} \right) = \rm const$

      则无量纲函数关系式可简化为

      $R = {m_{\rm e}^{1/3}}f\left( {\frac{{{m_{\rm{f}}}}}{m}} \right)$

      式(4)即为基于量纲分析得到的相遇距离的定性关系式,可知破片和冲击波的相遇距离取决于炸药总质量以及破片和炸药的质量比。

    • 在1.1节中使用量纲分析法得到的破片和冲击波相遇位置的定性关系式还不足以明确两载荷相遇位置和相关参数的具体关系,本节将结合爆炸驱动理论和相关经验公式对其具体形式进行推导。

      对破片,依据Gurney公式[9]可知最大破片速度

      ${v_0} = G\sqrt {{\alpha }/({1 + 0.5\alpha })} $

      式中:G为炸药的格林参数;α为战斗部破片总质量与装药量之比,对于同一战斗部的不同缩比模型,模型中的α应为常数。故不同模型中的破片速度应相等。

      考虑到破片的速度衰减,破片飞行时间tf与距离R的关系为[10]

      ${t_{\rm{f}}} = 250\frac{{m_{\rm{f}}^{1/3}}}{{{v_0}}}\left[ {\exp \left( {\frac{{0.004R}}{{m_{\rm{f}}^{1/3}}}} \right) - 1} \right]$

      对冲击波,由于破片的存在会在一定程度上削减冲击波的强度,依据能量法可以得到削减后的等效TNT当量mbe与装药时的TNT药量m的关系[10]

      ${m_{\rm{be}}} = \frac{{0.347 + 0.653\beta }}{{2 - \beta }}{m_{\rm{e}}}$

      式中:β为战斗部的装填系数,$\beta = {m_{\rm{e}}}/\left( {{m_{\rm{e}}} + {m_{\rm{f}}}} \right)$。式(7)可改写为

      ${m_{\rm{be}}} = \mu {m_{\rm{e}}}$

      对于指定的战斗部,由于各缩比模型的α相等,故µ为常数。

      冲击波波阵面传播到距离R处的时刻ts可由经验公式[9]改进得到

      ${t_{\rm{s}}} = 4.56 \times {10^{ - 4}} {R^{1.76}} m_{\rm{be}}^{ - 0.25}$

      要求得破片和冲击波同时到达的距离R,令tf = ts,可得

      $250\frac{{m_{\rm{f}}^{1/3}}}{{{v_0}}}\left[ {\exp \left( {\frac{{0.004R}}{{m_{\rm{f}}^{1/3}}}} \right) - 1} \right] = 4.56 \times {10^{ - 4}} {R^{1.76}} m_{{\rm{be}}}^{ - 0.25}$

      通过式(10)直接求出Rm的关系较为困难,由于近场环境下破片速度衰减较小,因此若不考虑破片速度衰减,式(10)可改写为

      ${R}/{{v_0}} = 4.56 \times {10^{ - 4}} {R^{1.76}} m_{\rm{be}}^{ - 0.25}$

      如前所述,v0为常数,则式(11)可改写为

      $R = \gamma m_{\rm{be}}^{0.33}$

      忽略破片速度衰减,视破片为匀速运动,可以得到破片和冲击波的相遇时间t

      $t = {R}/{{{v_0}}} ={{\gamma m_{\rm{be}}^{0.33}}}/{{{v_0}}}$

      前面已指出,同一原型战斗部的不同缩比模型中v0μ均为常数,而γv0的函数,故γ也为常数。那么由式(12)、式(13)可知,破片和冲击波的相遇距离和时间与等效TNT当量的0.33次方成正比。由此可推断出质量为m的缩比模型中破片冲击波相遇时间t1和距离R1与质量M的原模型中破片冲击波相遇时间t2和位置R2的关系为

      $\frac{{{t_1}}}{{{t_2}}} = \frac{{{R_1}}}{{{R_2}}} = \frac{{m_{\rm{be}}^{0.33}}}{{M_{\rm{be}}^{0.33}}} = {\left( {\frac{m}{M}} \right)^{0.33}}$

    • 为研究不同缩比尺寸下的工况,针对图1所示的战斗部,保持战斗部各部分尺寸(rhd)比例不变,对整体战斗部模型进行缩比。结合前文可知,破片和冲击波的相遇距离主要取决于炸药与破片的质量,故按照0.8∶1、0.4∶1、0.3∶1、0.2∶1、0.1∶1的总质量比设计包括原模型在内共6个模型,具体尺寸列于表2

      Model r/cm h/cm d/cm me/g mf/gMass shrinkage
      ratio
      Dimension shrinkage
      ratio
      11.966 5.0230.126 100 62.920.10.464
      22.476 6.3280.159 200125.830.20.585
      32.835 7.2440.181 300188.750.30.669
      43.120 7.9720.200 400251.670.40.737
      53.93110.0460.252 800503.340.80.928
      64.23710.8220.2711 000629.151.01.000

      表 2  缩比战斗部尺寸

      Table 2.  Scaled warhead size

    • 采用ANSYS/LS-DYNA非线性动力有限元分析程序,考虑到整体结构的对称性,建立不同缩比质量下的战斗部的1/8模型,如图2所示。

      图  2  有限元模型

      Figure 2.  Finite element analysis model

      数值模型由炸药、空气和破片3部分组成,均采用8节点Solid164三维实体单元建模,其中炸药、空气单元使用多物质ALE算法,破片采用Lagrange网格建模,破片与炸药和空气材料间采用流固耦合算法。在对称面上设立对称边界,对空气边界设定透射边界。模型采用cm-g-µs单位制。起爆方式为中心线起爆。

      炸药采用HIGH_EXPLOSIVE_BURN本构模型,对爆轰产物的膨胀采用JWL状态方程描述

      $p = A\left( {1 - \frac{\omega }{{{R_1}V}}} \right){{\rm e}^{ - {R_1}V}} + B\left( {1 - \frac{\omega }{{R_2V}}} \right){{\rm e}^{ - {R_2}V}} + \frac{{\omega E}}{V}$

      式中:p为爆轰压力,ABR1R2ω为试验确定的参数,e为指定函数,V为相对体积,E为初始内能。计算中所采用的炸药参数见表3[11]。其中,ρ为装药密度,D为爆轰速度,pCJ为爆轰波阵面的压力。

      ρ/(kg∙m–3)D/(m∙s–1)pCJ/GPaE/(GJ∙m–3)A1/GPaB1/GPaR1R2ωV
      1 6406 93019.46.23093.094.4850.790.301

      表 3  TNT炸药材料参数及JWL状态方程参数

      Table 3.  Parameters of TNT material and JWL equation of state

      空气采用NULL材料模型及LINEAR_POLY_NOMIAL状态方程描述

      $p = {C_0} + {C_1}\mu + {C_2}{\mu ^2} + {C_3}{\mu ^3} + \left( {{C_4} + {C_5}\mu + {C_6}{\mu ^2}} \right)E$

      式中:$\mu = {\rho }/{{{\rho _0}}} - 1$C0C1C2C3C4C5C6为试验定义的常数,E为单位体积的能量。空气材料参数见表4[11]

      ρ/(kg∙m–3)E/MPaC0/MPaC1C2C3C4C5C6
      1.250.25–0.100000.40

      表 4  空气材料参数及状态方程参数

      Table 4.  Equation of state parameters of air

      破片采用MAT _PLASTIC_KINEMATIC的钢材料模型,其应变率由Cowper-Symonds模型描述

      $ {\sigma _{\rm{d}}} = \left( {{\sigma _0} + \frac{{E{E_{\rm{h}}}}}{{E - {E_{\rm{h}}}}}} \right)\left[ {1 + {{\left( {\frac{\dot\varepsilon }{D}} \right)}^{1/n}}} \right] $

      式中:σd为动态屈服强度;σ0为静态屈服强度;E为弹性模量,取E = 210 GPa;Eh为硬化模量,取Eh = 319 MPa;εp为有效塑性应变;${\dot\varepsilon } $为等效塑性应变率;Dn为常数,对于低碳钢,通常取n = 5,D = 400 s–1。材料失效模型采用最大等效塑性应变失效准则。破片材料参数见表5[12],其中ν为泊松比。

      ρ/(kg∙m–3)νσ/MPa${\dot\varepsilon } $
      7.830.31 0750.9

      表 5  破片材料参数及状态方程参数

      Table 5.  Equation of state parameters of fragments

    • 从各工况的数值仿真结果中可以发现,各个模型中破片和冲击波的传播规律大致相同,图3图4给出了模型2中各阶段的破片和冲击波传播过程,主要分为4个阶段:第1阶段(Ⅰ)中破片领先于正后方的冲击波,此时冲击波在驱动破片加速的同时绕流过破片;第2阶段(Ⅱ)破片速度趋于平稳,冲击波汇集在破片前方,速度由快至慢衰减;第3阶段(Ⅲ)破片追赶上冲击波;第4阶段(Ⅳ)破片赶超冲击波后领先于冲击波。

      图  3  模型2中破片和冲击波的传播

      Figure 3.  Fragmentation and shock wave trajectory in Model 2

      图  4  模型2中破片和冲击波的传播距离与时间的关系

      Figure 4.  Propagation of blast wave and fragments as a function of time in air in Model 2

    • 表6给出了各个模型中破片和冲击波的相遇时间、相遇距离与${\left( {{m / M}} \right)^{0.33}}$的关系。由表6结果可知,随着模型的减小,破片和冲击波的相遇距离也会随之缩减,这是因为破片的速度保持不变,冲击波速度变小,导致破片追赶上冲击波的时间和距离也随之缩短。对比各结果中的相遇时间缩比、相遇距离缩比与${\left( {{m / M}} \right)^{0.33}}$的关系,可以看到,质量缩比在1.0~0.2范围内时,理论与仿真结果符合较好,但误差会随着缩比模型的缩小而增大。这主要是由于在有限元计算中采用的ALE算法考虑了实际的破片和空气的相互作用,而随着模型的缩小,破片迎风面积与质量之比却在增大,因此破片速度衰减带来的影响也越来越大。结合以上数据,可以认为该方法适用于战斗部质量缩比不小于0.2的模型。

      ModelMass reduction
      ratio
      ${\left( {\dfrac{m}{M}} \right)^{0.33}}$Meeting timeMeeting distance
      Value/µsReductionDeviation/%Value/cmReductionDeviation/%
      10.10.4641400.40313.1270.397 14.4
      20.20.5851880.531 9.2360.529 9.5
      30.30.6692290.642 4.0450.651 2.7
      40.40.7372700.763 3.5520.764 3.6
      50.80.9293340.941 1.3640.945 1.7
      61.01.000355 1.000 68 1.000

      表 6  理论结果与仿真结果对比

      Table 6.  Comparison of theoretical and simulation results

    • 缩比模型与原型战斗部爆炸产生的破片和冲击波的相遇位置之比和相遇时间之比取决于两模型的质量比,在不考虑破片速度衰减时,两模型中的相遇位置之比和相遇时间之比等于其质量比的0.33次方。通过有限元仿真验证了理论的有效性,同时由于破片速度衰减的影响,该方法适用于质量缩比不小于0.2的模型。

参考文献 (13)

目录

    /

    返回文章
    返回