A finite element formulation based on non-associated plasticity for sheet metal forming

In the present paper, a finite element formulation based on non-associated plasticity is developed. In the constitutive formulation, isotropic hardening is assumed and an evolution equation for the hardening parameter consistent with the principle of plastic work equivalence is introduced. The yield function and plastic potential function are considered as two different functions with functional form as the yield function of Hill [Hill, R., 1948. Theory of yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. A 193, 281–297] or Karafillis–Boyce associated model [Karafillis, A.P. Boyce, M., 1993. A general anisotropic yield criterion using bounds and a transformation weighting tensor. J. Mech. Phys. Solids 41, 1859–1886]. Algorithmic formulations of constitutive models that utilize associated or non-associated flow rule coupled with Hill or Karafillis–Boyce stress functions are derived by application of implicit return mapping procedure. Capabilities in predicting planar anisotropy of the Hill and Karafillis–Boyce stress functions are investigated considering material data of Al2008-T4 and Al2090-T3 sheet samples. The accuracy of the derived stress integration procedures is investigated by calculating iso-error maps.     

Evaluation of advanced anisotropic models with mixed hardening for general associated and non-associated flow metal plasticity

The main objective of this paper is to develop a generalized finite element formulation of stress integration method for non-quadratic yield functions and potentials with mixed nonlinear hardening under non-associated flow rule. Different approaches to analyze the anisotropic behavior of sheet materials were compared in this paper. The first model was based on a non-associated formulation with both quadratic yield and potential functions in the form of Hill’s (1948). The anisotropy coefficients in the yield and potential functions were determined from the yield stresses and r-values in different orientations, respectively. The second model was an associated non-quadratic model (Yld2000-2d) proposed by Barlat et al. (2003). The anisotropy in this model was introduced by using two linear transformations on the stress tensor. The third model was a non-quadratic non-associated model in which the yield function was defined based on Yld91 proposed by Barlat et al. (1991) and the potential function was defined based on Yld89 proposed by Barlat and Lian (1989). Anisotropy coefficients of Yld91 and Yld89 functions were determined by yield stresses and r-values, respectively. The formulations for the three models were derived for the mixed isotropic-nonlinear kinematic hardening framework that is more suitable for cyclic loadings (though it can easily be derived for pure isotropic hardening). After developing a general non-associated mixed hardening numerical stress integration algorithm based on backward-Euler method, all models were implemented in the commercial finite element code ABAQUS as user-defined material subroutines. Different sheet metal forming simulations were performed with these anisotropic models: cup drawing processes and springback of channel draw processes with different drawbead penetrations. The earing profiles and the springback results obtained from simulations with the three different models were compared with experimental results, while the computational costs were compared. Also, in-plane cyclic tension–compression tests for the extraction of the mixed hardening parameters used in the springback simulations were performed for two sheet materials.     

Model identification and FE simulations: Effect of different yield loci and hardening laws in sheet forming

The bi-axial experimental equipment [Flores, P., Rondia, E., Habraken, A.M., 2005a. Development of an experimental equipment for the identification of constitutive laws (Special Issue). International Journal of Forming Processes] developed by Flores enables to perform Bauschinger shear tests and successive or simultaneous simple shear tests and plane strain tests. Flores investigates the material behavior with the help of classical tensile tests and the ones performed in his bi-axial machine in order to identify the yield locus and the hardening model. With tests performed on one steel grade, the methods applied to identify classical yield surfaces such as [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic materials. Proceedings of the Royal Society of London A 193, 281–297; Hosford, W.F., 1979. On yield loci of anisotropic cubic metals. In: Proceedings of the 7th North American Metalworking Conf. (NMRC), SME, Dearborn, MI, pp. 191–197] ones as well as isotropic Swift type hardening, kinematic Armstrong–Frederick or Teodosiu and Hu hardening models are explained. Comparison with the Taylor–Bishop–Hill yield locus is also provided. The effect of both yield locus and hardening model choices is presented for two applications: plane strain tensile test and Single Point Incremental Forming (SPIF).     

Incorporation of yield surface distortion in finite deformation constitutive modeling of rigid-plastic hardening materials based on the Hencky logarithmic strain

In this paper a constitutive model for rigid-plastic hardening materials based on the Hencky logarithmic strain tensor and its corotational rates is introduced. The distortional hardening is incorporated in the model using a distortional yield function. The flow rule of this model relates the corotational rate of the logarithmic strain to the difference of the Cauchy stress and the back stress tensors employing deformation-induced anisotropy tensor. Based on the Armstrong–Fredrick evolution equation the kinematic hardening constitutive equation of the proposed model expresses the corotational rate of the back stress tensor in terms of the same corotational rate of the logarithmic strain. Using logarithmic, Green–Naghdi and Jaumann corotational rates in the proposed constitutive model, the Cauchy and back stress tensors as well as subsequent yield surfaces are determined for rigid-plastic kinematic, isotropic and distortional hardening materials in the simple shear deformation. The ability of the model to properly represent the sign and magnitude of the normal stress in the simple shear deformation as well as the flattening of yield surface at the loading point and its orientation towards the loading direction are investigated. It is shown that among the different cases of using corotational rates and plastic deformation parameters in the constitutive equations, the results of the model based on the logarithmic rate and accumulated logarithmic strain are in good agreement with anticipated response of the simple shear deformation.     

