• Structural Engineering and Mechanics
• /
• v.15 no.6
• /
• pp.717-733
• /
• 2003

A time domain method is presented for soil-structure interaction analysis under seismic excitations. It is based on the finite element formulation incorporating infinite elements for the far field soil region. Equivalent earthquake input forces are calculated based on the free field responses along the interface between the near and far field soil regions utilizing the fixed exterior boundary method in the frequency domain. Then, the input forces are transformed into the time domain by using inverse Fourier transform. The dynamic stiffness matrices of the far field soil region formulated using the analytical frequency-dependent infinite elements in the frequency domain can be easily transformed into the corresponding matrices in the time domain. Hence, the response can be analytically computed in the time domain. A recursive procedure is proposed to compute the interaction forces along the interface and the responses of the soil-structure system in the time domain. Earthquake response analyses have been carried out on a multi-layered half-space and a tunnel embedded in a layered half-space with the assumption of the linearity of the near and far field soil region, and results are compared with those obtained by the conventional method in the frequency domain.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.2
• /
• pp.231-244
• /
• 2000

This paper presents a dynamic analysis technique for a 2-D soil-structure interaction problem in the frequency domain, which can directly be applied as an analysis tool for seismic response analyses of underground structures, tunnels, embankments, and so on. In this method, the structure and near-field soil is modeled by the standard finite elements, while the unbounded far-field soil is represented using the dynamic infinite elements in the frequency domain. The earthquake-input motion is regarded as traveling P and SV waves which are incident vertically from the far-field of underlying half-space to the near-field of layered medium. The equivalent earthquake forces are then calculated utilizing so-called fixed-exterior-boundary-method and the free-field responses including displacements and tractions. For the verification of the present study, seismic response analyses are carried out for a multi-layered half-space free-field soil medium and a cylindrical cavity embedded in a homogeneous half-space. Comparisons of the present results with solutions by other approaches indicate that the proposed methodology gives accurate estimates. Finally, an application example of seismic response analysis for a subway station is presented, which demonstrates the applicability of the present study.

• The Journal of the Acoustical Society of Korea
• /
• v.27 no.2
• /
• pp.92-101
• /
• 2008

To dissolve the multi-layered model problems, and to complement 2-layered model's simplicity, assumed fluid-fluid-solid 3-layered model. Generally it is known that if the sediment thickness is more than 10 wavelength, the half space's influence to the in-water acoustic field could be disregarded. By tracking the maximum correlation coefficient of calculated results and experimental ones we confirmed that the requirement could be more realized. To calculate the maximum correlation coefficient we used single sensor transmission loss. On the assumption that the sediment sound velocity was 1813 m/s and frequency range 50 kHz to 120 kHz, the conversion condition was from 2.5 to 7.7 wavelength.

• Geomechanics and Engineering
• /
• v.16 no.1
• /
• pp.11-19
• /
• 2018

Propagation of the generalized Rayleigh waves in an elastic half-space covered by an elastic layer for different initial stress combinations and imperfect contact conditions is investigated. Three-dimensional linearized theory of elastic waves in initially stressed bodies in plane-strain state is employed, the corresponding dispersion equation is derived and an algorithm is developed for numerical solution to this equation. Numerical results on the influence of the initial stress patterns and on the influence of the contact conditions are presented and discussed. The case where the external forces are "follower forces" is considered as well. These investigations provide some theoretical foundations for the study of the near-surface waves propagating in layered mechanical systems and can be successfully used for estimation of the degree of the bonded defects between layers, fault characteristics and study of the behavior of seismic surface waves propagating under the bottom of the oceans.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.10a
• /
• pp.202-209
• /
• 1996

In this paper, two dimensional vertical and comer infinite elements which can include multiple wave components to model underlying half space are developed. These elements are natural and economical to model underlying stiff half space or rock. To verify the behavior of these infinite elements, vertical, horizontal, and rocking compliances of a rigid strip foundation on a viscoelastic soil profile are analyzed and compared with those of Tzong and Penzien who used the boundary solution method. Good agreements are noticed between the two methods. The influence of material properties like Poisson's ratio, material damping, and stiffness ratio of layers as well as the influence of geometrical properties such as layer thicknesses and depth of foundation embedment are studied. Example analysis is carried out for the shaking table which is located in KIMM(Korea Institute of Machinery and Materials), and the vertical and horizontal displacements of the analysis are compared with the measured, and show good results and demonstrate the efficiency of the proposed method.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1999.10a
• /
• pp.83-90
• /
• 1999

