A plane stress anisotropic plastic flow theory for orthotropic sheet metals

A phenomenological theory is presented for describing the anisotropic plastic flow of orthotropic polycrystalline aluminum sheet metals under plane stress. The theory uses a stress exponent, a rate-dependent effective flow strength function, and five anisotropic material functions to specify a flow potential, an associated flow rule of plastic strain rates, a flow rule of plastic spin, and an evolution law of isotropic hardening of a sheet metal. Each of the five anisotropic material functions may be represented by a truncated Fourier series based on the orthotropic symmetry of the sheet metal and their Fourier coefficients can be determined using experimental data obtained from uniaxial tension and equal biaxial tension tests. Depending on the number of uniaxial tension tests conducted, three models with various degrees of planar anisotropy are constructed based on the proposed plasticity theory for power-law strain hardening sheet metals. These models are applied successfully to describe the anisotropic plastic flow behavior of 10 commercial aluminum alloy sheet metals reported in the literature.     

Endochronic description of plastic anisotropy in sheet metal

It is shown in this paper that an extended form of Hills quadratic yield criterion for anisotropic sheet metal can be derived from an endochronic theory of plasticity. The extended form considers the combined isotropic–kinematic hardening and the anomalous behavior observed in the anisotropic plastic behavior of sheet metals can be accounted for by the concept of kinematic hardening.This form of anisotropic endochronic theory can accommodate the usual requirement of normality between the plastic strain rate and the yield function. In addition, the theory leads naturally to the expressions for back stresses. This work provides an additional example to show that the form of the intrinsic time is directly related to the form of the yield function.It is suggested that the coefficients of the quadratic yield function be determined from the yield stresses obtained from a set of tension tests.     

A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming

In this paper an anisotropic material model based on non-associated flow rule and mixed isotropic–kinematic hardening was developed and implemented into a user-defined material (UMAT) subroutine for the commercial finite element code ABAQUS. Both yield function and plastic potential were defined in the form of Hill’s [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. Lond. A 193, 281–297] quadratic anisotropic function, where the coefficients for the yield function were determined from the yield stresses in different material orientations, and those of the plastic potential were determined from the r-values in different directions. Isotropic hardening follows a nonlinear behavior, generally in the power law form for most grades of steel and the exponential law form for aluminum alloys. Also, a kinematic hardening law was implemented to account for cyclic loading effects. The evolution of the backstress tensor was modeled based on the nonlinear kinematic hardening theory (Armstrong–Frederick formulation). Computational plasticity equations were then formulated by using a return-mapping algorithm to integrate the stress over each time increment. Either explicit or implicit time integration schemes can be used for this model. Finally, the implemented material model was utilized to simulate two sheet metal forming processes: the cup drawing of AA2090-T3, and the springback of the channel drawing of two sheet materials (DP600 and AA6022-T43). Experimental cyclic shear tests were carried out in order to determine the cyclic stress–strain behavior and the Bauschinger ratio. The in-plane anisotropy (r-value and yield stress directionalities) of these sheet materials was also compared with the results of numerical simulations using the non-associated model. These results showed that this non-associated, mixed hardening model significantly improves the prediction of earing in the cup drawing process and the prediction of springback in the sidewall of drawn channel sections, even when a simple quadratic constitutive model is used.     

Comparison of the work-hardening of metallic sheets using tensile and shear strain paths

This work deals with the characterization of the kinematic work-hardening of a bake-hardening steel. A shear test device has been designed and its use for the characterization of the work-hardening of sheet metals is described. Two main results are presented. Firstly, a local strain measurement, based on the following of three dots drawn on the gauge area, gives the evolution of the strain tensor eigenvalues during the test. It is shown, by comparing the theoretical kinematics of simple shear with a slightly perturbated one, that the strain state is close to the ideal one in the center of the gauge area. Secondly, reversal of the shear direction is performed after several prestrain and the evolution of the kinematic work-hardening with the equivalent plastic strain has been identified using an anisotropic elasto-viscoplastic model of Hill 1948 type. Isotropic and kinematic contributions of the work-hardening are also calculated from loading–unloading tensile tests and are compared to those obtained from the simple shear tests. The results show a discrepancy between both identification for the isotropic and the kinematic hardening. However, they are in agreement concerning the evolution of the global work-hardening.     

