水下爆炸冲击波数值仿真研究

胡亮亮; 黄瑞源; 李世超; 秦健; 王金相; 荣光

doi:10.11858/gywlxb.20190773

水下爆炸冲击波数值仿真研究

doi: 10.11858/gywlxb.20190773

胡亮亮^1,,
黄瑞源^1,*, ,,
李世超¹,
秦健²,
王金相¹,
荣光¹

1.
南京理工大学瞬态物理国家重点实验室，江苏南京　210094
2.
海军研究院，北京　100161

基金项目: 国家自然科学基金（11802001，11402266，11672138）；装备预研基金（614260404021801）；中国空气动力研究与发展中心超高速碰撞研究中心开放基金（20190303）；中央高校基本科研业务费专项资金项目（30916011348）

详细信息

作者简介:
胡亮亮（1995－），男，硕士研究生，主要从事水下爆炸数值仿真研究.E-mail: 18362906302@163.com

通讯作者:
黄瑞源（1984－），男，博士，讲师，主要从事爆炸力学与冲击动力学研究.E-mail: ryhuang@njust.edu.cn

中图分类号: O383
计量
- 文章访问数: 10381
- HTML全文浏览量: 3638
- PDF下载量: 90
出版历程
- 收稿日期: 2019-05-10
- 修回日期: 2019-05-28

Shock Wave Simulation of Underwater Explosion

1.
National Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Nanjing 210094, Jiangsu, China
2.
Naval Research Academy, Beijing 100161, China

摘要

摘要: 由于在水下爆炸冲击波的数值仿真研究中，水的状态方程、人工黏性系数和网格尺寸对数值计算结果影响很大，采用常规TNT炸药的水下爆炸为例，以冲击波的峰值压力和比冲量为衡量指标，研究了这3个主要影响因素对数值仿真结果的影响。首先，通过采用常用的5种水的状态方程进行系列仿真，给出了各种状态方程的适用范围；其次，讨论了人工黏性系数对计算结果的影响，并给出了一次与二次人工黏性系数的建议取值范围；最后，通过对不同炸药当量及不同网格尺寸开展系列运算，从而得到不同炸药当量在满足工程计算精度要求下所对应的建议网格尺寸，并得到了不同炸药当量所对应的建议网格尺寸的表达式。
- 水下爆炸 /
- 数值仿真 /
- 水的状态方程 /
- 人工黏性系数 /
- 网格尺寸
Abstract: The state equation of water, artificial viscosity coefficient and mesh size have a great influence on the numerical results of underwater explosion shock wave. In order to improve the simulation accuracy of underwater explosion shock wave, the peak pressure and specific impulse of the conventional TNT explosive underwater explosion are taken as the measurement indicators, and the influence of these factors on the numerical simulation results is studied. For the five kinds commonly state equations of water, the specific values of the artificial viscosity coefficients under different working conditions and appropriate grid size for different explosive equivalents are given. These parameters can provide reference for improving simulation accuracy of underwater explosion shock wave under different working conditions. First, through a series of simulations of the commonly used five kinds of state equations of water, the calculation results of peak pressure and specific impulse are compared with the empirical formula, and the error analysis is carried out to give the applicable scope of each state equation. Secondly, the influence of the artificial viscosity coefficient on the calculation results is discussed, and a series of calculations are carried out for the primary and secondary artificial viscosity coefficients under different working conditions. The recommended range of values for the primary and secondary artificial viscosity coefficients under different working conditions is given. Finally, through a series of calculations on 0.1, 0.5, 1, 10, 50, 100, 500 and 1 000 kg equivalent explosives and different grid sizes, the recommended mesh sizes corresponding to different explosive equivalents under the requirement of engineering calculation accuracy are obtained by limiting the relative error of peak pressure less than 10%. The expressions of the recommended mesh sizes corresponding to different explosive equivalents are also given.
- underwater explosion /
- numerical simulation /
- equation of state /
- artificial viscosity /
- element density

HTML全文

图 1 采用不同状态方程得到的峰值压力计算结果

Figure 1. Calculated peak pressure results using different equations of states

下载: 全尺寸图片幻灯片

