• Transactions of Materials Processing
• /
• v.9 no.3
• /
• pp.274-283
• /
• 2000

A dynamic explicit finite element code for simulating sheet forming processes has been developed. The code utilizes the discrete Kirchhoff shell element and contact force is treated by a conventional penalty method. In order to reduce the computational cost, a new and robust contact searching algorithm has been developed and implemented into the code. In the method, a hierarchical structure of tool segments is built for each tool at the initial stage of the analysis. hierarchical structure is built in a way to divide a box to 8 sub-boxes, 2 in each direction, until the lowest level of the hierarchical structure contains exactly one segment of the tool or empty. Then at each time step, contact is checked from the box to sub-boxes hierarchically for each node. Comparisons of computational results of various examples with the existing ones show the validity of the method.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1996.04a
• /
• pp.546-550
• /
• 1996

A fully automatic computer simulation technique of axisymmetric multi-stage compound forging processes was presented in this paper. A penalty rigid-viscoplastic finite element method was employed together with an improved looping method for automatically remeshing with quadrilateral finite-elements only. An application example of six-stage axisymmetric forging processes involving one cold and two hot forging processes, two piercing processes and a sizing process was given with emphasis on automatically tracing the metal flow lines through the whole simulation.

• Transactions of Materials Processing
• /
• v.4 no.4
• /
• pp.365-374
• /
• 1995

Numerical simulation of blade forming is carried out as stretch forming by an elasto-plastic finite element method. The method adopts a Lagrangian formulation, which incorporates large deformation and rotation, with a penalty method to treat the contact boundary condition. Numerical integration is done with a directional reduced integration scheme to avoid shear locking. The numerical results demonstrates various final shapes of blades which depend on the variation of the stretching force. The strain distributions in deformed blades are also obtained with the variation of the stretching force.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.18 no.4
• /
• pp.337-350
• /
• 2014

In this paper, we consider discontinuous Galerkin finite element methods with interior penalty term to approximate the solution of nonlinear parabolic problems with mixed boundary conditions. We construct the finite element spaces of the piecewise polynomials on which we define fully discrete discontinuous Galerkin approximations using the Crank-Nicolson method. To analyze the error estimates, we construct an appropriate projection which allows us to obtain the optimal order of a priori ${\ell}^{\infty}(L^2)$ error estimates of discontinuous Galerkin approximations in both spatial and temporal directions.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2001.11a
• /
• pp.166-171
• /
• 2001

In order to improve the efficiency of the design process and the quality of the resulting design for the application-based exclusive rolling element bearings, this study propose design methodologies by using a genetic-based combinatorial optimization. By the presence of discrete variables such as the number of rolling element (standard component) and by the engineering point of views, the design problem of the rolling element bearing can be characterized by the combinatorial optimization problem as a fully discrete optimization. A genetic algorithm is used to efficiently find a set of the optimum discrete design values from the pre-defined variable sets. To effectively deal with the design constraints and the multi-objective problem, a ranking penalty method is suggested for constructing a fitness function in the genetic-based combinatorial optimization. To evaluate the proposed design method, a robust performance analyzer of ball bearing based on quasi-static analysis is developed and the computer program is applied to some design problems, 1) maximize fatigue life, 2) maximize stiffness, 3) maximize fatigue life and stiffness, of a angular contact ball bearing. Optimum design results are demonstrate the effectiveness of the design method suggested in this study. It believed that the proposed methodologies can be effectively applied to other multi-objective discrete optimization problems.

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

For the propagation of elastic waves in unbounded domains, absorbing boundary conditions at the fictitious numerical boundaries have been proposed. In this paper we focus on both first and second order paraxial boundary conditions(PBCs) in the framework of variational approximations which are based on paraxial approximations of the scalar and elastic wave equations. We propose a penalty function method for the treatment of PBCs and apply these into finite element analysis. The numerical verification of the efficiency is carried out through comparing PBCs with Lysmer-Kuhlemeyer's boundary conditions.

• Transactions of Materials Processing
• /
• v.12 no.1
• /
• pp.34-42
• /
• 2003

This paper presents a rigid-plastic finite element method to handle the frictional contact problem between two deformable bodies experiencing large deformation. The variational formulation combined with incremental quasi-static model is employed for treating the contact boundary condition. The frictional behavior of the model obeys Coulomb's law of friction. The proposed contact algorithms are classified into two categories, one for searching contacting nodes and the other for calculating contact forces at the contact surface. A slave node and master contact segment are defined using the geometric condition of finite elements on the contact interface. The penalty parameter is used to limit the penetration between contacting bodies, and the finite elements are coupled with contact boundary elements.us gates and cavity thicknesses. Through this study we have observed that the jetting is related to the die swell of material. This means that the jotting is strongly affected by the elastic flow property rather than the viscous flow property in viscoelastic characteristics of molten polymer. Different resins have different elastic properties, and elastic flow behavior depends on the shear rate of flow, i.e. injection speed. Large die swell would eliminate jetting however, the retardation of die swell would stimulate jetting. In the point of mole design, reducing the thickness ratio of cavity to gate can reduce or eliminate jetting regardless of amount of elasticity of polymer melt.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.3
• /
• pp.125-142
• /
• 2017

In this paper, development of the three-dimensional structural analysis is performed by applying FETI-local method. In the FETI-local method, the penalty term is added as a preconditioner. The OPT-DKT shell element is used in the present structural analysis. Newmark-${\beta}$ method is employed to conduct the dynamic analysis. The three-dimensional FETI-local static structural analysis is conducted. The contour and the displacement of the results are compared following the different number of sub-domains. The computational time and memory usage are compared with respect to the number of CPUs used. The three-dimensional dynamic structural analysis is conducted while applying FETI-local method. The present results show appropriate scalability in terms of the computational time and memory usage. It is expected to improve the computational efficiency by combining the advantages of the original FETI method, i.e., FETI-mixed using the mixed local-global Lagrange multiplier.

• Journal of applied mathematics & informatics
• /
• v.28 no.3_4
• /
• pp.953-966
• /
• 2010

In this paper, we develop discontinuous Galerkin methods with penalty terms, namaly symmetric interior penalty Galerkin methods to solve nonlinear parabolic equations. By introducing an appropriate projection of u onto finite element spaces, we prove the optimal convergence of the fully discrete discontinuous Galerkin approximations in ${\ell}^2(L^2)$ normed space.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1997.10a
• /
• pp.40-44
• /
• 1997

The governing functional in plastic deformation has to satisfy the incompressible condition. This incompressible condition imposed on the velocity fields can be removed by introducing either the Langrange multiplier or the penalty function into the functional. In the study two-dimensional rigid plastic FEM programs using by Lagrange multiplier and penalty function are developed. A compression of cylinder and a spike forging are simulated to compare the data of loads, local mean stresses and reductions of volume.