TATB基非均质炸药预冲击减敏的数值模拟

为了对一种TATB基非均质炸药的预冲击起爆现象展开数值模拟研究，将基于冲击温度及压力的AWSD反应率模型耦合进二维结构网格拉氏弹塑性流体力学程序。利用炸药及其产物的冲击雨贡纽实验数据校验了未反应炸药及产物的状态方程参数，通过一维冲击起爆的模拟，标定了反应速率模型参数。模拟了炸药在弱冲击0.45 μs后跟随的强冲击波的二次冲击实验，结果表明，受预压缩区域的炸药反应变慢，到爆轰距离增长了约1 mm，与该炸药二次冲击实验减敏现象相符。模拟拐角效应时，爆轰波经过拐角后，在拐角附近形成稳定的不起爆区域，与主要成分相同的LX-17炸药的拐角效应实验的死区特征相符。数值模拟结果表明，基于冲击温度及压力的AWSD反应率模型可以较好地模拟非均质炸药预冲击减敏问题。

Numerical simulation of pre-shock desensitization in TATB-based heterogeneous explosive

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.