图 2 采用不同状态方程得到的比冲量计算结果

Figure 2. Calculated specific impulse results obtained with different equations of states

下载: 全尺寸图片幻灯片

图 3 调整一次项黏性系数得到的计算结果

Figure 3. Calculated results obtained by adjusting the viscosity coefficient of the primary term

下载: 全尺寸图片幻灯片

图 4 调整二次项黏性系数得到的计算结果

Figure 4. Calculated results obtained by adjusting the viscosity coefficient of the quadratic term

下载: 全尺寸图片幻灯片

图 5 一次黏性系数对不同当量工况的影响

Figure 5. Effects of primary viscosity coefficient under different equivalent conditions

下载: 全尺寸图片幻灯片

图 6 采用不同网格下得到的峰值压力

Figure 6. Peak pressure obtained under different grids

下载: 全尺寸图片幻灯片

图 7 采用不同网格下得到的比冲量

Figure 7. Specific impulse obtained under different grids

下载: 全尺寸图片幻灯片

图 8 网格尺寸对不同炸药当量水下爆炸数值仿真的影响

Figure 8. Influence of grid size on different explosive equivalents numerical simulation of underwater explosion

下载: 全尺寸图片幻灯片

图 9 炸药当量与适宜网格尺寸的关系

Figure 9. Relationship between explosive equivalent and suitable grid size

下载: 全尺寸图片幻灯片

图 10 数值仿真得到的冲击波时程曲线与试验对比

Figure 10. Numerical and experimental comparison of shock wave time history curve

下载: 全尺寸图片幻灯片

表 1 TNT炸药的主要参数^[16]

Table 1. Main parameters of TNT explosives^[16]

A/GPa	B/GPa	R₁	R₂	ω	ρ₀/(kg∙m^–3)	e/(J∙kg^–1)
371.2	3.231	4.15	0.95	0.3	1630	4.19 ×10⁶

下载: 导出CSV

表 2 常用Mie-Grüneisen状态方程不同形式的参数

Table 2. Commonly used parameters of the Mie-Grüneisen equation of state

Equation of state	C₀/(km·s^–1)	S₁	S₂	S₃	${\gamma _0} $	a
SNL^[20]	1.647	1.92	0	0	0	0
HULL^[21]	1.483	1.75	0	0	0.28	0
Steiberg^[22]	1.480	2.56	–1.986	0.226 8	0.50	2.67

下载: 导出CSV

表 3 Autodyn程序提供的水多项式状态方程参数^{[12, 23]}

Table 3. Polynomial state equation parameters of water provided by the Autodyn program^{[12, 23]}

A₁/GPa	A₂ /GPa	A₃ /GPa	T₁ /GPa	T₂/GPa	B₀	B₁
2.2	9.54	14.57	2.2	0	0.28	0.28

下载: 导出CSV

表 4 Dytran中水的多项式状态方程的参数^[24–25]

Table 4. Polynomial state equations parameters of water in Dytran^[24–25]

a₁/GPa	a₂/GPa	a₃/GPa	b₀	b₁	b₂	b₃
2.002	9.224	8.767	0.493 4	1.393 7	0	0

下载: 导出CSV

表 5 1 kg炸药无限水域爆炸冲击波峰值压力

Table 5. Shock wave peak pressure of 1 kg explosive infinite water fields explosion MPa

State equation of water	R/R₀
State equation of water	7	10	15	20	25	30	35
Autodyn polynomial	165.7	95.6	56.7	40.6	31.3	24.6	20.3
Dytran polynomial	170.7	97.6	57.4	41.2	31.5	25.1	20.8
Steiberg	169.8	97.7	57.7	41.6	31.1	25.3	21.1
SNL	218.7	123.1	71.2	46.2	36.8	30.8	25.7
HULL	182.5	106.7	62.8	45.2	34.9	28.4	23.7
Empirical formula	196.7	115.2	68.3	49.4	38.4	31.2	26.2

下载: 导出CSV

表 6 1 kg炸药无限水域爆炸冲击波比冲量

Table 6. Shock wave impulse of 1 kg explosives infinite water fields explosion N·s·m^–2

