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

带隔板装药爆轰波马赫反射理论研究和数值模拟

潘建 张先锋 何勇 邓启斌

引用本文:
Citation:

带隔板装药爆轰波马赫反射理论研究和数值模拟

    作者简介: 潘建(1987—),男,博士研究生;
    通讯作者: 张先锋, lynx@mail.edu.cn
  • 基金项目:
    中央高校基本科研业务费专项项目 30916011305

  • 中图分类号: O381

Theoretical and numerical study on detonation wave Mach reflection in high explosive charge with waveshaper

    Corresponding author: Zhang Xianfeng, lynx@mail.edu.cn ;
  • CLC number: O381

  • 摘要: 基于三波理论和Whitham方法对带隔板装药爆轰波相互作用后发生的正规反射和非正规反射进行了理论分析,给出了爆轰波发生马赫反射时临界入射角和马赫杆增长角等参数的变化规律,提出了马赫杆高度的计算模型。基于凝聚炸药爆轰Jones-Wilkins-Lee(JWL)模型和冲击起爆的Lee-Tarver模型,利用有限元计算软件对带隔板装药爆轰波的传播过程进行了数值模拟。结果表明,发生马赫反射后,随着爆轰波的传播,马赫杆的高度不断增加。数值模拟结果与理论计算结果吻合较好,说明本文中采用的理论模型和数值模拟方法能够较准确地描述带隔板装药爆轰波马赫反射的传播过程。
  • 图 1  带隔板柱形装药结构及爆轰波示意图

    Figure 1.  Schematic diagram for cylindrical charge with waveshaper and detonation wave

    图 2  爆轰波马赫反射流场

    Figure 2.  Flow setup used to describe Mach reflection

    图 3  马赫反射临界入射角随多方指数的变化

    Figure 3.  Critical incident angle for the onset of Mach reflection varying with polytropic exponent

    图 4  马赫反射后爆轰波形位置与其垂直射线的关系简图

    Figure 4.  Relationship of detonation wave and its vertical line after Mach reflection

    图 5  马赫杆高度随入射角变化的计算结果

    Figure 5.  Calculated variation of Mach stem height with incident angle

    图 6  带隔板装药数值模拟的有限元模型

    Figure 6.  A finite element model for numerical simulation on charge with waveshaper

    图 7  采用不同的炸药状态方程模拟得到的爆轰波波形

    Figure 7.  Detonation waveform of explosive simulated by different equations of state

    图 8  不同时刻爆轰波的传播

    Figure 8.  Propagation of detonation wave at different times

    图 9  马赫反射发生后,随着爆轰波的传播,在轴线不同位置处压力的变化

    Figure 9.  Pressure at different positions of the axis after Mach reflection

    图 10  马赫反射发生后,随着爆轰波的传播,在轴线不同位置处密度的变化

    Figure 10.  Density at different positions of the axis after Mach reflection

    图 11  炸药Comp B中的马赫杆高度

    Figure 11.  Mach stem height in Comp B explosive

    图 12  炸药PBX9404中的马赫杆高度

    Figure 12.  Mach stem height in PBX9404 explosive

    图 13  炸药TNT中的马赫杆高度

    Figure 13.  Mach stem height in TNT explosive

    表 1  不同炸药装药中发生马赫反射时临界入射角的计算值与实验结果

    Table 1.  Calculated and experimental results of critical incident angles corresponding to Mach reflection in different explosive charges

    装药 ρ0/(g·cm-3) DCJ/(km·s-1) γ φⅠ, c/(°)
    实验 计算
    RDX 1.80 8.754 2.980 (44.5±2)[5] 44.70
    TNT 1.63 6.930 3.120 45.6[8] 44.12
    8321 1.70 8.212 2.838 43.72[3] 44.60
    PBX9501 1.83 8.802 2.097 56[6] 48.13
    下载: 导出CSV

    表 2  马赫参数的计算值与实验结果[7]的比较

    Table 2.  Calculated Mach parameters compared with experimental results[7]

    φ/(°) DCJ/(km·s-1) DM, exp/(km·s-1) DM, cal/(km·s-1) $\frac{{{D_{{\rm{M, exp}}}} - {D_{{\rm{M, cal}}}}}}{{{D_{{\rm{M, exp}}}}}}/\% $ z χexp/(°) χcal/(°) $\frac{{{\chi _{{\rm{exp}}}} - {\chi _{{\rm{cal}}}}}}{{{\chi _{{\rm{exp}}}}}}/\% $
    46.5 4.14 5.89 5.86 0.5 3.42 1.5 1.55 -3.2
    48.5 4.20 5.77 5.79 -0.3 3.22 2.27 2.13 6.6
    49.5 4.48 6.07 6.11 -0.7 3.13 2.80 2.45 14.3
    下载: 导出CSV
  • [1] Dowker J S. Quantum mechanics on group space and Huygens' principle[J]. Annals of Physics, 1971, 62(2):361-382. doi: 10.1016/0003-4916(71)90096-0
    [2] 汤明钧, 杨全中, 阎秀华.破甲弹中爆轰波相互作用的初步研究[J].南京理工大学学报:自然科学版, 1979(1):145-156.
    Tang Mingjun, Yang Quanzhong, Yan Xiuhua. Preliminary study of detonation wave interaction in HEAT[J]. Journal of Nanjing University of Science Technology:Nature Science, 1979(1):145-156.
    [3] 杨全中, 范宝春, 杨崇惠.带隔板装药中爆轰波的相互作用[J].南京理工大学学报:然科学版, 1982(2):23-30.
    Yang Quanzhong, Fan Baochun, Yang Chonghui. Interaction of the detonation waves in the explosive charge with the inserter[J]. Journal of Nanjing University of Science Technology:Nature Science, 1982(2):23-30.
    [4] Whitham G B. A new approach to problems of shock dynamics:Part Ⅰ:Two-dimension problem[J]. Journal of Fluid Mechanics, 1957, 2(2):145-171. doi: 10.1017/S002211205700004X
    [5] Dunne B B. Mach reflection of detonation waves in condensed high explosives:Ⅰ[J]. Physics of Fluids, 1961, 4(7):918-933. doi: 10.1063/1.1706425
    [6] Hull L M. Mach reflection of spherical detonation waves[C]//Short J M. Proceedings of the 10th International Symposium on Detonationm. Boston: Plenum Publishing Corporation, 1993: 12-16.
    [7] 张俊秀, 恽寿榕.凝聚炸药中爆轰波马赫反射的实验研究[J].爆炸与冲击, 1986, 6(3):208-213.
    Zhang Junxiu, Yun Shourong. Experimental investigations of the Mach relection in a condensed explosive[J]. Explosion and Shock Waves, 1986, 6(3):208-213.
    [8] Zhang Xianfeng, Huang Zhengxiang, Qiao Liang. Detonation wave propagation in double-layer cylindrical high explosive charge[J]. Propellants, Explosives, Pyrotechnics, 2011, 36(3):210-218. doi: 10.1002/prep.v36.3
    [9] 张先锋, 丁建宝, 赵晓宁.夹层装药作用过程的数字模拟[J].爆炸与冲击, 2009, 29(6):617-624. doi: 10.3321/j.issn:1001-1455.2009.06.011
    Zhang Xianfeng, Ding Jianbao, Zhao Xiaoning. Numberical simulation of duble layer shaped charge[J]. Explosion and Shock Waves, 2009, 29(6):617-624. doi: 10.3321/j.issn:1001-1455.2009.06.011
    [10] Murphy M J, Lee E L, Weston A M. Modeling shock initiation in composition B[C]//Short J M. Proceedings of the 10th International Symposium on Detonationm. Boston: Plenum Publishing Corporation, 1992: 965-970.
    [11] Grasso F, Paoli R. An analytical study of mach reflection in nonequilibrium steady flows[J]. Physics of Fluids, 1999, 11(10):3150-3167. doi: 10.1063/1.870172
    [12] Sternberg H M, Piacesi D. Interaction of oblique detonation waves with iron[J]. Physics of Fluids, 1966, 9(7):1307-1315. doi: 10.1063/1.1761845
    [13] Whithan G B. Liner and nonlinear waves[M]. New York:John Wiley & Sons Press, 1974:263-274.
    [14] Lambourn B D, Wright P W.Mach interaction of two plane detonation waves[C]//Proceeding of the Fourth International Symposium on Detonation. 1965: 142-152.
    [15] Lee E, Finger M, Collins W. JWL equation of state coefficients for high explosives[R]. Livermore: Lawrence Livermore Laboratory, 1973.
    [16] Davison L, Horie Y, Shahinpoor M. High-pressure shock compression of solids[M]. Springer Science & Business Media, 2012:123-143.
  • [1] 李志鹏曾贵玉刘兰蒋治海谯志强黄毅民 . TATB微观结构对炸药爆轰波传播性能的影响. 爆炸与冲击, 2009, 29(6): 665-668. doi: 10.11883/1001-1455(2009)06-0665-04
    [2] 陈军曾代朋孙承纬张震宇谭多望 . JB-9014炸药超压爆轰产物的状态方程. 爆炸与冲击, 2010, 30(6): 583-587. doi: 10.11883/1001-1455(2010)06-0583-05
    [3] 潘昊胡晓棉 . 钝感炸药的超压爆轰与冲击起爆过程数值模拟. 爆炸与冲击, 2006, 26(2): 174-178. doi: 10.11883/1001-1455(2006)02-0174-05
    [4] 王昌建徐胜利费立森 . 气相爆轰波绕射流场显示研究. 爆炸与冲击, 2006, 26(1): 27-32. doi: 10.11883/1001-1455(2006)01-0027-06
    [5] 张崇玉胡海波王翔刘宁文 . 铅飞层中斜冲击波对碰马赫反射行为实验研究. 爆炸与冲击, 2019, 39(4): 043102-1-043102-5. doi: 10.11883/bzycj-2017-0441
    [6] 缪玉松李晓杰王小红闫鸿浩陈翔 . 爆轰波碰撞的聚能效应. 爆炸与冲击, 2017, 37(3): 544-548. doi: 10.11883/1001-1455(2017)03-0544-05
    [7] 张先锋丁建宝赵晓宁 . 夹层聚能装药作用过程的数值模拟. 爆炸与冲击, 2009, 29(6): 617-624. doi: 10.11883/1001-1455(2009)06-0617-08
    [8] 任国武张崇玉张亚军洪仁楷陈永涛 . 金属铅在斜冲击波对碰加载下动态行为的数值模拟. 爆炸与冲击, 2018, 38(2): 397-403. doi: 10.11883/bzycj-2016-0177
    [9] 王昌建郭长铭徐胜利 . 气相爆轰波正反射激波加速研究. 爆炸与冲击, 2007, 27(2): 143-150. doi: 10.11883/1001-1455(2007)02-0143-08
    [10] 刘岩武丹王健平 . 低马赫数下斜爆轰波的结构. 爆炸与冲击, 2015, 35(2): 203-207. doi: 10.11883/1001-1455(2015)02-0203-05
    [11] 沈慧铭李伟兵王晓鸣李文彬董晓亮 . 圆柱壳体装药偏心多点起爆下破片速度的分布. 爆炸与冲击, 2017, 37(6): 1039-1045. doi: 10.11883/1001-1455(2017)06-1039-07
    [12] 张俊秀恽寿榕 . 凝聚炸药中爆轰波马赫反射的实验研究. 爆炸与冲击, 1986, 6(3): 208-213.
    [13] 易仰贤 . 空爆冲击波马赫反射近似计算. 爆炸与冲击, 1983, 3(2): 44-49.
    [14] 潘振华范宝春归明月张旭东 . 流动系统中爆轰波传播特性的数值模拟. 爆炸与冲击, 2010, 30(6): 593-597. doi: 10.11883/1001-1455(2010)06-0593-05
    [15] 刘玉存王建华陆春荣刘登程于雁武 . HMX粒度与爆轰波传播性能的关系研究. 爆炸与冲击, 2007, 27(6): 562-566. doi: 10.11883/1001-1455(2007)06-0562-05
    [16] 王保国张景林王作山刘玉存王建华 . 爆轰和爆炸冲击复合合成金刚石微粉得率的影响因素研究. 爆炸与冲击, 2006, 26(5): 429-433. doi: 10.11883/1001-1455(2006)05-0429-05
    [17] 张德良谢巍郭长铭胡湘渝 . 气相爆轰胞格结构和马赫反射数值模拟. 爆炸与冲击, 2001, 21(3): 161-167.
    [18] 夏昌敬周凯元 . 气相爆轰波在90矩形弯管中传播时胞格结构的演化. 爆炸与冲击, 2005, 25(2): 151-156. doi: 10.11883/1001-1455(2005)02-0151-06
    [19] 宁惠君王浩吴坛辉章猛华张成阮文俊 . 爆轰驱动异形杆在空气阻力作用下的动力学模拟. 爆炸与冲击, 2015, 35(4): 541-546. doi: 10.11883/1001-1455(2015)04-0541-06
    [20] 王爱峰赵伟姜宗林 . 斜爆轰的胞格结构及横波传播. 爆炸与冲击, 2010, 30(4): 349-354. doi: 10.11883/1001-1455(2010)04-0349-06
  • 加载中
