Abstract: |
The self similar solution of the strong explosion of spherically symmetrical charge is often used to give the initial wavefront parameters of a point source explosion in the air. The accuracy of the wavefront parameters seriously affects the subsequent shock wave propagation and the calculation of the impact load of the equipment and structure. In order to improve the calculation accuracy of the self similar solution for the strong explosion of spherically symmetrical charges, this paper uses the adaptive step Runge Kutta method instead of the original fixed step fourth order Runge Kutta method. The algorithm improves the accuracy of the self similar of the strong air explosion to the fifth order, and realizes the artificial control of the solution accuracy. In this paper, the improved algorithm is realized by Matlab software, and the self similar solution of spherical symmetrical charge in the air is solved. By comparing the original algorithm, the improved algorithm, and the reference value, it is considered that the improved algorithm can control the calculation error within the specified range. This research can provide theoretical basis and data support for the high precision solution of impact load on structures. |