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

圆柱形弹药空气中爆炸相似性规律

陈材 石全 尤志锋 郭驰名 戈洪宇

引用本文:
Citation:

圆柱形弹药空气中爆炸相似性规律

    作者简介: 陈 材(1990- ),男,博士研究生,caichen20091165@126.com;
  • 中图分类号: O381; TP391.9

Similarity law of cylindrical ammunition explosions in air

  • CLC number: O381; TP391.9

  • 摘要: 针对战场损伤试验中弹药费用高、使用风险系数大的问题,采用量纲分析的方法对圆柱形弹药相似性模型的建立方法进行分析,研究原型同相似性模型之间冲击波峰值压力的关系,利用有限元仿真软件AUTODYN对冲击波传播过程进行数值仿真,并结合实弹试验数据对相似模型的有效性进行验证。研究结果表明,AUTODYN仿真结果同量纲分析结果一致,同时与实弹试验结果也吻合较好,说明在实际试验中可以使用缩比模型替代原型进行试验。该研究为弹药缩比模型在战场损伤试验中的使用提供了理论基础,具有一定的工程价值。
  • 图 1  数值仿真模型

    Figure 1.  Numerical simulation model

    图 2  不同时刻弹药爆炸冲击波场

    Figure 2.  Shock wave field of ammunition explosionat different moments

    图 3  轴向冲击波峰值压力

    Figure 3.  Peak pressure of shock wave in the axial direction

    图 4  径向冲击波峰值压力

    Figure 4.  Peak pressure of shock wave in the radial direction

    图 5  弹药战斗部实体

    Figure 5.  Entity of ammunition warhead

    图 6  冲击波测量试验布置图

    Figure 6.  Test layout of shock wave measurement

    表 1  弹药空气中爆炸相关参数及其量纲

    Table 1.  Explosion-related parameters and their dimensions

    物理量符号量纲物理量符号量纲
    峰值压力pmML−1T−2空气初始状态压力p0ML−1T−2
    装药质量QM空气密度ρaML−3
    装药密度ρeML−3空气绝热指数γa1
    质量化学能eL2T−2弹壳密度ρsML−3
    膨胀指数γe1弹壳弹性模量EsML−1T−2
    长径比k1弹壳泊松比μs1
    爆距RL弹壳屈服极限YsML−1T−2
    赤道面夹角α1弹壳厚度hsL
    轴向夹角φ1
    下载: 导出CSV

    表 2  弹药原型及缩比模型结构参数值

    Table 2.  Ammunition prototype and scale model structural parameter values

    弹药模型TNT质量/kg长度/mm半径/mm外壳厚度/mm
    原型64.0004801606.00
    1/2缩比模型 8.000240 803.00
    1/3缩比模型 2.370162 542.00
    1/4缩比模型 1.000120 401.25
    1/5缩比模型 0.512 96 321.20
    下载: 导出CSV

    表 3  缩比模型中同原型具有相同冲击波峰值压力的爆距

    Table 3.  Detonation distance with same peak pressure of original model in the similar model

    弹药类型爆距/m
    原型1.0002.0003.0004.0005.0006.0007.000
    1/2缩比模型0.5001.0001.5002.0002.5003.0003.500
    1/3缩比模型0.3330.6671.0001.3331.6672.0002.333
    1/4缩比模型0.2500.5000.7501.0001.2501.5001.750
    1/5缩比模型0.2000.4000.6000.8001.0001.2001.400
    下载: 导出CSV

    表 4  轴向方向上冲击波峰值压力

    Table 4.  Peak pressure of shock wave in the axial direction

    弹药类型峰值压力/kPa
    原型12 744.14 763.61 894.9911.9552.1376.7277.7
    1/2缩比模型12 719.54 761.21 892.5907.7540.0366.3273.6
    1/3缩比模型12 764.24 819.81 922.4929.8560.6382.6280.3
    1/4缩比模型12 739.84 752.31 897.0910.3551.0374.6275.7
    1/5缩比模型12 709.54 714.01 891.2910.6543.1380.5279.4
    下载: 导出CSV

    表 5  径向方向上冲击波峰值压力

    Table 5.  Peak pressure of shock wave in the radial direction

    弹药类型峰值压力/kPa
    原型 9 927.24 357.52 160.91 145.5692.0466.4349.6
    1/2缩比模型 9 925.54 354.02 161.21 144.3691.9464.4345.4
    1/3缩比模型10 099.04 465.62 198.91 171.7701.6474.8353.6
    1/4缩比模型 9 921.24 356.22 158.81 138.9690.0463.3347.3
    1/5缩比模型 9 916.94 355.22 157.21 134.9685.9465.3343.5
    下载: 导出CSV

    表 6  缩比模型轴向方向峰值压力绝对值偏差

    Table 6.  Absolute deviation of peak pressure of scale model in the axial direction

    弹药类型缩比模型轴向峰值压力同原型偏差/%
    1/2缩比模型0.200.050.130.462.192.761.46
    1/3缩比模型0.161.181.451.961.551.582.77
    1/4缩比模型0.040.240.110.180.190.550.72
    1/5缩比模型0.271.040.210.141.631.020.63
    下载: 导出CSV

    表 7  缩比模型径向方向峰值压力绝对值偏差

    Table 7.  Absolute deviation of peak pressure of scale model in the radial direction

    弹药类型缩比模型径向峰值压力偏差/%
    1/2缩比模型0.020.080.010.100.020.431.19
    1/3缩比模型1.732.481.762.291.391.801.16
    1/4缩比模型0.060.030.100.580.290.670.66
    1/5缩比模型0.100.050.170.930.890.241.73
    下载: 导出CSV

    表 8  冲击波峰值压力测量距离

    Table 8.  Measure distance of peak pressure

    弹药类型测量距离/m
    原型14.0019.0024.00
    缩比模型A111.6115.7519.90
    缩比模型A2 7.31 9.9212.53
    缩比模型A3 4.61 6.25 7.90
    下载: 导出CSV

    表 9  现场试验与仿真试验冲击波峰值压力

    Table 9.  Peak pressure in field and simulation tests

    试验类型冲击波峰值压力/kPa
    原型实际试验257.3185.5145.3
    缩比模型A1仿真试验248.2176.5137.2
    缩比模型A2仿真试验246.0176.0137.2
    缩比模型A3仿真试验247.5177.0136.7
    下载: 导出CSV

    表 10  冲击波峰值压力误差

    Table 10.  Error of peak pressure

    试验
    类型
    误差/%
    缩比模型A1仿真试验3.54.95.6
    缩比模型A2仿真试验4.45.15.6
    缩比模型A3仿真试验3.84.65.9
    下载: 导出CSV
  • [1] 吴晓颖, 李帆, 张万君, 等. 装备毁伤模拟试验方案的设计与优化 [J]. 兵器装备工程学报, 2014, 35(10): 5–7. DOI: 10.11809/scbgxb2014.10.002.
    WU Xiaoying, LI Fan, ZHANG Wanjun, et al. Design and optimization of equipment damage simulation test [J]. Journal of Ordnance Equipment Engineering, 2014, 35(10): 5–7. DOI: 10.11809/scbgxb2014.10.002.
    [2] 石全. 装备战伤理论与技术 [M]. 北京: 国防工业出版社, 2007: 243−244.
    [3] 杨亚东, 李向东, 王晓鸣, 等. 密闭空间内爆炸缩比相似模型研究 [J]. 振动与冲击, 2014, 33(2): 128–133. DOI: 10.3969/j.issn.1000-3835.2014.02.025.
    YANG Yadong, LI Xiangdong, WANG Xiaoming, et al. Scale similarity model of internal explosion in closed field [J]. Journal of Vibration and Shock, 2014, 33(2): 128–133. DOI: 10.3969/j.issn.1000-3835.2014.02.025.
    [4] 杨亚东, 李向东, 王晓鸣, 等. 钢筋混凝土结构内爆炸相似模型试验研究 [J]. 南京理工大学学报, 2016, 40(2): 135–141.
    YANG Yadong, LI Xiangdong, WANG Xiaoming, et al. Experimental study on similarity model of reinforce concrete structure under internal explosion [J]. Journal of Nanjing University of Science and Technology, 2016, 40(2): 135–141.
    [5] 吕祥锋, 潘一山, 李忠华, 等. 冲击波作用下巷道破坏规律相似模拟研究 [J]. 振动与冲击, 2011, 30(10): 212–215. DOI: 10.3969/j.issn.1000-3835.2011.10.041.
    LV Xiangfeng, PAN Yishan, LI Zhonghua, et al. Similarity simulation tests for failure law of a rock burst roadway under shock wave [J]. Journal of Vibration and Shock, 2011, 30(10): 212–215. DOI: 10.3969/j.issn.1000-3835.2011.10.041.
    [6] 张玉磊, 王胜强, 袁建飞, 等. 不同量级TNT爆炸冲击波参数相似律实验研究 [J]. 弹箭与制导学报, 2016, 36(6): 53–56.
    ZHANG Yulei, WANG Shengqiang, YUAN Jianfei, et al. Experimental research on similarity law of explosive shock wave parameters with different orders of magnitude TNT [J]. Journal of Projectiles, Rocks, Missiles and Guidance, 2016, 36(6): 53–56.
    [7] 荣吉利, 何轩, 项大林, 等. 定向战斗部水下爆炸相似性研究 [J]. 兵工学报, 2014, 35(9): 1329–1334. DOI: 10.3969/j.issn.1000-1093.2014.09.001.
    RONG Jili, HE Xuan, XIANG Dalin, et al. The similarity research on underwater explosion of directional warhead [J]. Acta Armamentarii, 2014, 35(9): 1329–1334. DOI: 10.3969/j.issn.1000-1093.2014.09.001.
    [8] 刘文韬, 姚熊亮, 李帅, 等. 离心机水下爆炸缩比实验原理及数值研究 [J]. 爆炸与冲击, 2016, 36(6): 789–796. DOI: 10.11883/1001-1455(2016)06-0789-08.
    LIU Wentao, YAO Xiongliang, LI Shuai, et al. Experimental principle and numerical study of scaled-down underwater explosion model on a centrifuge apparatus [J]. Explosion and Shock Waves, 2016, 36(6): 789–796. DOI: 10.11883/1001-1455(2016)06-0789-08.
    [9] LI Jian, RONG Jili. Experimental and numerical investigation of the dynamic response of structures subjected to underwater explosion [J]. European Journal of Mechanics, 2012, 32(1): 59–69. DOI: 10.1016/j.euromechflu.2011.09.009.
    [10] 朱如意. 杀爆榴弹动静态爆炸仿真及毁伤威力差异性分析[D]. 石家庄: 军械工程学院, 2016: 31−34.
  • [1] 李雨成刘天奇周西华 . 基于量纲分析理论的煤尘爆炸能量预测模型. 爆炸与冲击, 2017, 37(3): 566-570. doi: 10.11883/1001-1455(2017)03-0566-05
    [2] 庞伟宾李永池何翔 . 化爆冲击波在T型通道内到时规律的实验研究. 爆炸与冲击, 2007, 27(1): 63-67. doi: 10.11883/1001-1455(2007)01-0063-05
    [3] 苏子星何继业 . 基于Cowper-Symonds方程的相似理论修正方法. 爆炸与冲击, 2018, 38(3): 654-658. doi: 10.11883/bzycj-2016-0308
    [4] 周文海梁瑞余建平杜超飞王敦繁楼晓明 . 边坡抛掷爆破峰值质点振动速度的无量纲分析. 爆炸与冲击, 2019, 39(5): 054201-1-054201-8. doi: 10.11883/bzycj-2017-0373
    [5] 宋梅利李文彬王晓鸣冯君刘志林 . 弹体高速侵彻效率的实验和量纲分析. 爆炸与冲击, 2016, 36(6): 752-758. doi: 10.11883/1001-1455(2016)06-0752-07
    [6] 周杰陶钢王健 . 爆炸冲击波对肺损伤的数值模拟. 爆炸与冲击, 2012, 32(4): 418-422. doi: 10.11883/1001-1455(2012)04-0418-05
    [7] 顾文彬郑向平刘建青李丹俊陆鸣 . 浅层水中爆炸冲击波对混凝土墩斜碰撞作用试验研究. 爆炸与冲击, 2006, 26(4): 361-366. doi: 10.11883/1001-1455(2006)04-0361-06
    [8] 陈文张庆明 , . 爆炸冲击波下导弹结构动力分析模型初步研究. 爆炸与冲击, 2009, 29(2): 199-204. doi: 10.11883/1001-1455(2009)02-0199-06
    [9] 彭常贤 , . 声振耦合分析在冲击波测量中的应用. 爆炸与冲击, 2008, 28(5): 427-432. doi: 10.11883/1001-1455(2008)05-0427-06
    [10] 沈雁鸣陈坚强 . 超高速碰撞相似律的数值模拟验证. 爆炸与冲击, 2011, 31(4): 343-348. doi: 10.11883/1001-1455(2011)04-0343-06
    [11] 赖华伟王占江杨黎明王礼立 . 线性黏弹性球面波的特征线分析. 爆炸与冲击, 2013, 33(1): 1-10. doi: 10.11883/1001-1455(2013)01-0001-10
    [12] 钟冬望何理操鹏张奎 . 爆破振动持时分析及微差爆破延期时间优选. 爆炸与冲击, 2016, 36(5): 703-709. doi: 10.11883/1001-1455(2016)05-0703-07
    [13] 卢强王占江王礼立赖华伟杨黎明 , . 基于ZWT方程的线黏弹性球面波分析. 爆炸与冲击, 2013, 33(5): 463-470. doi: 10.11883/1001-1455(2013)05-0463-08
    [14] 姚成宝李若田宙郭永辉 . 空气自由场中强爆炸冲击波传播二维数值模拟. 爆炸与冲击, 2015, 35(4): 585-590. doi: 10.11883/1001-1455(2015)04-0585-06
    [15] 张柱晋艳娟 . 反向起爆模型下的冲击波加载. 爆炸与冲击, 2014, 34(2): 223-228. doi: 10.11883/1001-1455(2014)02-0223-06
    [16] 师华强宗智贾敬蓓 . 水下爆炸冲击波的近场特性. 爆炸与冲击, 2009, 29(2): 125-130. doi: 10.11883/1001-1455(2009)02-0125-06
    [17] 周佩杰王坚陶钢周杰 . 泡沫材料对冲击波的衰减特性. 爆炸与冲击, 2015, 35(5): 675-681. doi: 10.11883/1001-1455(2015)05-0675-07
    [18] 刘贵兵侯海量朱锡张国栋 . 液滴对爆炸冲击波的衰减作用. 爆炸与冲击, 2017, 37(5): 844-852. doi: 10.11883/1001-1455(2017)05-0844-09
    [19] 李梅蒋建伟王昕 . 复合装药空气中爆炸冲击波传播特性. 爆炸与冲击, 2018, 38(2): 367-372. doi: 10.11883/bzycj-2016-0209
    [20] 姜涛由文立张可玉丁浩 . 水中爆炸表面空穴的理论研究. 爆炸与冲击, 2011, 31(1): 19-24. doi: 10.11883/1001-1455(2011)01-0019-06
  • 加载中
