Numerical simulation of pre-shock desensitization in TATB-based heterogeneous explosive
-
摘要: 为了对一种TATB基非均质炸药的预冲击起爆现象展开数值模拟研究,将基于冲击温度及压力的AWSD反应率模型耦合进二维结构网格拉氏弹塑性流体力学程序。利用炸药及其产物的冲击雨贡纽实验数据校验了未反应炸药及产物的状态方程参数,通过一维冲击起爆的模拟,标定了反应速率模型参数。模拟了炸药在弱冲击0.45 μs后跟随的强冲击波的二次冲击实验,结果表明,受预压缩区域的炸药反应变慢,到爆轰距离增长了约1 mm,与该炸药二次冲击实验减敏现象相符。模拟拐角效应时,爆轰波经过拐角后,在拐角附近形成稳定的不起爆区域,与主要成分相同的LX-17炸药的拐角效应实验的死区特征相符。数值模拟结果表明,基于冲击温度及压力的AWSD反应率模型可以较好地模拟非均质炸药预冲击减敏问题。Abstract: To study and simulate the pre-shock desensitization in the TATB-based heterogeneous explosive, the impact temperature and pressure based AWSD reaction flow model was implemented in a 2D structured mesh Lagrangian elastoplastic hydrodynamics program. The reactant and product EOS parameters were calibrated against the Hugoniot experimental data. To calibrate the parameters of the reaction flow model, one-dimensional numerical simulations of the shock initiation experiments were carried out. We simulated the double-shock experiments in which a first weak shock was followed by a second strong shock with a time interval of 0.45 μs. The results indicate the reaction becomes slower in the precompressed region and the run-to-detonation distance is about 1 mm longer than that in the uncompressed region, which is consistent with the desensitization in double-shock experiments. When simulating the corner-turning, the detonation wave passes through the corner and forms a stable non-initiation region near the corner, which is consistent with the dead zone characteristics of the corner-turning experiment of the LX-17 explosive with the same main composition. The numerical simulation results show that the AWSD reaction rate model based on the impact temperature and pressure can well simulate the pre-shock desensitization of heterogeneous explosives.
-
Key words:
- TATB-based /
- heterogeneous explosive /
- pre-shock /
- desensitization /
- AWSD model
-
表 1 未反应TATB基炸药的Davis状态方程参数
Table 1. Davis EOS parameters of unreacted TATB-based explosive
A/(mm·μs−1) B C Z Γ0s E0/(kJ·g−1) αst cVs/(kJ·g−1·K−1) 1.93 4.26 0.30 0 0.56 3.80 0.757 0 0.000 967 表 2 TATB基炸药产物的Davis状态方程参数
Table 2. Davis EOS parameters of reaction products of TATB-based explosive
A k vc/(cm3·g−1) pc/GPa n b cVg/(kJ·g−1·K−1) 0.835 603 1.30 0.926 85 1.485 1 4.242 66 0.85 0.001 072 表 3 TATB基炸药AWSD模型参数
Table 3. Parameters of the AWSD model for the TATB-based explosive
np ps/GPa k1/s−1 T1/K a1 b1 b2 k2/s−1 T2/K fs λc δλ 0.651 5 27.60 336 1 724 0.060 81 2.122 0.9 10 200 6 278 0.035 87 0.876 4 0.021 68 表 4 二次冲击起爆实验的到爆轰距离和到爆轰时间的实验和数值模拟结果
Table 4. Experimental and simulated results of distance and time of run to detonation by double-shock initiation
实验 方法 ρ0/(g·cm−3) v0/(km·s−1) p0/GPa pm/GPa L*/mm T*/μs 3 实验 1.881 1.669 6.135 6 12.649 9 9.845 1.716 数值模拟 10.141 1.679 2 实验 1.881 1.732 6.470 6 12.706 4 9.145 1.356 数值模拟 9.240 1.505 -
[1] CAMPBELL A W, DAVIS W C, RAMSAY J B, et al. Shock initiation of solid explosives [J]. Physics of Fluids, 1961, 4(4): 511. DOI: 10.1063/1.1706354. [2] CAMPBELL A W, TRAVIS J R. The shock desensitization of PBX-9404 and composition B-3 [C] // Proceedings of the 8th International Detonation Symposium. Albuquerque: Naval Surface Weapons Center, 1985: 1057−1068. [3] MULFORD R N, SHEFFIELD S A, ALCON R R. Preshock desensitization of PBX explosives [J]. AIP Conference Proceedings, 1994, 309(1): 1405–1408. DOI: 10.1063/1.46243. [4] TARVER C M, COOK T M, URTIEW P A, et al. Multiple shock initiation of LX-17 [C] // Proceedings of the 10th International Detonation Symposium. Boston: Office of Naval Research, 1993: 676−703. [5] BAT’KOV Y V, GLUSHAK B L, NOVIKOV S A. Desensitization of pressed explosive compositions based on TNT, RDX, and HMX under double shock-wave loading [J]. Combustion, Explosion and Shock Waves, 1995, 31(4): 482–485. DOI: 10.1007/BF00789372. [6] GUSTAVSEN R L, SHEFFIELD S A, ALCON R R, et al. Double shock initiation of the HMX based explosive EDC-37 [J]. AIP Conference Proceedings, 2002, 620(1): 999–1002. DOI: 10.1063/1.1483706. [7] MADER C L. Numerical modeling of explosives and propellants [M]. 3rd ed. Boca Raton: CRC Press, 2008. [8] HUSSAIN T, LIU Y, HUANG F L, et al. Modeling and simulation of preshock desensitization in heterogeneous explosives using a mesoscopic reaction rate model [J]. Simulation, 2015, 91(11): 980–988. DOI: 10.1177/0037549715608962. [9] BORDZILOVSKII S A, KARAKHANOV S M. Desensitization of pressed RDX/paraffin and HMX/paraffin compounds by multiple shock waves [J]. Combustion, Explosion and Shock Waves, 1995, 31(2): 227–235. DOI: 10.1007/BF00755754. [10] LEE E L, TARVER C M. Phenomenological model of shock initiation in heterogeneous explosives [J]. The Physics of Fluids, 1980, 23(12): 2362–2372. DOI: 10.1063/1.862940. [11] JOHNSON J N, TANG P K, FOREST C A. Shock-wave initiation of heterogeneous reactive solids [J]. Journal of Applied Physics, 1985, 57(9): 4323–4334. DOI: 10.1063/1.334591. [12] WESCOTT B L, STEWART D S, DAVIS W C. Equation of state and reaction rate for condensed-phase explosives [J]. Journal of Applied Physics, 2005, 98(5): 053514. DOI: 10.1063/1.2035310. [13] DEOLIVEIRA G, KAPILA A K, SCHWENDEMAN D W, et al. Detonation diffraction, dead zones and the ignition and growth model [C] // Proceedings of the 13th International Detonation Symposium. Norfolk: Office of Naval Research, 2006. [14] HUSSAIN T, LIU Y, HUANG F L, et al. Desensitization by pre-shocking in heterogeneous explosives and its numerical modelling [J]. Central European Journal of Energetic Materials, 2016, 13(2): 357–379. DOI: 10.22211/cejem/64990. [15] 郝鹏程, 冯其京, 洪滔, 等. 钝感炸药点火增长模型的欧拉数值模拟 [J]. 爆炸与冲击, 2012, 32(3): 243–250. DOI: 10.11883/1001-1455(2012)03-0243-08.HAO P C, FENG Q J, HONG T, et al. Eulerian simulation on insensitive explosives with the ignition-growth reactive model [J]. Explosion and Shock Waves, 2012, 32(3): 243–250. DOI: 10.11883/1001-1455(2012)03-0243-08. [16] STARKENBERG J. Shock-pressure and pseudo-entropic approaches to explosive initiation modeling [C] // Proceedings of the 15th International Detonation Symposium. San Francisco, 2014: 908−916. [17] STARKENBERG J. Modeling detonation propagation and failure using explosive initiation models in a conventional hydrocode [C] // Proceedings of the 12th International Detonation Symposium. Annapolis: Office of Naval Research, 2002. [18] HANDELEY C A. CREST reactive flow model [C] // Proceedings of the 13th International Detonation Symposium. Norfolk, 2006. [19] DESBIENS N, MATIGNON C, SORIN R. Temperature-based model for condensed-phase explosive detonation [J]. Journal of Physics: Conference Series, 2014, 500(15): 152004. DOI: 10.1088/1742-6596/500/15/152004. [20] ASLAM T D. Shock temperature dependent rate law for plastic bonded explosives [J]. Journal of Applied Physics, 2018, 123(14): 145901. DOI: 10.1063/1.5020172. [21] 裴红波, 刘俊明, 张旭, 等. 基于反向撞击法的JB-9014炸药Hugoniot关系测量 [J]. 爆炸与冲击, 2019, 39(5): 052301. DOI: 10.11883/bzycj-2017-0395.PEI H B, LIU J M, ZHANG X, et al. Measurement of Hugoniot relation for unreacted JB-9014 explosive with reverse-impact method [J]. Explosion and Shock Waves, 2019, 39(5): 052301. DOI: 10.11883/bzycj-2017-0395. [22] 刘俊明, 张旭, 裴红波, 等. JB-9014钝感炸药冲击Hugoniot关系测量 [J]. 高压物理学报, 2018, 32(3): 033202. DOI: 10.11858/gywlxb.20170669.LIU J M, ZHANG X, PEI H B, et al. Measurement of Hugoniot relation for JB-9014 insensitive explosive [J]. Chinese Journal of High Pressure Physics, 2018, 32(3): 033202. DOI: 10.11858/gywlxb.20170669. [23] 张旭, 池家春, 冯民贤. JB9014钝感炸药冲击绝热线测量 [J]. 高压物理学报, 2001, 15(4): 304–308. DOI: 10.11858/gywlxb.2001.04.011.ZHANG X, CHI J C, FENG M X. Hugoniot relation of JB9014 insensitive high explosive [J]. Chinese Journal of High Pressure Physics, 2001, 15(4): 304–308. DOI: 10.11858/gywlxb.2001.04.011. [24] 刘俊明. JB-9014炸药未反应状态方程研究[D]. 四川绵阳: 中国工程物理研究院, 2018. [25] 徐辉, 孙占峰. 钝感高能炸药JB-9014做功能力的实验研究 [J]. 高压物理学报, 2013, 27(4): 582–586. DOI: 10.11858/gywlxb.2013.04.018.XU H, SUN Z F. An experimental study on the capacity for work of insensitive high explosive [J]. Chinese Journal of High Pressure Physics, 2013, 27(4): 582–586. DOI: 10.11858/gywlxb.2013.04.018. [26] 曾代朋, 陈军, 谭多望. 超压爆轰产物冲击绝热线的实验研究 [J]. 高压物理学报, 2010, 24(1): 76–80. DOI: 10.11858/gywlxb.2010.01.014.ZENG D P, CHEN J, TAN D W. Experiment research on the Hugoniot of JB-9014 explosive over-driven detonation product [J]. Chinese Journal of High Pressure Physics, 2010, 24(1): 76–80. DOI: 10.11858/gywlxb.2010.01.014. [27] TANG P K, ANDERSON W W, FRITZ J N, et al. A study of the overdriven behaviors of PBX 9501 and PBX 9502 [C] // Proceedings of the 11th International Detonation Symposium. Snowmass, CO, 1998. [28] 张琪敏, 张旭, 赵康, 等. TATB基钝感炸药JB-9014的冲击起爆反应增长规律 [J]. 爆炸与冲击, 2019, 39(4): 041405. DOI: 10.11883/bzycj-2018-0050.ZHANG Q M, ZHANG X, ZHAO K, et al. Law of reaction growth of shock initiation on the TATB based insensitive explosive JB-9014 [J]. Explosion and Shock Waves, 2019, 39(4): 041405. DOI: 10.11883/bzycj-2018-0050. [29] 黄文斌, 文尚刚, 谭多望, 等. 在球面散心冲击波作用下JB-9014炸药冲击引爆过程实验研究 [J]. 爆炸与冲击, 2006, 26(4): 373–376. DOI: 10.11883/1001-1455(2006)04-0373-04.HUANG W B, WEN S G, TAN D W, et al. Experimental study on shock initiation process of JB-9014 under loading conditions of divergent spherical shock wave [J]. Explosion and Shock Waves, 2006, 26(4): 373–376. DOI: 10.11883/1001-1455(2006)04-0373-04. [30] 訾攀登, 陈军, 张蓉, 等. 二次压缩条件下JOB-9003炸药特性研究 [J]. 高压物理学报, 2017, 31(2): 155–161. DOI: 10.11858/gywlxb.2017.02.007.ZI P D, CHEN J, ZHANG R, et al. Characteristics of JOB-9003 in double shocks experiments [J]. Chinese Journal of High Pressure Physics, 2017, 31(2): 155–161. DOI: 10.11858/gywlxb.2017.02.007. [31] FERM E N, MORRIS C L, QUINTANA J P, et al. Proton radiography examination of unburned regions in PBX 9502 corner turning experiments [J]. AIP Conference Proceedings, 2002, 620(1): 966–969. DOI: 10.1063/1.1483699. [32] SOUERS P C, ANDRESKI H G, COOK Ⅲ C F, et al. LX-17 corner-turning [J]. Propellants, Explosives, Pyrotechnics, 2004, 29(6): 359–367. DOI: 10.1002/prep.200400067. [33] WHITWORTH N J. CREST modelling of PBX 9502 corner turning experiments at different initial temperatures [J]. Journal of Physics: Conference Series, 2014, 500(5): 052050. DOI: 10.1088/1742-6596/500/5/052050. -