Numerical implementation of orthotropic plasticity for sheet-metal forming analysis

Abstract Due to rolling and crystallographic texture development, thin sheet-metal for stamping exhibits an initial orthotropic plastic behaviour. The objective of the present work is to define a stress-computation algorithm taking into account this anisotropic plastic behaviour and the update of the orthotropic frame directions within a finite-element simulation of sheet-metal forming. Following Mandel's approach, a macroscopic behaviour model is considered which is based on the use of a rotational objective derivative defined from the initial micro-structure rotation governing the orthotropic direction update. The expression of the constitutive law in a frame rotated by this rotation permits a numerical scheme including an elastic prediction and a plastic correction for the stress calculation. The elastic prediction ensures the incremental objectivity and a second-order level of accuracy. The correction stage is based on an implicit Euler-backward scheme and leads to the solution of a non-linear equation using a local Newton method. In the proposed simulation approach, C 0 shell finite elements are used to model the blank. The proposed stress-calculation algorithm takes into account the mixed interpolation used in the element formulation in order to avoid shear locking. A set of numerical results is presented in cases of the simulation of deep-drawing processes. The accuracy of the solutions, and especially the good agreement with experimental results, shows the efficiency of the proposed approach.

[1]  R. Hill The mathematical theory of plasticity , 1950 .

[2]  O. C. Zienkiewicz,et al.  Analysis of thick and thin shell structures by curved finite elements , 1970 .

[3]  P. Boisse,et al.  A simple isoparametric three-node shell finite element , 1992 .

[4]  P. Boisse,et al.  A C0 three‐node shell element for non‐linear structural analysis , 1994 .

[5]  Siavouche Nemat-Nasser,et al.  Theoretical foundation for large-scale computations for nonlinear material behavior , 1984 .

[6]  R. Borst,et al.  Studies in anisotropic plasticity with reference to the Hill criterion , 1990 .

[7]  J. Mandel,et al.  Equations constitutives et directeurs dans les milieux plastiques et viscoplastiques , 1973 .

[8]  T. Hughes,et al.  Finite rotation effects in numerical integration of rate constitutive equations arising in large‐deformation analysis , 1980 .

[9]  B. Loret On the effects of plastic rotation in the finite deformation of anisotropic elastoplastic materials , 1983 .

[10]  J. Jonas,et al.  Axial stresses and texture development during the torsion testing of Al, Cu and α-Fe , 1984 .

[11]  M. Rashid,et al.  The effect of plastic spin on anisotropic material behavior , 1989 .

[12]  P. Pinsky,et al.  Numerical integration of rate constitutive equations in finite deformation analysis , 1983 .

[13]  F. Barlat,et al.  A six-component yield function for anisotropic materials , 1991 .

[14]  J. Jonas,et al.  Relation between axial stresses and texture development during torsion testing: A simplified theory , 1985 .