Next Article in Journal
Performance of Nano-Silica Modified Self-Compacting Glass Mortar at Normal and Elevated Temperatures
Next Article in Special Issue
A Theory of Elastic/Plastic Plane Strain Pure Bending of FGM Sheets at Large Strain
Previous Article in Journal
Investigation of Self-Assembly and Charge-Transport Property of One-dimensional PDI8-CN2 Nanowires by Solvent-Vapor Annealing
Previous Article in Special Issue
Stress Concentration and Damage Factor Due to Central Elliptical Hole in Functionally Graded Panels Subjected to Uniform Tensile Traction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Description of Residual Stress and Strain Fields in FGM Hollow Disc Subject to External Pressure

1
Department of Civil Engineering, Peoples’ Friendship University of Russia (RUDN University), Miklukho-Maklaya st. 6, 117198 Moscow, Russia
2
School of Mechanical Engineering and Automation, Beihang University, No. 37 Xueyuan Road, Beijing 100191, China
3
Ishlinsky Institute for Problems in Mechanics, 101-1 Prospect Vernadskogo, 119526 Moscow, Russia
*
Author to whom correspondence should be addressed.
Materials 2019, 12(3), 440; https://doi.org/10.3390/ma12030440
Submission received: 15 December 2018 / Revised: 21 January 2019 / Accepted: 25 January 2019 / Published: 31 January 2019

Abstract

:
Elastic/plastic stress and strain fields are obtained in a functionally graded annular disc of constant thickness subject to external pressure, followed by unloading. The elastic modulus and tensile yield stress of the disc are assumed to vary along the radius whereas the Poisson’s ratio is kept constant. The flow theory of plasticity is employed. However, it is shown that the equations of the associated flow rule, which are originally written in terms of plastic strain rate, can be integrated with respect to the time giving the corresponding equations in terms of plastic strain. This feature of the solution significantly facilitates the solution. The general solution is given for arbitrary variations of the elastic modulus and tensile yield stress along the radial coordinate. However, it is assumed that plastic yielding is initiated at the inner radius of the disc and that no other plastic region appears in the course of deformation. The solution in the plastic region at loading reduces to two ordinary differential equations. These equations are solved one by one. Unloading is assumed to be purely elastic. This assumption should be verified a posteriori. An illustrative example demonstrates the effect of the variation of the elastic modulus and tensile yield stress along the radius on the distribution of stresses and strains at the end of loading and after unloading. In this case, it is assumed that the material properties vary according to power-law functions.

1. Introduction

Stress and strain analyses of solid and hollow circular discs have long been an important topic in the mechanics of solids. The motivation of doing such analyses is that circular discs subject to mechanical, thermal, and inertial loading are used in many sectors of industry. The performance of such discs under service conditions can be improved by using functionally graded materials. The material may be continuously graded or be piecewise homogeneous. It is assumed in the present paper that the distribution of all material properties is axisymmetric. Discs made of homogeneous materials are not discussed. Discs made of functionally graded materials have been the subject of intense research. A linearly elastic solution under plane stress and plane strain conditions has been given in [1], assuming that the disc is loaded by external or internal pressure. It has been concluded that the stress response of the functionally graded disc is significantly different from that of the homogeneous disc. Another plane stress solution of this boundary value problem has been obtained in [2] and another plane strain solution in [3]. A thermoelastic stress solution for a disc of variable thickness has been found in [2]. A thermoelastic analysis of a disc subject to a steady-state temperature distribution together with external and internal pressures has been provided in [4]. It has been assumed in this work that the material properties are arbitrary smooth functions of the radial coordinate. A similar boundary value problem for a multilayered hollow cylinder has been solved in [5]. A design driven by the minimization of induced stresses in elastic multilayer cylinders under plane stress conditions has been proposed in [6]. All of the solutions above are purely elastic or thermoelastic. The process of autofrettage of a functionally graded cylinder has been studied in [7]. The analysis of this process requires the use of an elastic/plastic model. In [7], the deformation theory of plasticity together with the von Mises yield criterion has been employed. Another elastic/plastic plane strain solution for a functionally graded cylinder has been given in [8]. The solution is based on Tresca’s yield criterion, which significantly simplifies the analysis even in the case of the flow theory of plasticity.
There is a vast amount of literature on functionally graded rotating discs. The elastic response of an arbitrary functionally graded polar orthotropic disc has been investigated in [9]. Another purely elastic solution has been given in [10], using the finite difference method. Thermoelastic analyses have been presented in [11,12,13,14]. The effect of a non-uniform heat source on thermoelastic behavior of a functionally graded rotating disc has been investigated in [15]. The effect of viscosity on the response of a functionally graded rotating disc of variable thickness has been studied in [16]. The limit of elastic angular velocity has been determined in [17]. The effect of variable angular velocity on the elastic response of a functionally graded rotating disc has been analyzed in [18]. A design driven by weight optimization of a disc subject to thermomechanical loading has been proposed in [19]. Most of the available elastic/plastic solutions fall into three categories. A series of solutions is devoted to discs obeying Tresca’s yield criteria [20,21]. As it has been mentioned before, the use of Tresca’s yield criterion significantly simplifies the solution. Another category includes the solutions for the deformation theory of plasticity [22,23,24,25]. In some cases, using deformation theories of plasticity is justified since the stress path is nearly proportional. However, it has been shown in [26] that it may not be so in thin discs. The third category includes stress solutions [27,28]. In this case, no flow rule is necessary to find the solution.
An advantage of the present elastic/plastic solution is that the flow theory of plasticity in conjunction with the von Mises yield criterion is employed. It is assumed that a hollow disc is subject to external pressure, followed by unloading. First, the general solution is derived under plane stress conditions assuming that the elastic modulus and tensile yield stress are arbitrary smooth functions of the radial coordinate. It is, however, assumed that plastic yielding initiates at the inner radius of the disc and that there is one plastic region throughout the process of deformation. The Poisson’s ratio is supposed to be constant. This is a typical assumption for functionally graded discs [1,3,10,19]. Second, a numerical example is given assuming that the material properties vary according to power-law functions. This is also a typical assumption for functionally graded discs [1,9,10,11,12,13]. The solution found can be considered as an extension of the solution provided in [1] to the plastic range.

