• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

狭长直管约束条件水下电爆炸所产生的气泡运动和界面射流

张桂夫 朱雨建 李元超 杨基明

李如江, 韩宏伟, 孙素杰, 刘天生. 包覆板材料为陶瓷时平板装药的防护性能[J]. 爆炸与冲击, 2014, 34(1): 47-51. doi: 10.11883/1001-1455(2014)01-0047-05
引用本文: 张桂夫, 朱雨建, 李元超, 杨基明. 狭长直管约束条件水下电爆炸所产生的气泡运动和界面射流[J]. 爆炸与冲击, 2015, 35(5): 609-616. doi: 10.11883/1001-1455(2015)05-0609-08
Li Ru-jiang, Han Hong-wei, Sun Su-jie, Liu Tian-sheng. Ballistic resistance capabilities of explosive reactive armors encapsulated by ceramic layers[J]. Explosion And Shock Waves, 2014, 34(1): 47-51. doi: 10.11883/1001-1455(2014)01-0047-05
Citation: Zhang Gui-fu, Zhu Yu-jian, Li Yuan-chao, Yang Ji-ming. Bubble and jet induced by underwater wire explosion in a narrow tube[J]. Explosion And Shock Waves, 2015, 35(5): 609-616. doi: 10.11883/1001-1455(2015)05-0609-08

狭长直管约束条件水下电爆炸所产生的气泡运动和界面射流

doi: 10.11883/1001-1455(2015)05-0609-08
基金项目: 国家自然科学基金项目(11102204)
详细信息
    作者简介:

    张桂夫(1990—), 男, 博士研究生

    通讯作者:

    朱雨建, yujianrd@ustc.edu.cn

  • 中图分类号: O381

