Skip to main content
  • Original Paper
  • Open access
  • Published:

An analytical investigation of elastic-plastic deformation of FGM hollow rotors under a high centrifugal effect

Abstract

Functionally graded material shafts are the main part of many modern rotary machines such as turbines and electric motors. The purpose of this study is to present an analytical solution of the elastic-plastic deformation of functionally graded material hollow rotor under a high centrifugal effect and finally determine the maximum allowed angular velocity of a hollow functionally graded material rotating shaft. Introducing non-dimensional parameters, the equilibrium equation has been analytically solved. The results for variable material properties are compared with the homogeneous rotor and the case in which Young’s modulus is the only variable while density and yield stress are considered to be constant. It is shown that material variation has a considerable effect on the stress and strain components and radial displacement. Considering variable density and yield stress causes yielding onset from inner, outer, or simultaneously from both inner and outer rotor shaft radius in contrast to earlier researches that modulus of elasticity was the only variable. The effects of the density on the failure of a functionally graded material elastic fully plastic in a hollow rotating shaft are investigated for the first time in this study with regard to Tresca’s yield criterial. Numerical simulations are used to verify the derived formulations which are in satisfying agreement.

Introduction

Functionally graded materials (FGM) are finding vast applications in different rotary systems such as DC motors with a magnetic membrane and chemical resistant hydraulic motors (Mahamood & Akinlabi, 2017), gas turbine rotors (Bahaloo, Papadopolus, & Ghosha, 2016; Klocke, Klink, & Veselovac, 2014; Lal, Jagtap, & Singh, 2013), and modern vehicle drive train systems (Kaviprakash, Kannan, Lawrence, & Regan, 2014; Lee, Kim, Kim, & Kim, 2004; Moorthy, Mitiku, & Sridhar, 2013). Computing different stresses and the radial displacement of FGM rotors are required to determine the maximum allowed angular velocity (Nino, Hirai, & Watanabe, 1987). Timoshenko (Timoshenko & Goodier, 1970), Mendelson (Mendelson, 1968), Chakrabarty (Chakrabarty, 2006) and Mack (Mack, 1991) analyzed a homogenous rotor. You analyzed a rotating FGM disk (You, You, Zhang, & Li, 2007) and Dai considered the magnetic properties of the FGM disk (Dai & Dai, 2017). Fukui and Yamanaka studied elastic analysis for thick-walled FGM tubes subjected to internal pressure (Fukui & Yamanaka, 1991). Figueiredo studied FGM pipes (Figueiredo, Borges, & Rochinha, 2008). Tutunku and Ozturk determined solutions for stresses in FGM pressure vessels (Tutuncu & Ozturk, 2001). Jabbari (Jabbari, Sohrabpour, & Eslami, 2002) and Ansari (AnsariSadrabadi et al., 2017) investigated mechanical and thermal stresses in an FGM hollow cylinder under symmetric loads. You considered an FGM pressurized sphere with a nonlinear variable modulus of elasticity in a radial direction (You, Zhang, & You, 2005). Dai et al. studied a pressurized magneto elastic FGM tube (Dai, Fu, & Dong, 2006). Hosseini et al. analyzed the thermo-elastic behavior of an FG rotating disk (HosseiniKordkheili & Naghdabadi, 2006). Duc (Duc, Lee, Nguyen-Thoi, & Thang, 2017; Duc, Thang, Dao, & Vantac, 2015), Khoa (Khoa, Thiem, Thiem, & Duc, 2019), and El-Haina (El-Haina, Bakora, Bousahla, Tounsi, & Mahmoud, 2017) considered the buckling problem in their research. Thom (Thom, Kien, Duc, Duc, & Tinh, 2017) analyzed a two-dimensional analysis on an FGM plane by plane strain theories. Eraslan (Eraslan & Akis, 2006a) gave an analytical solution for rotating disks and tubes in plane stress and plane strain state, and studied stress solutions of FGM shafts and disks (Eraslan & Akis, 2006b; Eraslan & Akis, 2006c). Kargarnovin et al. (Kargarnovin, Faghidian, & Arghavani, 2007) investigated FGM circular plates with arbitrary rotational symmetric load. Akis (Akis & Eraslan, 2007) studied a rotating FGM shaft problem in the elastic-plastic state of stress with a variable modulus of elasticity. Tsiatas (Tsiatas & Babouskos, 2017) worked on torsional FGM bar. Akis studied the elasticity solution for thick-walled FG spherical pressure vessels with linearly and exponentially varying properties (Akis, 2009). ZamaniNejad and Rahimi studied the elasticity of an FGM rotating cylindrical pressure vessels (ZamaniNejad & Rahimi, 2010). Peng and Li investigated an orthotropic hollow rotating disk with a variable modulus of elasticity and density (Peng & Li, 2012). Some others studied creep for FGM material under thermal condition (Bose & Rattan, 2018; Khanna, Gupta, & Nigam, 2017; Zharfi & EkhteraeiToussi, 2018). Yildirim (Yildirim & Tutuncu, 2018), Seraj (Seraj & Ganesan, 2018), Bahaadini (Bahaadini & Saidi, 2018), Swaminathan (Swaminathan, Naveenkumar, Zenkour, & Carrera, 2015), Duc (Duc, 2013; Duc & Cong, 2018; Duc, Homayoun, Quan, & Khoa, 2019; Duc, Nguyen, & Khoa, 2017; Duc, Tran, & Cong, 2016), and Bouderba (Bouderba, Houari, Tounsi, & Mahmoud, 2016) worked on rotor instabilities and vibrations under different conditions. Burzyński (Burzyński, Chróścielewski, Daszkiewicz, & Witkowski, 2018) worked on a FEM method to understand elasto-plastic behaviors of FGM shells, and Mathew (Mathew, Natarajan, & Pañeda, 2018) considered size effects in his researches. Duc (Duc, 2016a; Duc, 2016b; Duc et al., 2015; Duc, Bich, & Cong, 2016; Duc, Khoa, & Thiem, 2018; Duc, Kim, & Chan, 2018; Duc, Thuy Anh, & Cong, 2014) specifically studied thermal effects such as buckling, thermal instability, and dynamic thermal loads circular sections. The authors (Torabnia, Hemati, & Aghajanib, 2019) considered the elastic behavior of a hollow FGM rotor.

