Cím: Anharmonikus rezgések periódusideje
Szerző(k):  Woynarovich Ferenc 
Füzet: 2020/szeptember, 365 - 373. oldal  PDF  |  MathML 
Témakör(ök): Szakmai cikkek, Egyéb nyújtás, összenyomás, Newton-törvények

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.

 
Anharmonikus rezgések periódusideje
 

Ha egy egyensúlyi helyzetéből kimozdított testet visszahúzó erő a Hooke-törvényt követi, azaz a kitéréssel arányos, a test harmonikus rezgőmozgást végez. Ennek egyik fontos jellegzetessége, hogy a rezgés ideje nem függ a kitérítés nagyságától. A valóságban azonban az erőtörvény mindig csak bizonyos határok között, valamilyen közelítésben tekinthető lineárisnak. Bizonyos esetekben ezek a határok lehetnek egészen tágak (pl. a direkt erre a célra készített rugalmas eszközök esetén), máskor el kell fogadnunk, hogy a vizsgált periodikus mozgás csak nagyon kicsiny amplitúdó esetén tekinthető harmonikusnak (ilyen pl. az inga mozgása, és általában az összetett rendszerek egyensúly körüli rezgései). Úgy is mondhatjuk, hogy ezekben az esetekben a harmonikus közelítés csak az első közelítés, ami tovább finomítható. Az anharmonikus rezgések periódusideje már nem független az amplitúdótól, de az első korrekció meglepően könnyen kiszámítható. A továbbiakban ezt fogjuk bemutatni egy egyszerű, ámde könnyen általánosítható példán.
Először elevenítsünk fel néhány dolgot, amit a harmonikus rezgésekről tudunk! Ha egy test mozgását az 
ma=-Dx(1)
mozgásegyenlet írja le (amelyben m a test tömege, a a gyorsulása, x egy adott ponthoz viszonyított kitérése és D egy pozitív állandó), akkor ez a test harmonikus rezgőmozgást végez. Erre jellemző, hogy a kitérést, a mozgás sebességét és a gyorsulást rendre az
x=Asin(ω0t+φ0),(2)v=Aω0cos(ω0t+φ0),(3)
és az
a=-Aω02sin(ω0t+φ0)(4)
függvények adják meg. Ezekben A a rezgés amplitúdója, az ω0 körfrekvencia ω0=D/m, a φ0 szöget pedig a test t=0 időpontban elfoglalt helyzete határozza meg. A rezgés periódusideje T0=2π/ω0. A későbbiekben fontos összefüggés lesz az energiamegmaradás törvénye, ezért érdemes felidéznünk, hogy ha egy testet F(x)=-Dx erő húz vissza, akkor ezen erővel szemben V(x)=Dx2/2 munkát kell végeznünk, hogy a testet az origóból az x pozícióba vigyük. (Ezt úgy szoktuk mondani, hogy az F(x)-nek a V(x) a potenciálja. Azzal, hogy ez mit is jelent pontosan, a Függelékben még foglalkozunk.) Ennek segítségével az energia mérlegegyenlete:
12mv2+V(x)=12DA2.(5)

