The calculation of buried large-sized structures for specified seismic impacts should take into account the field of gravity, contact interaction with the ground. The representation of the boundless half-space of the soil by a limited finite element domain requires taking measures to exclude the influence of boundary effects on the solution. The article proposes and analyzes a set of measures to eliminate edge effects: an algorithm for restoring kinematic load according to an experimental seismogram, the use of superimposed grids to separate the waves incident and radiated from the structure, and quasi-dimensional grids with viscosity to dampen waves at infinity. The algorithm for restoring the seismogram is based on comparing the discrete analogue of the experimental seismogram with the results of a one-dimensional problem of the test pulse run in a ground environment and considers the re-reflection of elastic waves in a multilayer soil with a horizontal upper surface and the boundaries of the media. The method of superimposed grids allows to bring the source of seismic impact closer to the foundation of the structure. The secondary interference waves emitted by the structure freely exit to the technical subdomain with a sparse grid, where calculations are carried out using linear viscosity to dampen the wave and noise associated with the numerical dispersion of the solution on the inhomogeneities of the grid. The proposed approach using superimposed grids makes it possible to reduce the time of numerical calculations by orders, since only in the subregion adjacent to the building small grids are set and nonlinear effects (contact interaction, elastic-plastic deformation) are taken into account. The use of quasi-uniform grids with the introduction of linear viscosity will reduce computational costs by two orders in two-dimensional problems and three orders in three-dimensional ones.