Effect of Single Crystal Orientation on Forming

Online since 05 April 2021

Article

Abstract

Among processes involving plastic deformation, sheet metal forming requires a most accurate description of plastic anisotropy. One of the main sources of mechanical anisotropy is the intrinsic anisotropy of the constituent crystals. In this paper, we present the single-crystal yield criterion recently developed by Cazacu et al. [1] and its application to the prediction of anisotropy in uniaxial tension of strongly textured polycrystalline sheets. Namely, it is shown that using this single crystal yield criterion the Lankford coefficients exist and have finite values for all loading orientations. Moreover, the variation of both the yield stress and Lankford coefficients with the crystallographic direction can be expressed analytically. An application of this criterion to forming a cylindrical cup from a single crystal of (100) orientation is presented. Finally, we show that using this single-crystal model, one can describe well the effect of the spread around an ideal texture component on the anisotropy in uniaxial tensile properties of a polycrystal.

Keywords

Single Crystal, Yield Criterion, Texture, Deep Drawing, Finite Element Analysis

Table of contents

Text

1 Introduction

Since the publication of Lankford landmark paper [2] in which it was shown that the material’s initial anisotropy may have a beneficial influence on forming performance, extensive experimental and theoretical efforts have been undertaken towards development of appropriate anisotropic yield criteria. For polycrystalline metallic materials, versatile orthotropic yield criteria have been developed (for example see [3], [4]). In comparison very few yield criteria have been proposed in the literature for the description of the yielding anisotropy of single crystals (e.g. see [5]). Recently, in [1] were demonstrated the benefits of using a single crystal yield criterion that is defined for any 3-D loadings. Specifically, it was shown that this criterion involves the correct number of independent anisotropy parameters it captures the differences in yield stress anisotropy between different crystals, e.g. the differences in the anisotropy in the tensile uniaxial yield stresses with the orientation between the loading direction and the crystallographic directions (e.g. see the examples provided in [1] for aluminum and copper single crystals). Moreover, it is possible to obtain analytically the expression of the Lankford coefficient and tensile uniaxial yield stress along any crystallographic orientation (see [6]). Another advantage of using this criterion is the fact that it is C2 differentiable. This in turn facilitates its implementation into finite element (FE) codes. Furthermore, using for the description of the individual grain behavior this single-crystal model, one can describe well the anisotropy of polycrystalline materials (e.g. see [7]).

In this paper, we present an application of this single crystal yield criterion to forming a cylindrical cup from a single-crystal sheet of (100) orientation. Moreover, we show that using this single-crystal model, one can describe well the effect of the spread of orientations around (100) on the anisotropy in uniaxial tensile properties.

2 Single-crystal yield criterion

The single-crystal yield criterion is defined for any stress-state. It is written in terms of cubic invariants that were deduced using rigorous theorems of representation of tensor functions. Therefore, the properties of invariance of the yield function with respect to the intrinsic symmetries associated to crystals belonging to the cubic system are automatically satisfied. Moreover, by using these representation theorems, it is ensured that the criterion involves the correct number of anisotropy coefficients such as to satisfy the crystal symmetries and the condition of yielding insensitivity to the hydrostatic pressure (for full mathematical proofs and further details, see [6] and Cazacu et al. [1]).

Let us define as Oxyz the Cartesian coordinate system associated with the crystal axes (i.e., the <100> crystal directions). In this coordinate system, the expression of the effective stress associated with the single-crystal yield criterion [1] is given by:

Image 100000000000056D000000B91B0120495DF1541E.png

In Eq. (1), σ´ denotes the Cauchy stress deviator, m1, m2, n1, n3, n4 are anisotropy coefficients while the parameter c describes the relative importance of the second-order and third-order cubic stress-invariants on yielding of the crystal. For example, all the parameters of the criterion can be determined using the experimental tensile uniaxial yield stresses in four crystallographic directions.

3 Application to forming a cup from a single crystal sheet

In the following we present the application of the single crustal yield criterion (1) to the prediction of the earing profile of a cup drawn from a single crystal sheet having the rolling (RD), transverse (TD), and normal direction (ND) of the sheet along the <100> axes of the crystal.

A schematic view of the cup-drawing set-up is shown in Fig. 1. The tool and blank dimensions are given in Table 1. As a result of material and geometrical symmetries of the problem, only one-quarter of the cup needs to be considered in the F.E. analysis. The blank is meshed with 4200 Abaqus C3D8H elements (8 node brick elements with constant pressure, see ABAQUS [8]). The circular blank has an initial thickness of 0.813 mm and an initial diameter of 79 mm. A blank force of 1000 N is used.

