• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.820-830
• /
• 1989

The tire inflation and contact problem has been solved by a finite element method. The finite element formulation is derived from the equilibrium equations by the principle of virtual work in the form of an updated Lagrangian formulation for incremental analysis. Then, a contact formulation is added to the finite element formulation to calculate stress state of tire in contact with flat rigid road under the load due to the self-weight of a vehicle. In the finite element analysis, equations of effective material properties are introduced to analyze a plane strain model of the shell-like tire by considering the bending effect of reinforced steel cords. The proposed equations of effective material properties produced stress concentration around the edge of belt layers, which does not appear when other well-known equations of material properties are adopted. The result from the above algorithm demonstrates the validity of the formulation and the proposed equations for the effective elastic constants. The result fully interprets the cause of separation between belt layers by showing the stress concentration.

• Coupled systems mechanics
• /
• v.7 no.5
• /
• pp.583-610
• /
• 2018

This paper discusses the issues associated with modeling frictional contact between solid bodies undergoing large deformations. The most common model for friction on contact interfaces in solid mechanics is the Coulomb friction model, in which two distinct responses are possible: stick and slip. Handling the transition between these two phases computationally has been a source of algorithmic instability, lack of convergence and non-unique solutions, particularly in the presence of large deformations. Most computational models for frictional contact have used penalty or updated Lagrangian approaches to enforce frictional contact conditions. These two approaches, however, present some computational challenges due to conditioning issues in penalty-type implementations and the iterative nature of the updated Lagrangian formulation, which, particularly in large simulations, may lead to relatively slow convergence. Alternatively, a plasticity-inspired implementation of frictional contact has been shown to handle the stick-slip conditions in a local, algorithmically efficient manner that substantially reduces computational cost and successfully avoids the issues of instability and lack of convergence often reported with other methods (Laursen and Simo 1993). The formulation of this approach, however, has been limited to the small deformations realm, a fact that severely limited its application to contact problems where large deformations are expected. In this paper, we present an algorithmically consistent formulation of this method that preserves its key advantages, while extending its application to the realm of large-deformation contact problems. We show that the method produces results similar to the augmented Lagrangian formulation at a reduced computational cost.

• Structural Engineering and Mechanics
• /
• v.12 no.1
• /
• pp.35-50
• /
• 2001

A numerical methodology is presented in this paper for the geometrically non-linear analysis of slender uni-dimensional structural elements under unilateral contact constraints. The finite element method together with an updated Lagrangian formulation is used to study the structural system. The unilateral constraints are imposed by tensionless supports or foundations. At each load step, in order to obtain the contact regions, the equilibrium equations are linearized and the contact problem is treated directly as a minimisation problem with inequality constraints, resulting in a linear complementarity problem (LCP). After the resulting LCP is solved by Lemke's pivoting algorithm, the contact regions are identified and the Newton-Raphson method is used together with path following methods to obtain the new contact forces and equilibrium configurations. The proposed methodology is illustrated by two examples and the results are compared with numerical and experimental results found in literature.

• Transactions of Materials Processing
• /
• v.14 no.1 s.73
• /
• pp.43-48
• /
• 2005

The production of metal strip with uniform thickness and flatness are two important problems associated thin strip rolling. The thickness and flatness of strip are affected by the flattening of contact surface between strip and roll, the elastic recovery and the bending of roll. Especially, the flatness of the strip is greatly affected by bending deflection of roll. The roll must be designed considered the elastic deformation of roll. This study describes the measurement of thickness and flatness of strip and shows the crown roll for producing flat strip. But it is difficult to produce the crown roller. The cross rolling that is a simple method which can produce the flat strip is introduced and it is found the optimal cross angle for improvement of flatness of plate. These problems are solved by the MARC code on the basis of elastic-plastic material and the updated Lagrangian formulation.

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

• Transactions of Materials Processing
• /
• v.4 no.3
• /
• pp.222-232
• /
• 1995

