Cím: Összetett rezgések I. rész
Szerző(k):  Vigh Máté 
Füzet: 2012/január, 38 - 47. oldal  PDF  |  MathML 
Témakör(ök): Csatolt rezgések

A szöveg csak Firefox böngészőben jelenik meg helyesen. Használja a fenti PDF file-ra mutató link-et a letöltésre.

 
Összetett rezgések
I. rész
 

 

 
Bevezetés

 
A középiskolában általában nincs lehetőség a rezgőmozgások mélyebb tárgyalására, legtöbbször csak az egyszerű harmonikus rezgés ismertetésére kerül sor. Ennek oka nem csak az idő szűkében rejlik, hanem abban is, hogy a bonyolultabb összetett rezgések szokásos tárgyalása a középiskolai tananyagon jócskán túlnyúló, lineáris algebrai ismereteket igényel (mátrixok sajátértékproblémája). A környezetünkben túlnyomórészt ezekkel az összetett rezgésekkel találkozunk: a pohár koccanása, a tortazselé remegése, a kristályokban a rácsionok nyüzsgése mind-mind összetett rezgés, amelyek nem érthetők meg pusztán az egyszerű harmonikus rezgőmozgásra vonatkozó ismeretekkel.
Ez a kétrészes cikk kísérletet tesz arra, hogy nehéz matematikai apparátus nélkül, egyszerű példákon keresztül, középiskolások számára is érthető módon tárgyalja az összetett rezgéseket és a témához kapcsolódó olyan érdekes jelenségeket, mint a csatolt ingák rezgése vagy a lebegés.
 
 
Az egyszerű harmonikus rezgőmozgás

 
Elevenítsük fel először a harmonikus rezgésre vonatkozó ismereteket a következő egyszerű példán keresztül!
 
1. példa. Függőleges falból vízszintesen egy egyenes, hosszú, merev huzal áll ki. A huzalra m tömegű gyöngy van felfűzve. A gyöngyöt D direkciós erejű húzó-nyomó rugó kapcsolja a falhoz (1. ábra).
 
A kezdetben álló gyöngyöt v0 kezdeti sebességgel elindítjuk. Írjuk le a mozgást! (A gyöngy súrlódásmentesen csúszhat a huzalon.)


 

1. ábra
 

 
Megoldás. A gyöngyre ható függőleges irányú erők kiegyenlítik egymást, vízszintesen pedig csak a rugóerő hat, így a mechanika alapegyenlete az
ma(t)=-Dx(t)
alakot ölti, ahol x(t) a gyöngy egyensúlyi helyzettől mért (időtől függő) kitérése, a(t) a gyorsulása, a negatív előjel pedig arra utal, hogy a gyorsulás mindig a kitéréssel ellentétes irányú. Az egyszerűség kedvéért a fizikai mennyiségek változási gyorsaságát1 jelöljük a fizikai mennyiség fölé írt ponttal, azaz legyen a gyöngy sebessége v(t)=x˙(t), a gyorsulása pedig a(t)=x¨(t)! A mozgásegyenletet átrendezve az
x¨(t)=-ω2x(t)(1)