Although the previous studies are valuable, none of them considered the plastic effects in an analytical model. All previous jobs used a numerical method such as FEM to solve the plastic model. In the present work, the analysis is based on small deformation theory. The shaft is assumed to be infinitely long (plane strain). The maximum allowed angular velocity has been defined as the angular velocity in which yielding initiates based on Tresca’s criterion. Non-dimensional parameters are introduced based on the geometry and material parameters. Stress components are derived using generalized Hook’s law. To identify the stress components ordering, non-dimensional stress components are plotted for the special case of equal exponent parameters with the variable radius ratio. The results show when the exponent parameters vary between − 2 and 2, hoop stress and radial stress components are the largest and the smallest stress components. The effect of variation of density and yield stress is investigated on the maximum allowed angular velocity and has a considerable effect on the stress distribution and yielding initiation and the maximum allowed angular velocity. For the first time, density variation is considered with variable density and radius ratio of a hollow rotor on elastic and plastic behavior and maximum allowed angular velocity are discussed (Fig. 1).

Fig. 1
figure 1

Schematic of the rotor

Methods/experimental

In this section, the aims and methodology of the study presented by an explanation of the governing equations of a hollow rotor with variable properties through its geometry. Material properties in an FGM may vary in any direction. Here, modulus of elasticity, density, and yield stress are functions of radial dimension:

$$ E(r)={E}_0{\left(r/b\right)}^{n_E},\rho (r)={\rho}_0{\left(r/b\right)}^{n_{\rho }},{\sigma}_{\mathrm{Y}}(r)={\sigma}_0{\left(r/b\right)}^{n_{\sigma }} $$
(1)

The material properties modeled with the power-law function. Different exponent parameters allow various shapes for material variation. By formulating in the cylindrical coordinate system (r, θ, z) for an infinitely long tube which rotates about longitude axis (Timoshenko & Goodier, 1970):

$$ \frac{d}{dr}\left(r{\sigma}_r\right)-{\sigma}_{\theta }=-\rho {r}^2{\omega}^2 $$
(2)

The strain and radial displacement relation is:

$$ {\varepsilon}_r= du(r)/ dr,{\varepsilon}_{\theta }=u\;(r)/r $$
(3)

Plane strain condition is due to a long tube which causes the zero value for the longitude strain. Manipulating stress-strain and radial displacement:

$$ {\sigma}_r=\frac{E}{\left(1+v\right)\left(1-2v\right)}\left(\left(1-v\right)\frac{du(r)}{dr}+v\frac{u(r)}{r}\right),{\sigma}_{\theta }=\frac{E}{\left(1+v\right)\left(1-2v\right)}\left(\left(1-v\right)\frac{u(r)}{r}+v\frac{du(r)}{dr}\right),{\sigma}_z=v\left({\sigma}_r+{\sigma}_{\theta}\right) $$
(4)

Substituting (4) into (2) in elastic region:

$$ {r}^2\frac{d^2}{d{r}^2}u(r)+\left(1+{n}_E\right)r\left(\frac{d}{dr}u(r)\right)-\frac{1-v\left(1+{n}_E\right)}{1-v}u(r)=-\frac{\left(1+v\right)\left(1-2v\right){\rho}_0{\omega}^2}{\left(1-v\right){E}_0}{b}^{\left({n}_E- n\rho \right)}{r}^{\left(3+ n\rho -{n}_E\right)} $$
(5)

