Next Article in Journal
Comparison of the Cracking Behavior of Powder Metallurgy and Ingot Metallurgy Ti-5Al-5Mo-5V-3Cr Alloys during Hot Deformation
Previous Article in Journal
Modification of Residual Stresses in Laser Additive Manufactured AlSi10Mg Specimens Using an Ultrasonic Peening Technique
Previous Article in Special Issue
Description of Residual Stress and Strain Fields in FGM Hollow Disc Subject to External Pressure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Theory of Elastic/Plastic Plane Strain Pure Bending of FGM Sheets at Large Strain

1
School of Mechanical Engineering and Automation, Beihang University, No. 37 Xueyuan Road, Beijing 100191, China
2
Ishlinsky Institute for Problems in Mechanics, 101-1 Prospect Vernadskogo, 119526 Moscow, Russia
3
Department of Civil Engineering, National Cheng Kung University, 1 University Road, Tainan 70101, Taiwan
*
Author to whom correspondence should be addressed.
Materials 2019, 12(3), 456; https://doi.org/10.3390/ma12030456
Submission received: 31 December 2018 / Revised: 24 January 2019 / Accepted: 28 January 2019 / Published: 1 February 2019

Abstract

:
An efficient analytical/numerical method has been developed and programmed to predict the distribution of residual stresses and springback in plane strain pure bending of functionally graded sheets at large strain, followed by unloading. The solution is facilitated by using a Lagrangian coordinate system. The study is concentrated on a power law through thickness distribution of material properties. However, the general method can be used in conjunction with any other through thickness distributions assuming that plastic yielding initiates at one of the surfaces of the sheet. Effects of material properties on the distribution of residual stresses are investigated.

1. Introduction

Structures made of functionally graded materials (FGM) are advantageous for many applications. A difficulty with theoretical analysis and design is that structures made of FGM are classified by a much greater number of parameters than similar structures made of homogeneous materials. For this reason, it is desirable to perform parametric studies by analytic or semi-analytic methods as much as possible. A review of results related to the analysis of FGM and published before 2007 is presented in [1]. This review focuses on structures with through-thickness variation of material properties. Analytic solutions derived in [1,2,3,4,5] belong to this class of FGM as well. In [2,3,4], elastic and elastic/plastic spherical vessels subjected to various loading conditions are considered. Thermo-elastic simply supported and clamped circular plates are studied in [5]. Many analytic and semi-analytic solutions are available for FGM discs and cylinders assuming that material properties vary in the radial direction but are independent of the circumferential and axial directions. Purely elastic solutions for a hollow disc or cylinder subjected to internal or/and external pressure are derived in [6,7,8]. An axisymmetric thermo-elastic solution for a hollow cylinder subjected quite a general system of thermo-mechanical loading is presented in [9]. It is assumed that the temperature varies along the radial coordinate. A plane strain analytic elastic/plastic solution for pressurized tubes is found in [10]. The solution is based on the Tresca yield criterion. Many solutions are proposed for functionally graded solid and hollow rotating discs. Purely elastic solutions for solid discs of constant thickness are given in [11,12], a purely elastic solution for a hollow disc of variable thickness in [13], a purely elastic solution for hollow polar orthotropic discs in [14], and a solution for hollow cylinders using the theory of electrothermoelasticity in [15]. An elastic perfectly plastic stress solution for hollow discs is derived in [16] using the von Mises yield criterion.
All of the aforementioned solutions deal with infinitesimal strain. A distinguished feature of the solution provided in the present paper is that strains are large. The process considered is pure bending of a FGM sheet under plane strain conditions. A review on bending of functionally graded sheets and beams at infinitesimal strains is given in [17]. The present solution is based on the approach proposed in [18]. It is shown in this paper that the use of Lagrangian coordinates facilitates the solution. Moreover, the equations describing kinematics can be solved independently of stress equations in the case of isotropic incompressible material. This is an advantage as compared to the classic approach developed in [19] where the stress equations are solved first. The classic approach is restricted to perfectly plastic materials, whereas the mapping in Equation (1) is valid for a large class of constitutive equations. The approach proposed in [18] has already been successfully extended to more general constitutive equations in [20,21,22,23]. It is shown in the present paper that the approach is also efficient for FGM sheets. It is worth noting that a rigid plastic solution for pure bending of laminated sheets (such sheets can also be referred to as functionally graded sheets) at large strain is given in [24].

2. Basic Equations

