Earthquake Motion Simulation with Multiscale Finite-Element Analysis on Hybrid Grid

The prediction of strong ground motion with high resolution is a challenging task. In this article, the authors propose a multiscale analysis based on a singular perturbation and a new finite-element method with a hybrid of structured and unstructured elements, for a full 3D numerical simulation of earthquake wave propagation from fault to surface, including soft surface layers. The multiscale analysis refines a solution of lower resolution by considering effects of ground structures on wave propagation, and the hybrid grid radically reduces the amount of numerical computation. Several numerical experiments are carried out to show the validity and usefulness of the finite-element method with the hybrid grids, by comparing the results with those obtained using existing methods. An exemplary problem of wave propagation is solved using the proposed method. The potential usefulness of the proposed method is discussed, with particular attention paid to the accuracy of computing strong ground motion with higher spatial resolution.

[1]  J. Bielak,et al.  Domain Reduction Method for Three-Dimensional Earthquake Modeling in Localized Regions, Part II: Verification and Applications , 2001 .

[2]  Ivo Oprsal,et al.  Three‐dimensional finite difference method and hybrid modeling of earthquake ground motion , 2002 .

[3]  J. Regan,et al.  Numerical modelling of SH L_g waves in and near continentalmargins , 1989 .

[4]  D. Komatitsch,et al.  The spectral element method: An efficient tool to simulate the seismic response of 2D and 3D geological structures , 1998, Bulletin of the Seismological Society of America.

[5]  Paul Christiano,et al.  On the effective seismic input for non-linear soil-structure interaction systems , 1984 .

[6]  N. Kikuchi,et al.  Homogenization theory and digital imaging: A basis for studying the mechanics and design principles of bone tissue , 1994, Biotechnology and bioengineering.

[7]  M. Dumbser,et al.  An arbitrary high-order discontinuous Galerkin method for elastic waves on unstructured meshes - I. The two-dimensional isotropic case with external source terms , 2006 .

[8]  Thomas J. R. Hughes,et al.  Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies , 1985 .

[9]  David R. O'Hallaron,et al.  Earthquake ground motion modeling on parallel computers , 1996, Supercomputing '96.

[10]  J. Bielak,et al.  Unified Symmetric BEM‐FEM for Site Effects on Ground Motion—SH Waves , 1991 .

[11]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[12]  J. Bielak,et al.  Domain Reduction Method for Three-Dimensional Earthquake Modeling in Localized Regions, Part I: Theory , 2003 .

[13]  Eduardo Kausel,et al.  The spring method for embedded foundations , 1978 .

[14]  J. Kristek,et al.  Seismic-Wave Propagation in Viscoelastic Media with Material Discontinuities: A 3D Fourth-Order Staggered-Grid Finite-Difference Modeling , 2003 .

[15]  Peter Moczo,et al.  Hybrid Seismic Modeling Based on Discrete-wave Number and Finite-difference Methods , 1996 .

[16]  Akira Mita,et al.  Dynamic response of embedded foundations: a hybrid approach , 1987 .

[17]  J. Kristek,et al.  3D Heterogeneous Staggered-grid Finite-difference Modeling of Seismic Motion with Volume Harmonic and Arithmetic Averaging of Elastic Moduli and Densities , 2002 .

[18]  Donald R. Smith,et al.  Singular-Perturbation Theory: An Introduction with Applications , 2009 .

[19]  Hiroyuki Fujiwara,et al.  The fast multipole method for solving integral equations ofthree‐dimensional topography and basin problems , 2000 .

[20]  Hiroyuki Fujiwara,et al.  3D finite-difference method using discontinuous grids , 1999, Bulletin of the Seismological Society of America.

[21]  J. Wolf.  Soil-structure-interaction analysis in time domain , 1988 .

[22]  R. Clough,et al.  Dynamics Of Structures , 1975 .

[23]  Arben Pitarka,et al.  3D Elastic Finite-Difference Modeling of Seismic Motion Using Staggered Grids with Nonuniform Spacing , 1999 .

[24]  José M. Carcione,et al.  Hybrid modeling of P-SV seismic motion at inhomogeneous viscoelastic topographic structures , 1997, Bulletin of the Seismological Society of America.

[25]  Sia Nemat-Nasser,et al.  On two micromechanics theories for determining micro–macro relations in heterogeneous solids , 1999 .

[26]  Francisco J. Sánchez-Sesma,et al.  Computing topographical 3D site effects using a fast IBEM/Conjugate gradient approach , 1998, Bulletin of the Seismological Society of America.

[27]  D. Komatitsch,et al.  Simulations of Ground Motion in the Los Angeles Basin Based upon the Spectral-Element Method , 2004 .

[28]  Ezio Faccioli,et al.  2d and 3D elastic wave propagation by a pseudo-spectral domain decomposition method , 1997 .

[29]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[30]  Tsuyoshi Ichimura,et al.  Macro–micro analysis method for wave propagation in stochastic media , 2006 .

[31]  K. Koketsu,et al.  Specific distribution of ground motion during the 1995 Kobe Earthquake and its generation mechanism , 1998 .

[32]  Giuliano F. Panza,et al.  A hybrid method for the estimation of ground motion in sedimentary basins: Quantitative modeling for Mexico city , 1994 .

[33]  M. Nuray Aydınoğlu.  Consistent formulation of direct and substructure methods in nonlinear soil-structure interaction , 1993 .

[34]  Hiroyuki Fujiwara,et al.  Finite-element Simulation of Seismic Ground Motion with a Voxel Mesh , 2004 .