• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1992.03a
• /
• pp.159-176
• /
• 1992

Limit analysis has been rendered versatile in many problems such as structural problems and metal forming problems. In metal forming analysis, a slip-line method and an upper bound method approach to limit solutions is considered as the most challenging areas. In the present work, a general algorithm for limit solutions of plastic flow is developed with the use of finite element limit analysis. The algorithm deals with a generalized Holder inequality, a duality theorem, and a combined smoothing and successive approximation in addition to a general procedure for finite element analysis. The algorithm is robust such that from any initial trial solution, the first iteration falls into a convex set which contains the exact solution(s) of the problem. The idea of the algorithm for limit solution is extended from rigid/perfectly-plastic materials to work-hardening materials by the nature of the limit formulation, which is also robust with numerically stable convergence and highly efficient computing time.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2667-2677
• /
• 1995

A new numerical algorithm using equal order linear finite element and fractional step method has been developed which is capable of analyzing unsteady fluid flow and heat transfer problems. Streamline Upwind Petrov-Galerkin (SUPG) method is used for the weighted residual formulation of the Navier-Stokes equations. It is shown that fractional step method, in which pressure term is splitted from the momentum equation, reduces computer memory and computing time. In addition, since pressure equation is derived without any approximation procedure unlike in the previously developed SIMPLE algorithm based FEM codes, the present numerical algorithm gives more accurate results than them. The present algorithm has been applied preferentially to the well known bench mark problems associated with steady flow and heat transfer, and proves to be more efficient and accurate.

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

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2007.05a
• /
• pp.340-343
• /
• 2007

This paper introduces a new approach to consider the bending history in finite element inverse analysis of the cylindrical cup drawing. A modified membrane element is adopted to add the bending-unbending energy to the total plastic energy on the bending-unbending region predicted from the geometry of the final shape and tools. The algorithm suggested was applied to a cylindrical cup deep drawing process. The blank shape and the distribution of the thickness strain are compared with those obtained from incremental finite element analysis. The comparison demonstrates the algorithm proposed reduces the difference between the results from inverse analysis and those from incremental analysis when the bending history is considered.

• Proceedings of the IEEK Conference
• /
• 2004.08c
• /
• pp.514-518
• /
• 2004

This paper presents an algorithm generating inverse element over finite fields $GF(3^{m})$, and constructing method of inverse element generator based on inverse element generating algorithm. A method computing inverse of an element over $GF(3^{m})$ which corresponds to a polynomial over $GF(3^{m})$ with order less than equal to m-l. Here, the computation is based on multiplication, square and cube method derived from the mathematics properties over finite fields.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Computational Structural Engineering
• /
• v.9 no.3
• /
• pp.145-155
• /
• 1996

This paper suggests a new algorithm of automatic mesh generation over planar domains with arbitrarily shaped boundaries and control curves. The algorithm is based on the method of recursively subdividing the domain by the path connecting, with minimum penalty value, two points on the super-loop, which consists of the boundaries and the control curves, The algorithm is not subject to any limitation on the shape of the domain, and its process can be fully automated. Therefore, this algorithm can be implemented into computer programs which require minimal user intervention while generating finite element meshes over complicated domains. This algorithm can also be easily extended for application to the generation of meshes over curved surfaces, or to the adaptive mesh generation.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.10
• /
• pp.38-45
• /
• 2004

In this paper, a parallel visualization algorithm is proposed for efficient visualization of the massive data generated from large-scale parallel finite element analysis through investigating the characteristics of parallel rendering methods. The proposed parallel visualization algorithm is designed to be highly compatible with the characteristics of domain-wise computation in parallel finite element analysis by using the sort-last-sparse approach. In the proposed algorithm, the binary tree communication pattern is utilized to reduce the network communication time in image composition routine. Several benchmarking tests are carried out by using the developed in-house software, and the performance of the proposed algorithm is investigated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.475-485
• /
• 2009

In this paper a two-step algorithm is proposed to efficiently generate rebar elements from 3D CAD data in the context of CAD/CAE data transfer. The first step is an algorithm to identify various type rebar objects and their attributes by analyzing 3D CAD data in STEP format, which is one of the international data standards. The second algorithmic step is a procedure to generate one-dimensional rebar elements from the object data made through the first step for finite element analysis or other CAE tasks. Successful rebar element data generation from real 3D CAD data for a reinforced concrete structure shows the efficacy of the proposed algorithm.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.14 no.5
• /
• pp.100-106
• /
• 2006

NiTi alloy known as its shape memory effect also has superelastic characteristic, which makes it possible to be elastic under large deformation. Since the tensile strength of the alloy is very high and density is low compared to carbon steel, it can be applied to lightweight structural design. In order to design structures with shape memory alloy, finite element analysis is used and a constitutive algorithm based on Aurrichio's model is added to LS-DYNA as a user subroutine. Explicit time integration and shell element formulation are used to simulate thin-walled structures. The algorithm uses Drucker-Prager type loading condition to calculate martensite volume fraction during the transformation. The implemented algorithm is verified in uni-axial loading condition and martensite phase transformation can be detected well with the algorithm. In this study, as a energy absorbing structure, thin-walled tube is modeled with finite elements and the deformation behavior is studied. Simulation results has shown that the martensite transformation was generated in loading condition. After plastic deformation reached, the load decreases linearly without reverse martensite transformation.