Tlumené kmitání

ikona
Tento článek není dostatečně ozdrojován, a může tedy obsahovat informace, které je třeba ověřit.
Jste-li s popisovaným předmětem seznámeni, pomozte doložit uvedená tvrzení doplněním referencí na věrohodné zdroje.
Tlumené kmity v případě tělesa na pružině

Tlumené kmitání je mechanické kmitání, které po určité době ustává. V reálném světě vždy existuje tření a různé jiné odporové síly, které způsobují, že oscilující systém postupně ztrácí energii a jeho amplituda se s časem zmenšuje. Jakkoliv bychom se snažili zamezovat těmto nepříznivým vlivům, omezuje nás druhý zákon termodynamiky, podle kterého se mechanická energie postupně přeměňuje na vnitřní tepelnou energii. Závislost odporových sil od výchylky a její časové derivace může být obecně velmi složitá. V nejjednodušším modelu je odporová síla přímo úměrná první časové derivaci výchylky a působí proti narůstání výchylky. Takovému případu tlumení říkáme lineární tlumení a popisujeme ho vztahem F t = b x ˙ {\displaystyle F_{t}=-b\cdot {\dot {x}}} . V případě tělesa na obrázku vpravo odporová síla vzduchu působí vždy proti směru pohybu tohoto tělesa.

Ve zmíněném ideálním případě lze zapsat pro časový vývoj výchylky následující diferenciální rovnici popisující silovou rovnováhu:

M x ¨ + b x ˙ + k x = 0 {\displaystyle M{\ddot {x}}+b{\dot {x}}+kx=0}

 

 

 

 

kde M {\displaystyle M} , k {\displaystyle k} a b {\displaystyle b} je v případě mechanického oscilátoru hmotnost, tuhost pružiny a koeficient lineárního tlumení.. Pokud obě strany rovnice vydělíme M {\displaystyle M} , dostaneme

x ¨ + b M x ˙ + k M x = 0 {\displaystyle {\ddot {x}}+{b \over M}{\dot {x}}+{k \over M}x=0}

 

 

 

 

Zavedeme substituce, které nám později ulehčí práci při popisu řešení:

ω 0 = k M {\displaystyle \omega _{0}={\sqrt {k \over M}}}

 

 

 

 

ζ = b 2 k M {\displaystyle \zeta ={b \over 2{\sqrt {k\,M}}}}

 

 

 

 

První parametr ω 0 {\displaystyle \omega _{0}} se nazývá vlastní úhlová rychlost a určuje frekvenci ( ω = 2 π f {\displaystyle \omega =2\pi f} ) v případě, kdyby na soustavu nepůsobily žádné tlumící vlivy. Druhý parametr ζ {\displaystyle \zeta } nazýváme relativní tlumení a podle jeho hodnoty pak určujeme typy tlumení.

Diferenciální rovnice nyní nabyla tvar:

x ¨ + 2 ζ ω 0 x ˙ + ω 0 2 x = 0 {\displaystyle {\ddot {x}}+2\zeta \omega _{0}{\dot {x}}+\omega _{0}^{2}x=0}

 

 

 

 


Tato homogenní diferenciální rovnice druhého stupně je řešitelná, předpokládáme, že hledané řešení x ( t ) {\displaystyle x(t)} má tvar

x ( t ) = C 1 e λ 1 t + C 2 e λ 2 t {\displaystyle x(t)=C_{1}e^{\lambda _{1}\,t}+C_{2}e^{\lambda _{2}\,t}}

 

 

 

 


kde λ jsou kořeny charakteristické rovnice.

λ 2 + 2 ζ ω 0 λ + ω 0 2 = 0 {\displaystyle \lambda ^{2}+2\zeta \omega _{0}\lambda +\omega _{0}^{2}=0}

 

 

 

 

λ 1 , 2 = ω 0 ζ ± ω 0 ζ 2 1 {\displaystyle \lambda _{1,2}=-\omega _{0}\zeta \pm \omega _{0}{\sqrt {\zeta ^{2}-1}}}

 

 

 

 

Chování systému při tlumených kmitech

Časový průběh výchylky oscilátoru při různých typech tlumení

Při řešení charakteristické rovnice mohou nastat tři případy v závislosti na parametru relativního tlumení, podle kterých rozlišujeme tři situace

Nadkritické tlumení (ζ > 1)

λ 1 λ 2 ;   I m ( λ 1 , 2 ) = 0 {\displaystyle \lambda _{1}\neq \lambda _{2};\ Im(\lambda _{1,2})=0} (dva různé reálné kořeny)

Při nadkritickém tlumení je systém přetlumený a při vychýlení se pomalu se vrací do rovnovážné polohy. Řešením diferenciální rovnice, která popisuje takový systém je funkce:


x ( t ) = C 1 e λ 1 t + C 2 e λ 2 t {\displaystyle x(t)=C_{1}e^{\lambda _{1}t}+C_{2}e^{\lambda _{2}t}}


