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

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.15 no.1
• /
• pp.1-8
• /
• 1990

In this paper, the transient electromagneti wave scattering at dielectric cylinder is studied by new numerical analysis method. Basic formulation of boundary integral equation (BIE) for numerical method is started weighted residual technique. BIE is made to two simultaneous equation at surface inner and outside point of dielectric cylinder in extended boundary condition (EBC) and surface boundary condition (SBC). Numerical method is used Boundary element method (BEM) that is two form, one is direct method and the other is indirect method, so that this method that transformes operator inversion martics is used numerical analysis. A good agreement of this numerical solution and the other results is obtained.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.5 no.5
• /
• pp.1-11
• /
• 2001

The purpose of this study is to develop the main program and pre/post processor for the geometric and plastic non-linear analysis of steel jacket structures subjected to wave and earthquake loadings. In this paper, steel jackets are modelled using geometric non-linear space frames and wave loadings re evaluated based on Morrison equation using the linear Airy theory and the fifth Stokes theory. Random wave is generated using JONSWAP spectrum. For earthquake analysis, dynamic analysis is performed using artificial earthquake time history. Also the plastic hinge method is presented for limit analysis of steel jacket. In the companion paper, the pre/post processor is developed and the numerical examples are presented for linear and non-linear dynamic analysis of steel jackets.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.321-327
• /
• 2007

This study presents a new gridless finite difference method for solving elastic crack problems. The method constructs the Taylor expansion based on the MLS(Moving Least Squares) method and effectively calculates the approximation and its derivatives without differentiation process. Since no connectivity between nodes is required, the modeling of discontinuity embedded in the domain is very convenient and discontinuity effect due to crack is naturally implemented in the construction of difference equations. Direct discretization of the governing partial differential equations makes solution process faster than other numerical schemes using numerical integration. Numerical results for mode I and II crack problems demonstrates that the proposed method accurately and efficiently evaluates the stress intensity factors.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.9
• /
• pp.674-684
• /
• 2002

In order to decrease remarkably the computation time and storage used in the direct integration method without the loss of accuracy, authors suggest a new transient analysis algorithm. This algorithm is derived from the combination of three techniques, that is, the transfer technique of the transfer stiffness coefficient method, the modeling technique of the finite element method, and the numerical integration technique of the Newmark method. In this paper, the transient analysis algorithm of a frame structure is formulated by the proposed method. The accuracy and computation efficiency of the proposed method are demonstrated through the comparing with the computation results by the direct integration method for three computation models under various excitations.

• Computational Structural Engineering
• /
• v.9 no.3
• /
• pp.187-197
• /
• 1996

An explicit direct time integration method based solution algorithm is presented to predict dynamic buckling response of Euler column. On the basis of large deflection beam theory, a plane frame finite element is formulated and implemented into the solution algorithm. The element formulation takes into account geometrical nonlinearity and overall buckling of steel structural frames. The solution algorithm employs the central difference method. Using the computer program developed by the author, dynamic instability behavior of Euler column under impact loading is investigated by considering the time variation of load, load magnitude, and load duration. The free vibration of Euler column caused by a short duration impact load is also studied. The validity and efficiency of the present formulation and solution algorithm are verified through illustrative numerical examples.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.3007-3019
• /
• 1993

For the same configuration, the stiffness of 6-node two dimensional isoparametric is stiffer than that of 8-node two dimensional isoparametric element. This phenomenon may be called 'Relative Stiffness Stiffening Phenomenon.' In this paper, the relative stiffness stiffening phenomenon was studied, and could be corrected by modifying the position of Gauss integral points used in the numerical integration of the stiffness matrix. For the same deformation (bending) energy of 6-node and 8-node two dimensional isoparametric elements, Gauss integral points of 6-node element have to move closer, in comparison with those of 8-node element, in the case of numerical integration along the thickness direction.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.82-98
• /
• 1994

A mathematical model for the hydrodynamic forces acting on the ship manoeuvring in waves is formulated and a numerical method for the problem is developed. The motion of a ship, which manoeuvres in waves, may be thought to have two components; one is a high frequency component due to encounter waves, and the other is a low frequency component due to manoeuvring motion. So the method of two time scale expansion is used to divide linear boundary value problem. For the effects of waves on the manoeuvring motion of a ship, only the second order drift forces are considered. The integral equation for the velocity potential is solved by 3 dimensional panel method and hydrodynamic forces are calculated by direct integral method.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.6 no.6
• /
• pp.542-546
• /
• 2005

The J-integral used as a ductile crack initiation criterion has been discussed for the impact loaded elastic-plastic 3PB specimens. The experimental method to measure or estimate the J-integral history has been investigated and its result has been compared to the obtained elastic-plastic values by the finite element model of this study. These numerical results and the experimental curves are found to agree closely. J-integral can be calculated by only numerical analysis with the finite element model. It is proved that simple calculation can be made in order to find the possible value of J-integral by crack mouth opening displacement(CMOD) in the dynamic nonlinear fracture experiment of 3-point bend(3PB) specimen. The property of elastic-plastic material is considered at different impact velocities. The J-integral may be estimated from the crack mouth opening displacement which can be measured directly kom photographs taken during impact experiments.

• Proceedings of the Acoustical Society of Korea Conference
• /
• autumn
• /
• pp.227-230
• /
• 2000

본 논문에서는 외부방사 문제에 널리 사용되는 음향경계 요소 법(Acoustic Boundary Element Method)을 확장하여 저주파수 대역은 기존의 방법을 이용하고 고주파수 대역에서는 경계요소의 특성길이에 의한 고주파수 상한을 극복할 수 있는 방법을 제시한다. 밴드를 대표하는 음압 레벨을 직접 계산하도록 상호상관도 가정을 통하여 간략화 하였다. 간단한 음향방사 모델인 배플이 있는 빔과 평판에 대하여 Rayleigh 적분식의 계산결과와 제안된 방법을 비교하여 적용효능성을 검토하였다. 적용을 위한 조건으로는 방사면 진동의 상호상관도가 작아야 하므로, 이를 만족시키기 위한 일반적인 조건식을 제시하여 이로부터 경계요소의 특성길이를 정하는 기준을 제시하였다.