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

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

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

• Journal of the korean Society of Automotive Engineers
• /
• v.12 no.6
• /
• pp.48-59
• /
• 1990

A finite element program was developed using line elements for simulating the stretch/draw forming operation of an arbitrarily-shaped plane-strain section. An implicit, incremental, updated Lagrangian formulation is employed, introducing a minimum plastic work path assumption for each time step. Geometric and material nonlinearities are also considered within each time step. The finite element equation is based on the mesh-normal, which compatibly describes arbitrary tool surfaces and FEM meshes without depending on the explicit spatial derivatives of tool surfaces. The membrane approximation is adopted under the plane stress assumption. The sheet material is assumed to obey a rigid-viscoplastic constitutive law. The developed program was tested in the die-tryout of typical automotive inner panels. In order to determine a single friction coefficient and boundary length, FEM results and measurements of thinning for a stretched section of final die were compared. After finding analysis parameters, the sheet forming operations of original and final die designs were simulated. Excellent agreement between measured and computed thickness strains was obtained and the developed program was able to identify die designs which were rejected during die tryout.

• Journal of the Korean Geotechnical Society
• /
• v.36 no.12
• /
• pp.125-132
• /
• 2020

Embedded joints in the rock mass are a major constituent influencing its mechanical behavior. Numerical analysis requires a rigorous modeling methodology for the rock mass with detailed information regarding joint properties, orientation, spacing, and persistence. This paper provides a mechanical model for a jointed rock mass based on the implicit joint-continuum approach. Stiffness tensors for rock mass are evaluated for an assemblage of intact rock separated by sets of joint planes. It is a linear summation of compliance of each joint sets and intact rock in the serial stiffness system. In the application example, kinematic analysis for a planar failure of rock slope is comparable with empirical daylight envelope and its lateral limits. Since the developed implicit joint-continuity model is formulated on a continuum basis, it will be a major tool for the numerical simulations adopting published plenteous thermal-hydro-chemical experimental results.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.245-251
• /
• 2009

Numerical algorithms for solving the incompressible Navier-Stokes equations using P2P1 finite element are compared regarding the eigenvalues of matrices. P2P1 element allocates pressure at vertex nodes and velocity at both vertex and mid nodes. Therefore, compared to the P1P1 element, the number of pressure variables in the P2P1 element decreases to 1/4 in the case of two-dimensional problems and to 1/8 in the three-dimensional problems. Fully-implicit-integrated, semi-implicit- integrated and semi-segregated finite element formulations using P2P1 element are compared in terms of elapsed time, accuracy and eigenvlue distribution (condition number). For the comparison,they have been applied to the well-known benchmark problems. That is, the two-dimensional unsteady flows around a fixed circular cylinder and decaying vortex flow are adopted to check spatial accuracy.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.107-116
• /
• 1990

The liquid sloshing in an elastic tank is a fluid-structure interaction problem. It requires nonlinear analysis to solve the complicated physics involved in the large interaction of fluid-structure, the variation of dynamic characteristics of structure due to hydrodynamic loading, and the distorsion of fluid flow due to structural vibration. In this paper a Lagrangian FEM is introduced to analyze the liquid sloshing in an elastic tank assuming that the elastic wall is one degree of freedom rigid wall. Numerical integration is performed using an implicit-explicit algorithm, which is formed by mixing the predictor-corrector method and the Runge-Kutta 4th order method. The influence of dynamic characteristics of the sloshing tank on the fluid flow is discussed. The numerical method is also applied for the simulation of the wall generated wave in the tank.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2000.10a
• /
• pp.136-140
• /
• 2000

It is desirable to predict springback quantitatively and accurately for the tool and process design in sheet stamping operations, however, it is blown very difficult. The result of springback analysis by the finite element method is sensitively influenced by numerical factors such as blank element size, number of integration point, punch velocity, contact algorithm etc. In the present work, a parametric study by Taguchi method is performed in order to evaluate the influence of numerical factors on springback Quantitatively and to obtain the combination of numerical factors which yields the best approximation to experimental data. Since springback is determined by the residual stress after forming process, it is important to evaluate stress distribution accurately. The oscillation in the time history curve of stress obtained by explicit FEM says that the stress solution at termination time is in very unstable state. Therefore, a variability study is also carried out in this study in order to assess the stability of implicit springback analysis starting from the stress solution by explicit forming simulation. The 2D draw bending process, one of the NUMISHEET '93 benchmark problems, is adopted as an application model.

• Computers and Concrete
• /
• v.25 no.1
• /
• pp.67-73
• /
• 2020

Traditionally used analytical approach to predict the fatigue failure of reinforced concrete (RC) structure is generally conservative and has certain limitations. The nonlinear finite element method (FEM) offers less expensive solution for fatigue analysis with sufficient accuracy. However, the conventional implicit dynamic analysis is very expensive for high level computation. Whereas, an explicit dynamic analysis approach offers a computationally operative modelling to predict true responses of a structural element under periodic loading and might be perfectly matched to accomplish long life fatigue computations. Hence, this study simulates the fatigue behaviour of RC beams with finite element (FE) assemblage presenting a simplified explicit dynamic numerical solution to show computer aided fatigue behaviour of RC beam. A commercial FEM package, ABAQUS has been chosen for this complex modelling. The concrete has been modelled as a 8-node solid element providing competent compression hardening and tension stiffening. The steel reinforcements are simulated as two-node truss elements comprising elasto-plastic stress-strain behaviour. All the possible nonlinearities are duly incorporated. Time domain analysis has been adopted through an automatic Newmark-β time incremental technique. The program consists of twelve RC beams to visualize the real behaviour during fatigue process and to obtain the reliability of the study. Both the numerical and experimental results indicate a redistribution of stresses along the time and damage accumulation of beam which severely affect the serviceability and ultimate capacity of RC beam. The output of the FEM analysis demonstrates good match with the experimental consequences which affirm the efficacy of the computer aided model. The controlled fatigue damage evolution at service fatigue load limits makes the FE model an efficient tool in predicting high cycle fatigue behaviour of RC structures.