• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.4
• /
• pp.221-230
• /
• 2021

This paper proposes full waveform inversion (FWI) for estimating the physical properties of a layered half-space. An FWI solution is obtained using a genetic algorithm (GA), which is a well-known global optimization approach. The dynamic responses of a layered half-space subjected to a harmonic vertical disk load are measured and compared with those calculated using the estimated physical properties. The responses are calculated using the thin-layer method, which is accurate and efficient for layered media. Subsequently, a numerical model is constructed for a layered half-space using mid-point integrated finite elements and perfectly matched discrete layers. An objective function of the global optimization problem is defined as the L2-norm of the difference between the observed and estimated responses. A GA is used to minimize the objective function and obtain a solution for the FWI. The accuracy of the proposed approach is applied to various problems involving layered half-spaces. The results verify that the proposed FWI based on a GA is suitable for estimating the material properties of a layered half-space, even when the measured responses include measurement noise.

• Earthquakes and Structures
• /
• v.3 no.3_4
• /
• pp.203-230
• /
• 2012

This paper presents results for impedance (and compliance) functions and input motions of foundations in a layered half-space computed on the basis of a procedure that combines a consistent transmitting boundary with continued-fraction absorbing boundary conditions which are accurate and effective in modeling wave propagation in various unbounded domains. The effects of obliquely incident seismic waves in a layered half-space are taken into account in the formulation of the transmitting boundary. Using the numerical model, impedance (and compliance) functions and input motions of rigid circular foundations on the surface of or embedded in a homogeneous half-space are computed and compared with available published results for verification of the procedure. Extrapolation methods are proposed to improve the performance in the very-low-frequency range and for the static condition. It is concluded from the applications that accurate analysis of foundation dynamics and soil-structure interaction in a layered half-space can be carried out using the enhanced consistent transmitting boundary and the proposed extrapolations.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2001.09a
• /
• pp.75-82
• /
• 2001

A dynamic interaction analysis of an adjacent surface fecundation on a layered half-space is performed in the frequency domain. A semi-analytical approach is employed to reduce the integration range of the wavenumber in the surface fundamental solution for a layered half-space in boundary element (BE) formulations. The present study then adopts a combined boundary and finite element method to analyze the dynamic behavior of a system of flexible surface foundations on an elastic homogeneous and layered half-space. Numerical examples are presented to demonstrate the accuracy of the developed method. The examples show the feasibility of an extended application fur the complicated dynamic interaction of foundations on layered media.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.3
• /
• pp.203-211
• /
• 2023

A Markov chain Monte Carlo (MCMC) simulation is proposed for probabilistic full waveform inversion (FWI) in a layered half-space. Dynamic responses on the half-space surface are estimated using the thin-layer method when a harmonic vertical force is applied. Subsequently, a posterior probability distribution function and the corresponding objective function are formulated to minimize the difference between estimations and observed data as well as that of model parameters from prior information. Based on the gradient of the objective function, a proposal distribution and an acceptance probability for MCMC samples are proposed. The proposed MCMC simulation is applied to several layered half-space examples. It is demonstrated that the proposed MCMC simulation for probabilistic FWI can estimate probabilistic material properties such as the shear-wave velocities of a layered half-space.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2006.03a
• /
• pp.628-636
• /
• 2006

This paper presents 3D infinite elements in Cartesian coordinates for the elastodynamic problem in multi-layered half-space. Five kinds of infinite elements are developed by using approximate expressions of multiple wave components for the wave function in exterior far-field soil region. They are horizontal, horizontal-corner, vertical, vertical-corner and vertical-horizontal-corner elements. The elements can be used for the multi-wave propagating problem. Numerical example analyses are presented for rigid disk, square footings and embedded footing on homogeneous and layered half-space. The numerical results obtained show the effectiveness of the proposed infinite elements.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.1
• /
• pp.39-50
• /
• 2007

