材料特性对爆轰驱动飞层对碰凸起影响的数值模拟

刘军; 王言金; 冯其京; 张树道

doi:10.11858/gywlxb.2015.01.011

材料特性对爆轰驱动飞层对碰凸起影响的数值模拟

doi: 10.11858/gywlxb.2015.01.011

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

基金项目: 国家自然科学基金(11072040)；中国工程物理研究院科学技术发展基金(2012B0202029, 2012A0202010)

详细信息

作者简介:
刘军(1981-), 男, 硕士, 助理研究员, 主要从事多介质计算方法和爆轰驱动数值模拟研究. E-mail:caepcfd@126.com

中图分类号: O358
计量
- 文章访问数: 6573
- HTML全文浏览量: 2013
- PDF下载量: 262
出版历程
- 收稿日期: 2013-06-13
- 修回日期: 2013-08-15

Simulation Study on Effects of Material Properties of the Bulging Flyers Driven by Colliding Detonation Waves

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

摘要

摘要: 应用多介质弹塑性L-R两步欧拉计算方法，对炸药两端起爆情况下柱面爆轰驱动飞层对碰凸起现象进行了数值模拟。在飞层为铅材料情况下，给出了铅飞层熔化区域及断裂图像，将对碰凸起计算结果与实验X光图像进行比较的基础上，进一步探索凸起现象的形成机制。从材料的熔点、弹塑性、密度、声速等方面对金属飞层对碰区凸起的影响进行了分析。综合分析几个计算模型的数值模拟结果，初步认为：材料强度大、密度大、熔点高能使飞层对碰凸起钝化，而材料声速较小会造成对碰区波系关系发生变化且正压作用时间增加，从而形成较严重的对碰凸起，材料声速大小对对碰凸起的影响较强度、密度等因素更明显。
- 爆轰驱动 /
- 熔化 /
- 断裂 /
- 欧拉方法
Abstract: The elastic-plastic hydrodynamic L-R two-step Euler method was used for studying metal flyers colliding behaviors driven by cylinder explosion.The simulation results just accorded with the radio graphs when the flyer metal was lead, and the melting and spall phenomenons were described by simulation.And then, the formation mechanism of bulging part was researched.The material properties including melting points, strength, density and sound velocity, were investigated in this study.By analyzing the simulation results, some phenomenons are found out:the strong strength brings on an blunt bulging area; and the material sound velocity affects the pressure working time on colliding area, which is much more sensitive than the density and strength.
- explosion drived /
- melt /
- spall /
- Euler method

HTML全文

图 1 实验模型(a)及计算模型(b)

Figure 1. The experimental setup (a) and simulation model (b)

下载: 全尺寸图片幻灯片

图 2 实验X光照相(a)及t=13 μs数值模拟密度图像(b)

Figure 2. Comparison between experimental result(a) and simulation result(b)

下载: 全尺寸图片幻灯片

图 3 t=13 μs时刻各飞层温度图像(a)及铅飞层熔化区域(b)

Figure 3. The temperature of metal flyers(a), and the melting part of Pb flyer(b)

下载: 全尺寸图片幻灯片

图 4 t=13 μs时刻模型(a)与模型(b)的计算结果密度图

Figure 4. Simulation results of Model(a) and Model(b) at t=13 μs

下载: 全尺寸图片幻灯片

图 5 t=13 μs时刻模型(c)与模型(d)的计算结果密度图及与铅模型对比

Figure 5. Simulation results of model(c) and model(d) at t=13 μs, and the comparison with Pb model

下载: 全尺寸图片幻灯片

图 6 t=6.0~8.5 μs时刻铜飞层、铅飞层在对碰区的密度、压力对比图像

Figure 6. Comparison of density and pressure between Cu flyer and Pb flyer at t=6.0-8.5 μs

下载: 全尺寸图片幻灯片

表 1 PETN炸药JWL参数

Table 1. The JWL EOS parameters of PETN

ρ/(kg/m³)	p_CJ/(GPa)	D/(m/s)	A/(GPa)	B/(GPa)	r₁	r₂	w
1 770	33.5	8 300	617.0	16.9	4.4	1.2	0.25

下载: 导出CSV

表 2 金属相关计算参数

Table 2. The simulation parameters of some metals

Metal	ρ₀/(g/cm³)	c₀/(km/s)	S₁	γ	α	G₀/(GPa)	Y₀/(GPa)	Y_max/(GPa)	T_m0/(K)
LY-12Al	2.79	5.328	1.34	2.00	0.47	28.6	0.26	0.76	1 222
Pb	11.34	2.100	1.46	2.74	0.54	8.6	0.01	0.10	760
304 steel	7.90	4.500	1.49	1.93	0.53	77.0	0.34	2.50	2 380
Cu	8.93	3.940	1.50	2.02	0.47	47.7	0.12	0.64	1 790
U	18.95	2.487	1.56	2.32	0.52	86.7	0.40	1.68	1 710

下载: 导出CSV

参考文献(13)