alakú egyenletre jutunk, ahol bevezettük az ω=D/m jelölést. Az (1) egyenlet egy másodrendű differenciálegyenlet a gyöngy kitérésének x(t) időfüggésére, melyet szerencsére nem kell tudnunk megoldani, hiszen a megoldása a jól ismert
x(t)=Asin(ωt+φ0)(2)
formát ölti, a gyöngy sebességét pedig a deriválással (vagy a mechanikai energiamegmaradásból) kapható
x˙(t)=Aωcos(ωt+φ0)(3)
kifejezés adja meg, ahol ω a rezgés fent bevezetett körfrekvenciája, A a rezgés amplitúdója (legnagyobb kitérése), φ0 pedig a rezgés kezdőfázisa. Utóbbi kettő a kezdeti feltételektől függ, vagyis a gyöngy kezdeti kitéréséből és sebességéből kiszámítható.
A konkrét példánkban kezdetben a gyöngy az egyensúlyi helyzetben volt (x(0)=0), ezért a φ0 kezdőfázis értéke zérus. A kezdősebesség nagysága
x˙(0)=Aω=v0 volt, így a mozgás amplitúdójának nagysága A=v0/ω. A feladatbeli kezdőfeltételek esetén tehát a gyöngyszem mozgását az
x(t)=v0ωsin(ωt)
összefüggés írja le.
 
 
Két, rugóval összekötött tömegpont mozgása

 
Az első példánál egy fokkal bonyolultabb a rugóval összekötött két tömegpont mozgása, amely jól ismert a fizikaversenyekre járók körében. A probléma szokásos megoldási módja a következő: A tömegközéppont egyenletesen haladó mozgását leválasztjuk úgy, hogy beleülünk a tömegközépponttal együtt mozgó koordinátarendszerbe. Itt a rugónak a rendszer tömegközéppontjával egybeeső pontja áll, ezért az egyes testek úgy mozognak, mintha a tömegközéppont által kettéosztott rugó hozzájuk közelebb eső darabja mozgatná őket. Ezzel a gondolattal meghatározható a rezgés körfrekvenciája. Itt azonban egy másik, a bonyolultabb rendszerek esetén is járható utat választunk a megoldáshoz.
 
2. példa. Egy egyenes, hosszú, vízszintes, merev huzalra m és M tömegű gyöngyök vannak felfűzve. A gyöngyöket D direkciós erejű húzó-nyomó rugó kapcsolja össze (2. ábra). A kezdetben álló rendszert az m tömegű gyöngynek adott v0 kezdeti sebességgel elindítjuk. Írjuk le a mozgást! (A gyöngyök súrlódásmentesen csúszhatnak a huzalon.)

 

2. ábra
 

 
Megoldás. Jelöljük az m tömegű gyöngy kezdeti helyzettől mért elmozdulását xm-mel, az M tömegű gyöngy elmozdulását pedig xM-mel! Az egyes testekre vízszintesen csak a rugóerő hat, így a gyöngyök mozgásegyenletei:
mx¨m=D(xM-xm),(4)Mx¨M=-D(xM-xm).
Ezek szembeötlően nem az (1) egyenlethez hasonló kifejezések, azaz a testek mozgása nem egyszerűen harmonikus rezgőmozgás. A további vizsgálódás kedvéért osszuk el az egyenletet a tömegekkel, majd adjuk hozzá az első egyenlethez a második α-szorosát (ahol α valós szám)! Az eredmény:
x¨m+αx¨M=-(Dm-αDM)(xm-xM).(5)
Vezessük be az eredeti koordináták lineáris kombinációjaként2 képzett
y(t)=xm(t)+αxM(t)
mennyiséget! A változási gyorsaság képzési szabályából következik, hogy az y(t) mennyiség ,,gyorsulása'' a gyöngyök gyorsulásaival hasonlóan egyszerű kapcsolatban van:
y¨(t)=x¨m(t)+αx¨M(t).
Az új y(t) változó bevezetésével észrevehetjük, hogy az (5) egyenlet α=M/m vagy -1 esetén
y¨(t)=-ω2y(t)
alakúra hozható, ahol
ω=Dm-αDM.
Tehát a két lehetséges lineárkombináció és a hozzájuk tartozó körfrekvencia:
y1(t)=xm(t)+MmxM(t),ω1=0,(6)y2(t)=xm(t)-xM(t),ω2=Dm+MmM.
Az y1(t) lineárkombináció egy zérus frekvenciájú rezgés egyenletét elégíti ki, vagyis y¨1(t)=0. Ez éppen az egyenletes (tehát nulla gyorsulású) mozgás egyenlete, aminek megoldása
y1(t)=y1(0)+Vt,
ahol y1(0)-t és V-t a kezdeti feltételek határozzák meg. Az y2(t) lineárkombináció egy ω2 körfrekvenciájú harmonikus rezgőmozgás egyenletét elégíti ki, így y2(t) és y˙2(t) időfüggése a (2) és (3) összefüggéseknek megfelelően
y2(t)=Asin(ω2t+φ0),y˙2(t)=Aω2cos(ω2t+φ0),
ahol az amplitúdót és a kezdőfázist most is a kezdeti feltételek határozzák meg.
Vajon milyen kapcsolatban vannak az y1(t) és y2(t) lineárkombinációk az m és M tömegű gyöngyök mozgásával?
 
