Tlumené kmitání
Author
Albert FloresTlumené 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 \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: {{vzorec|M \ddot{x} + b \dot{x} + k x = 0}}
kde M , k a 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, dostaneme {{vzorec|\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í: {{vzorec|\omega_0 = \sqrt{ k \over M }}} {{vzorec|\zeta = { b \over 2 \sqrt{k \, M} }}}
První parametr \omega_0
se nazývá vlastní úhlová rychlost a určuje frekvenci (\omega = 2\pi f ) v případě, kdyby na soustavu nepůsobily žádné tlumící vlivy. Druhý parametr \zeta nazýváme relativní tlumení a podle jeho hodnoty pak určujeme typy tlumení.
Diferenciální rovnice nyní nabyla tvar:
{{vzorec|\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) má tvar
{{vzorec| x(t) = C_1e^{\lambda_1 \, t} + C_2 e^{\lambda_2 \, t} }}
kde λ jsou kořeny charakteristické rovnice.
{{vzorec|\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)
\lambda_1 \ne \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^{\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 a počáteční rychlost 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) = \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]
=== Kritické tlumení (ζ = 1) === \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 \lambda_{1,2} = -\omega_0 a řešením diferenciální rovnice je pak funkce:
x(t) = (C_1 + C_2 \, t) \cdot e^{ -\omega_0 \, t}
Dosazením počátečních podmínek x_0 = C_1 a 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) = \left[x_0 + (v_0 + \omega_0 x_0) \, t \right] e^{-\omega_0 \, t}
Podkritické tlumení (ζ < 1)
Im(\lambda_{1,2}) \ne 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 \lambda_{1,2} = -b_r \omega_0 \pm i \Omega , kde \Omega = \sqrt{1 - b_r^2} je vlastní úhlová frekvence tlumených kmitů. +more Řešení diferenciální rovnice 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 \omega_0 t } \left[ A \ cos (\Omega \, t )+ B \ sin (\Omega \, t) \right]
Dosazením počátečních podmínek x_0 = A a v_0 = -A b_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 \omega_0 t } \left[ x_0 \ cos (\Omega \, t ) + \frac{x_0 b_r \omega_0 + v_0}{\Omega} \ sin (\Omega \, t) \right]