Simulation Study on Effects of Material Properties of the Bulging Flyers Driven by Colliding Detonation Waves
-
摘要: 应用多介质弹塑性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.
-
Key words:
- explosion drived /
- melt /
- spall /
- Euler method
-
表 1 PETN炸药JWL参数
Table 1. The JWL EOS parameters of PETN
ρ/(kg/m3) pCJ/(GPa) D/(m/s) A/(GPa) B/(GPa) r1 r2 w 1 770 33.5 8 300 617.0 16.9 4.4 1.2 0.25 表 2 金属相关计算参数
Table 2. The simulation parameters of some metals
Metal ρ0/(g/cm3) c0/(km/s) S1 γ α G0/(GPa) Y0/(GPa) Ymax/(GPa) Tm0/(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 -
[1] 陈军, 孙承纬, 蒲正美, 等.爆轰波对碰区产物驱动金属圆管的研究[J].爆炸与冲击, 2003, 23(5): 442-447. http://d.wanfangdata.com.cn/Periodical/bzycj200305010Chen 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=9991884Wang 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.htmYuan 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.htmHe 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-cg200909001He 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=27101393Feng 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.htmlZhang 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.htmLi 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