Nincs más dolgunk, mint y1(t) és y2(t) ismeretében megoldani a (6) egyenletrendszert xm(t)-re és xM(t)-re. Így a gyöngyök mozgását leíró egyenletek:
xm(t)=mm+My1(t)+Mm+My2(t),(7)xM(t)=mm+My1(t)-mm+My2(t).

További kérdés, hogy van-e fizikai jelentése az y1(t) és y2(t) lineárkombinációknak? A (7) egyenletekből leolvasható, hogy
ha y1(t)0 és y2(t)=0, a két test kitérése minden időpillanatban megegyezik, így a testek rezgés nélkül tisztán haladó mozgást végeznek. Az y1(t) tehát a testek egyenletes haladó mozgását írja le.
ha y2(t)0 és y1(t)=0, a két test kitérése minden időpillanatban ellentétes irányú, nagyságuk pedig fordítottan arányos a gyöngyök tömegével, így ebben az esetben a rendszer tömegközéppontja áll. y2(t) a gyöngyök egymáshoz viszonyított rezgőmozgását írja le.

Most térjünk vissza a konkrét példánk megoldásához! Az általános megoldásokat illesztenünk kell a kezdeti feltételekhez. Kezdetben mindkét gyöngy elmozdulása zérus, tehát xm(0)=xM(0)=0, azaz a rezgés kezdőfázisa nulla. Az indítás pillanatában az m tömegű test sebessége v0, az M tömegű test pedig áll, így x˙m(0)=v0 és x˙M(0)=0. (6) felhasználásával y˙1(0)=V=v0 és y˙2(0)=v0 adódik. Az y1(t) és y2(t) lineárkombinációk időfüggése tehát:
y1(t)=v0t,(8)y2(t)=v0ω2sin(ω2t).
Ezeket beírva a (7) egyenletekbe:
xm=mm+Mv0t+Mm+Mv0ω0sin(ω2t),xM=mm+Mv0t-mm+Mv0ω0sin(ω2t).

Végül tehát előzetes várakozásainknak megfelelő eredményt kaptunk: a testek olyan összetett mozgást végeznek, melyben a rendszer tömegközéppontja (az impulzusmegmaradás értelmében) v0m/(m+M) sebességgel egyenletesen halad, valamint a testek ω2 körfrekvenciával rezegnek a tömegközéppont körül. Ha az m/M tömegarány nagyon kicsi, akkor az M tömegű gyöngy lényegében állónak tekinthető és az m tömegű test mozgására visszakapjuk az 1. példa eredményeit.
 
 
Három, rugóval összekötött tömegpont mozgása

 
Az előbb megismert módszert alkalmazzuk most három testből és két rugóból álló rendszer rezgéseinek leírására. Legyen a konkrét feladat a következő:
 
3. példa. Egy egyenes, hosszú, vízszintes, merev huzalra három azonos m tömegű gyöngy van felfűzve. A középső gyöngyöt egy-egy D direkciós erejű, egyforma húzó-nyomó rugó kapcsolja a két szélsőhöz (3. ábra).
 
A kezdetben álló rendszert a bal szélső gyöngynek adott v0 kezdeti sebességgel elindítjuk. Írjuk le a mozgást! (A gyöngyök súrlódásmentesen csúszhatnak a huzalon.)


