Skip to main content

Thick cylinder stresses and failure criteria

Short answer

By Lamé equations: inner surface (r=80 mm) σθ=104 MPa, σr=−40 MPa; mid-wall (r=100 mm) σθ=78.1 MPa, σr=−14.1 MPa; outer surface (r=120 mm) σθ=64 MPa, σr=0; axial stress 32 MPa. Factor of safety: 2.43 (Tresca), 2.81 (von Mises). Hoop strain 5.32×10⁻⁴, axial strain 6.4×10⁻⁵, volumetric strain 1.92×10⁻⁴.

A thick steel cylinder has an inside diameter of 160 mm and an outside diameter of 240 mm. It is subjected to an internal pressure of 40 MPa, the outer surface is free. The ends are closed. Take E=200E = 200 GPa, ν=0.3\nu = 0.3, yield stress σy=350\sigma_y = 350 MPa. (a) Find the hoop and radial stresses at the inner and outer surfaces and at mid-wall radius of 100 mm, and the axial stress. (b) Find the factor of safety against yielding by the maximum shear stress (Tresca) theory and the von Mises theory. (c) Find the hoop strain, axial strain and volumetric strain at the inner surface.

Data

a=80a = 80 mm, b=120b = 120 mm, pi=40p_i = 40 MPa, po=0p_o = 0.

pia2b2−a2=40×640014400−6400=32 MPa\frac{p_ia^2}{b^2 - a^2} = \frac{40\times6400}{14400 - 6400} = 32\ \text{MPa}

(a) Lame's stresses

σθ=32(1+14400r2),σr=32(1−14400r2)\sigma_\theta = 32\left(1 + \frac{14400}{r^2}\right),\quad \sigma_r = 32\left(1 - \frac{14400}{r^2}\right)
Radius rr (mm)b2/r2b^2/r^2σθ\sigma_\theta (MPa)σr\sigma_r (MPa)
80 (inner)2.25104.00-40.00
1001.4478.08-14.08
120 (outer)1.0064.000.00

Axial stress for closed ends: σz=pia2/(b2−a2)=32.0\sigma_z = p_ia^2/(b^2 - a^2) = 32.0 MPa.

Check: σθ+σr=2A=64\sigma_\theta + \sigma_r = 2A = 64 MPa is the same at every radius, for example 104−40=64104 - 40 = 64 at the inner surface and 64+0=6464 + 0 = 64 at the outer surface. ✓\checkmark

At the inner surface: σθ=104\sigma_\theta = 104, σr=−40\sigma_r = -40, σz=32\sigma_z = 32 MPa.

(b) Failure criteria

Maximum shear stress (Tresca): σ1=σθ=104\sigma_1 = \sigma_\theta = 104, σ3=σr=−40\sigma_3 = \sigma_r = -40 MPa (the axial stress, 32 MPa, lies between them).

σ1−σ3=104−(−40)=144 MPa,τmax=72 MPa\sigma_1 - \sigma_3 = 104 - (-40) = 144\ \text{MPa},\qquad \tau_{max} = 72\ \text{MPa} FS=σyσ1−σ3=350144=2.43FS = \frac{\sigma_y}{\sigma_1 - \sigma_3} = \frac{350}{144} = 2.43

von Mises (distortion energy):

σvm=12[(σθ−σr)2+(σr−σz)2+(σz−σθ)2]=12[1442+722+722]=124.7 MPa\begin{aligned} \sigma_{vm} &= \sqrt{\tfrac12\left[(\sigma_\theta - \sigma_r)^2 + (\sigma_r - \sigma_z)^2 + (\sigma_z - \sigma_\theta)^2\right]} \\ &= \sqrt{\tfrac12\left[144^2 + 72^2 + 72^2\right]} = 124.7\ \text{MPa} \end{aligned} FS=350124.7=2.81FS = \frac{350}{124.7} = 2.81

Tresca is more conservative, which is why it is commonly used for thick cylinders.

(c) Strains at the inner surface

εθ=1E[σθ−ν(σr+σz)]=104−0.3(−40+32)200×103=5.32×10−4εz=1E[σz−ν(σθ+σr)]=32−0.3(64)200×103=6.4×10−5εr=1E[σr−ν(σθ+σz)]=−40−0.3(136)200×103=−4.04×10−4\begin{aligned} \varepsilon_\theta &= \frac{1}{E}\left[\sigma_\theta - \nu(\sigma_r + \sigma_z)\right] = \frac{104 - 0.3(-40 + 32)}{200\times10^3} = 5.32\times10^{-4} \\ \varepsilon_z &= \frac{1}{E}\left[\sigma_z - \nu(\sigma_\theta + \sigma_r)\right] = \frac{32 - 0.3(64)}{200\times10^3} = 6.4\times10^{-5} \\ \varepsilon_r &= \frac{1}{E}\left[\sigma_r - \nu(\sigma_\theta + \sigma_z)\right] = \frac{-40 - 0.3(136)}{200\times10^3} = -4.04\times10^{-4} \end{aligned}

Volumetric strain:

εv=εθ+εr+εz=1−2νE(σθ+σr+σz)=0.4×96200×103=1.92×10−4\varepsilon_v = \varepsilon_\theta + \varepsilon_r + \varepsilon_z = \frac{1-2\nu}{E}(\sigma_\theta + \sigma_r + \sigma_z) = \frac{0.4\times96}{200\times10^3} = 1.92\times10^{-4}

The inner radius increases by u=εθ a=5.32×10−4×80=0.0426u = \varepsilon_\theta\,a = 5.32\times10^{-4}\times80 = 0.0426 mm.

Answer: (a) At r=80r = 80: σθ=104\sigma_\theta = 104 MPa, σr=−40\sigma_r = -40 MPa; at r=100r = 100: σθ=78.1\sigma_\theta = 78.1 MPa, σr=−14.1\sigma_r = -14.1 MPa; at r=120r = 120: σθ=64\sigma_\theta = 64 MPa, σr=0\sigma_r = 0; σz=32\sigma_z = 32 MPa. (b) FS =2.43= 2.43 (Tresca), 2.812.81 (von Mises). (c) εθ=5.32×10−4\varepsilon_\theta = 5.32\times10^{-4}, εz=6.4×10−5\varepsilon_z = 6.4\times10^{-5}, εv=1.92×10−4\varepsilon_v = 1.92\times10^{-4}.

✓ Answer checked by Claude Sonnet 5.5

All Mechanics of Solids chapter 8 questions →