The planar anisotripic FEM analysis for predicting earing profiles and draw-in amounts in the deep-drawing process is introduced. An implicit, incremental, updated Lagrangian formulation with a rigid-viscoplastic constitutive equation is employed. Contact and friction are considered through the mesh-based unit vector and normal contact pressure. The consistent full set of governing relations, which is comprising euilbrium and geometric constraint equations, is appropriately linearized. Barlat's strain-rate potential is employed, whose in-plane anisotropic properties are taken into account with anisotropic coefficients and potential parameters. The linear triangular membrane elements are used for depicting the formed sheet. In the numerical simulations of deep drawing processes of a flat-top cylindrical cup for 2090-T3 aluminum alloy sheet show good agreement with experiments, although some discrepancies were observed in the directional trend of cup height and thickness strains.

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

The plane strain analysis for simulating the stretch/draw forming operation of arbitrarily-shaped tool profiles and arbitrarily draw-in conditions is introduced. An implicit, incremental, updated Lagrangian formulation is employed, introducing a rigid-viscoplastic constitutive equation. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshe without depending on the explicit spatial derivatives of tool surfaces. The FEM formulation is tested in the sections automotive inner panel and two-side draw-in. Not only the excellent agreement between measured and computed strains in the stretched section is obtained, but also the numerical stability of current formulation is verified in the two-side draw-in section.

• Transactions of Materials Processing
• /
• v.1 no.1
• /
• pp.95-103
• /
• 1992

The plane strain analysis for simulating the stretch/draw forming operation with an arbitrarily-shaped tool profile is introduced. An implicit, incremental, updated Lagrangian formulation with a rigid-viscoplastic constitutive equation is employed. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshes without depending on the explicit spatial derivatives of tool surfaces. The linear line elements are used for depicting the formed sheet, based on membrane approximation. The FEM formulation is tested in the sections of automotive inner panel and two-side draw-in. Not only the excellent agreement between measured and computed strains is obtained in the stretched section, but also the numerical stability of formulation is verified in the draw-in section.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1994.03a
• /
• pp.118-128
• /
• 1994

The planar anisotropic FEM analysis for predicting the earing profiles and draw-in amounts in the deep-drawing processes is introduced. An implicit, incremental, updated Lagrangian formulation with a rigid-viscoplastic constitutive equation is employed. Contact and friction are considered through the mesh-based unit vectors and the normal contact pressure. the consistent full set of governing relations, comprising equilibrium and geometric constraint equations, is appropriately linearized. Barlat's strain-rate potential is employed, whose in-plane anisotropic properties are taken into account with anisotropic coefficients and potential parameter. The linear triangular membrane elements are used for depicting the formed sheet. with the numerical simulations of deep drawing processes of flat-top cylindrical cup for the 2090-T3 aluminum effects on the earing behavior are examined. Earing predictions made for the 2090-T3 aluminum alloy sheet show good agreement with experiments, although some discrepancies were observed in the directional trend of cup height and thickness strains.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.8
• /
• pp.2113-2122
• /
• 1994

The 3-D FEM analysis for simulating the stamping operation of planar anisotropic sheet metals with arbitrarily-shaped tools is introduced. An implicit, incremental, updated Lagrangian formulation with a rigid-viscoplastic constitutive equation is employed. Contact and friction are considered through the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshes without depending on the explicit spatial derivatives of tool surfaces. The consistent full set of governing relations, comprising equilibrium equation and mesh-normal geometric constraints, is appropriately linearized. The linear triangular elements are used for depicting the formed sheet, based on membrane approximation. Barlat's non-quadratic anisotropic yield criterion(strain-rate potential) is employed, whose in-plane anisotropic properties are taken into account with anisotropic coefficients and non-quadratic function parameter. The planar anisotropic finite element formulation is tested with the numerical simulations of the stamping of an automotive hood inner panel and the drawing of a hemispherical punch. The in-plane anisotropic effects on the formability of both mild steel and aluminum alloy sheet metals are examined.