Image 1000000000000542000001B94A5A4668A97C9ACC.png

Fig. 1. (a) Schematic of the deep drawing experiment; (b) F.E. mesh for the blank.

The F.E. analysis was performed in ABAQUS/Standard using a user material subroutine (UMAT) developed by the authors for an elastic/plastic model with yielding described by the Cazacu et al. [1] single-crystal yield criterion (see Eq. (1)), associated flow rule, and a Swift-type hardening law, i.e.

Image 100000000000055A0000006A2E2B5FD5EFC7BC69.png

where K0, ε0 and n are parameters, and ε̅p is the equivalent plastic strain associated with the effective stress given by Eq. (1) using the work-equivalence principle. The numerical values considered for these parameters are K0= 478 MPa, ε0=0.0004 and n= 0.17. An implicit-time integration algorithm was used for the implementation of the model. In the elastic regime, the crystal was assumed isotropic, with Young’s modulus E = 69 GPa and Poisson’s ratio ν=0.3.

Table 1. Geometry of the tools for the deep drawing experiment

Image 1000000000000548000000A92700C96AAC31A0A3.png

The numerical values for the anisotropy coefficients considered are: m1=1, m2 =0.38, n1=0.98, n3=0.04, n4=0.08, c=2.3. The deformed shape of a fully drawn cup, namely the isocontours of the equivalent plastic strain are shown in Fig. 2(a). Note that the peaks in cup height are at 0° and 90° orientations (see Fig. 2(b) for the F.E. evolution of the cup height with the angle from RD ).

Image 100000000000053C000001BA05EFC5BD64A437C7.png

Fig. 2. (a) Deformed shape of a fully drawn cup of a {100}<001> single crystal sheet according to Cazacu et al. [1] single-crystal yield criterion showing the isocontours of the equivalent plastic strain; (b) Predicted cup height vs. the angle from the [100] direction (RD).

It is well known that even in the case of strongly textured polycrystalline sheets for which it can be considered that the texture has a single component, more a spread is generally observed around this component. Therefore, it is of great interest to estimate the effect of this spread on the plastic anisotropy in the uniaxial properties at the polycrystal level.

In the next section, using the single crystal yield criterion [1] for describing the plastic behavior of the constituent crystals, we describe the anisotropy in uniaxial tensile properties for sheets containing only the (100) component and being characterized by various scatter widths from this component (see Fig.3).

The effective stress of the polycrsytal, σ̅poly(σ) corresponding to the applied stress tensor, σ, is expressed in the loading frame is considered to be of the form:

Image 10000000000005A10000005035C228E523FAED4B.png

with N being the number of crystals considered in the polycrystal, σ̅grainj denotes the effective stress of the crystal j calculated using Eq.(1), Rj is the transformation matrix for passing from the frame associated to the crystal axes of the grain j to the loading frame. The plastic strain rate tensor of the polycrystal, dp, is expressed in the loading frame as:

Image 1000000000000539000000749FA691887FE1E46A.png

4 Influence of the spread from the ideal (100) cube texture on the evolution of the uniaxial tensile yield stress and Lankford coefficients

Simulations were carried out for samples of 400 crystals using the polycrystal model given by Eq. (3)-(4) for polycrystalline sheets with textures of scatter widths of 15º and 30º about the ideal (100) cube texture component (see Figure 3(b)-(c) for the (111) pole figures of the respective textures).

Image 10000000000005410000018152DF2F201E11B49D.png

Fig. 3. (111) pole figure for sheets with: (a) ideal cube texture (ω0=0°) and respectively a spread of (b) ω0=15° (c) ω0=30° about the ideal {100}<001> cube texture.

The predicted evolution of the anisotropy in uniaxial yield stresses and Lankford coefficients in the (RD, TD) plane are shown in Fig. 4.

Since in the case of a sheet with an ideal cube texture all constituent crystals have the same orientation, a single crystal yield criterion can be used. It is worth noting that with the Cazacu et al [1] yield criterion, it is possible to obtain analytically the variation of the tensile uniaxial yield stresses and Lankford coefficients with the loading direction. Indeed, for uniaxial tension in the (RD, TD) plane, with respect to the coordinate system associated with the <100> crystal axes, the only non-zero components of the stress tensor are: σxx=σ(θ)cos2θ, σxy=σ(θ)sinθcosθ, and σyy(θ)sin2θ, with σ(θ) denoting the yield stress in a direction at an angle θ with respect to the RD direction, [100]. Further substitution in Eq. (1) leads to:

