• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2000.10a
• /
• pp.143-150
• /
• 2000

A physical lumped parameter model is proposed for the time domain analysis of dam-reservoir system. The exact solution of transmitting boundary is derived for a semi-infinite 2-D reservoir of constant depth. The characteristics of the solution are examined in both frequency and the domains. Mass and damping coefficient are obtained from asymptotic behavior of the frequency domain solution. Further refinement to the lumped model is made by approximating the kernel function of the convolution integral in the exact solution. Finally a new physical lumped parameter model is proposed that consists of two masses, a spring and two dampers for each mode. It is demonstrated that new lumped parameter model of transmitting boundary can give excellent results.

• 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.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.37-43
• /
• 2006

Dynamic numerical analysis of geotechnical problems requires a way to simulate the decrease of energy as the domain of interest gets larger. This phenomenon is usually referred to as radiation damping or geometric attenuation and it is distinguished from material damping in which elastic energy is actually dissipated by viscous, hysteretic, or other mechanism. The fact that the domain of analysis in numerical modeling must be chosen, however, causes a need for special attention at the boundary. This observation leads directly to the idea of determining the dynamic response of the interior region from a finite model consisting of the interior region subjected to a boundary condition which ensures that all energy arriving at the boundary is absorbed. This paper presents a simple methodology to simulate transmitting boundaries condition using viscoelastic infinite elements within the recently developed "OpenSees" finite element code. The methodology used here provides that the level of absorption for traveling waves is efficient enough for practical purposes, but unsatisfactory for the case of sharp incident angles. The effectiveness of the infinite elements for the absorption of incident waves at boundaries is evaluated via example analysis.

• Structural Engineering and Mechanics
• /
• v.50 no.3
• /
• pp.349-364
• /
• 2014

Perfectly matched layers are employed in time harmonic analysis of dam-foundation systems. The Lysmer boundary condition at the truncation boundary of the PML region has been incorporated in the formulation of the dam-foundation FE model (including PML). The PML medium is defined in a way that the formulation of the system can be transformed into time domain. Numerical experiments show that applying Lysmer boundary conditions at the truncation boundary of the PML area reduces the computational cost and make the PML approach a more efficient technique for the analysis of dam-foundation systems.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2002.10a
• /
• pp.820-823
• /
• 2002

The asymptotic problem of a semi-infinite conducting crack parallel to the poling direction in ferroelectric ceramics subjected to electric fields is analyzed. The main mechanism for the conducting crack growth behavior is thought to be ferroelectric domain switching leading to the development of a process zone around the crack. The shape and size of the switching zone is shown to depend strongly on the relative magnitude on the ratio of the coercive electric field to the yield electric field. It is shown that the crack growth can be either enhanced or retarded depending on the ratio of the coercive electric field to yield electric field.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11a
• /
• pp.240-243
• /
• 2005

A design procedure using spatial Fourier transform is presented for a structural-acoustic coupled radiator that can emit sound in the desired direction with high power and low side lobe level. The design procedure consists of three steps. Firstly, the structural-acoustic coupled radiator is chosen to obtain strong coupling between structural vibration and acoustic pressure. The radiator is composed by two spaces which are separated by a wall. Spaces can be categorized as reverberant finite space and unbounded semi-infinite space, and the wall are composed of two plates and an opening. The velocities on the wall are predicted. Secondly, directivity and energy distribution of radiator are predicted in wave number domain using spatial Fourier transform. Finally, optimal design variables are calculated using a dual optimal algorithm. Its computational example is presented including the directivity and resulting pressure distribution using proposed procedure.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• Earthquakes and Structures
• /
• v.3 no.3_4
• /
• pp.533-548
• /
• 2012

Dynamic analysis of concrete gravity dam-reservoir systems is an important topic in the study of fluid-structure interaction problems. It is well-known that the rigorous approach for solving this problem relies heavily on employing a two-dimensional semi-infinite fluid element. The hyper-element is formulated in frequency domain and its application in this field has led to many especial purpose programs which were demanding from programming point of view. In this study, a technique is proposed for dynamic analysis of dam-reservoir systems in the context of pure finite element programming which is referred to as the wavenumber approach. In this technique, the wavenumber condition is imposed on the truncation boundary or the upstream face of the near-field water domain. The method is initially described. Subsequently, the response of an idealized triangular dam-reservoir system is obtained by this approach, and the results are compared against the exact response. Based on this investigation, it is concluded that this approach can be envisaged as a great substitute for the rigorous type of analysis.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.1 no.4
• /
• pp.59-68
• /
• 1997

In analysis of underground structures, the effects of artificial boundary conditions are considered as one of the major reasons for differences from experimental results. These phenomena can be overcome by using the boundary elements which satisfy the multi-layered half space conditions. The fundamental solutions of multi-layered half-space for boundary element method is formulated satisfying the transmission and reflection of waves at each layer interface and radiation conditions at bottom layer. The governing equations can be obtained from the displacements at each layer which are expressed in terms of harmonic functions. All types of waves can be included using the complete response from semi-infinite integrals with respect to horizontal wavenumbers using expansion of Fourier series and Hankel transformation. Two dimensional Green's functions are derived from cylindrical Navier equations and potentials performing infinite integration in y-direction. In this case, it is effective to transform into two dimensional problem using semi-analytical integration and sinusoidal Bessel function. Some verifications are given to show the accuracy and efficiency of the developed method, and numerical examples to demonstrate the dynamic behavior of underground with various properties.

• Structural Engineering and Mechanics
• /
• v.20 no.2
• /
• pp.225-239
• /
• 2005

Applications of the Trefftz boundary element method (BEM) to anti-plane electroelastic problems are presented in this paper. Both direct and indirect methods with domain decomposition are discussed in details. Each crack is treated as semi-infinite thin slit defined in a subregion, for which a particular solution of the anti-plane problem, satisfying exactly the crack-face condition, is derived. The stress intensity factors defined at each crack tip can be directly computed from the coefficients of the particular solution. The performance of the proposed formulation is assessed by two examples and comparison is made with results obtained by other approaches. The Trefftz boundary element approach is demonstrated to be suitable for the analysis of the anti-plane problem of piezoelectric materials.