It is important to investigate the dynamic behaviors of deep rocks near explosion cavity to reveal the mechanisms of deformations and fractures. Some improvements are carried out for Grigorian model with focuses on the dilation effects and the relaxation effects of deep rocks, and the high pressure equations of states with Mie-Grüneisen form are also established. Numerical calculations of free field parameters for deep underground explosions are carried out based on the user subroutines which are compiled by means of the secondary development functions of LS-DYNA9703 D software. The histories of radial stress, radial velocity and radial displacement of rock particles are obtained, and the calculation results are compared with those of U.S. Hardhat nuclear test. It is indicated that the dynamic responses of free field for deep underground explosions are well simulated based on improved Grigorian model, and the calculation results are in good agreement with the data of U.S. Hardhat nuclear test. The peak values of particle velocities are consistent with those of test, but the waveform widths and the rising times are obviously greater than those without dilation effects. The attenuation rates of particle velocities are greater than the calculation results with classic plastic model, and they are consistent with the results of Hardhat nuclear test. The attenuation behaviors and the rising times of stress waves are well shown by introducing dilation effects and relaxation effects into the calculation model. Therefore, the defects of Grigorian model are avoided. It is also indicated that the initial stress has obvious influences on the waveforms of radial stress and the radial displacements of rock particles.