2. Statement of the Problem

Consider a thin hollow disc of functionally graded material subject to uniform pressure p 0 over its outer radius b 0 , followed by unloading. The inner radius of the disc is denoted as a 0 . The thickness of the disc is constant. The mechanical properties of the disc are classified in terms of the yield stress tension σ Y , Poisson’s ratio ν , and Young’s modulus E. It is assumed that the Poisson’s ratio is constant, whereas the value of both σ Y and E vary with radius. It is convenient to use a cylindrical coordinate system ( r ,   θ ,   z ) whose z-axis coincides with the axis of symmetry of the disc. The normal stresses in this coordinate system are the principal stresses. The state of stress is plane (i.e., the axial stress in the cylindrical coordinate system vanishes). Therefore, Hooke’s law can be written as
ε r e = σ r ν σ θ E ,   ε θ e = σ θ ν σ r E ,   ε z e = ν ( σ r + σ θ ) E
here σ r is the radial stress, σ θ is the circumferential stress, ε r e , ε θ e , and ε z e are the elastic strains referred to the cylindrical coordinate system. Plastic yielding is controlled by the von Mises yield criterion. Under a plane stress condition, this criterion reads
σ r 2 + σ θ 2 σ r σ θ = σ Y 2
The flow theory of plasticity is adopted. The flow rule associated with the yield criterion (2) is
ξ r p = λ ( 2 σ r σ θ ) ,   ξ θ p = λ ( 2 σ θ σ r ) ,   ξ z p = λ ( σ r + σ θ )
here ξ r p , ξ θ p , and ξ z p are the plastic strain rates referred to the cylindrical coordinate system and λ is a non-negative multiplier. The total strain components in the cylindrical coordinate system are given by
ε r = ε r e + ε r p ,   ε θ = ε θ e + ε θ p ,   ε z = ε z e + ε z p
here ε r p , ε θ p , and ε z p are the plastic strains referred to the cylindrical coordinate system. The constitutive equations should be complemented with the equilibrium equation
σ r r + σ r σ θ r = 0
and the equation of strain compatibility of the form
r ε θ r = ε r ε θ
The boundary conditions at the stage of loading are
σ r = p 0
for r = b 0 and
σ r = 0
for r = a 0 . The boundary conditions at the stage of unloading will be formulated in Section 6.
It is convenient to introduce the following dimensionless quantities:
ρ = r b 0 ,   a = a 0 b 0 ,   k = σ 0 E 0 ,   p = p 0 σ 0
here σ 0 is the value of σ Y at r = b 0 and E 0 is the value of E at r = b 0 . Then, the variation of σ Y and E with ρ can be represented as
σ Y = σ 0 Φ ( ρ )   and   E = E 0 η ( ρ )
where Φ ( ρ ) and η ( ρ ) are arbitrary functions of ρ satisfying the conditions Φ ( ρ ) = 1 and η ( ρ ) = 1 at ρ = 1 . In what follows, it is assumed that these functions are such that plastic yielding initiates at the inner radius of the disc and no other plastic region appears in continued deformation.

3. Purely Elastic Solution

