Next Article in Journal
Cosmic-Ray Extremely Distributed Observatory
Next Article in Special Issue
Finite Pure Plane Strain Bending of Inhomogeneous Anisotropic Sheets
Previous Article in Journal
Discrepancy between Power Radiated and the Power Loss Due to Radiation Reaction for an Accelerated Charge
Previous Article in Special Issue
A Limit Load Solution for Anisotropic Welded Cracked Plates in Pure Bending
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of the Replacement of the Actual Plastic Orthotropy with Various Approximations of Normal Anisotropy on Residual Stresses and Strains in a Thin Disk Subjected to External Pressure

1
Metal Forming Department, Samara National Research University, 34 Moskovskoye shosse, 334086 Samara, Russia
2
Samara Federal Research Center of the Russian Academy of Sciences, 3A Studencheskiy Pereulok, 443001 Samara, Russia
3
School of Mechanical Engineering and Automation, Beihang University, 37 Xueyuan Road, Beijing 100191, China
4
Ishlinsky Institute for Problems in Mechanics, 101-1 Prospect Vernadskogo, Moscow 119526, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(11), 1834; https://doi.org/10.3390/sym12111834
Submission received: 1 October 2020 / Revised: 21 October 2020 / Accepted: 28 October 2020 / Published: 5 November 2020

Abstract

:
Plastic anisotropy is very common to metallic materials. This property may significantly affect the performance of structures. However, the actual orthotropic yield criterion is often replaced with a criterion based on the assumption of normal anisotropy. The present paper aims to reveal the influence of this replacement on the distribution of strains and residual strains in a thin hollow disk under plane stress conditions. The boundary-value problem is intentionally formulated such that it is possible to obtain an exact semi-analytical solution without relaxing the boundary conditions. It is assumed that the disk is loaded by external pressure, followed by elastic unloading. The comparative analysis of the distributions of residual strains shows a significant deviation of the distribution resulting from the solutions based on the assumption of normal anisotropy from the distribution found using the actual orthotropic yield criterion. This finding shows that replacing the actual orthotropic yield criterion with the assumption of normal anisotropy may result in very inaccurate predictions. The type of anisotropy accepted is of practical importance because it usually results from such processes as drawing end extrusion with an axis of symmetry.

1. Introduction

Many metallic materials are plastically anisotropic. The orthotropic form of anisotropy is most common, arising from such metal forming processes as rolling, drawing and extrusion [1]. This material property affects both subsequent sheet metal forming processes and the performance of structures and machine parts made of the products of these processes [2]. For example, the residual stresses and, hence, the elastic springing are very sensitive to plastic anisotropy [3,4]. The effect of plastic anisotropy on the solution behavior for thin rotating disks has been revealed in [5,6,7]. However, for simplifying theoretical calculations, the real orthotropic yield criterion is often replaced with a transversely isotropic yield criterion (i.e., it is assumed that the material properties are independent of the direction within the transverse plane). This type of anisotropy is called normal anisotropy. The effect of the replacement above is usually not evaluated. The primary objective of the present paper is to evaluate this effect for a particular structure.
The assumption of normal anisotropy is often used for the analysis of various sheet metal forming processes. In [8,9,10,11,12], theoretical analyses for the elastic/plastic bending of sheet metal have been carried out to predict springback, bendability, strain and stress distributions, and the maximum load on the punch. In [13], a numerical and analytical study on the stretch flanging of V-shaped sheet metal has been presented. A non-linear correlation between the hole expansion ratio and tensile properties, including the coefficient of normal anisotropy, has been revealed in [14,15]. Many instability (or failure) criteria used for predicting the limit strain and formability of sheet metal are applicable only for normal anisotropy [16,17,18]. In [19,20,21], the influence of plastic anisotropy on the plastic behavior and failure of porous ductile materials has been studied using different modifications of the Gurson yield criterion. Still, in all these cases, the matrix is supposed to be transversely isotropic. The assumption of normal anisotropy is also used to study the effect of anisotropy on the hardness of billets after upsetting [22], and to estimate the no-rivet connection quality [23].
Moreover, the assumption of normal anisotropy is often used not only to analyze technological processes but also for the estimation of plastic anisotropy of different materials. For example, in [24], Al/Cu bimetallic sheets have been studied, and the copper layer thickness’s influence on the normal anisotropy has been reported. The normal anisotropy coefficients have also been used. The authors of [25] studied the influence of continuous annealing temperature on battery shell steel’s mechanical properties. In [26], the normal anisotropy of crystallographic orientations has been calculated to evaluate the effect of crystallographic texture on the stretch-flangeability during the hole expansion test of hot-rolled multiphase steel.
There are various expressions for calculating the parameters involved in transversely isotropic yield criteria [2,27,28,29,30]. These expressions depend on the number of directions along which the yield stress is available. Therefore, it is of interest to reveal the effect of different transversely isotropic yield criteria on the distribution of residual stresses and strains in structures. For this purpose, it is not necessary to solve a complicated boundary value problem. It is more important to derive an analytic or semi-analytic solution without relaxing any boundary conditions. In [31], the effect of replacing the exact yield criterion by transversely isotropic yield criteria on the distribution of residual stresses in a thin hollow disk subjected to external pressure and subsequent unloading has been investigated. In the present study, this analysis is extended to the distribution of residual strains. A comprehensive overview of the boundary value problem solved has been provided in [32]. This overview does not include solutions for anisotropic materials.

2. Statement of the Problem