The general reformed solution of (5) is:

$$ \overline{u}\left(\overline{r}\right)={C}_1{\overline{r}}^{\left(\frac{-{n}_E-k}{2}\right)}+{C}_2{\overline{r}}^{\left(\frac{-{n}_E+k}{2}\right)}-{A}_1{\overline{\omega}}^2{\overline{r}}^{\left( n\rho -{n}_E+3\right)} $$
(6)

A solution of (6) is simplified and taken into a non-dimensional form to be independent of material properties. Non-dimensional quantities are presented:

$$ {\displaystyle \begin{array}{l}\overline{r}=r/b,\overline{h}=a/b,\overline{u}=u/b,{\overline{\omega}}^2={\rho}_0{\omega}^2{b}^2/{E}_0\\ {}k=\frac{\sqrt{{n_E}^2-2{n_E}^2v+{v}^2{n_E}^2+4-8v+4{v}^2-4v{n}_E+4{v}^2{n}_E}}{1-v}\\ {}{A}_1=\frac{2{v}^2+v-1}{\left(4+{n}_{\rho}\right)\left({n}_{\rho }-{n}_E+2\right)v-{n_{\rho}}^2+\left({n}_E-6\right){n}_{\rho }-8+3{n}_E}\end{array}} $$
(7)

ω̅ defined as the non-dimensional rotating velocity. Substituting radial displacement into (5):

$$ {\displaystyle \begin{array}{l}{\overline{\sigma}}_r=\frac{{\overline{r}}^{n_E}}{\left(1+v\right)\left(1-2v\right)}\left\{\begin{array}{l}v\left({C}_1{\overline{r}}^{m_1}+{C}_2{\overline{r}}^{m_2}-{A}_1{\overline{\omega}}^2{\overline{r}}^{\left(2+{n}_{\rho }-{n}_E\right)}\right)+\\ {}\left(1-v\right)\left(\begin{array}{l}{C}_1{\overline{r}}^{m_1}\frac{\left(-{n}_E-k\right)}{2}+{C}_2{\overline{r}}^{m_2}\frac{\left(-{n}_E+k\right)}{2}\\ {}-{A}_1{\overline{\omega}}^2\left(3+{n}_{\rho }-{n}_E\right)\overline{r}\left(2+{n}_{\rho }-{n}_E\right)\end{array}\right)\end{array}\right\}\\ {}{\overline{\sigma}}_{\theta }=\frac{{\overline{r}}^{n_E}}{\left(1+v\right)\left(1-2v\right)}\left\{\begin{array}{l}\left(1-v\right)\left({C}_1{\overline{r}}^{m_1}+{C}_2{\overline{r}}^{m_2}-{A}_1{\overline{\omega}}^2{\overline{r}}^{\left(2+{n}_{\rho }-{n}_E\right)}\right)+\\ {}v\left(\begin{array}{l}{C}_1{\overline{r}}^{m_1}\frac{\left(-{n}_E-k\right)}{2}+{C}_2{\overline{r}}^{m_2}\frac{\left(-{n}_E+k\right)}{2}\\ {}-{A}_1{\overline{\omega}}^2\left(3+{n}_{\rho }-{n}_E\right){\overline{r}}^{\left(2+{n}_{\rho }-{n}_E\right)}\end{array}\right)\end{array}\right\}\\ {}{\overline{\sigma}}_z=\frac{v\times {\overline{r}}^{n_E}}{\left(1+v\right)\left(1-2v\right)}\left\{{C}_1{\overline{r}}^{m_1}{m}_3+{C}_2{\overline{r}}^{m_2}{m}_4-{A}_1{\overline{\omega}}^2{\overline{r}}^{\left(2+{n}_{\rho }-{n}_E\right)}\left(4+{n}_{\rho }-{n}_E\right)\right\}\end{array}} $$
(8)

σ̅i stands for non-dimensional stress which is defined in the form σ̅i = σi/E0. The constants used in (8) are:

$$ {m}_1=\frac{-{n}_E-k-2}{2},{m}_2=\frac{-{n}_E+k-2}{2},{m}_3=\frac{-{n}_E-k+2}{2},{m}_4=\frac{-{n}_E+k+2}{2} $$
(9)

To obtain C1 and C2 in radial displacement (7), two boundary conditions are needed. Since no pressure is applied to the inner and outer surfaces of the rotor, the boundary conditions are considered as σ̅Y(r̅=h̅) = 0 & σ̅Y(r̅=1) = 0. Constants C1 and C2 are:

