Virtual topography: A fictitious domain approach for analyzing free‐surface irregularities in large‐scale earthquake ground motion simulation

SUMMARY This paper presents a numerical scheme based on a fictitious domain framework for the numerical modeling of earthquake-induced ground motion in the presence of realistic surface topography of the Earth's crust. We show that by adopting a non-conforming octree-based meshing approach associated with a virtual representation of the surficial irregularity, we can obtain accurate representations of ground motion. From the computational point of view, our methodology proves to be also efficient, and more importantly, it allows us to preserve the salient features of multi-resolution cubic-shaped finite elements for wave propagation applications. We implemented the non-conforming meshing scheme for the treatment of realistic topographies into Hercules, the octree-based finite-element earthquake simulator developed by the Quake Group at Carnegie Mellon University. We tested the benefits of the strategy by benchmarking its results against reference examples, and by means of numerical error estimate analyses. Our qualitative and quantitative comparisons showed a close agreement between our numerical results and the reference results, and also, that the order of convergence of the displacement field is preserved in the presence of surface topography. Moreover, this performance was obtained by using the same mesh refinement techniques with cubic elements as in traditional flat free-surface simulations. Copyright © 2014 John Wiley & Sons, Ltd.

[1]  Joseph E. Bishop,et al.  Rapid stress analysis of geometrically complex domains using implicit meshing , 2003 .

[2]  Gideon Juve,et al.  The ShakeOut earthquake scenario: Verification of three simulation sets , 2010 .

[3]  Hubert Maigre,et al.  An explicit dynamics extended finite element method. Part 1: Mass lumping for arbitrary enrichment functions , 2009 .

[4]  Peter Wriggers,et al.  Aspects of the computational testing of the mechanical properties of microheterogeneous material samples , 2001 .

[5]  Julien Réthoré,et al.  Efficient explicit time stepping for the eXtended Finite Element Method (X‐FEM) , 2006 .

[6]  T. Belytschko,et al.  MODELING HOLES AND INCLUSIONS BY LEVEL SETS IN THE EXTENDED FINITE-ELEMENT METHOD , 2001 .

[7]  David R. O'Hallaron,et al.  Large-scale simulation of elastic wave propagation in heterogeneous media on parallel computers , 1998 .

[8]  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.

[9]  T. Tu,et al.  From Mesh Generation to Scientific Visualization: An End-to-End Approach to Parallel Supercomputing , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[10]  Tsuyoshi Ichimura,et al.  A hybrid multiresolution meshing technique for finite element three‐dimensional earthquake ground motion modelling in basins including topography , 2009 .

[11]  T. Belytschko,et al.  The extended/generalized finite element method: An overview of the method and its applications , 2010 .

[12]  Ernst Rank,et al.  The finite cell method for three-dimensional problems of solid mechanics , 2008 .

[13]  Ted Belytschko,et al.  Elastic crack growth in finite elements with minimal remeshing , 1999 .

[14]  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.

[15]  Julien Yvonnet,et al.  A multiple level set approach to prevent numerical artefacts in complex microstructures with nearby inclusions within XFEM , 2011 .

[16]  Morgan T. Page,et al.  Effects of Large-Scale Surface Topography on Ground Motions, as Demonstrated by a Study of the San Gabriel Mountains, Los Angeles, California , 2007 .

[17]  Haim Waisman,et al.  A spline‐based enrichment function for arbitrary inclusions in extended finite element method with applications to finite deformations , 2013 .

[18]  Ted Belytschko,et al.  A finite element method for crack growth without remeshing , 1999 .

[19]  Francisco J. Sánchez-Sesma,et al.  DIFFRACTION OF P, SV AND RAYLEIGH WAVES BY TOPOGRAPHIC FEATURES: A BOUNDARY INTEGRAL FORMULATION , 1991 .

[20]  I. Babuska,et al.  The design and analysis of the Generalized Finite Element Method , 2000 .

[21]  Pierre-Yves Bard,et al.  THE EFFECT OF TOPOGRAPHY ON EARTHQUAKE GROUND MOTION: A REVIEW AND NEW RESULTS , 1988 .

[22]  Tsuyoshi Ichimura,et al.  Earthquake Motion Simulation with Multiscale Finite-Element Analysis on Hybrid Grid , 2007 .

[23]  G. Allaire,et al.  Structural optimization using sensitivity analysis and a level-set method , 2004 .

[24]  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 .

[25]  I. Babuska,et al.  The Partition of Unity Method , 1997 .

[26]  Jean-Pierre Vilotte,et al.  RegSEM: a versatile code based on the spectral element method to compute seismic wave propagation at the regional scale , 2012 .

[27]  Mark A Fleming,et al.  Smoothing and accelerated computations in the element free Galerkin method , 1996 .

[28]  Renaud Masson,et al.  Modeling the effective elastic behavior of composites: a mixed Finite Element and homogenisation approach , 2008 .

[29]  Ernst Rank,et al.  Finite cell method , 2007 .

[30]  I. Babuska,et al.  Generalized finite element method using mesh-based handbooks: application to problems in domains with many voids , 2003 .

[31]  Hiroyuki Goto,et al.  Galerkin boundary integral equation method for spontaneous rupture propagation problems: SH‐case , 2007 .

[32]  T. Strouboulis,et al.  The generalized finite element method: an example of its implementation and illustration of its performance , 2000 .

[33]  Pierre-Yves Bard,et al.  Diffracted waves and displacement field over two-dimensional elevated topographies , 1982 .

[34]  T. Belytschko,et al.  Arbitrary branched and intersecting cracks with the eXtended Finite Element Method , 2000 .

[35]  Jean-François Remacle,et al.  A computational approach to handle complex microstructure geometries , 2003 .

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

[37]  Colby C. Swan,et al.  INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING , 2022 .

[38]  Michel Bouchon,et al.  Effect of topography on surface motion , 1973, Bulletin of the Seismological Society of America.

[39]  N. Moës,et al.  X-FEM explicit dynamics for constant strain elements to alleviate mesh constraints on internal or external boundaries , 2008 .

[40]  Antonio Huerta,et al.  Imposing essential boundary conditions in mesh-free methods , 2004 .

[41]  Jinkyo F. Choo,et al.  Dynamic soil-structure interaction analysis for complex soil profiles using unaligned mesh generation and nonlinear modeling approach , 2013 .

[42]  Eduardo Reinoso,et al.  Three-dimensional scattering of seismic waves from topographical structures , 1997 .