The process of plane strain pure bending is illustrated in Figure 1. The approach proposed in [18] for solving the corresponding boundary value problem is based on the following transformation equations:
x H = ζ a + s a 2 cos 2 a η s a , y H = ζ a + s a 2 sin 2 a η .
where ( x , y ) is an Eulerian–Cartesian coordinate system and ( ζ , η ) is a Lagrangian coordinate system. Without loss of generality, it is possible to assume that the origin of the Cartesian coordinate system is located at the intersection of the axis of symmetry of the process and the outer surface A B and that the x-axis coincides with the axis of symmetry. The Lagrangian coordinate system is chosen such that
ζ = x / H and η = y / H
at the initial instant where H is the initial thickness of the sheet. It is evident from these relations and the geometry in Figure 1 that ζ = 0 on AB, ζ = 1 on CD, η = L / H on CB and η = L / H on AD throughout the process of deformation. Here, L is the initial width of the sheet. In Equation (1), a is a time-like variable. In particular, a = 0 at the initial instant. In Equation (1), s is a function of a. This function should be found from the stress solution and therefore depends on constitutive equations. The condition in Equation (2) is satisfied if
s = 1 4
at a = 0 . It is possible to verify by inspection that the mapping in Equation (1) satisfies the equation of incompressibility. Moreover, this mapping transforms initially straight lines A 1 B 1 and C 1 D 1 into circular arcs AB and CD and initially straight lines C 1 B 1 and A 1 D 1 into circular arcs CB and AD after any amount of deformation (Figure 1). Furthermore, coordinate curves of the Lagrangian coordinate system coincide with trajectories of the principal strain rates and, for coaxial models, with trajectories of the principal stresses. Thus, the shear stress vanishes in the Lagrangian coordinates. In particular, the contour ABCD is free of shear stresses. Let σ ζ and σ η be the physical stress components referred to the Lagrangian coordinates. The stress solution should satisfy the boundary conditions
σ ζ = 0
for ζ = 1 and ζ = 0 . The only non-trivial equilibrium equation in the Lagrangian coordinates has been derived in [18] as
σ ζ ζ + a σ ζ σ η 2 ζ a + s = 0 .
The initial plane strain yield criterion of the functionally graded sheet is supposed to be
σ ζ σ η = 2 3 σ 0 β x H .
where σ 0 is a material constant and β ( x / H ) is an arbitrary function of its argument. It is assumed that material properties are not affected by plastic deformation. Therefore, Equation (6) can be rewritten in the form
σ ζ σ η = 2 3 σ 0 β ζ .
In this case, the yield locus is invariant along the motion. The importance of this property of material models has been emphasized in [25]. Let τ ζ and τ η be the deviatoric portions of σ ζ and σ η , respectively. Since the material is incompressible, τ ζ + τ η = 0 under plane strain conditions. Then, the yield criterion in Equation (7) is equivalent to
τ ζ = τ η = σ 0 β ζ 3 .
Hooke’s law generalized on functionally graded materials reads
τ ζ = 2 G 0 g ζ ε ζ e , τ η = 2 G 0 g ζ ε η e .
It has been taken into account here that Poisson’s ratio is equal to 1/2 for incompressible materials. In addition, ε ζ e and ε η e are the total strain components in elastic regions and the elastic portions of the total strain components in plastic regions referred to the Lagrangian coordinate system, G 0 is a material constant and g ( ζ ) is an arbitrary function of its argument.
Geometric parameters shown in Figure 1 depend on a and are expressed as [18]
r AB H = s a , r CD H = s a a , θ 0 = 2 a L H , h H = s s a a .
Once s has been found as a function of a, these parameters are immediate from Equation (10).

3. Stress Solution at Loading

