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

钢箱内部爆炸破坏的SPH数值模拟

强洪夫 孙新亚 王广 黄拳章

引用本文:
Citation:

钢箱内部爆炸破坏的SPH数值模拟

    作者简介: 强洪夫(1963- ),男,博士,教授,博士生导师,qiang@263.net;
    通讯作者: 孙新亚, 1430167246@qq.com
  • 中图分类号: O385

Numerical simulation on steel box damage under internal explosion by smoothed particle hydrodynamics

    Corresponding author: SUN Xinya, 1430167246@qq.com ;
  • CLC number: O385

  • 摘要: 随着恐怖袭击的不断演化,船舶、桥梁等以钢箱为主要支撑的战略性结构逐渐成为恐怖袭击和敌方军事打击的重要目标。本文中采用光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)方法对钢箱内部爆炸过程进行了数值模拟,对比实验,分析了钢箱内部爆炸时钢箱表面的变形过程,得到了钢箱表面挠度值的变化趋势、不同时刻钢箱表面压强和von Mises应力的分布情况、钢箱表面中心点处速度和压强的变化趋势,验证了SPH方法在模拟钢箱内部爆炸问题上的有效性。通过进一步数值模拟,探讨了钢箱内部不同位置爆炸时钢箱的破坏形式和损伤程度,结果表明:炸药在钢箱内部角隅处爆炸时,钢箱的损伤程度最严重;炸药在钢箱内部正中心爆炸时,钢箱的损伤程度最轻。
  • 图 1  钢箱和TNT炸药试件模型图

    Figure 1.  Model diagrams for the steel box and TNT dynamite specimen

    图 2  钢箱内部爆炸算例模型

    Figure 2.  The case model for the explosion inside the steel box

    图 3  钢箱内部爆炸过程中不同时刻箱体内壁的变形情况

    Figure 3.  Its inner wall deformation at different times for the steel box subjected to inner explosion loading

    图 4  钢箱壁面破坏的数值模拟和实验对比

    Figure 4.  Comparison of steel box wall damage between numerical simulation and experiment

    图 5  钢箱壁面挠度变化趋势

    Figure 5.  Deflection in the center of the steel box wall varying with time

    图 6  钢箱内部爆炸过程中,壁面von Mises应力分布

    Figure 6.  Von Mises stress inside the steel wall during the explosion in the steel box

    图 7  钢箱内部爆炸过程中,壁面的压力分布

    Figure 7.  Pressure distribution of the wall during the explosion in the steel box

    图 8  钢箱壁面中心点处的速度变化趋势

    Figure 8.  Velocity varying with time at the wall center of the steel box

    图 9  钢箱壁面中心点处的压强变化趋势

    Figure 9.  Pressure varying with time at the wall center of the steel box

    图 10  钢箱中炸药位置示意图(a=60 mm)

    Figure 10.  The explosive location diagram in the steel box (a=60 mm)

    图 11  3种工况下,钢箱表面局部破坏示意图

    Figure 11.  Surface partial destruction diagrams of the steel box under three different conditions

  • [1] 吴保桦, 张尚根, 康正炎. 桁架式钢筋混凝土叠合板式墙抗爆性能实验 [J]. 爆炸与冲击, 2017, 37(1): 92–98. DOI: 10.11883/1001-1455(2017)01-0092-07.
    WU Baohua, ZHANG Shanggen, KANG Zhengyan. Anti-blast properties of RC superimposed slab shear wall [J]. Explosion and Shock Waves, 2017, 37(1): 92–98. DOI: 10.11883/1001-1455(2017)01-0092-07.
    [2] 姚术健, 张舵, 郑监, 等. 内部爆炸作用下钢箱结构变形规律性实验 [J]. 爆炸与冲击, 2017, 37(5): 964–968. DOI: 10.11883/1001-1455(2017)05-0964-05.
    YAO Shujian, ZHANG Duo, ZHENG Jian, et al. Experimental study of deformation of steel box subjected to internal blast loading [J]. Explosion and Shock Waves, 2017, 37(5): 964–968. DOI: 10.11883/1001-1455(2017)05-0964-05.
    [3] 耿少波, 刘亚玲, 薛建英. 钢箱梁缩尺模型爆炸冲击波作用下破坏实验研究 [J]. 工程力学, 2017, 34(S): 84–88. DOI: 1000-4750(2017)Suppl-0084-05.
    GENG Shaobo, LIU Yaling, XUE Jianying. Experimental studies on steel box girder scale model under the blast load [J]. Engineering Mechaincs, 2017, 34(S): 84–88. DOI: 1000-4750(2017)Suppl-0084-05.
    [4] 蒋志刚, 朱新明, 严波, 等. 钢箱梁爆炸冲击局部破坏的数值模拟 [J]. 振动与冲击, 2013, 32(13): 189–164. DOI: 10.13465/j.cnki.jvs.2013.13.008.
    JIANG Zhigang, ZHU Xinming, YAN Bo, et al. Numerical simulation for local failure of a steel box girder under blast loading [J]. Journal of Vibration and Shock, 2013, 32(13): 189–164. DOI: 10.13465/j.cnki.jvs.2013.13.008.
    [5] TANG E K C, HAO H. Numerical simulation of a cable-stayed bridge response to blast loads: part I: model development and response calculations [J]. Engineering Structures, 2010, 32(10): 3180–3192. doi: 10.1016/j.engstruct.2010.06.007
    [6] 姚术健, 蒋志刚, 卢芳云, 等. 汽车炸弹钢箱梁内部爆炸局部破坏效应分析 [J]. 振动与冲击, 2015, 34(7): 222–227. DOI: 10.13465/j.cnki.jvs.2015.07.035.
    YAO Shujian, JIANG Zhigang, LU Fangyun. Analysis on local damage of steel box girder under internal blast loading of vehicle bomb [J]. Journal of Vibration and Shock, 2015, 34(7): 222–227. DOI: 10.13465/j.cnki.jvs.2015.07.035.
    [7] MONAGHAN J J. Smoothed particle hydrodynamics [J]. Annual Review of Astronomical and Astrophysics, 1992(30): 543–574.
    [8] 强洪夫, 范树佳, 陈福振, 等. 基于拟流体模型的SPH新方法及其在弹丸超高速碰撞薄板中的应用 [J]. 爆炸与冲击, 2017, 37(6): 990–1000. DOI: 10.11883/1001-1455(2017)06-0990-11.
    QIANG Hongfu, FAN Shujia, CHEN Fuzhen, et al. A new smoothed particle hydrodynamics method based on the pseudo-fluid model and its application in hypervelocity impact of a projectile on a thin plate [J]. Explosion and Shock Waves, 2017, 37(6): 990–1000. DOI: 10.11883/1001-1455(2017)06-0990-11.
    [9] QIANG Hongfu, WANG Kunpeng, GAO Weiran. Numerical simulation of shaped charge jet using multiphase SPH method [J]. Transactions of Tianjin University, 2008, 14(1): 495–499.
    [10] LIU M B, LIU G R, LAM K Y. Meshfree particle simulation of the detonation process for high explosives in shaped charge unlined cavity configurations [J]. Shock Waves, 2003, 12(6): 509–520. doi: 10.1007/s00193-003-0185-2
    [11] JOHNSON G R, COOK W H. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures [C] // Proceedings of the Seventh International Symposium on Ballistics. The Hague, 1983.
    [12] 强洪夫. 光滑粒子流体动力学方法及应用 [M]. 北京: 科学出版社, 2017: 246−287.
    [13] BØRVIK T, LANGSETH M, HOPPERSTAD O S, et al. Ballistic penetration of steel plates [J]. International Journal of Impact Engineering, 1999, 22(9-10): 855–886. doi: 10.1016/S0734-743X(99)00011-1
    [14] NILSON A H. Bond stress-slip relationships in reinforced concrete: 345 [R]. New York: Department of Structure Engineering, Cornell University, 1971.
    [15] BONORCHIS D, NURICK G N. The analysis and simulation of welded stiffener plates subjected to localized blast loading [J]. International Journal of Impact Engineering, 2010, 37(3): 260–273. doi: 10.1016/j.ijimpeng.2009.08.004
    [16] 朱新明. 钢箱梁爆炸冲击局部破坏数值模拟研究 [D]. 长沙: 国防科技大学, 2011: 23−27.
  • [1] 崔云霄12胡永乐王春明胡昊陈鹏万 . 内部爆炸作用下多层钢筒的动态响应. 爆炸与冲击, 2015, 35(6): 820-824. doi: 10.11883/1001-1455(2015)06-0820-05
    [2] 郭伟国李玉龙黄福增 . 不同应变率下泡沫铝的形变和力学性能. 爆炸与冲击, 2008, 28(4): 289-292. doi: 10.11883/1001-1455(2008)04-0289-04
    [3] 宫翔飞刘文韬张树道杨基明 . 水下爆炸近场峰值压力的数值模拟. 爆炸与冲击, 2019, 39(4): 041409-1-041409-8. doi: 10.11883/bzycj-2017-0262
    [4] 张铱鈖赵隆茂 . 非均匀泡沫金属材料在冲击载荷下的变形模拟. 爆炸与冲击, 2006, 26(1): 33-38. doi: 10.11883/1001-1455(2006)01-0033-06
    [5] 徐志宏汤文辉罗永 . 光滑粒子模拟方法在超高速碰撞现象中的应用. 爆炸与冲击, 2006, 26(1): 53-58. doi: 10.11883/1001-1455(2006)01-0053-06
    [6] 张春波宋振飞谷卓伟卢纪赵士操 . 内爆压缩多层密绕螺线管的数值模拟. 爆炸与冲击, 2018, 38(5): 999-1005. doi: 10.11883/bzycj-2016-0052
    [7] 肖毅华胡德安韩旭杨刚 . 一种自适应轴对称FEM-SPH 耦合算法及其在高速冲击模拟中的应用. 爆炸与冲击, 2012, 32(4): 384-392. doi: 10.11883/1001-1455(2012)04-0384-09
    [8] 秦学军张德志杨军史国凯刘峻岭王惠熊琛 . 内部爆炸作用下钢筒变形过程的电探针测量技术. 爆炸与冲击, 2014, 34(1): 115-119. doi: 10.11883/1001-1455(2014)01-0115-05
    [9] 李艳超毕明树张大为高伟 . 热扩散不稳定和流体动力学不稳定对氢气/空气爆炸胞状火焰的影响. 爆炸与冲击, 2018, 38(5): 1064-1070. doi: 10.11883/bzycj-2017-0069
    [10] L.L.AltgilbersS.T.WuY.C.Xiao . 平面磁通量压缩流体动力学的数值模拟. 爆炸与冲击, 1991, 11(4): 297-308.
    [11] 柏劲松陈森华李平张展冀 . 多介质可压缩流体动力学界面捕捉方法. 爆炸与冲击, 2004, 24(1): 37-43.
    [12] 汤文辉张若棋赵国民 . 流体动力学数值计算中两种自由面格式的比较. 爆炸与冲击, 1995, 15(2): 155-159.
    [13] 邹广平孙杭其唱忠良熊海林 . 聚氨酯/钢夹芯结构爆炸载荷下动力学响应的数值模拟. 爆炸与冲击, 2015, 35(6): 907-912. doi: 10.11883/1001-1455(2015)06-0907-06
    [14] 周南程漱玉陈雨生丁升 . 弹-塑性流体动力学数值模拟计算的一种新自由面格式. 爆炸与冲击, 1997, 17(1): 10-14.
    [15] 梁霄王瑞利 . 混合不确定度量化方法及其在计算流体动力学迎风格式中的应用. 爆炸与冲击, 2016, 36(4): 509-515. doi: 10.11883/1001-1455(2016)04-0509-07
    [16] 崔溦宋慧芳张社荣 . 土中爆炸作用下箱涵动力响应的SPH-FE耦合分析. 爆炸与冲击, 2012, 32(5): 551-556. doi: 10.11883/1001-1455(2012)05-0551-06
    [17] 陆景德张万甲张玉松 . 2~#纯铁和不锈钢的变形硬化效应. 爆炸与冲击, 1984, 4(2): 75-79.
    [18] 李继承陈小伟陈刚 . 921A钢纯剪切帽状试件在SHPB实验中的动态变形. 爆炸与冲击, 2010, 30(3): 239-246. doi: 10.11883/1001-1455(2010)03-0239-08
    [19] 杨诗婷邢永明1赵燕茹郝贠洪李继军姜爱峰 . 单个弹丸撞击316L不锈钢引起的变形场. 爆炸与冲击, 2017, 37(1): 126-133. doi: 10.11883/1001-1455(2017)01-0126-08
    [20] 王吉王肖钧卞梁 . 光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用. 爆炸与冲击, 2007, 27(6): 522-528. doi: 10.11883/1001-1455(2007)06-0522-07
  • 加载中
