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

圆筒装置内爆炸压力载荷特性实验研究

李芝绒 王胜强 蒋海燕 张玉磊 袁建飞

引用本文:
Citation:

圆筒装置内爆炸压力载荷特性实验研究

    作者简介: 李芝绒(1968- ),女,研究员,lzr204@163.com;
    通讯作者: 王胜强, wangsheng1129@21cn.com
  • 中图分类号: 0382.4

Experimental studies on characteristics of explosion pressure loadin cylinder apparatus

    Corresponding author: WANG Shengqiang, wangsheng1129@21cn.com ;
  • CLC number: 0382.4

  • 摘要: 为了研究密闭环境压力载荷波形特征和分布规律,开展了长径比2∶1圆筒装置的TNT装药内爆炸实验,获取了圆筒壁面、端盖板表面的压力载荷数据,结合仿真计算结果,研究了圆筒壁面、端盖板上压力载荷波形特征和形成机理,建立了压力载荷峰值的表征模型,并开展了模型验证。结果表明,圆筒壁面与盖板压力载荷波形不完全相同,壁面压力载荷按照爆心距由近到远,从显著的单波峰向多峰值变化,在近场区域压力载荷衰减特性与一定炸高自由地面冲击波类似;盖板上压力载荷呈多波峰特点,中心区域压力载荷最大峰值大于首峰值的3倍,角隅区域大于首峰值6倍;圆筒壁面上的压力载荷最大峰呈现凹型分布规律。该研究结果可为密闭环境内爆载荷的特性分析和结构毁伤评估提供参考。
  • 图 1  圆筒爆炸装置及测点布设位置

    Figure 1.  Cylinder explosion apparatus and locations of measuring points

    图 2  120 g TNT装药圆筒壁面测点的压力载荷曲线

    Figure 2.  Pressure histories of 4 measuring points on cylinder wall at charge of 120 g TNT

    图 3  120 g TNT装药圆筒盖板测点的压力载荷曲线

    Figure 3.  Pressure histories of four measuring points on cover plate at charge of 120 g TNT

    图 4  40、80 g TNT装药在圆筒盖板P7测点的压力载荷曲线

    Figure 4.  Pressure histories of P7 measuring points on cover plate at charge of 40 g, 80 g TNT

    图 5  120 g TNT装药圆筒装置内爆炸冲击波压力云图

    Figure 5.  Pressure cloud diagram of cylinder internal explosion blast wave at charge of 120 g TNT

    图 6  120 g TNT装药圆筒装置测点P6、P7的冲击波压力曲线

    Figure 6.  Pressure histories of P6 and P7 measuring points at charge of 120 g TNT

    图 7  圆筒结构冲击波压力载荷曲线

    Figure 7.  Pressure histories of cylinder

    表 1  40、80、120 g TNT装药圆筒壁面、端盖的压力载荷峰值测量结果

    Table 1.  Measurement results of peak pressure of the cylinder wall and cover plate at charge of 40 g, 80 g and 120 g TNT

    测量点首峰值/MPa最大峰值/MPa首峰值/MPa最大峰值/MPa首峰值/MPa最大峰值/MPa
    40 g TNT80 g TNT120 g TNT
    P14.174.177.397.399.379.37
    P24.104.105.155.155.375.37
    P31.241.241.791.792.732.73
    P40.682.681.093.091.523.52
    P50.412.46 0.6772.341.203.34
    P6 0.3362.22 0.5362.661.153.85
    P7 0.3862.35 0.5552.791.183.67
    P8 0.3954.13 0.5175.091.197.63
    下载: 导出CSV
  • [1] ZYSKOWSKI A, SOCHET I, MAVROT G, et al. Study of the explosion process in a small scale experiment—structural loading [J]. Journal of Loss Prevention in the Process Industries, 2004, 17(4): 291–299. DOI: 10.1016/j.jlp.2004.05.003.
    [2] 黄雪峰, 王伟力, 姜颖资, 等. 穿爆战斗部装药对舰艇组合舱室内爆炸毁伤效应研究 [J]. 弹箭与制导学报, 2014, 34(6): 75–82. DOI: 10.3969/j.issn.1673-9728.2014.06.020.
    HUANG Xuefeng, WANG Weili, JIANG Yingzi, et al. Research on damage effect of semi-armor-piercing warhead explosion in combined cabin of vessel [J]. Journal of Projectiles Rockets Missiles and Guidance, 2014, 34(6): 75–82. DOI: 10.3969/j.issn.1673-9728.2014.06.020.
    [3] 连赟猛, 吴美云, 顾晓辉. 密闭装置内爆炸及数值模拟精度分析 [J]. 兵器装备工程学报, 2016, 37(8): 169–172. DOI: 10.11809/scbgxb2016.08.038.
    LIAN Yunmeng, WU Meiyun, GU Xiaohui. Reliability analysis on closed space due to internal explosion [J]. Journal of Ordnance Equipment Engineering, 2016, 37(8): 169–172. DOI: 10.11809/scbgxb2016.08.038.
    [4] 樊壮卿, 王伟力, 黄雪峰, 等. 典型舱室内爆炸仿真分析 [J]. 工程爆破, 2015, 21(3): 13–17. DOI: 10.3969/j.issn.1006-7051.2015.03.004.
    FAN Zhangqing, WANG Weili, HUANG Xuefeng, et al. Simulation analysis on typical cabin internal explosion [J]. Engineering Blasting, 2015, 21(3): 13–17. DOI: 10.3969/j.issn.1006-7051.2015.03.004.
    [5] 孔祥韶, 吴卫国, 李晓彬, 等. 舰船舱室内部爆炸的数值模拟研究 [J]. 中国舰船研究, 2009, 4(4): 7–11. DOI: 10.3969/j.issn.1673-3185.2009.04.002.
    KONG Xiangshao, WU Weiguo, LI Xiaobin, et al. Numerical simulation of cabin structure under inner explosion [J]. Chinese Journal of Ship Research, 2009, 4(4): 7–11. DOI: 10.3969/j.issn.1673-3185.2009.04.002.
    [6] 段雷琳, 高轩能, 江媛. 内爆炸下球面钢网壳结构的冲击波超压计算分析 [J]. 华侨大学学报(自然科学版), 2013, 34(5): 557–562. DOI: 10.11830/ISSN.1000-5013.2013.05.0557.
    DUAN Leilin, GAO Xuanning, JIANG Yan. Analysis for shock wave overpressure of spherical steel reticulated shell under internal explosion [J]. Journal of Huaqiao University (Natural Science), 2013, 34(5): 557–562. DOI: 10.11830/ISSN.1000-5013.2013.05.0557.
    [7] 姚熊亮, 屈子悦, 姜子飞, 等. 舰船舱内爆炸载荷特征与板架毁伤规律分析 [J]. 中国舰船研究, 2018, 13(3): 140–148. DOI: 10.19693/j.issn.1673-3185.01162.
    YAO Xiongliang, QU Ziyue, JIANG Zifei, et al. Analysis on characteristics of blast loading and stiffened plate damage due to internalblast in ship [J]. Chinese Journal of Ship Research, 2018, 13(3): 140–148. DOI: 10.19693/j.issn.1673-3185.01162.
    [8] 侯海量, 朱锡, 李伟, 等. 舱内爆炸冲击载荷特性实验研究 [J]. 船舶力学, 2010, 14(8): 902–908. DOI: 10.3969/j.issn.1007-7294.2010.08.011.
    HOU Hailiang, ZHU Xi, LI Wei, et al. Experimental studies on characteristics of blast loading when exploded inside ship cabin [J]. Journal of Ship Mechanics, 2010, 14(8): 902–908. DOI: 10.3969/j.issn.1007-7294.2010.08.011.
    [9] 李伟, 朱锡, 梅志远, 等. 战斗部舱内爆炸对舱室结构毁伤的实验研究 [J]. 舰船科学技术, 2009, 31(3): 34–37. DOI: 10.3404/j.issn.1672-7649.2009.03.005.
    LI Wei, ZHU Xi, MEI Zhiyuan, et al. Experimental studies on damage effect of missile warhead on cabin, structure under internal explosion [J]. Ship Science and Technology, 2009, 31(3): 34–37. DOI: 10.3404/j.issn.1672-7649.2009.03.005.
    [10] 北京工业学院八系. 爆炸及其作用[M]. 北京: 国防工业出版社, 1979: 400−420.
    [11] 易仰贤. 空爆冲击波马赫反射近似计算 [J]. 爆炸与冲击, 1983, 3(2): 44–49.
    YI Yangxian. Approximate calculation of mach, reflection of explosive shock waves in air [J]. Explosion and Shock Waves, 1983, 3(2): 44–49.
    [12] ANDERSON C E, BAKER W E, WAUTERS D K, et al. Quasi-static pressure, duration, and impulse for explosions (e.g. HE) in structures [J]. International Journal of Mechanical Sciences, 1983, 25(6): 45–464. DOI: 10.1016/0020-7403(83)90059-0.
  • [1] 谭成文王富耻李树奎王鲁 . 内爆炸加载条件下圆筒的膨胀、破裂规律研究. 爆炸与冲击, 2003, 23(4): 305-308.
    [2] 高洪泉卢芳云王少龙颜澎罗永锋袁伟胡建 . 具有爆炸反应特性装置在密闭容器内爆炸威力的评价方法. 爆炸与冲击, 2011, 31(3): 306-310. doi: 10.11883/1001-1455(2011)03-0306-05
    [3] 王宝兴马克辛张银花 . 亚麻尘在1米~3试验装置内的爆炸泄压研究. 爆炸与冲击, 1991, 11(1): 52-57.
    [4] 张玉磊苏健军李芝绒蒋海燕仲凯王胜强 . TNT内爆炸准静态压力特性. 爆炸与冲击, 2018, 38(6): 1429-1434. doi: 10.11883/bzycj-2017-0170
    [5] 范亚夫段祝平马岳陈捷 . 钛合金TC4材料热处理工艺对圆筒在内爆炸载荷下膨胀半径的影响. 爆炸与冲击, 2004, 24(4): 337-346.
    [6] 樊成飞王耀华王强 . 爆燃动力装置内弹道的优化设计. 爆炸与冲击, 2015, 35(2): 267-272. doi: 10.11883/1001-1455(2015)02-0267-06
    [7] 王志荣蒋军成周超 . 连通装置气体爆炸特性实验. 爆炸与冲击, 2011, 31(1): 69-74. doi: 10.11883/1001-1455(2011)01-0069-06
    [8] 尹晓春冯思艳龚育宁 . 利用 Laplace 变换确定双层厚壁长圆筒的轴对称界面碰撞压力. 爆炸与冲击, 1998, 18(4): 365-374.
    [9] 蔡正宇丁圆圆王士龙郑志军虞吉林 . 梯度多胞牺牲层的抗爆炸分析. 爆炸与冲击, 2017, 37(3): 396-404. doi: 10.11883/1001-1455(2017)03-0396-09
    [10] 路国运雷建平武勇忠张善元 . 多跨薄壁压力管道侧向冲击破坏的实验研究. 爆炸与冲击, 2003, 23(5): 454-459.
    [11] 张博孜王浩王珊珊 . 含扰流装置的中心燃气式抛撒内弹道过程的数值模拟. 爆炸与冲击, 2015, 35(2): 208-214. doi: 10.11883/1001-1455(2015)02-0208-07
    [12] 刘文祥张德志钟方平程帅张庆明 . 球形爆炸容器内炸药爆炸形成的准静态气体压力. 爆炸与冲击, 2018, 38(5): 1045-1050. doi: 10.11883/bzycj-2017-0056
    [13] 薛冰马宏昊沈兆武余勇 . 爆炸容器内小药量实验动态标定压力传感器. 爆炸与冲击, 2015, 35(3): 437-441. doi: 10.11883/1001-1455(2015)03-0437-05
    [14] 姚哲芳任辉启沈兆武 . 柱状装药爆炸条件下厚壁圆筒爆室内 流场及结构动力响应分析. 爆炸与冲击, 2012, 32(5): 449-456. doi: 10.11883/1001-1455(2012)05-0449-08
    [15] 李忠盛吴护林陈韵如易同斌李立潘复生 . 内爆炸载荷作用下7A55铝合金的动态性能及断裂行为. 爆炸与冲击, 2012, 32(2): 190-195. doi: 10.11883/1001-1455(2012)02-0190-06
    [16] 李永池董杰胡秀章范书群 . 内聚爆炸载荷下热塑性中厚球壳的变形和层裂. 爆炸与冲击, 2005, 25(4): 319-324. doi: 10.11883/1001-1455(2005)04-0319-06
    [17] 吴林杰侯海量朱锡陈鹏宇田万平 . 水下接触爆炸下防雷舱舷侧空舱的内压载荷特性. 爆炸与冲击, 2017, 37(4): 719-726. doi: 10.11883/1001-1455(2017)04-0719-08
    [18] 王海福冯顺山 . 爆炸载荷下聚氨酯泡沫材料中冲击波压力特性. 爆炸与冲击, 1999, 19(1): 78-83.
    [19] 徐胜利岳鹏涛彭金华 . 多爆源云雾爆炸波相互作用的三维数值研究. 爆炸与冲击, 2000, 20(1): 1-6.
    [20] 朱振海杨永琦 . 多火花式动态光弹性仪在爆炸力学实验中的初步应用. 爆炸与冲击, 1985, 5(3): 67-76.
  • 加载中