When p is small enough, the entire disc is elastic. In this case, the total strains are equal to the elastic strains. The system of equations comprises Hooke’ law, the equilibrium equation, and the equation of strain compatibility. Using (9) it is possible to rewrite Equations (5) and (6) as
σ r ρ + σ r σ θ ρ = 0
and
ρ ε θ e ρ = ε r e ε θ e
Eliminating the strains in this equation by means of (1) and using (9) results in
σ θ ρ + ν ρ ( σ r σ θ ) + η ρ ( σ θ ν σ r ) ( ρ / η ) ρ ( σ r ν σ θ ) ρ = 0
Equations (11) and (13) comprise the system for determining the distribution of stresses in the purely elastic disc. Then, the distribution of strains can be readily found from (1) and (9). However, the purely elastic solution is not of interest in the case under consideration. Therefore, the solution to Equations (11) and (13) is only necessary to determine the value of p at which plastic yielding is initiated. This value of p is denoted as p e . By assumption, plastic yielding is initiated at ρ = a . It follows from the boundary condition (8) and the yield criterion (2) that σ θ = σ Y at ρ = a at the initiation of plastic yielding. Using (10), this condition can be rewritten as σ θ = σ 0 Φ ( a ) at ρ = a . This is one of the boundary conditions of the boundary value problem to be solved. The other boundary condition is given by (8). Equations (11) and (13) should be solved together with these boundary conditions. The value of p e is readily found from this solution as p e = σ r ( 1 ) / σ 0 . In what follows, it is assumed that p > p e .

4. Elastic/Plastic Stress Solution

If p > p e , then the disc consists of two regions, elastic and plastic. The elastic region occupies the domain ρ c ρ 1 and the plastic region the domain a ρ ρ c . Here, ρ c is the elastic/plastic boundary. Consider the plastic region. The yield criterion (2) is satisfied by the following substitution:
σ r σ 0 = 2 Φ ( ρ ) sin ψ 3   and   σ θ σ 0 = Φ ( ρ ) ( sin ψ 3 + cos ψ )
here ψ is a new function of ρ . This function should be found from the solution. Substituting (14) into (11) yields
d ψ d ρ + tan ψ Φ d Φ d ρ + ( tan ψ 3 ) 2 ρ = 0
The boundary condition to this equation follows from (8) and (14). In particular, σ r = 0 if ψ = 0 or ψ = π . It is evident that σ θ < 0 at ρ = a . Then, the boundary condition to Equation (15) is
ψ = π
for ρ = a . Solving Equation (15) together with this boundary condition supplies the variation of ψ with ρ . This solution and (14) determine the distribution of the stresses in the plastic region. Let p p be the value of p at which the entire disc becomes plastic. Putting in the solution for the radial stress ρ = 1 gives the value of p p as p p = σ r ( 1 ) / σ 0 . In what follows, it is assumed that p < p p . In this case, a < ρ c < 1 . The value of σ r on the plastic side of the elastic/plastic boundary is denoted as σ r c and the value of σ θ on the plastic side of the elastic/plastic boundary as σ θ c . These values are readily found from the solution of (15) and (14). Equations (11) and (13) are valid in the elastic region. The radial and circumferential stresses must be continuous across the elastic/plastic boundary. Therefore, the boundary conditions to Equations (11) and (13) are
σ r = σ r c   and   σ θ = σ θ c
for ρ = ρ c . Solving Equations (11) and (13) together with these boundary conditions supplies the distribution of the radial and circumferential stresses in the elastic region. In particular, the value of p involved in (7) is determined from the equation p = σ r ( 1 ) / σ 0 . Therefore, the solution found connects p and ρ c . One of these parameters should be prescribed. Then, the other parameter is determined from the solution.

5. Elastic/Plastic Strain Solution