图(11)
计量
  • 文章访问数:  90
  • HTML全文浏览量:  437
  • PDF下载量:  27
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-12-06
  • 录用日期:  2018-01-26
  • 刊出日期:  2019-05-01

钢箱内部爆炸破坏的SPH数值模拟

    作者简介:强洪夫(1963- ),男,博士,教授,博士生导师,qiang@263.net
    通讯作者: 孙新亚, 1430167246@qq.com
  • 火箭军工程大学,陕西 西安 710025

摘要: 随着恐怖袭击的不断演化,船舶、桥梁等以钢箱为主要支撑的战略性结构逐渐成为恐怖袭击和敌方军事打击的重要目标。本文中采用光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)方法对钢箱内部爆炸过程进行了数值模拟,对比实验,分析了钢箱内部爆炸时钢箱表面的变形过程,得到了钢箱表面挠度值的变化趋势、不同时刻钢箱表面压强和von Mises应力的分布情况、钢箱表面中心点处速度和压强的变化趋势,验证了SPH方法在模拟钢箱内部爆炸问题上的有效性。通过进一步数值模拟,探讨了钢箱内部不同位置爆炸时钢箱的破坏形式和损伤程度,结果表明:炸药在钢箱内部角隅处爆炸时,钢箱的损伤程度最严重;炸药在钢箱内部正中心爆炸时,钢箱的损伤程度最轻。

