1D finite element artificial boundary method for transient response of ocean site under obliquely incident earthquake waves. Issue 126 (November 2019)

• Record Type: Journal Article
• Authors: Zhang, Guoliang; Zhao, Mi; Du, Xiuli; Zhang, Xiaoling; Wang, Piguang
• Is Part Of: Soil dynamics and earthquake engineering. Issue 126(2019)
• Publication Date: 2019-11
• Subjects: Ocean site response analysis -- Obliquely incident earthquake -- Poroelastic media -- Finite element method -- Artificial boundary condition; Soil dynamics -- Periodicals; Earthquake engineering -- Periodicals; Sols -- Dynamique -- Périodiques; Génie parasismique -- Périodiques; 624.176205
• DOI: 10.1016/j.soildyn.2019.105787
• Languages: English
• ISSNs: 0267-7261

Abstract

Abstract: Ocean site consisting of seawater-seabed-bedrock system is modeled as horizontally layered fluid-poroelastic-solid media. Earthquake wave is assumed as a plane P or SV wave of oblique incidence from the solid bedrock into the site. A time-domain numerical method called one-dimensional finite element artificial boundary (1D-FE-AB) method is proposed to effectively and accurately calculate the transient response of the ocean site under the obliquely incident earthquake waves. In this method, a horizontally line-shaped AB is firstly introduced in the solid bedrock to divide the whole site into the fluid-poroelastic-solid layer and the remaining homogenous solid half space. The former is semi-discretized by FE method along depth, and the latter with the incident earthquake wave is simulated by an exact AB condition. The resulting semi-discrete system with AB is then transformed into a spatially 1D dynamic equation along depth by using the Snell's law of plane body wave propagation in multilayer media. The 1D equation can be solved by the time integration algorithm such as Newmark method. Numerical examples are finally given to demonstrate the effectiveness of the proposed 1D-FE-AB method. Highlights: A time-domain numerical method is developed to calculate transient response of ocean site under obliquely incident earthquake. The ocean site of seawater-seabed-bedrock system is modeled as horizontally layered fluid-poroelastic-solid half space. The method can give accurate site response in time domain, and avoid analytical calculation for poroelastic wave motion theory.