Engineering › Materials ›

Rotating Discs and Cylinders

The Stresses and Strains generated in a rotating disc or cylinder

viewed 2548 times

View version details

Contents

Discs Of Uniform Thickness.
A Solid Disc
A Disc With A Central Hole
Long Cylinders
A Disc Of Uniform Strength
Page Comments

Discs Of Uniform Thickness.

For a "thin" disc it can be assumed that the Stress in the Axial direction is zero. Due to the rotation of the disc, Circumferential $\displaystyle f$ (Hoop) and radial $\displaystyle p$ Stresses will occur. Both these Stresses are Tensile. If the radial shift at a radius r is u, then the Strain Equations are:-

(1)

$E\times \frac{du}{dr}\;=\;p\;-\;\frac{f}{m}$

(2)

$And\;\;\;\;\;E\times \frac{u}{r}\;-\;\frac{f}{m}$

Differentiating equation (2) and equating it to equation (1) gives:-

(3)

$\left(f\;-\;p \right)\left(1\;-\;\frac{1}{m} \right)\;+\;r\times \frac{df}{dr}\;-\;\left(\frac{r}{m} \right)\left(\frac{dp}{dr} \right)\;=\;0$

If the angular velocity of rotation is $\displaystyle \omega$ and the density of the disc material is $\displaystyle \rho$ then the element shown in the diagram is subjected to a centrifugal force which:-

(4)

$\left(\frac{\rho r\,\delta \theta \times \delta r}{g} \right)\times r\;\omega ^2$

Which for unit thickness:-

(5)

$\left(\frac{\rho r\,\delta \theta \times \delta r}{g} \right)\times r\;\omega ^2$

(6)

$=\;\frac{\rho \,r^2\,\omega ^2\;\delta r\times \delta \theta }{g}$

The equilibrium equation in the radial direction is:-

(7)

$2f\,\delta r\,\sin\frac{1}{2}\,\delta \theta \;+\;\rho r\,\delta \theta \;-\;(p\;+\;\delta p)(r\;+\;\delta r)\delta \theta \;=\;\frac{\rho \,r^2\,\omega ^2\;\delta r\times \delta \theta }{g}$

In the limit this equation reduces to:-

(8)

$f\;-\;p\;-\;r\,\frac{dp}{dr}\;=\;\frac{\rho \,r^2\;\omega ^2}{g}$

From equation (8) obtain $\displaystyle f\;-\;p$ and substitute into equation (3)

(9)

$\left(r\,\frac{dp}{dr}\;+\;\frac{\rho \,r^2\,\omega ^2}{g} \right)\left(1\;-\;\frac{1}{m} \right)\;+\;r\times \frac{df}{dr}\;-\;\left(\frac{r}{m} \right)\frac{dp}{dr}\;=\;0$

Rearranging

(10)

$\frac{df}{dr}\;+\;\frac{dp}{dr}\;=\;-\;\left(\frac{\rho \.r\,\omega ^2}{g} \right)\left(1\;+\;\frac{1}{m} \right)$

Integrating:-

(11)

$f\;+\;p\;=\;-\;\left(\frac{\rho \,r^2\,\omega ^2}{2g} \right)\left(1\;+\;\frac{1}{m} \right)\;+\;2A$

Subtracting equation (8)

(12)

$2p\;+\;r\,\frac{dp}{dr}=\;-\;\left(\frac{\rho \,r^2\,\omega ^2}{2g} \right)\left(3\;+\;\frac{1}{m} \right)\;+\;2A$

(13)

$\therefore\;\;\;\;\;\frac{1}{r}\;\frac{d(p\,r)^2}{dr}\;=\;-\;\frac{\rho \,r^2\,\omega ^2(3\;+\;1/m)}{2g}\;+\;2A$

Integrating:-

(14)

$p\,r^2\;=\;-\;\frac{\rho \,r^4\,\omega ^2(3\;+\;1/m)}{8g}\;+\;2A\,r^2\;-\;B$

(15)

$Or\;\;\;\;\;\;\mathbf{p\;=\;A\;-\;\frac{B}{r^2}\;-\;\left(3\;+\;\frac{1}{m} \right)\left(\frac {\rho \,r^2\;\omega ^2}{8g} \right)}$

Substituting from equation (11)

(16)

$\mathbf{f\;=\;A\;+\;\frac{B}{r}\;-\;\left(1\;+\;\frac{3}{m} \right)\left(\frac{\rho \;r\;\omega ^2}{8g} \right)}$

A Solid Disc

Since the Stresses can no be infinite at the centre of a solid Disc B must be zero.

If R is the outside radius of the disc then rewriting equations (14) and (15) then:-

(17)

$p\;=\;0\;=\;A\;-\;\left(3\;+\;\frac{1}{m} \right)\left(\frac{\rho \,R^2\,\omega ^2}{8g} \right)$

