• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.04a
• /
• pp.128-136
• /
• 1994

In this study, dynamic boundary element method using mass matrix is derived, using fundamental solutions for the semi-infinite domain. In constituting boundary integral equations for the dynamic equilibrium condition, inertia term in the form of domain integral is transformed into boundary integral form. Corresponding system equations are derived, and a boundary element program is developed. In addition, equations for free vibration is formulated, and eigenvalue analysis is performed. The results from the dynamic boundary element analysis for a tunnel problem are compared with those from the finite element analysis. According to the comparison, boundary element method using mass matrix is consistent with the results of finite element method. Consequently, in tunnel vibration problems, it results in reasonable solution compared with other methods where relatively higher degree of freedoms are employed.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.5 no.4
• /
• pp.27-38
• /
• 2001

A mechanical lumped parameter model is proposed for the dynamic modeling of a semi-infinite reservoir. A semi-analytic 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 time domains. Mass, damping and spring coefficients of the mechanical model are obtained to preserve the major features of the solution such as eigenfrequencies and the shapes of Bessel functions that appear as kernels in the convolution integrals. The lumped parameter model in its final form consists of two masses, a spring and two dampers for each eigenfrequency. Application examples demonstrated that the new lumped parameter model could be used for the time domain analysis of dam-reservoir systems.

• Interaction and multiscale mechanics
• /
• v.6 no.4
• /
• pp.357-375
• /
• 2013

In this paper, the forced vibration problem of an Euler-Bernoulli beam that is joined with a semi-infinite field of a compressible fluid is considered as a boundary value problem (BVP). This BVP includes two partial differential equations (PDE) and some boundary conditions (BC), which are introduced comprehensively. After that, the closed-form solution of this fluid-structure interaction problem is obtained in the frequency domain. Some mathematical techniques are utilized, and two unknown functions of the BVP, including the beam displacement at each section and the fluid dynamic pressure at all points, are attained. These functions are expressed as an infinite series and evaluated quantitatively for a real example in the results section. In addition, finite element analysis is carried out for comparison.

• Coupled systems mechanics
• /
• v.3 no.4
• /
• pp.385-403
• /
• 2014

In this paper, the forced vibration problem of an Euler-Bernoulli beam that is joined with a semi-infinite field of a compressible fluid is considered as a boundary value problem (BVP). This BVP includes two partial differential equations (PDE) and some boundary conditions (BC), which are introduced comprehensively. After that, the closed-form solution of this fluid-structure interaction problem is obtained in the frequency domain. Some mathematical techniques are utilized, and two unknown functions of the BVP, including the beam displacement at each section and the fluid dynamic pressure at all points, are attained. These functions are expressed as an infinite series and evaluated quantitatively for a real example in the results section. In addition, finite element analysis is carried out for comparison.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.6
• /
• pp.427-435
• /
• 2021

An analysis considering the fluid-structure interaction is required to strictly evaluate the seismic behavior of facilities such as, environmental facilities and dams, that store fluids. Specifically, in the case of an infinite domain in the upstream direction, such as a dam-reservoir system, this should be carefully considered. In this study, we proposed a practical numerical model for both wave propagation and fluid-structure interaction analyses of an infinite domain, for a system with a semi-infinite domain such as a dam-reservoir system. This method was applicable to the time domain, and enabled accurate boundary analysis. For an infinite fluid domain, a small number of mid-point integrated acoustic finite elements were applied instead of a general acoustic finite element, and a viscous boundary was imposed on the outermost boundary. The validity and accuracy of the proposed method were secured by comparing analytic solutions of a reservoir having infinite domain, with the parametric analysis results, for the number of elements and the size of the modeling region. Furthermore, the proposed method was compared with other fluid-structure interaction methods using additional mass.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.157-175
• /
• 2008

A brief review of the research works on ground vibrations caused by trains moving in underground tunnels is first given. Then, the finite/infinite element approach for simulating the soil-tunnel interaction system with semi-infinite domain is summarized. The tunnel is assumed to be embedded in a homogeneous half-space or stratified soil medium. The train moving underground is modeled as an infinite harmonic line load. Factors considered in the parametric studies include the soil stratum depth, damping ratio and shear modulus of the soil with or without tunnel, and the thickness of the tunnel lining. As far as ground vibration is concerned, the existence of a concrete tunnel may somewhat compensate for the loss due to excavation of the tunnel. For a soil stratum resting on a bedrock, the resonance peak and frequency of the ground vibrations caused by the underground load can be rather accurately predicted by ignoring the existence of the tunnel. Other important findings drawn from the parametric studies are given in the conclusion.