Következő lépésként egy anharmonikus esetet próbálunk leírni. Bár nem minden egyensúlyi helyzet szimmetrikus (vagyis F(-x)=-F(x) tulajdonságú erőtörvénynek megfelelő), mi most csak ilyet vizsgálunk, ezek közül is a legegyszerűbbet vesszük. Ebben a kitérített testet visszahúzó erőt egy harmadfokú kifejezés adja meg, tehát a mozgásegyenlet
ma=F*(x)=-Dx-D*x3.(6)
Itt D* előjele nincs megkötve, és feltételezzük, hogy az anharmonikus D*x3 tag az egész mozgás során kicsiny az előtte állóhoz viszonyítva, azaz |D*|DA21. A továbbiakban azt vizsgáljuk meg, milyen mértékben változtatja meg ez a tag a rezgésidőt. Először két egyszerű, a harmonikushoz hasonló közelítéssel próbálkozunk.
1) A harmonikus rezgés esetén amax=Aω2. Ennek analógiájára a (6) egyenletből az
ω1=ω01+bA2(7)
érték adódik, ahol b=D*/D.
2) A harmonikus rezgés során vmax=Aω. Kiindulhatunk ebből az összefüggésből is! Ahogy azt a Függelékben bemutatjuk, a (6) egyenletben szereplő F*(x) erőhöz tartozó potenciál
V*(x)=12Dx2+14D*x4,
tehát az energiamegmaradás egyenlete most
12mv2+12Dx2+14D*x4=12DA2+14D*A4.(8)
A vmax sebességet ebből számolva az
ω2=ω01+bA22(9)
kifejezés adódik. Tanulságos megnéznünk azokat a harmonikus potenciálokat, amelyekben az m tömegű próbatestünk ω1, illetve ω2 körfrekvenciával rezegne. Ezek rendre
V1(x)=12D(1+bA2)x2-14D*A4
és
V2(x)=12D(1+bA22)x2.
A potenciálokhoz tetszőleges konstansokat hozzá lehet adni, ezek az energiamérlegből úgyis kiesnek. Mi most úgy választottuk meg a potenciális energiák nullpontját, hogy
V1(±A)=V*(±A)=V2(±A)(10)
legyen. A három potenciál jellegét D*>0 esetén az ábra mutatja.
Fontos észrevétel, hogy a szélső (±A) pontokban V1 érintője is megegyezik V*-éval, hiszen itt a potenciálfüggvény meredeksége az erővel, az pedig a test maximális gyorsulásával arányos, és V1 paraméterét éppen úgy választottuk meg, hogy a test legnagyobb gyorsulása ugyanannyi legyen a két esetben. Ez nem igaz V2-re, de V2 és V* érintik egymást az origóban, ami a maximális sebességek összehangolásának következménye.
 
 

Jól látszik, hogy a -AxA szakaszon a V1 és V2 közrefogja a V*-ot, pl. D*>0 mellett
V1(x)V*(x)V2(x),ha-AxA,
és a viszony pont fordított, ha D*<0. Ugyanakkor (10) miatt
12mv12+V1(x)=12mv2+V*(x)=12mv22+V2(x),
tehát a pálya bármely adott pontjában D*>0 mellett
|v1||v||v2|,
aminek egyenes következménye, hogy a tényleges rezgésidőre fennáll, hogy
T1=2πω1<T<T2=2πω2,haD*>0,
és a sorrend pont fordított az ellenkező esetben:
T1>T>T2,haD*<0.

Joggal vetődik fel, hogy most melyik érték áll közelebb a valódihoz. Az a meglepő válasz, hogy a kettő számtani közepe, vagyis
TT1+T22(11)
mindkettőnél sokkal jobb közelítés. A továbbiakban ezt fogjuk bemutatni, és becslést adunk (11) pontosságára is.
A rezgésidő kiszámítására kézenfekvőnek tűnik, hogy a (-A,+A) szakaszt kicsiny Δxi szakaszokra osztjuk, és felösszegezzük azokat a Δti időket, amelyek a kis szakaszok megtételéhez kellenek:
T=2iΔti=2iΔxivi.
(Az ehhez szükséges sebességeket az energiamegmaradásból számíthatjuk ki.) Ezzel azonban vigyáznunk kell, mert a szélső helyzetekhez közeledve a test nagyon lelassul, nagyon kicsiny sebességértékek kerülnek a nevezőbe, és ügyeskednünk kell, hogy az összeg értékét kellő pontossággal tudjuk meghatározni. Szerencsére ezt a problémát meg tudjuk kerülni, ha figyelembe vesszük, hogy a rezgés mégiscsak egy periodikus mozgás, azaz a kitérés megadható az
x(t)=Asinφ(t)
alakban. Itt φ(t) az idő monoton növekedő függvénye, és a mozgás periodikusságát az fejezi ki, hogy
φ(t+T)=φ(t)+2π,
de olyan szimmetrikus rezgés esetén, mint amit most vizsgálunk, az ennél szigorúbb
φ(t+T/2)=φ(t)+π
feltétel is teljesül. A φ(t) növekedési ütemét egy pillanatnyi szögsebességgel lehet megadni. Ennek definíciója a pillanatnyi sebesség analógiájára
ω(t)=ΔφΔt=φ(t+Δt)-φ(t)Δt,(Δt0).
(Ezt úgy kell értenünk, hogy miközben Δt-nek egyre kisebb értékeket választunk, aközben a számláló is egyre kisebb lesz, a kettő hányadosa viszont egy konkrét értékhez közeledik. Ahogy az út ‐ idő grafikon érintőjének a meredeksége a pillanatnyi sebesség, úgy a φ(t) grafikon érintőjének a meredeksége az ω(t) pillanatnyi szögsebesség.) A rezgés pillanatnyi sebességét annak definíciójából kiindulva határozzuk meg:
v(t)=ΔxΔt=Asin(φ+Δφ)-AsinφΔt=AsinΔφcosφ-(1-cosΔφ)sinφΔt.
Mivel kicsi Δφ esetén sinΔφΔφ, továbbá
1-cos(Δφ)=2sin2(Δφ/2)(Δφ)2/2,
a sebességet megadó képletben a számláló második tagja jóval kisebb, mint az első, ezért elhagyható, a maradék pedig a
v(t)=Aω(t)cosφ(t)(12)
kifejezést adja. Ez nagyon hasonlít a harmonikus rezgés sebességére, de vigyázzunk, az analógia nem teljes, mert a gyorsulás
a(t)-Aω2(t)sinφ(t).