A novel quadratic yield model to describe the feature of multi-yield-surface of rolled sheet metals

Since Hill’s quadratic yield model [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. Roy. Soc. Lond. A193, 281–297] cannot address enough experimental results for fairly describing the “anomalous” yield behavior as observed in some of rolled sheet metals, a new quadratic yield model is proposed. As the concept of multiple yielding systems is introduced into the new quadratic yield model, seven commonly used experimental results, three uniaxial tension stresses, one equibiaxial tension stress and three strain ratios, can all be taken into account for characterizing the anisotropy of rolled sheet metals. If more experimental results are extra needed for further improving the prediction, this yield model is still workable. As the experimental parameters are defined as functions of loading direction of corresponding test separately from the major part of yield model, the increase of experimental results regarding the same test does not vary the quadratic form of yield model. The representation of this yield model with axes of principal stresses demonstrates the similar form to Hill’s quadratic model. Therefore, many previous studies developed from Hill’s quadratic yield model can be directly upgraded by the new model to reach a higher accurate level.     

Analysis of sheet metal formability through isotropic and kinematic hardening models

The present paper aims at analysing the sheet metal formability through several isotropic and kinematic hardening models. Specifically, a special attention is paid to the physically-based hardening model of Teodosiu and Hu (1995), which accounts for the anisotropic work-hardening induced by the microstructural evolution at large strains, as well as to some more conventional hardening models, including the isotropic Swift strain-hardening power law, and the Voce saturation strain-hardening law, combined with a non-linear kinematic hardening described by the Armstrong–Frederick law. The onset of localized necking is simulated by an advanced sheet metal forming limit model which connects, through the Marciniak–Kuczinsky analysis, the hardening models with the anisotropic yield criterion Yld2000-2d (Barlat et al., 2003). Both linear and complex strain paths are taken into account. The selected material is a DC06 steel sheet. The validity of each model is assessed by comparing the predicted forming limits with experimental results carefully obtained on this steel. The origin of discrepancy in the predicted results using different hardening models is thoroughly analyzed.     

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

Modeling multi-axial deformation of planar anisotropic elasto-plastic materials,part II: Applications

In order to predict the deformations under multi-axial and multi-path loadings in a phenomenological framework, a new rotational-isotropic-kinematic (RIK) hardening model has been suggested in the theory part of the paper combining isotropic, kinematic and rotational hardening. Essential features of this material model are Armstrong–Frederick type backstress components for kinematic hardening and a plastic spin for the rotational hardening describing the evolution of the symmetry axes of the anisotropic yield function.The purpose of this article is to illustrate the significance of the RIK hardening model in sheet metal forming applications as well as in springback predictions. With the rotational hardening and a correction term related to the kinematic hardening, the flow stress in each orientation can be described with few material parameters. Several benchmark problems are considered to illustrate and assess the performance of the RIK hardening model in comparison with other hardening models and experimental results.     

Stress integration method for a nonlinear kinematic/isotropic hardening model and its characterization based on polycrystal plasticity

Sheet metal forming processes generally involve non-proportional strain paths including springback, leading to the Bauschinger effect, transient hardening, and permanent softening behavior, that can be possibly modeled by kinematic hardening laws. In this work, a stress integration procedure based on the backward-Euler method was newly derived for a nonlinear combined isotropic/kinematic hardening model based on the two-yield’s surfaces approach. The backward-Euler method can be combined with general non-quadratic anisotropic yield functions and thus it can predict accurately the behavior of aluminum alloy sheets for sheet metal forming processes. In order to characterize the material coefficients, including the Bauschinger ratio for the kinematic hardening model, one element tension–compression simulations were newly tried based on a polycrystal plasticity approach, which compensates extensive tension and compression experiments. The developed model was applied for a springback prediction of the NUMISHEET’93 2D draw bend benchmark example.     

Nonlocal gradient-dependent modeling of plasticity with anisotropic hardening