$$ {C}_1=-2{A}_1{A}_2{\overline{\omega}}^2\frac{\left[{\overline{h}}^{n_p+4}-{\overline{h}}^{-{m}_1}\right]}{\left({\overline{h}}^{-{m}_1}-{\overline{h}}^{-{m}_2}\right)},{C}_2=-2{A}_1{A}_3{\overline{\omega}}^2\frac{\left[{\overline{h}}^{n_{\rho }+4}-{\overline{h}}^{-{m}_2}\right]}{\left({\overline{h}}^{-{m}_1}-{\overline{h}}^{-{m}_2}\right)} $$
(10)

For A2 and A3:

$$ {A}_2=\frac{3-v{n}_{\rho }+v{n}_E-2v+{n}_{\rho }-{n}_E}{-{n}_E-k+v{n}_E+ vk+2v},{A}_3=\frac{3-v{n}_{\rho }+v{n}_E-2v+{n}_{\rho }-{n}_E}{-2v-v{n}_E-k+{n}_E+ vk} $$
(11)

Tresca’s criterion is considered to determine yield condition and allowed the angular velocity of the shaft. As the results show, yielding is a function of exponent parameters of material variables (nE, nρ, nσ). In this paper, the results are discussed on the equality of the exponents of material variables. According to Fig. 3 for the state of equal exponents in the range of − 2 ≤ ni ≤ 2 and 0.5 ≤ ≤ 1, the stress components have the order of σθ ≥ σz ≥ σr. The yield criterion is in the form of σθ-σr = σY. Rearranging into a non-dimensional form gives:

$$ {\overline{\sigma}}_{\theta }-{\overline{\sigma}}_{\overline{r}}=\frac{\sigma_0}{E_0}{\overline{r}}^{n_{\sigma }}={\sigma}_0{\overline{r}}^{n_{\sigma }}\Rightarrow {\overline{\sigma}}_{\mathrm{Tresca}}=\frac{\left({\overline{\sigma}}_{\theta }-{\overline{\sigma}}_{\overline{r}}\right)}{{\overline{\sigma}}_0}{-}_{\overline{r}}{n}_{\sigma } $$
(12)

By substituting hoop and radial stresses in the yield’s criterion following equation formed.

$$ {\overline{\sigma}}_{\mathrm{Tresca}}=\frac{{\overline{\omega}}^2{A}_1{\overline{r}}^{n_E}}{\left(1+v\right){\overline{\sigma}}_0}\left\{\frac{2\left[{}_{+{m}_2{A}_3\left({\overline{h}}^{4+{n}_{\rho }}-{\overline{h}}^{-{m}_2}\right){\overline{r}}^{m_2}}^{m_1{A}_2\left({\overline{h}}^{4+{n}_{\rho }}-{\overline{h}}^{-{m}_1}\right){\overline{r}}^{m_1}}\right]}{\begin{array}{l}{\overline{h}}^{-{m}_1}-{\overline{h}}^{-{m}_2}\\ {}+\left(3+{n}_{\rho }-{n}_E\right){\overline{r}}^{\left(2-{n}_E+{n}_{\rho}\right)}\end{array}}\right\}-{\overline{r}}^{n_{\sigma }} $$
(13)

Yielding occurs when above equation equals zero for the corresponding load parameters such as angular velocity ω̅, yield stress σ0, and modulus of elasticity E. These parameters are rearranged and defined together as non-dimensional loading parameter (NLP):

$$ NLP=\frac{\overline{\omega}}{\sqrt{{\overline{\sigma}}_0}} $$
(14)

NLP in the causes of which σ̅Tresca = 0 is called the maximum angular velocity of the shaft. As it is shown in Fig. 4 considering variable modulus of elasticity, density, and yield stress may cause yielding starts from the inner, outer, or simultaneously from the inner and outer surfaces of the shaft. The plastic region grows through the radial direction by increasing the angular velocity of the shaft results in raising the plastic elastic region ratio. Hence, determining the effect of radius ratio on stress ordonnance and then resumption of yielding initiation for the remaining elastic region is needed. Equilibrium equation of rotating tube (2) is independent of the elastic-plastic behavior of the material. Using Tresca’s yield criterion will give radial stress as follows:

$$ {\sigma}_r=\left({\sigma}_0/n\right){\overline{r}}^n-{\rho}_0{\overline{r}}^m\left[{r}^2{\omega}^2/\left(m+2\right)\right]+{C}_3 $$
(15)

To verify different conditions, it is supposed that yielding will initiate from both surfaces of the tube. For the yielding initiation from the inner radius of the shaft, the proposed boundary condition to determine C3 is (σ̅r(r̅=h̅) = 0). Substituting C3 into radial stress and using stress relations will result:

$$ {\displaystyle \begin{array}{l}{\overline{\sigma}}_{\overline{r}}=\frac{{\overline{\sigma}}_0}{n_{\sigma }}\left({\overline{r}}^{n_{\sigma }}-{\overline{h}}^{n_{\sigma }}\right)+\frac{{\overline{\omega}}^2}{n_{\rho }+2}\left({\overline{h}}^{n_{\rho }+2}-{\overline{r}}^{n_{\rho }+2}\right)\\ {}{\overline{\sigma}}_{\theta }=\frac{{\overline{\sigma}}_0}{n_{\sigma }}\left({\overline{r}}^{n_{\sigma }}\left(1+{n}_{\sigma}\right)-{\overline{h}}^{n_{\sigma }}\right)+\frac{{\overline{\omega}}^2}{n_{\rho }+2}\left({\overline{h}}^{n_{\rho }+2}-{\overline{r}}^{n_{\rho }+2}\right)\\ {}{\overline{\sigma}}_z=v\left(\frac{{\overline{\sigma}}_0}{n_{\sigma }}\left(-2{\overline{h}}^{n_{\sigma }}+{\overline{r}}^{n_{\sigma }}\left(2+{n}_{\sigma}\right)\right)+\frac{2{\overline{\omega}}^2}{n_{\rho }+2}\left({\overline{h}}^{n_{\rho }+2}-{\overline{r}}^{n_{\rho }+2}\right)\right)\end{array}} $$
(16)

Constants C1 and C2 related to the elastic region are obtained using stress continuity condition in different regions (σ̅r(r̅ = r̅ep)elastic = σ̅r(r̅ = r̅ep)plastic, σ̅r(r̅ = 1)elastic = 0, ep is the elastic-plastic border of the shaft).

In the case of yielding initiation from the outer radius of the shaft, the following boundary condition is used to determine C3 (σ̅r(r̅ = 1) = 0). The plastic stresses are:

$$ {\displaystyle \begin{array}{l}{\overline{\sigma}}_r=\frac{{\overline{\sigma}}_0}{n_{\overline{\sigma}}}\left[{\overline{r}}^{n_{\sigma }}-1\right]+\frac{{\overline{\omega}}^2}{n_{\rho }+2}\left[1-{\overline{r}}^{n_{\rho +2}}\right]\\ {}\left({n}_i\to 0\right)\Rightarrow {\overline{\sigma}}_r=\left(\frac{{\overline{\omega}}^2}{2}\right)\left(1-{\overline{r}}^2\right)+{\overline{\sigma}}_0\ln \left(r/a\right)\\ {}{\overline{\sigma}}_{\theta }=\frac{{\overline{\sigma}}_0}{n_{\sigma }}\left({\overline{r}}^{n_{\sigma }}\left(1+{n}_{\sigma}\right)-1\right)+\frac{{\overline{\omega}}^2}{n_{\rho }+2}\left(1-{\overline{r}}^{n_{\rho }+2}\right)\\ {}\left({n}_i\to 0\right)\Rightarrow {\overline{\sigma}}_{\theta }=\frac{{\overline{\omega}}^2}{2}\left(1-{\overline{r}}^2\right)+{\overline{\sigma}}_0\ln \left(r/a\right)+{\overline{\sigma}}_0\\ {}{\overline{\sigma}}_z=v\left(\frac{{\overline{\sigma}}_0}{n_{\sigma }}\left({\overline{r}}^{n_{\sigma }}\left(2+{n}_{\sigma}\right)-2\right)+\left[\frac{2{\overline{\omega}}^2}{n_{\rho }+2}\right]\left(1-{\overline{r}}^{n_{\rho +2}}\right)\right)\\ {}\left({n}_i\to 0\right)\Rightarrow {\overline{\sigma}}_z=v\left({\overline{\omega}}^2\left({1}^2-{\overline{r}}^2\right)+2{\overline{\sigma}}_0\ln \left(r/a\right)\right)\end{array}} $$
(17)

In the state of yielding initiation simultaneously from inner and outer radii of the shaft, the constants C1 and C2 related to radial elastic displacement and C3 and C4 related to radial plastic displacement and also rep1 and rep2 should be obtained simultaneously. Six equations are needed:

$$ {\displaystyle \begin{array}{l}{\overline{\sigma}}_{\overline{r}}{\left|{}^{elastic}\left(\overline{r}={\overline{r}}_{ep1}\right)={\overline{\sigma}}_{\overline{r}}\right|}^{elastic}\left(\overline{r}={\overline{r}}_{ep1}\right)\\ {}\overline{u}{\left|{}^{elastic}\left(\overline{r}={\overline{r}}_{ep1}\right)=\overline{u}\right|}^{plastic}\left(\overline{r}={\overline{r}}_{ep1}\right)\\ {}\overline{u}{\left|{}^{elastic}\left(\overline{r}={\overline{r}}_{ep1}\right)=\overline{u}\right|}^{plastic}\left(\overline{r}={\overline{r}}_{ep1}\right)\\ {}{\overline{\sigma}}_{\overline{r}}{\left|{}^{elastic}\left(\overline{r}={\overline{r}}_{ep2}\right)={\overline{\sigma}}_{\overline{r}}\right|}^{elastic}\left(\overline{r}={\overline{r}}_{ep2}\right)\\ {}\overline{u}{\left|{}^{elastic}\left(\overline{r}={\overline{r}}_{ep2}\right)=\overline{u}\right|}^{plastic}\left(\overline{r}={\overline{r}}_{ep2}\right)\\ {}{\overline{\sigma}}_{\theta }{\left|{}^{elastic}\left(\overline{r}={\overline{r}}_{ep2}\right)-{\overline{\sigma}}_{\overline{r}}\right|}^{elastic}\left(\overline{r}={\overline{r}}_{ep2}\right)={\overline{\sigma}}_Y\end{array}} $$
(18)

Associated flow rule for this state of stress order (Akis & Eraslan, 2007) is εθp = -εrp and εzp = 0. Superscripts e and p refer to elastic and plastic states (Fig. 2).

$$ {\displaystyle \begin{array}{l}{\varepsilon}^T={\varepsilon}^p+{\varepsilon}^e\\ {}{\varepsilon}^p={\varepsilon}_r^p+{\varepsilon}_{\theta}^p+{\varepsilon}_z^p\\ {}{\varepsilon}^e={\varepsilon}_r^e+{\varepsilon}_{\theta}^p+{\varepsilon}_z^e\end{array}} $$
(19)
Fig. 2
figure 2

Plastic deformations initiate from (a) the internal surface, (b) the outer surface of the rotor, and (c) both

The associated flow rule expresses that the total plastic strain equals zero (εp = 0). Hence, total elastic and plastic strains are as follows:

$$ {\varepsilon}^T={\varepsilon}^p+{\varepsilon}^e={\varepsilon}_r^e+{\varepsilon}_{\theta}^e=\frac{d\overline{u}}{d\overline{r}}+\frac{\overline{u}}{\overline{r}} $$
(20)

By knowing general stress-strain relations and using Hook’s general law and Tresca’s yield criterion, the stress-displacement equation becomes:

$$ {\displaystyle \begin{array}{l}{\varepsilon}_{ij}=\frac{1}{E}\left({\sigma}_{ij}-v\left({\sigma}_{kk}-{\sigma}_{ij}\right)\right);\sigma ij=\frac{\partial u}{\partial {\varepsilon}_{ij}};i,j,k=x,y,z\\ {}\frac{d\overline{u}}{d\overline{r}}+\frac{\overline{u}}{\overline{r}}=\frac{1}{{\overline{r}}^{n_E}}\left(\left(1-v-{v}^2\right)\left(2{\overline{\sigma}}_r+{\overline{\sigma}}_Y\right)\right)\end{array}} $$
(21)

Substituting obtained plastic stresses into the above relation and rearranging it, we have:

$$ \frac{d\overline{u}}{d\overline{r}}+\frac{\overline{u}}{\overline{r}}=\frac{2\left(1-v-2{v}^2\right)}{{\overline{r}}^{n_E}}\left\{\begin{array}{l}\frac{{\overline{\sigma}}_0}{n_{\sigma }}\left(\left(\frac{n_{\sigma }}{2}+1\right){\overline{r}}^{n_{\sigma }}-{\overline{h}}^{n_{\sigma }}\right)\\ {}+\frac{{\overline{\omega}}^2}{n_{\rho }+2}\left({\overline{h}}^{n_{\rho }+2}-{\overline{r}}^{n_{\rho }+2}\right)\end{array}\right\} $$
(22)

A non-dimensional solution of the above equation gives plastic radial displacement as follows:

$$ \frac{u(r)}{r}=\frac{1-v-2{v}^2}{n_{\rho }+2}\times \left\{\begin{array}{l}\frac{-2{\overline{r}}^{-{n}_E}}{n_{\sigma}\left(-2+{n}_E\right)}\left(\begin{array}{l}{\sigma}_0{\overline{h}}^{n_{\sigma }}\left({n}_{\rho }+2\right)\\ {}-{\overline{\omega}}^2{n}_{\sigma }{\overline{h}}^{n_{\rho }+2}\end{array}\right)-\frac{2{\overline{\omega}}^2{\overline{r}}^{-{n}_E+{n}_{\rho }+2}}{-4+{n}_E-{n}_{\rho }}\\ {}+\frac{{\overline{\sigma}}_0\left({n}_{\sigma }{n}_{\rho }+2{n}_{\sigma }+2{n}_{\rho }+4\right){\overline{r}}^{-{n}_E+{n}_{\sigma }}}{n_{\sigma}\left(-2+{n}_E-{n}_{\sigma}\right)}\end{array}\right\}+\frac{C_4}{r^2} $$
(23)