It is assumed that the functions β ( ζ ) and g ( ζ ) involved in Equations (7) and (9) are such that plastic yielding can only initiate at ζ = 0 or ζ = 1 . This assumption can be verified using the purely elastic solution with no difficulty. At the very beginning of the process, the entire sheet is elastic. As deformation proceeds, one of the following three cases arises: (i) plastic yielding initiates at the surface ζ = 1 ; (ii) plastic yielding initiates at the surface ζ = 0 ; and (iii) plastic yielding initiates simultaneously at the surfaces ζ = 1 and ζ = 0 . These cases should be treated separately. In the following, ζ 1 is the elastic/plastic boundary between the plastic region that propagates from the surface ζ = 0 and the elastic region and ζ 2 is the elastic/plastic boundary between the plastic region that propagates from the surface ζ = 1 and the elastic region. It is evident that both ζ 1 and ζ 2 depend on a. The general structure of the solution with two plastic regions is illustrated in Figure 2. Let M be the bending moment. Then, its dimensionless representation is in terms of the Lagrangian coordinates given by [18]
m = 2 3 M σ 0 H 2 = 3 a 1 0 σ η σ 0 d ζ .
In the elastic region, the whole strain is elastic. Therefore, it follows from Equation (1) that the principal logarithmic strains are
2 ε ζ e = 2 ε η e = ln 4 ζ a + s .
Since σ ζ σ η = τ ζ τ η , Equations (5) and (9) combine to give
σ ζ ζ + G 0 a g ζ ζ a + s ε ζ e ε η e = 0 .
Eliminating the strain components in Equation (13) by means of Equation (12) results in
σ ζ ζ G 0 a g ζ ζ a + s ln 4 ζ a + s = 0 .
Integrating this equation with respect to ζ and using the boundary condition in Equation (4) at ζ = 0 leads to
σ ζ σ 0 = a 3 k 0 ζ g χ ln 4 χ a + s χ a + s d χ , σ η σ 0 = σ ζ σ 0 + 2 3 k g ζ ln 4 ζ a + s ,
where k = σ 0 / ( 3 G 0 ) and χ is a dummy variable of integration. The expression for σ η in Equation (15) has been derived using the identity σ η = σ ζ τ ζ + τ η , and Equations (9) and (12). In the case of the purely elastic solution, Equation (15) must satisfy the boundary condition in Equation (4) at ζ = 1 . Then, the equation for the function s(a) is
1 0 g χ ln 4 χ a + s χ a + s d χ = 0 .
Using Equation (15), in which s should be eliminated by means of the solution of Equation (16), and the yield criterion in Equation (8), it is possible to determine which of the three cases mentioned above occurs for given material properties. Simultaneously, the value of a at which plastic yielding initiates is determined. This value of a is denoted as a e . In the following, it is assumed that a a e . It is now necessary to consider Cases (i), (ii) and (iii) separately.
Case (i). There are two regions. A plastic region occupies the domain 1 ζ ζ 2 and an elastic region the domain ζ 2 ζ 0 . Equation (15) is valid in the elastic region. However, the function s ( a ) is not determined from Equation (16). It is reasonable to assume that σ η < σ ζ in the plastic region. Therefore, the yield criterion in Equation (7) becomes
σ ζ σ η = 2 3 σ 0 β ζ .
Substituting Equation (17) into Equation (5) and integrating yields the dependence of the stress σ ζ on ζ . Using Equation (17) again provides the dependence of the stress σ η on ζ . As a result,
σ ζ σ 0 = a 3 1 ζ β χ χ a + s d χ , σ η σ 0 = σ ζ σ 0 2 3 β ζ .
It is evident that this solution satisfies the boundary condition in Equation (4) at ζ = 1 . Both σ ζ and σ η should be continuous across ζ = ζ 2 . Consequently, τ ζ is continuous across ζ = ζ 2 . The stress τ ζ on the elastic side of the elastic/plastic boundary is determined from Equation (15) and on the plastic side from Equation (8). Then, the condition of continuity of τ ζ across the surface ζ = ζ 2 is represented as
g ζ 2 ln 4 ζ 2 a + s = 3 k β ζ 2 .
Solving this equation for s yields
s = 1 4 exp 3 k β ζ 2 g ζ 2 ζ 2 a .
Using Equations (15) and (18), the condition of continuity of σ ζ across the surface ζ = ζ 2 is represented as
0 ζ 2 g χ ln 4 χ a + s χ a + s d χ = 3 k 1 ζ 2 β χ χ a + s d χ .
In this equation, s can be eliminated by means of Equation (20). The resulting equation should be solved numerically to find ζ 2 as a function of a. Then, s as a function of a is readily found from Equation (20). The yield criterion should be checked in the elastic region using the solution in Equation (15). The calculation should be stopped when the yield condition is satisfied at one point of the elastic region. Denote the corresponding value of a as a 2 .
In Case (i), Equation (11) becomes
m = 3 a 1 ζ 2 σ η σ 0 d ζ + 3 a ζ 2 0 σ η σ 0 d ζ .
In the first integrand, σ η / σ 0 should be eliminated by means of Equation (18) and in the second by means of Equation (15).
Case (ii). There are two regions. A plastic region occupies the domain ζ 1 ζ 0 and an elastic region the domain 1 ζ ζ 1 . The elastic solution in Equation (15) satisfies the boundary condition in Equation (4) at ζ = 0 . Therefore, it is convenient to rewrite this solution as
σ ζ σ 0 = a 3 k 1 ζ g χ ln 4 χ a + s χ a + s d χ , σ η σ 0 = σ ζ σ 0 + 2 3 k g ζ ln 4 ζ a + s .
The elastic solution in this form satisfies the boundary condition in Equation (4) at ζ = 1 . It is reasonable to assume that σ η > σ ζ in the plastic region. Therefore, the yield criterion in Equation (7) becomes
σ ζ σ η = 2 3 σ 0 β ζ .
Substituting Equation (24) into Equation (5) and integrating yields the dependence of the stress σ ζ on ζ . Using Equation (24) again provides the dependence of the stress σ η on ζ . As a result,
σ ζ σ 0 = a 3 0 ζ β χ χ a + s d χ , σ η σ 0 = σ ζ σ 0 + 2 3 β ζ .
It is evident that this solution satisfies the boundary condition in Equation (4) at ζ = 0 . Both σ ζ and σ η should be continuous across ζ = ζ 1 . Consequently, τ ζ is continuous across ζ = ζ 1 . The stress τ ζ on the elastic side of the elastic/plastic boundary is determined from Equation (23) and on the plastic side from Equation (8). Then, the condition of continuity of τ ζ across the surface ζ = ζ 1 is represented as
g ζ 1 ln 4 ζ 1 a + s = 3 k β ζ 1 .
Solving this equation for s yields
s = 1 4 exp 3 k β ζ 1 g ζ 1 ζ 1 a .
Using Equations (23) and (25), the condition of continuity of σ ζ across the surface ζ = ζ 1 is represented as
1 ζ 1 g χ ln 4 χ a + s χ a + s d χ = 3 k 0 ζ 1 β χ χ a + s d χ .
In this equation, s can be eliminated by means of Equation (27). The resulting equation should be solved numerically to find ζ 1 as a function of a. Then, s as a function of a is readily found from Equation (27). The yield criterion should be checked in the elastic region using the solution in Equation (23). The calculation should be stopped when the yield condition is satisfied at one point of the elastic region. Denote the corresponding value of a as a 1 .
In Case (ii), Equation (11) becomes
m = 3 a 1 ζ 1 σ η σ 0 d ζ + 3 a ζ 1 0 σ η σ 0 d ζ .
In the first integrand, σ η / σ 0 should be eliminated by means of Equation (23) and in the second by means of Equation (25).
Case (iii). In this case, there are two plastic regions, 1 ζ ζ 2 and ζ 1 ζ 0 , and one elastic region, ζ 1 ζ ζ 2 . At the beginning of this stage of the process, a = a 1 and ζ 2 = 1 or a = a 2 and ζ 1 = 0 . Let σ n 1 be the value of σ ζ at ζ = ζ 1 and σ n 2 be the value of σ ζ at ζ = ζ 2 . Then, the elastic solution in Equation (15) can be rewritten as
σ ζ σ 0 = a 3 k ζ 1 ζ g χ ln 4 χ a + s χ a + s d χ + σ n 1 σ 0 , σ η σ 0 = σ ζ σ 0 + 2 3 k g ζ ln 4 ζ a + s .
It follows from this solution that
σ n 2 σ 0 = a 3 k ζ 1 ζ 2 g χ ln 4 χ a + s χ a + s d χ + σ n 1 σ 0 .
The solution in Equation (18) is valid in the plastic region 1 ζ ζ 2 and the solution in Equation (25) in the plastic region ζ 1 ζ 0 . Then,
σ n 2 σ 0 = a 3 1 ζ 2 β χ χ a + s d χ ,
and
σ n 1 σ 0 = a 3 0 ζ 1 β χ χ a + s d χ .
Equations (20) and (27) are valid. Therefore,
exp 3 k β ζ 1 g ζ 1 4 ζ 1 a = exp 3 k β ζ 2 g ζ 2 4 ζ 2 a ,
and
a = 1 4 ζ 2 ζ 1 exp 3 k β ζ 2 g ζ 2 exp 3 k β ζ 1 g ζ 1 .
Equations (31)–(33) combine to give
1 ζ 2 β χ χ a + s d χ + 1 3 k ζ 1 ζ 2 g χ ln 4 χ a + s χ a + s d χ + 0 ζ 1 β χ χ a + s d χ = 0 .
Eliminating in this equation s by means of Equation (20) or Equation (27) and then a by means of Equation (35) supplies the equation to find ζ 1 as a function of ζ 2 (or ζ 2 as a function of ζ 1 ). Then, a as a function of ζ 1 (or ζ 2 ) is found from Equation (35) and s as a function of ζ 1 (or ζ 2 ) from Equation (20) or (27). The distribution of the stresses is determined from Equation (30) with the use of Equations (32) and (33) in the elastic region, from Equation (18) in the region 1 ζ ζ 2 and from Equation (25) in the region ζ 1 ζ 0 .
In Case (iii), Equation (11) becomes
m = 3 a 1 ζ 2 σ η σ 0 d ζ + 3 a ζ 2 ζ 1 σ η σ 0 d ζ + 3 a ζ 1 0 σ η σ 0 d ζ .
In the first integrand, σ η / σ 0 should be eliminated by means of Equation (18), in the second by means of Equation (30) and the third by means of Equation (25). As usual, it is necessary to verify that the yield criterion is not violated in the elastic region.