This work addresses the formulation of the thermodynamics of nonlocal plasticity using the gradient theory. The formulation is based on the nonlocality energy residual introduced by Eringen and Edelen (1972). Gradients are introduced for those variables associated with isotropic and kinematic hardening. The formulation applies to small strain gradient plasticity and makes use of the evanescent memory model for kinematic hardening. This is accomplished using the kinematic flux evolution as developed by Zbib and Aifantis (1988). Therefore, the present theory is a four nonlocal parameter-based theory that accounts for the influence of large variations in the plastic strain, accumulated plastic strain, accumulated plastic strain gradients, and the micromechanical evolution of the kinematic flux. Using the principle of virtual power and the laws of thermodynamics, thermodynamically-consistent equations are derived for the nonlocal plasticity yield criterion and associated flow rule. The presence of higher-order gradients in the plastic strain is shown to enhance a corresponding history variable which arises from the accumulation of the plastic strain gradients. Furthermore, anisotropy is introduced by plastic strain gradients in the form of kinematic hardening. Plastic strain gradients can be attributed to the net Burgers vector, while gradients in the accumulation of plastic strain are responsible for the introduction of isotropic hardening. The equilibrium between internal Cauchy stress and the microstresses conjugate to the higher-order gradients frames the yield criterion, which is obtained from the principle of virtual power. Microscopic boundary conditions, associated with plastic flow, are introduced to supplement the macroscopic boundary conditions of classical plasticity. The nonlocal formulation developed here preserves the classical assumption of local plasticity, wherein plastic flow direction is governed by the deviatoric Cauchy stress. The theory is applied to the problems of thin films on both soft and hard substrates. Numerical solutions are presented for bi-axial tension and simple shear loading of thin films on substrates.     

Modeling multi-axial deformation of planar anisotropic elasto-plastic materials,part I: Theory

In sheet metal forming processes local material points can experience multi-axial and multi-path loadings. Under such loading conditions, conventional phenomenological material formulations are not capable to predict the deformation behavior within satisfying accuracy. While micro-mechanical models have significantly improved the understanding of the deformation processes under such conditions, these models require large sets of material data to describe the micromechanical evolution and quite enormous computation expenses for industrial applications. To reduce the drawbacks of phenomenological material models under the multi-path loadings a new anisotropic elasto-plastic material formulation is suggested. The model enables the anisotropic yield surface to grow (isotropic hardening), translate (kinematic hardening) and rotate (rotation of the anisotropy axes) with respect to the deformation, while the shape of the yield surface remains essentially unchanged.Essentially, the model is formulated on the basis of an Armstrong–Frederick type kinematic hardening, the plastic spin theory for the reorientation of the symmetry axes of the anisotropic yield function, and additional terms coupling these expressions. The capability of the model is illustrated with multi-path loading simulations in ‘tension-shear’ and ‘reverse-shear’ to assess its performance with ‘cross’ hardening and ‘Bauschinger’ effects.     

A Finite Element Analysis of a Tapered Flat Sheet Tensile Specimen

A uniaxial tension sheet metal coupon with a tapered instead of a straight gage section has been used for centering the location of diffuse neck and for measuring sheet stretchability in a non-uniform strain field. A finite element analysis of such a tensile coupon made of automotive steel sheet metals has been carried out to assess the effect of the tapered gage section geometry and material plastic strain hardening characteristics on the development of local plastic deformation pattern and local stress state, especially beyond the onset of diffuse necking but before localized necking. In particular, the finite element analysis was used in this study to evaluate the accuracy and reliability of an experimental data analysis method for estimating the post-necking effective plastic stress-strain curve based on the direct local surface axial plastic strain measurements for base metal, heat-affected zone, and weld metals of a dual-phase steel DP600. It is concluded that the estimated lower and upper bounds of the effective stress-strain curve at large strains are not satisfactory for low strain-hardening materials such as heat-affected zone and weld metals with the tapered tension coupons. A simple correction method utilizing only the additional local surface strain measurement in the transverse direction is proposed and it is shown to be effective in correcting the estimated effective stress-strain curve of dual-phase steel weld metals obtained for two tapered gage section geometries.     

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.     