Consider the plastic region, a ρ ρ c . Eliminating λ between the equations in (3) gives
ξ r p ξ θ p = ( 2 σ r σ θ ) ( 2 σ θ σ r ) ,   ξ z p ξ θ p = ( σ r + σ θ ) ( 2 σ θ σ r )
Using (14), the stresses in these equations can be expressed in terms of ψ . Then, taking into account that ξ r p = ε r p / t , ξ θ p = ε θ p / t , and ξ z p = ε z p / t Equation (18) is transformed to
ε r p t = ( 3 tan ψ 1 ) 2 ε θ p t ,   ε z p t = ( 3 tan ψ + 1 ) 2 ε θ p t
here t is the time. It is seen from the structure of Equation (15) and the boundary condition (16) that ψ is independent of t. Therefore, the coefficients of ε θ p / t in (19) are also independent of t, and the equations in (19) can be immediately integrated with respect to the time to give
ε r p = ( 3 tan ψ 1 ) 2 ε θ p ,   ε z p = ( 3 tan ψ + 1 ) 2 ε θ p
It has been taken into account here that ε r p = ε θ p = ε z p = 0 at the elastic/plastic boundary. The elastic strains in the plastic region, ε r e p , ε θ e p and ε z e p , are determined from (1) and (14) with the use of (9) and (10). As a result,
ε r e p = k Λ [ ( 2 ν ) 3 sin ψ ν cos ψ ] ,   ε θ e p = k Λ [ ( 1 2 ν ) 3 sin ψ + cos ψ ] ,   ε z e p = ν k Λ ( 3 sin ψ + cos ψ )
here Λ = Φ / η . The total strains are found from (4), (20), and (21) as
ε r = k Λ [ ( 2 ν ) 3 sin ψ ν cos ψ ] + ( 3 tan ψ 1 ) 2 ε θ p , ε θ = k Λ [ ( 1 2 ν ) 3 sin ψ + cos ψ ] + ε θ p , ε z = ν k Λ ( 3 sin ψ + cos ψ ) ( 3 tan ψ + 1 ) 2 ε θ p
It follows from these equations that
ε r ε θ = ( tan ψ 3 ) [ k ( 1 + ν ) cos ψ 3 Λ + 3 2 ε θ p ] , ε θ ρ = k d Λ d ρ [ ( 1 2 ν ) 3 sin ψ + cos ψ ] + k Λ [ ( 1 2 ν ) 3 cos ψ sin ψ ] d ψ d ρ + ε θ p ρ .
.
These equations and Equation (6), in which r should be replaced with ρ by means of (9), combine to give
ε θ p ρ 3 2 ( tan ψ 3 ) ρ ε θ p + k d Λ d ρ [ ( 1 2 ν ) 3 sin ψ + cos ψ ] + k Λ [ ( 1 2 ν ) 3 cos ψ sin ψ ] d ψ d ρ k ( 1 + ν ) Λ ( sin ψ 3 cos ψ ) 3 ρ = 0 .
The derivative d ψ / d ρ in this equation can be eliminated by means of (15). Then,
ε θ p ρ 3 2 ( tan ψ 3 ) ρ ε θ p + k d Λ d ρ [ ( 1 2 ν ) 3 sin ψ + cos ψ ] k Λ [ ( 1 2 ν ) 3 cos ψ sin ψ ] [ tan ψ Φ d Φ d ρ + ( tan ψ 3 ) 2 ρ ] k ( 1 + ν ) Λ ( sin ψ 3 cos ψ ) 3 ρ = 0
Since ψ has already been determined as a function of ρ in Section 4, (23) is a linear differential equation for ε θ p . The boundary condition to this equation is
ε θ p = 0
for ρ = ρ c . Once Equation (23) has been solved, the distribution of ε r p and ε z p in the plastic region is found from (20) and the distribution of the total strains from (22).
The distribution of strains in the elastic region is determined from the solution for stress found in Section 4 and Hooke’s law.

6. Unloading

Let p f be the value of p at the end of loading. Then, the boundary conditions for the stage of unloading are
Δ σ r = 0
for ρ = a and
Δ σ r = σ 0 p f
for ρ = 1 . Here, Δ σ r is the increment of the radial stress after unloading ( Δ σ θ will stand for the increment of the circumferential stress). It is assumed that unloading is purely elastic. This assumption should be verified a posteriori. Equations (11) and (13), in which σ r should be replaced with Δ σ r and σ θ with Δ σ θ , are valid. An iterative procedure should be used for solving this system of equations together with the boundary conditions (25) and (26). Once this boundary value problem has been solved, the increment of strains is determined from Hooke’s law as
Δ ε r = Δ σ r ν Δ σ θ E ,   Δ ε θ = Δ σ θ ν Δ σ r E ,   Δ ε z = ν ( Δ σ r + Δ σ θ ) E
The distribution of residual stresses, σ r r e s and σ θ r e s , is given by
σ r r e s = σ r ( f ) + Δ σ r ,   σ θ r e s = σ θ ( f ) + Δ σ θ
Here, σ r ( f ) is the distribution of the radial stress and σ θ ( f ) is the distribution of the circumferential stress at the end of loading. These distributions have been found in Section 4. Substituting (28) into (2) provides the condition to verify that the process of unloading is purely elastic in the form
( σ r r e s ) 2 + ( σ θ r e s ) 2 σ r r e s σ θ r e s σ 0 2 Φ 2 0
Here, Equation (10) has been taken into account.
The distribution of residual strains, ε r r e s , ε θ r e s and ε z r e s , is given by
ε r r e s = ε r ( f ) + Δ ε r ,   ε θ r e s = ε θ ( f ) + Δ ε θ ,   ε z r e s = ε z ( f ) + Δ ε z
Here, ε r ( f ) is the distribution of the radial strain, ε θ ( f ) is the distribution of the circumferential strain and ε z ( f ) is the distribution of the axial strain at the end of loading. These distributions have been found in Section 5.

