3.6.7 · D1 · Physics › Spacecraft Structures & Systems Engineering › Shell buckling — thin-walled cylinder under axial load
Ek rocket tank ek giant thin can hai jo lengthwise squeeze ho raha hai, aur woh metal ke crush hone se pehle hi fold ho jaata hai (buckles ) — kyunki ek thin curved wall pehle apni shape stability kho deti hai. Parent page par jo bhi hai woh sirf do tareekon ke beech ek tug-of-war ka bookkeeping hai jisse wall us fold se ladti hai: bending stiffness aur membrane hoop stiffness .
Is page par assume kiya gaya hai ki parent page ke notation ke baare mein tumhe kuch nahi pata. Hum har letter, har ratio aur har picture ko zero se build karte hain, ek aisi order mein jahan har idea apne pehle wale idea par lean karta hai.
Kisi bhi formula se pehle, cheez ko picture karo. s01 mein, magenta arrow follow karo (woh R hai, round top ki spoke), left edge par orange band (woh wall thickness t hai), aur right par navy double-arrow (woh height L hai).
Definition Teen shape numbers
R , t , L
R — radius : centre line se wall tak ki distance. s01 mein yeh centre dot se rim tak ki magenta arrow hai.
t — wall thickness : metal skin kitni gehri hai. s01 mein yeh thin orange band hai — metal ka woh ring jo tube ko cross mein kaatne par dikhega.
L — length : tube apni axis ke saath kitni oonchi khadi hai, s01 mein navy double-arrow.
Topic ko teeno chahiye kyunki buckling ek geometry failure hai — shape decide karta hai kab fold hoga, na ki sirf material.
Definition "Thin-walled" aur ratio
R / t
Ek cylinder thin-walled hota hai jab wall radius ke comparison mein tiny ho: t ≪ R . Hum "kitna thin" ek single number ==R / t == (radius divided by thickness) se measure karte hain. Ek soda can ka R / t kuch hundred hota hai; rockets ~100–500 ke aas-paas hote hain.
Picture: bada R / t matlab bahut round, bahut flimsy skin — dimple karna aasaan. Yeh ek ratio, na ki R ya t akele, wahi hai jis par buckling formulas dhyan dete hain.
Ratio kyun, raw sizes kyun nahi?
Ek can ka radius aur thickness dono double karo toh woh utna hi "tinny" lagega — same R / t . Nature ko metres se koi matlab nahi; use proportion se matlab hai. Isliye parent page par har formula geometry ko t / R (ya uske flip R / t ) ke andar chhupa deta hai.
Tank ko upar stack kiye weight aur thrust ki wajah se apni axis ke along squeeze kiya jaata hai.
P
P axis ke along push hai — ek force, newtons (N ) ya meganewtons (MN = 1 0 6 N ) mein measure hota hai. Imagine karo ek giant haath tube ke upar seedha neeche daba raha hai.
σ (sigma)
Greek letter σ ka matlab stress hai: force us area par spread out jo use carry karti hai.
σ = area force
Units: pascals, Pa = N / m 2 ; hum usually megapascals use karte hain, MPa = 1 0 6 Pa .
Picture: ek moti wall par same push zyada metal par "dilute" ho jaata hai, toh stress kam hota hai. Stress squeeze ki intensity hai, total push nahi.
P nahi, stress kyun chahiye?
Materials ek fixed intensity of loading par toot te hain, fixed total force par nahi. Ek mota straw us push ko survive kar leta hai jo ek thin wale ko same P par crush kar deta hai, kyunki moota wala use spread out kar deta hai. Stress compare karne ka fair tarika hai — isliye parent pehle ek critical stress σ cr derive karta hai, phir force par wapas multiply karta hai.
Definition Load-bearing area
A = 2 π R t
Sirf metal ki thin ring axial push carry karti hai — hollow andar wala nahi. s02 mein hum us ring ko kaat ke flat unroll karte hain: woh ek orange strip ban jaati hai jis ki length 2 π R (circumference, magenta double-arrow) aur depth t (wall, orange double-arrow) hai. Uska area hai
A = 2 π R t .
Isliye P = σ A = σ ( 2 π R t ) — stress wapas total force mein convert ho jaata hai.
Do numbers metal khud ko describe karte hain, shape se independent.
Definition Young's modulus
E
E material ki stiffness in stretch hai: ek given fraction tak stretch karne ke liye kitna stress lagana padega. Bada E = stretch karna mushkil (steel, E ≈ 200 GPa ); chhota E = springier (aluminium, E ≈ 70 GPa ). Yahan GPa = 1 0 9 Pa .
Picture: stress-vs-stretch line ki slope — steep slope, stiff material.
Definition Poisson's ratio
ν (nu)
Jab tum ek material ko ek taraf stretch karte ho, toh woh crosswise directions mein thin ho jaata hai. ν (Greek letter "nu") hai ki lengthwise stretch ke liye woh sideways kitna thin hota hai . Zyaatar metals ke liye ν ≈ 0.3 : 1% long stretch karo, 0.3% wide shrink.
Picture: s03 ise directly dikhata hai — dashed block original hai; violet block ko magenta arrow se long kheencha gaya hai aur saath hi (orange arrows) thinner bhi shrink hua hai. ν "thinner" part ko "longer" ke per unit measure karta hai.
ν kyun chahiye?
Kyunki ek plate wide hai, slender rod nahi. Jab ek wide plate ek taraf bend hoti hai, uski Poisson sideways-shrink use doosri taraf bhi curve karna chahti hai, aur yeh fight use stiffen karti hai. Yahi reason hai ki plate rigidity mein ek factor ( 1 − ν 2 ) hota hai jo ek simple beam mein nahi hota — dekho Plate bending and flexural rigidity D .
Parent page ka poora point yahi hai ki yeh cylinder obvious tarike se fail nahi karta.
σ Y
σ Y woh stress hai jis par material permanently crush/deform ho jaata hai — ek strength limit. Imagine karo metal khud dab raha hai. Dekho Yield vs stability failure modes .
Definition Buckling (stability failure)
Buckling ek shape failure hai: ek critical load par straight shape stable nahi rehta aur wall suddenly ek folded pattern mein jump kar jaati hai — jab metal abhi bhi σ Y se bahut door hota hai. Imagine karo ek empty can par step karna: metal theek hai, shape collapse ho jaati hai.
Intuition Woh ek distinction jo sab kuch drive karti hai
Yield = material overwhelm ho gaya (strength ).
Buckling = geometry unstable ho gayi (stiffness/stability ).
Thin cylinders ke liye buckling limit lower hoti hai, toh wahi jeetti hai. Parent ka Example 3 dikhata hai ki steel ~108 MPa par buckle karta hai jabki yield 250 MPa hai. Isliye tum kabhi bhi sirf σ Y ke liye design nahi kar sakte. Yeh Euler column buckling ki classic story ka mirror hai, jahan ek slender strut crush hone se pehle bend ho jaata hai.
Jab wall thodi radial amount w andar ya bahar dimple karti hai, do "springs" wapas fight karti hain.
Definition Radial displacement
w
w hai ki wall apni perfect round position se in ya out kitna poke karti hai , har point par. Imagine karo ek single dimple ki gehraai.
Definition Buckle wavenumber
k
Ek dimple pattern ek wave hai jo wall ke around aur along run karti hai. Wavenumber k count karta hai ki woh folds kitne tightly packed hain: bada k = kai short folds close together, chhota k = few long gentle folds. Ise "waves per metre" sochho — fold spacing (wavelength λ = 2 π / k ) ka inverse. Hume k chahiye kyunki wall khud decide karti hai kitna tightly fold karna hai, aur yahi choice buckling load set karti hai.
s04 mein, top strip (small k ) mein kuch long humps hain; bottom strip (large k ) mein kai tight humps hain. Same wall, do possible fold patterns.
Definition Flexural (bending) rigidity
D
D measure karta hai ki wall curved (bent) hone ka kitna resist karti hai.
D = 12 ( 1 − ν 2 ) E t 3
t 3 ka matlab hai ki thickness enormously matter karti hai: wall ke middle se door metal hi resist karta hai. ( 1 − ν 2 ) upar se Poisson coupling wala plate stiffening hai. Puri tarah Plate bending and flexural rigidity D mein unpack kiya gaya hai.
Picture: woh springiness jo dimple ki curvature ki sharpness ko resist karti hai.
D k 4 kyun scale karti hai
w 0 depth aur wavenumber k ke fold mein ek curvature hoti hai — woh kitna sharply bend karta hai. Ek wave w = w 0 sin ( k x ) ke liye, curvature second derivative hai d 2 w / d x 2 = − k 2 w 0 sin ( k x ) : har derivative ek factor of k le aata hai, toh curvature ∝ k 2 — tighter folds (bada k ) zyada sharply bend karte hain.
Wall mein stored bending energy D × ( curvature ) 2 ke proportional hai (ek bent spring energy store karta hai ∝ deflection2 ). Curvature ko square karne par ( k 2 ) 2 = k 4 milta hai. Isliye
U bend ∝ D k 4 w 0 2 .
Words mein: short tight folds bend karne mein bahut expensive hote hain kyunki har fold metal ko ek sharp corner se ghumaata hai. Yahi hai k 4 .
Definition Membrane (hoop) spring constant
K = R 2 E t
Wall ko w radially push karna circumference (hoop) ko stretch ya squeeze karta hai . Stored stretch-energy dimple se ladti hai, aur uski stiffness ek spring constant hai
K = R 2 E t .
Notice karo: yahan koi ( 1 − ν 2 ) nahi. Hoop stretch circumferential fibres ka uniaxial stretch hai — ek single direction ka lengthen hona — toh yeh plain axial stiffness E t use karta hai, two-directional plate stiffness nahi. Yahi honest reason hai ki ν D mein appear karta hai lekin K mein nahi: D ek bending (two-way plate) effect hai, K ek stretch (one-way membrane) effect hai. (Fuller treatments mein ek mild ν correction rehta hai, lekin is level par K = E t / R 2 woh clean number hai jo 2 D K mein enter karta hai.)
Picture: round hoop smaller ya bigger banane se resist karna — yeh sirf isliye exist karta hai kyunki wall curved hai (1/ R ). Ek flat plate mein koi hoop nahi, toh aisa koi term nahi. Yeh wahi circumferential stretch idea hai jaise Hoop stress in pressurised cylinders .
K / k 4 kyun scale karti hai
Yahan subtle part hai. Jab wall ek jagah dimple karti hai, material ko kahin jaana hota hai. Ek long gentle fold (small k ) ke liye, wall sach mein radially in aur out push karti hai, toh poori circumference stretch honi chahiye — badi hoop energy. Ek tight short fold (bada k ) ke liye, metal hoop stretch karne ki jagah sideways wall ke along neighbouring dimple mein slosh kar sakta hai, toh almost koi circumference change nahi — tiny hoop energy.
Buckle pattern ki geometry ke through ise work karne par (radial motion jo folds ko compatible rakhne ke liye chahiye woh tighter hone par fall off hoti hai) hoop energy milti hai
U hoop ∝ k 4 K w 0 2 .
Words mein: long loose folds hoop ke liye expensive hain (woh real circumference stretch force karte hain); tight folds material ko sideways shuffle karke cheat karte hain. Yahi hai 1/ k 4 — bending term ka exact mirror image.
Intuition Do springs kyun, optimum
k ∗ , aur unka geometric mean
Dono costs side by side rakho (per unit dimple depth):
U ( k ) = D k 4 + k 4 K .
Bending tight folds ke liye worst hai (bada k ); hoop loose folds ke liye worst hai (chhota k ). Wall woh k choose karti hai jo total ko minimum banaye. Derivative zero set karo:
d k d U = 4 D k 3 − k 5 4 K = 0 ⟹ k 4 ⋅ k 4 = D K ⟹ k ∗ = ( D K ) 1/8 .
K = E t / R 2 aur D = E t 3 / [ 12 ( 1 − ν 2 )] substitute karne par natural fold size milti hai k ∗ ∼ [ 12 ( 1 − ν 2 ) / ( R 2 t 2 ) ] 1/8 ∼ 1/ R t , yani wavelength λ ∗ = 2 π / k ∗ ∼ R t — dimples lagbhag R t wide hote hain, R ya L dono se chhote.
k ∗ par dono costs equal hote hain (D k 4 = K / k 4 ), aur a + b ≥ 2 ab se (equality jab a = b ) minimum value geometric mean hai
min k ( D k 4 + k 4 K ) = 2 D K .
Isliye parent ka clean result padhta hai σ cr t = 2 D ⋅ R 2 E t — do springs exactly balance, aur buckle length L vanish ho gayi (wall ne apna natural fold size ∼ R t find kar liya, tube height se independent — provided tube bahut short na ho, dekho §6).
Humne L define kiya tha lekin woh σ cr se drop out ho gaya. Yeh sirf middle regime mein true hai — ends aur aspect ratio L / R decide karte hain ki kis tarah ki buckling hoti hai.
L / R aur do buckling families
Local shell buckling (short/medium tube): wall chhote diamond dimples mein fold hoti hai size ∼ R t ke, aur σ cr = 0.605 E t / R L ko ignore karta hai , kyunki natural fold tube se bahut chhoti hai. Yahi woh regime hai jo parent page solve karta hai. s05 mein, left panel.
Global (Euler) column buckling (bahut long, slender tube): poora cylinder ek bent straw ki tarah sideways bow karta hai, exactly Euler column buckling mode, aur ab L strongly matter karta hai — critical load L 2 badhne par fall karta hai. s05 mein, right panel.
Quantitative switch-over: local shell stress 0.605 E t / R ko ek tube ke global Euler stress se compare karo (σ Euler = π 2 E ( R 2 + ... ) /... ; ek thin ring ke liye Euler stress ≈ 2 π 2 E ( R / L ) 2 hai pinned ends ke liye). Dono equal set karne par, crossover roughly yahan hota hai
R L ≈ 0.605 π 2 /2 ⋅ t R ≈ 2.85 R / t .
Concrete rule of thumb: L / R ≲ 2.85 R / t ke liye (yani zyaatar tanks) local dimpling govern karta hai aur tum 0.605 E t / R use karte ho; sirf is se aage bahut slender tubes ke liye Euler formula par switch karo. Example: R / t = 600 crossover L / R ≈ 2.85 600 ≈ 70 deta hai — ek tube ko global bowing jeeetne se pehle ~70 radii long hona padega.
Definition Boundary conditions: clamped vs free ends
Tube apne ends par kaise held hai yeh σ cr change karta hai:
Clamped (rigidly welded/bolted, jaise ek rocket interstage joint) — ends rotate ya slide nahi kar sakte; yeh stiffest case hai aur highest σ cr deta hai.
Simply supported (round held lekin rotate karne ke liye free) — thoda softer.
Free at one end (cantilever tube) — weakest ; global mode ke liye Euler load clamped-clamped versus ~4 factor se drop kar sakta hai.
Picture: clamped ends rims ke paas "fold ko flat hold" karte hain, toh buckle ko middle mein squeeze hona padta hai — trigger karna mushkil. Free ends wall ko flap karne dete hain, toh woh jaldi dab jaati hai.
Intuition Classical formula quietly "long enough, ends far away" kyun assume karta hai?
0.605 E t / R result woh local answer hai jahan natural dimple wavelength ∼ R t ends ke beech comfortably fit ho jaati hai, toh boundary conditions use barely touch karti hain. Bahut short tubes ke liye (few dimples fit) ends dominate karte hain aur σ cr 0.605 se upar rise karta hai; bahut long tubes ke liye (L / R ≈ 2.85 R / t crossover ke past) global Euler mode take over karta hai aur σ cr iske neeche fall karta hai. Clean number middle mein rehta hai — ise trust karne se pehle hamesha check karo ki tumhara L / R kahan hai.
Definition Critical stress
σ cr
Woh stress jis par ek perfect cylinder buckle karta hai (local mode):
σ cr = 3 ( 1 − ν 2 ) E R t ≈ 0.605 E R t ( ν = 0.3 ) .
Subscript "cr" = critical = tipping point.
Definition Knockdown factor
γ (gamma) aur ϕ (phi)
Real cans mein dents hote hain, toh woh sirf σ cr ke ek fraction par buckle karte hain. Hum ==γ == se multiply karte hain, ek number lagbhag 0.2 aur 0.7 ke beech:
σ a l l o w = γ σ cr , γ = 1 − 0.901 ( 1 − e − ϕ ) , ϕ = 16 1 R / t .
Yahan ϕ ("phi") sirf ek helper number hai jo geometry R / t se bana hai, aur e − ϕ exponential decay function hai. Poori story Imperfection sensitivity and knockdown factors aur NASA SP-8007 buckling of thin-walled cylinders mein hai.
e − ϕ kyun?
Jaise R / t badhta hai (thinner, floppier shell), ϕ badhta hai, e − ϕ 0 ki taraf shrink karta hai, aur γ apni floor 1 − 0.901 = 0.099 ki taraf drop karta hai. Exponential woh empirical curve hai jo NASA ne hazaaron real buckling tests ko fit karke banayi — thinner shells ko harder punish kiya jaata hai, smoothly. Yeh directly Rocket tank and interstage structural design mein feed karta hai.
Optimal wavenumber k star
End conditions clamped or free
Buckling vs yield lower wins
Right side cover karo aur dekho kya tum har ek memory se bata sakte ho.
R / t ratio kya measure karta hai, aur ratio kyun raw sizes nahi?Wall kitni "thin/flimsy" hai; proportion, metres nahi, buckling decide karta hai — same R / t utna hi tinny lagta hai.
Stress σ ki definition aur uske units? Force per unit area, σ = P / A ; pascals (N / m 2 ) mein measure hota hai, usually MPa.
Axial load kaunsa area carry karta hai, aur uska formula? Sirf wall ring: A = 2 π R t .
Young's modulus E kya describe karta hai? Material stiffness in stretch — unit fractional stretch per stress kitna chahiye.
Poisson's ratio ν kya describe karta hai, aur topic ko yeh kyun chahiye? Lengthwise stretch per sideways thinning; yeh D mein ( 1 − ν 2 ) plate-stiffening factor produce karta hai.
Yield aur buckling failure mein difference? Yield = material crush (strength); buckling = shape unstable ho jaati hai (stiffness) — buckling pehle strike karta hai yahan.
Radial displacement w kya hai? Wall perfectly round se in/out kitna poke karti hai — ek dimple ki depth.
Buckle wavenumber k kya count karta hai? Folds kitne tightly packed hain — bada k kai short folds, chhota k few long folds; wavelength λ = 2 π / k .
Flexural rigidity D ka formula aur meaning? D = E t 3 / [ 12 ( 1 − ν 2 )] — wall ko curve karne ka resistance; tight folds ke liye D k 4 cost.
Bending energy D k 4 kyun scale karti hai? Curvature fold ki 2nd derivative hai, ∝ k 2 ; bending energy ∝ curvature2 ∝ ( k 2 ) 2 = k 4 .
Membrane spring constant K (exact) kya hai, aur uski cost kya hai? K = E t / R 2 (koi ν nahi — yeh ek one-way hoop stretch hai, two-way plate bending nahi); K / k 4 cost, loose folds ke liye worst.
Hoop energy 1/ k 4 kyun scale karti hai? Loose folds real circumference stretch force karte hain; tight folds material ko sideways shuffle karne dete hain, toh hoop cost folds ke tighter hone par fall off hoti hai.
ν do springs mein kahan enter karta hai, aur sirf wahan kyun?Sirf D mein (two-directional plate bending, Poisson coupling); K mein nahi (single-direction hoop stretch).
Optimal wavenumber k ∗ aur fold size kya hai? k ∗ = ( K / D ) 1/8 ; wavelength
λ ∗ ∼ R t — natural dimple size,
L se independent.
σ cr geometric mean 2 D K kyun hai?D k 4 + K / k 4 minimize karna
k ∗ par do springs ko equally balance karta hai;
a + b ≥ 2 ab geometric mean deta hai.
2 D K se boxed formula tak algebra dikhao.D , K sub karo,
E 2 t 4 / [ 12 ( 1 − ν 2 ) R 2 ] milega,
E t 2 / R bahar nikalo,
t se divide karo,
12 = 2 3 2 cancel karta hai →
E t / [ R 3 ( 1 − ν 2 ) ] .
Buckling local se global kab switch hoti hai L / R par? Roughly
L / R ≈ 2.85 R / t ; iske neeche local dimpling govern karta hai (
0.605 E t / R use karo), iske upar Euler bowing jeetta hai.
Clamped vs free ends σ cr ko kaise affect karte hain? Clamped = stiffest, highest σ cr ; free end = weakest (global Euler load ~4× drop kar sakta hai). Hamesha boundary condition check karo.
γ aur ϕ kya karte hain?γ perfect stress ko real dents ke liye knock down karta hai;
ϕ = 16 1 R / t empirical
γ = 1 − 0.901 ( 1 − e − ϕ ) mein feed karta hai.