To obtain C4, the continuity condition of radial displacement through the elastic and plastic border is considered.

$$ u{\left|{}^{elastic}\left(\overline{r}={\overline{r}}_{ep}\right)=u\right|}^{plastic}\left(\overline{r}={\overline{r}}_{ep}\right) $$
(24)

Results and discussion

Elastic results

Verification has been done by comparing results with articles discussed on homogenous materials and prior FGM articles which are shown on subsequent plots and considering modulus of elasticity as the only variable property of the material. nE = nρ = nσ = 0 (the homogenous material condition) for the limit of (21) and ni = 0 creates:

$$ \varOmega =\overline{\omega}/\sqrt{{\overline{\sigma}}_0}=2\;\left(\overline{h}\right)\;\sqrt{\left(1-v\right)/\left(1-2v\right)+\left(3-2v\right)\;{\left(\overline{h}\right)}^2} $$
(25)

The above equation is the maximum allowed angular velocity in a homogeneous tube (Nino et al., 1987). The results are discussed for − 2 ≤ ni ≤ 2 and = 0.5 and = 0.55 for ν = 0.3 (Akis & Eraslan, 2007; Dai et al., 2006). To form Tresca’s yield criterion, the stress collocation must be determined. Hoop, radial, and longitude stresses are plotted for different h̅ = a/b ratios and different exponent parameters for both constant and variable density (Fig. 3).

Fig. 3
figure 3

Effect of density and radius ratio on stress collocation for ni = 2 (a, b) and ni = − 2 (c, d)

When ni = 2, the Hoop and longitude stress rising up for higher r̅ = r/b, but radial stresses have a peak in the midrange of r̅. Hoop and longitude stresses have a higher value of dimensionless stress components for constant density in comparison with the variable density. Another fact is in > 0.9 all studied stresses remain constant for different . On the other hand, for ni = − 2, Hoop and longitude stresses descending, but there is no significant change in radial stress trends. The effects of the constant and variable density in values of dimensionless components are different and the higher values belong to variable density. These effects are identified in this study for the first time. Figure 3 reveals that the hoop stress is in maximum and radial stress is in the minimum value for − 2 ≤ ni ≤ 2 and 0.5 ≤ h̅ ≤ 1; hence, the Tresca’s yield criterion is as defined before. Also, higher ratios make the results more linear. In previous articles on the elastic-plastic behavior of FG rotating tube, the only variable of the material is defined as modulus of elasticity (Akis & Eraslan, 2007; Tsiatas & Babouskos, 2017). Circumferential stress has a smaller value when considering variable density. This trend is similar through a/b = 0.5 to 1 for ni = 2, but these results are reverse for ni = − 2. This phenomenon could be explained by considering power ni. The power value sign is the reverse of material distribution through rotor wall thickness. For instance, when power law is a positive value, material distributed at outer radius of the rotor has greater density. So negative value of power will reverse material distribution which will result in different stress orders.

Plastic region

To validate the model for plastic deformations, a variable module of elasticity material (nE = 1.3826, nρ = nσ) considered to compare with (Akis & Eraslan, 2007). Figure 4a and b show that the results are the same for different NLP and rep/b. As Fig. 4a and b show, considering variable modulus of elasticity and constant density and yield stress limit may cause yielding initiation from an inner and outer radius of the shaft, simultaneously. High angular velocities create high centrifugal forces that made plastic deformations in the rotor. As shown in Fig. 4c, considering variable modulus of elasticity, density, and yield stress limit with the equal exponential rates will cause yielding from the inner radius of the shaft (nE = nρ = nσ = ni and − 2 ≤ ni ≤ 2). Also, the homogenous behavior of the rotating shaft is obtained (nE = nρ = nσ = 0). As it is depicted in Fig. 4c, higher plastic growth happened for higher NLPs and higher rotational speed as expected.

Fig. 4
figure 4

NLP plotted in different non-dimensional elasto-plastic boundary radius for (a) different nis, (b, c) nE = 1.3826, nρ = nσ, comparing Akis results (Duc, 2016b) in different modulus of elasticity

Plastic growth through radial coordinate by increasing the angular velocity of the shaft is shown in Fig. 5 for the different exponent of parameters. Maximum elastic and plastic velocity are also shown. In both, increasing ni results in a reduction of non-dimensional loading parameter which causes yielding at lower speeds for the rotor. This is happening because of the lowering of the average of the material properties in higher exponential rates of material change. Neglecting yield and density changes make a considerable error not only in the calculation of non-dimensional loading parameters but also in the determination of yielding initiation point.

