Introduction..... Symbols..... Assumptions..... Basic Equations for a Rotating Disk..... Basic Equations for a Rotating Cylinder.....
Introduction These notes relate to the stresses and strains existing in thick walled cylinders when rotated at speed they are generally applicable to design of flywheels. The primary assumption is that the cylinders are not subject to internal or external prssure. A basic review of solid disks, rings and cylinders is carried out. Members of excelcalcs.com can upload a spreadsheet including all of the derived equations below at ExcelCalcs.com calculation Stresses in Rotating Disks & Rings ... Symbols / Units Tensile stresses are considered positive and compressive stresses are negative.
Initial Assumptions The following relationships are assumed for
the strains ε1,ε2,
ε3 associated with the stress σ 1,
σ 2 and σ 3.
ε1 = σ 1 /E - υσ 2 /E - υσ 3 /E Thick Disk basics Consider a "disk"/ "thin ring" subject internal stresses resulting from the inertial forces as a result of its rotational speed. Under the action of the inertial forces only, the three principal stress will be σ r tensile radial stress, σ t tensile tangential stress and σ a and axial stress which is generally also tensile. The stress conditions occur throughout the section and vary primarily relative to the radius r. It is assumed that the axial stress σ a is constant along the length of the section and because the disk is thin compared to its diameter the axial stress throughout the section is assumed to be zero.It is also assumed that the internal pressure (P1 ) and the external pressure ( P2 ) = 0. Consider a microscopically small area under stress as shown. u is the radial displacement at radius r . The circumferential (Hoop) strain due to the internal pressure is At the outer radius of the small section area (r + δr ) the radius will increase to (u + δ ). The resulting radial strain as δr -> 0 is Referring to the stress/strain relationships as stated above. The following equations are derived.
Basis of equations...
Eq. 1)......E ε a =0 - υσ t - υσ r Eq. 2)......E.ε t = E.u/r = σ t - υσ r Eq. 3)......E.εr = E.du/dr = σ r - υσ t - Multiplying 2) x r Eu = r ( σ t -- υσ r ) Differentiating Edu/dr = σ t - υσ r +
r. [ dσ t /dr - υ.( dσ r / dr ) ]
= σ r - υσ t .. from 3 above )
Eq. 4)........(σ t - σ r ). ( 1 + υ )
+ r.(dσ t/ dr ) ) - υ.r.(dσ r / dr) = 0 Now considering the radial equilibrium of the element of the section. Forces based on unit length of cylinder 2.σ t. δ r.sin(1/2.δθ) +
σrδθ - (σ r +
δσ r) (r + σ r )θδθ =
ρr2 ω2δr δθ In the limit this reduces to Eq 5 )..... σ t - σ r - r dσ r / d r = ρr2 ω2 Substitute for σ t - σ r in Equation 4 results in ( r dσ r / d r + ρr2 ω2 ). ( 1 + υ ) + r.(dσ t/ dr ) ) - υ.r.(dσ r / dr) = 0 Therefore dσ t/dr + dσ r/dr = - ρrω 2(1+υ) Integrating Eq 6 )..... σ t + σ r = - ρ r2 ω 2(1+υ)/2 + 2A Subtract equation 5... 2.σ r + r.dσ r/dr = - ρr2ω 2(3+υ)/2 + 2A This is the same as (1/r).d (σ r.r2)dr = - ρ r2 ω 2(3+υ)/2 + 2A Integrating σ r.r2 = - (ρ r4 ω 2(3+υ)/8 + Ar2 + B Eq 7 ... σ r = A + B/r2 - (3 + υ)ρ r2 ω 2/8 Eq 8 ......combining Eqn. 7 with Eqn. 6 σ t = A - B/r2 - ( 1+ 3υ)ρ r2 ω 2/8 1) Solid DiskAt the centre the B/r2 term implies infinite stresses which are clearly not credible and therefore B must equal 0. At r = R2 on the outside edge of the disk. the radial stress is equal to the surface stress which is equal to 0. Therefore at R2 σ r= 0 = A - (3 + υ)ρ R2 2 ω 2/8 σ t = (ρ ω 2/8)[ (3 + υ)R2 2 - (1 + 3υ) r2] σ r = (ρ ω 2/8)[ (3 + υ)( R2 2 - r2)] The maximum stress is at the centre σ t_max = σ r_max = (ρ ω 2/8).(3 + υ)R2 2 At the outside edge r = R2 and at the hole radius r = R1 the radial surface stress is asssumed to be 0 Therefore
Solving
2) Cylinders The primary difference between a long rotating cylinder and a thin one is that it the axial stress σ a is not equal to 0. The assumption in this case is that the longitudonal strain ε a is constant: cross sections remain plane.
Basis of equations... Refer to introductory notes at the top of webpage Multiplying 2) x r Eu = r ( σ t - υσ a - υσ r ) differentiating Edu/dr = σ t - υσ a - υσ r +
r. [ dσ t /dr - υ.( dσ a / dr ) - υ.( dσ r / dr ) ]
= σ r - υσ t - υσ a ..( from 3 above )
Eq. 4)........(σ t - σ r ). ( 1 + υ ) + r.(dσ t/ dr ) - υ.r.(dσ a / dr ) - υ.r.(dσ r / dr) = 0 substituting for dσ a/dr in Eq 4) (σ t - σ r ). ( 1 + υ ) + r.(dσ t/ dr ) - υ.r.(υ (dσ t/dr + dσ r /dr ) - υ.r.(dσ r / dr) = 0 (σ t - σ r ). ( 1 + υ )
+ r.(1-υ2)(dσ t/ dr )
- υ.r.(1 + υ ) (dσ r/dr)
= 0 Now considering the radial equilibrium of the element of the section as shown in the notes above Eq. 5 results Eq 5........ σ t - σ r - r dσ r / d r = ρr2 ω2 Substitute for (σ t - σ r )
r.(1-υ)(dσ t/ dr )
+ (1 - υ).r (dσ r/dr)
= - ρr2 ω2 (dσ t/ dr )
+ (dσ r/dr)
= - ρ r ω2 /(1-υ) Integrating Eq.6)... σt + σr = -ρr2ω2 /2(1-υ) + 2A This is similar to the equation 6 for the Rotating Disk analysis completed above. In fact the rotating disks equation can apply for the long cylinder if (1 + υ ) in the disk equations are replace by 1 / (1 - υ ) . Or if υ is replaced by (υ / (1- υ).. Now eliminating σ t by substituting into Eq 5) Integrating Eq.7).... σr = - ρr2ω2(3- 2υ ) /8(1-υ)
+ A - B / r2 Now substituting for σr in Equation 7 to obtain σt Eq.8).... σt = - ρr2ω2(1 + 2υ ) /8(1-υ)
+ A + B / r2 Solving for A and B for a solid cylinder At the centre of the cylinder R1 = 0 the stresses cannot be infinite so B is cleraly equal to 0. B = 0 At the outside diameter of the cylinder r = R2 . σ r = 0 = σr = - ρr2ω2(3- 2υ ) /8(1-υ)
+ A Eq.9).... .σr = ρω2(3- 2υ ) /8(1-υ ) (R2 2 - r 2) Eq.10).... σt = ρω2 /8(1-υ) [(3- 2υ )R22 -
(1 + 2υ ) r2 ]
The maximum radial stress and tangential stress are equal at r = 0 = Solving for A and B for a hollow cylinder At r = R1 and at r = R2 the radial stress σr = 0 Solving for A 0 = σr
= - ρR12ω2(3- 2υ ) /8(1-υ)
+ A -
R22.ρ ω2(3- 2υ ) /8(1-υ) Resulting in.. σr = ρω2(3- 2υ ) /8(1-υ ) [ - r2 + (R22 + R12) - (R12.R22)/r2 ] σr_max = ρω2(3- 2υ ) /8(1-υ ) ((R22 - R12)... Is located at at r = Sqrt (R1.R2 ) σt = ρω2 /8(1-υ ) [ - r2 (1 + 2υ ) + (3- 2υ ) {( R22 + R12) + (R12.R22)/r2} ] σt_max = ρω2 /4(1-υ ) [ (1 - 2υ )R12 + (3- 2υ )R22] ... Is located at r = R1. |