3. ábra

 

 
Megoldás. Jelöljük a bal szélső gyöngy elmozdulását x1-gyel, a középsőét x2-vel, a jobb szélsőét pedig x3-mal! Az egyes rugók megnyúlása x2-x1 és x3-x2, így a három gyöngyre a mozgásegyenletek:
mx¨1=D(x2-x1),(9)mx¨2=-D(x2-x1)+D(x3-x2),mx¨3=-D(x3-x2).
Bevezetve az ω0=D/m jelölést és rendezve:
x¨1=ω02(x2-x1),(10)x¨2=ω02(x1-2x2+x3),x¨3=ω02(x2-x3).
Ezek az egyenletek (akárcsak a két testből álló rendszernél) most sem hasonlítanak a szokásos harmonikus rezgés mozgásegyenletére, de látni fogjuk, hogy az x1, x2 és x3 koordinátákból megfelelő lineáris kombinációval most is kikeverhető olyan y(t) mennyiség (nem is egy, hanem három), melynek időbeli változását a harmonikus rezgés (1) egyenlete írja le. Keressük tehát azokat az α és β együtthatókat, melyekkel az
y(t)=x1(t)+αx2(t)+βx3(t)(11)
módon az eredeti koordinátákból képzett lineáris kombináció időfüggése harmonikus, azaz teljesíti az
y¨(t)=-ω2y(t)(12)
egyenletet. Ehhez szorozzuk meg a (10) összefüggések közül a másodikat α-val, a harmadikat pedig β-val, majd adjuk össze az egyenleteket!
x¨1+αx¨2+βx¨3=ω02((x2-x1)+α(x1-2x2+x3)+β(x2-x3)).
Rendezve:
x¨1+αx¨2+βx¨3=ω02((α-1)x1+(1-2α+β)x2+(α-β)x3).
Látható, hogy ha α=β=1, akkor az egyenlet jobb oldala eltűnik, azaz a koordináták
y1(t)=x1(t)+x2(t)+x3(t)
lineárkombinációja egy ω1=0 körfrekvenciájú y¨1(t)=0 alakú ,,rezgési egyenletet'' elégít ki. Ez azt jelenti, hogy az y1(t) mennyiség az idő függvényében egyenletesen változik, azaz
y1(t)=y1(0)+Vt,
ahol V egy sebesség dimenziójú állandó (lényegében a testek kezdősebességeinek összege). Hasonlóan a 2. példa megoldásához, látni fogjuk, hogy ez a lineárkombináció írja le a rendszer tömegközéppontjának mozgását.
A továbbiakban tegyük fel, hogy α1! Azért, hogy a (12) egyenlettel azonos alakú időfüggést kapjunk, végezzük el a következő kiemelést:
x¨1+αx¨2+βx¨3y¨(t)=-(1-α)ω02ω2(x1+1-2α+βα-1x2+α-βα-1x3)y(t).(13)
Innen leolvasható, hogy a harmonikus időfüggés teljesüléséhez az együtthatók között a következő összefüggéseknek kell fennállniuk:
α=1-2α+βα-1,β=α-βα-1,(14)
a körfrekvencia pedig (13) szerint
ω=ω01-α.(15)

