利用前沿推进法计算复杂炸药结构的起爆时间

王瑞利; 林忠; 魏兰; 刘学哲

doi:10.11858/gywlxb.2015.04.008

利用前沿推进法计算复杂炸药结构的起爆时间

doi: 10.11858/gywlxb.2015.04.008

北京应用物理与计算数学研究所，北京 100094

基金项目: 国家自然科学基金(11372051，11201035);中国工程物理研究院科学技术基金(2013A0101004)

详细信息

作者简介:
王瑞利(1964-)，男，研究员，主要从事计算流体力学与应用软件研究.E-mail: wang_ruili@iapcm.ac.cn

通讯作者:
刘学哲(1976-)，男，副研究员，主要从事计算流体力学与应用软件研究.E-mail: liu_xuezhe@iapcm.ac.cn

中图分类号: O382
计量
- 文章访问数: 5922
- HTML全文浏览量: 2394
- PDF下载量: 198
出版历程
- 收稿日期: 2014-03-19
- 修回日期: 2014-04-28

Calculating the Initiation Time of the Explosive with Complex Structure Using the Advancing Front Technique

Beijing Institute of Applied Physics and Computational Mathematics, Beijing 100094, China

摘要

摘要: 在各种工程定向爆破技术和爆轰数值模拟中，有时需要预先知道从起爆位置到装药内部某位置的起爆时间。为此，基于惠更斯原理，将网格生成前沿推进法的思想引入炸药爆轰计算中，给出了一种从炸药起爆位置到装药内部任意位置起爆时间的计算方法，通过几个数值算例验证了其有效性。该计算方法的特点是对于任何复杂装药结构均可给出从起爆位置到装药内部任意位置的起爆时间。
- 惠更斯原理 /
- 前沿推进法 /
- 复杂几何区域 /
- 起爆时间
Abstract: In the directional blasting in various engineering area and the numerical simulation of detonation, it is sometimes necessary to estimate the initiation time from the initiation position to a specific point of the explosive.Based on the Huygens' principle, the advancing front technique is introduced in the numerical mesh generation in order to calculate the detonation in the explosive.Furthermore, a computational method is proposed to evaluate the initiation time of the explosive from the initiation position to an arbitrary position of the explosive.Compared with traditional algorithms, this method is suitable for the explosive with complex geometric domain when coping with the initiation time.Several numerical examples are given to verify the effectiveness of this method.
- Huygens' principle /
- advancing front technique /
- complex geometric domain /
- initiation time

HTML全文

图 1 惠更斯原理示意图

Figure 1. Schematic of Huygens' principle

下载: 全尺寸图片幻灯片

图 2 炸药计算区域与初始网格

Figure 2. Computational domain and the initial meshes of explosives

下载: 全尺寸图片幻灯片

图 3 计算单元前沿的定义

Figure 3. Defining the front of calculation mesh

下载: 全尺寸图片幻灯片

图 4 起爆时间前沿最小时间单元的推进

Figure 4. Advancing front along the cell of minimum initiation time

下载: 全尺寸图片幻灯片

图 5 遇到其它物体时计算前沿的推进过程

Figure 5. Advancing front process when coming across an object

下载: 全尺寸图片幻灯片

图 6 3类网格体系

Figure 6. Three kinds of grid systems

下载: 全尺寸图片幻灯片

图 7 由3类网格计算得到的起爆时间等值线(单位:μs)

Figure 7. Contours of the initiation time corresponding to 3 kinds of grid systems (Unit:μs)

下载: 全尺寸图片幻灯片

图 8 绕爆模型及起爆时间等值线(单位:μs)

Figure 8. Explosion model around the substances and the contour of initiation time (Unit:μs)

下载: 全尺寸图片幻灯片

图 9(a) 计算模型结构和初始条件

Figure 9(a). Model structure and initial condition

下载: 全尺寸图片幻灯片

图 9(b) 起爆时间等值线(单位:μs)

