摘要

This study proposes a meshless symplectic algorithm for two-dimensional coupled elastic waves, which satisfy the equations of linear isotropic elasticity on a rectangular domain with appropriate initial and boundary conditions. This method is based on a combination of radial basis functions (RBFs) collocation and geometric integrators. The method involves two steps, first the spatial discretization by using RBFs interpolation transforms the seismic waves into a finite-dimensional Hamiltonian system, then evolves the semi-discretized system by using symplectic integrators. Furthermore, the conservations of energy and symplectic structure are investigated. The truncation error and global error of the method are estimated. Numerical experiments show that the scheme is efficient, easy to implement with the scattered knots and possesses a good behavior in the long time integration.