From Which:-

(18)

$p\;=\;\left(\frac{\rho \;\omega ^2}{8g} \right)\left(3\;+\;\frac{1}{m} \right)\left(R^2\;-\;r^2 \right)$

(19)

$And\;\;\;\;\;f\;=\;\left(\frac{\rho \;\omega ^2}{8g} \right)\left[ \left(3\;+\;\frac{1}{m} \right)R^2\:-\:\left(1\;+\;\frac{3}{m} \right)r^2 \right]$

At the centre r = 0 and so:-

(20)

$p\;=\;f\;=\;\left(3\;+\;\frac{1}{m} \right)\left(\frac{\rho \;\omega ^2\,R^2}{8g} \right)$

This is the Maximum Stress.

At the outside :-

(21)

$f\;=\;\left(1\;-\;\frac{1}{m} \right)\left(\frac{\rho \,\omega ^2\,r^2}{g} \right)$

For a value of $\displaystyle \frac{1}{m}\;=\;0.3$

(22)

$\hat{f}\;=\;\left(\frac{3.3}{8} \right)\left(\frac{\rho \,\omega ^2\,R^2}{g} \right)\;=\;0.41\times \frac{\rho \,\omega ^2\,R^2}{g}\;\;\;\;(At\;the\;centre)$

And at the outside:-

(23)

$f\;=\;\frac{0.7}{4}\times \frac{\rho \,\omega ^2\,R^2}{g}\;=\;0.425\;\hat{f}$

The variations of the Hoop and Radial Stresses with Radius are shown on the diagram.

A Disc With A Central Hole

The Radial Stress is Zero at both the inner and outer radii. If the value of these is $\displaystyle R_1\;\;and\;\;R_2$ then using equation (14)

(24)

$0\;=\;A\;-=\;\frac{B}{R_1^2}\;-\;\left(3\;+\;\frac{1}{m} \right)\left(\frac{\rho \;R_1^2\,\omega ^2}{ 8g }\right)$

(25)

$0\;=\;A\;-=\;\frac{B}{R_2^2}\;-\;\left(3\;+\;\frac{1}{m} \right)\left(\frac{\rho \;R_2^2\,\omega ^2}{ 8g }\right)$

Solving the above two equations for A and B:-

(26)

$B\;=\;\left(3\;+\;\frac{1}{m} \right)\left(\frac{\rho \,\omega ^2}{8g} \right)\left(R_1^2\,R_2^2 \right)$

(27)

$A\;=\;\left(3\;+\;\frac{1}{m} \right)\left(\frac{\rho \,\omega ^2}{8g} \right)\left(R_1^2\;+\;R_2^2 \right)$

(28)

$Thus\;\;\;\;\;\;p\;=\;\left(3\;+\;\frac{1}{m} \right)\left(\frac{\rho \,\omega ^2}{8g} \right)\left(R_1^2\;+\;R_2^2\;-\;\frac{R_1^2\,R_2^2}{r^2}\;-\;r^2\right)$

Substituting from equation(15)

(29)

$And\;\;\;\;\;f\;=\;\left(\frac{\rho \;\omega ^2}{8g} \right)\left[\left(3\;+\;\frac{1}{m} \right)\left(R_1^2\;+\;R_2^2\;+\;\frac{R_1^2\,R_2^2}{r^2} \right)\;-\;\left(1\;+\;\frac{3}{m} \right)r^2 \right]$

p is a maximum when $\displaystyle r\;=\;\sqrt{R_1\,R_2}$

(30)

$\hat{p}\;=\;\left(3\;+\;\frac{1}{m} \right)\left(\frac{\rho \;\omega ^2}{8g} \right)\left(R_2\;-\;R_1 \right)^2$

(31)

$\hat{f}\;=\;\left(\frac{\rho \;\omega ^2}{4g} \right)\left[\left(1\;-\;\frac{1}{m} \right)R_1^2\;+\;\left( 3\;+\;\frac{1}{m} \right)R_2^2 \right]$

Note that if $\displaystyle R_1$ is very small $\displaystyle \hat{f}\rightarrow \left(3\;+\;\frac{1 } {m} \right)\left(\frac{\rho \;\omega ^2\;R_2^2}{4g} \right)$ which is tweiew the value which would be found in a solid disc.

At the outside

(32)

$f\;=\;\left(\frac{\rho \;\psi ^2}{4g} \right)\left[\left(3\;+\;\frac{1}{m} \right)R_1^2\;+\;\left(1\;-\;\frac{1}{m} \right)R_2^2 \right]$

(33)

$If\;\;\;R_1\rightarrow R_2\;=\;R\;\;\;\;Then\;\;\;\;\hat{f}\rightarrow \frac{\rho \;\omega ^2\;R^2}{g}$

As in the case of a thin rotating cylinder. The variations in Stress are shown on the diagram.

