Numerical simulation is widely used in the study of wave fields in various media. One of the methods is to divide the domain of interest into elementary volumes and build a finite-difference scheme for numerical implementation. The work assumes that the domain can have a significant curvature of the surface, therefore, the technology of generating a mesh of curved cubes is used. This mesh provides good consistency between the discrete and physical models of the domain. A parallel algorithm is proposed for the numerical solution of a 3D linear system of elasticity theory, expressed via displacement velocities and stresses, using a curvilinear mesh and an explicit difference scheme based on the Lebedev scheme. The simulation results are presented. The calculations were carried out using the resources of the SSCC SB RAS.