This paper presents a method of seismic analysis for a cylindrical liquid storage structure on horizontally layered half-space considering the effects of the interior fluid and exterior soil medium in the frequency domain. the horizontal and rocking motions of the structures are included in this study. The fluid motion is expressed in terms of analytical velocity potential function which can be obtained by solving the boundary value problem including the sloshing behavior of the fluid as well as deformed configuration of the structure. The effect of the fluid is included in the equation of motion as the impulsive added mass and a frequency-dependent convective added mass along the nodes on the wetted boundary with structure. The soil medium is presented using the 3-D axisymmetric finite elements and dynamic infinite elements. The present method can be applied to the structures embedded in ground as well as on ground since it models the soil medium directly as well as the structure. For the purpose of vertification dynamci characteristics of a tank on homogeneous half-space is analyzed. Comparison of the present results with those by others shows good agreement.

• Earthquakes and Structures
• /
• v.26 no.2
• /
• pp.163-174
• /
• 2024

Dynamic equilibrium equations for finite element analysis were derived for the free field one-dimensional shear wave propagation through the horizontally layered soil deposits with the elastic half-space. We expressed Rayleigh's viscous damping consisting of mass and stiffness proportional terms. We considered two cases where damping matrices are defined in the total and relative displacement fields. Two forms of equilibrium equations are presented; one in terms of total motions and the other in terms of relative motions. To evaluate the performance of new equilibrium equations, we conducted two sets of site response analyses and directly compared them with the exact closed-form frequency domain solution. Results show that the base shear force as earthquake load represents the simpler form of equilibrium equation to be used for the finite element method. Conventional finite element procedure using base acceleration as earthquake load predicts exact solution reasonably well even in soil deposits with unrealistically high damping.

• Geophysics and Geophysical Exploration
• /
• v.24 no.3
• /
• pp.98-112
• /
• 2021

We followed and extended the algorithm originally made by Das (1995) to calculate the electric potential and field induced by electric current in arbitrary anisotropic layered structure. We confirmed all the theoretical contents and coded the corresponding program to acquire the electric potential and field. Further we extended to forward estimation of apparent resistivity to be attained by electrical resistivity survey on anisotropic layered structure with differing the electrode spacing and azimuth of anisotropy. The effects of anisotropy were reviewed by considering some examples.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1989.04a
• /
• pp.22-27
• /
• 1989

This paper presents a study of soil-structure interaction problems using infinite elements. The infinite elements are formulated for homogeneous and layered soil media, based on approximate expressions for three components of propagating waves, namely Rayleigh, compressive and shear waves. The integration scheme which was proposed for problems with single wave component by Zienkiewicz is expanded to the multi-wave problem. Verifications are carried out on rigid circular footings which are placed on and embedded in elastic half space. Numerical analysis is performed for a containment structure of a nuclear power plant subjected seismic excitation.

• Tunnel and Underground Space
• /
• v.12 no.3
• /
• pp.171-178
• /
• 2002

A dynamic analysis of a horseshoe_shaped tunnel near to cavity was performed to study the effect of the cavity on the dynamic behavior of the tunnel. In order to obtain the dynamic response of the tunnel embedded in a semi-infinite domain, a hybrid numerical technique was primarily developed. A dynamic fundamental solution in frequency domain for multi-layered half planes was derived and subsequently incorporated in the boundary element method. Coupling of the boundary element method for the far field with the finite element method for the near field is made by imposing compatibility condition of a displacement at the interface. The boundary element method is then coupled with the finite element method, which is utilized to model the near field including the tunnel and the cavity. In order to demonstrate the validity of the proposed technique, dynamic responses of single and multiply-layered semi-infinite structural systems are obtained by using the Kicker waveform and investigated in the limestone layer to find how the being and the location of the cavity affect the dynamic characteristics of the system.