• Geosystem Engineering
• /
• 제21권6호
• /
• pp.309-317
• /
• 2018

We employ a three-dimensional indirect boundary element method (BEM) to simulate temperature change around an underground liquefied natural gas storage cavern. The indirect BEM (IBEM) uses fictitious heat source strength on boundary elements as basic variables which are solved from equations of boundary conditions and then used to compute the temperature change at other points in the considered problem domain. The IBEM requires evaluation of singular integration for temperature change due to heat conduction from a constant heat source on a planar (triangular) region. The singularity can be eliminated by a semi-analytical integration scheme. However, it is found that the semi-analytical integration scheme yields sharp temperature gradient for points close to vertices of triangle. This affects the accuracy of heat flux, if they are evaluated by finite difference method at these points. This difficulty can be overcome by a combination of using a direct numerical integration for these points and the semi-analytical scheme for other points distance away from the vertices. The IBEM and the hybrid integration scheme have been verified with an analytic solution and then used to the application of the underground storage.

• Structural Engineering and Mechanics
• /
• 제17권6호
• /
• pp.735-749
• /
• 2004

This work presents a direct integration scheme, based on a fourth order finite difference approach, for elastodynamics. The proposed scheme was chosen as an alternative for attenuating the errors due to the use of the central difference method, mainly when the time-step length approaches the critical time-step. In addition to eliminating the spurious numerical oscillations, the fourth order finite difference scheme keeps the advantages of the central difference method: reduced computer storage and no requirement of factorisation of the effective stiffness matrix in the step-by-step solution. A study concerning the stability of the fourth order finite difference scheme is presented. The Finite Element Method and the Boundary Element Method are employed to solve elastodynamic problems. In order to verify the accuracy of the proposed scheme, two examples are presented and discussed at the end of this work.

• 한국지반공학회지:지반
• /
• 제14권4호
• /
• pp.129-140
• /
• 1998

탄소성 유한 요소 해석은 기지의 변형률 증분에 대한 응력 적분을 필요로 하며, 탄소성 구성 모델의 경우 특별한 경우를 제외하고는 해석적인 응력 적분이 불가능하고 수치적인 방법을 필요로 한다. 이때 응력 수치 적분의 정확도가 비선형 유한요소해의 전체적인 정확도에 상당히 큰 영향을 미치게 된다. 본 연구에서는 탄소성 구성 관계의 응력 적분을 위하여 외연적 방법중의 하나로서 Sloan이 제안한 단계분할 절차를 보완하여 안정적이고 정착한 응력 수치 적분법을 제시하고자 한다. 수정 오일러 절차에 따른 오차 조절의 기본 개념은 그대로 사용하고 오차를 평가하는 기준에 응력 수준이 영향을 미치는 단점을 보완하여 응력 수준에 관계없는 안정적이고 정확한 수치 적분법을 제시하였으며, 그 결과의 신뢰성을 삼축시험모사를 이용하여 검증하였다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집C
• /
• pp.402-407
• /
• 2001

Sensitivity analysis scheme is developed in the elasto-plastic finite element method with explicit time integration using direct differentiation method. The direct differentiation is concerned with the time integration, constitutive relation, shell element with reduced integration and the contact scheme. Sensitivity analysis results are mainly examined with the highly nonlinear and quasi-static problem with the complicated contact condition. The result shows stable sensitivity especially in the sheet metal forming analysis.

• 한국정밀공학회지
• /
• 제12권3호
• /
• pp.32-41
• /
• 1995

In the analysis of metal forming processes by the finite element method, there are many numerical instabilities such as element locking, hourglass mode and shear locking. These instabilities may have a bad effect upon accuracy and convergence. The present work is concerned with improvement of stability and efficiency in two-dimensional rigid-plastic finite element method using various type of elemenmts and numerical intergration schemes. As metal forming examples, upsetting and backward extrusion are taken for comparison among the methods: various element types and numerical integration schemes. Comparison is made in terms of stability and efficiency in element behavior and computational efficiency and a new scheme of adaptive directional reduced integration is introduced. As a result, the finite element computation has been stabilized from the viewpoint of computational time, convergency, and numerical instability.

• 동력기계공학회지
• /
• 제8권3호
• /
• pp.23-29
• /
• 2004

In this paper, a finite element formulation using dynamic-explicit time integration scheme is used for numerical analysis of Hydro-forming processes. The lumping scheme is employed for the diagonal mass matrix and dynamic explicit formulation. Hydro-forming process for auto-body panel forming is analyzed by using dynamic-explicit finite element method. Further, the simulated results of the Hydro-forming processes are shown and discussed. Its application is being increased especially in the automotive industrial area for the cost reduction, weight saving, and improvement of strength.

• 한국기계가공학회지
• /
• 제3권3호
• /
• pp.58-63
• /
• 2004

In the present work a finite element formulation using dynamic-explicit time integration scheme is used for numerical analysis of multi-stage stamping processes. The lumping scheme is employed for the diagonal mass matrix and dynamic explicit formulation Multi-Stage stamping is analyzed by using dynamic-explicit finite element method. Further, the simulated results for the panel stamping processes are shown and discussed. Its application is being increased especially in the stamping industrial area for the cost reduction, weight saving, and improvement of strength.

• 한국기계가공학회지
• /
• 제3권4호
• /
• pp.63-72
• /
• 2004

In this paper, a finite element formulation using dynamic-explicit time integration scheme is used for numerical analysis of auto-body panel stamping processes. The lumping scheme is employed for the diagonal mass matrix and dynamic explicit formulation. Auto-body panel forming is analyzed by using dynamic-explicit finite element method. Further, the simulated results of the auto-body panel stamping processes are shown and discussed. Its application is being increased especially in the stamping industrial area for the cost reduction, weight saving, and improvement of strength.

• 대한기계학회논문집
• /
• 제17권7호
• /
• pp.1719-1730
• /
• 1993

This paper presents the theoretical natural frequencies of out-of-plane deformable ring based on the variables such as out-of-plane deflection, torsional rotation and shear rotation. Based on the same variables, a finite element eigen analysis is carried out by using the $C^0$-continuous, isoparametric element which has three nodes per element and three degrees-of-freedom at each node. Numerical experiments are peformed to find the integration scheme which produces accurate natural frequencies, natural modes and correct rigid body motion. The uniformly reduced integration and the selective reduced integration give more accurate numerical frequencies than the uniformly full integration, but the uniformly reduced integration produces incorrect rigid body motion while selective reduced integration does correct one. Therefore, the ring element based on the three variables which employes selective reduced integration is recommended to avoid spurious modes, to alleviate the error due to shear locking and to produce correct rigid body motion, simultaneously.

• 한국산업융합학회 논문집
• /
• 제12권3호
• /
• pp.149-155
• /
• 2009

A new finite element equation is derived by applying quadratic and cubic time integration scheme to the variational formulation in time-integral for the analysis of the transient elastodynamic problems to increase the numerical accuracy and stability. Emphasis is focused on methodology for cubic time integration scheme procedure which are never presented before. In this semidiscrete approximations of the field variables, the time axis is divided equally and quadratic and cubic time variation is assumed in those intervals, and space is approximated by the usual finite element discretization technique. It is found that unconditionally stable numerical results are obtained in case of the cubic time variation. Some numerical examples are given to show the versatility of the presented formulation.