Barrowman equations — centre of pressure calculation for finned rockets: Visual walkthrough

3.4.8 · D2 · Physics › Rocket Flight Mechanics › Barrowman equations — centre of pressure calculation for fin

Step 1 — "Sideways air push" ka matlab kya hai

WHAT. Ek rocket normally apna pointed end seedha hawa ki taraf karke udta hai. Ab socho ise thoda dhakka lag jaata hai, toh uski axis airflow se thodi si angle par jhuk jaati hai. Us jhukav ka ek naam hai: angle of attack, jise $α$ (Greek letter "alpha") likhte hain. Dekho Angle of attack and restoring moment.

WHY. Jab rocket bilkul seedha hota hai ( $α = 0$ ), toh hawa symmetrically guzarti hai aur sabhi sides par barabar push karti hai — koi net sideways force nahi. Jaise hi yeh jhukta hai, har surface ka ek side hawa se zyada seedha takraata hai doosre se, aur ek sideways force paida hoti hai. Wahi sideways force rocket ko seedha karti hai (ya palat deti hai). Toh poora game $α$ se shuru hota hai.

PICTURE.

Sideways force, jo pink arrow ki tarah draw ki gayi hai, normal force kehlaati hai — yahan "normal" maths ka word hai perpendicular ke liye. Hum ise $N$ likhte hain: yeh rocket ki axis ke right angles par point karti hai.

Step 2 — Force ko ek pure number mein badalna

WHAT. Raw force $N$ kaafi messy cheezein depend karti hain: kitni tez ud rahe ho, hawa kitni thick hai, rocket kitna bada hai. Hum yeh sab strip out kar lete hain aur sirf shape ki apni personality rakhte hain divide karke: $C_{N} = \frac{N}{\frac{1}{2} ρ v ^{2} A _{r e f}}$

WHY. Do rockets jo shape mein identical hain lekin size mein alag hain, ya same rocket ek tez hawaadi din mein vs shaant din mein, inhe ek number se describe karna chahiye. $\frac{1}{2} ρ v^{2} A_{r e f}$ se divide karne par — jo natural "air pressure scale" hai — size, speed aur air-density ki dependency cancel ho jaati hai. Dekho Slender-body aerodynamic theory jisme bataya gaya hai ki yeh exact denominator kyun sahi hai.

PICTURE.

Step 3 — Hum force ko nahi, slope ko kyun track karte hain

WHAT. Hume ek fixed tilt par $C_{N}$ ki parwah nahi. Hume parwah hai ki jaise tilt badhti hai, yeh kitni tezi se badhta hai. Woh rate ek slope hai: $C_{N α} = \frac{\partial C _{N}}{\partial α}$

WHY this tool — the derivative. Symbol $\frac{\partial C _{N}}{\partial α}$ (padho "the derivative of $C_{N}$ with respect to $α$ ") bilkul usi sawaal ka jawaab deta hai "agar main thoda aur jhukoon, toh kitni extra sideways force milegi?" Derivative perfect tool hai kyunki stability ek chhoti disturbance ka response hai — aur "response per unit disturbance" bilkul wahi hai jo ek slope measure karta hai. $α = 0$ par force zero hai, toh number bekaar hai; uski rate of change sab kuch hai.

PICTURE.

Step 4 — Nose cone exactly "2" deta hai

WHAT. Kisi bhi slender nose ke liye, slender-body theory ek clean result deta hai: $C_{N α, nose} = 2$

WHY "2" aur kuch nahi. Slender-body theory kehti hai ek pointed body ka lift slope $2 \times \frac{base area}{A _{r e f}}$ hota hai. Agar hum nose ki apni base ko reference area choose karein, toh woh fraction exactly $1$ hai, aur $2$ bachta hai. Remarkably yeh is par depend nahi karta ki nose cone hai, ogive hai ya parabola — sirf nose ka peechla hissa kitna bada hai uss par depend karta hai.

PICTURE.

Recall Base area kyun, tip nahi?

Nose ki cross-section badhna kahan band hoti hai? ::: Base par — us ke aage body seedhi hai aur (is theory mein) koi naya lift nahi banati. Nose ki saari "growing" tip aur base ke beech hoti hai, isliye base area total set karta hai.

Step 5 — Nose ki force kahan act karti hai: "growing area" ka centroid

