• Geosystem Engineering
• /
• v.21 no.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
• /
• v.17 no.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.

• Geotechnical Engineering
• /
• v.14 no.4
• /
• pp.129-140
• /
• 1998

Elasto-plastic finite element analysis of geotechnical boundary value problems necessitate the stress integration for the known strain increments. For the elasto-plastic constitutive model, the stress integration is generally achieved by numerical schemes, because analytical integration is impossible for general strain path. In this case, the accuracy of numerical stress integration has an important role on the overall accuracy of nonlinear finite element solution. In this study, the Sloan's substepping method which is one of explicit integration methods has been adopted and iris applicability has been checked. The unstability and inaccuracy of ifs results initiated from initial stress level were revealed. So. a new modified numerical integration method which employs the basic concept of modified Euler scheme for error control is proposed and accuracy and stability of the solutions are confirmed by triaxial test simulation.

• Proceedings of the KSME Conference
• /
• 2001.06c
• /
• 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.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.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.

• Journal of Power System Engineering
• /
• v.8 no.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.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.3 no.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.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.3 no.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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.7 s.94
• /
• 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.

• Journal of the Korean Society of Industry Convergence
• /
• v.12 no.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.