Abstract:The development of viscoelastic boundary theory is systematically reviewed and the establishment of boundary stress equilibrium equations in different dimensional models is expounded. The applicability of viscoelastic theory under one-dimensional, two-dimensional and three-dimensional conditions is verified by using the advanced numerical simulation software ABAQUS. It is found that in the one-dimensional model, the simple viscous boundary condition can meet the required absorption accuracy. In the two-dimensional shear cylindrical wave model, the viscoelastic boundary conditions need to be introduced to achieve the same accuracy requirements. For the two-dimensional vertical incidence model, its characteristic is the unity of the input and absorption boundary. The equivalent nodal force method can be used to ensure that the incident wave is not affected by the viscoelastic boundary, and can still maintain a good absorption effect during reflection. Under three-dimensional condition, the theory of oblique incident wave field decomposition is introduced to angular incident, and the process of automatic application of viscoelastic boundary conditions using Python programming is described in detail. The studied results provide an effective technical means for the dynamic response analysis of engineering structures such as bridges and roadbeds under earthquake.