• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.2
• /
• pp.649-658
• /
• 1996

Small-amplitude harmonic oscillations of multi-cylinders are considered both experimentally and theoretically. For the theoretical model, the flow regime is separated into inner and outer regions. In the inner region, the flow is governed by the generalized Stokes boundary layer equation. In the outer region, the full Navier-Stokes equation for the steady streaming flow is solved numerically by using ADI scheme and FVM coupled with the boundary integral method. Flow visualization experiments are conducted by using the Laser Sheet Image Technique. The case of two circular cylinders and square cylinders with variable distances are chosen as a typical example. Although experimental results are based on the flow in the finite domain, both experimental and numerical results agree well qualitatively. As the separation of cylinders is increased, a numerical result shows the asymptotic convergence to a single cylinder case.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.23 no.7
• /
• pp.1764-1769
• /
• 1998

The FE-BI (finite element- boundary integral) method is used for the characterization of an arbitrary-shaped multiaperture-coupled directional coupler between two retangular waveguides crossing with an arbitrary angle. The method is shown to be much more effective than the ordinary FEM for the structures under consideration. the computation time is compared with that of the ordinary FEM. The results of other papers for Moreno coupler are compared with the theoretical prediction and a good agreement is obtained.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.1
• /
• pp.67-81
• /
• 2000

An accurate and efficient numerical method for two-dimensional nonlinear radiation problem has been developed. The wave motion due to a moving body is described by the assumption of ideal fluid flow, and the governing Laplace equation can be effectively solved by the higher-order boundary element method with the help of the GMRES (Generalized Minimal RESidual) algorithm. The intersection or corner problem is resolved by utilizing the so-called discontinuous elements. The implicit trapezoidal rule is used in updating solutions at new time steps by considering stability and accuracy. Traveling waves caused by the oscillating body are absorbed downstream by the damping zone technique. It is demonstrated that the present method for time marching and radiation condition works efficiently for nonlinear radiation problem. To avoid the numerical instability enhanced by the local gathering of grid points, the regriding technique is employed so that all the grids on the free surface may be distributed with an equal distance. This makes it possible to reduce time interval and improve numerical stability. Special attention is paid to the local flow around the body during time integration. The nonlinear radiation force is calculated by the "acceleration potential technique". Present results show good agreement with other numerical computations and experiments.

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

• Journal of Korean Society of Steel Construction
• /
• v.9 no.3 s.32
• /
• pp.333-340
• /
• 1997

This paper will expand the third-order shear deformation theory by the double-Fourier series and reduce to the solution of a system of ordinary differential equations in time, which are integrated numerically using Newmark's direct integration method and clarify the undamped dynamic responses for the cross-ply and antisymmetric angle-ply laminated composite plates and shells with simply supported boundary condition. Numerical results for deflections are presented showing the effect of side-to-thickness ratio, aspect ratio, material anisotropy, and lamination scheme.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.5
• /
• pp.1551-1562
• /
• 1991

최근까지 발표된 유한 요소법을 이용한 압연 해석 관련 주요 논문들을 정리해 보면 다음과 같다. Li와 Kobayashil는 강소성 유한 요소법(rigidplastic finite element method)을 여러가지 마찰조건에 대하여 해석하였다. 이때 압연롤은 강체 (rigid body)로 시편은 가공경화(workhardening)를 동반한 강소성체로 모델링하였다. Hwang과 Kobayashi는 강소성 유한 요소법을 이용한 평면 변형 압연에서 재료 손실을 최소화하는 예비 성형체(preform)의 설계에 대한 연구를 수행하였다. 이 경우에도 역시 압연롤은 강체로 시편은 가공 경화를 동반한 강소성체와 완전 소성체로 모델링 되었으나, 고착(sticking) 마찰 조건에 대해서만 해석을 수행하였다. Mori와 Osak- ada 그리고 Oda는 약간 압축성이 있는 재료의 평면 변형 압연에 대하여 연구하였다. 이때 압연롤은 강체로 시편은 가공 경화를 동반한 강소성체로 모델링 되었으며 경계 면에서는 Coulomb 마찰을 고려하였다. 이밖에도 오일러(Eulerian) 수식화를 이용한 Dawson과 Thompson, Berman의 해석 결과가 있으며, 또 폭 방향의 변형까지를 고려한 Li와 Kobayashi, Mori와 Osakada의 3차원 해석 결과가 있다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.1
• /
• pp.117-124
• /
• 2009

