A pressure vessel, and when the thin-wall formula stops being true

The thin-wall hoop stress is one of the most useful approximations in engineering and one of the most quietly misused. It is exact in the limit of a membrane and good to a few per cent while the wall is thin; it understates the stress at the bore of a thick one, and the error grows with the wall.

This worksheet works a cylinder both ways and puts the two answers side by side, so the approximation is checked rather than assumed.

The vessel: a 600 mm bore receiver at 2 MPa, 12 mm wall, S355.

use steel

Four significant figures.

digits 4
ri=300⁢mm
twall=12⁢mm
p=2⁢MPa
ro=ri+twall=300⁢mm+12⁢mm=312⁢mm

Is it thin?

The usual rule is that the wall is thin below a tenth of the radius. It is a rule about the answer rather than about the geometry: what it really says is that the through-thickness variation of hoop stress is small enough to ignore.

ratio=twallri=12⁢mm300⁢mm=0.04
checkratio≤0.1=0.04≤0.1pass

Thin wall

Hoop stress from a cut along the axis, longitudinal from a cut across it. The factor of two between them is why a cylinder that fails from internal pressure splits along its length rather than round its circumference.

σhoop_thin=p·ritwall=2⁢MPa·300⁢mm12⁢mm=50⁢MPa
σlong_thin=p·ri2·twall=2⁢MPa·300⁢mm2·12⁢mm=25⁢MPa

Thick wall

Lamé's solution, evaluated at the bore where the hoop stress is largest. It makes no assumption about the wall being thin and is what the thin-wall formula approaches as the wall gets thinner.

σhoop_bore=p·(ro2+ri2)ro2−ri2=2⁢MPa·((312⁢mm)2+(300⁢mm)2)(312⁢mm)2−(300⁢mm)2=51.02⁢MPa

The radial stress at the bore is not zero — it is the pressure itself, in compression — and the thin-wall treatment drops it. At this wall it is 4% of the hoop stress; at a quarter-radius wall it is not negligible at all.

σradial_bore=−p=−2⁢MPa=-2⁢MPa

How much the approximation is out by, which is the number this worksheet exists to produce.

errorthin=σhoop_thin−σhoop_boreσhoop_bore=50⁢MPa−51.02⁢MPa51.02⁢MPa=-0.01998
errorthin=-1.998⁢%

The check

Against the thick-wall stress, because it is the true one, and on the distortion-energy criterion with all three principal stresses — hoop, radial and longitudinal. Dropping the radial term is the same approximation again and this worksheet has just finished measuring it.

σlong=p·ri2·twall=2⁢MPa·300⁢mm2·12⁢mm=25⁢MPa
s1=σhoop_bore=51.02⁢MPa
s2=σlong=25⁢MPa
s3=σradial_bore=-2⁢MPa
σvm=(s1−s2)2+(s2−s3)2+(s3−s1)22=(51.02⁢MPa−25⁢MPa)2+(25⁢MPa−-2⁢MPa)2+(-2⁢MPa−51.02⁢MPa)22=45.92⁢MPa
ndesign=1.5
σallow=FyS355ndesign=355⁢MPa1.5=236.7⁢MPa
checkσvm≤σallow=45.92⁢MPa≤236.7⁢MPapass

Where the approximation goes wrong

Hoop stress at the bore against wall thickness, both ways, on one plot. They agree while the wall is thin and part company as it thickens — and the thin-wall line is the lower of the two, which is the dangerous direction.

fn thin defined
fn thick defined

The legend would otherwise name the curves after the functions that drew them, which here are thin and thick. label says what they are.

axis x "Wall thickness"
axis y "Hoop stress"
label "Thin-wall", "Thick-wall (Lamé)"
plot⁡(thin,thick,5,120)
50 100 0 50 100 150 Wall thickness Hoop stress Thin-wall Thick-wall (Lamé)

What this leaves out

A plain cylinder, away from its ends. The heads, the nozzles and the joints between them are where a real vessel is governed, the stress concentrations there are not small, and a vessel is signed off against a code — ASME VIII, EN 13445 — that fixes the allowable stress, the joint efficiency and the corrosion allowance rather than leaving them to a factor of safety chosen on the page.