Next Article in Journal
Intergranular Cracking in Mg-Gd-Y Alloy during Tension Test
Previous Article in Journal
Towards the Development of Novel Diclofenac Multicomponent Pharmaceutical Solids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Yield Function for Anisotropic Aggregate of FCC Crystallites

1
Department of Engineering Mechanics, Nanchang University, 999 Xuefu Avenue, Nanchang 330031, China
2
College of City Construction, Jiangxi Normal University, 99 Ziyang Avenue, Nanchang 330022, China
3
Design and Research Institute, Nanchang University, 235 East Nanjing Road, Nanchang 330096, China
*
Author to whom correspondence should be addressed.
Crystals 2022, 12(8), 1039; https://doi.org/10.3390/cryst12081039
Submission received: 22 June 2022 / Revised: 24 July 2022 / Accepted: 25 July 2022 / Published: 26 July 2022

Abstract

:
In this paper, we give an simple but approximate yield surface for single FCC crystals in Hill’s criterion form by Schmid’s law and nonlinear optimization theory. Assuming that all FCC crystallites in a polycrystal have the same (current) critical resolved shear stress τ c for slip, we derive two closed but approximate yield functions through the orientational averaging of all FCC crystallites’ yield surfaces in the polycrystal. The effect of crystallography on the two yield functions are described by the orientation distribution function.

1. Introduction

The studies of plasticity on polycrystals may be divided into two classes: the mathematical theory of plasticity and physical theory of plasticity [1,2,3,4]. The mathematical approach of plasticity is a formalization of known experimental results and does not penetrate deeply into the physical process of plastic yielding. An alternative approach to the mathematical theory of plasticity is the physical approach of plasticity. A polycrystalline material is an aggregate of tiny crystallites separated by grain boundaries. The microstructure of the polycrystal includes grain orientations and the grain boundary structure. Plastic properties of the polycrystal are related to slip mechanisms of crystallites and the microstructure of the polycrystal.
In order to study the properties of the polycrystal, we take a representative volume element (RVE) from the polycrystal. Although the volume of the RVE is small, the RVE contains information of the crystallographic texture on the polycrystal. The earlier study of polycrystalline plasticity was conducted by Sachs [5] and Taylor [6]. Sachs assumed that all crystallites in a RVE have the same stress. Taylor assumed that all crystallites in a RVE experience the same deformation. Bishop and Hill [7,8] adopted Taylor’s model to calculate a yield function for an isotropic aggregate of cubic crystallites. Let τ c denote a critical resolved shear stress. For an isotropic aggregate of FCC (face centered cubic) crystallites, assuming that all slip systems of cubic crystallites in a polycrystal have the same (current) critical resolved shear stress, the uniaxial tensile yield stresses should be 2.238 τ c under Sachs’ model and 3.06 τ c under Taylor’s model [9]. Maniatty et al. [10] employed the isotropic plasticity equivalent method, in which the actual stress tensor acting on an anisotropic material is transformed into a new stress tensor acting on an “isotropic plasticity equivalent” material, for obtaining an anisotropic yield function. Houtte [11] presented a method of obtaining the yield locus of textured polycrystals by constructing a data bank of the Taylor factors. Both Maniatty et al. and Houtte did not give closed expressions of yield functions with the effect of the texture coefficients.
The micromechanics of polycrystals is interesting only because one can measure the orientations of crystallites by orientation imaging microscopy or by X-ray diffraction. As we know, using the same manufacturing procedure, one can never make two polycrystalline samples with identical crystalline orientations. One can only make samples which seem macroscopically identical. Since the statistical information on the microstructure of a polycrystal can be described mathematically by the orientation distribution function (ODF) and since the plastic deformation of polycrystals involves slips inside crystallites, we will develop a statistical yield function with the effects of the orientation distribution function (ODF) and the the critical resolved shear stress τ c of crystallites.
In this paper, we provide a simple method to obtain two closed but approximate yield functions with the effect of the texture coefficients for the RVE, where the volume average (i.e., orientational averaging) of all crystallites’ yield surfaces is taken as a (macroscopic) yield function of the RVE. Using the physical approach to FCC crystallites, we derive a yield surface of single FCC crystals in Hill’s criterion form by Schmid’s law and nonlinear optimization theory, give a plastic anisotropy tensor M ( R ) of FCC crystallites with orientation R , and obtain the volume average of all crystallites’ yield surfaces in the RVE to obtain two closed but approximate yield functions for the RVE. The two yield functions are based on the assumption that all crystallites in the RVE have the same (current) critical resolved shear stress τ c for slip. The effect of crystallography in the two yield functions is described by the texture coefficients. The first yield function is based on Sachs’ model. The second yield function is based on the assumption that the effect of the ODF on the perturbation stress Δ σ i j is, up to the terms, linear in the texture coefficients. Man [12,13,14] found that for orthorhombic aggregates of cubic crystallites, a macroscopic yield function of the RVE, up to terms, is linear in the texture coefficients.Our two yield functions are different to Man’s yield function [12,13] in three aspects: (1) the two yield functions are for an anisotropic aggregate of cubic crystallites, but Man’s yield function is for an orthorhombic aggregate of cubic crystallites; (2) we can determine the relations between two material parameters and the (current) critical resolved shear stress τ c through some physical assumptions, but Man’s yield function contains two unspecified parameters ( Y 0 and β ) and Man’s yield function is independent of the physical model; (3) Man’s yield function includes the effect of the ODF being, up to terms, linear in the texture coefficients, while our second yield function is quadratic texture dependence.
The critical resolved shear stress τ c is related to the crystalline dislocation mechanism, the crystalline size, and the grain boundary structure. In this paper, however, we do not discuss the hardening process of the critical resolved shear stress τ c during plastic deformation. Herein, we only study the constitutive form of the yield function of the FCC polycrystal for the given ODF and the given τ c .

2. Approximate Yield Surface of FCC Crystal in Hill’s Criterion Form