7. Illustrative Example

It is often assumed that material properties vary according to a power law along the radius of the disc [1,9,10,11,12,13]. In the case under consideration one possible variant of this law reads
Φ ( ρ ) = ρ m   and   η ( ρ ) = ρ n
In all calculations, ν = 0.3 , a = 0.3 , and n = 0.3 . The value of m varies in the range 0 m 0.3 . It is worthy of note that there is no need to prescribe the values of σ 0 and E 0 for numerical analysis since the stress components are proportional to σ 0 , and the strain components are proportional to k. It is assumed that ρ c = 0.8 . Then, the value of p f has been found from the stress solution given in Section 4. The dependence of p f on m is presented in Table 1. The variation of the radial and circumferential stresses with ρ at p = p f is depicted in Figure 1 and Figure 2, respectively. It has been verified that the yield criterion (2) is not violated in the elastic region. The variation of the radial, circumferential, and axial strains with ρ at p = p f is depicted in Figure 3, Figure 4 and Figure 5, respectively.
Using the values of p f found (Table 1) the system of Equations (11) and (13) together with the boundary conditions (25) and (26) has been solved for Δ σ r / σ 0 and Δ σ θ / σ 0 . Having this solution and the radial distribution of the radial and circumferential stresses at the end of loading (Figure 1 and Figure 2) the radial distribution of the residual stresses is determined from (28). The variation of the residual radial and circumferential stresses with ρ is shown in Figure 6 and Figure 7, respectively. Then, it has been verified that the inequality (29) is satisfied in the range a ρ 1 . Having the solution for Δ σ r / σ 0 and Δ σ θ / σ 0 , the increment of the strains is determined from (27). It is evident from (10) and (27) that these increments are proportional to k introduced in (9). The radial distribution of the residual strains is determined from (30). The variation of the residual radial, circumferential, and axial strains with ρ is shown in Figure 8, Figure 9 and Figure 10, respectively.

8. Discussion

The paper presents a general solution for the distribution of stress and strain in an functionally graded hollow disc subject to external pressure, followed by unloading. The solution is valid for any variation of the yield stress and Young’s modulus with radius if the plastic region initiates at, and then propagates from, the inner radius of the disc. The purely elastic solution is valid if p p e . This solution is known, and therefore is not considered in the present paper. However, the general elastic solution is used in the elastic region, ρ c ρ 1 . The stress solution follows from the solution of Equation (13). Then, the distribution of strains is determined from Hooke’s law shown in (1).
The constitutive equations of the classical flow theory of plasticity are adopted. In particular, Hooke’s law shown in (1) is used to connect the stress components and the elastic strain components. This law is valid in the elastic region. The von Mises yield criterion (2) is adopted in the plastic region. In this region, the stress components are connected to the strain rate components rather than to the strain components. The corresponding constitutive equation is the associated flow rule (3). The total strain components in the plastic region are given by (4). A detailed description of this material model can be found in any textbook on plasticity theory (for example, [29,30]).
For any given functions Φ ( ρ ) and η ( ρ ) involved in (10), the distribution of stress and strain in the plastic region, a ρ ρ c , can be calculated as follows. Two ordinary differential equations, Equations (15) and (23), should be solved numerically. These equations can be solved one by one. In particular, the dependence of ψ on ρ is found from (15). Then, this numerical function is substituted into (23). As a result, a linear differential equation for the circumferential plastic strain, ε θ p , is obtained. The solution of Equation (15) supplies the distribution of stresses in the plastic region according to (14). A remarkable feature of the strain solution is that the equations in (18), which are derived from the associated flow rule written in terms of plastic strain rates, can be immediately integrated with respect to the time to result in the equations in terms of plastic strains (Equation (20)). This feature of the solution significantly facilitates the solution. Another remarkable feature of the strain solution is that all strain components are proportional to k introduced in (9). It is seen from Equations (22) and (23). Therefore, simple scaling of any strain solution provides the strain solution for similar discs of material with the same other properties and geometry, but any value of k. To calculate the distribution of strains, it is first necessary to solve Equation (23) for the circumferential plastic strain. Then, the radial and axial plastic strains are readily found from (20) and the solution of (15). The elastic portions of strain in the plastic region are determined from (14) and Hooke’s law shown in (1). Finally, the total strains follow from (4).
The general solution found is used to find a numerical solution, assuming that the power laws shown in (31) are valid. The effect on m-value on the radial and circumferential stresses in an a = 0.3 disc at ν = 0.3 , n = 0.3 , and ρ c = 0.8 is illustrated in Figure 1 and Figure 2, respectively. The corresponding values of the external pressure are depicted in Table 1. It is seen from these figures that this effect is quite significant, especially on the circumferential stress. The associated distribution of strains is illustrated in Figure 3, Figure 4 and Figure 5 (the radial strain is depicted in Figure 3, the circumferential strain in Figure 4 and the axial strain in Figure 5). The effect of m-value on these distributions is not so pronounced as compared to the stress distributions. This is associated with the imposed condition that ρ c = 0.8 in all cases. If the value of p f were fixed, then the effect of m-value on the strain distributions would be more significant as compared to its effect on the stress distribution. In particular, it is seen from Table 1 that the value of p f is quite sensitive to the value of m.
Residual stress and strain fields after purely elastic unloading are also obtained. These distributions are given by Equations (28) and (30). The validity of the pure elastic solution at unloading should be verified by means of Equation (29). The distribution of residual radial and circumferential stresses is illustrated in Figure 6 and Figure 7, respectively. As in the case of the stress distributions at the end of loading, the effect of m-value is most significant on the circumferential stress. The effect of this value on the distribution of the residual radial strain is small (Figure 8). It is moderate in the case of the residual circumferential strain (Figure 9) and large in the case of the residual axial strain (Figure 10).