Example 1:

A thin uniform disc of 10 in. diameter with a central hole of 2 in. , runs at 10,000 r.p.m. Calculate the maximum Principal Stress and the maximum Shear Stress in the disc. Take $\displaystyle \frac{1}{m}\;\;\;=\;\;0.3\;\;\;\;\;and\;\;\;\;\;Density\;=\;0.28\;lb.in.^{-3}$

The maximum Principal Stress will be at the inside and is given by equation (31)

(34)

$\hat{f}\;=\;\frac{0.28}{4\times 32.2\times 12}\left(\frac{10,000\times 2\pi }{60} \right)^2\;\left(0.7\times 1^2\;+\;3.3\times 5^2 \right)$

(35)

$=\;16,500\;lb.in^{-2}$

The maximum Shearing Stress at any radius is given by:-

(36)

$\frac{1}{2}(f\;-\;p)\;=\;\left(\frac{\rho \;\omega ^2}{8g} \right)\left[\left(3\;+\;\frac{1}{m} \right)\frac{R_1^2\;R_2^2}{r^2} \;+\;\left(1\;-\;\frac{1}{m} \right)r^2\right]$

It can be seen from the diagram that the greatest Stress difference occurs at

(37)

$r\;=\;R_1$

and the Maximum Shearing Stress is:-

(38)

$\frac{0.28}{8\times 32.2\times 12}\left(\frac{10.000\times 2\pi }{60} \right)^2\left(3.3\times\frac{1^2\times5^2}{1^2}\;+\;0.7\times 1^2 \right)\;=\;8250\;lb.in.^{-2}$

Long Cylinders

Assume that the Longitudinal Stress is $\displaystyle f_1$ and that the Longitudinal Strain e is constant. i.e. the cross sections remain plane which must be true away from the ends. The analysis is similar that used for a rotating disc and the Strain Equations are:-

(39)

$E\,e\;=\;f_1\;-\;\left(\frac{1}{m} \right)\left(f\;+\;p \right)$

(40)