WHAT. Nose ka single balance point hai: $\overset{ˉ}{X}_{nose} = k L_{nose}, k = \frac{2}{3} for a cone.$

WHY the centroid, aur why the integral. Sideways push ek jagah apply nahi hoti — yeh nose ke along phaili hui hai. Slender-body theory kehti hai ki unit length per local push $\frac{d S}{d x}$ ke proportional hai, jo ki rate hai jis se cross-section area $S$ badhti hai jab tum tip se $x$ distance peeche jaate ho. Ek phaile hue load ko ek equivalent point mein compress karne ke liye hum uska area-weighted average position lete hain — ek centroid — jo exactly wahi hai jo ek integral calculate karta hai: $\overset{ˉ}{X} = \frac{\int _{0}^{L} x \frac{d S}{d x} d x}{\int _{0}^{L} \frac{d S}{d x} d x}$

Cone ke liye radius seedha-line badhta hai, $r \propto x$ , toh area $S \propto x^{2}$ aur $\frac{d S}{d x} \propto x$ . Substitute karne par: $\overset{ˉ}{X} = \frac{\int _{0}^{L} x \cdot x d x}{\int _{0}^{L} x d x} = \frac{L ^{3} /3}{L ^{2} /2} = \frac{2}{3} L .$

PICTURE.

Step 6 — Fins: thin-airfoil flow se ek bada slope

WHAT. Har flat fin ek chhote wing ki tarah behave karta hai. Saare $n$ fins ko jod ke aur jo moti body par lage hain uska correction karke: $C_{N α, fins} = K_{f b} \frac{4 n ( s / d ) ^{2}}{1 + 1 + ( \frac{2 ℓ _{f}}{C _{r} + C _{t}} ) ^{2}}$

WHY is formula ki yeh shape. Chhote $α$ par ek thin flat plate apne span squared ke proportional lift banata hai (bada span = zyada airflow pakadta hai) — yahi upar $(s / d)^{2}$ hai, times $4 n$ for $n$ fins. Upar square-root wala ugly denominator ek aspect-ratio correction hai: chhote, mote fins lambe thin fins se kam efficient hote hain, aur $ℓ_{f}$ (woh slanted mid-chord length) woh sweep capture karta hai jo efficiency kam karta hai. Dekho Slender-body aerodynamic theory.

PICTURE.

Step 7 — Fins ki force kahan act karti hai (sweep use peeche kheenchta hai)

WHAT. Fin centre of pressure: $\overset{ˉ}{X}_{fins} = X_{f} + \frac{X _{r}}{3} \cdot \frac{C _{r} + 2 C _{t}}{C _{r} + C _{t}} + \frac{1}{6} [(C_{r} + C_{t}) - \frac{C _{r} C _{t}}{C _{r} + C _{t}}]$

WHY yeh teen pieces.

$X_{f}$ — fin ki leading edge body par kahan se shuru hoti hai. Har cheez yahan se measure hoti hai.
Beech wala piece $X_{r}$ ke saath (sweep: tip root se kitna peeche set hai) — ek taper-weighted centroid. Agar tip peeche swept hai, toh uska load aur peeche act karta hai, toh CP peeche jaata hai.
Bracket wala piece — ek chord-averaging correction jo root aur tip chords ko blend karta hai.

PICTURE.

Recall Degenerate check: ek rectangular, unswept fin

$C_{t} = C_{r} = C$ aur $X_{r} = 0$ rakho. Beech wala term vanish ho jaata hai; bracket $\frac{1}{6} [2 C - \frac{C}{2}] = \frac{C}{4}$ ban jaata hai. ::: Toh $\overset{ˉ}{X}_{fins} = X_{f} + \frac{C}{4}$ — CP quarter-chord par baithta hai, exactly classic thin-airfoil result. Accha, formula sahi tarike se reduce ho raha hai.

Step 8 — Sab kuch combine karna: weighted balance point

WHAT. Hamare paas, har part $i$ ke liye, ek slope $C_{N α, i}$ hai (yeh kitni kadi push karta hai) aur ek location $\overset{ˉ}{X}_{i}$ hai (yeh kahan push karta hai). Poore rocket ka CP hai: $\overset{ˉ}{X}_{c p} = \frac{\sum _{i} C _{N α, i} X ˉ _{i}}{\sum _{i} C _{N α, i}}$

