shell-buckling

name	shell-buckling
description	Shell buckling analysis — cylindrical shells under axial compression (Euler, ECCS knockdown factors), external pressure buckling (ASME UG-28, Windenburg-Trilling, Windenburg chart), spherical shells under external pressure, combined loading (axial + pressure + torsion interaction), imperfection sensitivity (Koiter theory, initial geometric imperfections), stiffened shells (ring-stiffened, stringer-stiffened), and NASA SP-8007/8019 design monographs.
metadata	{"priority":7,"promptSignals":{"phrases":["shell buckling","cylindrical shell buckling","external pressure buckling","shell stability","buckling of shells","thin-walled shell"],"minScore":3}}

Shell Buckling Analysis — Complete Skill

Shell Buckling Fundamentals

Classical (Perfect Shell) Critical Loads

Cylindrical shell axial compression (classical theory): N_cr,class = E × t² / (R × √3(1-ν²)) [N/m; force per unit circumference] σ_cr,class = E × t / (R × √3(1-ν²)) [MPa] For ν = 0.3: σ_cr,class ≈ 0.605 × E × t/R [classical theoretical buckling stress]

Cylindrical shell torsion (classical): N_θz,cr = 0.272 × E × (t/R)^(5/4) × (R/L)^(1/2) [approximate; Timoshenko] Better: τ_cr = k_s × π² × E × (t/L)² / (12(1-ν²)) [k_s depends on L/R, R/t]

Cylindrical shell bending (classical): Same σ_cr as axial compression for classical; but knockdown factors differ

Spherical shell external pressure (classical): p_cr = 2E × (t/R)² / √3(1-ν²) ≈ 1.21 × E × (t/R)² [for ν = 0.3]

Imperfection Sensitivity and Knockdown Factors

Critical Problem with Classical Theory

Imperfection sensitivity: shell buckling is acutely sensitive to geometric imperfections (waviness, out-of-roundness, eccentricities) Classical theory overestimates buckling load by factor of 2–10 for thin shells Test data scatter: σ_cr_test / σ_cr_class = 0.2–0.8 (wide scatter for axially loaded cylinders)

Empirical Knockdown Factor (γ)

ECCS (European Convention for Constructional Steelwork): For cylindrical shells under axial compression: σ_cr,design = γ × σ_cr,class [γ = knockdown factor] γ from ECCS interaction with ω = L²/(R×t) and r = R/t parameter

NASA SP-8007 knockdown factor α: α = 1 - 0.901 × (1 - exp(-φ)) [φ = (1/16) × √(R/t); dimensionless] Design stress: σ_cr = α × σ_cr,class For R/t = 100: φ = 1/16 × √100 = 0.625; α = 1 - 0.901×(1-0.535) = 1 - 0.418 = 0.582

Minimum knockdown factor for different loading: Axial compression: α ≈ 0.2–0.6 (most sensitive; most scatter) External pressure: α ≈ 0.6–0.8 (less sensitive) Torsion: α ≈ 0.65–0.80

Cylindrical Shell — External Pressure (Windenburg-Trilling / ASME)

ASME UG-28 Method

Procedure for external pressure design:

1. Calculate L/Do and Do/t
2. Find A (strain) from ASME Section II Part D Figure G: A = 0.125 / [(L/Do)/(Do/t)^0.333] [approximate; read from chart for exact value] [L = unsupported length; Do = outer diameter; t = shell thickness]
3. Find B (allowable stress) from material-specific charts: B = 2AE/(3(1 - ν²)) if below yield; B = σ_y if above (plastic correction) Alternatively: B = S_y × (A - 0.002) × E / (S_y + 200 × A × E) [simplified Johnson formula]
4. Allowable external pressure: P_a = B / (Do/t)
5. Require: P_design ≤ P_a

Example: Shell: Do = 600 mm; t = 8 mm; L = 2,000 mm; E = 200 GPa; σ_y = 250 MPa L/Do = 2,000/600 = 3.33; Do/t = 600/8 = 75 A = 0.125/(3.33/75^0.333) = 0.125/(3.33/4.22) = 0.125/0.789 = 0.158 → use chart value From ASME Figure G: A ≈ 0.003 (read directly; approximation above is simplified) B = 2×0.003×200,000/3 = 400 MPa → above σ_y → B = σ_y = 250 MPa (clamp at yield) P_a = 250/(600/8) = 250/75 = 3.33 MPa

Windenburg-Trilling Formula (Analytical)

Critical external pressure (unsupported cylinder): P_cr = 2.6 × E × (t/Do)^(5/2) / (L/Do - 0.45(t/Do)^(1/2)) [MPa; Windenburg 1934] Good approximation for medium-length cylinders (non-uniform pressure modes)

