Experimental Study and Numerical Simulation of Plane Flyer Driven by Detonation of JO-9159 and JB-9014 Composite Charge
-
摘要: 为了研究串联复合药柱(JO-9159/JB-9014)结构尺寸对能量输出的影响,采用有限元软件AUTODYN对标准平面飞片实验进行数值模拟,并进行了实验验证,结果表明,飞片速度实验值与计算值的相对误差为0.2%~3.0%,比动能相对误差为0.4%~6.0%,因此模型是可信的。利用该模型及材料物性参数,对不同高度比的串联复合装药结构进行数值计算,研究结构尺寸变化和复合装药能量释放的规律,得到高能炸药和钝感炸药尺寸比与飞片速度的指数关系公式。数值模拟研究表明, 随着高能炸药组成增加,爆轰驱动飞片的第1峰值速度和第2峰值速度越来越接近,而钝感炸药组成较大时,第1峰值速度较第2峰值速度较小,整个速度历史随着时间的推移有较大的跃升过程。Abstract: In order to research the characteristic of how the scale of series composite charge explosives (JO-9159/JB-9014) affect energy output, FEM software AUTODYN was used to simulate the plane flyer test, and relevant experiments were carried out to verify the result of analog.It is showed that the relative error of the velocity of the flyer between experimental value and calculated value is 0.2%-3.0% and the relative error of the ratio dynamic energy is 0.4%-6.0%.Therefore, the simulation model is reasonable.And then the model and the material parameters were applied to calculate the structures of series composite charge explosives in different height ratios.From the calculation, the relation of structure size and energy output was studied and the exponential function between the size of the explosives and the velocity of the flyer was obtained.In conclusion, the first peak velocity of the flyer grows closer to the second as the high explosive composition increases.In contrast, the first peak velocity is less than the second as the insensitive explosive composition increases, and the velocity has a major zoom during the whole process.
-
表 1 炸药的CJ爆轰参数
Table 1. CJ detonation parameters of explosive
Explosive ρ/(g/cm3) D/(km/s) uJ/(km/s) cJ/(km/s) JO-9159 1.863 8.87 2.22 6.65 JB-9014 1.900 7.71 1.93 5.78 表 2 炸药爆轰产物JWL状态方程参数
Table 2. JWL EOS parameters of detonation products
Explosive A/(GPa) B/(GPa) R1 R2 $\tilde{\omega}$ JO-9159 934.770 12.723 4.6 1.1 0.37 JB-9014 531.396 2.703 4.1 1.1 0.46 表 3 模拟计算结果
Table 3. Simulated results
No. JO-9159
h1/(mm)JB-9014
h2/(mm)Flyer velocity
uf/(km/s)Specific kinetic energy
e/(kJ/g)1 0 80 3.625 6.570 2 30 50 3.975 7.900 3 40 40 4.156 8.636 表 4 实验结果
Table 4. Experimental results
No. JO-9159
h1/(mm)JB-9014
h2/(mm)Flyer velocity
uf/(km/s)Specific kinetic energy
e/(kJ/g)1 0 80 3.691 6.812 2 30 50 3.854 7.427 3 40 40 4.164 8.669 表 5 飞片试验模拟结果
Table 5. Simulated results of flyer experiment
No. JO-9159
h1/(mm)JB-9014
h2/(mm)Flyer velocity
uf/(km/s)Specific kinetic energy
e/(kJ/g)1 20 60 3.851 7.415 2 25 55 3.905 7.625 3 30 50 3.975 7.900 4 35 45 4.078 8.315 5 40 40 4.156 8.636 6 45 35 4.228 8.938 7 50 30 4.280 9.159 8 55 25 4.352 9.470 9 60 20 4.444 9.875 -
[1] 梁增友, 黄风雷, 张震宇. PBX-9404的化学反应速率方程及起爆特性[J].爆炸与冲击, 2008, 28(1): 38-43. http://d.wanfangdata.com.cn/Periodical/bzycj200801007Liang Z Y, Huang F L, Zhang Z Y. Study on new reaction rate function model of PBX_9404 for damaged explosive initiation behavior[J]. Explosion and Shock Waves, 2008, 28(1): 38-43. (in Chinese) http://d.wanfangdata.com.cn/Periodical/bzycj200801007 [2] 于川, 刘文翰, 李良忠, 等.钝感炸药圆筒试验与爆轰产物JWL状态方程研究[J].高压物理学报, 1997, 11(3): 227-233. http://www.cnki.com.cn/article/cjfd1997-gywl703.011.htmYu C, Liu W H, Li L Z, et al. Studies on cylinder test and JWL equation of state of detonation product for insensitive high explosive[J]. Chinese Journal of High Pressure Physics, 1997, 11(3): 227-233. (in Chinese) http://www.cnki.com.cn/article/cjfd1997-gywl703.011.htm [3] Fried L E, Murphy M J, Souers P C, et al. Detonation modeling with an in-line thermo chemical equation of state[C]//Proceedings of 11th International Symposium on Detonation. Colorado: AFRL, 1998: 889-896. [4] Vitello P, Souers P C. Stability affects of artificial viscosity in detonation modeling[C]//Proceedings of 12th International Symposium on Detonation. San Diego, CA: NSWC, 2002: 649-656. [5] Bahl K L, Breithaupt R D, Tarver C M, et al. Fabry perot velocimetry on detonating LX-17 in planar and spherically divergent geometries[C]//Proceedings of 9th International Symposium on Detonation. Portland, Oregon, 1989: 133-141. [6] 虞德水, 赵锋, 谭多望, 等. JOB-9003和JB-9014炸药平面爆轰驱动飞片的对比研究[J].爆炸与冲击, 2006, 26(2): 140-144. http://www.cnki.com.cn/Article/CJFDTotal-BZCJ200602007.htmYu D S, Zhao F, Tan D W, et al. Experimental studies on detonation driving behavior of JOB-9003 and JB-9014 slab explosives[J]. Explosion and Shock Waves, 2006, 26(2): 140-144. (in Chinese) http://www.cnki.com.cn/Article/CJFDTotal-BZCJ200602007.htm [7] 张世文, 华劲松, 刘仓理, 等.金属圆管内爆轰波相互作用效应的数值模拟研究[J].爆炸与冲击, 2004, 24(3): 219-225. http://www.cnki.com.cn/Article/CJFDTotal-BZCJ200403005.htmZhang S W, Hua J S, Liu C L, et al. Research in numerical simulation on the effect under two detonation waves impacted in metallic tube[J]. Explosion and Shock Waves, 2004, 24(3): 219-225. (in Chinese) http://www.cnki.com.cn/Article/CJFDTotal-BZCJ200403005.htm [8] 孙承纬, 卫玉章, 周之奎.应用爆轰物理[M].北京: 国防工业出版社, 2000: 286-296.Sun C W, Wei Y Z, Zhou Z K. Application of Detonation Physics[M]. Beijing: National Defense Industry Press, 2000: 286-296. (in Chinese) [9] 张万甲, 吴强, 刘光祚.高-低爆速圆板炸药串联爆轰引起平面爆轰波的变凸现象[J].爆炸与冲击, 2007, 27(3): 259-265. http://www.cqvip.com/Main/Detail.aspx?id=24662448Zhang W J, Wu Q, Liu G Z. Analysis on a plane detonation wave changed into a convex wave induced by series-detonation of cylindrical explosives with low and high detonation velocities[J]. Explosion and Shock Waves, 2007, 27(3): 259-265. (in Chinese) http://www.cqvip.com/Main/Detail.aspx?id=24662448