Assume that a fixed spatial Cartesian coordinate system is chosen. To describe the orientation of a FCC crystallite in the RVE, we pick as a reference a FCC crystallite I , whose three four-fold axes of rotational symmetry coincide with the Cartesian coordinate axes. The orientation of any crystallite in the RVE is then specified by rotation R , which takes the reference crystallite I to the configuration of the crystallite.
By Schmid’s law, we know that yielding for the reference crystallite I would begin on a slip system n t (combination of a slip plane and a slip direction) when the resolved shear stress on this slip plane and in the slip direction reaches a critical value τ c , where the unit normal of the slip plane is denoted by n and the unit vector of the slip direction is denoted by t with t · n = t i n i = 0 . When volume average stress components in the reference crystallite I are denoted by σ i j ( I ) , the shear stress on the slip plane and in the slip direction should be τ = σ i j ( I ) t i n j . Because the shear stress of the reference crystallite on its slip planes and in its slip direction satisfies | τ | τ c , the stress components of the reference crystallite I are in space S :
S = ( σ 11 , σ 22 , σ 33 , σ 23 , σ 31 , σ 12 ) R 6 : | τ | = | σ i j t i n j | τ c ,
where σ i j = σ i j ( I ) . The yield surface S of the reference crystallite I should be
S = ( σ 11 , σ 22 , σ 33 , σ 23 , σ 31 , σ 12 ) R 6 : | τ | = | σ i j t i n j | = τ c .
The space S given in (1) is surrounded by a set of yield functions ( | σ i j t i n j | = τ c ) linear in stress components.
Assume that the space S in (1) can be surrounded approximately by a yield surface
f ( σ i j ( I ) ) = M i j k l ( I ) σ i j ( I ) σ k l ( I ) 1 = 0 ,
where I SO(3) is the second-order identity tensor for the reference crystallite I , M ( I ) is called the plastic anisotropy tensor [4] of the reference crystallite I with the minor and major symmetries, σ i j ( I ) is the volume average stress in the reference crystallite I . Any hydrostatic stress p δ i j does not change the shear stress value of the reference crystallite I because τ = σ i j + p δ i j t i n j = σ i j t i n j by t i n i = 0 for all p R . Hence, we have
f ( σ i j ( I ) ) = f ( σ i j ( I ) + p δ i j ) or 2 p M i j k l δ i j σ k l + p 2 M i j k l δ i j δ k l = 0 , p R
which leads to the traceless condition M i i k l ( I ) = 0 for any pair of indices k and l.
When the reference crystallite I and the traction acting on the reference crystallite I have a rotation R , the stress tensor in the reference crystallite I becomes σ i j R = R i s R j t σ s t ( I ) , and the plastic anisotropy tensor M ( I ) becomes M R with constitutive restriction:
f ( σ i j ( I ) ) = M p q n m ( I ) σ p q ( I ) σ n m ( I ) 1 = M i j k l R σ i j R σ k l R 1 = 0 .
where the three four-fold axes of rotational symmetry of the reference FCC crystallite I coincide with the Cartesian coordinate axes. Considering
M i j k l R σ i j R σ k l R = R i p R j q R k n R l m M p q n m ( I ) R i s R j t σ s t ( I ) R k u R l v σ u v ( I ) = δ s p δ t q δ u n δ v m M p q n m ( I ) σ s t ( I ) σ u v ( I ) = M p q n m ( I ) σ p q ( I ) σ n m ( I )
and R i p R i s = δ s p , we have M i j k l R = R i p R j q R k n R l m M p q n m ( I ) in (5). Assuming that all crystallites in the RVE have the same (current) critical resolved shear stress τ c for slip, we have the plastic anisotropy tensor M ( R ) of the crystallites with orientation R :
M i j k l ( R ) = M i j k l R = R i p R j q R k n R l m M p q n m ( I ) .
If we use σ i j ( R ) to denote the volume average stress of the crystallites with orientation R , the yield surface f ( σ i j ( R ) ) for the crystallites with orientation R can be expressed as
f ( σ i j ( R ) ) = M i j k l ( R ) σ i j ( R ) σ k l ( R ) 1 = 0
by (5) and (7). Since the reference FCC crystallite I has the octahedral symmetry O, the plastic anisotropy tensor should satisfy M i j k l R = M i j k l ( I ) for each rotation tensor R O , which with the minor and major symmetries of M ( I ) reads [15]
M I J ( I ) = m 11 m 12 m 12 0 0 0 m 12 m 11 m 12 0 0 0 m 12 m 12 m 11 0 0 0 0 0 0 m 44 0 0 0 0 0 0 m 44 0 0 0 0 0 0 m 44 ,
where M I J ( I ) is the matrix form of M ( I ) in Voigt’s notation (i.e., indices i j and k l in M i j k l are denoted by indices I and J of M I J with convention 11 1 , 22 2 , 33 3 , 23 4 , 13 5 , 12 6 ; for instances, M 1122 = M 12 , M 1123 = M 14 ), m 11 , m 12 , and m 44 are the material parameters of the reference crystallite I . The plastic anisotropy tensor M I J ( I ) in (9) can be decomposed into [16]
M I J ( I ) = m 12 B I J ( 1 ) + 2 m 44 B I J ( 2 ) + ( m 11 m 12 2 m 44 ) B I J ( 3 ) ( I )
where B I J ( 1 ) , B I J ( 2 ) , and B I J ( 3 ) ( I ) are the matrix form of the following tensors in Voigt’s notation, respectively
B i j k l ( 1 ) = δ i j δ k l , B i j k l ( 2 ) = 1 2 ( δ i k δ j l + δ i l δ j k ) , B i j k l ( 3 ) ( I ) = α = 1 3 δ α i δ α j δ α k δ α l ,
B ( 1 ) and B ( 2 ) are isotropic fourth-order tensor bases, B ( 3 ) ( I ) is a fourth-order tensor basis for cubic crystals. The traceless condition M i i k l ( I ) = 0 makes the plastic anisotropy tensor M ( I ) in (10) become
M i j k l ( I ) = m 12 B i j k l ( 1 ) + 2 m 44 B i j k l ( 2 ) + ( 3 m 12 2 m 44 ) B i j k l ( 3 ) ( I )
because of m 11 = 2 m 12 which makes that
M i i k l ( I ) = m 12 3 δ k l + 2 m 44 δ k l + ( m 11 m 12 2 m 44 ) δ k l 0
holds.
Putting (12) into (3), we find a yield surface of the reference crystallite I in Hill’s criterion form
f ( σ ( I ) ) = m 12 σ i i σ k k + 2 m 44 σ i j σ i j + ( 3 m 12 2 m 44 ) σ 11 2 + σ 22 2 + σ 33 2 = 1 η 2 ( σ 22 σ 33 ) 2 + ( σ 33 σ 11 ) 2 + ( σ 11 σ 22 ) 2 + 1 ς 2 σ 23 2 + σ 31 2 + σ 12 2 1 = 0
because
B i j k l ( 1 ) σ i j σ k l = σ i i σ k k , B i j k l ( 2 ) σ i j σ k l = σ i j σ i j , B i j k l ( 3 ) ( I ) σ i j σ k l = σ 11 2 + σ 22 2 + σ 33 2 ,
where σ i j = σ i j ( I ) and
η 2 = 1 m 12 , ς 2 = 1 4 m 44 .
When 3 m 12 2 m 44 = 0 , M ( I ) in (12) is an isotropic tensor whose corresponding yield surface is von Mises’ yield criterion.
For the reference FCC crystallite I , the slip occurs on {111} planes in <110> directions. The unit normal n of the slip planes and the unit vector t of the slip directions constitute the 24 slip systems n t [9]
{ n t } = { ± 1 3 ( 111 ) 1 2 [ 1 1 ¯ 0 ] , , ± 1 3 ( 11 1 ¯ ) 1 2 [ 011 ] } ,
in which there is t · n = 0 for each slip system. Substituting (17) into | τ | = | σ i j t i n j | τ c in (1), one can obtain 24 inequalities that are linear in stress components
| τ | = 1 6 | ( σ 22 σ 33 ) + σ 31 ± σ 12 | τ c , | τ | = 1 6 | ( σ 22 σ 33 ) σ 31 ± σ 12 | τ c , | τ | = 1 6 | ( σ 33 σ 11 ) + σ 12 ± σ 23 | τ c , | τ | = 1 6 | ( σ 33 σ 11 ) σ 12 ± σ 23 | τ c , | τ | = 1 6 | ( σ 11 σ 22 ) + σ 31 ± σ 23 | τ c , | τ | = 1 6 | ( σ 11 σ 22 ) σ 31 ± σ 23 | τ c
which constitute the the stress space S .
Now we assume that the stress space S in (1) can approximately be surrounded by a function in Hill’s criterion form (14). In order to determine η and ς in (14), we introduce two subspaces of S in (1)
S 1 = ( σ 11 , σ 22 , σ 33 , 0 , 0 , 0 ) R 6 : | τ | = | σ i j t i n j | τ c , S 2 = ( 0 , 0 , 0 , σ 23 , σ 31 , σ 12 ) R 6 : | τ | = | σ i j t i n j | τ c
and two subspaces in R 6
M ( η ) = { ( σ 11 , σ 22 , σ 33 , 0 , 0 , 0 ) R 6 : ( σ 22 σ 33 ) 2 + ( σ 33 σ 11 ) 2 + ( σ 11 σ 22 ) 2 η 2 0 } , N ( ς ) = ( 0 , 0 , 0 , σ 23 , σ 31 , σ 12 ) R 6 : σ 23 2 + σ 31 2 + σ 12 2 ς 2 0 .
The maximum and minimum values of η and ς in (20) are defined by
η max = inf { η : M ( η ) S 1 } , ς max = inf { η : N ( ς ) S 2 } , η min = sup { η : M ( η ) S 1 } , ς min = sup { ς : N ( ς ) S 2 } ,
where we determine η and ς in (14) through multi-axial tensions and pure shears, respectively. The proceeding maximum and minimum problems in (21) belong to nonlinear optimization problems. The relations (19) and (20) are used to obtain the best fitting of the FCC crystal yield surface | σ i j t i n j | τ c .
We use Microsoft Excel to solve the problems (21) with (20) and (18). Then we obtain the maximum and minimum values of η and ς
η min = 3 τ c , η max = 12 τ c , ς min = 3 τ c , ς max = 6 τ c .
Putting η = η max and ς = ς max into (14), we find a yield surface of FCC crystals for the reference FCC crystallite I
1 12 τ c 2 ( σ 22 σ 33 ) 2 + ( σ 33 σ 11 ) 2 + ( σ 11 σ 22 ) 2 + 1 6 τ c 2 σ 23 2 + σ 31 2 + σ 12 2 1 = 0
which is just the same as that of Maniatty et al. [10].
Since the goal of the paper is to estimate a macroscopic yield function of the RVE, let us choose η and ς in (14) as follows:
η = 1 2 ( η min + η max ) = 3.232 τ c , ς = 1 2 ( ς min + ς max ) = 2.091 τ c .
If all crystallites in the RVE have the same (current) critical resolved shear stress τ c for slip, we combine (7) with (12), (16) and (24) to obtain the plastic anisotropy tensor
M i j k l ( R ) = 1 η 2 B i j k l ( 1 ) + 1 2 ς 2 B i j k l ( 2 ) + ω B i j k l ( 3 ) ( R )
and the yield surface f ( σ i j ( R ) ) for FCC crystallites with orientation R
f ( σ i j ( R ) ) = 1 η 2 σ i i ( R ) 2 + 1 2 ς 2 σ i j ( R ) σ i j ( R ) + ω B i j k l ( 3 ) ( R ) σ i j ( R ) σ k l ( R ) 1 = 0
where ω = 3 η 2 1 2 ς 2 and
B i j k l ( α ) ( R ) = R i p R j q R k n R l m B p q n m ( α ) = B i j k l ( α ) , α = 1 , 2 , B i j k l ( 3 ) ( R ) = R i p R j q R k n R l m B p q n m ( 3 ) ( I ) = α = 1 3 R i α R j α R k α R l α .

