Time history analysis is widely applied in seismic analysis of underground structures

and the simulation of soil ground under dynamic loads has a significant influence on the result.Formulations of Rayleigh damping and Hardin-Drnevich model are introduced in the paper.The repeatable boundary condition is validated through comparison with far free boundary.On the basis of the repeatable boundary condition

the numerical simulation results with Shake91 model

Rayleigh damping model

Mohr coulomb model and Hardin-Drnevich model under 2 Hz and 8 Hz sine waves are compared from two aspects

that is

the ground surface acceleration responses and the peak acceleration profiles along depth.It is concluded that the Mohr coulomb model has limited applicability in dynamic time history analysis due to its constitutive model theory.It is recommended that the Rayleigh damping parameters be calculated with the fundamental frequency of soil ground and the natural frequency that approaching the dominant frequency of the imposed seismic wave.The viscoelastic model（e.g.Hardin-Drnevich model）is more suitable for seismic analysis of underground structures than the equivalent linear method

namely Shake91 model and Rayleigh damping model.