点火增长反应速率方程在LS-DYNA软件中嵌入及应用

An ignition and growth reaction rate equation with 10 parameters is implanted into LS-DYNA software to precisely describe the shock initiation of explosives and simplify the calibration of reaction rate equation. The reaction rate equation is used to simulate the shock initiation of a DNAN-based melt-cast explosive. The parameters of the reaction rate equation are calibrated by comparing the pressure-time curves in experiment and simulation, and the advantages of the reaction rate equation are analyzed. The results indicate that the reaction rate equation can be easily calibrated due to a few number of parameters. The difference between the simulated and experimental arrival times of shock wave does not exceed 0.2 μs, and thus the reaction rate equation can be used to precisely describe the shock initiation of explosives. The algorithm for implanting the reaction rate equation in LS-DYNA software is efficient because the iterations for every grid does not exceed 3 times.