Ahogy a v sebesség, úgy az ω pillanatnyi szögsebesség is kifejezhető az energiamegmaradás (8) egyenlete segítségével. Az ebből adódó
v=ω0(A2-x2)(1+b(A2+x2)2)
kifejezésből egyszerű behelyettesítéssel kapjuk, hogy
ω(t)=ω01+bA2(1+sin2φ(t))2.
Vegyük észre, hogy a szögsebesség az időtől csak a φ-n keresztül függ, tehát nyugodtan írhatjuk, hogy
ω(φ)=ω01+bA2(1+sin2φ)2.

Most már fel tudjuk írni a rezgésidőt olyan összeg formájában, amely esetében nem kell attól tartanunk, hogy a nevezőbe eltűnő mennyiség kerül. Ha a 2π tartományt kicsiny Δφi szakaszokra osztjuk, és φi-nek a megfelelő szakasz egy pontját (mondjuk a felezőpontját) választjuk, akkor az adott szakasz ,,megtételéhez'' szükséges idő ΔtiΔφi/ω(φi), azaz
T=iΔφiω(φi)=iΔφiω0(1+bA2(1+sin2φi)2)-1/2.
Az összegzés (integrál) kellő matematikai felkészültséggel egzakt módon is elvégezhető, de mi most sokkal kevesebbel is beérjük: megelégszünk azzal, hogy kiszámítsuk a T-hez a bA2 paraméterben legalacsonyabb rendű korrekciót. Ezt a megfelelő (1+y)1+y/2 és (1+y)-11-y közelítésekkel átalakított
TiΔφiω0(1-bA2(1+sin2φi)4)
formula adja meg, ami a sin2φ=1-cos2φ2 helyettesítés után
TiΔφiω0(1-38bA2+18bA2cos2φi).
Látható, hogy a zárójelben a harmadik tag az összegzésben nullát ad (mert ugyanannyi a pozitív és a negatív járulék), a többi pedig egyszerűen szorzódik az intervallum 2π hosszával, így
T2πω0(1-38bA2)=T0(1-38bA2).(13)
Összehasonlításképpen T1 és T2 ugyanilyen pontossággal (ugyanazon 1+y
1+y/2 és (1+y)-11-y közelítéseket alkalmazva):
T1=T01+bA2T0(1-12bA2),(14)
illetve
T2=T01+bA2/2T0(1-14bA2),(15)
tehát valóban fennáll, hogy nagyon jó közelítéssel
TT1+T22.(16)