9. Conclusions

Stress and strain fields in an elastic/plastic functionally graded annular disc of constant thickness subject to external pressure are obtained under plane stress conditions. Residual stress and strain fields after purely elastic unloading are also obtained. From this work, the following conclusions can be drawn.
  • A remarkable feature of the strain solution is that the equations in (18), which are derived from the associated flow rule written in terms of plastic strain rates, can be immediately integrated with respect to the time to result in the equations in terms of plastic strains (Equation (20)). This significantly facilitates the solution.
  • Another remarkable feature of the strain solution is that all strain components are proportional to k introduced in (9). Therefore, simple scaling of any strain solution provides the strain solution for similar discs of material with the same other properties and geometry, but any value of k.
  • In the case of the stress solution, the effect of m-value involved in (31) is most significant on the distribution of the circumferential stress and the residual circumferential stress.
  • In the case of the strain solution, the effect of m-value is most significant on the distribution of the axial strain and the residual axial strain.
The method used in the present paper is a generalization of the method developed in [31] for homogeneous discs. It is evident from the solutions provided in [31] that the method can be more successfully adopted for disc subject to other loading conditions than those used in the present paper. In particular, the basic equations derived are independent of boundary conditions. Therefore, the solution of these equations used in conjunction with any other boundary conditions (of course, the boundary value problem should be axisymmetric) supplies the distribution of stress and strain. This will be the subject of a subsequent investigation.

Author Contributions

The general method has been developed by S.A., analytical treatment has been carried out by S.S., numerical results have been obtained by L.L.

Funding

S.A. acknowledges support from the Russian Ministry of Science and Higher Education (Project AAAA-A17-117021310373-3).

Acknowledgments