State equation of water	R/R₀
State equation of water	7	10	15	20	25	30	35
Autodyn polynomial	12 522.9	9 116.8	6 355.1	4 919.6	4 033.4	3 429.3	2 989.7
Dytran polynomial	12 284.8	8 943.4	6 234.2	4 826.0	3 956.7	3 364.1	2 932.8
Steiberg	12 130.7	8 831.3	6 156.0	4 765.5	3 907.1	3 321.9	2 896.0
SNL	13 237.3	9 636.9	6 717.6	5 200.2	4 263.5	3 624.9	3 160.2
HULL	11 976.6	8 719.1	6 077.9	4 704.9	3 857.5	3 279.7	2 859.2
Empirical formula	14 007.7	10 197.8	7 108.6	5 502.9	4 511.7	3 835.9	3 344.1

下载: 导出CSV

表 7 5种常用水的状态方程的适用范围与特点

Table 7. Applicable scope and characteristics of common five kinds of state equations for water

State equation of water	Scope of application
Autodyn polynomial	The overall error of the peak pressure is too large, and the far field error is larger than the near field. The near field error of the specific impulse calculation is smaller than the far field.
Dytran polynomial	The near-field peak pressure error is small, the specific impulse error is large, the far-field peak pressure error is large, and the specific impulse error is small.
Steiberg	The near-field peak pressure error is small, the specific impulse error is large, the far-field peak pressure error is large, and the specific impulse error is small.
SNL	The mid-field attenuation of the pressure peak is fast, and the far-field error is small. The difference between the near-field and the far-field error is not large, which is suitable for far-field calculation.
Empirical formula	The pressure peak error is small overall, and the specific impulse error is gradually increased from near to far, suitable for near-field calculation.

下载: 导出CSV

表 8 人工黏性对计算结果的影响

Table 8. Artificial viscosity effects on calculation results

Artificial viscosity coefficient	Recommended range of values	Impact on calculation results
Primary viscosity coefficient	0.005–0.040	The peak pressure of the underwater shock wave has a great influence, the contrast impulse has little effect, the primary viscosity coefficient increases, and the peak pressure decreases.
Secondary viscosity coefficient	0.8–1.0	Less influence on peak pressure and specific impulse.

下载: 导出CSV

表 9 1 kg炸药采用不同网格尺寸得到的峰值压力

Table 9. Peak pressure of 1 kg explosive with different grid sizes MPa

Grid size/mm	R/R₀
Grid size/mm	7	10	15	20	25	30	35
2	188.8	110.6	65.6	47.4	36.8	30.0	25.2
5	179.0	104.8	62.1	44.9	34.9	28.4	23.9
8	161.3	94.5	56.0	40.5	31.4	25.6	21.5
12	135.7	81.5	46.1	35.7	26.5	21.5	18.1
15	112.1	65.7	38.9	27.1	21.9	17.8	14.9
20	85.6	52.4	28.7	20.7	16.1	13.1	11.0
Empirical formula	196.7	115.2	68.3	49.4	38.4	31.2	26.2

下载: 导出CSV

表 10 1 kg炸药采用不同网格尺寸得到的比冲量

Table 10. Specific impulse of 1 kg explosive with different grid sizes N·s·m^–2

Grid size/mm	R/R₀
Grid size/mm	7	10	15	20	25	30	35
2	12 228.8	10 197.7	6 205.8	4 803.9	3 938.6	3 348.7	2 919.4
5	12 004.6	8 902.6	6 092.0	4 715.9	3 866.5	3 287.3	2 865.9
8	12 172.7	8 739.4	6 177.3	4 781.9	3 920.6	3 333.3	2 906.0
12	11 752.4	8 861.8	5 964.1	4 616.8	3 785.3	3 218.3	2 805.7
15	11 472.3	8 555.9	5 821.9	4 506.8	3 695.0	3 141.6	2 738.8
20	11 822.5	8 351.9	5 999.6	4 644.4	3 807.8	3 237.4	2 822.4
Empirical formula	14 007.7	10 197.8	7 108.6	5 502.9	4 511.7	3 835.9	3 344.1