3. Yield Function and Effective Plastic Anisotropy Tensor of RVE

The crystallographic texture of a polycrystal is described by the orientation distribution function (ODF). The ODF is used to describe the probability density of finding a crystallite with orientation R in the RVE. The mechanical anisotropy of the RVE is caused mainly by its crystallographic texture and plastic anisotropy of crystallites. Let w ( R ) be the orientation distribution function. If L 2 ( SO ( 3 ) ) is the space of the square-integrable complex-valued functions, w L 2 ( SO ( 3 ) ) can be expanded as an infinite series in terms of the Wigner D -functions [17,18]
w ( R ) = w iso + l = 1 m = l l n = l l c m n l D m n l R , c m n l = ( 1 ) m n ( c m ¯ n ¯ l ) * ,
where w iso = 1 8 π 2 , c m n l ( l 1 ) are the texture coefficients, ( c m ¯ n ¯ l ) * denote the complex conjugate of the complex number c m ¯ n ¯ l , n ¯ = n . The Wigner D-functions D m n l constitute an orthogonal basis in L 2 ( SO ( 3 ) ) . The texture coefficients that we shall need can easily be measured by X -ray diffraction. The texture coefficients c m n l are related to Roe’s W l m n coefficients by formula [19]
W l m n = ( 1 ) n m 2 2 l + 1 c m n l .
The micromechanics of polycrystals becomes interesting only because one can obtain statistical information of the orientations by X-ray diffraction or measure the orientations of pcrystallites, pointwise, by orientation imaging microscopy. For the measurement of X-ray diffraction, one measures the specimen surface’s, 1 4 depth’s, and 1 2 depth’s texture coefficients, respectively, by sanding the specimen of sheets. The mean values of the texture coefficients on the surface, 1 4 depth, and 1 2 depth of the specimen are taken as the texture coefficients of the specimen. For the measurement of the orientation imaging microscopy, one obtains the orientation (Euler angles ψ p , θ p , ϕ p ) of each crystallite Ω p in the specimen. Since the FCC polycrystal can be taken as an aggregate of cubic crystallites (i.e., Ω = p = 1 N Ω p ), the texture coefficients of the specimen are given by [20]
c m n l = 1 | Ω | p = 1 N | Ω p | c m n 4 ( R p )
with
c m n l ( R p ) = 2 l + 1 8 π 2 · 1 24 i = 1 24 D m n l ( R p Q cr ( i ) ) * , Q cr ( i ) O ,
where R p = R ( ψ p , θ p , ϕ p ) and O is the group of cubic crystal symmetry.
Let σ ¯ i j be the volume average stress components of the RVE. By definition of the ODF, the orientational (volume) average stress can be obtained by
σ ¯ i j = SO ( 3 ) σ i j ( R ) w ( R ) d g ,
where g = 8 π 2 g H , g H is the Haar measure on SO ( 3 ) with g H ( SO ( 3 ) ) = 1 , σ i j ( R ) is the volume average stress tensor of crystallites with orientation R in the RVE. When R is defined by Roe’s notation [18], there is d g = sin θ d ψ d θ d ϕ .
When all crystallites in the RVE have the same (current) critical resolved shear stress τ c for slip, the yield surface f ( σ i j ( R ) ) of crystallites with orientation R is determined by (26). By the ODF definition, the volume average f ( σ i j ( R ) ) ¯ of all crystallites in the RVE can be expressed as
f ( σ i j ( R ) ) ¯ = SO ( 3 ) f σ i j ( R ) w ( R ) d g = 0 2 π 0 π 0 2 π f σ i j ( R ) w ( R ) sin θ d ψ d θ d ϕ = 0
which can be taken as the macroscopic yield function of the RVE.
If the macroscopic yield function f ( σ i j ( R ) ) ¯ can be expressed as
f ( σ i j ( R ) ) ¯ = σ ¯ i j M i j k l eff σ ¯ k l 1 = 0
in Hill’s criterion form, then we call M eff in (34) the effective plastic anisotropy tensor of the RVE. The effective plastic anisotropy tensor M eff gives the relationship between the macroscopic yield function and the volume average stress of the RVE. One recipe to compute the effective plastic anisotropy tensor M eff in (34) is by way of putting (8) into (33)
σ ¯ i j M i j k l eff σ ¯ k l = SO ( 3 ) σ i j ( R ) M i j k l ( R ) σ k l ( R ) w ( R ) d g .