图(13)表(2)
计量
  • 文章访问数:  355
  • HTML全文浏览量:  4
  • PDF下载量:  496
  • 被引次数: 0
出版历程
  • 收稿日期:  2014-11-17
  • 录用日期:  2015-03-09
  • 刊出日期:  2016-07-25

带隔板装药爆轰波马赫反射理论研究和数值模拟

    作者简介:潘建(1987—),男,博士研究生
    通讯作者: 张先锋, lynx@mail.edu.cn
  • 南京理工大学机械工程学院, 江苏 南京 210094
基金项目:  中央高校基本科研业务费专项项目 30916011305

摘要: 基于三波理论和Whitham方法对带隔板装药爆轰波相互作用后发生的正规反射和非正规反射进行了理论分析,给出了爆轰波发生马赫反射时临界入射角和马赫杆增长角等参数的变化规律,提出了马赫杆高度的计算模型。基于凝聚炸药爆轰Jones-Wilkins-Lee(JWL)模型和冲击起爆的Lee-Tarver模型,利用有限元计算软件对带隔板装药爆轰波的传播过程进行了数值模拟。结果表明,发生马赫反射后,随着爆轰波的传播,马赫杆的高度不断增加。数值模拟结果与理论计算结果吻合较好,说明本文中采用的理论模型和数值模拟方法能够较准确地描述带隔板装药爆轰波马赫反射的传播过程。

