Abstract
Recent ground motion simulations based on the spectral element method have mainly focused on source-internal problems. Reports on seismic responses due to plane wave incidence are rare, even though these are frequently encountered in site effect analysis and soil-structure interaction simulations. In this study, a seismic motion simulation method for 2D complex sites under plane SH wave incidence is proposed based on high-precision spectral element method (SEM) and multi-transmitting formula (MTF). The accuracy and stability of this method are validated through numerous numerical examples. Comparisons with analytical solutions and extended results demonstrate the high simulation accuracies for the ground motions in 2D complex site models when SH plane waves are vertically incident from the bottom. Compared with MTF in the traditional finite-element method (FEM), MTF integrated into SEM demonstrates better performance for both high- and low-frequency stabilities. First, until to the third-order MTF, this method shows a fairly good ability to eliminate high-frequency instability. This instability only occurs when the artificial wavespeed, ca, adopted in the MTF is larger than about five times the physical velocity, cs. For the most common case of ca = cs, long-term stable results can be achieved without any special treatments. In addition, under the condition of no interventions, the low-frequency drift instability for MTF in SEM is considerably lessened compared with that of MTF in FEM, and the drift instability can also be eliminated by a smaller perturbation parameter than in the latter. The method presented in this study combines the advantages of SEM and MTF, and exhibits better stability performance, indicating good application prospects for wave motion simulations of plane wave incidence problems.
Similar content being viewed by others
References
B’erenger JP (1994) A perfectly matched layer for the absorption of electromagnetic waves. J Comput Phys 114:185–200
Basabe D, Jonás D, Sen MK (2007) Grid dispersion and stability criteria of some common finite-element methods for acoustic and elastic wave equations. Geophysics 72(6):81–95
Chew WC, Wagner RL (1992) A modified form of Liao’s absorbing boundary condition // IEEE Antennas & Propagation Society International Symposium. IEEE: 536-539.
Clayton R, Engquist B (1977) Absorbing bonadary conditions for acoustic and elastic wave equations. Bull Seismol Soc Am 67(6):1529–1540
Dai ZJ, Li XJ, Hou CL (2015) A combination usage of transmitting formula and spectral element method and the study of its stability. Eng Mechan 32(11):40–50 (in Chinese)
Deeks AJ, Randolph MF (1994) Axisymmetric time-domain transmitting boundaries. J Eng Mech 120(1):25–42
Higdon RL (1987) Numerical absorbing boundary conditions for the wave equations. Math Comput 49(179):65–90
Huang JJ (2018) An incrementation-adaptive multi-transmitting boundary for seismic fracture analysis of concrete gravity dams. Soil Dyn Earthq Eng 110:145–158
Komatitsch D, Tromp J (2003) A perfectly matched layer absorbing boundary condition for the second-order seismic wave equation. Geophys J Int 154:146–153
Komatitsch D, Liu Q, Tromp J et al (2004) Simulations of ground motion in the Los Angeles basin based upon the spectral-element method. Bull Seismol Soc Am 94(1):187–206
Lee SJ, Chen HW, Liu QY et al (2008) Three-dimensional simulations of seismic-wave propagation in the Taipei basin with realistic topography based upon the spectral-element method. Bull Seismol Soc Am 98(1):253–264
Lee SJ, Chan YC, Komatitsch D, Huang BS, Tromp J (2009) Effects of realistic surface topography on seismic ground motion in the Yangminshan region of Taiwan based upon the spectral-element method and LiDAR DTM. Bull Seismol Soc Am 99(2A):681–693
Li XJ, Liao ZP (1996) The drift instability of local transmitting boundary in time domain. Acta Mech Sinica 28(5):627–632 (in Chinese)
Liao ZP (2002) Introduction to wave motion theories for engineering, 2nd edn. Science Press, Beijing, pp 236–237 (in Chinese)
Liao ZP, Liu JB (1992) Numerical instabilities of a local transmitting boundary. Earthq Eng Struct Dyn 21:65–77
Liao ZP, Wong HL (1984) A transmitting boundary for the numerical simulation of elastic wave propagation. Soil Dyn Earthq Eng 3:174–183
Liao ZP, Wong HL, Yang BP et al (1984) A transmitting boundary for transient wave analyses. Sci Sin (Ser A) 27(10):1063–1076
Liao ZP, Zhou ZH, Zhang YH (2002) Stable implementation of transmitting boundary in numerical simulation of wave motion. Chin J Geophys 45(4):554–568
Liu JB, Du YX, Du XL et al (2006) 3D viscous-spring artificial boundary in time domain. Earthq Eng Eng Vib 5(1):93–102
Lysmer J, Kuhlemeyer RL (1969) Finite dynamic model for infinite media. J Eng Mech Div 95(4):859–878
Martin R, Komatitsch D, Gedney SD (2008) A variational formulation of a stabilized unsplit convolutional perfectly matched layer for the isotropic or anisotropic seismic wave equation. Comput Model Eng Sci 37(3):274–304
Meza-Fajardo KC, Papageorgiou AS (2008) A nonconvolutional, split-field, perfectly matched layer for wave propagation in isotropic and anisotropic elastic media: stability analysis. Bull Seismol Soc Am 98(4):1811–1836
Patera AT (1984) A spectral element method for fluid dynamics: laminar flow in a channel expansion. J Comput Phys 54(3):468–488
Pilz M, Parolai S, Stupazzini M, Paolucci R, Zschau J (2011) Modelling basin effects on earthquake ground motion in the Santiago de Chile basin by a spectral element code. Geophys J Int 187(2):929–945
Pozrikidis C (2014) Introduction to finite and spectral element methods using MATLAB, 2nd edn. CRC Press, Baca Raton, pp 700–707
Seriani G (2004) Double-grid Chebyshev spectral elements for acoustic wave modeling. Wave Motion 39:351–360
Seriani G, Priolo E, Carcione JM, et al. (1992) High-order spectral element method for elastic wave modeling. Expanded Abstracts of 62nd SEG Annual Int Mtg 1285-1288
Shao XM, Lan ZL (1995) Numerical simulation of the seismic wave propagation in inhomogeneous isotropic elastic media. Chin J Geophys S1:39–55 (in Chinese)
Shi L, Wang P, Cai YQ et al (2016) Multi-transmitting formula for finite element modeling of wave propagation in a saturated poroelastic medium. Soil Dyn Earthq Eng 80:11–24
Tang H, Rong MS (2020) An improved wave motion input method for application of multi-transmitting boundary. Wave Motion 97:102600
Trifunac MD (1971) Surface motion of a semi-cylindrical alluvial valley for incident plane SH waves. Bull Seismol Soc Am 61(6):1755–1770
Wolf JP (1986) A comparison of time-domain transmitting boundaries. Earthq Eng Struct Dyn 14:655–673
Xie ZN, Liao ZP (2008) A note for the mechanism of high-frequency instability induced by aborbing boundary conditions. Acta Seimologica Sinica 3:302–306+328 (in Chinese)
Xie ZN, Zhang XB (2017) Analysis of high-frequency local coupling instability induced by multi-transmitting formula: P-SV wave simulation in a 2D waveguide. Earthq Eng Eng Vib 16(1):1–10
Xing HJ, Li HJ (2017) Implementation of Multi-transmitting boundary condition for wave motion simulation by spectral element method: 2-D case. Acta Mech Sinica 49(4):894–906 (in Chinese)
Yu YY, Ding HP, Liu QF (2017a) Three-dimensional simulations of strong ground motion in the Sichuan basin during the Wenchuan earthquake. Bull Earthq Eng 15:4661–4679
Yu YY, Ding HP, Liu QF (2017b) Integration of transmitting boundary and spectral element method and improvement on the accuracy of wave motion simulation. J Vibrat Shock 36(2):13–22 (in Chinese)
Zhang XB (2012) Several researches on numerical wave simulation in the interior computational domain. Institute of Engineering Mechanics, China Earthquake Administration, Harbin, pp 10–11 (in Chinese)
Zhang L, Yu T (2012) A Method of improving the stability of liao’s higher-order absorbing boundary condition. Prog Electromagn Research M 27:167–178
Zhang XB, Liao ZP, Xie ZN (2015) Mechanism of high frequency coupling instability and stable implementation for transmitting boundary: SH wave motion. Chin J Geophys 58(10):3639–3648 (in Chinese)
Acknowledgements
The authors would like to appreciate the constructive comments and suggestions of two anonymous reviewers for improving the manuscript.
Funding
This work was supported by the National Natural Science Foundation of China under Grant No. 51808371.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Yu, Y., Ding, H. & Zhang, X. Simulations of ground motions under plane wave incidence in 2D complex site based on the spectral element method (SEM) and multi-transmitting formula (MTF): SH problem. J Seismol 25, 967–985 (2021). https://doi.org/10.1007/s10950-021-09995-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10950-021-09995-y