Szólnunk kell még az eredményünk pontosságáról. Számolásainkban az egyhez képest kis paraméter a bA2 mennyiség. Ahogy azt a Függelékben bemutatjuk, az alkalmazott formulák az ott y-nal jelölt kis paraméterben első rendig pontosak, az elhanyagolt tagok közül a legnagyobb y2 nagyságrendű. Esetünkben ez azt jelenti, hogy az eredményünk hibája kb. (bA2)2. Ezt úgy szoktuk jelölni, hogy
T=T1+T22+O(A4),(17)
ahol az O(A4) (olvasd: ordó A a negyediken) olyan ‐ pontosabban, számszerűen nem részletezett ‐ tagokat jelent, amelyek legfeljebb A4 nagyságrendűek. (Az ordo latin szó, jelentése: rend.) Ezekben természetesen az A különböző hatványai mellett olyan (A-tól független) szorzótényezők állnak, amelyekkel együtt minden tag ugyanolyan dimenziójú, de ezek az együtthatók minden konkrét feladatban adottak, tehát a nagyságrendet maga az A hatványkitevője határozza meg. Összefüggésünk ebben a formában általánosnak tekinthető minden szimmetrikus anharmonikus rezgésre. Gondoljuk csak meg, ha az erőhöz hozzáadunk egy x5-nel (azaz a potenciálhoz egy x6-nal) arányos tagot, attól még a számolás ugyanígy elvégezhető, és arra az eredményre vezet, hogy a T, T1 és T2 csak A4 nagyságrendű tagokkal módosul. Nem nő a korrekciók nagyságrendje akkor sem, ha az erőhöz (potenciálhoz) további, még gyorsabban eltűnő tagokat adunk, tehát ezek a (17) összefüggés tartalmát már nem változtatják meg.
 
Megjegyzés. A leírt eljárás pontosságának érzékeltetésére tekintsünk egy konkrét példát! Legyen mondjuk D=0,1N/cm és D*=10-4N/cm3, vagyis b=10-3cm-2. Ha a rezgés amplitúdója A=10cm, akkor bA2=0,1. Ez annyit jelent, hogy a legnagyobb kitérésnél az erő képletében az anharmonikus tag 10-szer kisebb, mint a Hooke-törvénynek megfelelő harmonikus tag. A rezgő test m tömege meghatározza a nagyon kicsi kitérésekhez tartozó T0 rezgésidőt, ennek számszerű értékére a relatív eltérések számításánál nincs szükségünk.
A (14), (15) és (16) összefüggéseknek megfelelően a harmonikus mozgással közelített rezgések periódusideje:
T1=11+0,1T0=0,9535T0,T2=11+0,05T0=0,9759T0,
és ezek számtani közepe
Tátlag=0,9647T0.
Másrészt az ω(φ)=2πT01+0,05(1+sin2φ) függvény reciprokának integrálásával (pl. a www.wolframalpha.com segítségével) megkaphatjuk az anharmonikus rezgések ,,pontos'' periódusidejét:
T=02π1ω(φ)dφ=T02π02π11+0,05(1+sin2φ)dφ=0,9645T0.
Ezekből a számokból leolvasható, hogy T1 és a pontos érték relatív eltérése -1,14%, ugyanez T2-nél +1,18%, míg az átlagolt közelítő rezgésidőre a relatív eltérés mindössze 210-4.

 
Függelék

1. Erők és potenciálok egy dimenzióban
 