This paper presents 3D infinite elements for the elastodynamic problem with multi-layered half-space. Five different types of infinite elements are formulated by using approximate expressions of multiple wave components for the wave function in multi-layered soil media. They are horizontal, horizontal-corner, vortical, vertical-corner and vertical-horizontal-comer infinite elements. The elements can effectively be used for simulating wane radiation problems with multiple wave components. Numerical example analyses are presented for rigid disk, square footings and embedded footing on homogeneous and layered half-space. The numerical results show the effectiveness of the proposed infinite elements.

• Earthquakes and Structures
• /
• v.9 no.1
• /
• pp.173-194
• /
• 2015

Site response analysis is an important topic in earthquake engineering. A time-domain numerical method called as one-dimensional (1D) finite element artificial boundary method is proposed to simulate the homogeneous plane elastic wave propagation in a layered half space subjected to the obliquely incident plane body wave. In this method, an exact artificial boundary condition combining the absorbing boundary condition with the inputting boundary condition is developed to model the wave absorption and input effects of the truncated half space under layer system. The spatially two-dimensional (2D) problem consisting of the layer system with the artificial boundary condition is transformed equivalently into a 1D one along the vertical direction according to Snell's law. The resulting 1D problem is solved by the finite element method with a new explicit time integration algorithm. The 1D finite element artificial boundary method is verified by analyzing two engineering sites in time domain and by comparing with the frequency-domain transfer matrix method with fast Fourier transform.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.10a
• /
• pp.317-324
• /
• 1999

Dynamic interaction analysis of surface structure on layered half-space is performed in frequency domain under incident wave excitation. This present study adopts a coupling method that combines the finite element(FE) for the flexible structures and boundary element(BE) for the layered half-space. A semi-analytical approach is employed to reduce the integration range of wavenumbers in the BE formula. For the incident wave input, the response is decomposed and formulated after the impedance matrix for the structure system. Numerical examples are presented to demonstrate the accuracy of the method. The examples show the feasibility of an extended application to the complicated dynamic analysis of structures on layered media under incident wave excitation.

• Earthquakes and Structures
• /
• v.13 no.6
• /
• pp.573-588
• /
• 2017

The earthquake input is required when the soil-structure interaction (SSI) analysis is performed by the direct finite element method. In this paper, the earthquake is considered as the obliquely incident plane body wave arising from the truncated linearly elastic layered half space. An earthquake input method is developed for the time-domain three-dimensional SSI analysis. It consists of a new site response analysis method for free field and the viscous-spring artificial boundary condition for scattered field. The proposed earthquake input method can be implemented in the process of building finite element model of commercial software. It can result in the highly accurate solution by using a relatively small SSI model. The initial condition is considered for the nonlinear SSI analysis. The Daikai subway station is analyzed as an example. The effectiveness of the proposed earthquake input method is verified. The effect of the obliquely incident earthquake is studied.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.10a
• /
• pp.257-266
• /
• 1999

This paper presents a method of seismic analysis for a cylindrical liquid storage structure on/in horizontally layered half.space considering the effects of the interior fluid and exterior soil medium in the frequency domain. To capture the essence of fluid-structure-soil interaction effects effectively, a mixed finite element with two-field (u, p) approximation is employed to model the compressive inviscid fluid, while the structure and soil medium are presented by the 3-D axisymmetric finite elements and dynamic infinite elements. The present FE-based method can be applied to the system with complex geometry of fluid region as well as with inhomogeneous near-field soil medium, since it can directly model both the fluid and the soil. For the purpose of verification, dominant peak frequencies in transfer functions for horizontal motions of cylindrical fluid storage tanks with rigid massless foundation on a homogeneous viscoelastic half.space are compared with those by two different added mass approaches for the fluid motion. The comparison indicates that the Present FE-based methodology gives accurate solution for the fluid-structure-soil interaction problem. Finally, as a demonstration of versatility of the present study, a seismic analysis for a real-scale LNG storage tank embedded in layered half.space is carried out, and its member forces along the height of the structure are compared with those by an added mass approach developed by the present writers.