$E\;\frac{du}{dr}\;=\;p\;-\;\left(\frac{1}{m}\right)\left(f\;+\;f_1$

(41)

$E\;\frac{u}{r}\;=\;f\;-\;\left(\frac{1}{m}\right)\left(p\;+\;f_1 \right)$

Using Equations (40) and (41) eliminate $\displaystyle \frac{du}{dr}$

From Equation (41)

(42)

$E\,.\,\frac{du}{dr}\;=\;f\;-\;\left(\frac{1}{m} \right)\left(p\;+\;f_1 \right)\;+\;\left[\frac{df}{dr}\;-\;\left(\frac{1}{m} \right)\left(\frac{dp}{dr} \right)\;-\;\left(\frac{1}{m}\right)\left(\frac{df_1}{dr} \right) \right]$

From Equation (40)

(43)

$=\;p\;-\;\left(\frac{1}{m} \right)\left(f\;+\;f_1 \right)$

(44)

$\therefore\;\;\;\;\;\;(f\;-\;p)\left(1\;+\;\frac{1}{m} \right)\;+\;r\times \frac{df}{dr}\;-\;\left(\frac{r}{m} \right)\left(\frac{dp}{dr} \right)\;-\;\left(\frac{r}{m} \right)\left( \right)\;=\;0$

Substituting $\displaystyle \frac{df_1}{dr}\;=\;\left(\frac{1}{m} \right)\left(\frac{df}{dr}\;+\;\frac{dp}{dr} \right)$ which was obtained from equation (39)

(45)

$(f\;-\;p)\left(1\;+\;\frac{1}{m} \right)\;+\;r\times \frac{df}{dr}\left(1\;-\;\frac{1}{m^2} \right)\;-\;\left(\frac{r}{m} \right)\left(1\;+\;\frac{1}{m}\right)\left(\frac{dp}{dr} \right)\;=\;0$

(46)

$Or\;\;\;\;\;f\;-\;p\;-\;r\times \left(1\;-\;\frac{1}{m} \right)\left(\frac{dp}{dr} \right)\;=\;0$

The equilibrium equation is as equation (8) and subtracting equation (46) we get.

(47)

$-\;r\left(1\;-\;\frac{1}{m} \right)\left(\frac{df}{dr} \right)\;-\;r\left(1\;-\;\frac{1}{m} \right)\left(\frac{dp}{fdr} \right)\;=\;\frac{\rho \times r\times \omega ^2}{g}$

(48)

$Or\;\;\;\;\;\frac{df}{dr}\;+\;\frac{dp}{dr}\;=\;\frac{\rho \times r\times \omega ^2}{\left(1\;-\;\frac{1}{m} \right)g}$

Integrating;-

(49)

$Or\;\;\;\;\;f\;+\;p\;=\;\frac{\rho \times r^2\times \omega ^2}{2\left(1\;-\;\frac{1}{m} \right)g}\;+\;2A$

This can be compared to equation (11) which is for a thin disc. It can be seen that the result for a long cylinder can be obtained from those for a thin disc by writing $\displaystyle \frac{1}{1\;-\;\frac{1}{m}}\;instead\;of\;1\;+\;\frac{1}{m}$

Command [\B] not found.

{Solid Cylinder}

The maximum Stress occcurs at the centre where p and f are equal

(50)

$\hat{f}\;=\;\frac{3\;-\;\frac{2}{m}}{1\;-\;\frac{1}{m}}\times \frac{\rho \times\omega ^2\times R^2}{8g}$

(51)

$If\;\;\;\;\;\frac{1}{m}\;=\;0.3$

(52)

$Then\;\;\;\;\;\hat{f}\;=\;\frac{2.4}{5.6}\times \frac{\rho \times\omega ^2\times R^2}{8g}\;=\;0.43\;\frac{\rho \times \omega ^2\times R^2}{8g}$

( This compares with $\displaystyle 0.41\times\frac{\rho \times \omega ^2\times R^2}{8g}$ for a solid disc)

Hollow Cylinder

$\displaystyle \hat{p}\;=\;\left(\frac{3\;-\;\frac{2}{m}}{1\;-\;\frac{1}{m}} \right)\times \frac{\rho \times \omega ^2}{8g}\times\left(R_2\;-\;R_1 \right)^2$

$\displaystyle \hat{f}\;=\;\frac{\rho \times \omega ^2}{4g\left(1\;-\;\frac{1}{m} \right)}\left[\left(1\;-\;\frac{2}{m}\right)R_1^2\;+\;\left(3\;-\;\frac{2}{m} \right) R_2^2 \right]$

These values do not differ greatly from those for a thin disc.

A Disc Of Uniform Strength

Consider the condition of equal stress at all radii. i.e. $\displaystyle p\;=\;f\;=\;Constant$

Le $\displaystyle t$ be the thickness of the disc at a radius $\displaystyle r$ and

assume that the thickness will be $\displaystyle t\;+\;\delta t$ at radius $\displaystyle r\;+\;\delta r$

The mass of the element will be approximately

(53)

$\rho \times r\times \delta \theta \times \delta r\times \frac{t}{g}$

And the centrifugal Force will:-

(54)

$\rho \times r^2\times\omega ^2\times \delta \theta \times \delta r\times \frac{t}{g}$

Hence the equilibrium Equation is:-

(55)

$2\;f\;\delta r\times \sin\frac{1}{2}\delta \theta \;+\;f\;r\;\delta \theta \times t\;=\;f(r\;+\;\delta r)\;\delta \theta \;(t\;+\;\delta t)\;+\;\frac{\rho \;r^2\;\psi ^2\;t\;dr}{g}$

Which in the limit is:-

(56)

$f\;t\times dr\;=\;f\;r\times dt\;+\;f\;t\times dr\;+\;\frac{\rho\;\omega^2\;t\times dr}{g}$

Integrating:-

(57)

$\lnt\;=\;-\;\frac{\rho \;r^2\;\omega ^2}{2\;f\;g}\;+\;Constant$

(58)

$Or\;\;\;\;\;\;t\;=\;A\;e^{-\;\frac{p\;r^2\;\omega ^2}{2\;f\;g}}$

(59)

$Or\;\;\;\;\;\;t\;=\;t_0\;e^{-\;\frac{p\;r^2\;\omega ^2}{2\;f\;g}}$

Example 2:

A Turbine rotor Disc is 24 in. Diameter at the blade ring and is keyed to a 2 in diameter shaft. If the minimum thickness is 3/8 in. what should be the thickness at the shaft for a uniform stress of 30,000 lb./sq.in. aqt 10,000 r.p.m.? Density of material = 0.28 lb./cu.in.

(60)

$T\;=\;A\;e^{-\frac{\rho \,r^2\,\omega ^2}{2f\;g}}$

(61)

$At\;\;\;\;\;\;\;r\;=\;12\;in.$

(62)

$t\;=\;\frac{3}{8}\;=\;A\;e^{-\rho \,\omega ^2\times\frac{144}{2f\,g}}$

(63)

$At\;\;\;\;\;\;\;r\;=\;1\;in.$

(64)

$t\;=\;A\;e^{-\frac{\rho \,\omega ^2}{2f\,g}}\;=\;\frac{3}{8}\;e^{\rho \,\omega ^2\times \frac{143}{2f\,g}}$

(65)

$Where\;\;\;\;\;\rho \;\omega ^2\times \frac{143}{2f\;g}\;=\;0.28\left(\frac{10,000\pi }{30} \right)^2\times \frac{143}{2}\times30,000\times 32.2\times 12$