English Abstract

  • 在常规武器系统中,聚能装药战斗部扮演着越来越重要的角色,提高炸药能量利用率是战斗部设计者追求的目标之一。带隔板装药不仅结构简单而且可以在不增加装药量以及不改变装药长径比的条件下,提高其对药型罩的压垮做功能力。典型的带隔板装药结构如图 1所示,根据惠更斯原理[1],爆轰波经过隔板的绕射后将以O′点开始传播,以较小的入射角驱动药型罩形成聚能侵彻体。这类带隔板装药在驱动药型罩形成高速杆式射流、大长径比爆炸成型弹丸以及利用起爆控制技术形成多模毁伤元等方面应用非常广泛。

    图  1  带隔板柱形装药结构及爆轰波示意图

    Figure 1.  Schematic diagram for cylindrical charge with waveshaper and detonation wave

    早在20世纪70年代,汤明钧等[2]利用高速摄影测定了带隔板破甲弹中装药端面上爆轰波的爆速和爆压,得出马赫波的波速达到正常爆速的1.9倍,爆压达到正常CJ爆压的3倍。杨全中等[3]对带隔板装药爆轰波的相互作用进行了研究,并利用Whitham方法[4]计算了爆轰波传播过程中的波形和碰撞临界角等参数。B.B.Dunne等[5]指出了发生马赫反射的临界入射角取决于爆轰产物的多方指数,并给出了马赫杆压力等参数的计算方法。L.M.Hull等[6]计算了球面爆轰波的马赫反射参数,通过实验给出了PBX9501和PBXN-9等炸药的马赫反射参数。张俊秀等[7]进行了爆轰波马赫反射的实验研究,并根据修改的Whitham方法[4]计算了马赫反射的压力、爆速和增长角等参数。Zhang Xianfeng等[8]、张先锋等[9]分别利用理论计算和数值模拟分析了夹层装药爆轰波的传播过程,并基于质量守恒给出了爆轰波传播过程中马赫杆高度的计算模型。

    本文中,基于三波理论[3]和Whitham方法[4],对带隔板装药爆轰波相互作用的正规反射和马赫反射进行理论分析,并给出马赫杆高度的计算模型。在此基础上,利用AUTODYN有限元计算程序中的Lee-Tarver模型[10]对带隔板装药爆轰波传播过程进行数值模拟, 以期研究成果可为描述带隔板装药驱动典型介质提供理论基础和数据支持。

    • 2个等强度爆轰波的正规斜碰撞,可视为其中一侧爆轰波对刚性壁面的正规斜碰撞,当入射角φ较小时,流入界面I和界面R的粒子流的偏转角相等。入射角不断增大,达到马赫反射的临界入射角φⅠ, c时,偏转角θ2小于θ1。粒子流的不断积累迫使碰撞点P移动到Q点,便形成了图 2所示的马赫反射。基于F.Grasso等[11]关于爆轰波反射的描述,爆轰前反应区用带下标1的参数表示,反射波前反应区用带下标2的参数表示,超压爆轰反应区用带下标3的参数表示。

      图  2  爆轰波马赫反射流场

      Figure 2.  Flow setup used to describe Mach reflection

      假设炸药完全反应且满足CJ条件,区域(1)中的爆轰参数可由炸药的CJ参数、初始密度ρ0、多方指数γ和入射角获得。粒子流通过区域(1)发生偏转流入区域(2),偏转角和反射角满足下列关系:

      $ \tan {\theta _1} = \frac{{\left[ {Ma_1^2{{\sin }^2}\left( {{\varphi _Ⅱ} + \varepsilon } \right) - 1} \right]{{\tan }^{ - 1}}\left( {{\varphi _Ⅱ} + \varepsilon } \right)}}{{Ma_1^2\left[ {(\gamma + 1)/2 - {{\sin }^2}\left( {{\varphi _Ⅱ} + \varepsilon } \right)} \right] + 1}} $

      式中$M a_{1}=\sqrt{1+(\gamma+1)^{2} /\left(\gamma \tan \varphi_Ⅰ\right)^{2}}$为马赫数; 发生正规反射时,ε=θ1;发生马赫反射时,ε=θ1-θ2

      发生正规反射和非正规反射时, 式(1)中ε的表达式是不同的,若用发生正规反射时ε=θ1求解发生非正规反射时式子(1)的表达式,在0<φ+επ/2范围内必定无解。根据此条件可以获得发生马赫反射的临界入射角,不同炸药临界入射角的理论计算值和实验结果如表 1所示, 表中DCJ为炸药CJ爆轰速度。通过比较发现,两者比较接近,从而可以验证马赫反射临界入射角计算的正确性。通过计算还发现,发生马赫反射时的临界入射角与炸药多方指数成一一对应关系,随着多方指数的增大,发生马赫反射的临界入射角不断减小,如图 3所示,此结论与B.B.Dunne[5]给出的结果一致。

      装药 ρ0/(g·cm-3) DCJ/(km·s-1) γ φⅠ, c/(°)
      实验 计算
      RDX 1.80 8.754 2.980 (44.5±2)[5] 44.70
      TNT 1.63 6.930 3.120 45.6[8] 44.12
      8321 1.70 8.212 2.838 43.72[3] 44.60
      PBX9501 1.83 8.802 2.097 56[6] 48.13

      表 1  不同炸药装药中发生马赫反射时临界入射角的计算值与实验结果

      Table 1.  Calculated and experimental results of critical incident angles corresponding to Mach reflection in different explosive charges

      图  3  马赫反射临界入射角随多方指数的变化

      Figure 3.  Critical incident angle for the onset of Mach reflection varying with polytropic exponent

      大量理论和实验研究发现[5, 11-12],马赫杆是一个弯曲的曲面,交于对称轴P点,相切于爆轰波(I)Q点。在P点, 偏转角α=θ1-θ2=0°; 在Q点, α=θ1-θ2;马赫杆与对称轴切线的夹角从π/2变化到φ。由于碰撞点P处的压力和密度等都远大于CJ值,难以准确获取,本文中借助炸药爆轰产物的状态方程p=p(E, v)来求得此处的压力, 其中p为压力,E为炸药能量,v为比容。由H.M.Sternberg等[12]的定义,产生马赫反射时爆轰产物中将从三波点产生一条滑动线。此直线通过的地方,压力是连续的,但速度和密度非连续,因此,在p3=p2处存在一个三波点。

    • 基于Whitham方法[4, 13]获取马赫反射的增长角和高度等参数,并将马赫波简化为平面波。将爆轰波的每一元素沿每个微元射线管的传播当作固壁管道中的激波传播来处理,在二维情况下,爆轰波的位置和与其垂直的射线位置如图 4所示。爆轰波的传播过程可以看成是在变截面流管中的传播过程,当截面发生微小变化时,流动符合等熵流性质,得到马赫反射的A(z)关系式为:

      图  4  马赫反射后爆轰波形位置与其垂直射线的关系简图

      Figure 4.  Relationship of detonation wave and its vertical line after Mach reflection

      $ Af(z) = {A_0}f(1) = {\rm{const}} $

      式中:z=p/pCJ为爆轰波强度,且:

      $ f(z) = {z^{1/\gamma }}{(2z - 1)^{1/2}}B(z)C(z);\;\;\, C(z) = \exp \left\{ {{\gamma ^{ - 1/2}}{{\sin }^{ - 1}}\left[ {\frac{{(3 - \gamma )z - 2}}{{(\gamma + 1)z}}} \right]} \right\} $

      $ B(z) = {\left\{ {{{[2(\gamma - 1)(2z - 1)]}^{1/2}} + 2{{[(\gamma - 1)z + 1]}^{1/2}}} \right\}^{{{[2\gamma /(\gamma - 1)]}^{1/2}}}} $

      基于上述讨论,如图 4所示, 假设Δt时刻后,由爆轰波NS变成马赫波PQφ为爆轰波入射角,χ为三波点增长角,NQ点为不同时刻的三波点。设NS管道截面积为ACJPQ处管道截面积为AM,根据几何关系知:

      $ \tan {\varphi _Ⅰ} = \left( {\frac{{{D_{{\rm{CJ}}}}}}{{{D_{\rm{M}}}}} + \frac{{{A_{\rm{M}}}}}{{{A_{{\rm{CJ}}}}}}} \right){\left\{ {\left[ {1 - {{\left( {\frac{{{D_{{\rm{CJ}}}}}}{{{D_{\rm{M}}}}}} \right)}^2}} \right]\left[ {1 - {{\left( {\frac{{{A_{\rm{M}}}}}{{{A_{{\rm{CJ}}}}}}} \right)}^2}} \right]} \right\}^{ - 1/2}} $

      $ \frac{{\cos \chi }}{{\sin \left( {{\varphi _Ⅰ} - \chi } \right)}} = \frac{{{D_{\rm{M}}}}}{{{D_{{\rm{CJ}}}}}} $

      $ \frac{{\sin \chi }}{{\cos \left( {{\varphi _Ⅰ} - \chi } \right)}} = \frac{{{A_{\rm{M}}}}}{{{A_{{\rm{CJ}}}}}} $

      由式(5)知:

      $ \frac{{{A_{\rm{M}}}}}{{{A_{{\rm{CJ}}}}}} = \frac{{f(1)}}{{f(z)}} $

      参考B.D.Lambourn等[14]的研究,将爆速DMDCJ定义为:

      $ \frac{{{D_{\rm{M}}}}}{{{D_{{\rm{CJ}}}}}} = \frac{z}{{{{(2z - 1)}^{1/2}}}} $

      上述计算中未考虑爆轰波后稀疏波的影响,计算结果误差较大,对Whitham方法[4, 13]中的A(z)关系式进行修正:

      $ \frac{{{A_{{\rm{CJ}}}}}}{{{A_{\rm{M}}}}} = {(2z - 1)^{\frac{1}{{2k(z)}}}} $

      式中:k(z)是慢变函数,一般取常数, 通常由实验确定。

      根据张俊秀等[7]提供的入射角和爆速参数计算方法,通过方程(2)~(8),得到马赫波参数DMzχ的计算结果如表 2所示。利用Whitham方法[4, 13]计算获得的马赫波速度和三波点增长角与实验结果[7]吻合较好,从而验证此方法是可行的。

      φ/(°) DCJ/(km·s-1) DM, exp/(km·s-1) DM, cal/(km·s-1) $\frac{{{D_{{\rm{M, exp}}}} - {D_{{\rm{M, cal}}}}}}{{{D_{{\rm{M, exp}}}}}}/\% $ z χexp/(°) χcal/(°) $\frac{{{\chi _{{\rm{exp}}}} - {\chi _{{\rm{cal}}}}}}{{{\chi _{{\rm{exp}}}}}}/\% $
      46.5 4.14 5.89 5.86 0.5 3.42 1.5 1.55 -3.2
      48.5 4.20 5.77 5.79 -0.3 3.22 2.27 2.13 6.6
      49.5 4.48 6.07 6.11 -0.7 3.13 2.80 2.45 14.3

      表 2  马赫参数的计算值与实验结果[7]的比较

      Table 2.  Calculated Mach parameters compared with experimental results[7]

      同一装药中炸药的CJ爆速为常数,三波点增长角只与入射角和k(z)有关。根据带隔板装药爆轰波传播的特点,将入射角离散化,网格点(n, n)上的值可由网格点(n-1, n-1)上的值求得,马赫杆的高度可以定义为:

      $ \left\{ {\begin{array}{*{20}{l}} {{H_n} = \frac{{R\tan {\chi _{n - 1}}\tan {\varphi _{n - 1}} - {L_{n - 1}}\tan {\varphi _{n - 1}} + {H_{n - 1}}}}{{1 + \tan {\chi _{n - 1}}\tan {\varphi _{n - 1}}}}}\\ {{L_n} = \frac{{{L_{n - 1}}\tan {\chi _{n - 1}}\tan {\varphi _{n - 1}} + \left( {R - {H_{n - 1}}} \right)\tan {\varphi _{n - 1}}}}{{1 + \tan {\chi _{n - 1}}\tan {\varphi _{n - 1}}}}} \end{array}} \right. $

      式中:R为隔板半径;HnLn分别为马赫杆的高度和传播距离,n=1, 2, 3, …。

      n=1时,H0=1, L1=Rtanφ0。L.M.Hull[5]也认为,当入射角为临界入射角时,马赫杆高度为零或非常小。结合表 1的数据和方程(3)~(9),可以确定马赫杆高度。

      爆轰波绕射过隔板后,爆轰波不断向轴线方向汇聚,入射角不断增大。当入射角小于临界入射角时,发生正规反射,马赫杆高度为零。当入射角大于临界入射角时,发生马赫反射,马赫杆高度随着入射角的增大而增大,相同炸药中,随着k(z)的变化,马赫杆高度计算结果如图 5所示。马赫杆在不断增大的过程中,马赫杆上的压力和速度迅速降低,当马赫杆上的压力和速度降低到CJ值时,马赫杆消失,此时,马赫反射变成正常的CJ爆轰。

      图  5  马赫杆高度随入射角变化的计算结果

      Figure 5.  Calculated variation of Mach stem height with incident angle

    • 根据前述对带隔板装药爆轰波传播过程的分析,在保证隔板外侧爆轰波可以正常传播的条件下,隔板直径越大,调整爆轰波形状的能力越强,在装药结构中爆轰波马赫杆的形状就越明显。带隔板装药爆轰波的传播过程是一种多物质相互作用的大变形运动过程,因此本文中采用多物质耦合的Euler算法描述带隔板装药爆轰波的传播过程。图 6所示为装药的有限元模型。

      图  6  带隔板装药数值模拟的有限元模型

      Figure 6.  A finite element model for numerical simulation on charge with waveshaper

    • 凝聚炸药爆轰产物状态方程描述炸药爆轰达到CJ状态后爆轰产物中压力、温度和体积等物理量之间的关系,不仅体现炸药本身的做功能力,而且是研究爆炸力学各种问题的基础。本文中分别采用点火与增长模型的Lee-Tarver状态方程[10]和JWL状态方程[15]来描述爆轰波绕射过隔板后在炸药中的传播过程。炸药参数主要来自文献[10, 15]及AUTODYN中内置的参数。

      为了获取环形爆轰波形并保证透射波不起爆主装药,隔板厚度应足够大,这里隔板厚度选20 mm。隔板直径选42 mm,药柱直径选50 mm,使得隔板外侧装药厚度大于爆轰波正常传播时装药的临界直径。隔板材料为环氧树脂,其状态方程为Puff状态方程[16],主要参数分别为:ρ0=1.985 g/cm3A1=17.89 MPa,A2=201.2 MPa,A3=-6.053 MPa,Grüneisen系数为1.72,膨胀系数为0.25。

    • 针对上述讨论,对主装药分别采用JWL[15]和Lee-Travel[10]2种不同的状态方程进行描述,爆轰波波形的模拟结果如图 7所示。基于JWL模型[15]所形成的爆轰波的波阵面存在多段线组合形式,如图 7(a)所示。产生该结果的主要原因是:由于JWL模型为程序起爆模式,而炸药微元的爆轰采用预先设定时间起爆,JWL模型在描述带隔板装药爆轰(特别是隔板直径较大)时存在局限性,较难精确描述炸药微元的超压爆轰行为。基于炸药冲击起爆点火与增长的Lee-Travel模型[10]形成的爆轰波形较光滑,如图 7(b)所示。通过对比发现,Lee-Travel模型[10]能更真实地描述带隔板装药爆轰波的传播过程。

      图  7  采用不同的炸药状态方程模拟得到的爆轰波波形

      Figure 7.  Detonation waveform of explosive simulated by different equations of state

      对带隔板装药中的主副装药分别采用JWL[15]和Lee-Travel[10]2种不同的状态方程进行描述,爆轰波传播过程中不同时刻的波形如图 8所示。装药起爆后,爆轰波开始以球面波向前传播,经过隔板绕射后变成汇聚爆轰波。汇聚爆轰波在装药轴线处发生相互碰撞反射,当t=86 μs时,爆轰波已经发生马赫反射,随着爆轰波的传播, 马赫杆的高度不断增加,如图 8(b)(c)所示。随着爆轰波的传播,马赫杆的高度不断增加,但其压力和速度迅速降低,当马赫杆上的压力和爆轰速度等参数降低到CJ值时,马赫反射消失,马赫反射变成正常爆轰, 如图 8(d)所示。

      图  8  不同时刻爆轰波的传播

      Figure 8.  Propagation of detonation wave at different times

      为了进一步了解带隔板装药爆轰波传播过程中发生马赫发射后参数的变化,在装药轴线以及距轴线2.5、5.0、7.5和10.0 mm处分别设有一系列高斯点。以TNT炸药为例,通过设置高斯点可以获取爆轰产物的参数变化,当马赫反射发生后,随着爆轰波的传播,不同位置处压力和密度的变化见图 9~10。发生马赫反射时,装药轴线处(即马赫杆碰撞点处)的压力和密度均达到最大值,其中最大压力与CJ压力的比值达到4.3,最大密度与CJ密度的比为2.05,与Zhang Xianfeng等[8]对TNT装药发生马赫反射时最大压力和密度的理论计算值接近。随着爆轰波的传播,马赫杆碰撞点的压力和密度逐渐降低。通过对距轴线不同距离处设置高斯点发现,压力和密度先增大后减小,这说明马赫杆存在一定高度且是不断变化的。

      图  9  马赫反射发生后,随着爆轰波的传播,在轴线不同位置处压力的变化

      Figure 9.  Pressure at different positions of the axis after Mach reflection

      图  10  马赫反射发生后,随着爆轰波的传播,在轴线不同位置处密度的变化

      Figure 10.  Density at different positions of the axis after Mach reflection

      为了获取带隔板装药中马赫杆高度的变化,分别对带隔板装药中的Comp B、TNT和PBX9404炸药爆轰波的传播过程进行了数值模拟。随着爆轰波的传播,马赫杆高度的精确捕捉越来越难,其数值模拟结果和理论计算结果如图 11~13所示。

      图  11  炸药Comp B中的马赫杆高度

      Figure 11.  Mach stem height in Comp B explosive

      图  12  炸药PBX9404中的马赫杆高度

      Figure 12.  Mach stem height in PBX9404 explosive

      图  13  炸药TNT中的马赫杆高度

      Figure 13.  Mach stem height in TNT explosive

      当入射角较小时,马赫杆的高度为零,随着入射角不断增大并达到临界入射角时,马赫反射发生。此时,马赫杆碰撞点的压力和密度均达到最大值,随爆轰波的传播,马赫杆高度不断增加,马赫杆上的压力和密度不断减小。当马赫杆上的压力和密度等降到CJ值时,马赫反射消失,变成正常爆轰。

    • 通过对带隔板装药爆轰波马赫反射过程的理论分析以及数值模拟,得到的主要结论如下:

      (1) 基于三波理论和修正的Whitham方法,得到了马赫反射临界入射角和三波点增长角等参数,最终获得了马赫杆高度的计算模型。计算结果表明:爆轰波相互作用发生马赫反射时,马赫杆的高度为零;随着爆轰波的传播,马赫杆的高度不断增加,马赫杆上的压力和速度降低,当其降低到CJ值时,马赫反射消失。

      (2) 基于点火增长模型的Lee-Travel模型,带隔板装药爆轰波发生马赫反射后马赫杆高度的数值模拟结果与理论计算结果吻合较好,这说明利用Lee-Travel模型可以较准确地描述带隔板装药爆轰波的传播过程。

      (3) 带隔板装药技术是一种利用爆轰波相互作用产生超压爆轰的装药技术。马赫杆高度模型的建立不仅丰富和完善了带隔板装药爆轰波马赫反射的传播规律,而且可为描述带隔板装药驱动药型罩的过程提供理论基础。

参考文献 (16)

目录

    /

    返回文章
    返回