4. Volume Average of Plastic Anisotropy Tensor

The basic assumption for Sachs’ model [5] is that all crystallites in the RVE experience the same state of stress (i.e., σ i j ( R ) = σ ¯ i j , R SO(3)). Substituting σ i j ( R ) = σ ¯ i j into (35) reads σ ¯ i j M i j k l eff σ ¯ k l = SO ( 3 ) σ ¯ i j M i j k l ( R ) σ ¯ k l w ( R ) d g . Since the equation holds for each σ ¯ i j , the effective plastic anisotropy tensor has to be the volume average plastic anisotropy tensor M ¯ of the RVE
M ¯ i j k l = SO ( 3 ) M i j k l ( R ) w ( R ) d g = SO ( 3 ) 1 η 2 B i j k l ( 1 ) + 1 2 ς 2 B i j k l ( 2 ) + ω B i j k l ( 3 ) ( R ) w ( R ) d g
by (25).
From the results given by Huang and Man [21] and Huang [16,22], we know
SO ( 3 ) B i j k l ( 1 ) w ( R ) d g = B i j k l ( 1 ) , SO ( 3 ) B i j k l ( 2 ) w ( R ) d g = B i j k l ( 2 )
because of B i j k l ( 1 ) and B i j k l ( 2 ) being isotropic fourth order tensors and
SO ( 3 ) B i j k l ( 3 ) ( R ) w ( R ) d g = α = 1 3 0 2 π 0 π 0 2 π R i α R j α R k α R l α 1 8 π 2 sin θ d ψ d θ d ϕ + s = 1 m = s s n = s s c m n s α = 1 3 0 2 π 0 π 0 2 π R i α R j α R k α R l α D m n s R sin θ d ψ d θ d ϕ = 1 5 B i j k l ( 1 ) + 2 5 B i j k l ( 2 ) + Φ i j k l ,
because of the relation [19]
SO ( 3 ) B i j k l ( 3 ) ( R ) D m n s R d g = 0 , when s > 4
where B i j k l ( 3 ) ( R ) is a fourth-order tensor. For an anisotropic aggregate of cubic crystallites, Φ in (38) is totally symmetric (i.e., Φ i 1 i 2 i 3 i 4 = Φ i τ ( 1 ) i τ ( 2 ) i τ ( 3 ) i τ ( 4 ) for any permutation τ of { 1 , 2 , 3 , 4 } ) and traceless with non-trivial components [21]:
Φ 2233 = a 1 , Φ 1133 = a 2 , Φ 1122 = a 3 , Φ 1123 = a 5 a 8 , Φ 1113 = a 7 + 3 a 4 , Φ 1112 = a 6 + a 9 , Φ 3323 = 4 a 5 , Φ 3313 = 4 a 4 , Φ 3312 = 2 a 6
and
a 1 = 32 π 2 105 ( c 00 4 + 5 2 Re c 20 4 ) , a 2 = 32 π 2 105 ( c 00 4 5 2 Re c 20 4 ) , a 3 = 8 π 2 105 ( c 00 4 70 Re c 40 4 ) , a 4 = 8 5 π 2 105 Re c 10 4 , a 5 = 8 5 π 2 105 Im c 10 4 , a 6 = 8 10 π 2 105 Im c 20 4 , a 7 = 8 35 π 2 105 Re c 30 4 , a 8 = 8 35 π 2 105 Im c 30 4 , a 9 = 8 70 π 2 105 Im c 40 4 .
Here Re c m 0 4 and Im c m 0 4 denote the real and imaginary parts of the complex number c m 0 4 , respectively. c 00 4 is a real number because of (28)2.
Putting (37) and (38) into (36), we obtain an effective plastic anisotropy tensor M eff of the RVE under the Sachs’ model
M i j k l eff = M ¯ i j k l = ρ B i j k l ( 1 ) + 3 ρ B i j k l ( 2 ) + ω Φ i j k l
with
ρ = 1 5 ( 2 η 2 + 1 2 ς 2 ) = 1 16.348 τ c 2 , ω = 3 η 2 1 2 ς 2 = 1 5.786 τ c 2 .
The texture coefficients c m 0 4 above can be measured by the X -ray diffraction and inversion of pole figures (Roe, 1965 [18]). Obviously, the effective plastic anisotropy tensor satisfies the traceless condition M i i k l eff = 0 because of Φ i i k l = 0 . Substituting (42) and (15)1,2 into (34), we obtain a macroscopic yield function of the RVE
f ( σ i j ( R ) ) ¯ = ρ σ ¯ i i σ ¯ j j + 3 ρ σ ¯ i j σ ¯ i j + ω Φ i j k l σ ¯ i j σ ¯ k l = 1 Y 0 2 3 2 σ ¯ i j σ ¯ i j 1 2 σ ¯ i i σ ¯ j j + β Φ i j k l σ ¯ i j σ ¯ k l 1 = 0
where
Y 0 = 1 2 ρ = 2.859 τ c , β = ω 2 ρ = 1.413 .
The macroscopic yield function f ( σ i j ( R ) ) ¯ of the RVE with a group of texture symmetry G tex can be obtained from the constraints imposed by G tex on the texture coefficients. For instance, if the coordinate axes are the two-fold axes of orthorhombic symmetry of the texture, then c m 0 4 = 0 for odd m and Im c m 0 4 = 0 for all m. Substituting these restriction into (44) and (41), we obtain the macroscopic yield function of the RVE in the Hill’s criterion form
f ( σ i j ( R ) ) ¯ = F ( σ ¯ 22 σ ¯ 33 ) 2 + G ( σ ¯ 33 σ ¯ 11 ) 2 + H ( σ ¯ 11 σ ¯ 22 ) 2 + 2 L σ ¯ 23 2 + 2 M σ ¯ 31 2 + 2 N σ ¯ 12 2 1 = 0
for an orthorhombic aggregate of FCC crystallites, where
( F , G , H ) = 1 Y 0 2 1 2 ( 1 , 1 , 1 ) β ( a 1 , a 2 , a 3 ) , ( L , M , N ) = 1 Y 0 2 3 2 ( 1 , 1 , 1 ) + 2 β ( a 1 , a 2 , a 3 ) .
a 1 , a 2 , and a 3 in (41) are determined by the texture coefficients. Y 0 and β in (45) are dependent on the (current) critical resolved shear stress τ c . The formula (46) is the same as that of Man [12,13] in form. However, Man’s yield function is independent of the physical model; hence Y 0 and β for Man’s yield function are two unspecified material parameters.