准脆性材料的弹塑性各向异性损伤耦合本构模型研究

针对准脆性材料的非线性特征：强度软化和刚度退化、单边效应、侧限强化和拉压软化、不可恢复变形、剪胀及非弹性体胀,在热动力学框架内,建立了准脆性材料的弹塑性与各向异性损伤耦合的本构关系。对准脆性材料的变形机理和损伤诱发的各向异性进行了诠释,并给出了损伤构形和有效构形中各物理量之间的关系。在有效应力空间内,建立了塑性屈服准则、拉压不同的塑性随动强化法则和各向同性强化法则。在损伤构形中,采用应变能释放率,建立了拉压损伤准则、拉压不同的损伤随动强化法则和各向同性强化法则。基于塑性屈服准则和损伤准则,构建了塑性势泛函和损伤势泛函,并由正交性法则,给出了塑性和损伤强化效应内变量的演化规律,同时,联立塑性屈服面和损伤加载面,给出了塑性流动和损伤演化内变量的演化法则。将损伤力学和塑性力学结合起来,建立了应变驱动的应力-应变增量本构关系,给出了本构数值积分的要点。以单轴加载-卸载往复试验识别和校准了本构材料常数,并对单轴单调试验、单轴加载-卸载往复试验、二轴受压、二轴拉压试验和三轴受压试验进行了预测,并与试验结果作了比较,结果表明,所建本构模型对准脆性材料的非线性材料性能有良好的预测能力。     

Directional distortional hardening in metal plasticity within thermodynamics

This paper presents a complete theory for metal plasticity that includes isotropic, kinematic, and directional distortional hardening, within the framework of thermodynamics. Directional distortion is defined here as the formation of a region of high curvature on the yield surface, approximately in the direction of loading, and a region of flattening approximately in the opposite direction, as observed in experiments on various types of metals. The distinguishing features of this theory are the introduction of a fourth order tensor-valued internal variable, whose evolution in conjunction with a directional scalar multiplier describes the evolving directional distortion, and the fact that the hardening laws for all internal variables are derived on the basis of sufficient conditions to satisfy the thermodynamic requirement of positive dissipation. The applicability of the theory is illustrated by fitting experimental data on distorted yield surfaces in the course of plastic deformation.     

Kinematic hardening in large strain plasticity

A finite strain hyper elasto-plastic constitutive model capable to describe non-linear kinematic hardening as well as non-linear isotropic hardening is presented. In addition to the intermediate configuration and in order to model kinematic hardening, an additional configuration is introduced – the center configuration; both configurations are chosen to be isoclinic. The yield condition is formulated in terms of the Mandel stress and a back-stress with a structure similar to the Mandel stress.It is shown that the non-dissipative part of the plastic velocity gradient not governed by the thermodynamical framework and the corresponding quantity associated with the kinematic hardening influence the material behaviour to a large extent when kinematic hardening is present. However, for isotropic elasticity and isotropic hardening plasticity it is shown that the non-dissipative quantities have no influence upon the stress–strain relation.As an example, kinematic hardening von Mises plasticity is considered, which fulfils the plastic incompressibility condition and is independent of the hydrostatic pressure. To evaluate the response and to examine the influence of the non-dissipative quantities, simple shear is considered; no stress oscillations occur.     

On the evolution of isotropic and kinematic hardening with finite plastic deformation Part I: compression/tension loading of OFHC copper cylinders

Compression tests followed by tension tests after re-machining were performed on annealed oxygen-free-high-conductivity copper cylinders. These tests were conducted at nine levels of maximum strain ranging from 5 to 50%. From this data, isotropic and kinematic hardening were calculated using 50, 1000 and 2000 microstrain offset definitions. Both isotropic and kinematic hardening were found to depend on the yield definition. Isotropic hardening, which increased with plastic strain with no signs of saturation, also increased with larger offset definition of yield. Kinematic hardening, which increased to 40% strain and appeared to saturate thereafter, decreased with higher offset definitions of yield.