3.4.5 · D3 · Physics › Rocket Flight Mechanics › 6DOF equations — translational (Newton), rotational (Euler's
Intuition Yeh page kis liye hai
Parent note ne do laws build ki thi: Newton sliding (translation) ke liye aur Euler twisting
(rotation) ke liye. Yahan hum unhe har us situation ke against stress-test karte hain jo ek rocket
tumhe de sakta hai —
har sign, har zero, har degenerate shape, plus ek word problem aur ek exam twist.
Pehle scenario matrix padho, phir har example solve karo. End tak tumhe koi aisa case nahi milna chahiye
jo tumne pehle solve hote nahi dekha ho.
Kuch bhi shuru karne se pehle: woh symbols jo hum baar baar use karenge, plain words mein.
Definition Symbols ki cast (sab ek baar, yahan define)
m — rocket ki mass (kitna saamaan hai), kilograms mein.
v = ( u , v , w ) — center of mass (balance point) ki velocity body
axes mein likhi: u = forward speed, v = sideways speed, w = up/down speed. Units m/s.
ω = ( p , q , r ) — rocket kitni tezi se twist karta hai: p = roll (ek log ki tarah barrel-roll),
q = pitch (nose up/down), r = yaw (nose left/right). Units rad/s.
F = ( F x , F y , F z ) — total external force (push), newtons mein.
M = ( M x , M y , M z ) — total external moment/torque (twisting effort), N·m mein.
I x , I y , I z — teen body axes ke baare mein moments of inertia : "is axis ko spin karwana kitna mushkil hai". Units kg·m². Inertia tensor and principal axes mein build kiya gaya hai.
I t — ek axisymmetric body ka transverse moment of inertia: jab do "side" axes identical hote hain to hum I y = I z = I t likhte hain (do ki jagah ek shared number). Units kg·m².
ℓ — ek lever arm : center of mass se kisi force ki line tak ki perpendicular distance.
Force ko uske lever arm se multiply karne par moment milta hai, M = F ℓ . Units metres.
Symbol × cross product hai — yeh do arrows leta hai aur ek teesra arrow return karta hai
jo dono ke perpendicular ho, jiska length "yeh dono kitne parallel nahi hain" yeh batata hai.
Do master equations jo hum baar baar apply karenge (parent se):
Har row ek class of situation hai. Last column us worked example ka naam batata hai jo use cover karta hai.
#
Case class
Kya special hai
Covered by
A
Zero rotation (ω = 0 )
sab cross terms vanish → plain 1D Newton
Ex 1
B
Pure turn, constant body speed
v ˙ = 0 lekin real acceleration exist karta hai
Ex 2
C
Sign of a coupling term (q r > 0 vs < 0 )
gyroscopic torque kis direction mein point karta hai
Ex 3
D
Symmetric spinner (I y = I z )
steady coning, oscillation frequency
Ex 4
E
Fully asymmetric, torque-free (I x = I y = I z )
intermediate-axis instability
Ex 5
F
Degenerate inertia (sphere I x = I y = I z )
all coupling dies, ω ˙ from torque only
Ex 6
G
Decoupled single-axis maneuver
one moment, clean angular accel
Ex 7
H
Word problem (real thruster + gravity)
build F , M from a story
Ex 8
I
Limiting / exam twist (mass → variable)
why naive d ( m v ) / d t is wrong
Ex 9
Worked example Vertical climb, no spin
Ek rocket seedha upar jaata hai. Body-x upar point karta hai. Thrust T = 30 , 000 N upar, weight
m g with m = 2000 kg, g = 9.8 m/s². Koi rotation nahi: p = q = r = 0 , aur v = w = 0 .
Nikalo forward acceleration u ˙ .
Forecast: padhne se pehle guess karo — kya answer sirf ( T − m g ) / m hona chahiye? Kyun ya kyun nahi?
Newton-x likho: F x = m ( u ˙ + q w − r v ) .
Yeh step kyun? Isme hi u ˙ hai, jo hamara unknown hai.
q = r = 0 set karo: terms q w aur r v khatam ho jaate hain, bacha F x = m u ˙ .
Yeh step kyun? Rotation na hone par koi transport correction nahi hoti — body frame aur
inertial frame acceleration par agree karte hain.
Net force hai F x = T − m g = 30000 − 2000 ( 9.8 ) = 10400 N.
Yeh step kyun? Dono forces body-x ke along act karte hain; upar positive hai, weight negative hai, isliye
hum dono vertical forces ko unke signs ke saath add kar dete hain.
Solve karo: u ˙ = F x / m = 10400/2000 = 5.2 m/s².
Yeh step kyun? Net force ko mass se divide karo — F x = m u ˙ ka final rearrangement.
Verify karo: units N / kg = m/s 2 ✓. Agar T = m g = 19600 N to u ˙ = 0
(hovering) — intuition se match karta hai. Yeh woh "floor" hai jis par har richer case reduce hona chahiye.
Worked example Coordinated turn — acceleration saamne chupi hui hai
Ek rocket constant body velocity u = 200 m/s, v = w = 0 se cruise karta hai, jabki steady
r = 0.1 rad/s se yaw kar raha hai (p = q = 0 ). Body speed feel hoti hai constant. Mass m = 1500 kg.
Nikalo zaruri sideways force F y , aur true inertial acceleration.
Forecast: v ˙ = 0 ke saath, kya F y zero hai? (Nahi hai — yahi toh pura lesson hai.)
Newton-y : F y = m ( v ˙ + r u − pw ) .
Yeh step kyun? Sideways force yahan hai; term r u transport correction hai.
Constant body velocity → v ˙ = 0 . p = 0 ke saath: F y = m r u .
Yeh step kyun? Bhale hi numbers u , v , w na badlein, velocity arrow yaw ki wajah se
swing ho raha hai — uski inertial direction badal rahi hai, isliye real acceleration exist karta hai.
Plug in karo: F y = 1500 × 0.1 × 200 = 30000 N.
Yeh step kyun? Required force nikaalane ke liye F y = m r u mein given numbers substitute karo.
Inertial acceleration ki magnitude hai a = F y / m = r u = 20 m/s².
Yeh step kyun? Yeh exactly centripetal acceleration v 2 / R hai jahan R = u / r = 2000 m:
u 2 / R = 20 0 2 /2000 = 20 m/s² ✓.
Neeche figure dekho. Gray dashed curve R = 2000 m radius ki flight arc hai. Blue arrow
velocity v hai, arc ke tangent — uski length (200 m/s) kabhi nahi badlti, lekin
green curved arrow (yaw r ) uski direction ko swing karte rehta hai. Orange arrow force
F y = 30000 N hai jo turn center ki taraf point kar raha hai: wahi force seedhi
flight ko circle mein bend karti hai. Toh "constant speed" turn ko bhi real sideways push chahiye.
Verify karo: a tak pahunchne ke dono raaste (transport term r u aur centripetal u 2 / R ) 20 m/s² par agree karte hain.
F y ke units: kg·(rad/s)·(m/s) = kg·m/s² = N ✓ (radians dimensionless hote hain).
Worked example Gyroscopic torque kis direction mein push karta hai?
Ek rocket roll p = 5 rad/s aur yaw r = 2 rad/s ke saath spin karta hai, is instant pitch q = 0 hai.
Inertias: I x = 10 , I y = 20 , I z = 40 kg·m². Koi applied moments nahi (M = 0 ).
Nikalo pitch acceleration q ˙ aur uska sign .
Forecast: zero pitch aur zero torque ke saath, kya nose usi jagah rahega? q ˙ ke sign ka guess karo.
Euler-y : M y = I y q ˙ + ( I x − I z ) r p . M y = 0 set karo.
Yeh step kyun? Hume q ˙ chahiye; coupling term ( I x − I z ) r p uska sign decide karta hai.
Solve karo: q ˙ = − I y ( I x − I z ) r p = − 20 ( 10 − 40 ) ( 2 ) ( 5 ) .
Yeh step kyun? Euler-y ko rearrange karna; ( I x − I z ) ka sign times r p jo hum padh rahe hain.
Sign pieces evaluate karo: ( I x − I z ) = − 30 (negative), r p = 10 (positive). Saamne wale
minus aur negative bracket ka product ek positive q ˙ deta hai.
Yeh step kyun? Number crunch karne se pehle signs track karo, taaki physics (nose kis taraf pitch hoga)
arithmetic mein dabbi na rahe.
Number: q ˙ = − 20 − 30 ( 10 ) = + 15 rad/s².
Yeh step kyun? Do axes par ek saath spin karne wala body ek teesra twist develop karta hai khud hi —
yahi gyroscopic coupling hai.
Verify karo: yaw sign ko r = − 2 kar do: to r p = − 10 aur q ˙ = − 20 − 30 ( − 10 ) = − 15
rad/s² — same size, opposite sign, jaisa formula ki symmetry demand karti hai ✓.
Worked example Sounding rocket coning frequency
Ek spin-stabilized rocket axisymmetric hai: I y = I z = I t = 30 kg·m², spin inertia
I x = 6 kg·m². Yeh p 0 = 20 rad/s par spin karta hai, torque-free. Ek chhoti si disturbance use
transverse rates de deti hai. Nikalo coning (nutation) frequency.
Forecast: kya transverse wobble grow karega, decay karega, ya hamesha oscillate karta rahega?
Euler-x : I x p ˙ = 0 ⇒ p = p 0 constant.
Yeh step kyun? I y = I z ke saath x -coupling term ( I z − I y ) q r = 0 ; spin rate lock ho jaati hai.
Euler-y , z : I t q ˙ = ( I x − I t ) r p 0 aur I t r ˙ = − ( I x − I t ) q p 0 .
Yeh step kyun? Symmetric inertias substitute karna; note karo ki q aur r ek doosre ko feed karte hain.
Define karo λ = I t I x − I t p 0 = 30 6 − 30 ( 20 ) = − 16 rad/s. Phir
q ˙ = λ r , r ˙ = − λ q ⇒ q ¨ = − λ 2 q .
Yeh step kyun? Do first-order coupled equations ek oscillator equation mein combine ho jaate hain —
yeh vector ( q , r ) ki circular motion ki pehchaan hai.
Coning frequency hai ∣ λ ∣ = 16 rad/s (lagbhag 2.5 Hz).
Yeh step kyun? q ¨ = − λ 2 q ek spring ka equation hai; uski frequency ∣ λ ∣ hai.
Neeche figure dekho. Blue aur orange curves do transverse rates q ( t ) aur
r ( t ) hain. Notice karo ki yeh ek cosine aur ek (negative) sine hain — quarter-cycle apart — isliye jab ek badhta hai
doosra girta hai: pair apne aap ko circle mein chase karta rahta hai. Khaas baat yeh hai ki curves na badhte hain na ghatte;
yeh ± 1 ke beech bounded rehte hain. Green dashed line ek full period 2 π /∣ λ ∣ ≈ 0.39 s mark karti hai. Woh bounded, repeating wobble exactly wahi hai jo spin stabilization ko safe banata hai.
Verify karo: wobble na grow karta hai na decay — yeh bounded oscillation hai, exactly isliye
spin stabilization kaam karta hai. Units: (kg·m²/kg·m²)(rad/s)
= rad/s ✓.
Worked example Middle axis ke baare mein tumbling
Ek body ke teen alag inertias hain I x = 2 < I y = 5 < I z = 9 kg·m² (sab distinct).
Yeh almost purely intermediate axis y ke baare mein spin karta hai: q = 10 rad/s, tiny
p = r = 0.01 rad/s. Torque-free. Nikalo kya chhote p , r grow karte hain ya chhote rehte hain.
Forecast: kaun sa axis unstable hai — sabse chhota, middle, ya sabse bada?
Euler-x : I x p ˙ = ( I y − I z ) q r , Euler-z : I z r ˙ = ( I x − I y ) pq (q ≈ const ke saath).
Yeh step kyun? Hum bade spin q ke baare mein linearize karte hain; p , r chhote deviations hain.
Euler-x differentiate karo aur Euler-z substitute karo:
p ¨ = I x ( I y − I z ) q r ˙ = I x I z ( I y − I z ) ( I x − I y ) q 2 p .
Yeh step kyun? Hum p mein ek single equation chahte hain taaki growth vs oscillation padh sakein.
Coefficient compute karo k = I x I z ( I y − I z ) ( I x − I y ) q 2 = 2 ⋅ 9 ( 5 − 9 ) ( 2 − 5 ) ( 100 ) = 18 ( − 4 ) ( − 3 ) ( 100 ) = 18 12 ( 100 ) = 3 200 ≈ 66.7 .
Yeh step kyun? k ka sign sab kuch decide karta hai: k > 0 matlab p ¨ = + k p →
exponential blow-up.
Kyunki k > 0 , p e k t ki tarah grow karta hai — spin unstable hai. Body flip ho jaata hai.
Yeh step kyun? Middle axis ke baare mein spin karne par hamesha ek negative aur ek positive
bracket milta hai, jinका product positive hota hai → instability.
Verify karo: largest axis z ko spin axis lekar repeat karo (r large): do brackets
( I z − I x ) ( I y − I z ) -type ban jaate hain opposite signs ke saath, k < 0 dete hain → oscillation → stable. Yahi
Dzhanibekov effect / intermediate axis theorem hai. Growth rate k = 200/3 ≈ 8.16 /s ✓.
Worked example Ball-shaped satellite
Ek satellite uniform sphere hai: I x = I y = I z = I = 12 kg·m². Ek torque M z = 3 N·m
z ke baare mein apply hota hai jabki woh pehle se p = 4 , q = 5 rad/s spin kar raha hai. Nikalo r ˙ .
Forecast: kya p , q spins z ke baare mein koi coupling torque produce karte hain? Yes/no guess karo.
Euler-z : M z = I z r ˙ + ( I y − I x ) pq .
Yeh step kyun? Teeno inertias aate hain; coupling bracket dekho.
Kyunki I x = I y , bracket ( I y − I x ) = 0 ; coupling term bilkul khatam ho jaata hai.
Yeh step kyun? Sphere ka koi "hard" ya "easy" spin axis nahi hota, isliye kabhi koi gyroscopic coupling nahi hoti.
Isliye M z = I r ˙ ⇒ r ˙ = M z / I = 3/12 = 0.25 rad/s².
Yeh step kyun? Rotation ab bilkul translation ki tarah behave karta hai: torque = inertia × angular accel.
Verify karo: sphere ke liye har Euler equation M i = I ω ˙ i tak reduce ho jaata hai —
teen independent scalar laws, koi cross-talk nahi. Units: N·m / kg·m² = 1/s² = rad/s² ✓.
Worked example Nose thruster se clean pitch-up
Ek rocket rotation mein rest par hai (p = q = r = 0 ), ek side thruster F = 800 N fire karta hai
CoM se ℓ = 3 m aage lever arm par, ek pure pitch moment produce karta hai. I y = 500 kg·m².
Nikalo pitch acceleration q ˙ aur t = 2 s burn ke baad pitch rate.
Forecast: kya yahan roll/yaw state matter karta hai? (Nahi karega — isliye humne p = r = 0 choose kiya.)
Pitch axis ke baare mein moment: M y = F ℓ = 800 × 3 = 2400 N·m.
Yeh step kyun? Force times lever arm = torque; yeh drive term hai.
Euler-y : M y = I y q ˙ + ( I x − I z ) r p . r = p = 0 ke saath coupling khatam ho jaati hai.
Yeh step kyun? Non-spinning shuru karne ka matlab maneuver ek clean axis mein decouple ho jaata hai.
q ˙ = M y / I y = 2400/500 = 4.8 rad/s².
Yeh step kyun? Torque ko inertia se divide karo — M y = I y q ˙ ka final rearrangement.
Burn ke baad: q = q ˙ t = 4.8 × 2 = 9.6 rad/s.
Yeh step kyun? Kyunki thruster force, lever arm aur inertia sab burn ke dauran constant hain,
q ˙ constant hai, aur ek constant rate of change simply
q = q ˙ t (rate × time) tak integrate ho jaata hai, exactly jaise distance = speed × time.
Verify karo: units N·m / kg·m² = rad/s² ✓ q ˙ ke liye, phir rad/s² · s = rad/s ✓ q = 9.6
rad/s ke liye. Sanity check: agar burn sirf t = 0 s chala to pitch rate 0 hogi, aur yeh
burn time ke saath linearly badhti hai — jaise constant acceleration mein expect hota hai. Agar rocket roll kar raha hota
(p = 0 ), Case C dikhata hai ki ek cross term aa jaata — isliye "start from rest" trick kaam aati hai.
Worked example Off-axis thruster ke saath gravity turn
Ek rocket (m = 1000 kg, I z = 400 kg·m²) horizontally fly kar raha hai: body-x forward u = 150
m/s, v = w = 0 , aur yeh yaw kar raha hai r = 0.2 rad/s par (p = q = 0 ). Ek gimballed engine
thrust T = 12 , 000 N body-x ke along deta hai, plus ek chhota yaw-control side force F s = 500 N
CoM se ℓ = 4 m peeche lever arm par. Gravity neeche act karti hai (body-z yahan neeche point karta hai) m g par, g = 9.8 .
Nikalo (a) forward accel u ˙ , (b) sideways accel v ˙ , (c) yaw accel r ˙ .
Forecast: teen answers mein se kaun sa simply force/mass nahi hai? (Step 1 se pehle guess karo.)
Forward: F x = T = 12000 N. Newton-x : F x = m ( u ˙ + q w − r v ) , aur q = v = 0 isliye
u ˙ = F x / m = 12000/1000 = 12 m/s².
Yeh step kyun? Yahan x par koi transport term survive nahi karta (use q ya v chahiye).
Sideways: side force F y = F s = 500 N. Newton-y : F y = m ( v ˙ + r u − pw ) ,
p = 0 ke saath: v ˙ = F y / m − r u = 500/1000 − ( 0.2 ) ( 150 ) = 0.5 − 30 = − 29.5 m/s².
Yeh step kyun? Yahi non-trivial wala hai — yaw velocity arrow ko swing karta hai, ek
bada transport term − r u add karta hai jo thruster ke direct contribution ko peeche chhod deta hai. Hum v ˙ isolate karte hain
us transport term ko right-hand side par move karke.
Yaw: moment M z = F s ℓ = 500 × 4 = 2000 N·m. Euler-z :
M z = I z r ˙ + ( I y − I x ) pq , p = q = 0 ke saath: r ˙ = M z / I z = 2000/400 = 5 rad/s².
Yeh step kyun? Rotation decouple ho jaata hai kyunki p = q = 0 ; sirf drive term M z survive karta hai, isliye
promised yaw acceleration seedha moment over inertia se aata hai.
Verify karo: (a) u ˙ = 12 m/s², (b) v ˙ = − 29.5 m/s², (c) r ˙ = 5 rad/s². v ˙ ka
sign negative hai: bhale hi thruster + y push karta ho, turn ki centripetal demand body ko doosri taraf zyada
tezi se sideways pull karti hai — ek real coordinated-turn effect. Units sab m/s² aur rad/s² ✓.
d t d ( m v ) kyun fail ho jaata hai
Ek rocket propellant burn karta hai: m ˙ = − 20 kg/s, current m = 800 kg, rocket ke relative
exhaust speed u e = 2500 m/s. Yeh seedha fly karta hai, ω = 0 , no gravity, u = 400 m/s.
Ek student likhta hai F = d t d ( m u ) = m u ˙ + m ˙ u aur only force ko zero set karta hai.
Nikalo sahi u ˙ aur error expose karo.
Forecast: student conclude karta hai u ˙ = − m ˙ u / m . Kya woh real acceleration hai?
Sahi physics (Meshchersky ):
F e x t + m ˙ u e , rel = m u ˙ . Thrust T = − m ˙ u e ko external force treat karo.
Yeh step kyun? Ejected gas momentum carry karta hai; m ˙ u ko d t d ( m v ) mein lump karke
Newton nahi keh sakte — reference momentum gas ka hai, rocket ka nahi.
Thrust: T = − m ˙ u e = − ( − 20 ) ( 2500 ) = 50 , 000 N (forward).
Yeh step kyun? Mass peeche se jaata hai (m ˙ < 0 ), isliye thrust forward aur positive hai.
Koi aur external force nahi, isliye m u ˙ = T ⇒ u ˙ = 50000/800 = 62.5 m/s².
Yeh step kyun? Thrust ko external force maante hue clean constant-mass Newton form.
Student ka galat answer: u ˙ wrong = − m ˙ u / m = − ( − 20 ) ( 400 ) /800 = 10 m/s²
— six se zyada factor se galat, aur exhaust speed u e ki jagah flight speed u use karta hai.
Yeh step kyun? Trap vehicle momentum se exhaust momentum replace kar deta hai.
Verify karo: sahi u ˙ = 62.5 m/s² vs naive 10 m/s². Missing physics ka ratio
u e / u = 2500/400 = 6.25 hai, exactly discrepancy — proof hai ki error swapped velocity hai ✓.
Units: N/kg = m/s² ✓.
Recall Har case ka ek-line summary
Case A ::: no spin → plain F = ma .
Case B ::: constant body speed lekin turning → ω × v term se real accel.
Case C ::: ( I x − I z ) r p ka sign sets karta hai ki gyroscopic torque kis direction mein point karega.
Case D ::: symmetric spinner → bounded coning frequency ∣ λ ∣ = I t ∣ I x − I t ∣ p 0 par.
Case E ::: intermediate-axis spin → coefficient k > 0 → exponential flip.
Case F ::: sphere → sab coupling zero → teen independent M i = I ω ˙ i .
Case G ::: rest se shuru → single-axis maneuver cleanly decouple ho jaata hai.
Case H ::: word problem → F , M build karo, dekho transport term v ˙ par dominate karta hai.
Case I ::: variable mass → thrust ek external force hai − m ˙ u e , kabhi d t d ( m v ) nahi.
"Zero kills the cross, spin wakes the coupling." Agar koi rate zero hai, uske cross terms vanish ho jaate hain
(Cases A, G). Agar inertias tie kar jaati hain, coupling vanish ho jaati hai (Cases D on spin axis, F). Warna
ω × terms alive hain aur unhe carry karna padega.
In Ex 2, F y nonzero kyun hai bhale hi v ˙ = 0 ho? Yaw rate r velocity arrow ko swing karta hai, isliye transport term r u real inertial (centripetal) acceleration r u = u 2 / R deta hai.
Fully asymmetric body ke liye kaun sa spin axis unstable hota hai? Intermediate (middle) axis — coupling coefficient k positive ho jaata hai, giving exponential growth (Dzhanibekov effect).
Perfect sphere ke liye Euler ke coupling terms ka kya hota hai? Sab vanish ho jaate hain kyunki I x = I y = I z , teen independent scalar laws M i = I ω ˙ i bacha jaate hain.
Burning rocket ke liye naive d t d ( m v ) kyun galat hai? Yeh exhaust momentum ki jagah vehicle momentum use karta hai; sahi thrust external force − m ˙ u e hai (Meshchersky).