Bubble and jet induced by underwater wire explosion in a narrow tube

  • 摘要: 以高速摄影为主要手段,揭示直管中爆炸诱导气泡和射流的典型演变过程,并测试爆炸深度和爆炸能量对该现象的影响。研究发现直管中爆炸诱导的表面射流分为光滑和粗糙的两段,这区别于自由表面射流的形态;爆炸气泡的发展经历一个先膨胀再坍缩的过程,其中封闭坍缩以气泡顶部形成内向射流为特征。表面射流速度主要来自爆炸早期短时间内气泡膨胀赋予水体的冲量,且整体上与起爆能量成正相关,而与爆炸深度成反相关;用准一维的简化模型能够很好地描述它们之间的依赖关系,计算结果不仅在趋势上与实验结果一致,数值上也能很好吻合。
  • 平板装药, 又称爆炸反应装甲(explosive reactive armour, ERA), 首先由M.Held[1]发明, 由于具有防护效益高、重量轻、装卸方便、成本低等优点, 而广泛应用于现代装甲车辆的防护。传统的爆炸反应装甲结构为两层对称或非对称的钢板中间夹层炸药组成的三明治结构。当射流高速撞击具有一定倾角的爆炸反应装甲时, 会引爆中间炸药层, 钢板在高压、高速爆轰产物作用下, 沿其法线方向高速运动, 对射流形成连续或断续切割, 使射流受到严重干扰, 从而降低射流对主装甲的侵彻能力。

    对传统平板装药的大量研究[2-5]表明, 它爆炸后会产生大量的高速破片(通常达500~1 500m/s), 有可能对战场己方人员和车辆造成损伤, 因此需寻找新的材料替代传统的钢板。H.Kaufmamnn等[6]对低密度包覆板材料(包括铝、陶瓷、玻璃钢等)进行了实验研究, 发现低密度材料具有和钢一样的防护效能, 且附带损伤明显降低。P.J.Hazell等[7]运用X射线研究了陶瓷和硼硅玻璃等包覆板材料与射流的作用, 发现脆性包覆板在爆炸产物作用下会很快破碎成碎末, 可有效降低对周围环境的损伤。王玲等[8]进行了纤维增强复合材料为包覆板材料时反应装甲的抗弹性能的实验研究, 发现其具有较好的防弹效果, 且可以消除爆炸反应装甲的二次杀伤效应。

    本文中, 使用口径36mm的精密成型装药, 进行包覆板为氧化铝和碳化硅陶瓷时平板装药的抗弹性能的实验研究, 并与等重厚钢板平板装药进行对比。运用LS-DYNA3D软件对作用过程进行三维数值模拟, 拟为平板装药包覆板新材料的选用提供参考。

    用穿深检验陶瓷包覆板对射流的防护性能。图 1为实验装置示意图, 包括聚能装药、陶瓷平板装药和后效靶, 通过测量残余射流在后效靶上的穿深衡量其防护性能。选用氧化铝陶瓷和碳化硅陶瓷作为平板装药的面背板材料, 并与45钢板进行对比。其中, 氧化铝陶瓷板(w(Al2O3)=99.3%)的厚度为5.1mm, 密度为4.05~4.10g/cm3; 碳化硅陶瓷板的厚度为6.4mm, 密度为3.8~3.9g/cm3; 钢板厚2.5mm, 密度为7.85g/cm3。夹层炸药为B炸药, 密度为1.715g/cm3, 药厚为4mm。后效靶为603装甲钢, 平板装药距后效靶69mm, 倾角均为30°。精密聚能装药口径为36mm, 当炸高为26mm时, 3发平均穿深为100mm。

    图  1  实验装置示意图
    Figure  1.  Experimental setup

    侵彻实验后, 测量后效靶上最大残余穿深衡量平板装药的防护效能, 通过观察后效靶表面的损伤情况衡量平板装药的损伤效应。对不同包覆板材料的平板装药进行单发实验, 结果如表 1图 2所示, 其中abc为弹坑的长、宽和穿深, l为压痕深度。可见, 碳化硅和氧化铝陶瓷平板装药的防护水平优于45钢, 且碳化硅陶瓷比氧化铝陶瓷更优。

    表  1  不同包覆板材料平板装药的防护性能
    Table  1.  Protection performance of ERA with different plate materials
    包覆板材料 a/mm b/mm c/mm l/mm
    SiC 15 10 6 麻点
    45 steel 25 10 10 2~3
    Al2O3 25 11 9 麻点
    下载: 导出CSV 
    | 显示表格
    图  2  实验结果
    Figure  2.  Comparisons of the experimental results

    采用非线性有限元LS-DYNA3D软件对聚能射流和平板装药的作用过程进行数值模拟计算。采用ALE算法, 此种算法中包括拉格朗日与欧拉两种算法。由于在欧拉算法中, 网格固定不动, 材料和能量在网格间运动, 所以在建立模型时不仅包括金属射流部分, 还需加入空气。根据结构的对称性, 考虑计算时间与周期, 建立1/2模型, 模型包括聚能装药、空气、平板装药三部分。其中聚能装药、空气采用欧拉法, 平板装药采用拉格朗日算法。

    聚能药柱使用8701炸药压制, 采用高能燃烧材料模型和JWL状态方程描述。其主要参数为:ρ0=1.685g/cm3, D=8 130m/s, A=625.3GPa, B=23.29GPa, R1=5.25, R2=1.6, ω=0.28。

    紫铜药型罩、45钢聚能装药壳体以及平板装药包覆板, 使用Johnson-Cook模型和Grüneisen状态方程描述:

    σ=(A+Bεp)(1+Clnεp)[1(T)m]
    (1)

    式中:ABCnm为材料常数; εp*为有效塑性应变; ˙εp为有效塑性应变率; Tm为材料熔化温度; Tr为室温。紫铜和45钢材料本构参数见表 2

    表  2  紫铜和45钢的计算参数
    Table  2.  Calculating parameters for the copper and 45 steel
    材料 ρ/(g·cm-3) A/GPa B/GPa n C m Tm/K Tr/K
    Fe 7.85 0.792 0.510 0.26 0.014 1.03 1 793 294
    Cu 8.96 0.090 0.292 0.31 0.025 1.09 1 356 294
    下载: 导出CSV 
    | 显示表格

    空气采用空物质材料模型, 状态方程采用多线性状态方程描述, ρ=1.25mg/cm3

    夹层炸药为B炸药, 采用Lee-Tarver反应模型[9]和JWL状态方程描述, Lee-Tarver反应模型为:

    dλ/d<italic>t</italic>=I(1λ)b(ρ/ρ01a)x+G1(1λ)cλdpy+G2(1λ)eλgpz
    (2)

    式中:λ为反应速率分数; t为时间, p为压力。B炸药的计算参数分别为:I=4.4×1017s-1, b=0.667, a=0, x=20, G1=310GPa, c=0.667, d=0.111, y=1.0, G2=400GPa, e=0.333, g=1.0, z=2.0。

    陶瓷材料采用氧化铝和碳化硅, 使用JH-2模型描述材料在高压下的力学行为。模型主要包括了对材料的强度、压力和损伤的变化关系的描述。通过这3部份关系的耦合, 可以描述完整材料和破碎材料的在载荷作用下的响应。陶瓷材料的本构参数见表 3

    表  3  陶瓷的计算参数[10]
    Table  3.  Calculating parameters for SiC and AL2O3
    材料 G/GPa A B C M N T σc/GPa σHEL/GPa pHEL/GPa K1 K2 K3 β D1 D2 σf, max/GPa
    SiC 183.00 0.96 0.35 0 1.0 0.65 0.37 14.57 13 5.90 204.8 0 0 1.0 0.480 0.480 1.0
    Al2O3 90.16 0.35 0.31 0 0.6 0.60 0.20 2.79 2 2.79 131.0 0 0 1.0 0.005 1.00 1.0
    下载: 导出CSV 
    | 显示表格

    包覆板的等效钢板厚度δequ=ρphp/ρFe, ρphp为包覆板的密度和厚度。根据实验条件可以计算, 氧化铝包覆板的等效钢板厚度为2.50mm, 碳化硅陶瓷的等效钢板厚度为2.51mm, 两者和45钢板的板厚几乎相同。根据Gurney公式, 炸药爆炸后驱动飞板的运动速度为:

    v=2E(2mm/me+0.333)0.5
    (3)

    式中: 2E为Gurney能, mm/me为金属和炸药质量之比。可以看出, 碳化硅陶瓷、氧化铝陶瓷以及45钢板在药厚都为4mm时, 爆炸后的驱动速度应近似相同, 为约1 000m/s。由实验结果可以看出, 氧化铝和碳化硅陶瓷平板装药使射流的垂直侵彻能力分别下降了82%和88%, 而45钢平板装药使射流的侵彻能力下降了80%, 陶瓷板的防护性能不低于等重厚的钢板, 趋势与文献[7]基本相符。

    图 2后效靶表面的损伤结果可以看出:45钢板在高速碰撞后效靶后, 压痕比较完整、深度达到了2~3mm; 而氧化铝和碳化硅陶瓷碰撞后的后效靶表面无完整痕迹, 仅有大量麻点存在, 凹坑深度较浅。说明陶瓷背板在此距离内已完全破碎, 后效损伤与钢板相比已明显降低。

    图 3分别为碳化硅陶瓷、氧化铝陶瓷以及45钢包覆板不同时刻与射流作用的模拟结果。对比图 3可以看出:当t=35μs时, 碳化硅陶瓷板经爆炸驱动后保持完整, 而氧化铝陶瓷板内表面呈现破坏, 原因是由于氧化铝陶瓷抗压强度较低所致。当t=45μs时, 出现了比较有趣的现象, 从陶瓷板周边开始出现破碎和断裂, 原因可能是陶瓷材料为典型的抗压不抗拉材料, 在边界处由于爆轰产物的强烈剪切和稀疏波拉伸作用下, 先从陶瓷板边缘出现破坏。当t=60μs时, 氧化铝陶瓷板已完全破碎, 而碳化硅中间部分仍保持有一定的完整性, 但是表面出现了裂纹。在t=45μs时, 45钢板由于稀疏波的作用在周边出现了层裂, 但整个过程保持完好, 这也正是陶瓷材料在碰撞后效靶表面压痕不明显而钢板压痕比较深的原因。

    图  3  碳化硅、氧化铝和钢反应装甲与射流的作用过程的模拟结果
    Figure  3.  Comparison of the interaction process between shaped charge jets and ERA with SiC, Al2O3 and steel plates

    对比图 3还可以看出, 钢板和陶瓷板与射流的作用机理不一样:当t=45μs时, 氧化铝陶瓷背板后出现了一定长度的逃逸射流, 为整个射流长度的约7%, 随着射流的延伸, 逃逸射流很快断裂成更小的碎片, 而碳化硅陶瓷背板后的逃逸射流不明显。由于氧化铝陶瓷抗压强度较低, 在冲击波作用下内表面破碎成更小的碎末, 对射流的干扰作用减弱, 造成射流头部逃逸。当t=60μs时, 氧化铝陶瓷由于破碎完全, 其与射流尾部的作用减弱, 而碳化硅仍保持一定的完整性, 仍可与充分射流作用, 此为碳化硅陶瓷板抗弹性能优于氧化铝的原因。45钢板与射流作用后造成射流前部出现断裂, 但无明显逃逸射流, 与射流后部作用后使射流局部出现了大的失稳。模拟结果说明, 钢板与射流尾部作用为断续干扰, 而陶瓷板为连续干扰, 其中钢板对射流的干扰周期为约125kHz, 氧化铝为约140kHz, 碳化硅为约160kHz。连续干扰后的射流射流在飞行过程中会很快断裂成破片, 失去侵彻能力, 而断续干扰后的射流在飞行过程中会进一步延伸, 对后效靶仍具有一定的侵彻能力, 造成钢反应装甲的防护性能低于等重厚氧化铝和碳化硅陶瓷反应装甲。

    实验结果显示, 碳化硅和氧化铝陶瓷平板装药具有优良的防护性能, 对于此结构的平板装药, 可使垂直穿深100mm的聚能装药侵彻能力分别下降88%和82%, 优于钢板的防护性能。钢包覆板与后效靶作用后有明显的压痕, 而陶瓷包覆板仅在后效靶表面产生麻点, 压痕不明显。

    模拟结果显示, 陶瓷包覆板先从边缘出现断裂和粉碎现象, 逐渐向中心发展, 当t=60μs时, 已基本破碎。钢板与射流尾部作用为断续干扰, 而陶瓷板为连续干扰。

  • 图  1  实验装置

    Figure  1.  Experimental setup

    图  2  Uc=300 V,h=15 mm,直管中气泡及界面射流发展高速摄影图像

    Figure  2.  Sequential images of the interfacial flow in a straight tube with Uc=300 V and h=15 mm

    图  3  斜压效应诱导射流示意图

    Figure  3.  Schematic of jet formation by baroclinic effect

    4(a)  Uc=300 V,h=5 mm, 直管中气泡及界面射流发展高速摄影图像

    4(a).  Sequential images of the interfacial flow in a straight tube with Uc=300 V and h=5 mm

    4(b)  Uc=300 V, h=25 mm时, 直管中气泡及界面射流发展高速摄影图像

    4(b).  Sequential images of the interfacial flow in a straight tube with Uc=300 V and h=25 mm

    图  5  气泡最大长度和表面射流速度随爆炸深度的变化

    Figure  5.  Maximum bubble length and surface jet velocity varied with explosion depth

    图  6  直管中不同深度爆炸波与自由界面作用示意图

    Figure  6.  Schematic of interaction patterns of a free surface and an explosion at varied depths in a straight tube

    图  7  h=20 mm, Uc=150, 350 V时,直管中气泡及界面射流发展高速摄影图像

    Figure  7.  Sequential images of the interfacial flow in a straight tube with h=20 mm and Uc=150, 350 V

    图  8  气泡最大长度和表面射流速度随充电电压(能量)的变化

    Figure  8.  Maximum bubble length and surface jet velocity varied with charging voltage (energy)

    图  9  理论模型设置示意图

    Figure  9.  Schematic of theoretical model

    图  10  理论计算得到的气泡长度随时间的变化

    Figure  10.  Evolution history of bubble length by theoretical calculation

    表  1  气柱初始条件设置

    Table  1.   Initial conditions for gas column

    初始条件L0/mmpb0/MPa
    10.001367.00
    20.01038.20
    30.1004.23
    下载: 导出CSV
  • [1] Hirano T, Komatsu M, Takahashi A, et al. Enhancement of fibrinolytics with a laser-induced liquid jet[J]. Lasers in Surgery and Medicine, 2001, 29(4): 360-368. doi: 10.1002/lsm.1129
    [2] Nakagawa A, Hirano T, Komatsu M, et al. Holmium: YAG laser-induced liquid jet knife: Possible novel method for dissection[J]. Lasers in Surgery and Medicine, 2002, 31(2): 129-135. doi: 10.1002/lsm.10055
    [3] Blake J R, Taib B B, Doherty G. Transient cavities near boundaries: Part 1: Rigid boundary[J]. Journal of Fluid Mechanics, 1986, 170: 479-497. doi: 10.1017/S0022112086000988
    [4] Blake J R, Taib B B, Doherty G. Transient cavities near boundaries: Part 2: Free surface[J]. Journal of Fluid Mechanics, 1987, 181: 197-212. doi: 10.1017/S0022112087002052
    [5] Blake J R, Gibson D C. Cavitation bubbles near boundaries[J]. Annual Review of Fluid Mechanics, 1987, 19: 99-123. doi: 10.1146/annurev.fl.19.010187.000531
    [6] Klaseboer E, Hung K C, Wang C, et al. Experimental and numerical investigation of the dynamics of an underwater explosion bubble near a resilient/rigid structure[J]. Journal of Fluid Mechanics, 2005, 537: 387-413. doi: 10.1017/S0022112005005306
    [7] Klaseboer E, Khoo B C. An oscillating bubble near an elastic material[J]. Journal of Applied Physics, 2004, 96(10): 5808. doi: 10.1063/1.1803925
    [8] Klaseboer E, Turangan C K, Khoo B C. Dynamic behavior of a bubble near an elastic infinite interface[J]. International Journal of Multiphase Flow, 2006, 32(9): 1110-1122. doi: 10.1016/j.ijmultiphaseflow.2006.05.005
    [9] Ong G P, Khoo B C, Turangan C, et al. Behavior of oscillating bubbles near elastic membranes: An experimental and numerical study[J]. Modern Physics Letters: B, 2005, 19(28/29): 1579-1582. doi: 10.1142/S021798490500995X
    [10] Turangan C K, Ong G P, Klaseboeret E, et al. Experimental and numerical study of transient bubble-elastic membrane interaction[J]. Journal of Applied Physics, 2006, 100(5): 054910. doi: 10.1063/1.2338125
    [11] Miao H, Gracewski S M. Coupled FEM and BEM code for simulating acoustically excited bubbles near deformable structures[J]. Computational Mechanics, 2008, 42(1): 95-106. doi: 10.1007/s00466-007-0238-y
    [12] Xie W F, Young Y L, Liu T G, et al. Multiphase modeling of dynamic fluid-structure interaction during close-in explosion[J]. International Journal for Numerical Methods in Engineering, 2008, 74(6): 1019-1043. doi: 10.1002/nme.2207
    [13] Dadvand A, Khoo B C, Shervani-Tabar M T. A collapsing bubble induced micro injector: An experimental study[J]. Experiments in Fluids, 2009, 46(3): 419-434. doi: 10.1007/s00348-008-0568-3
    [14] 逄春京.水下爆炸气泡动过程的数值模拟研究[D].长沙: 国防科技大学, 2008: 28-59.
    [15] Zhang A M, Cui P, Wang Y. Experiments on bubble dynamics between a free surface and a rigid wall[J]. Experiments in Fluids, 2013, 54: 1602. doi: 10.1007/s00348-013-1602-7
  • 期刊类型引用(8)

    1. 史安顺,洪晓文,秦超,樊雪飞,夏开朗,孙素杰. 爆炸载荷下拼接陶瓷板膨胀破碎行为及飞散特性研究. 火工品. 2024(01): 65-70 . 百度学术
    2. 洪晓文,史安顺,秦超,樊雪飞,夏开朗,尹飞. 爆轰驱动陶瓷面背板碎裂行为研究. 兵器装备工程学报. 2024(02): 188-196 . 百度学术
    3. 吴立林,钱晨晨,陆静. 陶瓷防护结构对引信抗破片冲击性能影响研究. 兵器装备工程学报. 2024(06): 159-165 . 百度学术
    4. 张华申,李如江,梁宇. 橡胶叠层装药反应装甲干扰聚能射流作用研究. 火工品. 2023(03): 27-31 . 百度学术
    5. 梁宇,李如江,张华申,王云鹏. 射流斜侵彻薄层装药结构数值模拟和试验研究. 火工品. 2023(04): 43-49 . 百度学术
    6. 周杰,王凤英,原诗瑶,吴鹏. 楔形装药对射流干扰的数值模拟. 高压物理学报. 2018(02): 135-142 . 百度学术
    7. 聂鹏松,刘天生,阮光光,吴鹏,刘金彪. 铝合金平板装药对射流干扰的数值模拟. 火炮发射与控制学报. 2018(01): 43-47 . 百度学术
    8. 宋思维,刘天生,王凤英,晏蜜. 碳化硅及氧化铝陶瓷反应装甲附带损伤研究. 科学技术与工程. 2016(13): 151-155 . 百度学术

    其他类型引用(3)

  • 加载中
图(11) / 表(1)
计量
  • 文章访问数:  2951
  • HTML全文浏览量:  410
  • PDF下载量:  256
  • 被引次数: 11
出版历程
  • 收稿日期:  2014-02-19
  • 修回日期:  2015-01-20
  • 刊出日期:  2015-10-10

目录

/

返回文章
返回