4. Unloading

It is assumed that unloading is purely elastic. This assumption should be verified a posteriori. At this stage of the process, the strains can be considered as infinitesimal. Let a f and s f be the values of a and s, respectively, at the end of loading. These values are known from the solution given in the previous section. Using Equation (10), the values of r CD and r AB at the end of loading, r CD f and r AB f , are determined as
r CD f H = s f a f 2 1 a f R f , r AB f H = s f a f = r f .
It is convenient to introduce a polar coordinate system ( r , θ ) with the origin at x = H s f / a f and y = 0 (point O 1 in Figure 1). The coordinate curves of this coordinate system coincide with the coordinate curves of the ( ζ , η ) -coordinate system. Therefore, σ ζ = σ r and σ η = σ θ where σ r and σ θ are the normal stresses in the polar coordinate system. Moreover, r = R f H at ζ = 0 and r = r f H at ζ = 1 . The equilibrium equation for the increment of the stresses, Δ σ ζ and Δ σ η , in the polar coordinate system can be written as
Δ σ ζ ρ = Δ σ η Δ σ ζ ρ ,
where ρ = r / H . Since σ ζ = 0 at ζ = 0 and ζ = 1 at any stage of the process, the increment of this stress should satisfy the conditions
Δ σ ζ = 0 ,
for ζ = 0 and ζ = 1 .
The displacement components from the configuration corresponding to the end of loading in the polar coordinate system are supposed to be
u r = H U 0 R f 2 ρ ρ V 0 2 and u θ = H ρ θ V 0 ,
where U 0 and V 0 are dimensionless constants. Using Equation (41), the increment of the normal strains in the polar coordinate system is determined as
Δ ε r = V 0 2 U 0 R f 2 ρ 2 , Δ ε θ = V 0 2 + U 0 R f 2 ρ 2 .
The increment of the deviatoric stresses is found from Equation (42) and the Hooke’s law (Equation (9)) where the stresses and strains should be replaced with the corresponding increments. Then,
Δ τ r = G 0 g ζ V 0 + 2 U 0 R f 2 ρ 2 , Δ τ θ = G 0 g ζ V 0 + 2 U 0 R f 2 ρ 2 .
Using this solution, the right hand side of Equation (36) can be rewritten as
Δ σ η Δ σ ζ ρ = Δ σ θ Δ σ r ρ = Δ τ θ Δ τ r ρ = 2 G 0 g ζ V 0 + 2 U 0 R f 2 ρ 2 .
The Lagrangian coordinate ζ at the end of loading is expressed in terms of ρ as [18]
ζ = ρ 2 a f s f a f .
Using this equation, it is possible to eliminate ζ in Equation (44). Then, substituting Equation (44) into Equation (39) and integrating gives
Δ σ ζ σ 0 = 2 3 k r f ρ g ζ χ V 0 + 2 U 0 R f 2 χ 2 d χ .
It is evident that this solution satisfies the boundary condition in Equation (40) at ζ = 1 (or ρ = r f ). The other boundary conditions in Equations (40) and (46) combine to yield
V 0 r f R f g ζ ρ d ρ + 2 U 0 R f 2 r f R f g ζ ρ 3 d ρ = 0 .
Solving this equation for V 0 results in
V 0 = 2 U 0 R f 2 r f R f g ζ ρ 3 d ρ r f R f g ζ ρ d ρ 1 .
Using Equations (43) and (46), it is possible to represent the distribution of Δ σ η as
Δ σ η σ 0 = Δ σ ζ σ 0 2 Δ τ r σ 0 = 2 3 k r f ρ g ζ χ V 0 + 2 U 0 R f 2 χ 2 d χ + 2 3 k g ζ V 0 + 2 U 0 R f 2 ρ 2 .
The constant V 0 can be eliminated in Equations (46) and (49) by means of Equation (48). It is then obvious that both Δ σ ζ and Δ σ η are proportional to U 0 . The distribution of the residual stresses follows from Equations (46) and (49) in the form
σ ζ r e s σ 0 = σ ζ f σ 0 + 2 3 k r f ρ g ζ μ V 0 + 2 U 0 R f 2 χ 2 d χ , σ η r e s σ 0 = σ η f σ 0 + 2 3 k r f ρ g ζ χ V 0 + 2 U 0 R f 2 χ 2 d χ + 2 3 k g ζ V 0 + 2 U 0 R f 2 ρ 2 .
As before, ζ should be eliminated by means of Equation (45) and V 0 by means of Equation (48). The constant U 0 remains to be found. To this end, it is necessary to use the condition that the bending moment vanishes at the end of unloading. Using Equations (11) and (45), this condition can be represented as
R f r f σ η r e s σ 0 ρ d ρ = 0 .
This equation should be solved for U 0 numerically. Then, Equation (50) supplies the distribution of the residual stresses. To verify that the solution given in Section 4 is valid, this distribution should be substituted into the yield criterion in Equation (7) where σ ζ and σ η should be replaced with σ ζ r e s and σ η r e s , respectively. The left-hand side of Equation (7) should be less than or equal to 2 / 3 σ 0 β ( ζ ) in the range 1 ζ 0 .