• Korean Chemical Engineering Research
• /
• v.57 no.5
• /
• pp.723-729
• /
• 2019

A theoretical analysis was conducted of convective instability driven by buoyancy forces under transient temperature fields in an annular porous medium bounded by coaxial vertical cylinders. Darcy's law and Boussinesq approximation are used to explain the characteristics of fluid motion and linear stability theory is employed to predict the onset of buoyancy-driven motion. The linear stability equations are derived in a global domain, and then cast into in a self-similar domain. Using a spectral expansion method, the stability equations are reformed as a system of ordinary differential equations and solved analytically and numerically. The critical Darcy-Rayleigh number is founded as a function of the radius ratio. Also, the onset time and corresponding wavelength are obtained for the various cases. The critical time becomes smaller with increasing the Darcy-Rayleigh number and follows the asymptotic relation derived in the infinite horizontal porous layer.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.213-221
• /
• 2013

This paper introduces a mesh continuation scheme for a one-dimensional inverse medium problem to reconstruct the spatial distribution of elastic wave velocities in heterogeneous semi-infinite solid domains. To formulate the inverse problem, perfectly-matched-layers(PMLs) are introduced as wave-absorbing boundaries that surround the finite computational domain truncated from the originally semi-infinite extent. To tackle the inverse problem in the PML-truncated domain, a partial-differential-equations(PDE)-constrained optimization approach is utilized, where a least-squares misfit between calculated and measured surface responses is minimized under the constraint of PML-endowed wave equations. The optimization problem iteratively solves for the unknown wave velocities with their updates calculated by Fletcher-Reeves conjugate gradient algorithms. The optimization is performed using a mesh continuation scheme through which the wave velocity profile is reconstructed in successively denser mesh conditions. Numerical results showed the robust performance of the mesh continuation scheme in reconstructing target wave velocity profile in a layered heterogeneous solid domain.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.1
• /
• pp.45-53
• /
• 1991

The underground structure, which has infinite or semi-infinite boundary conditions, is subjected by body forces and in-situ stresses. It also has stress concentration, which causes material nonlinear behavior, in the vicinity of the excavated surface. In this paper, some methods which can be used to transform domain integrals into boundary integrals are reviewed in order to analyze the effect of the body forces and the in-situ stresses. First, the domain integral of the body force is transformed into boundary integral by using the Galerkin tensor and divergence theorem. Second, it is transformed by writing the domain integral in cylindrical coordinates and using direct integration. The domain integral of the in-situ stress is transformed into boundary integral applying the direct integral method in cylindrical coordinates. The methodology is verified by comparing the results from the boundary element analysis with those of the finite element analysis. Coupling the above boundary elements with finite elements, the nonlinear behavior that occurs locally in the vicinity of the excavation is analyzed and the results are verified. Thus, it is concluded that the domain integrals of body forces and in-situ stresses could be performed effectively by transforming them into the boundary integrals, and the nonlinear behavior can be reasonably analyzed by coupled nonlinear finite element and boundary element method. The result of this research is expected to he used for the analysis of the underground structures in the effective manner.

• Coupled systems mechanics
• /
• v.2 no.4
• /
• pp.375-387
• /
• 2013

Soils consist of an assemblage of particles with different sizes and shapes which form a skeleton whose voids are filled with water and air. Hence, soil behaviour must be analyzed by incorporating the effects of the transient flow of the pore-fluid through the voids, and therefore requires a two-phase continuum formulation for saturated porous media. The present paper presents briefly the Biot's basic theory of dynamics of saturated porous media with u-P formulation to determine the responses of pore fluid and soil skeleton during cyclic loading. Kelvin elements are attached to transmitting boundary. The Pastor-Zienkiewicz-Chan model has been used to describe the inelastic behavior of soils under isotropic cyclic loadings. Newmark-Beta method is employed to discretize the time domain. The response of fluid-saturated porous media which are subjected to time dependent loads has been simulated numerically to predict the liquefaction potential of a semi-infinite saturated sandy layer using finite-infinite elements. A settlement of 17.1 cm is observed at top surface. It is also noticed that liquefaction occurs at shallow depth. The mathematical advantage of the coupled finite element analysis is that the excess pore pressure and displacement can be evaluated simultaneously without using any empirical relationship.