This work was performed while S.V. Strashnov was a visiting research fellow at Beihang University, Beijing, China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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]
  2. 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]
  3. Tutuncu, N. Stresses in thick-walled FGM cylinders with exponentially-varying properties. Eng. Struct. 2007, 29, 2032–2035. [Google Scholar] [CrossRef]
  4. Gönczi, D.; Ecsedi, I. Thermoelastic Analysis of Functionally Graded Hollow Circular Disc. Arch. Mech. Eng. 2015, 62, 5–18. [Google Scholar] [CrossRef] [Green Version]
  5. Yeo, W.H.; Purbolaksono, J.; Aliabadi, M.H.; Ramesh, S.; Liew, H.L. Exact solution for stresses/displacements in a multilayered hollow cylinder under thermo-mechanical loading. Int. J. Press. Vessels Pip. 2017, 151, 45–53. [Google Scholar] [CrossRef]
  6. Shabana, Y.M.; Elsawaf, A.; Khalaf, H.; Khalil, Y. Stresses minimization in functionally graded cylinders using particle swarm optimization technique. Int. J. Press. Vessels Pip. 2017, 154, 1–10. [Google Scholar] [CrossRef]
  7. Haghpanah Jahromi, B.; Farrahi, G.H.; Maleki, M.; Nayeb-Hashemi, H.; Vaziri, A. Residual stresses in autofrettaged vessel made of functionally graded material. Eng. Struct. 2009, 31, 2930–2935. [Google Scholar] [CrossRef]
  8. 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]
  9. Yildirim, 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]
  10. Jalali, M.H.; Shahriari, B. Elastic Stress Analysis of Rotating Functionally Graded Annular Disk of Variable Thickness Using Finite Difference Method. Math. Probl. Eng. 2018, 2018, 1–11. [Google Scholar] [CrossRef]
  11. Hosseini Kordkheili, S.A.; Naghdabadi, R. Thermoelastic analysis of a functionally graded rotating disk. Compos. Struct. 2007, 79, 508–516. [Google Scholar] [CrossRef]
  12. Hassani, A.; Hojjati, M.H.; Farrahi, G.; Alashti, R.A. Semi-exact elastic solutions for thermo-mechanical analysis of functionally graded rotating disks. Compos. Struct. 2011, 93, 3239–3251. [Google Scholar] [CrossRef]
  13. Dai, T.; Dai, H.-L. Thermo-elastic analysis of a functionally graded rotating hollow circular disk with variable thickness and angular speed. Appl. Math. Model. 2016, 40, 7689–7707. [Google Scholar] [CrossRef]
  14. Allam, M.N.M.; Tantawy, R.; Zenkour, A.M. Thermoelastic stresses in functionally graded rotating annular disks with variable thickness. J. Theor. Appl. Mech. 2018, 56, 1029–1041. [Google Scholar] [CrossRef]
  15. Leu, S.-Y.; Chien, L.-C. Thermoelastic analysis of functionally graded rotating disks with variable thickness involving non-uniform heat source. J. Therm. Stress. 2015, 38, 415–426. [Google Scholar] [CrossRef]
  16. Allam, M.N.M.; Tantawy, R.; Yousof, A.; Zenkour, A.M. Elastic and Viscoelastic Stresses of Nonlinear Rotating Functionally Graded Solid and Annular Disks with Gradually Varying Thickness. Arch. Mech. Eng. 2017, 64, 423–440. [Google Scholar] [CrossRef] [Green Version]
  17. Sondhi, L.; Sanyal, S.; Saha, K.; Bhowmick, S. Limit Elastic Speeds of Functionally Graded Annular Disks. FME Trans. 2018, 46, 603–611. [Google Scholar]
  18. Zheng, Y.; Bahaloo, H.; Mousanezhad, D.; Vaziri, A.; Nayeb-Hashemi, H. Displacement and Stress Fields in a Functionally Graded Fiber-Reinforced Rotating Disk with Nonuniform Thickness and Variable Angular Velocity. J. Eng. Mater. Technol. 2017, 139, 1–10. [Google Scholar] [CrossRef]
  19. Khorsand, M.; Tang, Y. Design functionally graded rotating disks under thermoelastic loads: Weight optimization. Int. J. Press. Vessels Pip. 2018, 161, 33–40. [Google Scholar] [CrossRef]
  20. Nejad, M.Z.; Rastgoo, A.; Hadi, A. Exact elasto-plastic analysis of rotating disks made of functionally graded materials. Int. J. Eng. Sci. 2014, 85, 47–57. [Google Scholar] [CrossRef]
  21. Nejad, M.Z.; Fatehi, P. Exact elasto-plastic analysis of rotating thick-walled cylindrical pressure vessels made of functionally graded materials. Int. J. Eng. Sci. 2015, 86, 26–43. [Google Scholar] [CrossRef]
  22. Hassani, A.; Hojjati, M.H.; Farrahi, G.H.; Alashti, R.A. Semi-exact solution for thermo-mechanical analysis of functionally graded elastic-strain hardening rotating disks. Commun. Nonlinear Sci. Numer. Simulat. 2012, 17, 3747–3762. [Google Scholar] [CrossRef]
  23. Haghpanah Jahromi, B.; Nayeb-Hashemi, H.; Vaziri, A. Elasto-Plastic Stresses in a Functionally Graded Rotating Disk. J. Eng. Mater. Technol. 2012, 134, 1–11. [Google Scholar]
  24. Mahdavi, E.; Ghasemi, A.; Alashti, R.A. Elastic–plastic analysis of functionally graded rotating disks with variable thickness and temperature-dependent material properties under mechanical loading and unloading. Aerosp. Sci. Technol. 2016, 59, 57–68. [Google Scholar] [CrossRef]
  25. Kalali, A.T.; Hadidi-Moud, S.; Hassani, B. Elasto-Plastic Stress Analysis in Rotating Disks and Pressure Vessels Made of Functionally Graded Materials. Latin Am. J. Solids Struct. 2016, 13, 819–834. [Google Scholar] [CrossRef] [Green Version]
  26. Pirumov, A.; Alexandrov, S.; Jeng, Y.-R. Enlargement of a Circular Hole in a Disc of Plastically Compressible Material. Acta Mech. 2013, 224, 2965–2976. [Google Scholar] [CrossRef]
  27. Çallıoglu, H.; Sayer, M.; Demir, E. Elastic–plastic stress analysis of rotating functionally graded discs. Thin-Walled Struct. 2015, 94, 38–44. [Google Scholar] [CrossRef]
  28. Zheng, Y.; Bahaloo, H.; Mousanezhad, D.; Mahdi, E.; Vaziri, A.; Nayeb-Hashemi, H. Stress analysis in functionally graded rotating disks with non-uniform thickness and variable angular velocity. Int. J. Mech. Sci. 2016, 116, 283–293. [Google Scholar] [CrossRef]
  29. Chakrabarty, J. Theory of Plasticity; McGraw-Hill Book Company: New York, NY, USA, 1987. [Google Scholar]
  30. Rees, D.W.A. Basic Engineering Plasticity; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  31. Alexandrov, S. Elastic/Plastic Discs under Plane Stress Conditions; Springer: Berlin, Germany, 2015. [Google Scholar]