金属板料成形的快速有限元分析

研究了基于形变理论的金属板形成快速有限元分析的方法－反向方法,并实现了计算程序。通过实例的计算结果和实验以及增量方法进行了比较,表明此方法能够定性地分析成形工件中的变形情况。由于计算速度快、建立分析模模型简单,此方法可用于设计早期估计零件的可成形性,以及部分工艺参数对形的影响。     

Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheets: Application to sheet springback

The constitutive model for the unusual asymmetric hardening behavior of magnesium alloy sheet presented in a companion paper (Lee, M.G., Wagoner, R.H., Lee, J.K., Chung, K., Kim, H.Y., 2008. Constitutive modeling for anisotropic/asymmetric hardening behavior of magnesium alloy sheet, Int. J. Plasticity 24(4), 545–582) was applied to the springback prediction in sheet metal forming. The implicit finite element program ABAQUS was utilized to implement the developed constitutive equations via user material subroutine. For the verification purpose, the springback of AZ31B magnesium alloy sheet was measured using the unconstrained cylindrical bending test of Numisheet (Numisheet ’2002 Benchmark Problem, 2002. In: Yang, D.Y., Oh, S.I., Huh, H., Kim, Y.H. (Eds.), Proceedings of 5th International Conference and Workshop on Numerical Simulation of 3D Sheet Forming Processes, Jeju, Korea) and 2D draw bend test. With the specially designed draw bend test the direct restraining force and long drawn distance were attainable, thus the measurement of the springback could be made with improved accuracy comparable with conventional U channel draw bend test. Besides the developed constitutive models, other models based on isotropic constitutive equations and the Chaboche type kinematic hardening model were also considered. Comparisons were made between simulated results by the finite element analysis and corresponding experiments and the newly proposed model showed enhanced prediction capability, which was also supported by the simple bending analysis adopting asymmetric stress–strain response.     

FE-analysis of sheet metal forming processes using continuous contact treatment

Discrete meshes cause stepwise propagation of the contact nodes of a sheet despite the fact that the contact region in the actual forming process is altered very smoothly. This can cause problems of convergence and accuracy in contact-sensitive processes, such as a bending process. In this study, a scheme for a continuous contact treatment is proposed in order to consider the more realistic behavior of the contact phenomena during the forming process. For verification of the proposed method, the contact pressures and forming load are evaluated during the compression forming of a tube. The analysis of a hemispherical dome formed without a blank holder is also presented in order to investigate the effects of the proposed algorithm. The results show that the precise deformation mode is predicted by the utilization of the proposed method.     

Advances in experiments on metal sheets and tubes in support of constitutive modeling and forming simulations

This work is a review of experimental methods for observing and modeling the anisotropic plastic behavior of metal sheets and tubes under a variety of loading paths, such as biaxial compression tests; biaxial tension tests on metal sheets and tubes using closed-loop electrohydraulic testing machines; the abrupt strain path change method for detecting a yield vertex and subsequent yield loci without unloading; in-plane stress reversal tests on metal sheets; and multistage tension tests. Observed material responses are compared with the predictions of phenomenological plasticity models. Special attention is paid to the plastic deformation behavior of materials commonly used in industry, and to verifying the validity of conventional anisotropic yield criteria for those materials and associated flow rules at large plastic strains. The effects of using appropriate anisotropic yield criteria on the accuracy of simulations of forming defects, such as large springback and fracture, are also presented to highlight the importance of accurate material testing and modeling.     

An Eulerian multiplicative constitutive model of finite elastoplasticity

An Eulerian rate-independent constitutive model for isotropic materials undergoing finite elastoplastic deformation is formulated. Entirely fulfilling the multiplicative decomposition of the deformation gradient, a constitutive equation and the coupled elastoplastic spin of the objective corotational rate therein are explicitly derived. For the purely elastic deformation, the model degenerates into a hypoelastic-type equation with the Green–Naghdi rate. For the small elastic- and rigid-plastic deformations, the model converges to the widely-used additive model where the Jaumann rate is used. Finally, as an illustration, using a combined exponential isotropic-nonlinear kinematic hardening pattern, the finite simple shear deformation is analyzed and a comparison is made with the experimental findings in the literature.     

Anisotropic finite elastoplasticity with nonlinear kinematic and isotropic hardening and application to sheet metal forming

