A modified domain reduction method for numerical simulation of wave propagation in localized regions

A modified domain reduction method (MDRM) that introduces damping terms to the original DRM is presented in this paper. To verify the proposed MDRM and compare the computational accuracy of these two methods, a numerical test is designed. The numerical results of the MDRM and DRM are compared using an extended meshed model. The results show that the MDRM significantly improved the computational accuracy of the DRM. Then, the MDRM is compared with two existing conventional methods, namely Liao’s transmitting boundary and viscous-spring boundary with Liu’s method. The MDRM shows its great advancement in computational accuracy, stability and range of applications. This paper also discusses the influence of boundary location on computational accuracy. It can be concluded that smaller models tend to have larger errors. By introducing two dimensionless parameters, φ1 and φ2, the rational distance between the observation point and the MDRM boundary is suggested. When φ1>2 or φ2>13, the relative PGA error can be limited to 5%. In practice, the appropriate model size can be chosen based on these two parameters to achieve desired computational accuracy.

[1]  Thomas J. R. Hughes,et al.  Improved numerical dissipation for time integration algorithms in structural dynamics , 1977 .

[2]  Liu Jing-bo,et al.  A direct method for analysis of dynamic soil-structure interaction based on interface idea , 1998 .

[3]  E. Lindman.  “Free-space” boundary conditions for the time dependent wave equation , 1975 .

[4]  Roberto Paolucci,et al.  Comparison of 3D, 2D and 1D numerical approaches to predict long period earthquake ground motion in the Gubbio plain, Central Italy , 2011 .

[5]  Zou Jing-xiang.  A high -frequency instability mechanism in numerical realization of multi-transmitting formula , 2002 .

[6]  Boris Jeremić,et al.  Time domain simulation of soil–foundation–structure interaction in non‐uniform soils , 2009 .

[7]  D. Potts,et al.  An assessment of the domain reduction method as an advanced boundary condition and some pitfalls in the use of conventional absorbing boundaries , 2009 .

[8]  Liu Jin,et al.  Seismic stability of jointed rock slopes under obliquely incident earthquake waves , 2018, Earthquake Engineering and Engineering Vibration.

[9]  关慧敏,et al.  A METHOD FOR THE STABILITY ANALYSIS OF LOCAL ARTIFICIAL BOUNDARIES , 1996 .

[10]  Zhenpeng Liao,et al.  A TRANSMITTING BOUNDARY FOR TRANSIENT WAVE ANALYSES , 1984 .

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

[12]  Eduardo Kausel,et al.  Local transmitting boundaries , 1988 .

[13]  W. B. Joyner,et al.  Calculation of nonlinear ground response in earthquakes , 1975 .

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

[15]  I. M. Idriss,et al.  Seismic Response of Horizontal Soil Layers , 1968 .

[16]  John Lysmer,et al.  Analytical Procedures in Soil Dynamics , 1978 .

[17]  Z. Liao,et al.  Numerical instabilities of a local transmitting boundary , 1992 .

[18]  L. Kellezi.  Local transmitting boundaries for transient elastic analysis , 2000 .

[19]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[20]  杜义欣,et al.  3D viscous-spring artificial boundary in time domain , 2006 .

[21]  Lidija Zdravković,et al.  Seismic response and interaction of complex soil retaining systems , 2012 .

[22]  Haibing Chen,et al.  Accuracy of three-dimensional seismic ground response analysis in time domain using nonlinear numerical simulations , 2017, Earthquake Engineering and Engineering Vibration.

[23]  J. Lysmer,et al.  Finite Dynamic Model for Infinite Media , 1969 .

[24]  Zhinan Xie,et al.  Analysis of high-frequency local coupling instability induced by multi-transmitting formula–P-SV wave simulation in a 2D waveguide , 2017, Earthquake Engineering and Engineering Vibration.

[25]  Lidija Zdravković,et al.  The domain reduction method for dynamic coupled consolidation problems in geotechnical engineering , 2008 .

[26]  Jerome Solberg,et al.  Nonlinear time-domain soil–structure interaction analysis of embedded reactor structures subjected to earthquake loads , 2016 .

[27]  George P. Mavroeidis,et al.  A Mathematical Representation of Near-Fault Ground Motions , 2003 .

[28]  Mark Randolph,et al.  Axisymmetric Time‐Domain Transmitting Boundaries , 1994 .

[29]  Zhen‐Peng Liao,et al.  Stable Implementation of Transmitting Boundary in Numerical Simulation of Wave Motion , 2002 .

[30]  Loukas F. Kallivokas,et al.  Seismic wave amplification by topographic features: A parametric study , 2017 .

[31]  Jacobo Bielak,et al.  Coupled Soil-Structure Interaction Effects of Building Clusters during Earthquakes , 2015 .

[32]  B. Jeremić,et al.  Seismic behavior of NPP structures subjected to realistic 3D, inclined seismic motions, in variable layered soil/rock, on surface or embedded foundations , 2013 .

[33]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[34]  Z. Liao,et al.  A transmitting boundary for the numerical simulation of elastic wave propagation , 1984 .

[35]  Roberto Paolucci,et al.  Seismic analysis of deep tunnels in near fault conditions: a case study in Southern Italy , 2011 .

[36]  Ronald Y. S. Pak,et al.  PML solution of longitudinal wave propagation in heterogeneous media , 2016, Earthquake Engineering and Engineering Vibration.