Image 1000000000000559000000A93F924ED897D193FA.png

It is worth noting that the symmetry of the ideal cube-textured sheet is correctly described by the Cazacu et al. [1] single crystal model. Indeed, Eq.(6) is invariant by the transformation θ (90°-θ), so the yield stress remains identical under rotations of π/2 about the normal to the sheet. In particular, irrespective of the values of the parameters m1, m2, n1, n3, n4 , and c, we have σ()=σ(9) i.e. the yield stresses in uniaxial tension along the [100] and [010] directions are the same (see Fig. 4(a)). By the flow rule, the Lankford coefficient r(θ) is given by

Image 100000000000054E000000A7E4F669E9BC7FD927.png

with the effective stress σ̅ given by Eq.(1). It can be shown that irrespective of the values of the parameters of the criterion, r(θ)=r(90°-θ), and in particular, r(0°)=r(90°)=1 (for full proof, see [5]). In Fig. 4, the evolution of Lankford coefficients r(θ) and yield stress ratios σ(θ)/σ(0) with the loading direction θ for the sheets with textures given in Fig. 3, were calculated with the same set of values for the parameters describing the crystal level yielding behavior , namely m1=1.0, m2 =0.38, n1=0.98, n3=0.04, n4=0.08, c=2.3.

Image 100000000000054400000222E1F6F69B77F2D076.png

Fig. 4. Numerical simulation using the polycrystal model of the anisotropy in (a) yield stress ratio σ(θ)/σ(0) and (b) Lankford coefficients r(θ) in the plane of the cube textured polycrystalline sheets. The textures for different scatter width are shown in Fig.3. The results for the ideal texture ω0=0 were obtained with the analytical formulas (Eq.(5)-Eq.(6)), respectively.

Note that even for the ideal cube texture both the predicted yield stresses and Lankford coefficients vary smoothly with the loading orientation (Fig.4), for the values of the parameters considered, there are extrema only at 0º , 45º, and 90º to the [100] direction in the (100) plane. It is worth recalling that for an ideal cube texture, the yield stress variation with the loading orientation according to the Taylor-Bishop-Hill model displays two cusps and infinite values are predicted for the Lankford coefficients along the 0º (i.e. RD) and 90º (i.e. TD) tensile loadings (see [9]). Additionally, the TBH model predicts a smooth variation only if the texture is characterized by a very large spread (see simulation results presented in [9]).

5 Conclusions

The capabilities of the single-crystal yield criterion recently developed by Cazacu et al [1] was illustrated. It was shown that using this single crystal yield criterion, the Lankford coefficients exist and have finite values for all loading orientations. Moreover, the variation of both the yield stress and Lankford coefficients with the crystallographic direction can be expressed analytically. Furthermore, it was shown that using Cazacu et al.[1] single-crystal model to describe the anisotropy of individual constituent crystals, the effect of texture on the Lankford coefficients and yield stresses of a polycrystalline sheet can be reasonably described. With regards to deep drawing simulation, according to Cazacu et al [1] single-crystal model, the number and location of the extrema in the height of the drawn cup is dictated by the crystal orientation, for e.g. a minimum of four ears for the (100) crystal orientation. Whether at the respective locations, the criterion predicts minima or peaks for the height of the cup depends on the parametrization used for the criterion. For the parameterization considered in this papers, peaks are at 0º and 90º. The same general conclusion concerning the influence of the parameterization on the earing predictions also holds true for all the yield criteria for polycrystalline materials (e.g. Hill 1948 [10], Yld 91 [11], Yld 2004-18p [12], BBC 2010 [13], Cazacu et al. [14], etc.) as has been clearly demonstrated over the years in the literature and Numisheet benchmarks.

Acknowledgements

Partial support for this work provided by the AFOSR grant # FA9550-18-1-0517 is gratefully acknowledged.

Bibliography

[1] Cazacu, O. Revil-Baudard, B. Chandola, N. A yield criterion for cubic single crystals. International Journal of Solids Structures, 2018, 151, 9-19.