English Abstract

  • 桥梁是公路、铁路等交通网的咽喉。为了减小恒载,对斜拉桥、悬索桥等大跨度桥广泛采用钢箱梁作其结构支撑的主梁。随着恐怖袭击的不断演化,战略性桥梁成为恐怖主义爆炸破坏的重点目标,因此,有必要对钢箱结构内部爆炸过程中钢箱的变形规律进行研究[1-4]。Tang等[5]运用LS-DYNA软件研究了汽车炸弹爆炸冲击作用下桥塔、桥墩、桥面结构的局部破坏模式;姚术健等[2,6]通过实验,研究了钢箱内部爆炸过程中内壁失效变形的模式, 并利用ANSYS软件分析了箱体的损伤特性。钢箱内部爆炸相比外部爆炸而言,爆炸冲击波在约束空间内多次反射、叠加、汇聚,该物理过程相当复杂,目前的研究大多局限于爆炸过程中钢箱变形的宏观现象。光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)作为一种无网格的粒子方法,因其在计算过程中可避免网格重分和算法耦合,被广泛用于研究爆炸冲击问题。自Monaghan等[7]修正的人工黏度解决了SPH计算结果的非物理震荡问题后,强洪夫等[8]、Qiang等[9]和Liu等[10]利用该方法对爆炸冲击波冲击钢板的过程成功地进行了数值模拟。本文中,拟利用自编的SPH程序对钢箱内部爆炸过程进行数值模拟,并与实验对比,验证SPH算法在模拟钢箱内部爆炸问题上的有效性;分析爆炸过程中钢箱表面压力、应力波的传播规律以及特定点处速度变化趋势;得到钢箱表面的von Mises应力分布;进一步对钢箱内不同位置爆炸的过程进行数值模拟,得到不同工况下钢箱的失效变形模式,以期通过研究为钢箱类结构的爆炸毁伤提供一种评估方法。

    • 钢箱的内部爆炸是绝热过程,外力相对于内部驱动力可忽略不计,即可不考虑重力的影响,在流体控制方程中忽略黏性项。因此,其控制方程组表述为:

      $ \left\{ \begin{aligned} &{\frac{\rm{d} \rho}{{\rm{d}} t}=-\rho \nabla \cdot {v}} \\ &{\frac{{\rm{d}} {{v}}}{{\rm{d}} t}={-\frac{\nabla p}{\rho}}}\\ &{\frac{{\rm{d}} e}{{\rm{d}} t}=-\frac{p}{\rho} \nabla \cdot {v}}\\ & {\frac{{\rm{d}} {x}}{{\rm{d}} t}={v}} \end{aligned} \right. $

      式中:ρvpex分别为密度、速度、压强、能量、位移。

    • 数值模拟过程中,TNT炸药的爆轰气体压力通过标准的JWL状态方程[11]获得:

      $ p(e,v) = A\left( {1 - \frac{w}{{{R_1}v}}} \right)\exp ( - {R_1}v) + B\left( {1 - \frac{w}{{{R_2}v}}} \right)\exp ( - {R_2}v) + \frac{{w{\rho _0}e}}{v} $

      式中:$v = \displaystyle\frac{1}{\rho }$$\rho $为爆轰气体的密度;e为爆轰气体的比内能;$A\text{、}B\text{、}{R_1}\text{、}{R_2}\text{、}w$为实验拟合得到的系数[11],且$A = 371.2\;{\rm{ GPa}}$, $B = 3.231\;{\rm{ GPa}}$, ${R_1} = 4.15$, ${R_2} = 0.95$, $w = 0.3$${\rho _0} = 1\; 500\;{\rm{ kg/}}{{\rm{m}}^{\rm{3}}}$为TNT炸药初始密度。

      对于钢箱材料,选用Q235钢。爆炸过程中钢箱的各向同性压力p利用Mie-Grüneisen状态方程[12]求得:

      $ p(\rho ,e) = \left( {1 - \frac{1}{2}\gamma \eta } \right){p_{\rm{H}}}(\rho ) + \gamma \rho e $

      $ {p_{\rm{H}}} = \left\{ {\begin{array}{*{20}{l}} {{a_0}\eta + {b_0}{\eta ^2} + {c_0}{\eta ^3}}&{\eta \text{≥} 0}\\ {{a_0}\eta }&{\eta \text{<} 0} \end{array}} \right. $

      式中:$\gamma = 1.99$为Grüneisen参量;$\eta = \displaystyle\frac{\rho }{{{\rho _0}}} - 1$${\rho _0} = 7\;800\;{\rm{kg}}\rm{/}{{\rm{m}}^{\rm{3}}}$为Q235钢的初始密度;${a_0} = {\rho _0}c_{\rm{s}}^2$${c_{\rm{s}}} = 4\;569\;{\rm{m}}\rm{/}{\rm{s}}$为冲击速度与粒子速度的线性相关常量;${b_0} = {a_0}[1 + 2({S_{\rm{s}}} - 1)]$${S_{\rm{s}}} = 1.49$为斜率[12]${c_0} = {a_0}[2({S_{\rm{s}}} - 1) + $$ 3{({S_{\rm{s}}} - 1)^2}]$

    • 钢箱材料模型采用Børvik等[13]、Nilson[14]提出的修正Johnson-Cook本构模型,该模型通过流动应力和失效应变描述金属材料的大变形、高应变率等动态力学特性。其流动应力的表达式为:

      $ {\sigma _{\rm{eq}}} = (1 - D)(A + B{r^n})(1 + C\ln {\dot r^ * })(1 - {T^{ * m}}) $

      式中:$r = (1 - D)g$为累积损伤塑性应变,${\dot r^ * } = (1 - D){\dot g^ * }$${\dot g^ * } = \displaystyle\frac{g}{{{{\dot g}_0}}}$g为累积塑性应变,${\dot g_0}$为自定义参考应变率,${T^ * } = (T - {T_0})/({T_{\rm{m}}} - {T_0})$T0为室温,Tm为材料熔点,$A{\text{、}}B{\text{、}}C{\text{、}}n{\text{、}}m$为材料常数。${D = }\sum {{\rm{\Delta }}{g}} {/}{{g}_{\rm{f}}}$$D = 0$表示材料没有损伤,$D = 1$表示材料完全失效, ${g_{\rm{f}}} = [{D_1} + {D_2}\exp ({D_3}{\sigma ^*})](1 + {D_4}\ln {\dot p^*})(1 + {D_5}{T^*})$为断裂塑性应变,${\dot p^*} = \displaystyle\frac{p}{{p_0^*}}$$p_0^*$为初始压力,${D_1} \sim {D_5}$为材料常数, ${\sigma ^ * } = {\sigma _{\rm{mean}}}/{\sigma _{\rm{eq}}}$为应力三轴度,${\sigma _{\rm{mean}}} = ({\sigma _x} + {\sigma _y} + {\sigma _{\textit z}})/3$为应力平均值。

    • SPH算法是模拟流体流动的一种拉格朗日型粒子方法,通过使用一系列任意分布的粒子来求解具有各种边界条件的积分方程或偏微分方程。SPH方法通常通过核函数插值实现场变量或场变量梯度的插值,通过粒子近似实现对核函数估计积分表达式的粒子离散[9]

      本文中利用SPH方法模拟钢箱内部爆炸过程属于该方法在具有大变形和高度非均匀的动力学极端情形下的应用。为了解决爆炸域内求解结果的非物理振荡,且将冲击面内的动能以热能的形式耗散,引入人工黏性[7],同时为了更好地处理爆炸过程中钢箱大变形、大扭曲等引起的密度和光滑长度变化剧烈的问题,采用Qiang等[9]提出的完全变光滑长度SPH方法对钢箱内部爆炸过程进行数值模拟,结合后的方程组如下:

      $ \left\{ {\begin{aligned} &{\frac{{{\rm{d}}{\rho _i}}}{{{\rm{d}}t}} = \sum\limits_{j = 1}^{{N}} {{m_j}} \left[ {{{{v}}_{ij}} \cdot {\nabla _i}{W_{ij}} + \frac{1}{2}\left( {\frac{{{\rm{d}}{h_i}}}{{{\rm{d}}t}} + \frac{{{\rm{d}}{h_j}}}{{{\rm{d}}t}}} \right)\frac{{\partial {W_{ij}}}}{{\partial h}}} \right]}\\ &{\frac{{{\rm{d}}v_i^\alpha }}{{{\rm{d}}t}} = \sum\limits_{j = 1}^{{N}} {{m_j}} \left( {{f_i}\frac{{p_i^{\alpha \beta }}}{{\rho _i^2}} + {f_j}\frac{{p_j^{\alpha \beta }}}{{\rho _j^2}} + {\displaystyle\varPi _{ij}}} \right)\frac{{\partial {W_{ij}}}}{{\partial x_i^\beta }}}\\ &{\frac{{{\rm{d}}{e_i}}}{{{\rm{d}}t}} = \left( {\frac{{{p_i} + {p_j}}}{{{\rho _i} + {\rho _j}}} + {\displaystyle\varPi _{ij}}} \right)\frac{{{\rm{d}}{\rho _i}}}{{{\rm{d}}t}}}\\ &{\frac{{{\rm{d}}{{{x}}_i}}}{{{\rm{d}}t}} = {{{v}}_i}} \end{aligned}} \right. $

      $ {f_i} = {\left( {1 + \frac{{{h_i}}}{{d{\rho _i}}}\sum\limits_j {{m_j}\frac{{\partial {W_{ij}}}}{{\partial {h_i}}}} } \right)^{ - 1}} $

      式中:${f_i}$为修正系数;${W_{ij}} = W({x_i} - {x_j},h)$为核函数,一般选用三次样条函数;Πij为人工黏度;N为求解区域内粒子总数;αβ表示空间坐标轴方向; h为光滑长度,在连续性方程中将其对时间求导,得:

      $ \frac{{\rm{d}} h}{{\rm{d}} t}=-\frac{1}{ d} \frac{h_{i}}{\rho_{i}} \frac{{\rm d} \rho_{i}}{{\rm d} t} $

      式中:d为空间维数。

    • 为了证明SPH方法在数值模拟钢箱结构内部爆炸过程时的有效性,本文中设定了与姚术健等[2]的实验一致的工况条件,如图1(a)所示,钢箱模型为立方体形状,边长为600 mm,且采用了120 mm的延伸结构,Q235钢板厚度为4 mm;TNT炸药为立方体形状,边长为40 mm。具体粒子配置如图1(b)所示,粒子间距为4 mm,炸药粒子数为1 331,钢箱粒子数为265 600。

      图  1  钢箱和TNT炸药试件模型图

      Figure 1.  Model diagrams for the steel box and TNT dynamite specimen

      为保证箱体各壁面受力均衡,TNT炸药位于钢箱内部正中心,且各边与箱体对应边平行,如图2所示。计算中,TNT炸药密度为1 500 kg/m3, 采用点起爆方式,起爆点设置在炸药上表面中心点(−73.9,5.6,50.7),爆速为6 930 m/s,初始内能为4 290 kJ/kg,状态方程选用JWL状态方程;钢箱密度为7 800 kg/m3, 使用Johnson-Cook(J-C)本构模型,具体材料参数A=369 MPa, B=684 MPa, n=0.73, C=0.083, m=1.7, Tm=775 K, ${\dot g_0}=0.1$ s−1, CV=875 J/(kg·K), D1=0.13, D2=0.13, D3=−1.5, D4=0.011, D5=0。在组成钢箱的各板之间施加焊接锚固作用,由文献[15]可知焊缝金属强度不应低于母材。为简化计算,忽略焊缝附近区域母板材料性能的变化,将焊缝与母材视为同一结构,即认为钢箱是一体的。

      图  2  钢箱内部爆炸算例模型

      Figure 2.  The case model for the explosion inside the steel box

    • 数值模拟中,核函数选用三次样条函数,光滑长度取1.5倍粒子间距,时间积分采用蛙跳格式,时间步长取0.1 μs。钢箱内部爆炸的数值模拟中,钢箱表面的变形情况如图3所示。起始时刻,钢箱表面无变形,1.97 ms时中心处产生明显的向外鼓起,并向四周传播变形;3.34 ms时变形首次到达表面边沿处,且表面角隅处发生变形;12.03 ms时,表面变形达到最大,并保持不变。

      图  3  钢箱内部爆炸过程中不同时刻箱体内壁的变形情况

      Figure 3.  Its inner wall deformation at different times for the steel box subjected to inner explosion loading

      图4给出了钢箱内部爆炸实验和数值模拟的对比,图中A处钢板发生卷边破坏,B处钢板向外鼓起,C处钢板发生扭曲变形,通过对比发现,数值模拟和实验所得钢箱壁面失效变形模式吻合较好。图5给出了钢箱内部爆炸过程中箱体壁面挠度的变化情况。从图5可以看出,85 μs时钢箱壁面开始变形,挠度值随时间逐渐增大,12.03 ms后基本保持不变。数值模拟得到的钢箱壁面最终挠度值为31.80 mm,而实验值为32.93 mm,实验和数值模拟的相对误差为3.5%。产生这一误差的主要原因可能是:(1)炸药理论设计的质量为96 g,由于密度不均匀等原因实验中炸药的真实质量却为98.443 g;(2)本文数值模拟过程中Q235钢材料密度是均匀的,没有考虑制作中的工艺误差。综上所述,可得SPH方法能够数值模拟钢箱内部爆炸过程中钢箱壁面的变形情况。

      图  4  钢箱壁面破坏的数值模拟和实验对比

      Figure 4.  Comparison of steel box wall damage between numerical simulation and experiment

      图  5  钢箱壁面挠度变化趋势

      Figure 5.  Deflection in the center of the steel box wall varying with time

    • 装药在完全封闭的钢箱内部爆炸时,爆炸冲击波在钢箱四周壁面上发生多次反射,冲击波强度逐渐降低。图6给出了钢箱内部爆炸时,钢箱表面的von Mises应力分布;图7给出了钢箱内部爆炸过程中,钢箱表面压力分布情况。从图6图7可得:钢箱内部炸药爆炸后,很快产生爆炸冲击波并迅速(85.0 μs)到达钢箱表面,爆心投影点附近出现压强和应力最大值;随后冲击波沿着钢箱的内表面向板与板的连接处(角隅位置)运动,同时爆心投影点处超压峰值逐渐降低,436.1 μs时各板向角隅运动的冲击波在角隅处汇聚,形成较强的冲击波压力,并出现应力集中和应力增强现象;随后各板反射的冲击波与角隅处汇聚加强了的冲击波向箱室中心汇聚,使箱室中心形成较高超压,箱室中心超压又会向四周传播,此过程与起爆后传播过程类似,钢箱表面承受多次类似的冲击载荷作用后,钢箱表面压强和von Mises应力逐渐消失,且箱内维持较长时间的余压。

      图  6  钢箱内部爆炸过程中,壁面von Mises应力分布

      Figure 6.  Von Mises stress inside the steel wall during the explosion in the steel box

      图  7  钢箱内部爆炸过程中,壁面的压力分布

      Figure 7.  Pressure distribution of the wall during the explosion in the steel box

    • 图8给出了钢箱内部爆炸过程中,钢箱表面中心点处速度的变化趋势;图9给出了钢箱内部爆炸过程中,钢箱表面中心点处压强的变化趋势。从图8可以看出:101.2 μs时(如点A所示),表面中心点的速度跳跃式增加至22.3 m/s,其后速度值又迅速减小,最终减小至0,且速度值在B区域内出现上下微小浮动。从图9可以看出:85.0 μs时,钢箱表面中心点处的压强迅速增大,达到峰值487 MPa; 131.6 μs时表面压强达到第二次峰值394 MPa;随后压强值上下扰动并逐渐减弱,出现扰动的主要原因是爆炸冲击波在钢箱表面上的来回反射;2 000.0 μs时压强值逐渐减弱直至消失。由文献[16]可得:压力和速度时程曲线形状与实验及经典冲击波理论相吻合,即初次峰值较大,然后发生多次反射,速度和压强逐渐减弱。

      图  8  钢箱壁面中心点处的速度变化趋势

      Figure 8.  Velocity varying with time at the wall center of the steel box

      图  9  钢箱壁面中心点处的压强变化趋势

      Figure 9.  Pressure varying with time at the wall center of the steel box

    • 为了进一步探讨同一爆炸当量下,钢箱内部不同位置爆炸时,钢箱结构的破坏程度。本文设定了以下3种计算工况(如图10所示):工况1,炸药位于钢箱底部正中心;工况2,炸药在底部中心正上方60 mm处;工况3,炸药位于底板与侧板之间的角隅处,距离各板60 mm。

      图  10  钢箱中炸药位置示意图(a=60 mm)

      Figure 10.  The explosive location diagram in the steel box (a=60 mm)

      图11给出了3种工况下钢箱表面的失效变形模式。图11(a)为工况1钢箱的变形情况,钢箱底板局部特征破坏明显,产生圆形的冲切破坏,断裂口处相对比较平滑,且底板产生裂纹(A);同时底板与侧板的角隅处产生了长条枣核形状的裂口(B)。图11(b)为工况2钢箱的变形情况,钢箱底板剧烈破坏,产生不规则的撕裂口,断裂口呈锯齿状,且底板没有产生比较规则的裂纹(A);同时底板与侧板的角隅处也产生了长条状枣核形状的裂口,裂口周围向外凸起(B)。图11(c)为工况3钢箱的变形情况,距离炸药较近的角隅处发生严重破坏,断裂口周围钢板发生扭曲变形(A);同时发现,在远离炸药的3个钢板的角隅处发生了撕裂破坏(B)。结合各工况下钢箱失效变形的情况,可得炸药在钢箱角隅处爆炸时破坏最为严重。

      图  11  3种工况下,钢箱表面局部破坏示意图

      Figure 11.  Surface partial destruction diagrams of the steel box under three different conditions

    • 采用光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)方法对钢箱内部爆炸的过程进行数值模拟,研究了爆炸过程中钢箱表面的变形,分析了相同当量炸药在钢箱内部不同位置爆炸时钢箱表面的失效模式。

      (1)通过比较实验与SPH方法数值模拟获得的结果,数值模拟与实验所得的失效变形模式吻合较好,实验和数值模拟得到的钢箱表面挠度值相对误差为3.5%,验证了本文SPH方法在模拟钢箱内部爆炸问题上的有效性;

      (2)采用SPH方法数值模拟钢箱内部爆炸的过程,得出与实验一致的变形结果。同时也发现,在钢箱内部爆炸过程中,钢箱角隅处出现超压和应力集中现象,这与实验吻合较好。

      (3)进一步研究了相同当量炸药在钢箱内部不同位置爆炸时钢箱失效破坏情况,发现炸药在钢箱角隅处爆炸时钢箱破坏最严重。

      以上研究结果可为工事结构的防护设计提供参考。

参考文献 (16)

目录

    /

    返回文章
    返回