WHY ek weighted average — torques balance karne se. Torque = force × lever arm. Saare parts ka total torque ek kalpanaa ki gayi force ke barabar hona chahiye jo ek point $\overset{ˉ}{X}_{c p}$ par hai: $\sum_{i} N_{i} \overset{ˉ}{X}_{i} = \overset{ˉ}{X}_{c p} \sum_{i} N_{i} .$ Dono sides ko $\frac{1}{2} ρ v^{2} A_{r e f}$ se divide karo aur har $N_{i}$ $C_{N α, i}$ ban jaata hai. $\overset{ˉ}{X}_{c p}$ ke liye solve karo — aur boxed formula nikal aata hai. Jis part ka bada slope hai (fins!) woh balance point ko apni taraf kheenchta hai.

PICTURE.

Phir CG se compare karo static margin paane ke liye: CG $38$ cm par aur $d = 4$ cm ke saath, margin $= \frac{45.3 - 38}{4} = 1.8$ calibres — comfortably stable.

Step 9 — Edge cases jinpar kabhi mat phislo

WHAT / WHY / PICTURE saath mein. Formula mein kuch quiet corners hain:

Boat-tail (body narrow hoti hai), $d_{2} < d_{1}$ . Transition slope $2 [(d_{2} / d)^{2} - (d_{1} / d)^{2}]$ negative ho jaata hai — ek destabilising forward pull. Yeh ek real effect hai; sign rakho.
Zero angle of attack, $α = 0$ . Har $N_{i} = 0$ hai, lekin slopes $C_{N α, i}$ aur locations $\overset{ˉ}{X}_{i}$ unchanged hain. CP location $α$ par depend nahi karta — isliye hum ise ek fixed number ki tarah quote kar sakte hain.
Saare slopes zero (ek finless, noseless tube). Denominator $\sum C_{N α, i} \to 0$ aur $\overset{ˉ}{X}_{c p}$ undefined hai — physically koi restoring force nahi aur koi CP nahi. Maths honestly jawaab dene se mana kar deta hai.
Chhote $α$ se zyada ya Mach 1 ke paas. Step 3 ki straight-line assumption toot jaati hai; asli $C_{N}$ curve mein bend aata hai. Barrowman sirf $α ≲ 1 0^{\circ}$ aur subsonic ke liye trustworthy hai.

Ek-picture summary

Ek board par sab kuch: har part ek push strength contribute karta hai (uska slope $C_{N α, i}$ , arrow thickness ki tarah drawn) ek position ( $\overset{ˉ}{X}_{i}$ ) par. Rocket ka CP unka strength-weighted balance point hai. Ise CG ke peeche rakho → restoring moment flight ko seedha karta hai, ek stable trajectory kisi bhi 6-DOF simulation mein feed karta hai.

Recall Poore walkthrough ki Feynman retelling

Rocket ko thoda tilt karo — woh tilt $α$ hai. Kyunki ab yeh hawa ke thoda sideways hai, har surface ko thodi sideways shove milti hai, $N$ . Hum speed aur air-density scrub kar dete hain ek clean number rakhne ke liye, $C_{N}$ , aur phir woh ek hi sawaal poochte hain jis ki stability ko parwah hai: jaise main aur jhukoon woh shove kitni tezi se badhti hai? — wahi slope $C_{N α}$ hai. Nose, ek neat theorem se, hamesha exactly 2 score karta hai aur apni length ke do-teesre par push karta hai (apni growing width ka centre). Har fin ek tiny wing ki tarah act karta hai; bade spans kadi push karte hain, sweep unke push spot ko peeche kheenchti hai, aur ek moti body par baithna unhe aur bhi harder push karta hai ( $K_{f b}$ ). Aakhir mein, woh single spot dhundhne ke liye jahan yeh saari shoves balance karti hain, hum wahi karte hain jo tum karte ho ek seesaw par kai bachchon ka balance point dhundhne ke liye: unki positions ka strength-weighted average. Fins sabse kadi push karte hain, isliye woh tug-of-war jeette hain aur balance point ko tail ke paas kheench lete hain — bilkul wahan jo tum chahte ho, weight balance ke safely peeche, taaki rocket ek acchi tarah se phenke gaye dart ki tarah ude.