Abstract: Using logarithmic function as a skeleton curve, the concept of "modified dynamic skeleton curve" and "damping ratio degradation coefficient" is introduced. Considering the influence of test damping on the hysteresis curve, the nonlinear dynamic constitutive of test damping and hysteresis damping is realized. This constitutive method is only suitable for nonlinear seismic response analysis of one-dimensional soil layers. In this paper, based on the logarithmic dynamic skeleton constitutive function relation, the method of solving the logarithmic dynamic skeleton parameters is given, and the time-varying tangent shear modulus of the loading and unloading curve is derived. Based on the operating platform of ABAQUS software and the equivalent shear strain algorithm, a display subroutine module based on a logarithmic dynamic skeleton constitutive model is developed, which is suitable for nonlinear seismic response analysis of two-dimensional and three-dimensional soil layers. The numerical simulation results of EI centro ground motion with different acceleration peaks input from bedrock were compared to verify the log-dynamic skeleton constitutive considering soil damping effect and the rationality and availability of the subprogram developed.