Figure 1. Variation of the radial stress, σ r , with the dimensionless radius, ρ , at p = p f .
Figure 1. Variation of the radial stress, σ r , with the dimensionless radius, ρ , at p = p f .
Materials 12 00440 g001
Figure 2. Variation of the circumferential stress, σ θ , with the dimensionless radius, ρ , at p = p f .
Figure 2. Variation of the circumferential stress, σ θ , with the dimensionless radius, ρ , at p = p f .
Materials 12 00440 g002
Figure 3. Variation of the radial strain, ε r , with the dimensionless radius, ρ , at p = p f .
Figure 3. Variation of the radial strain, ε r , with the dimensionless radius, ρ , at p = p f .
Materials 12 00440 g003
Figure 4. Variation of the circumferential strain, ε θ , with the dimensionless radius, ρ , at p = p f .
Figure 4. Variation of the circumferential strain, ε θ , with the dimensionless radius, ρ , at p = p f .
Materials 12 00440 g004
Figure 5. Variation of the axial strain, ε z , with the dimensionless radius, ρ , at p = p f .
Figure 5. Variation of the axial strain, ε z , with the dimensionless radius, ρ , at p = p f .
Materials 12 00440 g005
Figure 6. Variation of the residual radial stress, σ r r e s , with the dimensionless radius, ρ .
Figure 6. Variation of the residual radial stress, σ r r e s , with the dimensionless radius, ρ .
Materials 12 00440 g006
Figure 7. Variation of the residual circumferential stress, σ θ r e s , with the dimensionless radius, ρ .
Figure 7. Variation of the residual circumferential stress, σ θ r e s , with the dimensionless radius, ρ .
Materials 12 00440 g007
Figure 8. Variation of the residual radial strain, ε r r e s , with the dimensionless radius, ρ .
Figure 8. Variation of the residual radial strain, ε r r e s , with the dimensionless radius, ρ .
Materials 12 00440 g008
Figure 9. Variation of the residual circumferential strain, ε θ r e s , with the dimensionless radius, ρ .
Figure 9. Variation of the residual circumferential strain, ε θ r e s , with the dimensionless radius, ρ .
Materials 12 00440 g009
Figure 10. Variation of the residual axial strain, ε z r e s , with the dimensionless radius, ρ .
Figure 10. Variation of the residual axial strain, ε z r e s , with the dimensionless radius, ρ .
Materials 12 00440 g010
Table 1. Dependence of the value of pressure at the end of loading on the value of m introduced in (31).
Table 1. Dependence of the value of pressure at the end of loading on the value of m introduced in (31).
m p f
00.79
0.10.75
0.20.71
0.30.68

Share and Cite

MDPI and ACS Style

Strashnov, S.; Alexandrov, S.; Lang, L. Description of Residual Stress and Strain Fields in FGM Hollow Disc Subject to External Pressure. Materials 2019, 12, 440. https://doi.org/10.3390/ma12030440

AMA Style

Strashnov S, Alexandrov S, Lang L. Description of Residual Stress and Strain Fields in FGM Hollow Disc Subject to External Pressure. Materials. 2019; 12(3):440. https://doi.org/10.3390/ma12030440

Chicago/Turabian Style

Strashnov, Stanislav, Sergei Alexandrov, and Lihui Lang. 2019. "Description of Residual Stress and Strain Fields in FGM Hollow Disc Subject to External Pressure" Materials 12, no. 3: 440. https://doi.org/10.3390/ma12030440

APA Style

Strashnov, S., Alexandrov, S., & Lang, L. (2019). Description of Residual Stress and Strain Fields in FGM Hollow Disc Subject to External Pressure. Materials, 12(3), 440. https://doi.org/10.3390/ma12030440

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