Anisotropic hardening model based on non-associated flow rule and combined nonlinear kinematic hardening for sheet materials

 A material model for more effective analysis of plastic deformation of sheet materials is presented in this paper. The model is capable of considering the following aspects of plastic deformation behavior of sheet materials: the anisotropy in yielding stresses in different directions by using a quadratic yield function (based on Hill’s 1948 model and stress ratios), the anisotropy in work hardening by introducing non-constant flow stress hardening in different directions, the anisotropy in plastic strains in different directions by using a quadratic plastic potential function and non-associated flow rule (based on Hill’s 1948 model and plastic strain ratios, r-values), and finally some of the cyclic hardening phenomena such as Bauschinger’s effect and transient behavior for reverse loading by using a coupled nonlinear kinematic hardening (so-called Armstrong-Frederick-Chaboche model). Basic fundamentals of the plasticity of the model are presented in a general framework. Then, the model adjustment procedure is derived for the plasticity formulations. Also, a generic numerical stress integration procedure is developed based on backward-Euler method (so-called multistage return mapping algorithm). Different aspects of the model are verified for DP600 steel sheet. Results show that the new model is able to predict the sheet material behavior in both anisotropic hardening and cyclic hardening regimes more accurately. By featuring the above-mentioned facts in the presented constitutive model, it is expected that more accurate results can be obtained by implementing this model in computational simulations of sheet material forming processes. For instance, more precise results of springback prediction of the parts formed from highly anisotropic hardened materials or that of determining the forming limit diagrams is highly expected by using the developed material model.

A study on compressive anisotropy and nonassociated flow plasticity of the AZ31 Magnesium Alloy in hot rolling

Effect of anisotropy in compression is studied on hot rolling of AZ31 magnesium alloy with a three-dimensional constitutive model based on the quadratic Hill48 yield criterion and nonassociated flow rule (non-AFR). The constitutive model is characterized by compressive tests of AZ31 billets since plastic deformations of materials are mostly caused by compression during rolling processes. The characterized plasticity model is implemented into ABAQUS/Explicit as a user-defined material subroutine (VUMAT) based on semi-implicit backward Euler's method. The subroutine is employed to simulate square-bar rolling processes. The simulation results are compared with rolled specimens and those predicted by the von Mises and the Hill48 yield function under AFR. Moreover, strip rolling is also simulated for AZ31 with the Hill48 yield function under non-AFR. The strip rolling simulation demonstrates that the lateral spread generated by the non-AFR model is in good agreement with experimental data. These comparisons between simulation and experiments validate that the proposed Hill48 yield function under non-AFR provides satisfactory description of plastic deformation behavior in hot rolling for AZ31 alloys in case that the anisotropic parameters in the Hill48 yield function and the non-associated flow rule are calibrated by the compressive experimental results.

A non-associated plasticity model with anisotropic and nonlinear kinematic hardening for simulation of sheet metal forming

A material model for more thorough analysis of plastic deformation of sheet materials is presented in this paper. This model considers the following aspects of plastic deformation behavior of sheet materials: (1) the anisotropy in yield stresses and in work hardening by using Hill's 1948 quadratic yield function and non-constant stress ratios which leads to different flow stress hardening in different directions, (2) the anisotropy in plastic strains by using a quadratic plastic potential function and non-associated flow rule, also based on Hill's 1948 model and r-values, and (3) the cyclic hardening phenomena such as the Bauschinger effect, permanent softening and transient behavior for reverse loading by using a coupled nonlinear kinematic hardening model. Plasticity fundamentals of the model were derived in a general framework and the model calibration procedure was presented for the plasticity formulations. Also, a generic numerical stress integration procedure was developed based on backward-Euler method, so-called multi-stage return mapping algorithm. The model was implemented in the framework of the finite element method to evaluate the simulation results of sheet metal forming processes. Different aspects of the model were verified for two sheet metals, namely DP600 steel and AA6022 aluminum alloy. Results show that the new model is able to accurately predict the sheet material behavior for both anisotropic hardening and cyclic hardening conditions. The drawing of channel sections and the subsequent springback were also simulated with this model for different drawbead configurations. Simulation results show that the current non-associated anisotropic hardening model is able to accurately predict the sidewall curl in the drawn channel sections.

Anisotropic behavior in plasticity and ductile fracture of an aluminum alloy

Anisotropic mechanical behavior is investigated for an aluminum alloy of 6K21-IH T4 both in plastic deformation and ductile fracture. Anisotropic plastic deformation is characterized by uniaxial tensile tests of dog-bone specimens, while anisotropy in ductile fracture is illustrated with specimens with a central hole, notched specimens and shear specimens. All these specimens are cut off at every 15º from the rolling direction. The r-values and uniaxial tensile yield stresses are measured from the tensile tests of dog-bone specimens. Then the anisotropic plasticity is modeled by a newly proposed J2-J3 criterion under non-associate flow rule (non-AFR). The testing processes of specimens for ductile fracture analysis are simulated to extract the maximum plastic strain at fracture strokes as well as the evolution of the stress triaxiality and the Lode parameter in different testing directions. The measured fracture behavior is described by a shear-controlled ductile fracture criterion proposed by Lou et al. (2014. Modeling of shear ductile fracture considering a changeable cut-off value for stress triaxiality. Int. J. Plasticity 54, 56-80) for different loading directions. It is demonstrated that the anisotropic plastic deformation is described by the J2-J3 criterion with high accuracy in various loading conditions including shear, uniaxial tension and plane strain tension. Moreover, the anisotropy in ductile fracture is not negligible and cannot be modeled by isotropic ductile fracture criteria. Thus, an anisotropic model must be proposed to accurately illustrate the directionality in ductile fracture.

Study on yield function and plastic potential under non-associated flow for accurate earing prediction in cup drawing

Prediction of more than four ears in a cup drawing process can be successfully achieved by considering r-value and stress directionalities. Yld2004-18p based on associated flow rule and Yld2000-2D based on non-associated flow rule are the examples. The former, however, is more costly in terms of computational efficiency than the latter. In this work, an anisotropic constitutive model based on non-associated flow rule which combines two different functions, Hill (1948) and Yld2000-2d, is implemented to a user defined material model. The accuracy of the anisotropic directionalities (yield stresses and plastic strain ratios) is evaluated. Simulation of a mini-die cup drawing with a body stock alloy predicted eight ears, in good agreement with the experimental results. The use of Hill (1948) model for the yield function and Yld2000-2d for plastic potential under the framework of non-associated flow rule led to accurate prediction of up to eight ears at the lower computational cost.