DynaLab dynamika ramene

Tohle je statická kopie kapitoly pro vyhledávače. Interaktivní verze má animované obrázky, kontrolní otázky a tlačítka, která příklad načtou do kalkulačky.

Odstředivé a Coriolisovy členy

Prostřední člen rovnice je ten, kterému lidé nejmíň věří, protože v něm nevystupuje ani zrychlení, ani gravitace — a přesto tam moment je. Rameno, které se otáčí naprosto rovnoměrně, tedy s nulovým , potřebuje v ostatních kloubech nenulový moment. To není chyba, to je celá podstata pohybu v otáčející se soustavě.

Členy jsou dvojího druhu a rozeznají se podle indexů:

DruhTvarKdy vzniká
odstředivýj²otáčí-li se jeden kloub, tahá hmotu ven
Coriolisůvi·q̇jotáčejí-li se dva klouby zároveň

Coriolisův člen je zajímavější, protože zmizí, jakmile se hýbe jen jeden kloub. To je zrada pro zkoušení: rameno se otestuje po jedné ose, všechno vychází, a pak se pustí složený pohyb a moment je jinde, než čekal.

Nejmenší stroj, na kterém to jde vidět

tyč se točí, korálek na ní klouže ven a zpátkyjede-li ven, tyč se brzdíjede-li dovnitř, roztáčí semoment hybnosti m·r²·θ̇se sám od sebe neměníroste-li r, musí θ̇ klesnout
Korálek na roztočené tyči. Nic než hmotný bod, jedna rotace a jeden posuv — a obsahuje to oba členy naráz.

Kinetická energie korálku je T = ½m(ṙ² + r²θ̇²). Z Lagrangeovy rovnice vyjdou dvě rovnice a v každé je jeden z těch členů:

τθ = m·r²·θ̈  +  2m·r·ṙ·θ̇Coriolisův — brzdí otáčení, když korálek jede ven
Fr = m·r̈  −  m·r·θ̇²odstředivý — tahá korálek ven, ať se točí kterýmkoli směrem

Ten první má překvapivě jasný fyzikální smysl: moment hybnosti m·r²·θ̇ se sám od sebe nemění, takže když r roste, musí θ̇ klesat. Aby otáčení nekleslo, musí pohon dodat moment — a to je přesně ten Coriolisův člen. Je to piruetový efekt, jen zapsaný rovnicí.

m = 2,4 kgθ̇ = 2.0 rad/sr = 0.30 mṙ = 0.8 m/sCoriolis na otáčení: 2m·r·ṙ·θ̇2.30 N·modstředivá na posuv: −m·r·θ̇²-2.88 N
Coriolisův moment: 2,304 N·m — otáčení brzdí (korálek jede ven a odnáší si moment hybnosti)
Odstředivá síla: -2,880 N — vždy záporná, tedy vždy tahá ven, ať se točí kterýmkoli směrem
Zkus si to. Otoč znaménko rychlosti posuvu a Coriolisův moment se otočí taky — korálek jedoucí dovnitř roztáčí. Odstředivá síla se naproti tomu neotočí nikdy, protože v ní je θ̇ na druhou: točit se doleva nebo doprava je pro ni totéž.

Christoffelovy symboly a proč není C jednoznačná

Vektor c(q,q̇) se dá zapsat jako matice krát rychlost, c = C·q̇. Jenže takových matic je nekonečně mnoho — rychlost se dá vytknout mnoha způsoby. Jedna z nich je ale lepší než ostatní a dostane se z derivací matice setrvačnosti:

ckij = ½·(∂mki/∂qj + ∂mkj/∂qi − ∂mij/∂qk)Christoffelovy symboly prvního druhu

Rychlostní členy tedy nejsou nic nového — celé jsou schované v tom, jak se M(q) mění s polohou. Kdyby byla matice setrvačnosti konstantní, nebyly by žádné odstředivé ani Coriolisovy členy. To je asi nejužitečnější věta v téhle kapitole.

Šikmá symetrie, a proč na ní záleží

Právě s touhle volbou C platí, že matice Ṁ − 2C je šikmo symetrická. Formálně to znamená xᵀ(Ṁ − 2C)x = 0 pro každé x; fyzikálně to znamená, že rychlostní členy energii přeskupují, ale nikdy ji nevytvářejí ani neubírají. Přesně proto se odstředivé síly nedají použít jako zdroj energie, i když vypadají, že něco tahají.

Kde se to používá. Šikmá symetrie je základní stavební kámen důkazů stability regulátorů ramen — Lyapunovova funkce se postaví z celkové energie a ten člen se v ní vyruší. Je to zároveň nejcitlivější test dynamického řešiče, jaký existuje: prochází skrz derivace M, skrz Christoffelovy symboly i skrz rekurzi naráz, a nedá se ohnout jen v jedné z těch cest, aby vyšel.

Vzorce v této kapitole

c_kij — Christoffelovy symboly
c_kij = ½(∂m_ki/∂q_j + ∂m_kj/∂q_i − ∂m_ij/∂q_k) [kg·m²] De Luca, EL, slide 21
S — šikmá symetrie
xᵀ(Ṁ − 2C)x = 0 pro všechna x [—] De Luca, EL, slide 28
τ_cor — Coriolisův moment korálku
τ = 2m·r·ṙ·θ̇ [N·m] Lagrange
F_cf — odstředivá síla korálku
F = m·r·θ̇² [N] Lagrange