A (14) egyenletrendszer így is írható:
(α+2)(α-1)=β-1,α(β-1)=0.
Ennek két megoldása van:
α=0,β=-1;α=-2,β=1.
A korábban talált y1(t) mellett megkaptuk tehát az x1, x2, x3 koordináták további két olyan lineárkombinációját, amelyek az időben (12) szerint harmonikusan változnak. Tehát a három független lineárkombináció (és a nekik megfelelő rezgési frekvencia) a következő:
y1(t)=x1(t)+x2(t)+x3(t),ω1=0,(16)y2(t)=x1(t)-x3(t),ω2=ω0,y3(t)=x1(t)-2x2(t)+x3(t),ω3=3ω0.
Ezek időfüggése pedig:
y1(t)=y1(0)+Vt,(17)y2(t)=A2sin(ω2t+φ2),y3(t)=A3sin(ω3t+φ3).
Az y1(0) és V értéke, az A2, A3 amplitúdók és a φ2, φ3 kezdőfázisok a kezdeti feltételekből (azaz a kiinduló koordinátákból és sebességekből) határozhatók meg. A feladatunk most is az x1(t), x2(t) és x3(t) kitérés-idő függvények meghatározása. A (16) egyenlet alapján ezek kifejezhetők az ismert időfüggésű y1(t), y2(t), y3(t) mennyiségekkel:

x1(t)=13y1(t)+12y2(t)+16y3(t),(18)x2(t)=13y1(t)-13y3(t),x3(t)=13y1(t)-12y2(t)+16y3(t).

A két testből álló rendszerhez hasonlóan könnyen kideríthetjük, hogy mi az yi(t) (i=1,2,3) lineárkombinációk fizikai jelentése (4. ábra). A (18) összefüggésből látható, hogy ha
y1(t)0 és y2(t)=y3(t)=0, akkor x1(t)=x2(t)=x3(t), azaz mindhárom test kitérése megegyezik: ilyenkor a rendszer nem rezeg, ez a lineárkombináció a tisztán haladó mozgást írja le.
y2(t)0 és y1(t)=y3(t)=0, akkor x1(t)=-x3(t) és x2(t)=0, vagyis az y2(t) lineárkombináció olyan ,,lélegző'' mozgásformát ír le, melyben a középső test áll, a két szélső test kitérése pedig minden időpillanatban ellentétes, egyenlő nagyságú. A rendszer tömegközéppontja ilyenkor áll.
y3(t)0 és y1(t)=y2(t)=0, akkor x1(t)=x3(t)=-12x2(t), ami olyan ,,tili-toli'' mozgásformát jelent, hogy a két szélső test kitérése minden időpillanatban megegyezik, a középső pedig velük ellentétesen, kétszer akkora kitéréssel mozog. A tömegközéppont ebben az esetben sem mozog.


4. ábra.Az yi(t) (i=1,2,3) lineárkombinációk fizikai jelentésének szemléltetése: y1(t) az egyenletes mozgást (a), y2(t) a lüktető mozgást (b), y3(t) pedig a ,,tili-toli'' mozgást (c) írja le. A rendszer mozgása e három mozgásforma kezdőfeltételektől függő, megfelelő arányú szuperpozíciójaként (keverékeként) állítható elő

 

Hátra van még a konkrét példánk megoldása, azaz a (18) általános megoldások illesztése a kezdeti feltételekhez. Kezdetben mindhárom test egyensúlyi helyzetben volt, így x1(0)=x2(0)=x3(0)=0. Csak a bal oldali testet löktük meg v0 kezdősebességgel, így x˙1(0)=v0, x˙2(0)=x˙3(0)=0. Ezekből és a (16) egyenletekből következik, hogy y1(0)=y2(0)=y3(0)=0 (azaz φ2=φ3=0) és y˙1(0)=y˙2(0)=y˙3(0)=v0. A (17) összefüggéseket figyelembe véve V=v0, A2ω2=v0, A3ω3=v0. A testek mozgását leíró egyenletek tehát:

x1(t)=13v0t+12v0ω0sin(ω0t)+16v03ω0sin(3ω0t),(19)x2(t)=13v0t-13v03ω0sin(3ω0t),x3(t)=13v0t-12v0ω0sin(ω0t)+16v03ω0sin(3ω0t).

Az egyes gyöngyök elmozdulását ábrázolhatjuk az idő függvényében, ahogy azt az 5. ábra mutatja az yi(t) lineárkombinációk időfüggésével együtt.

