Vizualizace dynamických modelů ve Virtual RealityToolboxu v Matlabu
Vizualizace dynamických modelů ve Virtual RealityToolboxu v Matlabu
Vizualizace dynamických modelů ve Virtual RealityToolboxu v Matlabu
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
5.Wattův odstředivý regulátor 34<br />
W kr<br />
t= 1 2 J 2 t sin 2 t= 1 2 m l 2 2 tsin 2 t .<br />
Dále budeme uvažovat síly působící proti odstředivé síle. Kinetickou energii jednoho<br />
závaží při změně úhlu natočení ramene<br />
W ks<br />
t= 1 2 J ˙ 2 r1t = 1 2 m l 2 ˙ 2 r1<br />
t ,<br />
kde m je hmotnost závaží a ˙ [rad s −1 ] úhlová rychlost. Potenciální energie závaží<br />
W p<br />
t[ J ] je závislá na jeho výšce h=l cos a na gravitačním zrychlení g [m s −2 ]<br />
W p<br />
t= 1 2 m g h= 1 2 m g l cos .<br />
Lagrangeova funkce pro složený pohyb obou ramen<br />
Po dosazení energií do 5.12<br />
L t=2W ks<br />
t2W kr<br />
t2W p<br />
t . (5.12)<br />
L t=m l 2 ˙ 2 t m l 2 2 t sin 2 t 2 m l g cost .<br />
Přidáme tlumení způsobené třením závislém na rychlosti otáčení ramene ˙<br />
a vypočteme Lagrangerovu rovnici<br />
d<br />
dt ∂ L<br />
∂ ˙ − ∂ L =−2 δ J ˙t <br />
∂ <br />
2 m l 2 ¨t −4 m l 2 2 sin t cost2l m gsin t=2δ m l 2 ˙t<br />
Po úpravě dostaneme pohybovou rovnici úhlu natočení ramene v závislosti na otáčkách<br />
stroje <br />
¨t−2 2 sin 2t g l sint=δ ˙t . (5.13)<br />
Z rovnice nám vypadla hmotnost m . Ta by se projevila, pokud bychom uvažovali silové<br />
působení závaží na otáčky stroje. To ale lze zanedbat z důvodu nepoměru hmotnosti závaží<br />
regulátoru a hmotnosti stroje většinou opatřeného těžkým setrvačníkem.<br />
Nakonec je třeba převést úhel natočení ramene na <strong>ve</strong>likost průtoku páry nebo její tlak.<br />
Lze si vybrat mezi rotačními nebo posuvnými <strong>ve</strong>ntily, pákovými převody a kladkami.<br />
Vybereme posuvný <strong>ve</strong>ntil přímo připojený na nejspodnější část Wattova regulátoru. To lze<br />
vyjádřit jako x=2l cos s tím že se maximální otevření a zavření přepočítá z omezení<br />
natočení ramene . Pro nulové otáčky je úhel natočení ramene =0.2[rad ] a <strong>ve</strong>ntil je