下载: 导出CSV

表 11 不同炸药当量下适宜网格尺寸

Table 11. Appropriate grid size for different explosive equivalents

Explosive equivalent/kg	Suitable grid size/mm
0.1	1.5
0.5	3
1	5
10	10
50	16
100	20
500	35
1 000	50

下载: 导出CSV

参考文献(36)

[1]	杨振, 沈晓乐. 爆破战斗部水中兵器爆炸威力评定方法研究 [J]. 爆破, 2015, 32(2): 51–53. doi: 10.3963/j.issn.1001-487X.2015.02.009 YANG Z, SHEN X L. Research on evaluation method of underwater blast brisance of weapon’s explosive [J]. Blasting, 2015, 32(2): 51–53. doi: 10.3963/j.issn.1001-487X.2015.02.009
[2]	张志华, 钟强晖, 李庆民. 小药量水下爆炸对水下目标的毁伤有效值评估 [J]. 兵工学报, 2009, 30(10): 1344–1348. ZHANG Z H, ZHONG Q H, LI Q M. Evaluation of damage valid value of underwater target detonated by small charge [J]. Acta Armamentarii, 2009, 30(10): 1344–1348.
[3]	SRINIVAS K A, UMAPATHI GOKUL K, PYDISETTY V K R, et al. Blast loading of underwater targets – a study through explosion bulge test experiments [J]. International Journal of Impact Engineering, 2015, 76: 189–195. doi: 10.1016/j.ijimpeng.2014.09.007
[4]	ZHANG J, SHI X H, SOARES C G. Experimental study on the response of multi-layered protective structure subjected to underwater contact explosions [J]. International Journal of Impact Engineering, 2016, 100: 23–34.
[5]	罗松林, 叶序双, 顾文彬, 等. 水下爆炸研究现状 [J]. 工程爆破, 1999, 5(1): 84–87. doi: 10.3969/j.issn.1006-7051.1999.01.023
[6]	梁向前. 水下爆破技术[M]. 化学工业出版社, 2013.
[7]	QIAN K J, GANG Y D. A finite element analysis of ship sections subjected to underwater explosion [J]. International Journal of Impact Engineering, 2011, 38(7): 558–566. doi: 10.1016/j.ijimpeng.2010.11.005
[8]	BRETT J M, YIANNAKOPOLOUS G. A study of explosive effects in close proximity to a submerged cylinder [J]. International Journal of Impact Engineering, 2008, 35(4): 206–225. doi: 10.1016/j.ijimpeng.2007.01.007
[9]	王军, 孙丰, 陈舸, 等. 水下爆炸载荷作用下MK46鱼雷结构动态响应分析 [J]. 鱼雷技术, 2013(4): 293–298.
[10]	北京工业学院八系《爆炸及其运用》编写组. 爆炸及其运用[M]. 北京: 国防工业出版社, 1979.
[11]	徐豫新, 王树山, 李园. 水下爆炸数值仿真研究 [J]. 弹箭与制导学报, 2009, 29(6): 95–97. doi: 10.3969/j.issn.1673-9728.2009.06.026 XU Y X, WANG S S, LI Y. Study on numerical simulation of the underwater explosive [J]. Journal of Projectiles Rockets, Missiles and Guidance, 2009, 29(6): 95–97. doi: 10.3969/j.issn.1673-9728.2009.06.026
[12]	刘科种, 徐更光, 辛春亮, 等. AUTODYN水下爆炸数值模拟研究 [J]. 爆破, 2009, 26(3): 18–21. doi: 10.3963/j.issn.1001-487X.2009.03.005 LlU K Z, XU G G, XlN C L, et al. Research on numerical simulation in underwater explosion by AUTODYN [J]. Blasting, 2009, 26(3): 18–21. doi: 10.3963/j.issn.1001-487X.2009.03.005
[13]	方斌, 朱锡, 张振华, 等. 水下爆炸冲击波数值模拟中的参数影响 [J]. 哈尔滨工程大学学报, 2005, 26(4): 416–424. FANG B, ZHU X, ZHANG Z H, et al. Effect of parameters in numerical simulation of underwater shock wave [J]. Journal of Harbin Engineering University, 2005, 26(4): 416–424.
[14]	张振华, 朱锡, 白雪飞. 水下爆炸冲击波的数值模拟研究 [J]. 爆炸与冲击, 2004, 24(2): 182–188. doi: 10.3321/j.issn:1001-1455.2004.02.014 ZHANG Z H, ZHU X, BAI X F. The study on numerical simulation of underwater blast wave [J]. Explosion and Shock Waves, 2004, 24(2): 182–188. doi: 10.3321/j.issn:1001-1455.2004.02.014
[15]	梁龙河, 曹菊珍, 王元书. 水下爆炸特性的一维球对称数值研究 [J]. 高压物理学报, 2002, 16(3): 199–203. doi: 10.3969/j.issn.1000-5773.2002.03.007 LIANG L H, CAO J Z, WANG Y S. One-dimensional numerical simulations of underwater spherical explosions [J]. Chinese Journal of High Pressure Physics, 2002, 16(3): 199–203. doi: 10.3969/j.issn.1000-5773.2002.03.007
[16]	DOBRATZ B M. Properties of chemical explosives and explosive simultants [J]. International Journal of Neuroscience, 1981, 51(3–4): 339–340.
[17]	盛振新, 刘荣忠, 郭锐. 水下爆炸冲击波相互作用的仿真分析 [J]. 火工品, 2012(3): 25–29. SHENU Z X, LlU R Z, GUO R. Study on the shock wave interaction of underwater explosons [J]. Initiators and Pyrotechnics, 2012(3): 25–29.
[18]	荣吉利, 李健, 杨荣杰, 等. 水下爆炸气泡脉动的实验及数值模拟 [J]. 北京理工大学学报, 2008, 28(12): 1035–1038. RONG J L, LI J, YANG R J, et al. Experiment and numerical simulation for the bubble impulse in underwater explosion [J]. Transactions of Beijing Institute of Technology, 2008, 28(12): 1035–1038.
[19]	李晓杰, 张程娇, 闫鸿浩, 等. 水下爆炸近场非均熵流的特征线差分解法 [J]. 爆炸与冲击, 2012, 32(6): 604–608. doi: 10.3969/j.issn.1001-1455.2012.06.008 LI X J, ZHANG C J, YAN H H, et al. Difference method of characteristics in isentropic flow of underwater explosion in near-field region [J]. Explosion and Shock Waves, 2012, 32(6): 604–608. doi: 10.3969/j.issn.1001-1455.2012.06.008
[20]	辛春亮, 秦健, 刘科种, 等. 基于LS-DYNA软件的水下爆炸数值模拟研究 [J]. 弹箭与制导学报, 2008, 28(3): 156–158. doi: 10.3969/j.issn.1673-9728.2008.03.047 XlN C L, QlN J, LlU K Z, et al. Research on UNDEX numerical simulation hased on LS-DYNA [J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2008, 28(3): 156–158. doi: 10.3969/j.issn.1673-9728.2008.03.047
[21]	STENBERG D J. Spherical explosions and the equation of state of water, UCID-2097 [R]. Livermore, CA: Lawrence Livermore National Laboratory,1987.
[22]	SHIN Y S, LEE M, LAM K Y, et al. Modeling mitigation effects of watershield on shock waves [J]. Shock and Vibration, 1998, 5(4): 225–234. doi: 10.1155/1998/782032
[23]	LIU M B, LIU G R, LAM K Y, et al. Smoothed particle hydrodynamics for numerical simulation of underwater explosion [J]. Computational Mechanics, 2003, 30(2): 106–118. doi: 10.1007/s00466-002-0371-6
[24]	梅群, 侯中华, 朱俊锋, 等. 水下爆炸冲击波压力时程的数值模拟 [J]. 河南科技大学学报(自然科学版), 2010, 31(4): 57–59. MEI Q, HOU Z H, ZHU J F, et al. Numerical simulation of underwater explosion shock wave [J]. Journal of Henan University of Science and Technology: Natural Science, 2010, 31(4): 57–59.
[25]	尹群. 水面舰船设备冲击环境与结构抗冲击性能研究 [D]. 南京: 南京航空航天大学, 2006. YIN Q. Studies on shock environment for equipments on surface ship and anti-shock characteristics of structures [D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2006. DOI: 10.7666/d.d037570.
[26]	COLE R H. Underwater explosions [M]. New York: Dover Publications, 1965.
[27]	ZAMYSHLYAEV B V, YAKOVLEV Y S. Dynamic loads in underwater explosion [J]. Dynamic Loads in Underwater Explosion, 1973.
[28]	辛春亮. 高能炸药爆炸能量输出结构的数值分析[D]. 北京: 北京理工大学, 2008.
[29]	杨坤, 陈朗, 伍俊英, 等. 计算网格与人工黏性系数对炸药水中爆炸数值模拟计算的影响分析 [J]. 兵工学报, 2014(Suppl 2): 237–243. YANG K, CHEN L, WU J Y, et al. The effects of computing grid and artificial viscosity coefficient on underwater explosion numerical simulation [J]. Acta Armamentarii, 2014(Suppl 2): 237–243.
[30]	KIM J H, SHIN H C. Application of the ALE technique for underwater explosion analysis of a submarine liquefied oxygen tank [J]. Ocean Engineering, 2008, 35(8/9): 812–822. doi: 10.1016/j.oceaneng.2008.01.019
[31]	张社荣, 李宏璧, 王高辉, 等. 水下爆炸冲击波数值模拟的网格尺寸确定方法 [J]. 振动与冲击, 2015, 34(8): 93–100. ZHANG S R, LI H B, WANG G H, et al. A method to determine mesh size in numerical simulation of shock wave of underwater explosion [J]. Journal of Vibration and Shock, 2015, 34(8): 93–100.
[32]	余晓菲, 刘土光, 张涛. 水下爆炸冲击波的载荷强度计算 [J]. 舰船科学技术, 2006, 28(5): 22–28. YU X F, LIU T G, ZHANG T. Computation of the blast loading strength of underwater explosion shock waves [J]. Ship Science and Technology, 2006, 28(5): 22–28.
[33]	方斌, 朱锡, 张振华. 水下爆炸冲击波载荷作用下船底板架的塑性动力响应 [J]. 哈尔滨工程大学学报, 2008, 29(4): 326–331. doi: 10.3969/j.issn.1006-7043.2008.04.002 FANG B, ZHU X, ZHANG Z H. Plastic dynanic response of ship hull grillage to underwater blast loading [J]. Journal of Harbin Engineering University, 2008, 29(4): 326–331. doi: 10.3969/j.issn.1006-7043.2008.04.002
[34]	吴国民, 周心桃, 肖汉林, 等. 水下爆炸数值仿真 [J]. 舰船科学技术, 2012, 34(9): 20–26. doi: 10.3404/j.issn.1672-7649.2012.09.004 WU G M, ZHOU X T, XIAO H L, et al. Numerical simulation of underwater explosion [J]. Ship Science & Technology, 2012, 34(9): 20–26. doi: 10.3404/j.issn.1672-7649.2012.09.004
[35]	肖秋平, 陈网桦, 贾宪振, 等. 基于AUTODYN的水下爆炸冲击波模拟研究 [J]. 舰船科学技术, 2009, 31(2): 38–43. doi: 10.3404/j.issn.1672-7649.2009.02.005 XIAO Q P, CHEN W H, JIA X Z, et al. Numerical study of underwater explosion shock wave based on AUTODYN [J]. Ship Science and Technology, 2009, 31(2): 38–43. doi: 10.3404/j.issn.1672-7649.2009.02.005
[36]	胡毅亭, 贾宪振, 饶国宁, 等. 水下爆炸冲击波和气泡脉动的数值模拟研究 [J]. 舰船科学技术, 2009, 31(2): 134–140. doi: 10.3404/j.issn.1672-7649.2009.02.027 HU Y T, JIA X Z, RAO G N, et al. Numerical study of underwater explosion shock wave and bubble pulse [J]. Ship Science and Technology, 2009, 31(2): 134–140. doi: 10.3404/j.issn.1672-7649.2009.02.027