A time-domain approach for the simulation of three-dimensional seismic wave propagation using the scaled boundary finite element method.

A direct time-domain approach to simulate seismic wave propagation in three-dimensional unbounded media is proposed based on the Scaled Boundary Finite Element Method (SBFEM). A domain of interest is commonly partitioned into a far field and a near field. The far field is modelled by the semi-analytical SBFEM satisfying rigorously the radiation conditions at infinity. Separate scaled boundary finite elements are employed to reach a balance between computational efficiency and accuracy. The near field is discretized into arbitrarily-shaped scaled boundary finite elements without the occurrence of hanging nodes. This advantage of the SBFEM in mesh generation is leveraged by incorporating the automatic octree-based meshing technique. By exploiting the geometrical similarity of both bounded and unbounded SBFE subdomains the computational cost is reduced. Inspired by the Domain Reduction Method (DRM), seismic waves are introduced to the system via a single layer of elements in the near field. This formulation of seismic input is mathematically convenient as it avoids the direct participation of the formulation of the far field. The proposed approach is attractive in a reliable simulation of the far field, flexible mesh generation of the near field and simple formulation of the seismic excitations. These merits are demonstrated through numerical simulations of seismic wave propagation in a free field and different examples featuring complex geometries in the near fields. • A direct time-domain approach to simulate seismic wave propagation in three-dimensional media is proposed. • Seismic wave fields in unbounded domains can be rigorously represented by SBFEM. • A bounded domain is be automatically discretized into polyhedra without hanging node problems. • Geometrical similarities of both bounded and unbounded domains are exploited to reduce computational costs. • A simple model for seismic wave inputs is presented base on the domain reduction method. [ABSTRACT FROM AUTHOR]