5. ábra. Első oszlop: A három test xi (i=1,2,3) kitérése az idő függvényében.
Második oszlop: az xi(t) kitérésekből képzett yi(t) lineárkombinációk.
Az időt T0=2πω0 egységekben mértük

 

 
 
Általános következtetések

 
Az előzőekben láttuk, hogy a két és három testből álló rendszerek egyenes mentén történő (egydimenziós) rezgéseinél az egyes testek xi(t) koordinátáinak időfüggését leíró egyenletek általában nem hasonlítanak az egyszerű harmonikus rezgőmozgás egyenletére, de a koordinátákból megfelelő lineárkombináció képzésével előállíthatók olyan yi(t) mennyiségek, melyek már az időben harmonikusan (egy jellemző körfrekvenciával tisztán szinuszosan) változnak. Az yi(t) lineárkombinációk pontos időfüggése a rendszer kezdőfeltételeinek (kezdeti kitérések és sebességek) ismeretében meghatározható, majd abból az egyes testek elmozdulásainak xi(t) értéke is megállapítható.
Az eljárás általánosítható nagyobb rendszerekre is. Álljon a rendszerünk N tetszőleges tömegű pontszerű testből, melyek egy egyenes mentén mozoghatnak, és a kössék össze a testeket (akár különböző erősségű) rugók! Ha az i-edik test kitérését xi(t)-vel jelöljük, a rendszer mozgása a következő recepttel írható le:
1.Felírjuk a rendszer összes elemének mozgásegyenletét, azaz minden i=1,2,...,N-re egy
x¨i(t)=j=1Nηjxj(t)
alakú egyenletet (itt az ηj együtthatókat a rugók direkciós ereje és a testek tömegei határozzák meg).
2.Megkeressük azokat az α1,α2,...,αN együtthatókat, melyekkel a testek kitéréseiből képzett
y(t)=j=1Nαjxj(t)
lineárkombinációk a harmonikus rezgőmozgás
y¨(t)=-ω2y(t)
egyenletét elégítik ki alkalmas ω körfrekvenciával. Általában N különböző ilyen yi(t) lineárkombinációt fogunk találni (ahol i=1,2,...,N). Ezekből egy kombinációhoz mindig zérus körfrekvencia tartozik, ez írja le az egyenletesen haladó mozgást.
3.Az yi(t) lineárkombinációkra vonatkozó rezgési egyenleteket megoldjuk és a megoldásokat illesztjük a kezdeti feltételekhez (kezdősebességek és kezdeti koordináták).
4.Végül az egyes testek xi(t) koordinátáit kifejezzük az yi(t) mennyiségekkel.

Az yi(t) mennyiségeket, vagy általánosabban fogalmazva a rendszer eredeti koordinátáinak olyan lineárkombinációit, amelyek időfüggése tisztán harmonikus (azaz az időnek egyetlen, adott körfrekvenciájú szinuszos függvényével írható le), normálkoordinátáknak, a nekik megfelelő ωi körfrekvenciákat pedig a rendszer saját(kör)frekvenciáinak nevezzük. Az olyan mozgásformákat, melyekben egyetlen normálkoordináta kivételével mindegyik zérus, sajátrezgésnek vagy normálrezgésnek nevezzük. Például a 3. példában szereplő gyöngyök az egyenletes, a lüktető és a tili-toli mozgásban sajátrezgéseket végeznek. Minden sajátrezgéshez tartozik egy amplitúdó-konfiguráció (a 3. példában a három sajátrezgéshez az amplitúdók 1:1:1, 1:0:1 és 1:(-2):1 aránya tartozott), ezeket normálmódusoknak nevezzük.
 
Vigh Máté


1Idő szerinti deriváltját.

2Az x1,x2,x3,... mennyiségek lineáris kombinációján a λ1x1+λ2x2+λ3x3+... alakú összeget értjük, ahol λ1,λ2,λ3,... (most) valós számok.