图(7)表(1)
计量
  • 文章访问数:  67
  • HTML全文浏览量:  298
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-08-31
  • 录用日期:  2019-02-17
  • 刊出日期:  2019-10-01

圆筒装置内爆炸压力载荷特性实验研究

    作者简介:李芝绒(1968- ),女,研究员,lzr204@163.com
    通讯作者: 王胜强, wangsheng1129@21cn.com
  • 西安近代化学研究所,陕西 西安 710065

摘要: 为了研究密闭环境压力载荷波形特征和分布规律,开展了长径比2∶1圆筒装置的TNT装药内爆炸实验,获取了圆筒壁面、端盖板表面的压力载荷数据,结合仿真计算结果,研究了圆筒壁面、端盖板上压力载荷波形特征和形成机理,建立了压力载荷峰值的表征模型,并开展了模型验证。结果表明,圆筒壁面与盖板压力载荷波形不完全相同,壁面压力载荷按照爆心距由近到远,从显著的单波峰向多峰值变化,在近场区域压力载荷衰减特性与一定炸高自由地面冲击波类似;盖板上压力载荷呈多波峰特点,中心区域压力载荷最大峰值大于首峰值的3倍,角隅区域大于首峰值6倍;圆筒壁面上的压力载荷最大峰呈现凹型分布规律。该研究结果可为密闭环境内爆载荷的特性分析和结构毁伤评估提供参考。

