Abstract
Summary Numerical simulation of seismic wave propagation is utilized increasingly, and it can provide realistic and useful information about complex structures. The finite element method (FEM) is a powerful method for seismic forward modeling, and it has been well established in 2-D case with quadrilateral elements due to its computational accuracy and efficiency. However, the implementation of FEM with regular quadrilateral elements can not adapt to mesh complex structures. In this paper, the stiffness matrix of arbitrary quadrilateral meshes is derived, and numerical solution of FEM with arbitrary quadrilateral meshes is given too. In order better to eliminate the boundary reflection, we introduce improved Sarma boundary condition containing a transition region which can suppress the spurious reflections arising on the interfaces between the purely elastic region and the damping region. To improve computational efficiency and reduce memory occupation, we employ compact storage format to store global stiffness matrix. After these, a numerical experiment is presented to compare FEM with arbitrary quadrilateral meshes and FEM with rectangular meshes. The results confirm that the arbitrary quadrilateral meshes technique in this paper is more effective.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call