5. Numerical Examples

Several numerical examples are presented in this section, based on the analytical solutions developed in the previous sections. Our chosen modulus gradient function is g ( ζ ) = 1 + ( G 1 / G 0 1 ) ( ζ ) N , and yield stress gradient function β ( ζ ) = 1 + ( σ 1 / σ 0 1 ) ( ζ ) N . The power law exponent N controls the functional distribution of material properties along the thickness coordinate ζ . The power law distributions in modulus and yield stress with the same N have been proposed in the literature [26,27]. The material parameters used in our numerical calculations are listed in Table 1.

5.1. Homogeneous Sheet under Bending

When the homogenous sheet is under bending, both edges will simultaneously develop plastic zones. Figure 3a shows the movement of the two elastic-plastic boundaries toward the centerline of the sheet, as deformation magnitude increases. The deformation magnitude is measured by parameter a. The applied bending moment is a function of a, as shown in Figure 3b. As an illustration of the developed analytical solutions in previous sections, Figure 4 shows the stress distributions along the sheet under two different deformation magnitudes. As can be seen, the plastic zones increase with a for σ η , while σ ζ remains in elastic regime. The reason for σ η is not perfectly horizontal in the plastic zone is because our numerical codes do not allow N set equal to zero, hence a very small N is chosen, as shown in Table 1. After unloading, Figure 5a,b shows the residual stress distributions under two different as. Larger a increases the magnitude of residual stresses after unloading. Moreover, the residual stress σ ζ is zero at the left and right edges, as indicated by the red short-dashed line.

5.2. FGM Sheet Belonged to Case (i) under Bending

In Case (i), the left edge ( ζ = 1 ) of the sheet has smaller yield stress, hence a plastic zone will start on the left edge first. Figure 6 shows the relationship between the applied bending moment and deformation magnitude a with the gradient function exponent N = 1 and 3. As can be seen, larger m is required for N = 3 than that for N = 1 , as a increases. Under given deformation magnitudes, Figure 7 shows the stress distributions in the Case (i) FGM under plastic deformation. Larger plastic zone is developed at the left edge as deformation increases. After unloading, residual stress distributions are shown in Figure 8 and Figure 9 for the N = 1 and N = 3 FGM, respectively. Residual stresses are more predominant at the left edge. In addition, the residual stress σ ζ is zero at the left and right edges, as indicated by the red short-dashed line.

5.3. FGM Sheet Belonged to Case (ii) under Bending

In Case (ii), the right edge ( ζ = 0 ) of the sheet has smaller yield stress, hence a plastic zone will start on the right edge first. Figure 10 shows the relationship between the applied bending moment and deformation magnitude a with the gradient function exponent N = 1 and 3. As can be seen, larger m is required for N = 1 than that for N = 3 , as a increases. Under given deformation magnitudes, Figure 11 shows the stress distributions in the Case (ii) FGM under plastic deformation. Larger plastic zone is developed at the right edge as deformation increases. After unloading, residual stress distributions are shown in Figure 12 and Figure 13 for the N = 1 and N = 3 FGM, respectively. Residual stresses are more predominant at the right edge. The residual stress σ ζ is zero at the left and right edges, as indicated by the red short-dashed line. The results from the illustrative examples solved here may serve as benchmark solutions for data obtained from numerical or experimental methods.

6. Conclusions

