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

装药动爆冲击波特性研究

陈龙明 李志斌 陈荣

引用本文:
Citation:

装药动爆冲击波特性研究

    作者简介: 陈龙明(1996- ),男,硕士研究生,chenlongming14@nudt.edu.cn;
    通讯作者: 李志斌, lizhibin@nudt.edu.cn
  • 中图分类号: O382.1

Characteristics of dynamic explosive shock wave of moving charge

    Corresponding author: LI Zhibin, lizhibin@nudt.edu.cn ;
  • CLC number: O382.1

  • 摘要: 导弹、炮弹等战斗部爆炸时具有一定的速度,较大的运动速度会使爆炸冲击波场分布发生变化,进而对弹药的毁伤威力产生影响。本文中采用AUTODYN软件对速度分别为0、272、340、680、1 020和1 700 m/s的TNT球形裸装药在空气中爆炸的冲击波场进行了仿真计算,定量研究装药在动爆条件下的峰值超压、比冲量和正压作用时间等威力参数特性。结果表明,方位角小于90°时装药速度与冲击波超压、比冲量成正相关,与正压作用时间成负相关;方位角大于90°时装药速度与冲击波超压、比冲量成负相关,与正压作用时间成正相关。超压峰值大小沿方位角成正弦变化。最后,分析了冲击波峰值超压数据,建立了动爆冲击波超压的计算模型,该模型计算结果与仿真和实验结果吻合较好。
  • 图 1  参数设置

    Figure 1.  Parameter setting

    图 2  t=0.4 ms时不同速度下的冲击波压力云图

    Figure 2.  Pressure contour of shock waves under different velocities at t=0.4 ms

    图 3  不同比例距离时超压峰值随方位角的变化

    Figure 3.  Peak overpressure Δp changes with θ under different $\overline R $

    图 4  不同方向的增大系数

    Figure 4.  Increasing ratios r in different directions

    图 5  不同比例距离时比冲量随方位角的变化

    Figure 5.  Specific impluse I changes with θ under different $\overline R $

    图 6  不同比例距离时正压时间随方位角的变化

    Figure 6.  Position pressure time τ changes with θ under different $\overline R $

    图 7  拟合曲线f (v)

    Figure 7.  Fitting curve f (v)

    图 8  拟合曲线f ($\overline R $)

    Figure 8.  Fitting curve f ($\overline R $)

    图 9  拟合曲线f (θ)

    Figure 9.  Fitting curves f (θ)

    表 1  超压数据与公式计算结果比较

    Table 1.  Comparison of peak overpressure results with theoretical values

    Δp/MPaη/%
    模拟 Henrych公式
    2.492.384.41
    0.970.943.78
    0.470.471.07
    0.280.272.22
    下载: 导出CSV

    表 2  仿真计算结果对比

    Table 2.  Comparison of simulation results

    v/(m·s−1)$\overline R $/(m·kg−1/3)θ/(°)Δp/kPaη/%
    公式静爆模拟
    1000.72 151 609.031 516.351 663.803.29
    1001.81 15214.91219.84236.429.09
    1001.26 73394.31480.98493.1320.04
    1000.541641 155.982 552.412 322.7950.23
    5000.72 152 324.811 516.352 367.871.82
    5001.81 15310.52219.84317.832.30
    5001.26 73569.71480.98533.516.79
    5000.541641 670.222 552.411 579.855.72
    1 5000.72 155 551.671 516.355 249.825.75
    1 5001.81 15741.52219.84690.907.33
    下载: 导出CSV

    表 3  实验结果对比

    Table 3.  Comparison of experimental results

    v/(m·s−1)$\overline R $/(m·kg−1/3)θ/(°)Δp/kPaη/%
    实验公式
    534.311.37 15627.42638.32 1.71
    534.311.37 45565.37572.21 1.20
    534.311.37105361.97375.33 3.56
    587.651.19 15751.53926.0318.84
    587.651.19 45599.84830.1327.74
    587.651.19105379.21544.5030.36
    下载: 导出CSV
  • [1] 杜红棉, 王燕, 祖静, 等. 导弹动爆对模拟船舱毁伤效果试验研究 [J]. 高压物理学报, 2011, 25(3): 261–267. DOI: 10.11858/gywlxb.2011.03.011.
    DU H M, WANG Y, ZU J, et al. Experimental study on damage effect of simulated ship cabin subjected to missile dynamic explosion [J]. Chinese Journal of High Pressure Physics, 2011, 25(3): 261–267. DOI: 10.11858/gywlxb.2011.03.011.
    [2] 高元浩. 榴弹战斗部动爆条件下毁伤评估方法研究[D]. 太原: 中北大学, 2016.
    [3] 卢芳云, 李翔宇, 林玉亮. 战斗部结构与原理[M]. 北京: 科学出版社, 2009.
    [4] 赵岩. 动爆冲击波超压测试方法研究[D]. 太原: 中北大学, 2012.
    [5] 张光莹, 周旭, 黄咏政, 等. 动爆冲击波特性分析方法研究[C] // 第四届全国计算爆炸力学会议论文集. 2008: 282−287.
    [6] 蒋海燕, 李芝绒, 张玉磊, 等. 运动装药空中爆炸冲击波特性研究 [J]. 高压物理学报, 2017, 31(3): 286–294. DOI: 10.11858/gywlxb.2017.03.010.
    JIANG H Y, LI Z R, ZHANG Y L, et al. Characteristics of air blast wave field for explosive charge moving at different velocities [J]. Chinese Journal of High Pressure Physics, 2017, 31(3): 286–294. DOI: 10.11858/gywlxb.2017.03.010.
    [7] 聂源, 蒋建伟, 李梅. 球形装药动态爆炸冲击波超压场计算模型 [J]. 爆炸与冲击, 2017, 37(5): 951–956. DOI: 10.11883/1001-1455(2017)05-0951-06.
    NIE Y, JIANG J W, LI M. Overpressure calculation model of sphere charge blasting with moving velocity [J]. Explosion and Shock Waves, 2017, 37(5): 951–956. DOI: 10.11883/1001-1455(2017)05-0951-06.
    [8] 杨鑫, 石少卿, 程鹏飞, 等. 爆炸冲击波在空气中传播规律的经验公式对比及数值模拟 [J]. 四川建筑, 2007, 27(5): 71–73. DOI: 10.3969/j.issn.1007-8983.2007.05.033.
    [9] ARMENDT B F, SPERRAZZA J. Air blast measurements around moving explosive charges, Part Ⅲ: AD0114950 [R]. Aberdeen: Army Ballistics Research Laboratory, 1956.
  • [1] 丁永红尤文斌马铁华 . 舰用动爆冲击波记录系统的设计与应用. 爆炸与冲击, 2013, 33(2): 194-199. doi: 10.11883/1001-1455(2013)02-0194-06
    [2] 林谋金崔晓荣马宏昊郑炳旭贾虎 . RDX基铝纤维炸药空中爆炸性能. 爆炸与冲击, 2016, 36(2): 230-235. doi: 10.11883/1001-1455(2016)02-0230-06
    [3] 李凌锋韦灼彬唐廷董琪刘靖晗邱艳宇 . 爆炸荷载下沉箱重力式码头模型毁伤效应. 爆炸与冲击, 2019, 39(1): 012202-1-012202-7. doi: 10.11883/bzycj-2017-0406
    [4] 薛冰马宏昊沈兆武余勇 . 爆炸容器内小药量实验动态标定压力传感器. 爆炸与冲击, 2015, 35(3): 437-441. doi: 10.11883/1001-1455(2015)03-0437-05
    [5] 王力韩峰陈放马晓青 . 偏心对称起爆战斗部破片初速的增益. 爆炸与冲击, 2016, 36(1): 69-74. doi: 10.11883/1001-1455(2016)01-0069-06
    [6] 程祥杨明郭亚丽韩焱 . 修正的Friedlander方程指数衰减因子. 爆炸与冲击, 2009, 29(4): 425-428. doi: 10.11883/1001-1455(2009)04-0425-04
    [7] 廖学燕沈兆武姚保学董树男马宏昊蒋耀港 . 含铝纤维复合炸药的能量输出和力学强度. 爆炸与冲击, 2010, 30(4): 424-428. doi: 10.11883/1001-1455(2010)04-0424-05
    [8] 林俊德 . 空中爆炸冲击波(Ⅱ)测量技术. 爆炸与冲击, 1986, 6(1): 89-96.
    [9] 乔登江 . 空中爆炸冲击波(Ⅰ)基本理论. 爆炸与冲击, 1985, 5(4): 78-85.
    [10] 张洪武何扬张昌权 . 空中爆炸冲击波地面荷载的数值模拟. 爆炸与冲击, 1992, 12(2): 156-165.
    [11] 王新颖王树山卢熹王建民 . 空中爆炸冲击波对生物目标的超压-冲量准则. 爆炸与冲击, 2018, 38(1): 106-111. doi: 10.11883/bzycj-2017-0031
    [12] 项大林荣吉利李健杨荣杰 . 金属壳体装药水下爆炸的冲击波特性. 爆炸与冲击, 2012, 32(1): 67-72. doi: 10.11883/1001-1455(2012)01-0067-06
    [13] 李润珊易嘉钰陶崇铸沈国光 . 触、近水面爆炸自由场冲击波特性的试验研究. 爆炸与冲击, 1991, 11(4): 331-338.
    [14] 姚成宝王宏亮浦锡锋寿列枫王智环 . 空中强爆炸冲击波地面反射规律数值模拟研究. 爆炸与冲击, 2019, 39(11): 112201-1-112201-8. doi: 10.11883/bzycj-2018-0287
    [15] 林俊德周忠海 . 一种空中核爆炸近区冲击波测量用的压传感器. 爆炸与冲击, 1982, 2(2): 64-70.
    [16] 钱胜国张伟林徐光耀 . 近自由水面水下爆炸时水中激波特性. 爆炸与冲击, 1983, 3(4): 53-63.
    [17] 宗琦 . 爆炸冲击波的动态破岩特性研究. 爆炸与冲击, 1997, 17(4): 369-374.
    [18] 师华强宗智贾敬蓓 . 水下爆炸冲击波的近场特性. 爆炸与冲击, 2009, 29(2): 125-130. doi: 10.11883/1001-1455(2009)02-0125-06
    [19] 李梅蒋建伟王昕 . 复合装药空气中爆炸冲击波传播特性. 爆炸与冲击, 2018, 38(2): 367-372. doi: 10.11883/bzycj-2016-0209
    [20] 王海福冯顺山 . 爆炸载荷下聚氨酯泡沫材料中冲击波压力特性. 爆炸与冲击, 1999, 19(1): 78-83.
  • 加载中