Short cylinder (Euler ring mode dominates): P_cr = 0.855 × E × (t/Do)^3 / (1 - ν²) [ring buckling under radial pressure]

Long cylinder (column mode): P_cr = 2E × (t/R)^3 / (3(1-ν²)) × 1/(L/R)^2 [Euler column analog]

Spherical Shells — External Pressure

Classical critical pressure: p_cr = 2E(t/R)² / √(3(1-ν²)) ≈ 1.21E(t/R)² [for ν = 0.3]

Knockdown factor for spheres: α_sphere = 0.124 (from NASA SP-8019; based on test data scatter) Design: p_allow = 0.124 × 2E(t/R)² / √3(1-ν²)

ASME UG-32 dished heads under external pressure: Use Do = 0.8 × crown radius × 2 = 1.6 R_crown as equivalent sphere (2:1 ellipsoidal heads) Apply UG-28 procedure with this equivalent diameter

Stiffened Shells

Ring-Stiffened Cylindrical Shell

Ring stiffeners prevent overall collapse; shell bays collapse between rings: Bay buckling: use unsupported length L = ring spacing for shell between rings Ring design: prevent overall ring buckling (tripping): I_ring_min = P_cr_overall × L × R³ × (n²-1) / E [I = second moment; n = number of lobes]

Optimal ring spacing: Minimize total weight: balance bay shell thickness t vs. ring weight Rule of thumb: L_ring_spacing ≈ 1.7 × √(R × t) (approximate for minimum weight)

Smeared stiffener theory (for closely spaced stiffeners): Replace stiffened shell with equivalent orthotropic shell (effective moduli D_xx, D_θθ, D_xθ) PANDA2, STAGS, ABAQUS can model this directly

Stringer-Stiffened Shell (Axial Load)

Stringer stringers increase axial buckling load: N_cr = (n² × π² × E_x_eff × t_eff × R²) / L² [n = number of lobes; E_x_eff = smeared modulus] For column collapse: use Johnson-Euler buckling with effective stringer area

Combined Loading Interaction

Axial + External Pressure Interaction

Conservative linear interaction (ECCS/API): (σ_x / σ_cr,x) + (P / P_cr,P) ≤ 1.0 [bilinear interaction; conservative]

Better interaction (Teng, 2004): (σ_x / σ_cr,x)^(1.2) + (P / P_cr,P)^(0.5) ≤ 1.0 [curved interaction; less conservative]

Axial + Bending + External Pressure: Check maximum compressive stress: σ_max = σ_axial + σ_bending Use combined N_x = (max compressive + P×R/2) × t → apply to axial buckling check Plus: external pressure alone check

Torsion + External Pressure: (τ / τ_cr)^2 + (P / P_cr)^(3/2) ≤ 1.0 [Galletly interaction formula]

Numerical Methods and Software

Classical/empirical: NASA SP-8007 (cylinders under compression), SP-8008 (cylinders under pressure), SP-8019 (spheres) FEM: ABAQUS BUCKLE step → linear eigenvalue; then NLGEOM with imperfection → nonlinear collapse Imperfection seeding: scale first linear buckling mode by δ_imp = 0.1t to 0.5t (sensitivity study) STAGS (Structural Analysis of General Shells): specialized for shell stability; used by NASA/aerospace

Standards and References

Standard	Scope
ASME Sec. VIII UG-28	External pressure vessel design
NASA SP-8007	Buckling of thin-walled circular cylinders
NASA SP-8019	Buckling of thin-walled doubly curved shells
ECCS Publication No. 125	Shell buckling recommendations
EN 1993-1-6 (Eurocode 3)	Strength and stability of shell structures
DNV-RP-C202	Buckling strength of shells (offshore)

Output

Provide: shell geometry (R [mm]; t [mm]; L [mm]; material: E [GPa]; σ_y [MPa]; ν), loading (axial N_x [kN/m]; external pressure P [MPa]; bending M [kN·m]; torsion T [kN·m]; or combinations), classical critical load (σ_cr,class or P_cr,class [MPa] from formula, with source), knockdown factor (γ or α; source: NASA SP-8007/ECCS/empirical; basis for selection), design critical load (σ_cr = α × σ_cr,class [MPa] or P_cr [MPa]), safety factor (load/resistance γ_R ≥ 2.0 for ECCS; ≥ 3.0 for external pressure ASME; verify), combined loading interaction (if multiple loads: interaction equation; result ≤ 1.0), stiffener design (if needed: ring spacing [mm]; I_ring [mm⁴]; stringer dimensions), FEM recommendation (ABAQUS BUCKLE + nonlinear verification with imperfection δ_imp = 0.1t), and applicable standard (ASME UG-28 for pressure vessels; NASA SP-8007/8019 for aerospace; ECCS/EN 1993-1-6 for civil/offshore).