Efficient analytical and numerical methods and procedures have been developed and programmed to predict the distribution of stresses in a sheet of incompressible material subject to plane strain pure bending at large strain and then the distribution of residual stresses after unloading. Springback is also predicted. It has been assumed that the sheet is made of functionally graded material. The general theory has been developed for an arbitrary through thickness distribution of material properties assuming that the initiation of plastic yielding occurs at one of the surfaces of the sheet. This assumption can be verified using the purely elastic solution (Equation (15)) and the yield criterion (Equation (7)). It is possible to use the general solutions (Equations (15), (18) and (25)) even if the assumption is not satisfied but constants of integration should be added. Then, these general solutions should be combined to satisfy the boundary conditions and the conditions at elastic/plastic boundaries. An illustrative example is concentrated on power law distributions of material properties. Using the numerical code developed in this work enables the effect of parameters involved in these laws to be predicted effectively. The calculated examples show the analytical solutions derived here can systematically treat the plastic problems of the homogenous or functionally graded sheet. The magnitudes of applied moment may be strongly influenced by the power law exponent as deformation increases, which provides an effective way to design the functionally graded sheets.
The method employed to derive the solution in this paper can be extended to cyclic loading. This new solution may be useful for the interpretation of experimental data from the reverse bending fatigue test (for example, [28]).

Author Contributions

Conceptualization, S.A., Y.-C.W. and L.L.; methodology, S.A. and Y.-C.W.; software, S.A. and Y.-C.W.; validation, S.A., Y.-C.W. and L.L.; formal analysis, S.A. and Y.-C.W.; investigation, S.A., Y.-C.W. and L.L.; resources, S.A., Y.-C.W. and L.L.; data curation, S.A. and Y.-C.W.; writing–original draft preparation, S.A., Y.-C.W. and L.L.; writing–review and editing, S.A., Y.-C.W. and L.L.; visualization, S.A. and Y.-C.W.; supervision, S.A., Y.-C.W. and L.L.; project administration, S.A., Y.-C.W. and L.L.; funding acquisition, S.A. and Y.-C.W.

Funding

S.A. acknowledges support from the Russian Foundation for Basic Research (Project 17-58-560005). This research was funded, in part, by Ministry of Science and Technology, Taiwan, grant number MOST 107-2221-E-006-028. This research was, in part, supported by the Ministry of Education, Taiwan, and the Aim for the Top University Project to the National Cheng Kung University (NCKU).

Acknowledgments

