Tlumené kmitání

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\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: Šablona:Vzorec

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 Šablona:Vzorec

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

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:

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

kde λ jsou kořeny charakteristické rovnice.

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 nadkriticky 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}[({\lambda }_2{x}_0-v_0)\cdot e^{{\lambda }_{1}t}+(v_0-{\lambda }_1{x}_0)\cdot e^{{\lambda }_{2}t}]$

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)=[x_0+(v_0+{\omega }_0{x}_0)t]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\mathrm{\Omega }$, kde $\mathrm{\Omega }={\omega }_0\sqrt{1-b_r^2}$ je vlastní úhlová frekvence tlumených kmitů. Ř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}[A\cos (\mathrm{\Omega }t)+B\sin (\mathrm{\Omega }t)]$

Dosazením počátečních podmínek $x_0=A$ a $v_0=-A{b}_r{\omega }_0+B\mathrm{\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}[x_0\cos \left(\mathrm{\Omega }t\right)+\frac{x_0{b}_r{\omega }_0+v_0}{\mathrm{\Omega }}\sin \left(\mathrm{\Omega }t\right)]$

Externí odkazy

• Šablona:Commonscat