The paper discusses the derivation and the numerical implementation of a finite strain material model for plastic anisotropy and nonlinear kinematic and isotropic hardening. The model is derived from a thermodynamic framework and is based on the multiplicative split of the deformation gradient in the context of hyperelasticity. The kinematic hardening component represents a continuum extension of the classical rheological model of Armstrong–Frederick kinematic hardening. Introducing the so-called structure tensors as additional tensor-valued arguments, plastic anisotropy can be modelled by representing the yield surface and the plastic flow rule as functions of the structure tensors. The evolution equations are integrated by a new form of the exponential map that preserves plastic incompressibility and uses the spectral decomposition to evaluate the exponential tensor functions in closed form. Finally, the applicability of the model is demonstrated by means of simulations of several deep drawing processes and comparisons with experiments.     

Effects of spatial grain orientation distribution and initial surface topography on sheet metal necking

The finite element method is used to numerically simulate localized necking in AA6111-T4 under stretching. The measured EBSD data (grain orientations and their spatial distributions) are directly incorporated into the finite element model and the constitutive response at an integration point is described by the single crystal plasticity theory. We assume that localized necking is associated with surface instability, the onset of unstable growth in surface roughening. It is demonstrated that such a surface instability/necking is the natural outcome of the present approach, and the artificial initial imperfection necessitated by the macroscopic M–K approach [Marciniak and Kuczynski (1967). Int. J. Mech. Sci. 9, 609–620] is not relevant in the present analysis. The effects of spatial orientation distribution, material strain rate sensitivity, texture evolution, and initial surface topography on necking are discussed. It is found that localized necking depends strongly on both the initial texture and its spatial orientation distribution. It is also demonstrated that the initial surface topography has only a small influence on necking.     

Micromechanical analysis of kinematic hardening in natural clay

This paper presents a micromechanical analysis of the macroscopic behaviour of natural clay. A microstructural stress–strain model for clayey material has been developed which considers clay as a collection of clusters. The deformation of a representative volume of the material is generated by mobilizing and compressing all the clusters along their contact planes. Numerical simulations of multistage drained triaxial stress paths on Otaniemi clay have been performed and compared the numerical results to the experimental ones in order to validate the modelling approach. Then, the numerical results obtained at the microscopic level were analysed in order to explain the induced anisotropy observed in the clay behaviour at the macroscopic level. The evolution of the state variables at each contact plane during loading can explain the changes in shape and position in the stress space of the yield surface at the macroscopic level, as well as the rotation of the axes of anisotropy of the material.     

Investigation of advanced strain-path dependent material models for sheet metal forming simulations

Sheet metal forming processes often involve complex loading sequences. To improve the prediction of some undesirable phenomena, such as springback, physical behavior models should be considered. This paper investigates springback behavior predicted by advanced elastoplastic hardening models which combine isotropic and kinematic hardening and take strain-path changes into account. A dislocation-based microstructural hardening model formulated from physical observations and the more classical cyclic model of Chaboche have been considered in this work. Numerical implementation was carried out in the ABAQUS software using a return mapping algorithm with a combined backward Euler and semi-analytical integration scheme of the constitutive equations. The capability of each model to reproduce transient hardening phenomena at abrupt strain-path changes has been shown via simulations of sequential rheological tests. A springback analysis of strip drawing tests was performed in order to emphasize the impact of several influential parameters, namely: process, numerical and behavior parameters. The effect of the two hardening models with respect to the process parameters has been specifically highlighted.     

Finite element modeling of elasto-plastic contact between rough surfaces

This paper presents a finite element calculation of frictionless, non-adhesive, contact between a rigid plane and an elasto-plastic solid with a self-affine fractal surface. The calculations are conducted within an explicit dynamic Lagrangian framework. The elasto-plastic response of the material is described by a J₂ isotropic plasticity law. Parametric studies are used to establish general relations between contact properties and key material parameters. In all cases, the contact area A rises linearly with the applied load. The rate of increase grows as the yield stress σ_y decreases, scaling as a power of σ_y over the range typical of real materials. Results for A from different plasticity laws and surface morphologies can all be described by a simple scaling formula. Plasticity produces qualitative changes in the distributions of local pressures in the contact and of the size of connected contact regions. The probability of large local pressures is decreased, while large clusters become more likely. Loading-unloading cycles are considered and the total plastic work is found to be nearly constant over a wide range of yield stresses.     

Solution existence conditions for elastoplastic constitutive models of granular materials

In this paper, the conditions of solution existence for stress rates under given strain rates are investigated. The focus of the solution existence investigation is on the non-associated flow rule and elastic stress–strain relationship. Granular materials characterized with strong non-associated plastic flows are used as a particular example for analysis. Various flow rules for granular materials are analyzed, including Rowe’s, Roscoe’s flow rules and their modified versions. In the elastic stress–strain relationships of materials, the effects of Poisson’s ratio on solution existence are investigated. Both isotropic and anisotropic elasticity are considered. Given a granular material and its states, it is found that there exists a critical Poisson’s ratio for a particular non-associated flow rule. When the Poisson’s ratio of a material is above this critical Poisson’s ratio, its constitutive model is susceptible to solution non-existence. It is suggested that special attentions should be paid to the selection of material Poisson’s ratio and non-associated flow rule to ensure the existence of elastoplastic solutions.