Fig. 5
figure 5

Ωfp and Ωy to the loading parameter exponent nj ( = 0.5)

In Fig. 6, the elastic-plastic stresses are plotted for the state of plastic growth. To compare obtained results, plastic radial displacement is plotted for two conditions: considering variable modulus elasticity as the only variable of the material as Eraslan and Akis (Akis & Eraslan, 2007) (Fig. 6a) and considering variable modulus of elasticity, density, and yield stress limit as discussed above (Fig. 6b–d). Radial and hoop strains are plotted to verify and compare results. The plotted results are similar to the results from reference (Akis & Eraslan, 2007). According to Fig. 6, considering variable density and yield stress will change elastic-plastic radial displacement significantly. Similar to Fig. 3, the effects of the exponent rate are presented in Fig. 6c and d.

Fig. 6
figure 6

Stress distribution for (a) nE = 1.3826, np = no = 0. Yields from inside and outside ( = 0.55, Ω = 1.307), (b) nE = np = no = 1.3826, = 0.55 and ep = 0.65(Ω = 1.1999) yields from inside, (c) nE = np = no = − 2, and (d) nE = np = no = 2

Plastic strains for the case of constant yield and density (Fig. 7a) and equal exponent rate for density, yield, and elastic modulus (Fig. 7b) investigated. The results are quite different in two cases for the equivalent non-dimensional loading parameter. Non-dimensional loading parameter (Ω) calculated as 1.307 according to Fig. 5 for fully plastic behavior (Fig. 7a) and 1.14 for a yield initiated case (Fig. 7b). The analogy between two cases again reveals that the yield initiates from inside if density, yield, and elastic modulus are variable. To have a better comparison, radial displacements for two cases (Fig. 7c and d) presented.

Fig. 7
figure 7

Stress distribution for (a) nE = 1.3826, np = no = 0. Yields from inside and outside ( = 0.55, Ω = 1.307), (b) nE = np = no = 1.3826, = 0.55 and ep = 0.65 (Ω = 1.19999) yields from inside, (c) nE = np = no = − 2, and (d) nE = np = no = 2

Finally, to get an analogy in different cases, the effect of variable density on radial displacement for a/b = 0.5 graphed in Fig. 8. It shows that for the same elastic modulus and yield stress exponent, lower density exponent reveals lower rotor displacement. Which is expected based on the authors’ experience in different gas turbine rotor design and maintenance.

Fig. 8
figure 8

Effect of variable density on radial displacement for a/b = 0.5

Conclusion

In the present article, elastic-plastic behavior of a rotating shaft made of FGM under high centrifugal forces is investigated for the first time. Modulus of elasticity, density, and yield stress is assumed to have a power-law function of the cylindrical coordinate system and all parameters concluded in an analytical model which is an improvement regarding previous jobs. The analytical equations derived based on different studies, and non-dimensional parameters defined to create comprehensive and analogical outcomes. The results are compared and validated with homogenous materials and previously published articles which considered modulus of elasticity as the only variable of the material (Akis & Eraslan, 2007).

According to the presented research, the shaft’s deformations and strength have a great dependency on the material property definition. The results show that neglecting the variety of density and yield stress causes a considerable difference in stress and strain and yielding initiation behavior may change from the inner surface of the shaft to its outer. It is essential to put great care to determine material properties for high-speed components such as hollow shaft to prevent design flaws in such sensible parts of the machine.

Due to the experience of the authors in gas turbine design industries, there is a need to have some robust formulas to check the yield point of hollow shafts during turbine maintenance. Different Rolls Royce gas turbine series such as Trent and AVON, Siemens SGT 800, and many other midsized turbines are using hollow shafts in their compressor and turbine parts. This model will help to control the yield start point for a hollow shaft measured during maintenance. This paper will pave a reliable way to design many high-speed rotary components. Although this research is about FGM materials, it could be outstretched for orthogonal and non-isotropic materials as well. The results will help designers to get a better perception of hollow shafts possible weaknesses and failures to design more efficient rotary machines.

Availability of data and materials

Data are available by the request.

Abbreviations

FGM:

Functionally graded materials

References

Download references

Acknowledgements

The paper is theoretical work of the authors.

Funding

Not applicable

Author information

Authors and Affiliations

Authors

Contributions

ST coordinated, validated, and formulated the data of the study. SA and MH contributed to literature review, calculations, and numerical analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Shams Torabnia.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Torabnia, S., Aghajani, S. & Hemati, M. An analytical investigation of elastic-plastic deformation of FGM hollow rotors under a high centrifugal effect. Int J Mech Mater Eng 14, 16 (2019). https://doi.org/10.1186/s40712-019-0112-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40712-019-0112-7

Keywords