图(9)表(3)
计量
  • 文章访问数:  1904
  • HTML全文浏览量:  721
  • PDF下载量:  33
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-22
  • 录用日期:  2019-05-08
  • 网络出版日期:  2020-01-07
  • 刊出日期:  2020-01-01

装药动爆冲击波特性研究

    作者简介:陈龙明(1996- ),男,硕士研究生,chenlongming14@nudt.edu.cn
    通讯作者: 李志斌, lizhibin@nudt.edu.cn
  • 国防科技大学文理学院,湖南 长沙 410073

摘要: 导弹、炮弹等战斗部爆炸时具有一定的速度,较大的运动速度会使爆炸冲击波场分布发生变化,进而对弹药的毁伤威力产生影响。本文中采用AUTODYN软件对速度分别为0、272、340、680、1 020和1 700 m/s的TNT球形裸装药在空气中爆炸的冲击波场进行了仿真计算,定量研究装药在动爆条件下的峰值超压、比冲量和正压作用时间等威力参数特性。结果表明,方位角小于90°时装药速度与冲击波超压、比冲量成正相关,与正压作用时间成负相关;方位角大于90°时装药速度与冲击波超压、比冲量成负相关,与正压作用时间成正相关。超压峰值大小沿方位角成正弦变化。最后,分析了冲击波峰值超压数据,建立了动爆冲击波超压的计算模型,该模型计算结果与仿真和实验结果吻合较好。

