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

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.5 s.176
• /
• pp.1075-1083
• /
• 2000

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method. The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.7
• /
• pp.1042-1049
• /
• 1997

LMFBR(Liquid Metal Fast Breeder Reactor) vessel is operated under the high temperatures of 500-550.deg. C. Thus, transient thermal loads were severe enough to cause inelastic deformation due to creep-fatigue and plasticity. For reduction of such inelastic deformations, Y-piece structure in the form of a thermal sleeve is used in LMFBR vessel under repeated start-up, service and shut-down conditions. Therefore, a systematic method for inelastic analysis is needed for design of the Y-piece structure subjected to such loading conditions. In the present investigation, finite element analysis of heat transfer and inelastic thermal stress were carried out for the Y-piece structure in LMFBR vessel under service conditions. For such analysis, ABAQUS program was employed based on the elasto-plastic and Chaboche viscoplastic constitutive equations. Based on numerical data obtained from the analysis, creep-fatigue damage estimation according to ASME Code Case N-47 was made and compared to each other. Finally, it was found out that the numerical predictio of damage level due to creep based on Chaboche unified viscoplastic constitutive equation was relatively better compared to elasto-plastic constitutive formulation.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.11
• /
• pp.3516-3524
• /
• 1996

Service conditions for structures at elevated temperatures in nuclear power plant involve transient thermal and mechanical load levels that are severe enough to caeuse inelastic deformations due to creep and plasticity. Therefore, a systematic mehtod of inelastic analysis is needed for the design of structural components in nuclear poser plants subjected to such loading conditions. In the present investigation, the Chabodhe model, one of the unified viscoplastic constitutive equations, was selected for systematic inelastic analysis. The material response was integrated based on GMR ( generallized mid-point rule) time integral scheme and provided to ABAQUS as a material subroutine, UMAT program. By comparing results obtaned from uniaxial analysis using the developed UMAT program with those from Runge-Kutta solutions and experimentaiton, the validity of the adopted Chaboche model and the numerical stability and accuracy of the developed UMAT program were verified. In addition, the developed material subroutine was applied for uniaxial creep and tension analyses for the plate with a hole in the center. The application further demonstrates usefulness of the developed program.

• Geomechanics and Engineering
• /
• v.25 no.3
• /
• pp.195-205
• /
• 2021

Time effect on the deformation and strength characteristics of geogrid reinforced sand retaining wall has become an important issue in geotechnical and transportation engineering. Three physical model tests on geogrid reinforced sand retaining walls performed under various loading conditions were simulated to study their rate-dependent behaviors, using the presented nonlinear finite element method (FEM) analysis procedure. This FEM was based on the dynamic relaxation method and return mapping scheme, in which the combined effects of the rate-dependent behaviors of both the backfill soil and the geosynthetic reinforcement have been included. The rate-dependent behaviors of sands and geogrids should be attributed to the viscous property of materials, which can be described by the unified three-component elasto-viscoplastic constitutive model. By comparing the FEM simulations and the test results, it can be found that the present FEM was able to be successfully extended to the boundary value problems of geosynthetic reinforced soil retaining walls. The deformation and strength characteristics of the geogrid reinforced sand retaining walls can be well reproduced. Loading rate effect, the trends of jump in footing pressure upon the step-changes in the loading rate, occurred not only on sands and geogrids but also on geogrid reinforced sands retaining walls. The lateral earth pressure distributions against the back of retaining wall, the local tensile force in the geogrid arranged in the retaining wall and the local stresses beneath the footing under various loading conditions can also be predicted well in the FEM simulations.

• Transactions of Materials Processing
• /
• v.26 no.5
• /
• pp.286-292
• /
• 2017

Magnesium alloy shows strong anisotropy and asymmetric behavior in tension and compression curve, especially at room temperature. These characteristics limit the application of finite element method (FEM) which is based on conventional continuum mechanics. To accurately predict the material behavior of magnesium alloy at microstructural level, a methodology of fully coupled multiscale simulation is presented and a crystal plasticity model as a constitutive equation in the simulation of metal forming process is introduced in this study. The existing constitutive equation for rigid plastic FEM is modified to accommodate deviatoric stress component and its derivatives with respect to strain rate components. Viscoplastic self-consistent (VPSC) polycrystal model was selected as a constitutive model because it was regarded as the most robust model compared to Taylor model or Sachs model. Stiffness matrix and load vector were derived based on the new approach and implemented into $DEFORM^{TM}-3D$ via a user subroutine handling stiffness matrix at an elemental level. The application to extrusion and rolling process of pure magnesium is presented in this study to assess the validity of the proposed multiscale process.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.6
• /
• pp.1411-1421
• /
• 1995

A rigid visco-plastic finite element method has been developed for modeling superplastic forming processes. The optimum pressure-time relationship for a target strain rate and thickness distributions was predicted using two-node line element based on membrane approximation for plane strain and axisymmetric condition. Analysis of superplastic forming was carried out using the developed program and the numerical results were compared to the values available in the literature for plane strain problems. For description of the contact between the dies and sheet, the direct projection method was applied to the complicated problem and the validity of the scheme was tested. Experiments for the various geometries such as hemisphere and cone were performed with the developed forming machine using the calculated optimum pressure-time curves. Comparison between analysis and experiments showed good agreement.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.46 no.4
• /
• pp.75-84
• /
• 2004

The evaluation technique of seismic performance on agricultural infrastructure based on dynamic numerical simulations, which Included a cyclic elasto-plastic and a viscoelastic-viscoplastic constitutive model to actual multi-layered ground conditions during large earthquake were performed by a liquefaction analysis in the present study. From the liquefaction analysis, it was verified that the models can give a good description of the damping characteristics and liquefaction phenomena of ground accurately during large event which induces plastic deformation in large strain range.

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

A full three-dimensional thermo-coupled rigid-viscoplastic finite element method and the currently developed microstructural evolution system which includes semi-empirical mathematical equations suggested by different research groups were used together to form an integrated system of process and microstructure simulation of hot rolling. The distribution and time history of thermomechanical variables such as temperature, strain, strain rate, and time during pass and between passes were obtained FEM analysis of multipass hot rolling processes. Then distribution of metallurgical variables were calculated successfully on the basis of instantaneous thermomechanical data. For the verification of this method the evolution of microstructure in plate rolling and shape rolling was simulated and their results were compared with the data available in literature. Consequently, this approach makes it passible to describe the realistic evolution of microstructure by avoiding the use of erroneous average value and can be used in CAE of multipass hot rolling.

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

The purpose of this research is for the development of a new type forging equipment H.C.G(Hyundai Continuous Grain-flow) by using two virtual build-up tools rigid viscoplastic FEM and downsized plasticine experiment. This forging equipment consists of consecutive horizontal and vertical pressure while the traditional forging method consists of only vertical pressure. Using this method high quality crankshafts can be forged as it can maintain a continuous grain flow. The factors considered in the development of equipment are die geometry for flawless deformed shape die reaction forces stress/strain distributions and continuous material flow. We carried out several numerical simulations and downsized plasticine experiments for the proper design of the forging equipment. The validity of those simulation results is confirmed by checking with the actual test results. Based on these simulation results the proper design of the H.C.G for ging equipment is enabled.