English Abstract

  • 伴随着精确制导技术的发展,导弹战斗部的毁伤效能成为了战斗部研制和使用部门极为关切和重点关注的问题,特别是反舰型、侵爆型战斗部。了解和掌握战斗部的内爆炸压力载荷特征及分布规律,是检验和评价战斗部毁伤效能的关键。

    在密闭环境内,由于壁面结构的约束,爆炸冲击波与壁面结构的耦合效应显著,产生的反射冲击波在壁面间来回反射,形成了比开放环境下更加复杂的压力波流场。目前,由于对内爆炸压力载荷特性的认识尚不完善,在理论上很难建立一个合理的内爆炸载荷解析模型。学者们主要通过数值模拟仿真和模拟实验的方法开展研究:Zyskowski等[1]开展了小比例模型内爆炸试验和仿真计算,研究内爆炸压力载荷波形特征;黄雪峰等[2]、连赟猛等[3]、樊壮卿等[4]、孔祥韶等[5]开展舱室内爆炸仿真,分析舰船舱室内冲击波流场时空变化和壁面应力变化规律;段雷琳等[6]开展大跨球面钢网壳结构压力载荷的数值模拟计算研究;姚熊亮等[7]开展的舰船舱内爆炸仿真计算,研究冲击波压力载荷波形特征与板架结构的毁伤规律。由于在数值仿真计算中缺乏准确的炸药爆轰产物状态方程,仿真结果只能够定性地表征压力载荷特性,定量结果有一定的误差。模拟内爆炸试验是一种最直观的方法,通过测试手段,获取爆炸压力载荷数据:侯海量等[8]、李伟等[9]开展了舱室内爆炸冲击载荷特性实验研究,分析了爆炸冲击波流场和压力载荷的作用过程以及舱室结构毁伤特性。在内爆炸环境下,特别是钢质立方体模型结构,可能产生大变形或破坏,通常一次使用,试验性价比较高,且由于结构的破坏,常常难以获取完整的测量数据。因此,开展实体模型的内爆炸试验研究较少,可参考的内爆炸载荷特性研究少见文献报道。

    本文中,设计一种能够多次使用的圆筒型爆炸装置,开展质量分别为40、80、120 g TNT装药内爆炸实验,获取壁面、盖板上的压力载荷数据,并结合仿真计算结果,分析研究内爆炸压力载荷的波形特征和分布规律,拟为内爆炸压力载荷特性研究提供参考。

    • 圆筒爆炸装置如图1(a)所示,主要由圆筒、圆筒两端的法兰盘、前端盖、后端盖、压环、密封圈组成。圆筒长1 670 mm,内径800 mm,壁厚12 mm;前端盖厚度12 mm,后端盖厚度1.5 mm。法兰盘上均布16个螺栓孔,通过螺杆将圆筒、前后端盖连接,形成了密封装置。装置材料均采用Q235钢。

      图  1  圆筒爆炸装置及测点布设位置

      Figure 1.  Cylinder explosion apparatus and locations of measuring points

      在圆筒的经线和端盖的中心线上,分别设置了压力传感器安装孔,圆筒壁面间距200 mm,前端盖间距100 mm,如图1(b)所示。压力传感器选用PCB公司的113B24、113B22型压力传感器,量程分别为6.9、34.5 MPa。113B24型传感器安装于前端盖,113B22型传感器安装于圆筒壁面,传感器敏感面距圆筒壁面或盖板表面约0.8 mm。传感器安装完成后,在传感器敏感面与安装孔形成的凹槽内涂满油脂,并使油脂表面与圆筒壁面的弧度一致,或与盖板表面平齐。设置油脂的目的是延缓爆炸热向传感器敏感头传递的时间,减小压力传感器测压过程中产生的热冲击信号对压力测量信号的干扰。

      实验装药为柱形TNT,质量分别为40、80、120 g,长径比近似1∶1。实验时,TNT装药放置于圆筒装置几何中心的木质弹架上,药柱几何中心的径向平面与圆筒壁面测点连线在同一平面,由电雷管起爆装药。

    • 表1是质量分别为40、80、120 g TNT装药在圆筒壁面、端盖的压力载荷峰值测量结果。

      测量点首峰值/MPa最大峰值/MPa首峰值/MPa最大峰值/MPa首峰值/MPa最大峰值/MPa
      40 g TNT80 g TNT120 g TNT
      P14.174.177.397.399.379.37
      P24.104.105.155.155.375.37
      P31.241.241.791.792.732.73
      P40.682.681.093.091.523.52
      P50.412.46 0.6772.341.203.34
      P6 0.3362.22 0.5362.661.153.85
      P7 0.3862.35 0.5552.791.183.67
      P8 0.3954.13 0.5175.091.197.63

      表 1  40、80、120 g TNT装药圆筒壁面、端盖的压力载荷峰值测量结果

      Table 1.  Measurement results of peak pressure of the cylinder wall and cover plate at charge of 40 g, 80 g and 120 g TNT

      图2为120 g TNT装药圆筒壁面4个测点获取的冲击波压力载荷曲线。

      图  2  120 g TNT装药圆筒壁面测点的压力载荷曲线

      Figure 2.  Pressure histories of 4 measuring points on cylinder wall at charge of 120 g TNT

      图2表1可以看出,在圆筒壁面上,按照爆心距由近到远,压力载荷曲线呈现由显著单波峰向多波峰变化的特点,在爆心近场区域,冲击波首峰值为最大峰值。随着爆心距增大,首峰值逐渐减小,第2波峰逐渐出现,在角隅区域,第2波峰甚至超过第1峰值成为最大峰。从壁面中心到角隅,压力载荷的最大峰值逐渐衰减,然后又增大,呈现了凹形的变化规律。

      图3为120 g TNT装药在前端盖板4个测点获取的冲击波压力载荷曲线。

      图  3  120 g TNT装药圆筒盖板测点的压力载荷曲线

      Figure 3.  Pressure histories of four measuring points on cover plate at charge of 120 g TNT

      图3表1可以看出,圆筒壁面和端盖上的压力载荷曲线特征不完全相同,端盖板上各测点压力载荷曲线呈现多波峰的特点,首峰值较小,第2峰值最大。各装药4个测点的首峰值基本相当。除了P8测点外,各装药P5、P6、P7测点的最大峰值也基本相当。测点P8最大峰值比其他3个测点的最大峰值大,约等于其他3个测点平均值的1.7~2.1倍。从端盖中心到角隅,压力载荷最大峰值基本相当,当接近角隅区域时最大峰值突跃增大。

      图4是40、80 g TNT装药在端盖P7测点的压力载荷曲线,结合图3(c)的120 g TNT装药P7测点的压力载荷曲线,可观察到端盖上压力载荷曲线呈现周期性的压力波动,40 g TNT前3个压力波动周期分别为2.48、2.64、2.98 ms,80 g TNT分别为1.95、2.45、2.58 ms,120 g TNT分别为1.75、2.30、2.41 ms。各装药的压力波动周期随时间的增长而逐渐变大,后一周期的最大峰值小于前一周期的最大峰值,大约3~4个周期后,压力波趋于匀化。压力波动周期也随炸药质量的变化而变化,炸药质量越大,压力波动周期越短,质量越小,压力波动周期越长。

      图  4  40、80 g TNT装药在圆筒盖板P7测点的压力载荷曲线

      Figure 4.  Pressure histories of P7 measuring points on cover plate at charge of 40 g, 80 g TNT

    • 为了进一步分析圆筒装置内压力载荷的产生机理,开展了该工况的内爆炸仿真计算。应用有限元分析软件AUTODYN的Euler算法,计算炸药装药爆炸冲击波的形成、传播过程。按照试验测点位置,在圆筒壁面、盖板设置观察点,获取观察点冲击波压力载荷随时间的变化曲线。

      在仿真计算中,依据圆筒装置几何尺寸及对称性,建立了1/2模型。考虑到炸药质量小,装置各组成的厚度及强度,圆筒壁面、前端盖采用了实体单元模型;后端盖为薄钢板,采用shell壳单元模型,圆筒壁与盖板之间采用jion连接。炸药和空气场采用Euler网格划分单元。炸药采用JWL状态方程,空气场采用理想气体状态方程,在装置结构与炸药、空气之间采用流固耦合的计算方法。

      图5是120 g TNT装药圆筒装置内冲击波压力云图的变化过程。

      图  5  120 g TNT装药圆筒装置内爆炸冲击波压力云图

      Figure 5.  Pressure cloud diagram of cylinder internal explosion blast wave at charge of 120 g TNT

      图5可以看出,炸药装药爆炸输出的爆炸冲击波先与圆筒壁面耦合反射(见图5(c)),产生的反射波部分向爆炸中心传播,部分沿圆筒壁面向盖板传播。向爆心传播的壁面反射波在炸点位置区域汇聚(见图5(d)),形成极度压缩的汇聚波,然后膨胀,向周围扩散,形成了两头强、中间弱的梭形二次冲击波向四周传播(见图5(e))。向端盖方向传播的二次冲击波,跟随于爆炸冲击波之后。端盖的中心区域先受到爆炸冲击波的作用(见图5(e)),然后是壁面反射波、二次冲击波、盖板反射冲击波的叠加,形成了复杂压力波与盖板耦合反射,形成的复杂反射波向圆筒中心传播(见图5(f))。在两角角隅区域,先是沿壁面传播的马赫反射波在角隅耦合反射,然后是角隅的反射波、二次冲击波、壁面二次反射波的叠加,在角隅耦合反射,形成复杂反射波向圆筒中心传播。在圆筒壁面、端盖上产生的复杂反射波再次在炸点位置区域再次汇聚(见图5(h)(i))、压缩,然后膨胀,形成三次冲击波向四周扩散(图5(j)),再次与圆筒壁、端盖耦合反射(见图5(k)(l))。如此反复多次,冲击波能逐渐转化为气体内能,压力峰值逐渐减小,准静态压力逐渐增大,直至形成压力稳定、均匀的准静态压力场。

      图5还可以看出,圆筒左右部分的压力云图不对称,主要是由于后端盖为1.5 mm的薄钢板,在冲击波压力载荷的作用下,板面先从应力最大的两角角隅产生塑性变形,并逐渐向端盖中心扩展。由于板面的外凸,使前、后端盖反射波传播的距离、时间不同步,形成了冲击波流场左右不对称现象。

      图6是仿真计算获取的40 g TNT在盖板P6、P7测点位置的冲击波压力载荷曲线。P6测点压力载荷最大峰值为2.15 MPa,P7测点压力载荷最大峰值为1.87 MPa。与实验对应点的压力载荷对比,测点计算值比实验值小,并且压力波峰的上升前沿比实验曲线缓慢,但是压力波峰出现的时刻与实验基本相同。产生这些差异的主要原因,一是仿真计算中材料模型及参数数值不够准确,二是建模网格密度不高。

      图  6  120 g TNT装药圆筒装置测点P6、P7的冲击波压力曲线

      Figure 6.  Pressure histories of P6 and P7 measuring points at charge of 120 g TNT

    • 从圆筒壁面所受到的压力载荷曲线可以看出,距爆心较近的近场区域,压力载荷曲线呈现单波峰特点。在受到二次冲击波作用之前,爆炸冲击波在壁面传播特性,与地面环境中一定炸高炸药装药爆炸冲击波在地面传播规律类似。因此,圆筒壁面压力载荷的首峰值,可依据爆炸冲击波在地面传播模型计算[10]。当入射角小于39.14°时,冲击波产生规则反射。在规则反射区内,当冲击波产生正反射时,冲击波压力首峰值角度为:

      $\Delta {p_{\rm{2}}} = 2\Delta {p_{\rm{1}}} + \frac{{6\Delta {p^2_{1}}}}{{7{p_0} + \Delta {p_{\rm{1}}}}}$

      式中:$\Delta {p_{\rm{1}}}$为炸药在无限空气中超过空气初始压力的冲击波压力峰值,可按照萨道夫斯基公式或工程设计规范中的经验公式得到;${p_{\rm{0}}}$为空气初始压力。

      当冲击波产生斜反射时,冲击波压力首峰$\Delta {p_{\rm{2}}}$为:

      $\Delta {p_2} = (1 + \cos \alpha )\Delta {p_{\rm{1}}} + \frac{{6\Delta {p_{\rm{1}}^{\rm{2}}}}}{{7{p_0} + \Delta {p_{\rm{1}}}}}{\cos ^2}\alpha $

      当入射角大于39.14°时,在地面环境,冲击波产生马赫反射现象。马赫反射区冲击波参数的计算相对复杂,参考易仰贤[11]的马赫反射压近似计算公式(3)~(4),得到:

      $\Delta {p_3} = \frac{1}{6}{p_0}\left( {{A^2}\left( {\frac{6}{{{p_0}}}\Delta {p_1} + 7} \right) - 7} \right)$

      $A = \frac{1}{{\sin \alpha }} + \left( {\left( {1 - \frac{{r^*}}{r}} \right)\frac{{{{\cos }^2}\alpha }}{{\cos (\alpha - \alpha^ *)}} + \frac{{r^*}}{r}\tan(\alpha - \alpha^ *)\cos \alpha } \right)$

      式中:$\alpha $为冲击波入射角;r*、α*分别为马赫反射起点处的爆心距、入射角;r为测点处的爆心距。

      将模型计算结果与试验测量结果比对,对于规则反射区,正反射点冲击波压力载荷首峰值的偏差小于5%,斜反射点压力载荷首峰值的偏差小于14%;对于马赫反射区,除了40 g TNT偏差较大外,其他测点压力载荷首峰值的偏差小于13%。

    • 从圆筒端盖所受到的压力载荷曲线和仿真计算压力云图可以看出,端盖上先受到了爆炸冲击波的作用,与圆筒壁面类似。在受到其他冲击波作用之前,爆炸冲击波的传播特性,与地面环境中一定炸高炸药装药爆炸冲击波在地面的传播规律类似。由于端盖各测点产生规则反射,因此,各测点压力载荷的首峰值可依据入射角用式(1)或式(2)计算。比对计算结果与试验测量结果,测点P5首峰的偏差小于10%;斜反射测点P6、P7、P8首峰的偏差的小于16%。

      表1的实验结果可知,40、80、120 g TNT装药在端盖中心区域的最大峰值等于各自首峰值的6.23、4.48、3.07倍,在角隅区域的最大峰值等于各首峰值的10.46、9.85、6.52倍。由图5可以看出,由于压力载荷的最大压力峰值是多波耦合叠加结果,形成机理复杂,目前还无法建立一个合理的解析模型。因此,依据爆炸相似定律,按照实验工况和实验炸药装药质量,可以认为:在长径比为2∶1的圆筒端盖上,在冲击波入射角小于13.5°的中心区域,对比距离为1.74~2.50时,压力载荷的最大峰值大于首峰值的3倍;在冲击波入射角为19.7°~25.6°的角隅区域,对比距离为1.80~2.70时,压力载荷的最大峰值大于首峰值的6倍。

    • 从获取的圆筒装置所受到的压力载荷曲线可以看出,冲击波经过3~4个周期性波动后逐渐匀化,在圆筒内形成了稳定、均匀的准静态压力,如图7所示。冲击波最大峰值的周期性衰减特性,大致符合Anderson等[12]得到的反射冲击波峰值为前一个反射冲击波峰值的一半的规律。主要差异表现为,在第2周期以后,后一个冲击波最大峰值比前一个冲击波最大峰值的一半略大。主要原因为:一方面是圆筒内的准静态压力提高了圆筒的整体压力,增大了压力载荷峰值的绝对值;另一方面可能是由测量系统的测量误差引起。如采用的PCB压力传感器低频响应特性(−5%)为0.005 Hz,用于测量准静态压力时会产生测量误差;如爆炸热持续作用于压力传感器的敏感面上,传感器输出的热响应信号叠加到压力测量信号中,造成了准静态压力的测量误差。

      图  7  圆筒结构冲击波压力载荷曲线

      Figure 7.  Pressure histories of cylinder

    • 通过在圆筒型装置内的爆炸试验,分析了长径比2∶1圆筒壁面、端盖上压力载荷特征,结论如下。

      (1)圆筒壁面的压力载荷按照爆心距由近到远,从显著的单波峰向多峰值变化。在近场区域压力载荷衰减特性与一定炸高自由地面冲击波传播特性类似,压力载荷的首峰值可按照规则反射区和马赫反射区的压力峰值计算公式得到。

      (2)端盖上压力载荷呈多波峰特点。在中心区域,对比距离为1.74~2.50时,压力载荷的最大峰值大于首峰值的3倍;在角隅区域,对比距离为1.80~2.70时,压力载荷的最大峰值大于首峰值的6倍。

      (3)在圆筒壁面和盖板上,从几何中心到角隅区域,压力载荷峰值呈凹型分布特征;在盖板上,从几何中心到角隅区域,压力载荷最大峰值基本相当,当接近角隅区域时最大峰值突跃增大。

      (4)圆筒内冲击波经过3~4个周期性波动后,冲击波压力逐渐匀化,形成了稳定、均匀的准静态压力。

      虽然用试验和数值模拟的方法,验证了TNT装药在该工况和有限质量范围内压力载荷的特征和分布规律。但对于其他比例模型的圆筒装置,本文的结论是否适合,还有待于进一步验证。由于内爆炸冲击波的反射及叠加效应复杂,今后还需要进一步研究。

参考文献 (12)

目录

    /

    返回文章
    返回