The p-convergent boundary element method has been proposed to analyze two-dimensional potential problem on the basis of high order Legendre shape functions that have different property comparing with the shape functions in conventional boundary element method. The location of nodes corresponding to high order shape function are not defined along the boundary, called by nodeless node, similar to the p-convergent finite element method. As the order of shape function increases, the collocation point method is used to solve linear simultaneous equations. The collocation patterns of p-convergent boundary element method consist of non-symmetric hierarchial or symmetric non-hierarchical. As the order of shape function increases, the number of collocation point increases. The singular integral that appears in p-convergent boundary element has been calculated by special numeric quadrature technique and semi-analytical integration technique. The L-shape domain problem including singularity in the vicinity of reentrant comer is analyzed and the numerical results show that the relative error is smaller than $10^{-2}%$ range as compared with other results in literatures. In case of same condition, the symmetric p-collocation point pattern shows high accuracy of solution.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1634-1638
• /
• 2008

하천의 2차원 흐름 해석, 유사이동 해석, 오염확산 해석을 위한 유체의 수치해석법에는 유한요소법, 유한차분법, 유한차분법의 변형인 유한체적법, 경계적분법 등이 있다. 유체에 대한 수치해석 기법으로 전통적으로 가장 많이 사용되고 있는 방법은 유한차분법이지만, 비구조적 요소망(unstructured mesh)을 이용하여 복잡한 형상을 표현하기가 상대적으로 용이한 유한요소법이 다양한 형태의 하천 해석에는 더욱 적합할 것이다. 본 연구에서는 비구조적 요소망을 advanced front method를 이용하여 생성해 보았다. Advanced front method는 해석하고자 하는 영역에 적절한 절점들을 생성한 후 삼각 요소망을 구성하는 grid based advanced front method와 절점들을 생성하지 않고 해석 영역에 삼각 요소를 바로 구성하는 element based advanced front method로 나눌 수 있다. Grid based advanced front method에서 해석 영역에 적절한 절점을 생성하는 방법으로는 일반적인 격자 구조의 절점 생성 방법을 적용하였으며 경계와의 거리가 가까운 절점은 생성되지 않으며, 삼각 요소를 구성할 때에는 두 개의 인접 절점을 비교하여 최적의 삼각 요소를 구성하게 된다. 단 두 개의 인접 절점만을 비교함으로서 비교적 빠른 시간 안에 최적의 삼각 요소망을 구성할 수 있다. 삼각 요소망을 생성한 후에는 Laplacian smoothing을 이용하여 삼각 요소망의 형질을 개선하였다. Element based advanced front method는 외부 경계에서부터 시작된 Front가 내부 영역으로 확대되어지며 각 Front에서 적절한 절점을 직접 생성하여 바로 삼각 요소를 구성하게 된다. Front에서 생성된 절점은 인접 절점들이 있는지 검색하여 인접 절점이 있다면 생성된 절점은 삭제되어지며 인접 절점이 삼각 요소를 위한 나머지 한 점으로 채택되어진다. Front는 외부 경계와 교차되어지지 않아야 하며 또한 연속된 Front를 효율적으로 관리하기 위해 list 자료 구조를 활용하였다.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.34D no.4
• /
• pp.1-7
• /
• 1997

This paper presents an analysis of the wave propagation in cross-type tunnels using the finite volume time domain (FVTD) method. Because the FVTD method is based on the volume intergrations of themaxwell's equations with respected to arbitrary shaped small polyhedron cells and the fields at every center point of the cells ar eassigned in a average fashion, the method can handle arbitrary boundary problems with inhomogeneous media. In this paper, the wave propagation in cross-type tunnels has been analyzed using the fVTD method with the PML (perfectly matched layer) absorbing bundary conditons, and the numerical results are verified with a set of experimental data.

• Geophysics and Geophysical Exploration
• /
• v.15 no.2
• /
• pp.57-65
• /
• 2012

Various numerical methods in simulation of seismic wave propagation have been developed. Recently an innovative numerical method called as the Spectral Element Method (SEM) has been developed and used in wave propagation in 3-D elastic media. The SEM that easily implements the free surface of topography combines the flexibility of a finite element method with the accuracy of a spectral method. It is generally used a weak formulation of the equation of motion which are solved on a mesh of hexahedral elements based on the Gauss-Lobatto-Legendre integration rule. Variational formulations of velocity-stress motion are newly modified in order to implement ADE-PML (Auxiliary Differential Equation of Perfectly Matched Layer) in wave propagation in 3-D elastic media, because a general weak formulation has a difficulty in adapting CFS (Complex Frequency Shifted) PML (Perfectly Matched Layer). SEM of Velocity-Stress motion having ADE-PML that is very efficient in absorbing waves reflected from finite boundary is verified with simulation of 1-D and 3-D wave propagation.