图(7)表(10)
计量
  • 文章访问数:  135
  • HTML全文浏览量:  607
  • PDF下载量:  33
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-07-14
  • 录用日期:  2018-08-25
  • 刊出日期:  2019-09-01

圆柱形弹药空气中爆炸相似性规律

    作者简介:陈 材(1990- ),男,博士研究生,caichen20091165@126.com
  • 1. 陆军工程大学,河北 石家庄 050003
  • 2. 中国白城兵器试验中心,吉林 白城 137001

摘要: 针对战场损伤试验中弹药费用高、使用风险系数大的问题,采用量纲分析的方法对圆柱形弹药相似性模型的建立方法进行分析,研究原型同相似性模型之间冲击波峰值压力的关系,利用有限元仿真软件AUTODYN对冲击波传播过程进行数值仿真,并结合实弹试验数据对相似模型的有效性进行验证。研究结果表明,AUTODYN仿真结果同量纲分析结果一致,同时与实弹试验结果也吻合较好,说明在实际试验中可以使用缩比模型替代原型进行试验。该研究为弹药缩比模型在战场损伤试验中的使用提供了理论基础,具有一定的工程价值。

English Abstract

  • 为打赢高新技术条件下的现代战争,提升部队实战化水平,要求部队的作战训练更贴近实战[1]。因此在战场损伤试验中,使用制式弹药对实际装备进行打击毁伤试验,成为了研究装备损伤规律、提高装备战场环境下生存力的重要途径[2]。然而对于制式弹药而言:一方面其造价高昂,造成试验成本过高;另一方面其毁伤威力大,造成试验风险系数过高,同时对装备的损伤程度不易把控,无法高效获取损伤数据。这两方面因素共同造成了目前装备战场损伤试验费效比过高、无法广泛开展的问题。因此,研究弹药空气中爆炸相似性问题,探寻制式弹药替代物的选择与建立方法,具有十分重要的现实意义。

    目前,采用相似理论对炸药空气中爆炸进行研究,目的多为研究冲击波的毁伤机理及优化其理论公式。杨亚东等[3-4]以相似理论为基础,针对毁伤对象研究了利用相似模型对原模型毁伤规律进行分析的方法,吕祥锋等[5]通过建立相似模拟试验模型,研究了在不同装药量下巷道破坏规律,张玉磊等[6]研究了不同量级TNT裸装药在空气中爆炸的冲击波参数相似规律。也有学者对炸药在水下爆炸的相似性规律进行了研究:荣吉利等[7]和刘文韬等[8]利用数值仿真方法,研究了水下爆炸试验中缩比模型的建立方法及冲击波的毁伤规律。

    本文以相似理论为理论基础,采用量纲分析的方法,研究圆柱形弹药战斗部的相似性规律,分析不同缩比模型同原型之间冲击波峰值压力的关系,并利用有限元仿真软件AUTODYN进行仿真计算,在此基础上,使用实弹试验所采集数据同缩比模型仿真数据进行对比分析,以此验证相似性模型的正确性。

    • 弹药在空气中爆炸时,所产生的冲击波强度主要受到装药、空气和弹药壳体性质的影响。因此冲击波的强度特征主要受到以下几个方面的控制因素影响。

      (1)装药相关参数:装药质量Q;单位质量装药所释放化学能e;装药密度ρe;爆炸产物膨胀指数γe;弹药长径比k

      (2)距弹药几何中心距离R;测试点同弹药几何中心连线同弹药赤道面夹角α;测试点同弹药几何中心连线同弹药轴向夹角φ

      (3)空气相关参数:初始状态压力p0;空气密度ρa;空气绝热指数γa。忽略空气温度、传热性和空气黏性等次要因素。

      (4)弹药壳体相关参数:密度ρs;弹性模量Es;泊松比μs;屈服极限Ys;厚度hs

    • 以时间T,质量M,长度L为基本量纲,各个物理量量纲见表1

      物理量符号量纲物理量符号量纲
      峰值压力pmML−1T−2空气初始状态压力p0ML−1T−2
      装药质量QM空气密度ρaML−3
      装药密度ρeML−3空气绝热指数γa1
      质量化学能eL2T−2弹壳密度ρsML−3
      膨胀指数γe1弹壳弹性模量EsML−1T−2
      长径比k1弹壳泊松比μs1
      爆距RL弹壳屈服极限YsML−1T−2
      赤道面夹角α1弹壳厚度hsL
      轴向夹角φ1

      表 1  弹药空气中爆炸相关参数及其量纲

      Table 1.  Explosion-related parameters and their dimensions

      因此,由相似第二定理(Π定理)可知,弹药在空气中爆炸的峰值压力pm与以上控制参数的关系为:

      $ {p_{\rm{m}}} = f(Q,{\rho _{\rm{e}}},e,{\gamma _{\rm{e}}},k,R,\alpha ,\phi ,{p_0},{\rho _{\rm{a}}},{\gamma _{\rm{a}}},{Y_{\rm{s}}},{\mu _{\rm{s}}},{h_{\rm{s}}},{\rho _{\rm{s}}},{E_{\rm{s}}}) $

      Qρee作为基本量,则式(1)可转化为以下的无量纲形式:

      $ \frac{{{p_{\rm{m}}}}}{{{\rho _{\rm{e}}}e}} = f\left( {{\gamma _{\rm{e}}},k,\frac{R}{{{{(Q \cdot \rho _{\rm{e}}^{ - 1})}^{1/3}}}},\alpha ,\phi ,\frac{{{p_0}}}{{{\rho _{\rm{e}}}e}},\frac{{{\rho _{\rm{a}}}}}{{{\rho _{\rm{e}}}}},{\gamma _{\rm{a}}},{\kern 1pt} \frac{{{Y_{\rm{s}}}}}{{{\rho _{\rm{e}}}e}},{\mu _{\rm{s}}},{\kern 1pt} \frac{{{h_{\rm{s}}}}}{{{{(Q \cdot \rho _{\rm{e}}^{ - 1})}^{1/3}}}},\frac{{{\rho _{\rm{s}}}}}{{{\rho _{\rm{e}}}}},\frac{{{E_{\rm{s}}}}}{{{\rho _{\rm{e}}}e}}} \right) $

      使用相同炸药在相同环境下做相似模型试验,则弹药峰值压力影响参数中的10个有关参数同原型保持一致,即:

      $ ({\rho _{\rm{e}}},e,{\gamma _{\rm{e}}},{p_0},{\rho _{\rm{a}}},{\gamma _{\rm{a}}},{\rho _{\rm{s}}},{E_{\rm{s}}},{\mu _{\rm{s}}},{Y_{\rm{s}}}) = {\rm{const}} $

      因此式(2)可简化为:

      $ \frac{{{p_{\rm{m}}}}}{{{\rho _{\rm{e}}}e}} = f\left( {k,\frac{R}{{{{(Q \cdot \rho _{\rm{e}}^{ - 1})}^{1/3}}}},\alpha ,\phi ,\frac{{{h_{\rm{s}}}}}{{{{(Q \cdot \rho _{\rm{e}}^{ - 1})}^{1/3}}}}} \right) $

      由相似定律可知,若要使模型同原型峰值压力相似,则需要满足以下条件:

      $ \left\{ \begin{array}{l} {k_{\rm{m}}} = {k_{\rm{p}}}\\ {\alpha _{\rm{m}}} = {\alpha _{\rm{p}}}\\ {\phi _{\rm{m}}} = {\phi _{\rm{p}}}\\ {R_{\rm{m}}}/{({Q_{\rm{m}}} \cdot \rho {_{\rm{e,m}}^{ - 1})^{1/3} }}= {R_{\rm{p}}}/{({Q_{\rm{p}}} \cdot \rho {{_{\rm{e,p}}^{ - 1})}^{1/3}}}\\ {h_{\rm{m}}}/{({Q_{\rm{m}}} \cdot \rho {_{\rm{e,m}}^{ - 1})}^{1/3}} = {h_{\rm{p}}}/{({Q_{\rm{p}}} \cdot \rho {{_{\rm{e,p}}^{ - 1}}})^{1/3}} \end{array} \right. $

      即:

      $ \left\{ \begin{array}{l} {k_{\rm{m}}} = {k_{\rm{p}}}\\ {\alpha _{\rm{m}}} = {\alpha _{\rm{p}}}\\ {\phi _{\rm{m}}} = {\phi _{\rm{p}}}\\ {R_{\rm{m}}}/{R_{\rm{p}}} = {\left( {{Q_{\rm{m}}}/{Q_{\rm{p}}}} \right)^{1/3}}\\ {h_{\rm{m}}}/{h_{\rm{p}}} = {\left( {{Q_{\rm{m}}}/{Q_{\rm{p}}}} \right)^{1/3}} \end{array} \right. $

      式中:下标m代表相似模型参数,下标p代表原型参数。

      设相似模型同原型的装药质量缩比比例Qm/Qp=θ3,对于相似模型而言,只需要满足缩比后的弹药长径比k不变,弹药壳体厚度hm=θhp,则当爆距Rm=θRp,赤道面夹角αmp,轴向夹角φm=φp时,其爆炸产生的冲击波峰值压力pm同原型中爆距为Rp、赤道面夹角为αp、轴向夹角为φp时所产生的冲击波峰值压力相同。

      通过以上分析可知,在缩比模型满足相似率的条件下,若其装药质量为原型质量的θ3倍,则其在θ倍爆距处所产生的冲击波峰值压力与原模型相同。以27 kg炸药为例,其在爆距为3 m处同1 kg炸药在爆距为1 m处所产生的冲击波峰值压力相同。

    • 采用显式有限元分析软件AUTODYN对带壳装药爆炸冲击波场的形成和传播过程进行数值模拟。炸药采用TNT装药,形状为圆柱形,由文献[9]可知,当圆柱形炸药的长径比为1.5时,其爆炸产生的峰值压力较大,故将装药长径比定为1.5,起爆方式采用中心单点起爆。因此建立60 kg的带壳装药,其半径为160 mm,长度为480 mm,外壳厚度为6 mm,网格单元尺寸为1 mm×1 mm,将其置于边长为20 m×10 m的长方形空气域当中,为了得到冲击波传播过程中峰值压力的变化规律,以装药中心为圆心,在半径为10 m的范围内,在轴向及径向每隔1 m设置一个观测点,以此测量冲击波参数。由于冲击波问题具有对称性特点,为提高运算效率,模型采用2D轴对称的方式进行构建,单位制定为mm-mg-ms。空气域选用Euler算法进行计算,网格单元尺寸为1 mm×1 mm,采用填充的方式在空气域中建立TNT装药模型。为实现无限空气域的仿真效果,在空气域模型边界节点处设置无反射边界,从而避免冲击波在边界处产生折射和反射的现象。建立好的数值仿真模型如图1所示。在此基础上,分别建立缩比比例为1/2、1/3、1/4和1/5的缩比模型,模型数据如表2所示。

      图  1  数值仿真模型

      Figure 1.  Numerical simulation model

      弹药模型TNT质量/kg长度/mm半径/mm外壳厚度/mm
      原型64.0004801606.00
      1/2缩比模型 8.000240 803.00
      1/3缩比模型 2.370162 542.00
      1/4缩比模型 1.000120 401.25
      1/5缩比模型 0.512 96 321.20

      表 2  弹药原型及缩比模型结构参数值

      Table 2.  Ammunition prototype and scale model structural parameter values

    • 弹药外壳材料选用V250钢,采用Shock状态方程和Steinberg-Guinan强度方程对其进行描述。

      Shock状态方程表达式为${u_{\rm{s}}} = c + {s_1}{u_{\rm{p}}} + {s_2}u_{\rm{p}}^{\rm{2}}$,在强冲击条件下,对于大多数材料而言,其表达式可以满足us=c+sup。其中需要设置的部分参数:密度,8.129×103 kg/m3;Grüneisen系数,1.6;c,3.98 km/s;s1,1.58;参考温度,291.15 K;比热,408 J/(kg·K)。

      Steinberg-Guinan强度方程其表达式为:

      $ G = {G_0}\left[ {1 + \left( {\frac{{G_p'}}{{{G_0}}}} \right)\frac{p}{{{\eta ^{1/3}}}} + \left( {\frac{{G_T'}}{{{G_0}}}} \right)\left( {T - 300} \right)} \right] $

      $ Y = {Y_0}\left[ {1 + \left( {\frac{{Y_p'}}{{{Y_0}}}} \right)\frac{p}{{{\eta ^{1/3}}}} + \left( {\frac{{G_T'}}{{{G_0}}}} \right)\left( {T - 300} \right)} \right]{\left( {1{\rm{ + }}\beta \varepsilon } \right)^n} $

      $ {Y_0}{\left( {1 + \beta \varepsilon } \right)^n} \text{≤} {Y_{\max }} $

      式中: G为剪切模量,G0为初始状态的剪切模量,$G_p' = {\rm{d}}G/{\rm{d}}p$p为压力,η为压缩比,$G_T' = {\rm{d}}G/{\rm{d}}T$T为温度,Y为屈服应力,Y0为初始状态的屈服应力,$Y_p' = {\rm{d}}Y/{\rm{d}}p$β为硬度常数,ε为等效塑性应变,n为硬度幂数,Ymax为最大屈服应力, T为温度(K)。具体参数值:G0,71.8 GPa;Y0,1.56 GPa;Ymax,GPa;β,2;n,0.5;dG/dp,1.479;dG/dT,−22.62 MPa/K,dY/dp,0.032 14;熔化温度,2 310 K。材料失效模型选择主应变失效。失效主应变取值为0.2,最大切应变取值0.9。

    • 采用理想气体状态方程对空气进行描述:

      $ {p_1} = (\gamma - 1){\rho _1}{e_{\rm{a}}} + {p_{{\rm{shift}}}} $

      式中:p1为压力;ρ1为初始密度;ea为比热力学能;γ为理想气体常数,pshift为压力偏移量。在空气模型中,取γ为1.4,ρ1为1.225 kg/m3ea则取为206.8 J/g。

      TNT采用JWL状态方程进行描述:

      $ {p_2} = A\left( {1 - \frac{\omega }{{{R_1}V}}} \right){{\rm{e}}^{ - {R_1}V}} + B\left( {1 - \frac{\omega }{{{R_2}V}}} \right){{\rm{e}}^{ - {R_2}V}} + \frac{{\omega E}}{V} $

      式中:p2为爆轰压力,V为相对体积,E为初始比内能,ωABR1R2均为表征炸药特性的常数。TNT参数:ρTNT,1.63 g/cm3A,374 GPa;B,3.75 GPa;R1,4.15;R2,0.09;ω,0.35。

    • 通过对弹药爆炸及冲击波传播过程进行数值计算,可知弹药原型及各个缩比模型冲击波传播规律基本相同,以1/2缩比模型为例,可得到其不同时刻冲击波场空间状态信息如图2所示。由图2可知,在爆炸初始时刻,弹药内部装药急剧反应产生高压,并迅速向四周扩散。由于弹药是圆柱形结构,因此冲击波场在传播过程中将出现局部超高压区,在爆炸初期,当冲击波场撕裂弹药外壳后,首先在轴线方向出现超高压区,随着冲击波的不断传播,其在径向方向产生叠加,从而产生局部超高压区。由此分析可知,弹药爆炸后冲击波峰值压力将出现在径向及轴线方向,因此重点对径向及轴线方向上的峰值压力进行测量及分析。

      图  2  不同时刻弹药爆炸冲击波场

      Figure 2.  Shock wave field of ammunition explosionat different moments

    • 测量弹药轴向及径向方向不同距离处的冲击波峰值压力。对于弹药原型,分别测量爆距为1.000、2.000、3.000、4.000、5.000、6.000和7.000 m处的冲击波峰值压力。由量纲分析结果可知,对于不同缩比比例η的缩比模型,其与原型具有相同冲击波峰值压力的爆距如表3所示,由此分别测量轴向及径向对应各点的冲击波峰值压力,结果分别如表45所示。

      弹药类型爆距/m
      原型1.0002.0003.0004.0005.0006.0007.000
      1/2缩比模型0.5001.0001.5002.0002.5003.0003.500
      1/3缩比模型0.3330.6671.0001.3331.6672.0002.333
      1/4缩比模型0.2500.5000.7501.0001.2501.5001.750
      1/5缩比模型0.2000.4000.6000.8001.0001.2001.400

      表 3  缩比模型中同原型具有相同冲击波峰值压力的爆距

      Table 3.  Detonation distance with same peak pressure of original model in the similar model

      弹药类型峰值压力/kPa
      原型12 744.14 763.61 894.9911.9552.1376.7277.7
      1/2缩比模型12 719.54 761.21 892.5907.7540.0366.3273.6
      1/3缩比模型12 764.24 819.81 922.4929.8560.6382.6280.3
      1/4缩比模型12 739.84 752.31 897.0910.3551.0374.6275.7
      1/5缩比模型12 709.54 714.01 891.2910.6543.1380.5279.4

      表 4  轴向方向上冲击波峰值压力

      Table 4.  Peak pressure of shock wave in the axial direction

      弹药类型峰值压力/kPa
      原型 9 927.24 357.52 160.91 145.5692.0466.4349.6
      1/2缩比模型 9 925.54 354.02 161.21 144.3691.9464.4345.4
      1/3缩比模型10 099.04 465.62 198.91 171.7701.6474.8353.6
      1/4缩比模型 9 921.24 356.22 158.81 138.9690.0463.3347.3
      1/5缩比模型 9 916.94 355.22 157.21 134.9685.9465.3343.5

      表 5  径向方向上冲击波峰值压力

      Table 5.  Peak pressure of shock wave in the radial direction

      表45分别可得不同缩比比例下轴向及径向方向上冲击波峰值压力变化趋势图,如图34所示。

      图  3  轴向冲击波峰值压力

      Figure 3.  Peak pressure of shock wave in the axial direction

      图  4  径向冲击波峰值压力

      Figure 4.  Peak pressure of shock wave in the radial direction

      图34可知,弹药爆炸后冲击波峰值压力随着传播距离的增大而降低,且初期降低速率较快,后期逐渐趋于平缓;缩比模型爆炸冲击波峰值压力在表3所示的相应爆距处,能够与原型峰值压力基本保持一致。但是,随着缩比比例的逐渐加大,冲击波在单位距离上的变化率迅速加大,由于传感器测量精度的限制,造成实际试验中测量成本加大、测量精度降低,因此在实际试验中应根据试验具体情况合理选择缩比比例。

      以原型为参考,计算轴向上缩比模型冲击波峰值压力同原型对应点峰值压力的绝对值偏差,计算结果如表6所示,计算径向上缩比模型冲击波峰值压力同原型对应点峰值压力的绝对值偏差,结果如表7所示。

      弹药类型缩比模型轴向峰值压力同原型偏差/%
      1/2缩比模型0.200.050.130.462.192.761.46
      1/3缩比模型0.161.181.451.961.551.582.77
      1/4缩比模型0.040.240.110.180.190.550.72
      1/5缩比模型0.271.040.210.141.631.020.63

      表 6  缩比模型轴向方向峰值压力绝对值偏差

      Table 6.  Absolute deviation of peak pressure of scale model in the axial direction

      弹药类型缩比模型径向峰值压力偏差/%
      1/2缩比模型0.020.080.010.100.020.431.19
      1/3缩比模型1.732.481.762.291.391.801.16
      1/4缩比模型0.060.030.100.580.290.670.66
      1/5缩比模型0.100.050.170.930.890.241.73

      表 7  缩比模型径向方向峰值压力绝对值偏差

      Table 7.  Absolute deviation of peak pressure of scale model in the radial direction

      表67可知,对于不同缩比模型,轴向方向上峰值压力绝对值偏差最大为2.77%,径向方向上峰值压力绝对值偏差最大为2.48%,属于数值模型建立及网格划分所造成的误差范围之内[7],因此可以说明通过量纲分析所建立的相似性模型是正确的。

    • 为进一步验证通过量纲分析所建立的弹药相似性模型,在进行数值模拟的基础上,开展某型弹药战斗部冲击波测试试验。战斗部实体如图5所示,其装药TNT当量为112.5 kg。分别建立缩比模型A1、A2和A3,其对应的TNT质量m1=64 kg,m2=16 kg,m3=4 kg,则由量纲分析可知缩比比例分别为θ1=0.829,θ2=0.522,θ3=0.329,由此可根据实体战斗部参数建立相应缩比模型。

      图  5  弹药战斗部实体

      Figure 5.  Entity of ammunition warhead

      根据试验方案,在距离爆心14、19、24 m的3个点处分别放置3组(每组2个)超压测试传感器,传感器距离地面1.8 m,与爆心同高。在赤道面夹角α和轴向夹角ϕ同实际试验保持一致的前提下,利用缩比比例θ进行计算,可得缩比模型中同实际测量点处具有相同冲击波峰值压力的测量距离如表8所示。

      弹药类型测量距离/m
      原型14.0019.0024.00
      缩比模型A111.6115.7519.90
      缩比模型A2 7.31 9.9212.53
      缩比模型A3 4.61 6.25 7.90

      表 8  冲击波峰值压力测量距离

      Table 8.  Measure distance of peak pressure

      试验布置示意图及现场图如图6所示。试验共使用2发弹药战斗部,分别记录每发爆炸后各个传感器冲击波超压值,并将其换算为冲击波峰值压力值,计算各次测量平均值作为不同距离处的冲击波峰值压力值。试验结果与仿真结果如表9所示,由此可得仿真试验冲击波峰值压力同实际试验冲击波峰值压力误差值如表10所示。

      图  6  冲击波测量试验布置图

      Figure 6.  Test layout of shock wave measurement

      试验类型冲击波峰值压力/kPa
      原型实际试验257.3185.5145.3
      缩比模型A1仿真试验248.2176.5137.2
      缩比模型A2仿真试验246.0176.0137.2
      缩比模型A3仿真试验247.5177.0136.7

      表 9  现场试验与仿真试验冲击波峰值压力

      Table 9.  Peak pressure in field and simulation tests

      试验
      类型
      误差/%
      缩比模型A1仿真试验3.54.95.6
      缩比模型A2仿真试验4.45.15.6
      缩比模型A3仿真试验3.84.65.9

      表 10  冲击波峰值压力误差

      Table 10.  Error of peak pressure

      由试验结果可知,通过缩比模型仿真试验所获得的冲击波峰值压力与实际试验所得峰值压力误差最小为3.5%,最大为5.9%,随着测量距离的增大,两者间误差有所增加。考虑误差产生原因,仿真试验是在完全理想化的条件下进行,而实际试验则受到天气、地形以及人为因素等诸多因素的影响,故产生的误差属于试验误差的合理范围之内。因此通过缩比模型仿真试验与原型实际试验的对比分析,证明了弹药相似性模型的正确性。

    • (1)当缩比模型装药量为原型的θ3倍时,其在θ倍爆距处所产生的冲击波峰值压力与弹药原型所产生的冲击波峰值压力相同。因此在满足相似律的条件下,冲击波峰值压力的测试可以用弹药缩比模型替代原型。

      (2)随着缩比比例的增大,冲击波在单位尺度上的变化率迅速增大,这将给实际测量工作带来困难,降低试验测量精度。同时,缩比比例越大,缩比模型的外壳越薄,这不仅会对模型的加工制造带来困难,也会对材料的强度属性造成一定的影响。因此在实际应用中需要合理选择缩比比例。

      (3)通过对弹药原型试验同缩比模型仿真试验结果进行对比分析,得到缩比模型仿真试验同弹药原型试验冲击波峰值压力误差在试验允许范围之内,能够较好地反映原型试验的试验效果,证明了弹药相似性模型的正确性。

参考文献 (10)

目录

    /

    返回文章
    返回