In many cases, the actual orthotropic yield criterion (i.e., this criterion describes experimental data with sufficient accuracy) is replaced with the assumption of normal anisotropy [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23]. The effect of this replacement on the final result is unknown. It may depend on the actual yield criterion and its approximations. In the present paper, the quadratic yield criterion proposed by Hill [33] is chosen as the actual yield criterion. However, any other criterion can play this role. Having Hill’s yield criterion, one can find parameters that are necessary for formulating any yield criterion based on the assumption of normal anisotropy. These different yield criteria are used for finding the distributions of residual stresses and strains under the same system of loading. The difference in these distributions is considered as the effect of the replacement of the actual yield criterion with approximate yield criteria.
It is seen from the above that the choice of the boundary value problem is not so crucial for the primary objective of this study. For the reason of transparency of the final result, it is desirable to have semi-analytic solutions without relaxing the equations and boundary conditions. The boundary value problem below satisfies this requirement.
Consider the deformation of a thin hollow disk of inner radius a0 and outer radius b0 by a uniform external pressure p 0 , followed by unloading (Figure 1). The internal pressure is zero. It is natural to solve this boundary value problem in a cylindrical coordinate system ( r , θ , z ) whose z axis coincides with the axis of symmetry of the disk. Let σ r , σ θ and σ z be the normal stresses referred to this coordinate system. The state of stress is plane, σ z = 0 . The stress boundary conditions are
σ r = 0
for r = a 0 and
σ r = p 0
for r = b 0 .
The constitutive equations comprise Hooke’s law, the yield criterion, and its associated flow rule. It is assumed that the material of the disk is orthotropic and obeys the Hill’s quadratic yield criterion [33]. Moreover, the principal axes of anisotropy coincide with the coordinate curves of the cylindrical coordinate system. Then, since σ z = 0 , the yield criterion is represented as
( G + H ) σ r 2 2 H σ r σ θ + ( H + F ) σ θ 2 = 1 ,
where
2 F = 1 Y 2 + 1 Z 2 1 X 2 , 2 G = 1 Z 2 + 1 X 2 1 Y 2 , 2 H = 1 X 2 + 1 Y 2 1 Z 2
and X, Y and Z are the yield stresses in tension in the radial, circumferential, and axial directions, respectively. One can also express the coefficients F, G, and H in terms of the Lankford’s coefficients R 0 and R 90 as [2]
F = 1 X 2 R 0 R 90 ( R 0 + 1 ) , G = 1 X 2 1 ( R 0 + 1 ) , H = 1 X 2 R 0 ( R 0 + 1 ) .
here, the coefficient R 0 corresponds to the radial direction and the coefficient R 90 to the circumferential direction. Actual orthotropic yield criteria are often approximated by the normally anisotropic yield criterion of the form
σ r 2 2 R ¯ ( 1 + R ¯ ) σ r σ θ + σ θ 2 = σ ¯ s 2 ,
where σ ¯ s is the average yield stress in tension in any direction in planes z = constant. The value of R ¯ is expressed through the Lankford’s coefficients in different directions. Several expressions are available in the literature [2,27,28,29,30]. All of them can be represented as [30]
R ¯ = { R 0 + 2 i = 1 n 2 R Δ α i + R 90 2 ( n 1 ) if n > 2 , R 0 + R 90 2 if n = 2 .
Similar expressions are used for the average yield stress in tension:
σ ¯ s = { X + 2 i = 1 n 2 σ s Δ α i + Y 2 ( n 1 ) if n > 2 , X + Y 2 if n = 2 .
The difference between the expressions in [2,27,28,29,30] is in the value of n. In Equations (7) and (8), n is the number of directions in which the Lankford’s coefficients and yield stresses are measured (including the radial and circumferential directions), R Δ α i is the Lankford’s coefficient in the direction γ = Δ α i , σ s Δ α i is the yield stress in tension in the direction γ = Δ α i , γ is the angle measured from the radial direction anticlockwise, and Δ α = π / [ 2 ( n 1 ) ] .
It is evident from Equations (1)–(3) and (6) that the solution for any yield criterion is independent of θ . Moreover, the stresses σ r , σ θ and σ z are the principal stresses. Therefore, the only non-trivial equilibrium equation is
σ r r + σ r σ θ r = 0 .
Let ε r , ε θ and ε z be the normal strains referred to the cylindrical coordinate system. These strains are the principal strains. Therefore, the Hooke’s law reduces to
E ε r e = σ r ν σ θ , E ε θ e = σ θ ν σ r , E ε z e = ν ( σ r + σ θ ) .
here, the superscript e denotes the elastic portion of the strain components, E is the Young’s modulus, and ν is the Poisson’s ratio. The plastic flow rule associated with Yield Criterion (3) is
ε ˙ r p = λ [ ( G + H ) σ r H σ θ ] , ε ˙ θ p = λ [ ( F + H ) σ θ H σ r ] , ε ˙ z p = λ ( G σ r + F σ θ )
and with Yield Criterion (6) is
ε ˙ r p = λ 1 [ σ r R ¯ σ θ ( 1 + R ¯ ) ] , ε ˙ θ p = λ 1 [ σ θ R ¯ σ r ( 1 + R ¯ ) ] , ε ˙ z p = λ 1 ( 1 + R ¯ ) ( σ r + σ θ ) .
here, λ and λ 1 are non-negative multipliers, and the superscript p denotes the plastic portion of the strain components, and the superimposed dot denotes the time derivative. The total strain components are determined as
ε r = ε r e + ε r p , ε θ = ε θ e + ε θ p , ε z = ε z e + ε z p .
In what follows, it is convenient to use the dimensionless quantities:
ρ = r b 0 , a = a 0 b 0 , k = X E , s = σ ¯ s X , p = p 0 X .
The present paper aims to show the effect of replacing Yield Criterion (3) with (6) on the distribution of stresses and strain at the end of loading and residual stresses and strains.

3. Stress Solution

A stress solution has been found in [31]. For the reason of readability, the main results of this solution are presented in this section.
In the case of Yield Criterion (3), the distribution of the radial and circumferential stresses in the plastic region a 0 r r c (or a ρ ρ c ) is given in parametric form as
σ r X = 2 sin ψ 4 η 2 , σ θ X = η 1 ( η sin ψ 4 η 2 + cos ψ ) , ln ( ρ a ) = η 1 4 η 2 2 [ 1 + η 1 ( η 1 η ) ] ψ + + ( 2 η η 1 ) 2 [ 1 + η 1 ( η 1 η ) ] ln [ η 1 4 η 2 η 1 4 η 2 cos ψ ( 2 η η 1 ) sin ψ ] ,
where η = 2 H X Y , η 1 = Y / X , ψ is an auxiliary variable, r c is the radius of the elastic/plastic boundary and ρ c is its dimensionless representation, ρ c = r c / b 0 . The distribution of the radial and circumferential stresses in the elastic region r c r b 0 (or ρ c ρ 1 ) is given as
σ r X = A ( 1 ρ 2 ) ρ 2 p , σ θ X = A ( 1 + ρ 2 ) ρ 2 p , A = ρ c 2 [ η 1 2 ( η sin ψ c 4 η 2 + cos ψ c ) sin ψ c 4 η 2 ] .
Let ψ c be the value of ψ at ρ = ρ c . One can find ψ c and ρ c as functions of p from the following equations:
ln ( ρ c a ) = η 1 4 η 2 2 [ 1 + η 1 ( η 1 η ) ] ψ c + + ( 2 η η 1 ) 2 [ 1 + η 1 ( η 1 η ) ] ln [ η 1 4 η 2 η 1 4 η 2 cos ψ c ( 2 η η 1 ) sin ψ c ] , sin ψ c 4 η 2 ( 1 ρ c 2 ) + η 1 2 ( η sin ψ c 4 η 2 + cos ψ c ) ( 1 + ρ c 2 ) = p .
Then, the radial distribution of the stresses σ r and σ θ at a given value of p immediately follows from Equations (15) and (16).
In the case of Yield Criterion (6), the distribution of the radial and circumferential stresses in the plastic region a ρ ρ c is given in parametric form as
σ r σ s = q 1 sin φ , σ θ σ s = q 2 sin φ cos φ , ln ( ρ a ) = 1 + 2 R ¯ 2 { φ + ( q 1 q 2 ) ln [ cos φ + ( q 1 q 2 ) sin φ ] } , q 1 = 1 + R ¯ 1 + 2 R ¯ , q 2 = R ¯ 1 + 2 R ¯ ,
where φ is an auxiliary variable. The distribution of the radial and circumferential stresses in the elastic region ρ c ρ 1 is given as
σ r σ s = A 1 ρ 2 B 1 , σ θ σ s = A 1 ρ 2 B 1 , B 1 = 1 2 [ ( q 1 + q 2 ) sin φ c cos φ c ] , A 1 = ρ c 2 2 [ ( q 1 q 2 ) sin φ c + cos φ c ] ,
where φ c is the value of φ at ρ = ρ c . It is seen from Equation (18) that
ln ( ρ c a ) = 1 + 2 R ¯ 2 { φ c + ( q 1 q 2 ) ln [ cos φ c + ( q 1 q 2 ) sin φ c ] } .
Moreover, it follows from Equations (2), (14) and (19) that
B 1 = A 1 + p s .
Substituting Equation (19) into Equation (21) leads to
( q 1 + q 2 ) sin φ c cos φ c = ρ c 2 [ ( q 1 q 2 ) sin φ c + cos φ c ] + 2 p s .
The solution of Equations (20) and (22) supplies the dependencies of φ c and ρ c on p. Then, the radial distribution of the stresses σ r and σ θ at a given value of p immediately follows from Equations (18) and (19).
The radial stress increment at the end of unloading should satisfy the boundary conditions:
Δ σ r = 0
for ρ = a and
Δ σ r = p 0
for ρ = 1 . Assuming that the process of unloading is purely elastic, one can find from the equilibrium equations and Hooke’s law that the solution for the radial and circumferential stress increments satisfying the boundary conditions (23) and (24) is
Δ σ r = p 0 ( 1 a 2 ) ( 1 a 2 ρ 2 ) , Δ σ θ = p 0 ( 1 a 2 ) ( 1 + a 2 ρ 2 ) .
This solution is independent of the yield criterion accepted for the stage of loading.
Let σ r f and σ θ f be the radial and circumferential stresses at the end of loading, respectively. The residual stresses are determined from the equations σ r r e s = σ r f + Δ σ r and σ θ r e s = σ θ f + Δ σ θ . Then, using Equations (14) and (25) one can find the radial distribution of residual stresses as
σ r r e s X = σ r f X + p ( 1 a 2 ) ( 1 a 2 ρ 2 ) , σ θ r e s X = σ θ f X + p ( 1 a 2 ) ( 1 + a 2 ρ 2 )
in the case of Yield Criterion (3) and
σ r r e s σ s = σ r f σ s + p s ( 1 a 2 ) ( 1 a 2 ρ 2 ) , σ θ r e s σ s = σ θ f σ s + p s ( 1 a 2 ) ( 1 + a 2 ρ 2 )
in the case of Yield Criterion (6). The radial and circumferential stresses on the right-hand side of Equation (26) should be eliminated using Equations (15) and (16), on the right-hand side of Equation (27) using Equations (18) and (19).
The solutions (26) and (27) are valid if the respective yield criterion is not violated in the range a ρ 1 . This condition is represented as
( G + H ) ( σ r r e s ) 2 2 H σ r r e s σ θ r e s + ( H + F ) ( σ θ r e s ) 2 1
in the case of Yield Criterion (3) and
( σ r r e s σ s ) 2 2 R ¯ ( 1 + R ¯ ) ( σ r σ s ) ( σ θ σ s ) + ( σ θ r e s σ s ) 2 1
in the case of Yield Criterion (6).

4. Strain Solution

The solution in the elastic region immediately follows from the stress solution and Hooke’s law. In particular, substituting Equation (16) into Equation (10) and using Equation (14) leads to
ε r e k = ε r k = A [ ( 1 ρ 2 ) + ν ( 1 + ρ 2 ) ] ρ 2 p ( 1 ν ) , ε θ e k = ε θ k = A [ ( 1 + ρ 2 ) + ν ( 1 ρ 2 ) ] ρ 2 p ( 1 ν ) , ε z e k = ε z k = 2 ν ( A + p )
in the range ρ c ρ 1 . One can eliminate A in Equation (30) using the third equation in (16). Equation (30) is valid for Yield Criterion (3). Substituting Equation (19) into Equation (10) and using Equation (14) leads to
ε r e k = ε r k = s [ A 1 ( 1 + ν ) ρ 2 B 1 ( 1 ν ) ] , ε θ e k = ε θ k = s [ A 1 ( 1 + ν ) ρ 2 + B 1 ( 1 ν ) ] , ε z e k = ε z k = 2 s ν B 1
in the range ρ c ρ 1 . One can eliminate A 1 and B 1 in Equation (31) using the third and fourth equations in (19). Equation (31) is valid for Yield Criterion (6).
One can eliminate λ in Equation (11) to arrive at
ε ˙ r p = [ ( G + H ) σ r H σ θ ] [ ( F + H ) σ θ H σ r ] ε ˙ θ p , ε ˙ z p = ( G σ r + F σ θ ) [ ( F + H ) σ θ H σ r ] ε ˙ θ p .
Using the parameters η and η 1 introduced after Equation (15), these equations transform to
ε ˙ r p = η 1 ( 2 η 1 σ r η σ θ ) ( 2 σ θ η η 1 σ r ) ε ˙ θ p , ε ˙ z p = [ ( 2 η η 1 ) σ θ + η 1 ( 2 η 1 η ) σ r ] ( 2 σ θ η η 1 σ r ) ε ˙ θ p .
Equations (15) and (32) combine to give
ε ˙ r p = η 1 2 ( 4 η 2 tan ψ η ) ε ˙ θ p , ε ˙ z p = 1 2 ( 2 η η 1 + η 1 4 η 2 tan ψ ) ε ˙ θ p .
The third equation in (15) shows that ψ is independent of the time. Therefore, Equation (33) can be immediately integrated to give
ε r p = η 1 2 ( 4 η 2 tan ψ η ) ε θ p , ε z p = 1 2 ( 2 η η 1 + η 1 4 η 2 tan ψ ) ε θ p .
The elastic strains in the plastic region are determined from Equations (10), (14) and (15) as
ε r e k = ν η 1 cos ψ + ( ν η η 1 2 ) 4 η 2 sin ψ , ε θ e k = ( 2 ν η η 1 ) 4 η 2 sin ψ η 1 cos ψ , ε z e k = ν [ ( 2 + η η 1 ) 4 η 2 sin ψ + η 1 cos ψ ] .
In the case under consideration, the equation of strain compatibility reduces to
ρ ε θ ρ + ε θ ε r = 0 .
Using Equation (15), one can replace differentiation with respect to ρ with differentiation with respect to ψ . As a result, Equation (36) becomes
[ η 1 4 η 2 ( 2 η η 1 ) tan ψ ] ε θ ψ + 2 ( ε θ ε r ) = 0 .
Substituting Equations (34) and (35) into Equation (13) yields
ε r k = ν η 1 cos ψ + ( ν η η 1 2 ) 4 η 2 sin ψ + η 1 2 ( 4 η 2 tan ψ η ) ε θ p k , ε θ k = ( 2 ν η η 1 ) 4 η 2 sin ψ η 1 cos ψ + ε θ p k , ε z k = ν [ ( 2 + η η 1 ) 4 η 2 sin ψ + η 1 cos ψ ] 1 2 ( 2 η η 1 + η 1 4 η 2 tan ψ ) ε θ p k .
Equations (37) and (38) combine to give
( ε θ p / k ) ψ + ( 2 + η η 1 η 1 4 η 2 tan ψ ) [ η 1 4 η 2 ( 2 η η 1 ) tan ψ ] ε θ p k cos ψ ( 2 + η η 1 η 1 4 η 2 tan ψ ) 4 η 2 = 0 .
This is a linear differential equation for ε θ p . The plastic strains vanish at the elastic/plastic boundary. Therefore, the boundary condition to Equation (39) is
ε θ p = 0
for ψ = ψ c . The solution of Equation (39) satisfying this boundary condition is
ε θ p k = 1 4 η 2 exp [ ψ c ψ ( 2 + η η 1 η 1 4 η 2 tan ω ) [ η 1 4 η 2 ( 2 η η 1 ) tan ω ] d ω ] × ψ c ψ cos μ ( 2 + η η 1 η 1 4 η 2 tan μ ) exp [ ψ c μ ( 2 + η η 1 η 1 4 η 2 tan ω ) [ η 1 4 η 2 ( 2 η η 1 ) tan ω ] d ω ] d μ .
Having found the dependence of ε θ p on ψ , one can determine the other plastic strains from Equation (34) and the total strains from Equation (38) as functions of ψ . These functions and the third equation in (15) combine to supply the dependence of all these strains on ρ in parametric form.
Turning to Yield Criterion (6), one can eliminate λ 1 in Equation (12) to arrive at
ε ˙ r p = [ σ r + R ¯ ( σ r σ θ ) ] [ σ θ + R ¯ ( σ θ σ r ) ] ε ˙ θ p , ε ˙ z p = ( σ r + σ θ ) [ σ θ + R ¯ ( σ θ σ r ) ] ε ˙ θ p .
Equations (18) and (42) combine to give
ε ˙ r p = ( R ¯ cos φ 1 + 2 R ¯ sin φ ) ( 1 + R ¯ ) cos φ ε ˙ θ p , ε ˙ z p = ( cos φ + 1 + 2 R ¯ sin φ ) ( 1 + R ¯ ) cos φ ε ˙ θ p .
The third equation in (18) shows that φ is independent of the time. Therefore, Equation (43) can be immediately integrated to give
ε r p = ( R ¯ cos φ 1 + 2 R ¯ sin φ ) ( 1 + R ¯ ) cos φ ε θ p , ε z p = ( cos φ + 1 + 2 R ¯ sin φ ) ( 1 + R ¯ ) cos φ ε θ p .
The elastic strains in the plastic region are determined from Equations (10), (14) and (18) as
ε r e k s = ( ν q 2 q 1 ) sin φ + ν cos φ , ε θ e k s = ( ν q 1 q 2 ) sin φ cos φ , ε z e k s = ν ( 1 + 2 R ¯ sin φ + cos φ ) .
Equation (36) is valid. Using Equation (18), one can replace differentiation with respect to ρ with differentiation with respect to φ . As a result, Equation (36) becomes
( 1 + 2 R ¯ tan φ ) ( 1 + R ¯ ) ε θ φ + ε θ ε r = 0 .
Substituting Equations (44) and (45) into Equation (13) yields
ε r k s = ( ν q 2 q 1 ) sin φ + ν cos φ ( R ¯ cos φ 1 + 2 R ¯ sin φ ) ( 1 + R ¯ ) cos φ ε θ p k s , ε θ k s = ( ν q 1 q 2 ) sin φ cos φ + ε θ p k s , ε z k s = ν ( 1 + 2 R ¯ sin φ + cos φ ) ( cos φ + 1 + 2 R ¯ sin φ ) ( 1 + R ¯ ) cos φ ε θ p k s .
Equations (46) and (47) combine to give
ε θ p k φ + 1 + 2 R ¯ ε θ p k + s ( sin φ 1 + 2 R ¯ cos φ ) = 0 .
This is a linear differential equation for ε θ p . The plastic strains vanish at the elastic/plastic boundary. Therefore, the boundary condition to Equation (48) is
ε θ p = 0
for φ = φ c . The solution of Equation (48) satisfying this boundary condition is
ε θ p k s = cos φ exp [ 1 + 2 R ¯ ( φ c φ ) ] cos φ c .
Having found the dependence of ε θ p on φ , one can determine the other plastic strains from Equation (44) and the total strains from Equation (47) as functions of φ . These functions and the third equation in (18) combine to supply the dependence of all these strains on ρ in parametric form.
In the course of unloading, the increments of the strain components follow Hooke’s law:
Δ ε r = Δ σ r ν Δ σ θ E , Δ ε θ = Δ σ θ ν Δ σ r E , Δ ε z = ν ( Δ σ r + Δ σ θ ) E .
Using Equations (14) and (25), one can transform this equation to
Δ ε r k = p ( 1 a 2 ) [ 1 ν ( 1 + ν ) a 2 ρ 2 ] , Δ ε θ k = p ( 1 a 2 ) [ 1 ν + ( 1 + ν ) a 2 ρ 2 ] , Δ ε z k = 2 ν p ( 1 a 2 ) .
Using the strain solutions above and Equation (52), one can find the distribution of the residual strains as
ε r r e s k = ε r f k + Δ ε r k = ε r f k + p ( 1 a 2 ) [ 1 ν ( 1 + ν ) a 2 ρ 2 ] , ε θ r e s k = ε θ f k + Δ ε θ k = ε θ f k + p ( 1 a 2 ) [ 1 ν + ( 1 + ν ) a 2 ρ 2 ] , ε z r e s k = ε z f k + Δ ε z k = ε z f k 2 ν p ( 1 a 2 ) .
Here, the superscript f denotes the radial distribution of the strains at the end of loading.

5. Numerical Example and Discussion

The numerical example provided in this section evaluates the effect of replacing Hill’s quadratic yield criterion with several transversely isotropic yield criteria for a disk of inner radius a = 0.35 made of aluminum alloy 5352. This alloy’s mechanical properties involved in Yield Criterion (3) are available in the literature [34]. One of the principal axes of anisotropy coincides with the axial direction independently of the orientation of the two other axes. The left-hand side of Yield Criterion (3) is known from [34]. Then, Z = 220 MPa and the yield stresses in the directions of the other two principal axes of anisotropy are 198 and 279 MPa. The corresponding Lankford’s coefficients are 0.53 and 2.27.
Using these data and the formulation of the boundary value problem in Section 2, one can arrive at two different Hill’s yield criteria depending on the orientation of the principal axes of anisotropy relative to the radial direction. In particular, X = 198 MPa, R0 = 0.53, and R90 = 2.27 in Case (i) and X = 279 MPa, R0 = 2.27, and R90 = 0.53 in Case (ii). Using these parameters and Equation (5), one can find Hill’s locus from Yield Criterion (3) for each case. Then, the Lankford’s coefficient and the yield stresses in any direction can be found by means of any of these loci. Case (i) has been used. As a result, if Δα = 15° in Equations (7)-(8), R15 = 0.46, R30 = 0.65, R45 = 1.10, R60 = 1.41, and R75 = 1.59; σS15 = 199 MPa, σS30 = 203 MPa, σS45 = 215 MPa, σS60 = 236 MPa, and σS75 = 264 MPa. The elastic properties are E = 68.5 GPa and ν = 0.33. The mechanical properties responsible for plastic anisotropy are summarized in Table 1.
It is worthy of note that it is assumed in Equation (13) that X = 198 MPa independently of the yield criterion adopted. This choice allows one to compare the solutions at the same value of p0 rather than p.
The solution for Yield Criterion (3) is valid for Cases (i) and (ii). The input data for the solution for Yield Criterion (6) have been derived from Equations (7) and (8) for three n-values (Table 1). In particular,
R ¯ = R ¯ I = R 0 + R 90 2 , R ¯ = R ¯ I I = R 0 + 2 R 45 + R 90 4 , R ¯ = R ¯ I I I = R 0 + 2 ( R 15 + R 30 + R 45 + R 60 + R 75 ) + R 90 12 ,
σ ¯ s = σ ¯ s I = X + Y 2 , σ ¯ s = σ ¯ s I I = X + 2 σ s 45 + Y 4 , σ ¯ s = σ ¯ s I I I = X + 2 ( σ s 15 + σ s 30 + σ s 45 + σ s 60 + σ s 75 ) + Y 12 .
The radius of the elastic/plastic boundary at the end of loading is an essential parameter for autofrettage [35]. Even though the system of loading in the case of hydraulic autofrettage is different from that considered in the present paper, it is of interest to understand the dependence of this radius on the yield criterion adopted. The variation of ρ c with p found with the use of Equations (16) and (19) is depicted in Figure 2. It is seen from this figure that the effect of the yield criterion on ρ c is significant and increases as p increases. Thus, the difference between the plastic region’s predicted sizes found using the actual orthotropic yield criterion and the assumption of normal anisotropy becomes larger as the pressure increases.
In what follows, it is assumed that the disk is loaded with the pressure p0 = 140 MPa. The effect of the yield criterion on the distribution of stresses has been discussed in [31]. Therefore, the present paper is restricted to evaluating this effect on the distribution of strains.
Figure 3 shows the variation of the principal strains with the dimensionless radius at the end of loading found using the solutions given in Section 4. The distribution of the radial strain is the least sensitive to the choice of the yield criterion (Figure 3c). The distribution of the circumferential strain is more sensitive because of the solution for case (ii) (Figure 3b). The choice of the yield criterion has the most significant effect on the distribution of the axial strain (Figure 3a). The effect is pronounced in the vicinity of the inner radius of the disk.
The effect of the yield criterion on the distribution of the residual strains found from Equation (52) is shown in Figure 4. It is seen from Figure 3 and Figure 4 that the residual strains are more strongly affected than the strains at the end of loading. In particular, there is a significant deviation of all residual strains for Case (ii) from four other cases. However, the difference in the residual axial strain for the different yield criteria is negligible in the common elastic region (Figure 4a).

6. Conclusions

Many metallic materials have different yield stresses in the directions of three principal axes of anisotropy. In the case of a sheet, the different yield stresses in its plane are often replaced with average yield stress independent of the direction. It has been shown in the present paper that this simplification in the theoretical description of material behavior can result in a very significant difference in the distribution of strains and residual strains even in the case of infinitesimal strains. It is expected that this difference will be even more considerable in the case of using the assumption of transversely isotropic yield criteria in the modeling of metal forming processes at large strains.

Author Contributions

Data curation, S.S.; Formal analysis, S.A.; Investigation, Y.E.; Supervision, L.L.; Validation, Y.E.; Visualization, S.S.; Writing—original draft, S.A.; Writing—review & editing, L.L. All authors have participated in the research and in the writing of this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation, grant number 20-79-10340.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a0, b0inner and outer radius of hollow disk
p 0 external pressure
( r , θ , z ) cylindrical coordinate system
σ r , σ θ , σ z normal stresses
X, Y, Zyield stresses in tension in the radial, circumferential and axial directions
F, G, HHill’s coefficients
R 0 , R 90 Lankford’s coefficients in the radial and circumferential directions
σ ¯ s , R ¯ average yield stress in tension and Lankford’s coefficient
R Δ α i , σ s Δ α i Lankford’s coefficient and yield stress in tension in the direction γ = Δ α i
nnumber of directions in which the Lankford’s coefficients and yield stresses are measured
ε r , ε θ , ε z
ε r e , ε θ e , ε z e
ε r p , ε θ p , ε z p
normal strains (the superscript e denotes the elastic portion of the strain components, the superscript p denotes the plastic portion of the strain components)
ε ˙ r p , ε ˙ θ p , ε ˙ z p normal plastic strain rates
EYoung’s modulus
ν Poisson’s ratio
λ , λ 1 non-negative multipliers in Equations (10) and (11)
ρ , a , k , s , p dimensionless quantities introduced in Equation (13)
η , η 1 parameters introduced in Equation (14)
ψ , φ auxiliary variables
r c , ρ c radius of the elastic/plastic boundary and its dimensionless representation
ψ c , φ c values of ψ and φ at ρ = ρ c
A parameter introduced in Equation (15)
q 1 , q 2 parameters introduced in Equation (17)
A 1 , B 1 parameters introduced in Equation (18)
Δ σ r , Δ σ θ stress increments at the end of unloading
σ r f , σ θ f stresses at the end of loading
σ r r e s , σ θ r e s residual stresses
μ , ω integration variables
Δ ε r , Δ ε θ , Δ ε θ strain increments at the end of unloading
ε r f , ε θ f , ε z f strains at the end of loading
ε r r e s , ε θ r e s , ε z r e s residual strains

References

  1. Rees, D. Basic Engineering Plasticity; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar] [CrossRef]
  2. Banabic, D. Sheet Metal Forming Processes: Constitutive Modeling and Numerical Simulation; Springer: Heidelberg, Germany, 2010. [Google Scholar] [CrossRef]
  3. Wagoner, R.; Li, M. Simulation of springback: Through-thickness integration. Int. J. Plast. 2007, 23, 345–360. [Google Scholar] [CrossRef]
  4. Prime, M.B. Amplified effect of mild plastic anisotropy on residual stress and strain anisotropy. Int. J. Solids Struct. 2017, 70–77. [Google Scholar] [CrossRef]
  5. Essa, S.; Argeso, H. Elastic analysis of variable profile and polar orthotropic FGM rotating disks for a variation function with three parameters. Acta Mech. 2017, 228, 3877–3899. [Google Scholar] [CrossRef]
  6. Yıldırım, V. Numerical/analytical solutions to the elastic response of arbitrarily functionally graded polar orthotropic rotating discs. J. Braz. Soc. Mech. Sci. Eng. 2018, 40, 320. [Google Scholar] [CrossRef]
  7. Jeong, W.; Alexandrov, S.; Lang, L. Effect of plastic anisotropy on the distribution of residual stresses and strains in rotating annular disks. Symmetry 2018, 10, 420. [Google Scholar] [CrossRef] [Green Version]
  8. Wang, C.; Kinzel, G.; Altan, T. Mathematical modeling of plane-strain bending of sheet and plate. J. Mater. Process. Technol. 1993, 39, 279–304. [Google Scholar] [CrossRef]
  9. Pourboghrat, F.; Chung, K.; Richmond, O. A hybrid membrane/shell method for rapid estimation of Springback in anisotropic sheet metals. J. Appl. Mech. 1998, 65, 671–684. [Google Scholar] [CrossRef]
  10. Chakrabarty, J.; Bin Lee, R.; Chan, K. An exact solution for the elastic/plastic bending of anisotropic sheet metal under conditions of plane strain. Int. J. Mech. Sci. 2001, 43, 1871–1880. [Google Scholar] [CrossRef]
  11. Leu, D.-K. Relationship between mechanical properties and geometric parameters to limitation condition of springback based on springback–radius concept in V-die bending process. Int. J. Adv. Manuf. Technol. 2018, 101, 913–926. [Google Scholar] [CrossRef]
  12. Salimi, M.; Jamshidian, M.; Beheshti, A.; Sadeghi Dolatabadi, A. Bending-unbending analysis of anisotropic sheet under plane strain condition. J. Adv. Mat. Eng. 2008, 26, 77–86. [Google Scholar]
  13. Li, D.; Luo, Y.; Peng, Y.; Hu, P. The numerical and analytical study on stretch flanging of V-shaped sheet metal. J. Mater. Process. Technol. 2007, 189, 262–267. [Google Scholar] [CrossRef]
  14. Paul, S.K. Non-linear correlation between uniaxial tensile properties and shear-edge hole expansion ratio. J. Mater. Eng. Perform. 2014, 23, 3610–3619. [Google Scholar] [CrossRef]
  15. Kim, J.H.; Kwon, Y.J.; Lee, T.; Lee, K.-A.; Kim, H.S.; Lee, C.S. Prediction of hole expansion ratio for various steel sheets based on uniaxial tensile properties. Met. Mater. Int. 2018, 24, 187–194. [Google Scholar] [CrossRef]
  16. Marciniak, Z.; Kuczyński, K. Limit strains in the processes of stretch-forming sheet metal. Int. J. Mech. Sci. 1967, 9, 609–620. [Google Scholar] [CrossRef]
  17. Bressan, J.; Williams, J. The use of a shear instability criterion to predict local necking in sheet metal deformation. Int. J. Mech. Sci. 1983, 25, 155–168. [Google Scholar] [CrossRef]
  18. Bressan, J.D.; Wang, Q.; Simonetto, E.; Ghiotti, A.; Bruschi, S. Formability prediction of Ti6Al4V titanium alloy sheet deformed at room temperature and 600 °C. Int. J. Mater. Form. 2020, 1–15. [Google Scholar] [CrossRef]
  19. Liao, K.-C.; Pan, J.; Tang, S. Approximate yield criteria for anisotropic porous ductile sheet metals. Mech. Mater. 1997, 26, 213–226. [Google Scholar] [CrossRef]
  20. Huang, H.-M.; Pan, J.; Tang, S. Failure prediction in anisotropic sheet metals under forming operations with consideration of rotating principal stretch directions. Int. J. Plast. 2000, 16, 611–633. [Google Scholar] [CrossRef]
  21. Chien, W.Y.; Pan, J.; Tang, S.C. Modified anisotropic Gurson yield criterion for porous ductile sheet metals. J. Eng. Mater. Technol. 2000, 123, 409–416. [Google Scholar] [CrossRef]
  22. Harikrishna, C.; Davidson, M.J.; Nagaraju, C.; RatnaPrasad, A.V. Effect of lubrication and anisotropy on hardness in the upsetting test. Trans. Indian Inst. Met. 2015, 69, 1449–1457. [Google Scholar] [CrossRef]
  23. Yang, X.-Y.; Yang, C.; Liang, F.; Zhang, W.-X.; Niu, Y. Influence of normal anisotropy coefficient on no-rivet connection quality. Suxing Gongcheng Xuebao/J. Plast. Eng. 2019, 26, 245–250. [Google Scholar] [CrossRef]
  24. Uscinowicz, R. Characterization of directional Elastoplastic properties of Al/Cu bimetallic sheet. J. Mater. Eng. Perform. 2019, 28, 1350–1359. [Google Scholar] [CrossRef]
  25. Ning, B.; Zhao, Z.; Mo, Z.; Wu, H.; Peng, C.; Gong, H. Influence of continuous annealing temperature on mechanical properties and texture of battery shell steel. Metals 2019, 10, 52. [Google Scholar] [CrossRef] [Green Version]
  26. Moura, A.N.; Ferreira, J.L.; Martins, J.B.R.; De Souza, M.V.; Castro, N.A.; Orlando, M.T.D. Microstructure, crystallographic texture, and stretch-Flangeability of hot-rolled multiphase steel. Steel Res. Int. 2020, 91. [Google Scholar] [CrossRef]
  27. Lequeu, P.; Jonas, J.J. Modeling of the plastic anisotropy of textured sheet. Met. Mater. Trans. A 1988, 19, 105–120. [Google Scholar] [CrossRef]
  28. Wang, X.; Cao, L.; Peng, X.; Zhang, J.; Zhuang, L.; Guo, M. Relationship among mechanical properties anisotropy, microstructure and texture in AA 6111 alloy sheets. J. Wuhan Univ. Technol. Sci. Ed. 2016, 31, 648–653. [Google Scholar] [CrossRef]
  29. Inoue, H.; Takasugi, T. Texture control for improving deep Drawability in rolled and annealed aluminum alloy sheets. Mater. Trans. 2007, 48, 2014–2022. [Google Scholar] [CrossRef] [Green Version]
  30. Serebryany, V.N.; Pozdnyakova, N.N. Evaluation of the normal anisotropy coefficient in AZ31 alloy sheets. Russ. Met. (Metally) 2009, 2009, 58–64. [Google Scholar] [CrossRef]
  31. Grechnikov, F.V.; Erisov, Y.A.; Alexandrov, S.E. Effect of the anisotropic yield condition on the predicted distribution of residual stresses in a thin disk. Dokl. Phys. 2019, 64, 233–237. [Google Scholar] [CrossRef]
  32. Cohen, T.; Masri, R.; Durban, D. Analysis of circular hole expansion with generalized yield criteria. Int. J. Solids Struct. 2009, 46, 3643–3650. [Google Scholar] [CrossRef] [Green Version]
  33. Hill, R. The Mathematical Theory of Plasticity; Oxford University Press: New York, NY, USA, 1950. [Google Scholar]
  34. Watson, M.; Dick, R.; Huang, Y.H.; Lockley, A.; Cardoso, R.; Santos, A. Benchmark 1—Failure prediction after cup drawing, reverse redrawing and expansion part A: Benchmark description. J. Phys. Conf. Ser. 2016, 734, 022001. [Google Scholar] [CrossRef] [Green Version]
  35. Dixit, U.S.; Kamal, S.M.; Shufen, R. Autofrettage Processes: Technology and Modelling; CRC Press: Boka Raton, FL, USA, 2019. [Google Scholar]
Figure 1. Schematic of the disk subjected to external pressure p 0 .
Figure 1. Schematic of the disk subjected to external pressure p 0 .
Symmetry 12 01834 g001
Figure 2. Dependence of ρ c on p for different yield criteria.
Figure 2. Dependence of ρ c on p for different yield criteria.
Symmetry 12 01834 g002
Figure 3. Radial distribution of axial (a), circumferential (b) and radial (c) strains.
Figure 3. Radial distribution of axial (a), circumferential (b) and radial (c) strains.
Symmetry 12 01834 g003
Figure 4. Radial distribution of residual axial (a), circumferential (b) and radial (c) strains.
Figure 4. Radial distribution of residual axial (a), circumferential (b) and radial (c) strains.
Symmetry 12 01834 g004
Table 1. Input data for all the cases considered.
Table 1. Input data for all the cases considered.
Hill’s Yield Criterion (Equation (3))Normal Anisotropy (Equation (6))
Case (i)Case (ii)n = 2n = 3n = 7
X, MPa (R0)198 (0.53)279 (2.27)238 (1.40)226 (1.25)225 (1.11)
Y, MPa (R90)279 (2.27)198 (0.53)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Erisov, Y.; Surudin, S.; Alexandrov, S.; Lang, L. Influence of the Replacement of the Actual Plastic Orthotropy with Various Approximations of Normal Anisotropy on Residual Stresses and Strains in a Thin Disk Subjected to External Pressure. Symmetry 2020, 12, 1834. https://doi.org/10.3390/sym12111834

AMA Style

Erisov Y, Surudin S, Alexandrov S, Lang L. Influence of the Replacement of the Actual Plastic Orthotropy with Various Approximations of Normal Anisotropy on Residual Stresses and Strains in a Thin Disk Subjected to External Pressure. Symmetry. 2020; 12(11):1834. https://doi.org/10.3390/sym12111834

Chicago/Turabian Style

Erisov, Yaroslav, Sergei Surudin, Sergei Alexandrov, and Lihui Lang. 2020. "Influence of the Replacement of the Actual Plastic Orthotropy with Various Approximations of Normal Anisotropy on Residual Stresses and Strains in a Thin Disk Subjected to External Pressure" Symmetry 12, no. 11: 1834. https://doi.org/10.3390/sym12111834

APA Style

Erisov, Y., Surudin, S., Alexandrov, S., & Lang, L. (2020). Influence of the Replacement of the Actual Plastic Orthotropy with Various Approximations of Normal Anisotropy on Residual Stresses and Strains in a Thin Disk Subjected to External Pressure. Symmetry, 12(11), 1834. https://doi.org/10.3390/sym12111834

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