This work was initiated while Y.-C.W. was a visiting researcher at Beihang University, Beijing, China. Y.-C.W. is grateful to the National Center for High-performance Computing, Taiwan, for computer time and facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Birman, V.; Byrd, L.W. Modeling and analysis of functionally graded materials and structures. Appl. Mech. Rev. 2007, 60, 195–216. [Google Scholar] [CrossRef]
  2. Tutuncu, N.; Ozturk, M. Exact solutions for stresses in functionally graded pressure vessels. Compos. Part B Eng. 2001, 32, 683–686. [Google Scholar] [CrossRef]
  3. Akis, T. Elastoplastic analysis of functionally graded spherical pressure vessels. Comp. Mater. Sci. 2009, 46, 545–554. [Google Scholar] [CrossRef]
  4. Sadeghian, M.; Toussi, H.E. Axisymmetric yielding of functionally graded spherical vessel under thermos-mechanical loading. Comp. Mater. Sci. 2011, 50, 975–981. [Google Scholar] [CrossRef]
  5. Li, X.-Y.; Li, P.-D.; Kang, G.-Z. Axisymmetric thermo-elasticity field in a functionally graded circular plate of transversely isotropic material. Math. Mech. Solids 2013, 18, 464–475. [Google Scholar] [CrossRef]
  6. Horgan, C.O.; Chan, A.M. The pressurized hollow cylinder or disk problem for functionally graded isotropic linearly elastic materials. J. Elast. 1999, 55, 43–59. [Google Scholar] [CrossRef]
  7. Tutuncu, N. Stresses in thick-walled FGM cylinders with exponentially-varying properties. Eng. Struct. 2007, 29, 2032–2035. [Google Scholar] [CrossRef]
  8. You, L.H.; Wang, J.X.; Tang, B.P. Deformations and stresses in annular disks made of functionally graded materials subjected to internal and/or external pressure. Meccanica 2009, 44, 283–292. [Google Scholar] [CrossRef]
  9. Jabbari, M.; Sohrabpour, S.; Eslami, M.R. Mechanical and thermal stresses in a functionally graded hollow cylinder due to radially symmetric loads. Int. J. Press. Vessels Pip. 2002, 79, 493–497. [Google Scholar] [CrossRef]
  10. Eraslan, A.N.; Akis, T. Plane strain analytical solutions for a functionally graded elastic–plastic pressurized tube. Int. J. Press. Vessels Pip. 2006, 83, 635–644. [Google Scholar] [CrossRef]
  11. Horgan, C.O.; Chan, A.M. The stress response of functionally graded isotropic linearly elastic rotating disks. J. Elast. 1999, 55, 219–230. [Google Scholar] [CrossRef]
  12. Peng, X.-L.; Li, X.-F. Effects of gradient on stress distribution in rotating functionally graded solid disks. J. Mech. Sci. Technol. 2012, 26, 1483–1492. [Google Scholar] [CrossRef]
  13. Bayat, M.; Saleem, M.; Sahari, B.B.; Hamouda, A.M.S.; Mahdi, E. Analysis of functionally graded rotating disks with variable thickness. Mech. Res. Commun. 2008, 35, 283–309. [Google Scholar] [CrossRef] [Green Version]
  14. Peng, X.-L.; Li, X.-F. Elastic analysis of rotating functionally graded polar orthotropic disks. Int. J. Mech. Sci. 2012, 60, 84–91. [Google Scholar] [CrossRef]
  15. Dai, H.-L.; Dai, T.; Zheng, H.-Y. Stresses distributions in a rotating functionally graded piezoelectric hollow cylinder. Meccanica 2012, 47, 423–436. [Google Scholar] [CrossRef]
  16. Callioglu, H.; Sayer, M.; Demir, E. Elastic–plastic stress analysis of rotating functionally graded discs. Thin-Walled Struct. 2015, 94, 38–44. [Google Scholar] [CrossRef]
  17. Elishakoff, I.; Pentaras, D.; Gentilini, C. Mechanics of Functionally Graded Material Structures; World Scientific: Singapore, 2016; ISBN 978-981-4656-58-0. [Google Scholar]
  18. Alexandrov, S.; Kim, J.-H.; Chung, K.; Kang, T.-J. An alternative approach to analysis of plane-strain pure bending at large strains. J. Strain Anal. Eng. Des. 2006, 41, 397–410. [Google Scholar] [CrossRef]
  19. Hill, R. The Mathematical Theory of Plasticity; Clarendon Press: Oxford, UK, 1950; ISBN 0-19-850367-9. [Google Scholar]
  20. Alexandrov, S.; Hwang, Y.-M. The bending moment and springback in pure bending of anisotropic sheets. Int. J. Solids Struct. 2009, 46, 4361–4368. [Google Scholar] [CrossRef] [Green Version]
  21. Alexandrov, S.; Hwang, Y.-M. Plane strain bending with isotropic strain hardening at large strains. Trans. ASME J. Appl. Mech. 2010, 77, 064502. [Google Scholar] [CrossRef]
  22. Alexandrov, S.; Gelin, J.-C. Plane strain pure bending of sheets with damage evolution at large strains. Int. J. Solids Struct. 2011, 48, 1637–1643. [Google Scholar] [CrossRef] [Green Version]
  23. Alexandrov, S.; Hwang, Y.-M. Influence of Bauschinger effect on springback and residual stresses in plane strain pure bending. Acta Mech. 2011, 220, 47–59. [Google Scholar] [CrossRef]
  24. Verguts, H.; Sowerby, R. The pure plastic bending of laminated sheet metals. Int. J. Mech. Sci. 1975, 17, 31–51. [Google Scholar] [CrossRef]
  25. Romano, G.; Barretta, R.; Diaco, M. Geometric continuum mechanics. Meccanica 2014, 49, 111–133. [Google Scholar] [CrossRef]
  26. Nakumura, T.; Wang, T.; Sampath, S. Determination of properties of graded materials by inverse analysis and instrumented indentation. Acta Mater. 2000, 48, 4293–4306. [Google Scholar] [CrossRef]
  27. Huang, H.; Chen, B.; Han, Q. Investigation on buckling behaviors of elastoplastic functionally graded cylindrical shells subjected to torsional loads. Comput. Struct. 2014, 118, 234–240. [Google Scholar] [CrossRef]
  28. Faghidian, S.A.; Jozie, A.; Sheykhloo, M.J.; Shamsi, A. A novel method for analysis of fatigue life measurements based on modified Shepard method. Int. J. Fatigue 2014, 68, 144–149. [Google Scholar] [CrossRef]