Dosadíme-li počáteční podmínky pro vychýlení x 0 = C 1 + C 2 {\displaystyle x_{0}=C_{1}+C_{2}} a počáteční rychlost v 0 = C 1 λ 1 + C 2 λ 2 {\displaystyle v_{0}=C_{1}\lambda _{1}+C_{2}\lambda _{2}} , můžeme potom odvodit vztah pro výchylku nakriticky tlumeného oscilátoru v závislosti na čase a počátečních podmínkách:


x ( t ) = 1 λ 2 λ 1 [ ( λ 2 x 0 v 0 ) e λ 1 t + ( v 0 λ 1 x 0 ) e λ 2 t ] {\displaystyle x(t)={\frac {1}{\lambda _{2}-\lambda _{1}}}\left[(\lambda _{2}x_{0}-v_{0})\cdot e^{\lambda _{1}t}+(v_{0}-\lambda _{1}x_{0})\cdot e^{\lambda _{2}t}\right]} [1]

Kritické tlumení (ζ = 1)

λ 1 = λ 2 ;   I m ( λ 1 , 2 ) = 0 {\displaystyle \lambda _{1}=\lambda _{2};\ Im(\lambda _{1,2})=0} (jeden dvojnásobný reálný kořen)

Při kritickém tlumení se oscilátor ustálí v rovnovážné poloze nejrychleji. Získáme řešení charakteristické rovnice λ 1 , 2 = ω 0 {\displaystyle \lambda _{1,2}=-\omega _{0}} a řešením diferenciální rovnice je pak funkce:


x ( t ) = ( C 1 + C 2 t ) e ω 0 t {\displaystyle x(t)=(C_{1}+C_{2}\,t)\cdot e^{-\omega _{0}\,t}}


Dosazením počátečních podmínek x 0 = C 1 {\displaystyle x_{0}=C_{1}} a v 0 = ω 0 C 1 + C 2 {\textstyle v_{0}=-\omega _{0}C_{1}+C_{2}} a následnou úpravou pak dostaneme vztah pro závislost výchylky na počátečních podmínkách a čase:


x ( t ) = [ x 0 + ( v 0 + ω 0 x 0 ) t ] e ω 0 t {\displaystyle x(t)=\left[x_{0}+(v_{0}+\omega _{0}x_{0})\,t\right]e^{-\omega _{0}\,t}} [1]

Podkritické tlumení (ζ < 1)

I m ( λ 1 , 2 ) 0 {\displaystyle Im(\lambda _{1,2})\neq 0} (komplexně sdružené kořeny)

Při podkritickém tlumení zůstává v průběhu oscilátoru kmitavá složka, jejíž amplituda exponenciálně klesá.

Součinitel tlumení způsobí, že je pod odmocninou v řešení charakteristické rovnice nula a jejími kořeny jsou komplexní čísla. Ty pak přepíšeme do tvaru λ 1 , 2 = b r ω 0 ± i Ω {\displaystyle \lambda _{1,2}=-b_{r}\omega _{0}\pm i\Omega } , kde Ω = 1 b r 2 {\displaystyle \Omega ={\sqrt {1-b_{r}^{2}}}} je vlastní úhlová frekvence tlumených kmitů. Řešení diferenciální rovnice x ( t ) = C 1 e λ 1 t + C 2 e λ 2 t {\displaystyle x(t)=C_{1}e^{\lambda _{1}t}+C_{2}e^{\lambda _{2}t}} se upraví podle Eulerova vztahu a dostaneme tak funkci:


x ( t ) = e b r ω 0 t [ A   c o s ( Ω t ) + B   s i n ( Ω t ) ] {\displaystyle x(t)=e^{b_{r}\omega _{0}t}\left[A\ cos(\Omega \,t)+B\ sin(\Omega \,t)\right]}


Dosazením počátečních podmínek x 0 = A {\displaystyle x_{0}=A} a v 0 = A b r ω 0 + B Ω {\displaystyle v_{0}=-Ab_{r}\omega _{0}+B\Omega } a úpravou získáme výsledný vztah pro výchylku v závislosti na čase a počátečních podmínkách:


x ( t ) = e b r ω 0 t [ x 0   c o s ( Ω t ) + x 0 b r ω 0 + v 0 Ω   s i n ( Ω t ) ] {\displaystyle x(t)=e^{-b_{r}\omega _{0}t}\left[x_{0}\ cos(\Omega \,t)+{\frac {x_{0}b_{r}\omega _{0}+v_{0}}{\Omega }}\ sin(\Omega \,t)\right]} [1]

Reference

  1. a b c VALÁŠEK, Michael, Václav BAUMA a Zbyněk ŠIKA. Mechanika B. Vyd. 1. Praha: ČVUT, 2004. 121 s. .

Externí odkazy

  • Logo Wikimedia Commons Obrázky, zvuky či videa k tématu Tlumené kmitání na Wikimedia Commons