Dynamic stability — pitch - yaw damping derivatives

3.4.11 · Physics › Rocket Flight Mechanics

WHAT: hum kya describe kar rahe hain?

Jab ek rocket pitch karta hai (naak upar/neeche) ya yaw karta hai (naak left/right), center of gravity (CG) ke peeche ki har surface sideways hawa mein sweep karti hai. Yeh sideways sweep sirf rotation ki wajah se ek local angle of attack create karti hai, jo ek restoring aerodynamic force produce karti hai — aur isliye ek torque jo rotation ko oppose karta hai.

HOW: first principles se damping moment derive karo

Step 1 — Ek rotating rocket har point ko ek local vertical velocity deta hai. CG ko origin maano, $x$ ko aft measure karo. CG ke peeche $x$ distance par ek point, jab body $q$ rate se pitch karti hai, airstream ke relative downward speed se move karta hai $w_{local} = q x .$ Yeh step kyun? Rigid-body rotation: kisi point ki velocity = $ω \times r$ . CG ke baare mein ek pitch rate tail ko utha deta hai jab naak neeche jaati hai, toh $+ x$ par ek point (CG ke peeche) transverse velocity $q x$ gain karta hai.

Step 2 — Woh transverse velocity ek local angle of attack hai. Airflow $V$ speed se axially aati hai, plus yeh naya transverse $q x$ , toh us element ko ek extra angle of attack milta hai $Δ α (x) = \frac{w _{local}}{V} = \frac{q x}{V} .$ Yeh step kyun? Angle of attack bas $arctan (transverse / axial) \approx$ transverse/axial hai small angles ke liye.

Step 3 — Extra angle of attack → extra local normal force. Local lift-curve slope $c_{n_{α}} (x)$ aur local area contribution wala ek element ek incremental normal force generate karta hai

=\tfrac12\rho V^2\,c_{n_\alpha}(x)\,\frac{qx}{V}\,dA.$$ *Yeh step kyun?* Force = dynamic pressure $\times$ effective area $\times$ (slope $\times$ angle). Standard linear aero. **Step 4 — Har force moment arm $x$ par ek torque banata hai; integrate karo.** Restoring torque (CG ke baare mein moment) hai $$M=-\int x\,dN=-\tfrac12\rho V^2\frac{q}{V}\int c_{n_\alpha}(x)\,x^2\,dA .$$ *Minus sign kyun?* $\Delta\alpha$ se force motion ko *reduce* karne ki taraf point karta hai (neeche jaata tail upar push hota hai), toh resulting moment $q$ ko oppose karta hai. **Integral mein $x^2$ aata hai**: yeh damping ka geometric dil hai — CG ke door peeche ki surfaces dominate karti hain, kyunki woh zyada tezi se sweep karti hain *aur* longer lever arm bhi hai. **Step 5 — Non-dimensionalize karo.** $C_m=M/(\tfrac12\rho V^2 S d)$ aur $\hat q = qd/(2V)$ define karo. Substitute karne par, $$\boxed{\,C_{m_q}=\frac{\partial C_m}{\partial \hat q} = -\frac{2}{S\,d^2}\int c_{n_\alpha}(x)\,x^2\,dA\, }$$ > [!formula] Yeh result yaad rakhna zaroori hai > $$C_{m_q}=-\frac{2}{S d^2}\int c_{n_\alpha}\,x^2\,dA \;<\;0 \quad\text{(damping)}$$ > **Sign negative hai → yeh hamesha rotation ko oppose karta hai.** Magnitude CG ke baare mein aerodynamic area ke **second moment** ke saath scale karta hai. Fins ko zyada peeche le jaao → $x$ badhta hai → damping badhti hai (quickly). --- ## HOW: damping equation of motion mein kaise enter karta hai *Static* restoring moment ($C_{m_\alpha}\alpha$, spring) ko *damping* moment ($C_{m_q}\hat q$, dashpot) ke saath combine karo. Linearised pitch dynamics: $$I_y\ddot\theta - \tfrac12\rho V^2 S d\left(C_{m_\alpha}\,\theta + C_{m_q}\frac{d}{2V}\dot\theta\right)=0.$$ Yeh ek **damped harmonic oscillator** $\ddot\theta + 2\zeta\omega_n\dot\theta + \omega_n^2\theta=0$ hai jisme $$\omega_n=\sqrt{\frac{-\,\tfrac12\rho V^2 S d\,C_{m_\alpha}}{I_y}},\qquad \zeta=\frac{-\tfrac14\rho V S d^2\,C_{m_q}}{2 I_y\,\omega_n}.$$ *Yeh kyun matter karta hai:* $\omega_n$ ko $C_{m_\alpha}<0$ chahiye (static stability); $\zeta>0$ ko $C_{m_q}<0$ chahiye (damping). Sach mein calm down karne ke liye **dono chahiye**. ![[3.4.11-Dynamic-stability-—-pitch-yaw-damping-derivatives.png]] --- ## Worked examples > [!example] 1. Single tail-heavy fin set > Ek rocket ki saari aerodynamic area fins par concentrated hai, $x_f = 0.6\,$m CG ke peeche, jisme $c_{n_\alpha}A_f = 8\text{ m}^2$ (per rad) hai. Reference $S=0.008\text{ m}^2$, $d=0.1$ m. $C_{m_q}$ nikalo. > > $C_{m_q}=-\dfrac{2}{Sd^2}\,(c_{n_\alpha}A_f)\,x_f^2 = -\dfrac{2}{0.008\times 0.01}\times 8 \times 0.36$ > **Yeh step kyun?** Jab area ek station par lumped ho toh integral bas value $\times x_f^2$ hai. > $=-\dfrac{2\times8\times0.36}{8\times10^{-5}} = -7.2\times10^{4}.$ Strongly damped — achha hai. > [!example] 2. Fin distance double karna > Unhi fins ko $x_f=1.2$ m par le jaao. Kyunki $C_{m_q}\propto x_f^2$ hai, damping **chaar guna** ho jaati hai: $\approx -2.9\times10^{5}$. > **Yeh step kyun?** Integral mein $x^2$ ka matlab hai ki fins ko double distance peeche le jaane se damping do nahi, chaar guna ho jaati hai. Yeh sabse important design lever hai. > [!example] 3. Kya oscillation sach mein decay kar rahi hai? > Diya hai $\zeta=0.06>0$, response $\propto e^{-\zeta\omega_n t}\cos(\omega_d t)$ decay karta hai. Agar ek design error ne fins ko CG ke *aage* flip kar diya, toh $x^2$-weighting abhi bhi ek moment deta hai lekin ab $C_{m_\alpha}>0$ hai: $\omega_n^2<0$, roots real & positive → pure **divergence**, koi oscillation hi nahi. > **Yeh step kyun?** Damping sirf ek *statically stable* body ki help karta hai; ek unstable body par yeh tumhe bacha nahi sakta. --- ## Common mistakes > [!mistake] "Static stability guarantee karta hai ki rocket safe hai." > **Kyun sahi lagta hai:** Ek statically stable rocket *sach mein* wind ki taraf wapas turn karna shuru karta hai, jo stabilizing lagta hai. **Flaw yeh hai:** woh overshoot karta hai aur oscillate karta hai; agar $C_{m_q}$ bahut weak ho (ya, high altitude par, air density $\rho$ bahut kam ho), toh oscillation barely damp hoti hai aur tumble mein excite ho sakti hai. **Fix:** hamesha check karo $\zeta>0$ *aur* $\zeta$ itna bada ho (typically $\zeta \gtrsim 0.05$–$0.1$ pure flight envelope mein). > [!mistake] "Damping fin distance $x$ ke saath scale karta hai." > **Kyun sahi lagta hai:** Longer lever arm = zyada torque, toh surely linear hai. **Flaw yeh hai:** har surface jo sweep karta hai uski *velocity* bhi $x$ ke saath badhti hai (Step 1), toh torque $\propto x \times x = x^2$ hai. **Fix:** **second moment** yaad rakho — damping distance ke *square* ke saath jaati hai. > [!mistake] "$C_{m_q}$ rocket ka ek constant hai." > **Kyun sahi lagta hai:** Yeh wind-tunnel tables mein ek fixed geometry number hai. **Flaw yeh hai:** *physical damping moment* $=\tfrac12\rho V S d^2 C_{m_q}\cdot(\text{rate})$ $\rho$ aur $V$ par depend karta hai. Jaise rocket climb karta hai, $\rho\downarrow$ aur damping dramatically weak ho jaati hai. **Fix:** $\zeta$ ko har altitude par evaluate karo, sirf launch par nahi. --- ## #flashcards/physics $C_{m_q}$ kaunsa physical effect quantify karta hai? ::: Woh aerodynamic moment jo pitch *rate* ko oppose karta hai — rotational shock absorber jo oscillations ko damp karta hai. Damping $x^2$ ke saath (area ka second moment) kyun scale karta hai? ::: $x$ distance par ek surface transversely $qx$ speed se sweep karti hai (∝x) AUR $x$ (∝x) lever arm par kaam karti hai; torque ∝ $x·x = x^2$. Stability ke liye $C_{m_q}$ ka kya sign hona chahiye, aur iska kya matlab hai? ::: Negative; moment hamesha rotation rate ko oppose karta hai, oscillation se energy remove karta hai. Static stability ke liye kaunsa derivative negative hona chahiye? ::: $C_{m_\alpha}<0$ (angle of attack ke against restoring moment). Ek decaying oscillation ke liye kaunsi do conditions milke zaroori hain? ::: $C_{m_\alpha}<0$ (deta hai $\omega_n^2>0$) AUR $C_{m_q}<0$ (deta hai $\zeta>0$). Jab rocket altitude gain karta hai toh damping ratio kaise badalta hai, aur kyun? ::: Yeh decrease hota hai, kyunki physical damping moment $\propto \rho V$ hai aur air density $\rho$ altitude ke saath girti hai. Ek rocket ke liye $C_{n_r}=C_{m_q}$ kyun hota hai? ::: Axial symmetry — pitch aur yaw geometrically identical planes hain. Non-dimensional pitch rate $\hat q$ ki definition kya hai? ::: $\hat q = qd/(2V)$, pitch rate ko chord-time se scale kiya gaya. Agar fins ko CG ke 0.6 m se 1.2 m peeche le jaao, toh damping kitne factor se badlega? ::: ×4 (kyunki ∝ $x^2$). --- > [!recall]- Feynman: ek 12-saal ke bachche ko samjhao > Socho ki tum ek office chair par arms failaake ghoom rahe ho. Hawa tumhare haathon par push back karti hai aur tumhe slow karti hai — jitna tezi ghoomoge, utna zyada push back hoga. Ek rocket ke fins unhi arms ki tarah hain: jab rocket wobble karta hai, fins hawa mein sweep karte hain aur us direction mein push hote hain jo *wobble ko rokta hai*. Jo fins tail mein bahut peeche hain woh hawa mein bahut tezi se sweep karte hain AUR unka lamba "handle" bhi hota hai, toh woh wobbles rokne mein double-good hain. Agar fins bahut chhote hon ya hawa bahut patli ho (bahut upar aasman mein), toh wobble rukne mein bahut time lagta hai — aur yeh dangerous hai. > [!mnemonic] Yaad rakho > **"Far Fins Fight Faster, Force ∝ x²."** Aur do rules ke liye: **"Spring kehta hai kaunsi taraf (α), Dashpot kehta hai kitni tezi se (q)."** Tumhe spring ($C_{m_\alpha}$) aur dashpot ($C_{m_q}$) dono chahiye. ## Connections - [[Static Stability — Center of Pressure & Margin]] — $C_{m_\alpha}$ provide karta hai; $C_{m_q}$ ke saath milke zaroori hai. - [[Damped Harmonic Oscillator]] — ODE template ($\zeta$, $\omega_n$) jisme yeh sab reduce hota hai. - [[Fin Design & Sizing]] — geometry kaise $x^2$ integral set karta hai. - [[Atmospheric Density vs Altitude]] — kyun $\rho\downarrow$ flight mein damping erod karta hai. - [[Barrowman Equations]] — standard tarika jisme area distribution $c_{n_\alpha}(x)$ compute hoti hai. - [[Moments of Inertia of a Rocket]] — $I_y$ jo $\omega_n$ set karta hai. ## 🖼️ Concept Map ```mermaid flowchart TD SS[Static stability] -->|causes overshoot| OSC[Oscillation after snap-back] OSC -->|question| DS[Dynamic stability] DS -->|needs| DM[Damping moment] DM -->|quantified by| DD[Damping derivatives Cmq Cnr] DD -->|defined as| DEF[dCm / d q-hat] Q[Pitch rate q] -->|rigid-body rotation| WL[Local velocity w = q x] WL -->|divide by V| DA[Local angle of attack q x / V] DA -->|linear aero| DN[Incremental normal force dN] DN -->|times arm x, integrate| M[Restoring torque M] M -->|opposes rotation| DM DD -->|axisymmetry| SYM[Cnr equals Cmq] DD -->|stable if| NEG[Cmq less than 0] NEG -->|oscillations| DECAY[Oscillations die out] ```