English Abstract

  • 从一般的火炮到集各种先进技术于一身的导弹,大多数现代武器的战斗部速度均大于音速。战斗部在高速运动中爆炸所带来的各种特殊效应与战斗部静爆时相比有着特殊之处[1]。因此,研究战斗部动爆条件下的爆炸冲击波特性,对于战斗部的威力评价、武器毁伤效能评估以及动爆试验测试方案设计均有重要意义[2]

    冲击波静爆超压计算公式大多是建立在爆炸相似律的基础上,根据实测的多组数据拟合出来的经验公式,进行修正得到[3]。比较著名的经验公式有Henrych公式、萨多夫斯基公式、Brode公式等[4]。在动爆超压计算模型方面,现有的动爆理论中对于运动装药的处理,一般采取等效药量的方法。从能量相似原理出发,可以将运动TNT装药导致的威力增加与静止TNT装药量的增加相等效。根据等效药量公式结合上述静爆超压经验公式,可得到相应的动爆超压计算公式[5];杜红棉等[1]得到了运动装药在空中爆炸时形成的冲击波正方向和反方向的超压计算公式。蒋海燕等[6]利用AUTODYN软件对装药动爆冲击波场进行仿真,分析了动爆冲击波场的分布规律,并对数据进行回归分析,建立了与动爆试验结果符合较好的工程计算模型;聂源等[7]采用高精度显示欧拉流体力学软件SPEED,模拟了球形装药在空气介质中的爆炸过程,得到了相似的结果。

    本文中,通过AUTODYN软件进行建模计算,研究装药速度对爆炸冲击波流场演化的影响。分析研究动爆冲击波的冲击波超压、正压作用时间、比冲量的特性。并且从物理原理出发,利用模拟数据进行拟合分析,研究动爆冲击波对角度的依赖关系。最后建立修正因子δ的表达式,据此出发可以得到动爆冲击波的超压计算公式,与已有的结果相比应具备更好的准确度。

    • 为准确分析爆炸瞬时装药运动速度对冲击波威力场的影响规律,选用TNT球形裸装药进行计算。采用中心起爆,装药质量W=0.170 1 kg,装药密度ρ=1.63 g/cm3。采用二维轴对称模型,计算空气域尺寸为2 400 mm×1 200 mm,网格尺寸为2 mm×2 mm。采用多物质欧拉算法,材料直接从AUTODYN材料库中选取,空气采用理想气体状态方程,炸药采用JWL状态方程。在空气域的x轴中间填充质量为0.170 1 kg的球形TNT炸药。边界条件设定为压力流出边界条件,以模拟无限空气域。如图1所示,将TNT炸药中心定义为坐标原点,各监测点成环形围绕着坐标原点排列于空气域中,距爆心300 mm起间隔200 mm放置4个监测点,每列监测点之间夹角为30°,共7列28个监测点。定义方位角θ为爆心连线与装药运动方向的夹角,即爆心连线与x轴正方向的夹角。

      图  1  参数设置

      Figure 1.  Parameter setting

      为能较好地得到贴合实际的动爆冲击波特性,分别计算v为0、272、340、680、1 020和1 700 m/s等6种速度条件下的冲击波场,分别对应0Ma、0.8Ma、1.0Ma、2.0Ma、3.0Ma和5.0Ma的情况。

    • 分析装药动爆冲击波流场演化图(见图2),可以发现:不同运动速度的装药在无限空中爆炸形成的冲击波均以球面波形式扩展;当v=0 m/s(静爆)时,压力分布呈现规则的球形,以x=0为对称轴呈对称结构,超压峰值约为0.37 MPa,因此在静爆条件下球形裸装药的爆炸冲击波场是规则的球对称结构;当v=272,340,680,1 020 m/s(动爆)时,冲击波压力分布随着装药运动速度的增大而严重畸变,与速度方向夹角越小的地方压力峰值越大,相同运动速度时的冲击波超压峰值出现在装药运动速度方向上。对比爆轰产物分布状态:当v=0 m/s时,TNT材料分布以x=0为对称轴呈左右对称状态,形成蝴蝶状;随着速度的增加,爆轰产物分布逐渐变形,形状变得不规则;当v=1 700 m/s时,超压峰值达到了0.96 MPa,是v=0 m/s(静爆)时超压峰值的2.6倍,此时的爆轰产物分布则不呈现蝴蝶状,形状不规则。

      图  2  当t=0.4 ms时不同速度下的冲击波压力云图

      Figure 2.  Pressure contour of shock waves under different velocities at t=0.4 ms

      为了验证数值模拟的有效性,将速度v=0 m/s、方位角θ=0°的数据与Henrych公式计算结果[8]进行比较(见表1)。数值计算的结果与经验公式相比较,误差都小于5%,这验证了数值计算结果的有效性和准确性。

      Δp/MPaη/%
      模拟 Henrych公式
      2.492.384.41
      0.970.943.78
      0.470.471.07
      0.280.272.22

      表 1  超压数据与公式计算结果比较

      Table 1.  Comparison of peak overpressure results with theoretical values

      本文中,主要研究起毁伤破坏作用的冲击波超压峰值Δp、冲击波正压作用时间τ、冲击波比冲量I等3个威力参数。从图3可以发现,当比例距离$\overline R $变大时,不同装药速度的超压越来越小,冲击波超压随$\overline R $增大衰减得较快。从超压的角分布来看:当θ=90°时,不同装药运动速度的动爆冲击波超压曲线汇聚于一点,数值大小与静爆时相近;当θ>90°时,超压曲线随装药速度大小由下至上排列,可以得到装药速度越大冲击波超压越小的结论;当θ<90°时,超压曲线随装药速度大小由上至下排列,可以得到装药速度越大冲击波超压越大的结论。因此可以得到结论:当90°<θ<180°、比例距离一定时,装药速度越大冲击波超压峰值越小;当0°<θ<90°、比例距离一定时,装药速度越大冲击波超压峰值越大。这个结果与冲击波流场演化图中观察到的结果一致。

      图  3  不同比例距离时超压峰值随方位角的变化

      Figure 3.  Peak overpressure Δp changes with θ under different $\overline R $

      为了描述对称分布的监测点测得的冲击波超压之间的大小关系,引入一个增大系数r,定义为以中轴线为对称轴对称分布、相同比例距离处的监测点测得的超压峰值之间的比。

      增大系数描述了装药运动速度导致θ<90°方向与θ>90°方向的冲击波超压峰值变化程度。由图4可以看出,当速度值小于2.0Ma时,运动导致的超压峰值增幅小于2倍,然而当速度大于2.0Ma时,增幅愈加明显。从2.0Ma开始,速度增幅分别为50%和150%,然而增大系数的增幅达到78%和418%,超压增幅程度远大于速度增幅程度。后两个对称方向的增大系数均小于θ=0°方向与θ=180°方向,说明装药速度一定时,超压峰值沿θ=0°至θ=180°呈现一个由大到小的分布状态。

      图  4  不同方向的增大系数

      Figure 4.  Increasing ratios r in different directions

      图5所示,比冲量也明显呈现出随方向角分布的特点:相同比例距离下,比冲量随角度的增大而减小;大于90°,装药速度越大比冲量越大;小于90°,装药速度越小比冲量越小。

      图  5  不同比例距离时比冲量随方位角的变化

      Figure 5.  Specific impluse I changes with θ under different $\overline R $

      总体上,有装药速度越大、正压作用时间越大的趋势。由图6可见:当θ<90°时,有装药速度越小、正压作用时间越大的趋势;当θ>90°时,有装药速度越大、正压作用时间越大的趋势。此外,当比例距离较小时,正压作用时间的变化规律不是很明显,可以知道距离爆炸中心近区的物理参数变化情况比较复杂,变化的规律性不强。

      图  6  不同比例距离时正压时间随方位角的变化

      Figure 6.  Position pressure time τ changes with θ under different $\overline R $

    • 对动爆冲击波场进行简单分析,刚形成的冲击波可以近似为以爆炸中心为球心的球形冲击波。球面上的冲击波速度可以通过静爆冲击波速度${{{v}}_{\rm s}}$与装药运动速度${{v}}_0$的矢量相加,近似求出该处的动爆冲击波速度:${{{v}}_{\rm s}} = {{{v}}_{{\rm s} }}+{{{v}}_0}$。显然,随着方位角增大(0°<θ<180°),冲击波速度逐渐减小,因此冲击波超压大小变化与角度相关。这与前文的仿真数据分析得到的结论一致。

      在描述动爆超压时可以引入修正因子δ[7],定义δ为:

      $\delta = \frac{{\Delta {p_{\rm m}}}}{{\Delta {p_0}}}$

      修正因子δ是与装药运动速度v、方位角θ和比例距离$\overline R $等3个因子相关的量,即:

      $\delta = f(v,\overline R,\theta )$

      式中:f (v)为在特定的方位角θ*和特定的比例距离$\overline R ^*$处(θ*$\overline R ^*$均为任取)、不同装药速度v时的δf ($\overline R $)为在特定的方位角θ*、取一系列装药速度v时、不同比例距离$\overline R ^*$处的δ$\overline R $ = $\overline R ^*$时的δ之比;f (θ)为取一系列装药速度v和比例距离$\overline R $时、不同方位角θ处的δθ=θ*时的δ之比。

      θ*=0°、$\overline R ^*$ = 0.903 m/kg1/3,得到不同装药运动速度下的修正因子δ,可拟合得到该平均值随装药运动速度v的关系曲线(见图7),即函数f(v)为:

      图  7  拟合曲线f (v)

      Figure 7.  Fitting curve f (v)

      $f(v) = 1 + 0.27\frac{v}{{{c_0}}} + 0.08{\left( {\frac{v}{{{c_0}}}} \right)^2}$

      式中:装药速度v的单位为m/s,c0为声速,取340 m/s。

      θ* = 0°,得到一系列不同装药运动速度v、不同比例距离$\overline R $处的修正因子δ。在各种不同装药运动速度下,对不同比例距离$\overline R $时的修正因子δ取平均值,可拟合得到该平均值随比例距离$\overline R $的关系曲线(见图8),即函数f($\overline R $)为:

      图  8  拟合曲线f ($\overline R $)

      Figure 8.  Fitting curve f ($\overline R $)

      $f(\overline R) = 0.75 + 0.47 \overline R - 0.21 {\overline R^2}$

      式中:比例距离$\overline R $的单位为m/kg1/3

      取一系列不同装药运动速度v和比例距离$\overline R $时,画出不同比例距离$\overline R $处的修正因子δ与方位角θ之间的关系曲线。在各种不同装药运动速度和比例距离下,对不同方位角θ时的修正因子δ取平均值,可拟合得到该平均值随方位角θ的关系曲线(见图9),即函数f (θ)为:

      图  9  拟合曲线f (θ)

      Figure 9.  Fitting curves f (θ)

      $f(\theta ) = 0.73 + 0.29\sin \left( {\frac{{\theta + 1.61}}{{2.95}}{{\text{π}}} } \right)$

      式中:方位角θ的单位为rad。

      得到修正因子δ的表达式为:

      $\delta = \left[ {1 + 0.27\frac{v}{{{c_0}}} + 0.08{{\left( {\frac{v}{{{c_0}}}} \right)}^2}} \right] {\left( {0.75 + 0.47\overline R - 0.21{{\overline R}^2}} \right)}\left[ {0.73 + 0.29\sin \left( {\frac{{\theta + 1.61}}{{2.95}}{{\text{π}}}} \right)} \right]$

      式中:修正因子的适用范围为$0.5\;{\rm{m}}/{\rm{k}}{{\rm{g}}^{1/3}} {\text{≤}} \overline R {\text{≤}} 1.5\;{\rm{m}}/{\rm{k}}{{\rm{g}}^{1/3}}$。将δ与各种经验公式组合便能得到不同公式的动爆冲击波计算模型。

    • 为了验证更高的可靠度,仍选用TNT球形裸装药进行计算,TNT装药质量W=0.170 1 kg,装药密度ρ=1.63 g/cm3。计算得到装药运动速度v为100、500、1 500 m/s时的动爆超压场的压力分布。分别选取比例距离$\overline R $为0.72、1.81、1.26、0.54 m/kg1/3,方位角为θ为15°、73°、164°处的10组超压数据,见表2

      v/(m·s−1)$\overline R $/(m·kg−1/3)θ/(°)Δp/kPaη/%
      公式静爆模拟
      1000.72 151 609.031 516.351 663.803.29
      1001.81 15214.91219.84236.429.09
      1001.26 73394.31480.98493.1320.04
      1000.541641 155.982 552.412 322.7950.23
      5000.72 152 324.811 516.352 367.871.82
      5001.81 15310.52219.84317.832.30
      5001.26 73569.71480.98533.516.79
      5000.541641 670.222 552.411 579.855.72
      1 5000.72 155 551.671 516.355 249.825.75
      1 5001.81 15741.52219.84690.907.33

      表 2  仿真计算结果对比

      Table 2.  Comparison of simulation results

      采用运动装药爆炸冲击波的测试数据[9],对公式进行验证,见表3。实验所采用的是质量为0.17 kg的B炸药,其与TNT当量换算系数取1.3。

      v/(m·s−1)$\overline R $/(m·kg−1/3)θ/(°)Δp/kPaη/%
      实验公式
      534.311.37 15627.42638.32 1.71
      534.311.37 45565.37572.21 1.20
      534.311.37105361.97375.33 3.56
      587.651.19 15751.53926.0318.84
      587.651.19 45599.84830.1327.74
      587.651.19105379.21544.5030.36

      表 3  实验结果对比

      Table 3.  Comparison of experimental results

      静爆冲击波由Henrych公式计算[8]

      $\Delta {p_{\rm m}} = \left\{ \begin{aligned} &\frac{{{\rm{1}}{\rm{.407\;17}}}}{{\overline R }} + \frac{{0.553\;97}}{{{{\overline R }^2}}} + \frac{{0.035\;72}}{{{{\overline R }^3}}} + \frac{{0.000\;625}}{{{{\overline R }^4}}}\quad\quad\quad 0.05 {\text{≤}} \overline R {\text{≤}} 0.3\\ &\frac{{{\rm{0}}{\rm{.619\;38}}}}{{\overline R }} - \frac{{0.032\;62}}{{{{\overline R }^2}}} + \frac{{0.213\;24}}{{{{\overline R }^3}}}\quad\quad\quad\quad\quad\quad\quad\quad 0.3 {\text{<}} \overline R {\text{≤}} 1.0\\ &\frac{{{\rm{0}}{\rm{.066\;2}}}}{{\overline R }} + \frac{{0.405}}{{{{\overline R }^2}}} + \frac{{0.328\;8}}{{{{\overline R }^3}}}\quad\quad\quad\quad\;\,\quad\quad\quad\quad\quad\quad 1.0 {\text{<}} \overline R {\text{≤}} 10 \end{aligned} \right.$

      式中:Δpm的单位为MPa,$ {\overline R } $的单位为m/kg1/3

      通过误差分析,得到比例距离、角度、装药速度不全相同的10组数据的公式验证,其中8组的误差小于10%,模型计算结果基本达到预期效果。在与实验数据的检验中,误差的平均值小于14%,考虑到当时实验器材精度状况以及实验误差,这个结果符合得比较好。

    • 通过使用AUTODYN软件进行仿真计算,分析动爆冲击波场的变化,研究了冲击波超压峰值、冲击波正压作用时间、冲击波比冲量3个重要参量,得出了相应的变化规律。主要有:装药速度对小于90°的方向的冲击波超压、比冲量起增强的作用,对90°到180°的方向的冲击波超压、比冲量起减弱的作用;而装药速度对小于90°的方向的冲击波正压作用时间起减小作用,对于大于90°的方向的正压作用时间起增大作用;并且基于仿真结果进行分析,得到了动爆冲击波的一个工程计算模型。将本文得出的计算模型与软件仿真算出的不同装药速度下的结果作对比。读取不同比例距离、角度处测得的冲击波峰值超压,并将数据与根据模型计算得出的动爆超压峰值以及与实验测试数据对比分析其误差,结果表明计算模型能得到较好的结果,因此本次研究得出的计算模型具备一定的应用价值。

参考文献 (9)

目录

    /

    返回文章
    返回