• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.6 s.249
• /
• pp.705-712
• /
• 2006

A finite element method based on the penalized subdomain-interface framework is proposed for fully-coupled, nonlinear thermomechanical analyses with thermal contact anuor radiation boundaries. In the variational formulation, a well-known penalty functional scheme is adopted for connecting subdomains and interfaces that satisfy various continuity requirements. As a logical consequence, the whole domain can be arbitrarily divided into independently-modeled subdomains without considering the conformity of meshes along their interfaces. Since the nonlinearities due to the contact and radiation boundaries can be localized within a few subdomains, the computational efficiency of the present method is greatly increased with appropriate solution algorithms. By solving some numerical problems, these advantageous features are confirmed carefully.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.03b
• /
• pp.253-256
• /
• 1999

The Streamline Upwind Petrov-Galerkin(SUPG) finite element method is used to solve the two-dimensional laminar and turbulent flow. The flow is simulated by averaged Navier-Stokes equations with a penalty function approach and the lograithmic(k-$\varepsilon$) turbulent model is employed to take into account its turbulent behavior. The near-wall viscous sub-layer model is employed to approach the dominant viscous effects in the near wall zones. To find a good-enough initial guess of the Newton-Raphson iteration solving Nonlinear Matrix the Incremental method is considered for momentum and the Incomplete logarithmic turbu-lent equations for Turbulence. The validation of our method is investigated in comparision with published experimental data.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1995.10a
• /
• pp.102-108
• /
• 1995

Internal contact scheme between two free surfaces on one deforming body has been proposed by using the penalty method. It has been imposed to be internal boundary condition on two-dimensional thermo-viscoplastic finite element method so as to analyze one deforming body, which has two free surfaces penetrating each others. Analysis of side pressing with a circular void and a inclined elliptic hole have been carried out in order to verity the proposed contact scheme. A finite element code imposed internal boundary condition has been applied to two-dimensional analysis of free forging of large ingot with a void. Through the analysis, effects of working parameters in order to consolidate voids have been investigated.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.03b
• /
• pp.22-25
• /
• 1999

A dynamic explicit finite element code for simulating sheet forming processes has been developed The code utilises 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 called a tree structure is built for each tool at the initial stage of the analysis Tree is built in a way to divide a trunk to 8 sub-trunk 2 in each direction until the lowest level of the tree(leaf) contains exactly one segment of the tool. In order to have a well-balanced tree each box on each sub level contains one eighth of the segments. Then at each time step contact line from a node comes out of the surface of the tool. Simulation of various sheet forming processes were performed to verify the validity of the developed code with main focus on he usefulness of the developed contact searching algorithm.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.03b
• /
• pp.103-107
• /
• 1999

In this paper, a computer simulation technique for the forging process having a floating die is presented. The penalty rigid-plastic finite element method is employed together with an iteratively force-balancing method, in which the convergence is achieved when the floating die part is in force equilibrium within the user-specified tolerance. The force balance is controled by adjusting the velocity of the floating die in an automatic manner. An application example of a three-stage cold forging process is given.

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.39-46
• /
• 1999

The FVM(Finite Volume Method) have been used mainly for the flow analyses in the piston-cylinder. The objective of the present study is to analyze numerically the piston-driven intake flows using the FEM(Finite Element Method). The FEM algorithm used in this study is 4-step time-splitting method which requires much less execution time and computer storage than the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the momentum equations to prevent checkerboard pressure oscillations. Also, the ALE(arbitrary Lagrangian Eulerian) method is adopted for the moving grids. The calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.

• Korean Journal of Computational Design and Engineering
• /
• v.2 no.2
• /
• pp.71-76
• /
• 1997

The deformations of polycrystalline materials are modelled by linking a constitutive equation for the crystallographic slip of a single crystal to the macroscopic behavior of the aggregate. In this study, anisotropic elasticity (lattice stretching) of a cubic crystal is incoporated into the anisotropic plasticity from crystallographic slip. The constitutive description for the aggregate, derived from a crystal plasticity theory, is used to formulate a Consistent Penalty Finite Element Method for the anisotropic elasto-viscoplastic deformation of polycrystalline materials. As an application, a plane-strain forging process is simulated and the effects of the initial textures on the deformation behavior of the workpiece are examined.

• Korean Journal of Mathematics
• /
• v.31 no.4
• /
• pp.419-431
• /
• 2023

In this paper, we propose an iterative method for a symmetric interior penalty Galerkin method for heterogeneous elliptic problems. The iterative method consists mainly of two parts based on a non-overlapping domain decomposition approach. One is an intermediate preconditioner constructed by understanding the properties of the discontinuous finite element functions and the other is a preconditioning related to the dual-primal finite element tearing and interconnecting (FETI-DP) methodology. Numerical results for the proposed method are presented, which demonstrate the performance of the iterative method in terms of various parameters associated with the elliptic model problem, the finite element discretization, and non-overlapping subdomain decomposition.

• Bulletin of the Korean Mathematical Society
• /
• v.49 no.4
• /
• pp.829-850
• /
• 2012

In this paper, we adopt symmetric interior penalty discontinuous Galerkin (SIPG) methods to approximate the solution of nonlinear viscoelasticity-type equations. We construct finite element space which consists of piecewise continuous polynomials. We introduce an appropriate elliptic-type projection and prove its approximation properties. We construct semidiscrete discontinuous Galerkin approximations and prove the optimal convergence in $L^2$ normed space.

• Transactions of Materials Processing
• /
• v.6 no.5
• /
• pp.386-394
• /
• 1997

The analysis of the springback is done based on the stress of sheet after forming. Therfore, it is important to get the accurate stress from forming analysis. In this study, some parameters that influence on the accuracy of the springback estimation are investigated. Discretization of sheet and tools, choice of penalty constant and damping in contact treatment, and tool speed scaling are chosen as parameters. As a numerical example, the 2D draw bending benchmark problem of the NUMISHEET'93 is used. Also, the springback results of the s-rail benchmark problem of the NUMISHEET'96 are presented.