Figure 9(b). Contour of initiation time (Unit:μs)

下载: 全尺寸图片幻灯片

图 10 不同时刻的压力等值线

Figure 10. Contours of pressure at different times

下载: 全尺寸图片幻灯片

图 11 平面绕爆模型计算网格和压力等值线

Figure 11. Mesh and contours of pressure at different times

下载: 全尺寸图片幻灯片

参考文献(11)

[1]	章冠人, 陈大年.凝聚炸药起爆动力学[M].北京: 国防工业出版社, 1991. Zhang G R, Chen D N. Condensed Explosives Detonation Dynamics[M]. Beijing: National Defense Industry Press, 1991. (in Chinese)
[2]	Sun J S, Zhu J S. Theory of Detonation Physics[M]. Beijing: National Defense Industry Press, 1995. (in Chinese) Sun J S, Zhu J S. Theory of Detonation Physics[M]. Beijing: National Defense Industry Press, 1995. (in Chinese)
[3]	孙锦山.凝聚炸药非理想爆轰的数值模拟[J].力学进展, 1995, 25(1): 127-133. http://qikan.cqvip.com/Qikan/Article/Detail?id=1824896 Sun J S. Numerical modeling of non-ideal detonation in condensed explosives[J]. Advances in Mechanics, 1995, 25(1): 127-133. (in Chinese) http://qikan.cqvip.com/Qikan/Article/Detail?id=1824896
[4]	Wilkins M L. Computer Simulation of Dynamic Phenomena[M]. Heidelberg, Germany: Springer-Verlag, 1999.
[5]	刘尔岩, 王元书.高能炸药散心爆轰波绕射传播的数值模拟[J].爆炸与冲击, 2000, 20(1): 64-67. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bzycj200001011 Liu E Y, Wang Y S. Numerical simulation of diffraction of a divergent detonation[J]. Explosive and Shock Waves, 2000, 20(1): 64-67. (in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bzycj200001011
[6]	尧德中.波动方程积分解的简化与惠更斯原理[J].电子科技大学学报, 1996, 25(3): 245-250. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199600042446 Yao D Z. Simplification of the integral solution of wave equation and Huygens' principle[J]. Journal of University of Electronic Science and Technology of China, 1996, 25(3): 245-250. (in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199600042446
[7]	朱振和.惠更斯原理的修正[J].光学技术, 1996(2): 36-38. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199600147369 Zhu Z H. Correction of the Huygen's principle[J]. Optical Technology, 1996(2): 36-38. (in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199600147369
[8]	Lo S H. Volume discretization into tetrahedra Ⅱ. 3D triangulation by advancing front approach[J]. Comput Struct, 1991, 39(5): 501-511. doi: 10.1016/0045-7949(91)90059-U
[9]	Lohner R, Parikh P, Gurnbert C. Interactive generation of unstructured grids for three dimensional problems[C]//Numerical Grid Generation in Computational Fluid Mechanics '88. Swansea UK: Pineridge Press, 1988: 687-697.
[10]	王瑞利, 林忠, 闫伟.多介质问题分析的前处理程序PreGenGrid[J].应用数学和力学, 2013, 34(5): 535-545. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yysxhlx201305012 Wang R L, Lin Z, Yan W. Simulations of the multi-medium problem pre-process: PreGenGrid[J]. Applied Mathematics and Mechanics, 2013, 34(5): 535-545. (in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yysxhlx201305012
[11]	王瑞利, 林忠, 倪国喜.基于任意多边形拉氏网格的有限体积方法研究及应用[J].数值计算与计算机应用, 2006, 27(1): 31-38. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=szjsyjsjyy200601004 Wang R L, Lin Z, Ni G X. Base on arbitrary n-polygon Lagrange grids finite volume method and applications[J]. Journal on Numerical Methods and Computer Applications, 2006, 27(1): 31-38. (in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=szjsyjsjyy200601004