Figure 1. Geometric configuration of the bending problem: (a) before deformation; and (b) after deformation.
Figure 1. Geometric configuration of the bending problem: (a) before deformation; and (b) after deformation.
Materials 12 00456 g001
Figure 2. Schematics of elastic and plastic zones.
Figure 2. Schematics of elastic and plastic zones.
Materials 12 00456 g002
Figure 3. For the homogeneous sheet under bending: (a) ζ 1 or ζ 2 vs. a; and (b) applied bending moment m vs. a.
Figure 3. For the homogeneous sheet under bending: (a) ζ 1 or ζ 2 vs. a; and (b) applied bending moment m vs. a.
Materials 12 00456 g003
Figure 4. Stress distributions for the bent homogenous sheet with two different loading conditions a = a 1 h = 0.0121 and a = a 2 h = 0.0775 . Two plastic zones, one developed from the left edge and the other from the right edge, increase their size as a increases for σ η .
Figure 4. Stress distributions for the bent homogenous sheet with two different loading conditions a = a 1 h = 0.0121 and a = a 2 h = 0.0775 . Two plastic zones, one developed from the left edge and the other from the right edge, increase their size as a increases for σ η .
Materials 12 00456 g004
Figure 5. Residual stress distribution along the homogeneous sheet for: (a) a = a r 1 h = 0.0096 ; and (b) a = a r 1 h = 0.0295 . Solid line is for residual σ η and dashed line for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Figure 5. Residual stress distribution along the homogeneous sheet for: (a) a = a r 1 h = 0.0096 ; and (b) a = a r 1 h = 0.0295 . Solid line is for residual σ η and dashed line for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Materials 12 00456 g005
Figure 6. Bending moment m vs. a for Case (i) with N = 1 and N = 3 . With sufficiently large deformation, i.e., large a, Case (iii) is automatically developed, hence both edges are plastically deformed.
Figure 6. Bending moment m vs. a for Case (i) with N = 1 and N = 3 . With sufficiently large deformation, i.e., large a, Case (iii) is automatically developed, hence both edges are plastically deformed.
Materials 12 00456 g006
Figure 7. Stress distributions for Case (i) with N = 1 with two different deformation magnitudes, i.e., two different as, for: (a) N = 1 ; and (b) N = 3 . Two plastic zones, one developed from the left edge and the other from the right edge, increase their size as a increases for σ η .
Figure 7. Stress distributions for Case (i) with N = 1 with two different deformation magnitudes, i.e., two different as, for: (a) N = 1 ; and (b) N = 3 . Two plastic zones, one developed from the left edge and the other from the right edge, increase their size as a increases for σ η .
Materials 12 00456 g007
Figure 8. Residual stress distributions for Case (i) with N = 1 under deformation: (a) a = a r 1 ( i ) = 0.011785065 ; and (b) a = a r 1 ( i ) = 0.074421129 . Solid line is for residual σ η and dashed line for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Figure 8. Residual stress distributions for Case (i) with N = 1 under deformation: (a) a = a r 1 ( i ) = 0.011785065 ; and (b) a = a r 1 ( i ) = 0.074421129 . Solid line is for residual σ η and dashed line for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Materials 12 00456 g008
Figure 9. Residual stress distributions for Case (i) with N = 3 under deformation: (a) a = a r 3 ( i ) = 0.011023776 ; and (b) a = a r 4 ( i ) = 0.051782217 . Solid line is for residual σ η and dashed line for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Figure 9. Residual stress distributions for Case (i) with N = 3 under deformation: (a) a = a r 3 ( i ) = 0.011023776 ; and (b) a = a r 4 ( i ) = 0.051782217 . Solid line is for residual σ η and dashed line for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Materials 12 00456 g009
Figure 10. Bending moment m vs. a for Case (ii) with N = 1 and N = 3 . With sufficiently large deformation, i.e., large a, Case (iii) is developed, hence both edges are plastically deformed.
Figure 10. Bending moment m vs. a for Case (ii) with N = 1 and N = 3 . With sufficiently large deformation, i.e., large a, Case (iii) is developed, hence both edges are plastically deformed.
Materials 12 00456 g010
Figure 11. Stress distributions for Case (ii) with two different deformation magnitudes, i.e., two different as, for: (a) N = 1 ; and (b) N = 3 . Two plastic zones, one developed from the left edge and the other from the right edge, increase their size as a increases for σ η .
Figure 11. Stress distributions for Case (ii) with two different deformation magnitudes, i.e., two different as, for: (a) N = 1 ; and (b) N = 3 . Two plastic zones, one developed from the left edge and the other from the right edge, increase their size as a increases for σ η .
Materials 12 00456 g011
Figure 12. Residual stress distributions for Case (ii) with N = 1 under deformation: (a) a = a r 1 ( i i ) = 0.004915193 ; and (b) a = a r 1 ( i i ) = 0.019238382 . Solid curve is for residual σ η and dashed curve for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Figure 12. Residual stress distributions for Case (ii) with N = 1 under deformation: (a) a = a r 1 ( i i ) = 0.004915193 ; and (b) a = a r 1 ( i i ) = 0.019238382 . Solid curve is for residual σ η and dashed curve for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Materials 12 00456 g012
Figure 13. Residual stress distributions for Case (ii) with N = 3 under deformation: (a) a = a r 3 ( i i ) = 0.005565273 ; and (b) a = a r 4 ( i i ) = 0.016428139 . Solid curve is for residual σ η and dashed curve for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Figure 13. Residual stress distributions for Case (ii) with N = 3 under deformation: (a) a = a r 3 ( i i ) = 0.005565273 ; and (b) a = a r 4 ( i i ) = 0.016428139 . Solid curve is for residual σ η and dashed curve for residual σ ζ . Zeros of σ η r e s and σ ζ r e s are indicated by orange solid line and red dashed line, respectively.
Materials 12 00456 g013
Table 1. Material parameters used in the numerical examples.
Table 1. Material parameters used in the numerical examples.
G 0 , GPa G 1 , GPa σ 0 , GPa σ 1 , GPaN
Homogeneous3030110.0001
FGM Case (i)301010.11
FGM Case (i)301010.13
FGM Case (ii)10300.111
FGM Case (ii)10300.113

Share and Cite

MDPI and ACS Style

Alexandrov, S.; Wang, Y.-C.; Lang, L. A Theory of Elastic/Plastic Plane Strain Pure Bending of FGM Sheets at Large Strain. Materials 2019, 12, 456. https://doi.org/10.3390/ma12030456

AMA Style

Alexandrov S, Wang Y-C, Lang L. A Theory of Elastic/Plastic Plane Strain Pure Bending of FGM Sheets at Large Strain. Materials. 2019; 12(3):456. https://doi.org/10.3390/ma12030456

Chicago/Turabian Style

Alexandrov, Sergey, Yun-Che Wang, and Lihui Lang. 2019. "A Theory of Elastic/Plastic Plane Strain Pure Bending of FGM Sheets at Large Strain" Materials 12, no. 3: 456. https://doi.org/10.3390/ma12030456

APA Style

Alexandrov, S., Wang, Y. -C., & Lang, L. (2019). A Theory of Elastic/Plastic Plane Strain Pure Bending of FGM Sheets at Large Strain. Materials, 12(3), 456. https://doi.org/10.3390/ma12030456

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