[2] Lankford, W. T. Snyder, S. C. Bausch, J. New Criteria for Predicting the Press Performance of Deep Drawing Steels. Transactions of the American Society of Metals, 1950, 42, 1197-1232.

[3] Cazacu, O. Revil-Baudard, B. Plasticity of Metallic Materials: Modeling and Applications to Forming. Elsevier, ISBN 978-0-12-817984, 535 pp, 2021.

[4] Banabic, D. Barlat, B. Cazacu, O. Kuwabara, T. Advances in anisotropy of plastic behaviour and formability of sheet metals. International Journal of Material Forming, 2020, 13(5), 749-787.

[5] Cazacu, O., Revil-Baudard, B., Chandola, N. Plasticity-Damage Couplings: From Single Crystal to Polycrystalline Materials. Springer, ISBN 978-3-319-92921-7 2019.

[6] Cazacu, O. Chandola, N. Revil-Baudard, B. Analytical expressions for the yield stress and Lankford coefficients of polycrystalline sheets based on a new single crystal model. International Journal of Material Forming, 2018, 11(4), 571-581.

[7] Cazacu, O. Chandola, N. Revil-Baudard, B. Frodal, B.H. Børvik, T. Hopperstad, O.S. Modeling the effect of notch geometry on the deformation of a strongly anisotropic aluminum alloy. European Journal of Mechanics-A/Solids, 2020, 104004.

[8] ABAQUS, Abaqus Documentation. Dassault systemes. Providence, RI, USA, 2011.

[9] Lequeu, P. Gilormini, P. Montheillet, F. Bacroix, B. Jonas, J. Yield surfaces for textured polycrystals—I. Crystallographic approach. Acta Metallurgica, 1987, 35:439–451

[10] Hill, R. A theory of the yielding and plastic flow of anisotropic metals. In: Proceedings of the Royal Society of London A: mathematical, physical and engineering sciences, The Royal Society, 1948, 281–297.

[11] Barlat, F. Panchanadeeswaran, S., Richmond, O. Prediction of earing in cup drawing fcc materials. Texture, Stress, and Microstructure, 1991, 14, 507-512.

[12] Barlat F, Aretz H, Yoon J, Karabin M, Brem J, Dick R. Linear transfomation-based anisotropic yield functions. International Journal of Plasticity, 2005, 21:1009–1039.

[13] Banabic, D. Barlat, F. Cazacu, O. Kuwabara, T. Advances in anisotropy and formability. International Journal of Material Forming, 2010, 3, 165-189.

[14] Cazacu, O. Plunkett, B. Barlat, F. Orthotropic yield criterion for hexagonal closed packed metals. International Journal of Plasticity, 2006, 22:1171–1194.

Illustrations

Table 1. Geometry of the tools for the deep drawing experiment

Table 1. Geometry of the tools for the deep drawing experiment

References

Electronic reference

Nitin Chandola, Oana Cazacu and Benoit Revil-Baudard, « Effect of Single Crystal Orientation on Forming  », ESAFORM 2021 [Online], Online since 05 April 2021, connection on 28 March 2024. URL : https://popups.uliege.be/esaform21/index.php?id=1621

Authors

Nitin Chandola

Department of Mechanical and Aerospace Engineering, University of Florida/REEF, Shalimar, FL 32579-1163, USA

Oana Cazacu

Department of Mechanical and Aerospace Engineering, University of Florida/REEF, Shalimar, FL 32579-1163, USA

Corresponding author: Oana Cazacu. E-mail address: cazacu@reef.ufl.edu

Benoit Revil-Baudard

Department of Mechanical and Aerospace Engineering, University of Florida/REEF, Shalimar, FL 32579-1163, USA

Details

Title
Effect of Single Crystal Orientation on Forming
Language
en
Author, co-author
Nitin Chandola, Oana Cazacu and Benoit Revil-Baudard,
Publication date
14 April 2021
Journal title
ESAFORM 2021
Copyright
CC-BY
DOI
10.25518/esaform21.1621
Permanent URL
https://popups.uliege.be/esaform21/index.php?id=1621

Statistics

Views : 288 (0 ULiège)
Downloads : 0 (0 ULiège)

Cite

Chandola, N., Cazacu, O., & Revil-Baudard, B. (2021). Effect of Single Crystal Orientation on Forming . Paper presented at ESAFORM 2021. 24th International Conference on Material Forming, Liège, Belgique. doi: 10.25518/esaform21.1621