Abstract: The logarithmic dynamic skeleton constitutive model uses the logarithmic function as the skeleton curve and introduces the concepts of "modified dynamic skeleton curve" and "damping ratio degradation coefficient". Considering the influences of the test damping on the hysteresis curve, the nonlinear dynamic constitutive model for the test damping and hysteresis damping is realized. This constitutive method is only suitable for the nonlinear seismic response analysis of one-dimensional soil layers. In this study, 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 the 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 the input EI centro ground motion with different acceleration peaks from bedrock are compared to verify the logarithmic dynamic skeleton constitutive model considering damping effects of soils and the rationality and availability of the developed subprogram.