Ha egy esetben az erő kizárólag a helytől függ (azaz nem függ pl. a mozgásállapottól), tehát F=F(x), akkor ebben a rendszerben definiálható a V(x) helyzeti energia, más néven potenciális energia, vagy egyszerűen potenciál. Ez a fizikai mennyiség megegyezik azzal a munkával, amit az F(x) ellenében kell végeznünk, ha a testet egy előre kiválasztott x0 pontból az x pontba visszük. A kiszámításához az (x,x0) szakaszt xi osztópontokkal olyan kicsiny szakaszokra vágjuk, hogy azokban az erő egy-egy szakaszon belül már állandónak vehető legyen, és a potenciált a
V(x)=i-F(x'i)Δxi(F. 1)
összeg adja meg. Itt az xi'Δxi=(xi+1-xi) előjeles (irányított) szakasz egy pontja. Nyilván, ha az a koordinátájú pontból visszük a testet a b koordinátájú pontba, akkor W(b,a)=V(b)-V(a) munkát kell végeznünk. Jól látszik, hogy a V(x)-hez hozzáadhatunk egy tetszőleges x-független értéket, mert az az energiaváltozásból úgyis kiesik. Azt is könnyen beláthatjuk, hogy ha az eredő erőt több erő összege adja meg, az eredő potenciál az egyes erőkhöz tartozó potenciálok összege lesz. A potenciál konstrukciójából következik, hogy
F(x)=-V(x+Δx)-V(x)Δx,(F. 2)
ha Δx nagyon kicsi. (A differenciálszámításban jártasak felismerik: F(x)-V(x)  függvény deriváltja.)
Ennek alapján pl. könnyen belátható, hogy ha az erő F=-kx3 alakú, akkor a hozzá tartozó potenciál V=(kx4)/4 szerint függ az x-től:
k4(x+Δx)4-x4Δx=k4(4x3+6x2Δx+4x(Δx)2+(Δx)3).
(Ezt a kifejezést a binomiális tétel alkalmazásával kaptuk.) Ha Δx-et egyre kisebbnek (,,végtelenül kicsinek'') választjuk, a zárójelben csak az első tag marad véges, tehát az erő valóban -kx3.
 
Megjegyzés. A fenti megfontolás (miszerint az erőfüggvény mindig egy potenciálfüggvényből származtatható deriválással) csak egydimenziós erőtereknél érvényes. Két- vagy háromdimenziós erőtereknél (ha azok ,,nem konzervatívak'', más náven: örvényesek) előfordulhat, hogy potenciálfüggvény nem létezik.

 
2. Két közelítő formula és ezek pontossága

2/1. A geometriai sor összegképlete szerint |y|<1 esetén
11+y=n=0(-1)nyn.
Ha |y|1, az egymást követő tagok lényegesen kisebbek, mint a megelőző tag. Akkor, ha valahol ,,levágjuk'' az összeget, az ezzel elkövetett hiba nagyságrendje azonos az első elhanyagolt tag nagyságrendjével. Így az
11+y1-y(F. 3)
közelítés hibája y2 nagyságrendű. Fontos megjegyeznünk, hogy itt a nagyságrend nem egyszerűen azt jelenti, hogy ,,nagyjából akkora (jelen esetben: olyan kicsi), mint'', hanem ‐ függvényről lévén szó ‐ azt is, hogy ,,körülbelül úgy viselkedik, mint''. Esetünkben a hiba pontosan y2/(1+y), de mivel kicsi y esetén az 1/(1+y) lassan változik, a hiba viselkedésének a jellegét is az y hatványkitevője határozza meg.
 

2/2. Egy másik közelítő képlet, amit gyakran használunk, a 1+y-re vonatkozik. Itt a következő átalakítással élünk:
1+y=(1+y2)2-(y2)2=(1+y2)1-(y/2)2(1+y/2)2.(F. 4)
Látni való, hogy gyök alatt egy ugyanolyan típusú kifejezés (az egy mellett egy kicsi tag) szerepel, mint a kiinduló kifejezésben, ráadásul ez közelebb van az egyhez, mint az eredeti. Így az eljárás akárhányszor megismételhető, ezzel egyre több tényező hozható ki a gyök alól úgy, hogy közben az ott maradó kifejezés egyre jobban megközelíti az egyet. Mi most nem ezt az utat választjuk, hanem megelégszünk azzal, hogy megbecsüljük a
1+y1+y2(F. 5)
első közelítés hibáját. Ez
(1+y2)-1+y=(y/2)2(1+y/2)+1+y.
Ha most |y|1, a jobb oldal nevezője egy 2-höz közeli szám, vagyis (F. 5) hibája is y2 nagyságrendű.
 
Woynarovich Ferenc  
Budapest