[1]	陈军, 孙承纬, 蒲正美, 等.爆轰波对碰区产物驱动金属圆管的研究[J].爆炸与冲击, 2003, 23(5): 442-447. http://d.wanfangdata.com.cn/Periodical/bzycj200305010 Chen J, Sun C W, Pu Z M, et al. Expansion of metallic tubes driven by detonation product behind two head-on colliding detonation waves[J]. Explosion and Shock Waves, 2003, 23(5): 442-447. (in Chinese) http://d.wanfangdata.com.cn/Periodical/bzycj200305010
[2]	王裴, 秦承森, 张树道, 等. SPH方法对金属表面微射流的数值模拟[J].高压物理学报, 2004, 18(2): 149-156. http://www.cqvip.com/Main/Detail.aspx?id=9991884 Wang P, Qin C S, Zhang S D, et al. Simulated microjet from free surface of aluminum using smoothed particle hydrodynamics[J]. Chinese Journal of High Pressure Physics, 2004, 18(2): 149-156. (in Chinese) http://www.cqvip.com/Main/Detail.aspx?id=9991884
[3]	袁帅, 胡海波, 张崇玉, 等.应用SPH方法对平板铅飞层对碰凸起动载行为的数值模拟研究[J].高压物理学报, 2010, 24(5): 383-387. http://www.cnki.com.cn/Article/CJFDTotal-GYWL201005011.htm Yuan S, Hu H B, Zhang C Y, et al. Simulation on head-on bulging dynamic behavior of lead flyer with SPH method[J]. Chinese Journal of High Pressure Physics, 2010, 24(5): 383-387. (in Chinese) http://www.cnki.com.cn/Article/CJFDTotal-GYWL201005011.htm
[4]	何长江, 于志鲁, 冯其京.高速碰撞的三维欧拉数值模拟方法[J].爆炸与冲击, 1999, 19(3): 216-221. http://www.cnki.com.cn/Article/CJFDTotal-BZCJ903.004.htm He C J, Yu Z L, Feng Q J. 3D Eulerian numerical simulation method of high speed impact[J]. Explosion and Shock Waves, 1999, 19(3): 216-221. (in Chinese) http://www.cnki.com.cn/Article/CJFDTotal-BZCJ903.004.htm
[5]	何长江, 周海兵, 杭义洪.爆轰驱动金属铝界面不稳定性的数值分析[J].中国科学, 2009, 39(9): 1170-1173. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkx-cg200909001 He C J, Zhou H B, Hang Y H. Numerical research of expansion driving interface instability of Al[J]. Science in China, 2009, 39(9): 1170-1173. (in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgkx-cg200909001
[6]	冯其京, 郝鹏程, 杭义洪, 等.聚能装药的欧拉数值模拟[J].爆炸与冲击, 2008, 28(2): 138-143. http://www.cqvip.com/Main/Detail.aspx?id=27101393 Feng Q J, Hao P C, Hang Y H, et al. Eulerian numerical simulation of a shaped charge[J]. Explosion and Shock Waves, 2008, 28(2): 138-143. (in Chinese) http://www.cqvip.com/Main/Detail.aspx?id=27101393
[7]	刘军, 何长江, 梁仙红.三维弹塑性流体力学自适应欧拉方法研究[J].高压物理学报, 2008, 22(1): 72-78. Liu J, He C J, Liang X H. An Eulerian adaptive mesh refinement method for three dimensional elastic-plastic hydrodynamic simulations[J]. Chinese Journal of High Pressure Physics, 2008, 22(1): 72-78. (in Chinese)
[8]	Chandle E A, Egan P O, Stokes J, et al. Application of Pegasus Ⅱ to study inertial instability and fracture of cylindrical tubes of solid aluminum, UCRL-ID-134657[R]. Livermore, CA: Lawrence Livermore National Laboratory, 1999.
[9]	张崇玉, 胡海波, 李庆忠, 等.爆轰波对碰驱动下平面铅飞层对碰区动载行为实验研究[J].高压物理学报, 2009, 23(4): 283-287. http://www.cqvip.com/QK/96553X/200904/32713244.html Zhang C Y, Hu H B, Li Q Z, et al. Experimental study on dynamic behavior of lead plate driven by two head-on colliding detonation waves[J]. Chinese Journal of High Pressure Physics, 2009, 23(4): 283-287. (in Chinese) http://www.cqvip.com/QK/96553X/200904/32713244.html
[10]	张崇玉, 胡海波, 李庆忠, 等.柱面内爆对碰驱动下铅飞层对碰区凸起行为研究[C]//第十届全国冲击动力学学术会议.太原, 2011. Zhang C Y, Hu H B, Li Q Z, et al. Bulging dynamic behavior of lead flyer driven by cylinder explosion[C]//The 10th Impact Dynamics Conference. Taiyuan, 2011. (in Chinese)
[11]	张崇玉.对碰爆轰波驱动下几种金属平面飞层的动载行为[R].绵阳: 中国工程物理研究院, 2009. Zhang C Y. Dynamic behavior of some metal flyers driven by two head-on colliding detonation waves[R]. Mianyang: CAEP, 2009. (in Chinese)
[12]	李德元, 徐国荣, 水鸿寿, 等.二维非定常流体力学数值方法[M].北京: 科学出版社, 1998. Li D Y, Xu G R, Shui H S, et al. Numerical Simulation Method of 2D-Unsteady Fluid Flow[M]. Beijing: Science Press, 1998. (in Chinese)
[13]	李茂生, 陈栋泉.高温高压下材料的本构模型[J].高压物理学报, 2001, 15(1): 24-31. http://www.cnki.com.cn/Article/CJFDTotal-GYWL200101003.htm Li M S, Chen D Q. A constitutive model for materials under high temperature and pressure[J]. Chinese Journal of High Pressure Physics, 2001, 15(1): 24-31. (in Chinese) http://www.cnki.com.cn/Article/CJFDTotal-GYWL200101003.htm