5. Yield Function with Quadratic Texture Dependence

Under Sachs’ model, the effective plastic anisotropy tensor M eff is the volume average plastic anisotropy tensor of the RVE. Man’s yield function [12,13] and the yield function (46) contain the effect of the ODF only being, up to terms, linear in the texture coefficients, which may not suffice for strongly textured polycrystals [20]. In this section, we will establish a new physical model for deriving a macroscopic yield function of the RVE. The macroscopic yield function delineates the effect of the crystallographic texture on the plastic anisotropy up to terms quadratic in the texture coefficients.
Considering the collection of crystallites with orientation R in the RVE and using σ i j ( R ) to denote the volume average stress in these crystallites, we have the average perturbation stress of these crystallites
Δ σ i j = σ i j ( R ) σ ¯ i j ,
where σ ¯ i j is given in (32). We assume that the average perturbation stress Δ σ i j of crystallites with orientation R are governed by a constitutive relation of the form
Δ σ i j = Δ σ i j ( R , w , σ ¯ k l )
with constraint
Δ σ i j ( R , w , 0 ) = 0 .
From (32) and (48), we have
SO ( 3 ) Δ σ i j ( R , w , σ ¯ k l ) w ( R ) d g = 0 .
Similar to the process of deriving the HV-V model given by Huang and Man [21], the average perturbation stress Δ σ i j of the crystallites with orientation R can approximately be expressed as
Δ σ i j = κ Ψ i j k l ( R ) Φ i j k l σ ¯ k l ,
where we consider the effect of the ODF to Δ σ i j being, up to terms, linear in the texture coefficients,
Ψ i j k l ( R ) = 1 5 B i j k l ( 1 ) 2 5 B i j k l ( 2 ) + B i j k l ( 3 ) ( R ) ,
and κ is an unspecified material parameter. From (37) and (38), we know that Δ σ i j in (52) satisfy the relation (51).
Through (35) and (48), we obtain an effective plastic anisotropy tensor by relation
σ ¯ i j M i j k l eff σ ¯ k l = SO ( 3 ) σ i j ( R ) M i j k l ( R ) σ k l ( R ) w ( R ) d g = SO ( 3 ) σ ¯ i j M i j k l ( R ) σ ¯ k l w ( R ) d g + SO ( 3 ) Δ σ i j M i j k l ( R ) Δ σ k l + Δ σ i j M i j k l ( R ) σ ¯ k l + σ ¯ i j M i j k l ( R ) Δ σ k l w ( R ) d g .
From (52) and (25), we have
σ ¯ i j M i j k l ( R ) Δ σ k l = κ 2 ς 2 σ ¯ i j Ψ i j k l ( R ) Φ i j k l σ ¯ k l + κ ω σ ¯ i j 1 5 B i j k l ( 1 ) + 3 5 B i j k l ( 3 ) ( R ) B i j m n ( 3 ) ( R ) Φ m n k l σ ¯ k l .
by relations B i j m n ( 1 ) Ψ m n k l ( R ) = B i j m n ( 1 ) Φ m n k l = 0 and
B i j m n ( 1 ) B m n k l ( 1 ) = 3 B i j k l ( 1 ) , B i j m n ( 1 ) B m n k l ( 3 ) ( R ) = B i j k l ( 1 ) , B i j m n ( 3 ) ( R ) B m n k l ( 3 ) ( R ) = B i j k l ( 3 ) ( R ) , B i j m n ( 2 ) B m n k l ( α ) ( R ) = B i j k l ( α ) ( R ) , α = 1 , 2 , 3 .
By means of the relations (37) and (38), we can complete the following integrations and obtain the volume average value on (55)
SO ( 3 ) Δ σ i j M i j k l ( R ) σ ¯ k l w ( R ) d g = SO ( 3 ) σ ¯ i j M i j k l ( R ) Δ σ k l w ( R ) d g = κ ω σ ¯ i j 2 25 B i j k l ( 1 ) + 6 25 B i j k l ( 2 ) + 1 5 Φ i j k l Φ i j m n Φ m n k l σ ¯ k l .
Similarly, from (52) and (25) we have
Δ σ i j M i j k l ( R ) Δ σ k l = κ 2 σ ¯ i j ( Ψ i j p q ( R ) Φ i j p q ) [ 1 η 2 B p q m n ( 1 ) + 1 2 ς 2 B p q m n ( 2 ) + ω B p q m n ( 3 ) ( R ) ] ( Ψ m n k l ( R ) Φ m n k l ) σ ¯ k l
whose volume average value is
SO ( 3 ) Δ σ i j M i j k l ( R ) Δ σ k l w ( R ) d g = κ 2 2 ς 2 σ ¯ i j 2 25 B i j k l ( 1 ) + 6 25 B i j k l ( 2 ) + 1 5 Φ i j k l Φ i j m n Φ m n k l σ ¯ k l + κ 2 ω σ ¯ i j 6 125 B i j k l ( 1 ) + 18 125 B i j k l ( 2 ) 3 25 Φ i j k l 4 5 Φ i j m n Φ m n k l + o ( | Φ 2 | σ ¯ k l .
Substituting (36), (42), (57), and (58) into (54), we obtain the effective plastic anisotropy tensor for anisotropic aggregates of FCC crystallites with quadratic texture dependence
M i j k l eff = ρ eff B i j k l ( 1 ) + 3 ρ eff B i j k l ( 2 ) + ω eff Φ i j k l + ϑ eff Φ i j m n Φ m n k l
with
ρ eff = ρ + 4 25 κ ω + 6 25 κ 2 ρ + 2 125 κ 2 ω = ( 0.0612 + 0.0277 κ + 0.0175 κ 2 ) / τ c 2 , ω eff = ω + 2 5 κ ω + 3 5 κ 2 ρ 1 5 κ 2 ω = 0.00213 ( κ + 29.573 ) ( κ + 2.731 ) / τ c 2 , ϑ eff = 2 κ ω 3 κ 2 ρ 2 5 κ 2 ω = 0.253 ( κ + 1.368 ) κ / τ c 2
where 1 ς 2 = 2 5 ( 15 ρ 2 ω ) because of (43). Substituting (59) into (34), we obtain a macroscopic yield function of the RVE
f ( σ i j ( R ) ) ¯ = 1 Y 0 eff 2 3 2 σ ¯ i j σ ¯ i j 1 2 σ ¯ i i σ ¯ j j + β eff Θ i j k l σ ¯ i j σ ¯ k l + γ eff Θ i j m n Θ m n k l σ ¯ i j σ ¯ k l 1 = 0
where we discard terms o ( | Φ 2 | ) in (59) and
Y 0 eff = 2 ρ eff 1 / 2 = 0.707 τ c ( 0.0612 + 0.0277 κ + 0.0175 κ 2 ) 1 / 2 , β eff = ω eff 2 ρ eff = 0.00107 ( κ + 29.573 ) ( κ + 2.731 ) 0.0612 + 0.0277 κ + 0.0175 κ 2 , γ eff = ϑ eff 2 ρ eff = 0.126 ( κ + 1.368 ) κ 0.0612 + 0.0276 κ + 0.0174 κ 2 .
The effective plastic anisotropy tensor in (59) satisfies the traceless condition M i i k l eff = 0 .
In metallurgical practice, because of the effects of alloying elements, it is difficult to obtain accurate estimates of the material parameters τ c and κ in (61). The macroscopic yield function (61) is then more properly looked upon as a representation formula with three undetermined parameters Y 0 eff ,   β eff and γ eff .

6. Discussion and Example

When the yield function is in Hill’s criterion form, the Taylor factor M ( 0.5 , 0 ) for an orthorhombic weakly textured aggregate of FCC crystallites can be expressed as [14]
M ( 0.5 , 0 ) = σ ¯ 11 τ c = 3.067 10.095 W 400 2 10 3 W 420 + 70 3 W 440 ,
where σ ¯ 11 is a uniaxial tensile yield stress.
Under Sachs’ model, we obtained the macroscopic yield function (46) for an aggregate of FCC crystallites. For uniaxial tensile problems, the yield function in (46) becomes G + H ( σ ¯ 11 ) 2 1 = 0 or
σ ¯ 11 = ( G + H ) 1 / 2 = Y 0 ( 1 β ( a 2 + a 3 ) ) 1 / 2 Y 0 1 + β 2 ( a 2 + a 3 )
from (47). When the RVE is an orthorhombic aggregate of FCC crystallites, there are Re c 20 4 = c 20 4 and Re c 40 4 = c 40 4 (Roe, 1965 [18]) in (41). From (41), (45) and (29), we rewrite (64) into
σ ¯ 11 = 2.859 τ c 9.666 τ c W 400 2 10 3 W 420 + 70 3 W 440
which is close to the result (63). However, the derivation of (65) is much simpler than that of (63) because the computation of the Taylor factor for a textured polycrystal is very complicated and time consuming.
The parameter κ in (62) can be determined by an experiment. For instance, from the relation of an isotropic polycrystal (aluminum) stress–strain curve and a single crystal stress–strain curve, Taylor [6] found that Y 0 eff = 3.06 τ c agrees with the experimental data very well. Taking Y 0 eff = 3.06 τ c , we obtain κ = 0.365 by solving (62)1 and then get β eff = 1.384 and γ eff = 0.866 from (62)2,3. For an orthorhombic weakly textured aggregate of FCC crystallites, if we discard the o ( c m 0 4 ) terms of the texture coefficients in (61), similar to the procedure of deriving (46) and (64) from (44), we obtain the yield stress of uniaxial tension
σ ¯ 11 = Y 0 eff 1 + β eff 2 ( a 2 + a 3 ) = 3.06 τ c 10.133 τ c W 400 2 10 3 W 420 + 70 3 W 440 .
The formula above is almost the same as (63). To check the form of the expression (66), we give one example as follows.
Example of yield stress under uniaxial tension of sample S θ :
Take an example S θ from a metal sheet with the rolling direction θ as shown in Figure 1. To obtain the expression of the uniaxial yield stress σ θ and the rolling direction θ , we give the stress tensor of a uniaxial stress σ θ on as follows
σ ( θ ) = σ θ n n = σ θ cos 2 θ cos θ sin θ 0 cos θ sin θ sin 2 θ 0 0 0 0
where n = [ cos θ , sin θ , 0 ] T . Similar to the derivation of (66), for an orthorhombic weakly textured aggregate of FCC crystallites, putting σ ij ( θ ) in (67) and W 4 m 0 = 2 9 c m 0 4 in (29) into (61) and after discarding the o ( W m 0 4 ) terms of the texture coefficients, we obtain the yield stress σ θ of the uniaxial tension sample S θ as follows:
σ θ = 3 . 06 τ c 10 . 133 τ c W 400 2 10 3 W 420 cos 2 θ + 70 3 W 440 cos 4 θ .
Yu Xiang [23] in the doctoral dissertation gave the experimental data on two aluminum alloys (AA5754 hot band, AA3105 O-temper). As shown in Figure 1, uniaxial test samples were cut with tensile axis with the rolling direction (RD) θ = 0 , 15 , 30 , , 90 . For each tensile direction, measurements were made on two samples to obtain the flow stress at 4% longitudinal plastic strain. The experimental results are listed in Table 1, where the data Y i ( 1 ) and Y i ( 2 ) of Xiang’s Experiment 1 and Experiment 2 are labeled with superscripts (1) and (2), respectively, and the subscript 1 i 7 indicates a sample with tensile axis making an angle ϕ i = ( i 1 ) × 15 with the rolling direction (RD) shown in Figure 1. The texture coefficients of the two alloys were measured by X-ray diffraction as shown in Table 2.
To determine τ c in (68) for the aluminum alloys from Xiang’s experimental data, we introduce an objective function as follows:
Π ( τ c ) = i = 1 7 [ 1 2 ( Y i ( 1 ) + Y i ( 2 ) ) σ θ | θ = θ i ] 2 , θ i = ( i 1 ) × 15 .
By the method of least squares, we solve the equation Π ( τ c ) τ c = 0 and obtain the critical resolved shear stress τ c of the FCC crystals for the two aluminum alloys
τ c = 86.93 MPa for AA 5754 hot band ; τ c = 41.66 MPa for AA 3105 O - temper .
Putting the fitted values τ c in (70) and the texture coefficients W lmn in Table 2 for the two alloys into (68), we plot the fitting results of the yield stress σ θ on the uniaxial tension sample S θ as shown in Figure 2. The expression (68) of the yield stress σ θ on the uniaxial tension sample S θ can fit Xiang’s experimental data on two aluminum alloys (AA5754 and AA3105) very well.

7. Conclusions

In this paper, we provide a simple method to obtain two closed but approximate yield functions with the effect of the texture coefficients for the RVE, where the volume average (i.e., orientational averaging) of all crystallites’ yield surfaces is taken as a (macroscopic) yield function of the RVE. Through theoretical analysis, we draw the following important conclusions:
  • Assuming that all crystallites in the RVE have the same (current) critical resolved shear stress τ c for slip and employing Schmid’s law and the nonlinear optimization theory, we give a plastic anisotropy tensor M ( R ) of the FCC crystals with orientation R and obtain the volume average of all crystallites’ yield surfaces in the RVE.
  • By the Microsoft Excel, we solve nonlinear optimization problems and obtain an approximate plastic anisotropy tensor and a yield surface of FCC crystallites with orientation.
  • For an orthorhombic weakly textured aggregates of FCC crystallites, we obtain two yield stresses of uniaxial tension, respectively.
  • Through an example, we obtain the yield stress under the uniaxial tension on sample S θ . We find that the expression of the yield stress can fit Xiang’s experimental data very well.

Author Contributions

Conceptualization, L.Z.; software, T.Z.; formal analysis, L.Z. and G.M.; investigation, G.M.; resources, Z.L.; data curation, M.H.; writing—original draft preparation, L.Z. and M.H.; writing—review and editing, T.Z. and Z.L.; project administration, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Jiangxi Graduate Education and Teaching Reform Research Project (Awards Nos. JXYJG-2021-057) and the National Natural Science Foundation of China (Awards Nos. 51568046).

Data Availability Statement

The data of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hill, R. The Mathematical Theory of Plasticity; Clarendon Press: Oxford, UK, 1950. [Google Scholar]
  2. Lin, T.H. Physical theory of plasticity. Adv. Appl. Mech. 1971, 11, 255–311. [Google Scholar]
  3. Khan, A.S.; Huang, S. Continuum Theory of Plasticity; John Wiley &Sons: New York, NY, USA, 1995. [Google Scholar]
  4. Gambin, W. Plasticity and Textures; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2001. [Google Scholar]
  5. Sachs, G. Zur ableilung einer fleissbedingung. Z. Verein Deut. Ing. 1928, 72, 734–736. [Google Scholar]
  6. Taylor, G.I. Plastic strain in metals. J. Inst. Met. 1938, 62, 307–324. [Google Scholar]
  7. Bishop, J.F.W.; Hill, R. XLVI. A theory of the plastic distortion of a polycrystalline aggregate under combined stresses. London Edinb. Dublin Philos. Mag. J. Sci. 1951, 42, 414–427. [Google Scholar] [CrossRef]
  8. Bishop, J.F.W.; Hill, R. CXXVIII. A theoretical derivation of the plastic properties of a polycrystalline face-centred metal. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1951, 42, 1298–1307. [Google Scholar] [CrossRef]
  9. Hosford, W.F. The Mechanics of Crystals and Textured Polycrystals; Oxford University: New York, NY, USA, 1993. [Google Scholar]
  10. Maniatty, A.M.; Yu, J.S.; Keane, T. Anisotropic yield criterion for polycrystalline metals using texture and crystal symmetries. Int. J. Solids Struct. 1999, 36, 5331–5355. [Google Scholar] [CrossRef]
  11. Houtte, P.V. Calculation of the yield locus of textured polycrystals using the Taylor and the relaxed Taylor theory. Textures Microstruct. 1970, 7, 29–72. [Google Scholar] [CrossRef]
  12. Man, C.S. Elastic compliance and Hill’s quadratic yield function for weakly orthotropic sheets of cubic metals. Metall. Mater. Trans. A 1994, 25, 2835–2837. [Google Scholar] [CrossRef]
  13. Man, C.S. On the correlation of elastic and plastic anisotropy in sheet metals. J. Elast. 1995, 39, 165–173. [Google Scholar] [CrossRef]
  14. Man, C.S.; Huang, M. Identification of material parameters in yield functions and flow rules for weakly textured sheets of cubic metals. Int. J. Non-Linear Mech. 2001, 36, 501–514. [Google Scholar] [CrossRef]
  15. Mura, T. Micromechanics of Defects in Solids; Martinus Nijhoff Publishers: Dordrecht, The Netherlands, 1987. [Google Scholar]
  16. Huang, M. Perturbation approach to elastic constitutive relations of polycrystals. J. Mech. Phys. Solids 2004, 52, 1827–1853. [Google Scholar] [CrossRef]
  17. Bunge, H.J. Texture Analysis in Material Science: Mathematical Methods; Butterworths: London, UK, 1982. [Google Scholar]
  18. Roe, R.J. Description of crystallite orientation in polycrystalline materials: III, General solution to pole figure inversion. J. Appl. Phys. 1965, 36, 2024–2031. [Google Scholar] [CrossRef]
  19. Man, C.S. On the Constitutive Equations of Some Weakly-Textured Materials. Arch. Ration. Mech. Anal. 1998, 143, 77–103. [Google Scholar] [CrossRef]
  20. Huang, M.; Man, C.S. A finite-element study on constitutive relation HM-V for elastic polycrystals. Comput. Mater. Sci. 2005, 32, 378–386. [Google Scholar] [CrossRef]
  21. Huang, M.; Man, C.S. Constitutive relation of elastic polycrystal with quadratic texture dependence. J. Elast. 2003, 72, 183–212. [Google Scholar] [CrossRef]
  22. Huang, M. Elastic constants of a polycrystal with an orthorhombic texture. Mech. Mater. 2004, 36, 623–632. [Google Scholar] [CrossRef]
  23. Xiang, Y. Effects of Grain Shape and Crystallographic Texture on Plastic Anisotropy of Aluminum Alloy Sheets. Ph.D. Thesis, University of Kentucky, Lexington, KY, USA, 2004. [Google Scholar]
Figure 1. Samples S θ of metal sheet.
Figure 1. Samples S θ of metal sheet.
Crystals 12 01039 g001
Figure 2. Experimental data of Xiang and fitting results of (68).
Figure 2. Experimental data of Xiang and fitting results of (68).
Crystals 12 01039 g002
Table 1. Flow stresses Y i ( j ) (MPa) of two aluminum alloys.
Table 1. Flow stresses Y i ( j ) (MPa) of two aluminum alloys.
Aluminum Alloy15°30°45°60°75°90°
Hot Band Y i ( 1 ) 282.31279.14274.78258.58259.41270.97279.58
AA5754 Y i ( 2 ) 282.00278.82270.15257.99257.53271.77280.87
O-temper Y i ( 1 ) 126.72125.13124.11122.15124.19124.72126.10
AA3105 Y i ( 2 ) 125.19125.43124.82124.50123.91124.65126.04
Table 2. Texture coefficients W lmn of two aluminum alloys.
Table 2. Texture coefficients W lmn of two aluminum alloys.
AA5754 (Hot band) W 400 = 0.0033 W 420 = 0.0025 W 440 = 0.0071
AA3105 (O-temper) W 400 = 0 . 0072 W 420 = 0 . 0019 W 440 = 0 . 0022
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, L.; Min, G.; Zhao, T.; Huang, M.; Lan, Z. Estimation of Yield Function for Anisotropic Aggregate of FCC Crystallites. Crystals 2022, 12, 1039. https://doi.org/10.3390/cryst12081039

AMA Style

Zhou L, Min G, Zhao T, Huang M, Lan Z. Estimation of Yield Function for Anisotropic Aggregate of FCC Crystallites. Crystals. 2022; 12(8):1039. https://doi.org/10.3390/cryst12081039

Chicago/Turabian Style

Zhou, Lu, Guanhua Min, Tengfei Zhao, Mojia Huang, and Zhiwen Lan. 2022. "Estimation of Yield Function for Anisotropic Aggregate of FCC Crystallites" Crystals 12, no. 8: 1039. https://doi.org/10.3390/cryst12081039

APA Style

Zhou, L., Min, G., Zhao, T., Huang, M., & Lan, Z. (2022). Estimation of Yield Function for Anisotropic Aggregate of FCC Crystallites. Crystals, 12(8), 1039. https://doi.org/10.3390/cryst12081039

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop