2.1.20 · D3 · Physics › Analytical Mechanics › Normal modes — coupled oscillators, normal coordinates
Yeh page parent note on normal modes ke liye practice arena hai. Wahan humne machinery derive ki thi. Yahan hum use har tarah ki situation ke khilaf run karte hain jo yeh topic dे sakta hai — clean modes, energy sloshing, weird masses, zero coupling, infinite coupling, ek real-world word problem, aur ek exam trap.
Numbers chhune se pehle, puri landscape lay out karte hain taaki koi bhi case surprise na kare.
Intuition "Har scenario" se yahan kya matlab hai
Ek coupled-oscillator problem do cheezein fix karta hai: the system (masses, spring constants) aur how you kick it (initial positions & velocities). Neeche har example (system type) × (kick type) ke grid mein ek point hai. Grid cover karo = topic cover karo.
Matrix mein do symbols aayenge, toh pehle inhe plain words mein name karte hain:
Definition Mass matrix aur identity matrix
M mass matrix hai — ek chhota grid jo har mass ko apne diagonal par rakhta hai, jaise ( m 1 0 0 m 2 ) . Yeh equations ko batata hai ki har coordinate kitna heavy hai. Symbol 1 identity matrix hai, ek grid jisme diagonal par 1 's hain aur baaki jagah 0 's hain (isse multiply karne se kuch nahi badalta). Jab dono masses equal hote hain, M = m 1 — ek single number m times woh plain identity. "Unequal masses" ka matlab hai M = m 1 : diagonal ab ek repeated number nahi hai.
#
Case class
Kya special hai
Example
A
Pure symmetric kick
initial condition is eigenvector 1 → ek frequency
Ex 1
B
Pure antisymmetric kick
initial condition is eigenvector 2 → ek frequency
Ex 2
C
Mixed kick (beats)
general start → superposition, energy sloshes
Ex 3
D
Weak-coupling limit k c ≪ k
modes nearly degenerate → bahut slow beats
Ex 4
E
Zero coupling k c = 0
degenerate case: coupling vanish ho jaata hai, masses independent
Ex 5
F
Strong / stiff-link limit k c ≫ k
ek mode race karta hai, link rigid act karta hai
Ex 6
G
Unequal masses M = m 1
generalized eigenproblem, M-orthogonality bites
Ex 7
H
Real-world word problem (units!)
do atoms / do carts, Hz mein real frequency nikalo
Ex 8
I
Velocity kick + timing twist
struck (pull nahi) start, aur full-transfer time
Ex 9
Nau examples, nau cells. Note karo ki "kick" ek displacement ho sakta hai (mass ko pull karke release karo) ya ek velocity (rest mein mass ko strike karo) — Cell I velocity case cover karta hai taaki initial conditions ka poora space exhaust ho jaye. Chalo shuru karte hain.
Throughout, system parent ka two-mass chain hai (jab tak cell kuch aur na kahe), jinke results hum reuse karenge:
Definition Potential energy
V aur K kahan se aata hai
V potential energy hai — stretched/compressed springs mein stored energy. Hamare chain ke liye V = 2 1 k x 1 2 + 2 1 k x 2 2 + 2 1 k c ( x 2 − x 1 ) 2 . Stiffness matrix K yeh dekh ke banti hai ki har coordinate kitni strongly push back karta hai: entry K ij = ∂ 2 V / ∂ x i ∂ x j , equivalently coordinate i par restoring force − ∂ V / ∂ x i hai. Toh K bas spring stiffnesses ko ek grid mein arrange karna hai.
Worked example Dono masses ek hi direction mein pull kiye
System: m = 1 kg , k = 4 N/m , k c = 3 N/m .
Kick: x 1 ( 0 ) = x 2 ( 0 ) = 0.05 m , dono rest mein.
Find: motion, aur woh frequency jo tum measure karoge.
Forecast: Padhne se pehle guess karo — kya tum beats dekhoge, ya ek clean tone? (Ek clean tone: kick exactly eigenvector 1 hai.)
t = 0 par normal coordinates compute karo.
η 1 ( 0 ) = 2 1 ( 0.05 + 0.05 ) = 2 0.10 , η 2 ( 0 ) = 2 1 ( 0.05 − 0.05 ) = 0 .
Yeh step kyun? Normal coordinates hi woh ek maatra variables hain jo independently move karte hain. t = 0 par inhe padhne se pata chalta hai ki kaun se modes switched on hain.
Note karo ki η 2 ≡ 0 hamesha rahega. Uska equation η ¨ 2 = − ω 2 2 η 2 hai jisme η 2 ( 0 ) = 0 aur η ˙ 2 ( 0 ) = 0 → zero hi rahega.
Yeh step kyun? Ek oscillator jo origin par rest se start ho woh kabhi wahan se nahi jaata. Mode 2 dead hai.
Sirf mode 1 alive hai: η 1 ( t ) = 2 0.10 cos ω 1 t jisme ω 1 = k / m = 4 = 2 rad/s .
Yeh step kyun? Ek surviving normal coordinate η ¨ 1 = − ω 1 2 η 1 obey karta hai, jiska zero-velocity solution usi mode ki apni frequency par ek pure cosine hota hai.
Back-transform: x 1 = x 2 = 2 1 η 1 = 0.05 cos ( 2 t ) m .
Yeh step kyun? x j = 2 1 ( η 1 ± η 2 ) ; η 2 = 0 ke saath dono masses ek motion share karte hain.
Verify: Dono masses identically move karte hain, toh middle spring kabhi stretch nahi hoti — usse ω affect nahi honi chahiye. Indeed ω = k / m mein koi k c nahi hai. ✓ Units: ( N/m ) / kg = s − 2 = s − 1 . ✓ (Cell A)
Worked example Masses opposite directions mein pull kiye
Same system (m = 1 , k = 4 , k c = 3 ).
Kick: x 1 ( 0 ) = + 0.05 , x 2 ( 0 ) = − 0.05 , rest mein.
Find: frequency.
Forecast: Case A se faster ya slower? (Faster — coupling spring ab fight kar rahi hai.)
Normal coordinates: η 1 ( 0 ) = 2 1 ( 0.05 − 0.05 ) = 0 , η 2 ( 0 ) = 2 1 ( 0.05 + 0.05 ) = 2 0.10 .
Yeh step kyun? Same reason as A: find karo kaun sa mode lit hai. Ab woh mode 2 hai.
Mode 1 dead, mode 2 alive: ω 2 = ( k + 2 k c ) / m = ( 4 + 6 ) /1 = 10 ≈ 3.162 rad/s .
Yeh step kyun? Antisymmetric stretch coupling spring ka effect double kar deta hai: k + 2 k c , na ki k + k c .
Motion: x 1 = − x 2 = 0.05 cos ( 10 t ) m .
Yeh step kyun? Back-transform x j = 2 1 ( η 1 ± η 2 ) jisme η 1 = 0 single frequency ω 2 par equal-and-opposite motion deta hai.
Verify: ω 2 > ω 1 (3.162 > 2 ) ✓ — mnemonic SLOW=Same, FAST=Fight se match karta hai. Factor-2 trap k + 2 k c use karke dodge kiya. ✓ (Cell B)
Worked example Ek hi mass pull karo
Same system (m = 1 , k = 4 , k c = 3 ).
Kick: x 1 ( 0 ) = a = 0.05 , x 2 ( 0 ) = 0 , rest mein.
Find: x 1 ( t ) , aur motion describe karo.
Forecast: Kya mass 2, jo still start karta hai, still rahega? (Nahi — energy flood hoti hai; beats milte hain.)
Normal coords: η 1 ( 0 ) = 2 a , η 2 ( 0 ) = 2 a — dono modes lit, equally.
Yeh step kyun? Ek one-sided pull symmetric aur antisymmetric ka 50/50 mix hai; wahi mix sloshing produce karta hai.
Har ek evolve karo: η 1 = 2 a cos ω 1 t , η 2 = 2 a cos ω 2 t , jisme ω 1 = 2 , ω 2 = 10 .
Yeh step kyun? Har normal coordinate ek independent SHM hai; rest se start kiya toh woh bas apni frequency par cosines karta hai, aur ab dono frequencies alag hain.
Back-transform karo & sum-to-product identity use karo:
x 1 = 2 a ( cos ω 1 t + cos ω 2 t ) = a cos ( 2 ω 2 − ω 1 t ) cos ( 2 ω 1 + ω 2 t )
Yeh step kyun? Identity isse split karta hai ek fast carrier (average frequency) ke andar ek slow envelope (half the difference) mein — yahi exactly beats hote hain.
Envelope frequency: 2 ω 2 − ω 1 = 2 3.162 − 2 ≈ 0.581 rad/s .
Yeh step kyun? Yahi slow rate hai jo tumhari aankh ko actually energy "sloshing" back and forth dikhti hai; yeh set karta hai ki ek full transfer kitna time leta hai (Case I mein phir use hoga).
Neeche ka figure yeh exact motion plot karta hai. Horizontal axis par time t (s) hai aur vertical axis par x 1 (units of a mein) hai, legend mein har curve ka naam hai. Magenta curve dekho jisko "x1(t) full motion" label kiya hai: yeh ek fast wiggle hai jo violet dashed envelope ke andar trapped hai jisko "slow envelope" label kiya hai, equal to ± a cos ( 2 ω 2 − ω 1 t ) . Jahan envelope zero par pinch hota hai (orange dotted line ), mass 1 quiet ho jaata hai aur saari energy mass 2 par cross ho jaati hai — woh beat hai.
Verify: t = 0 par, x 1 = a ⋅ 1 ⋅ 1 = a ✓ aur x 2 = 2 a ( cos ω 1 t − cos ω 2 t ) = 0 ✓. Energy periodically mass 1 se bilkul leave ho jaati hai (envelope zero hit karta hai). ✓ (Cell C)
Worked example Loosely tied swings
System: m = 1 , k = 100 , k c = 1 . Case C jaisa one-sided pull.
Find: beat period ek single swing period se kaise compare karta hai.
Forecast: Energy exchange fast hai ya slow? (Bahut slow — weak link ka matlab hai modes almost same pitch par hain.)
Frequencies: ω 1 = 100 = 10 , ω 2 = 102 ≈ 10.0995 .
Yeh step kyun? Humein dono eigenfrequencies chahiye inki difference form karne ke liye; tiny k c ke saath woh almost coincide karte hain, yahi is limit ka poora point hai.
Envelope frequency: 2 ω 2 − ω 1 ≈ 2 0.0995 ≈ 0.0498 rad/s .
Yeh step kyun? Beat rate frequency difference hai; tiny difference ⇒ tiny beat rate ⇒ lamba slosh time.
Beat period T beat = ∣ ω 2 − ω 1 ∣ 2 π = 0.0995 2 π ≈ 63.1 s , versus ek single swing period ≈ 10 2 π ≈ 0.628 s .
Yeh step kyun? Ek full slosh over and back envelope ka ek period leta hai, aur envelope difference frequency par oscillate karta hai — toh ∣ ω 2 − ω 1 ∣ se 2 π divide karna us rate ko ek actual time mein convert karta hai jo ek single swing se compare ho sake.
Verify: Beat period ≫ swing period (63 ≫ 0.63 , ratio ≈ 100 ) — weak coupling ⇒ har slow slosh mein bahut saare fast swings, loosely coupled oscillators ki hallmark. ✓ (Cell D)
Worked example Middle spring kaat do
System: m = 1 , k = 4 , k c = 0 . Kick: x 1 ( 0 ) = 0.05 , x 2 ( 0 ) = 0 , rest mein.
Find: har mass kya karta hai.
Forecast: Koi link nahi, toh kya energy mass 2 par cross hogi? (Nahi — ab woh independent hain.)
Frequencies collapse together: ω 1 2 = k / m = 4 , ω 2 2 = ( k + 0 ) / m = 4 . Dono = 2 rad/s — degenerate .
Yeh step kyun? k c = 0 ke saath dono modes ki same frequency hai; system ab truly coupled nahi hai.
Beat envelope frequency = 2 ω 2 − ω 1 = 0 → envelope kabhi nahi dipti.
Yeh step kyun? Zero frequency difference ka matlab hai koi sloshing nahi; mass 2 ko kabhi energy nahi milti.
Motion: x 1 = 0.05 cos ( 2 t ) , x 2 = 0 har time ke liye.
Yeh step kyun? Ek baar decouple hone par, mass 2 m x ¨ 2 = − k x 2 obey karta hai jisme x 2 ( 0 ) = x ˙ 2 ( 0 ) = 0 , jiska ek hi solution x 2 ≡ 0 hai; mass 1 ek lone cosine oscillator ki tarah ω 1 par reh jaata hai.
Verify: x 2 kabhi 0 se depart nahi karta — equations instant mein decouple ho gaye jab k c → 0 : m x ¨ 2 = − k x 2 jisme x 2 ( 0 ) = x ˙ 2 ( 0 ) = 0 . ✓ Yeh sanity anchor hai: zero coupling zaroori hai independent SHM dene ke liye. (Cell E)
Worked example Masses ko weld kar do
System: m = 1 , k = 4 , k c = 1000 .
Find: do frequencies aur fast mode ka kya "matlab" hai.
Forecast: Jab k c → ∞ , ω 2 ka kya hoga? (Yeh infinity tak run karta hai — antisymmetric motion impossibly stiff ho jaata hai.)
Symmetric mode untouched rehta hai: ω 1 = k / m = 2 rad/s , chahe k c kitna bhi large ho.
Yeh step kyun? Symmetric mode mein coupling spring kabhi stretch nahi hoti, toh uski stiffness irrelevant hai.
Antisymmetric mode: ω 2 = ( k + 2 k c ) / m = 2004 ≈ 44.77 rad/s .
Yeh step kyun? Masses ki koi bhi relative motion bahut stiff link ko stretch karti hai, enormous energy cost karta hai — hence bahut high frequency.
Physical reading: slow motions ke liye sirf ω 1 excitable hai; pair ek rigid block ki tarah outer springs par 2 m mass ke saath move karta hai. Stiff mode "freeze out" ho jaata hai.
Yeh step kyun? Kyunki ω 2 itna far above ω 1 hai, koi bhi gentle, slow driving sirf symmetric mode stir kar sakta hai — antisymmetric ek ke liye frequencies itni high chahiye jo reach nahi ho sakti, toh effectively dono masses ek single rigid object mein lock ho jaate hain aur sirf mode 1 observable rehta hai.
Verify: ω 1 2 par raha (k c se independent) ✓; ω 2 ≫ ω 1 (44.8 ≫ 2 ) ✓. Yeh freeze-out exactly waise hai jaise high-frequency optical branches acoustic ones se alag hote hain Phonons and Lattice Vibrations mein. (Cell F)
Worked example Heavy mass, light mass
System: m 1 = 1 , m 2 = 2 (kg), har ek wall se spring k = 6 se tied, coupled by k c = 6 (N/m). Ab M = ( 1 0 0 2 ) nahi hai m 1 .
Find: do frequencies.
Forecast: Kya eigenvectors phir bhi ( 1 , 1 ) aur ( 1 , − 1 ) honge? (Nahi — unequal masses modes ko tilt kar dete hain; simple symmetric/antisymmetric shapes break ho jaati hain.)
Matrices likho. Stiffness (potential energy V se, − ∂ V / ∂ x ke through):
K = ( k + k c − k c − k c k + k c ) = ( 12 − 6 − 6 12 )
Yeh step kyun? K potential energy se aata hai aur masses par depend nahi karta.
Generalized secular equation solve karo det ( K − ω 2 M ) = 0 :
det ( 12 − ω 2 − 6 − 6 12 − 2 ω 2 ) = ( 12 − ω 2 ) ( 12 − 2 ω 2 ) − 36 = 0.
Yeh step kyun? Yeh generalized hai — ω 2 M se multiply hota hai, identity se nahi. Yeh woh trap hai jo parent ke mistake box mein flag kiya gaya tha.
Expand karo: 2 ω 4 − 36 ω 2 + 108 = 0 ⇒ ω 4 − 18 ω 2 + 54 = 0 . ω 2 mein quadratic:
ω 2 = 2 18 ± 324 − 216 = 9 ± 27 = 9 ± 3 3 .
Toh ω 2 ≈ 14.196 ya 3.804 , yaani ω ≈ 3.768 aur ω ≈ 1.950 rad/s .
Yeh step kyun? Determinant multiply out karne se matrix condition ek ordinary polynomial ban jaati hai; ω 2 ko ek unknown treat karne se yeh ek plain quadratic ban jaata hai jo hum standard formula se solve kar sakte hain, aur har positive root ek squared mode frequency hai.
Verify: Dono ω 2 > 0 (stable) ✓. Roots ka product = ( 9 ) 2 − ( 3 3 ) 2 = 81 − 27 = 54 = det K / det M = 108/2 ✓. Eigenvectors ( 1 , ± 1 ) nahi hain aur orthogonal sirf M ke saathe hain — dekho Eigenvalues and Eigenvectors . (Cell G)
Worked example Do coupled pendulum-carts
Do carts, har ek m = 0.50 kg , har ek wall se spring k = 8.0 N/m se anchored aur k c = 2.0 N/m se coupled. Hertz mein do natural frequencies nikalo.
Forecast: Dono kuch Hz? (Haan — light carts, moderate springs.)
Symmetric: ω 1 = k / m = 8.0/0.50 = 16 = 4.0 rad/s .
Yeh step kyun? Woh mode use karo jahan link so rahi ho.
Antisymmetric: ω 2 = ( k + 2 k c ) / m = ( 8.0 + 4.0 ) /0.50 = 24 ≈ 4.899 rad/s .
Yeh step kyun? Opposite-way motion coupling spring ko work karta hai, effective stiffness mein 2 k c add karta hai.
Hz mein convert karo f = ω /2 π use karke: f 1 = 4.0/ ( 2 π ) ≈ 0.637 Hz , f 2 ≈ 4.899/ ( 2 π ) ≈ 0.780 Hz .
Yeh step kyun? Real instruments cycles per second padhte hain, rad/s nahi; angular frequency ko 2 π (radians per cycle) se divide karna hertz mein convert karta hai.
Verify: Units: ( N/m ) / kg = s − 1 ✓; 2 π (rad/cycle) se divide karna Hz deta hai ✓. f 2 > f 1 ✓. (Cell H)
Worked example Ek mass strike karo, phir hand-off time karo
System (m = 1 , k = 4 , k c = 3 ). Ab kick ek velocity hai, pull nahi: dono masses equilibrium par baithe hain, aur mass 1 ko strike kiya jaata hai taaki x ˙ 1 ( 0 ) = v 0 = 0.30 m/s , jisme x 1 ( 0 ) = x 2 ( 0 ) = 0 aur x ˙ 2 ( 0 ) = 0 .
Find: (a) x 1 ( t ) ; (b) kitna time lagega jab mass 1 momentarily quiet ho aur mass 2 full swing carry kare.
Forecast: Zero initial displacement ke saath, kya phir bhi beats milte hain? (Haan — velocity kick bhi dono modes excite karta hai; ab sine beats milte hain cosine ki jagah.)
Normal coordinates aur unke velocities t = 0 par.
η 1 ( 0 ) = η 2 ( 0 ) = 0 ; η ˙ 1 ( 0 ) = 2 1 ( x ˙ 1 + x ˙ 2 ) = 2 v 0 , η ˙ 2 ( 0 ) = 2 1 ( x ˙ 1 − x ˙ 2 ) = 2 v 0 .
Yeh step kyun? Ek SHM jo origin par speed ke saath start hota hai woh ek sine hai, cosine nahi — toh humein positions nahi, initial velocities padhni chahiye har mode ka amplitude set karne ke liye.
Har mode ko sine ke roop mein evolve karo. η ¨ a = − ω a 2 η a ke liye jisme η a ( 0 ) = 0 , η ˙ a ( 0 ) = u : η a ( t ) = ω a u sin ω a t . Toh
η 1 = 2 ω 1 v 0 sin ω 1 t , η 2 = 2 ω 2 v 0 sin ω 2 t , jisme ω 1 = 2 , ω 2 = 10 .
Yeh step kyun? Sine solution woh hai jo "origin par, move karte hue" match karta hai; ω a se divide karna units correct rakhta hai (ek velocity ek frequency se divide = ek length).
Back-transform karo: x 1 = 2 1 ( η 1 + η 2 ) = 2 v 0 ( ω 1 s i n ω 1 t + ω 2 s i n ω 2 t ) .
Yeh step kyun? Same reconstruction rule jaise hamesha; do nearly-equal frequencies phir se ek slow beat banate hain, ab sine waves par.
Hand-off time. Is near-beat ka envelope pehli baar t = ω 2 − ω 1 π = 10 − 2 π ≈ 2.703 s par pinch karta hai.
Yeh step kyun? Full energy transfer tab hota hai jab dono modes half a cycle out of step drift kar jaate hain — exactly wahi ω 2 − ω 1 π timing jaise displacement beat mein, kyunki timing sirf frequency difference par depend karti hai, iss par nahi ki tumne kick kaisa lagaya.
Verify: t = 0 par: x 1 = 0 ✓ (struck-from-rest position match karta hai) aur x ˙ 1 ( 0 ) = 2 v 0 ( cos 0 + cos 0 ) = v 0 ✓. Transfer time = π / ( 10 − 2 ) ≈ 2.703 s , aur yeh half beat period T beat /2 ke equal hai jisme T beat = 2 π / ( 10 − 2 ) ≈ 5.406 s ✓. Toh displacement kicks aur velocity kicks same slosh clock share karte hain. (Cell I)
Recall Matrix par quick self-test
Kaun se cell mein ω 1 = ω 2 hai? ::: Cell E — zero coupling, modes degenerate hain.
Kaun se cell mein ω 2 → ∞ hota hai jab link tighten karo? ::: Cell F — strong coupling, antisymmetric mode freeze out ho jaata hai.
Cell G mein mode shapes ( 1 , ± 1 ) kyun NAHI hain? ::: Kyunki unequal masses ise ek generalized eigenproblem banate hain; modes eigenvectors hain jo M ke saathe orthogonal hain, identity ke nahi, toh tidy symmetric/antisymmetric shapes ab nahi chalti.
Kya velocity kick phir bhi beats produce karta hai? ::: Haan — Cell I; yeh dono modes ko sines ki tarah excite karta hai, same slosh timing π / ( ω 2 − ω 1 ) ke saath sine beats deta hai.
Mnemonic Matrix ek saansmein
Kick decide karta hai kaun se modes light up hote hain; system decide karta hai har ek kitna fast burns karta hai. Same-way kick → slow lone tone; opposite kick → fast lone tone; one-sided kick → beats; struck kick → sine beats; cut link → independence; weld link → freeze-out.
Deep-dive built on Lagrangian Mechanics (the EOMs), Small Oscillations (kyun linearizing legal hai), aur Eigenvalues and Eigenvectors (the modes). Intuitive story ke liye dekho Hinglish companion .