• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.121-128
• /
• 1999

This paper presents a time domain method for soil-structure interaction analysis for seismic loadings. It is based on the finite element formulation incorporating analytical frequency-dependent infinite elements for the far field soil. The dynamic stiffness matrices of the far field region formulated using the present method in frequency domain can be easily transformed into the corresponding matrices in time domain. At first, the equivalent earthquake forces are evaluated along the interface between the near and the far fields from the free-field response analysis carried out in frequency domain, and the results are transformed into the time domain. An efficient procedure is developed for the convolution integrals to evaluate the interaction force along the interface, which depends on the response on the interface at the past time instances as well as the concurrent instance. Then, the dynamic responses are obtained for the equivalent earthquake force and the interaction force using Newmark direct integration technique. Since the response analysis is carried out in time domain, it can be easily extended to the nonlinear analysis. Example analysis has been carried out to verify the present method in a multi-layered half-space.

• Structural Engineering and Mechanics
• /
• v.13 no.6
• /
• pp.615-638
• /
• 2002

This paper presents a method of seismic analysis for a cylindrical liquid storage structure considering the effects of the interior fluid and exterior soil medium in the frequency domain. The horizontal and rocking motions of the structure are included in this study. The fluid motion is expressed in terms of analytical velocity potential functions, which can be obtained by solving the boundary value problem including the deformed configuration of the structure as well as the sloshing behavior of the fluid. The effect of the fluid is included in the equation of motion as the impulsive added mass and the frequency-dependent convective added mass along the nodes on the wetted boundary of the structure. The structure and the near-field soil medium are represented using the axisymmetric finite elements, while the far-field soil is modeled using dynamic infinite elements. The present method can be applied to the structure embedded in ground as well as on ground, since it models both the soil medium and the structure directly. For the purpose of verification, earthquake response analyses are performed on several cases of liquid tanks on a rigid ground and on a homogeneous elastic half-space. Comparison of the present results with those by other methods shows good agreement. Finally, an application example of a reinforced concrete tank on a horizontally layered soil with a rigid bedrock is presented to demonstrate the importance of the soil-structure interaction effects in the seismic analysis for large liquid storage tanks.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.7 no.6
• /
• pp.1-7
• /
• 2003

This paper presents a novel hybrid time-frequency-domain method for nonlinear soil-structure interaction(SSI) analysis. It employs, in a practical manner, a computer code for equivalent linear SSI analysis and a general-purpose nonlinear finite element program. The proposed method first (calculates dynamic responses on a truncated finite element boundary utilizing an equivalent linear SSI program in the frequency domain. Then, a general purpose nonlinear finite element program is employed to analyze the nonlinear SSI problem in the time domain, in which boundary conditions at the truncated boundary are imposed with the responses calculated in the previous frequency domain SSI analysis, In order to validate the proposed method, seismic response analyses are carried out for a 2-D underground subway station in a multi-layered half-space, For the analyses, a equivalent linear SSI code KIESSI-2D is coupled to ANSYS program. The numerical results indicate that the proposed methodology can be a viable solution for nonlinear SSI problems.

• Geophysics and Geophysical Exploration
• /
• v.9 no.2
• /
• pp.148-154
• /
• 2006

The effect of small-scale near surface inhomogeneity on Rayleigh wave propagation and dispersion has been investigated in this study using two-dimensional FEM elastic modeling. Various inhomogeneity models with a variety of geometrical shape and embedment depth which exist in homogeneous half-space and two-layered media are considered. Results show that any near surface inhomogeneity greater than one wavelength in terms of minimum wavelength of Rayleigh wave shows dispersion characteristics. Such dispersion effect become stronger as the dimensions of the inhomogeneity increase. The effect of horizontal dimension is more dominant factor governing the dispersion characteristics than vertical dimension. However, the dispersion effect can not be identifiable in seismogram if the horizontal dimension is not wide enough. Nonetheless, even in this case, the existence of inhomogeneity can be inferred by the reflection or transmission event of Rayleigh wave. The results can be expected to provide insights on the behavior of Rayleigh wave which may be helpful for designating field work or new processing scheme to detect near surface inhomogeneity by surface wave method.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.361-366
• /
• 2009

Layered-earth Green's functions in electormagnetic (EM) surveys play a key role in modeling the response of exploration targets. They are computed through the Hankel transforms of analytic kernels. Computational precision depends upon the choice of algebraically equivalent forms by which these kemels are expressed. Since three-dimensional (3D) modeling can require a huge number of Green's function evaluations, total computational time can be influenced by computational time for the Hankel transform evaluations. Linear digital filters have proven to be a fast and accurate method of computing these Hankel transforms. In EM modeling for 3D inversion, electric fields are generally evaluated by the secondary field formulation to avoid the singularity problem. In this study, three components of electric fields for five different sources on the surface of homogeneous half-space were derived as primary field solutions. Moreover, reflection coefficients in TE and TM modes were produced to calculate EM responses accurately for a two-layered model having a sea layer. Accurate primary fields should substantially improve accuracy and decrease computation times for Green's function-based problems like MT problems and marine EM surveys.

• Geophysics and Geophysical Exploration
• /
• v.17 no.3
• /
• pp.129-136
• /
• 2014

Electromagnetic (EM) methods are generally divided into frequency-domain EM (FDEM) and time-domain EM (TDEM) methods, depending on the source waveform. The FDEM and TDEM fields are mathematically related by the Fourier transformation, and the TDEM field can thus be obtained as the Fourier transformation of FDEM data. For modeling in time-domain, we can use fast frequency-domain modeling codes and then convert the results to the time domain with a suitable numerical method. Thus, frequency-to-time transformations are of interest to EM methods, which is generally attained through fast Fourier transform. However, faster frequency-to-time transformation is required for the 3D inversion of TDEM data or for the processing of vast air-borne TDEM data. The diffusion expansion method (DEM) is one of smart frequency-to-time transformation methods. In DEM, the EM field is expanded into a sequence of diffusion functions with a known frequency dependence, but with unknown diffusion-times that must be chosen based on the data to be transformed. Especially, accuracy of DEM is sensitive to the diffusion-time. In this study, we developed a method to determine the optimum range of diffusion-time values, minimizing the RMS error of the frequency-domain data approximated by the diffusion expansion. We confirmed that this method produces accurate results over a wider time range for a homogeneous half-space and two-layered model.

• Geophysics and Geophysical Exploration
• /
• v.5 no.2
• /
• pp.108-117
• /
• 2002

We have introduced a new approach to obtain the conductivity information of subsurface using Cagniard impedance over two-dimensional (2-D) model in the presence of horizontal magnetic dipole source with the frequency range of $1\;kHz\~1\;MHz$. Firstly, we designed the method to calculate the apparent resistivity from the ratio between horizontal electric and magnetic fields, Cagniard impedance, considering the source effects when the plane wave assumption is failed in finite source EM problem, and applied it to several numerical models such as homogeneous half-space or layered-earth model. It successfully provided subsurface information even though it is still rough, while the one with plane wave assumption is hard to give useful information. Next, through analyzing Cagniard impedance and apparent resistivity considering source effect over 2-D models containing conductive- or resistive-block, we showed that the possibility of obtaining conductivities of background media and anomaly using this approach. In addition, the apparent resistivity considering source effect and phase pseudosections constructed from Cagniard impedance over the isolated conductive- and resistive block model well demonstrated outlines of anomalies and conductivity distribution even though there were some distortions came from sidelobes caused by 2-D body.