Matematické Fórum

Nevíte-li si rady s jakýmkoliv matematickým problémem, toto místo je pro vás jako dělané.

Nástěnka
22. 8. 2021 (L) Přecházíme zpět na doménu forum.matweb.cz!
04.11.2016 (Jel.) Čtete, prosím, před vložení dotazu, děkuji!
23.10.2013 (Jel.) Zkuste před zadáním dotazu použít některý z online-nástrojů, konzultovat použití můžete v sekci CAS.

Nejste přihlášen(a). Přihlásit

#1 18. 04. 2014 00:43 — Editoval jelena (01. 05. 2014 19:32)

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Odvození dif. rovnice z modelu

Ahoj,


poradil by někdo jak z model:

//forum.matweb.cz/upload3/img/2014-04/74557_model.jpg

popsat diferenciální rovnicí?

popis modelu v příspěvku 3.

Offline

 

#2 19. 04. 2014 11:56

jelena
Jelena
Místo: Opava
Příspěvky: 30020
Škola: MITHT (abs. 1986)
Pozice: plním požadavky ostatních
Reputace:   100 
 

Re: Odvození dif. rovnice z modelu

Zdravím,

to dost skromně popsaný model. Předpokládám, že je tvořen z pružin (zadána tuhost pro každou pružinu) a systémem závaží (systém se pohybuje v uzavřeném pístu?). Tak by prospělo všechno ještě doplnit a podrobně popsat. Tučné černé čáry jsou nehmotné?

Začala bych náhradou paralelních pružin za 1 pružinu s tuhosti $k_{34}$. Potom budeš mít systém pružin "za sebou". A pak se uvidí. Rozhodně tématu prospěje podrobnější úvod a co k tomu všemu umíš a máš. Děkuji.

Offline

 

#3 01. 05. 2014 10:54 — Editoval karlB (01. 05. 2014 10:55)

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

Zdravím Vás,

dolňuji některé informace:

Ano je dná se o pružiny.Zadána tuhost (k), tření/tumení  (b), na vrchol konstrukce bude později položeno závaží.
X je směr působení síly (vyvolané položením závaží). Půjde tedy o tlumené kmity, které je možné popsast rovnicí:

$m\frac{\mathrm{d} x}{\mathrm{d} t}^{2} + b\frac{\mathrm{d} x}{\mathrm{d} t}+kx = F(t)$

O čarách se nic nepíše, takže předpokládám, že v modelu nemají význam. Pouze vím, že se jedná se o konstrukci.

Podařilo se mi dohledat jak popstat soustavu dvou pružin v sérii a paralelně.

Sériově:

$F = -kx=-\frac{k_{1}k_{2}}{k_{1}+k_{2}}x$

Paralelně:

$F = -kx=-({k_{1}+k_{2}})x$

Dále že pro model musí platit z jistého fyz. zákona rovnost sil

$F=F_{g}+F_{b}+F_{p}$ 

z čehož vychází výše zmíněná dif. rovnice pro tlumené kmity.

Jen teď nevím jak to všechno dát do hromady..

Zatím nastíním moji úvahu:

1)
Kdybych měl jen závaží o hmotnosti m , které položím na "desku" pod kterou bude jedna pružina, potom přiložením závaží na desku půdobí gravitační síla směrem dolů a síly tlumící, pružiny, a akcelerační (m*a) směrem nahoru..

Rovnice pro tento zjednodušený případ by tedy moha mít tvar:

$m\frac{\mathrm{d} x}{\mathrm{d} t}^{2} + b\frac{\mathrm{d} x}{\mathrm{d} t}+kx = m*g$

2) Rovnice popisující celý model bude mít podobný tvar, akorát bude složitější. Např síla pružin by mohla být definována podle toho jestli jsou seriově nebo paralelně:

$F=-kx=-(\frac{k_{1}k_{2}}{k_{1}+k_{2}}+k_{3}+k_{4})x$

Celková třecí síla - tady už nevím jak by se mohla definovat..

Děkuji za jakoukoliv radu:-)

Offline

 

#4 03. 05. 2014 18:31

jelena
Jelena
Místo: Opava
Příspěvky: 30020
Škola: MITHT (abs. 1986)
Pozice: plním požadavky ostatních
Reputace:   100 
 

Re: Odvození dif. rovnice z modelu

Zdravím a děkuji za doplnění modelu, doufala jsem, že někoho z kolegů již doplněný model osloví a téma se pohne. Tak v to pořád doufám :-)

Ohledně spojení pružin do série a paralelně, to bych měla stejně. Ohledně tření/tlumení - řekla bych, že zadání $b$ je pro celý systém a bude zahrnovat jak tření v horním válci, tak i odpor vzduchu, případně další vlivy směrem k tlumení. Jelikož osobně nejsem si jistá s náhradou části "série", když v tom jsou ještě hmotnosti, tak poprosím kolegy o další náměty. Děkuji.

Offline

 

#5 29. 05. 2014 19:04

FliegenderZirkus
Příspěvky: 544
Škola: RWTH Aachen
Reputace:   25 
 

Re: Odvození dif. rovnice z modelu

↑ karlB:

Ozvi se prosím, pokud je téma ještě aktuální.

Offline

 

#6 31. 05. 2014 00:27

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

↑ FliegenderZirkus:

aktuální to je, něco jsem už vymyslel, ale nevím jestli je to 100 pro dobře. Můžu ti to někam poslat, nechce se mi to celý přepisovat..

Offline

 

#7 31. 05. 2014 00:56

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

ale vlastně by to mohlo jít takhle:

//forum.matweb.cz/upload3/img/2014-05/90544_img.jpg

a = klidová délka pružin..

Nevím jestli je to na 100 pro dobře...Můžete někdo potvrdit?? :-)

Offline

 

#8 31. 05. 2014 01:31

FliegenderZirkus
Příspěvky: 544
Škola: RWTH Aachen
Reputace:   25 
 

Re: Odvození dif. rovnice z modelu

↑ karlB:

Ahoj, druhý Newtonův zákon je určitě správná cesta. Jen je potřeba si ujasnit zadání.

1) Má ta spodní deska (x3) nějakou hmotnost? Ani na jednom z obrázků není označena.
2) Používáš vzoreček pro sílu tlumení úměrnou relativní rychlosti, jenže žádné tlumiče na obrázku nakreslené nejsou. Pod pojmem "friction" bych si představil spíše coulombovské tření, tedy konstantní sílu proti směru pohybu. Je to možné? Např. ve třetí rovnici je člen $\frac{b_1}{m} (\dot{x}_2-\dot{x}_1)$, kterému úplně nerozumím. Taky tuhosti pružin jsou značeny různě.
3) Volné délky pružin a tíhovou sílu je výhodné z rovnic eliminovat měřením souřadnic od staticky rovnovážné polohy. V dynamice nás většinou zajímají vlastní frekvence atd., které na statické výchylce (v lineárním případě) nezávisí. Kdyby nás zajímala celkový výchylka, tak je možné k výsledku přičíst konstantní výchylku danou vlastní tíhou...

Online budu zase až v neděli večer.

Offline

 

#9 31. 05. 2014 10:54

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

Ahoj,

díky za reakci, tady jsou odpovědi:

1) deska má hmotnost 10 kg, ostatní závaží m1,m2 mají hmotnost 2kg, ale to potřebuju až když budu tento model modelovat. Hmotnost tedy má.

2) tlumiče tam prý mají být. Jeden si lze představit mezi deskou a m2 a druhý mezi deskou a m1. (tření b1 a b2 působí na desku)  Za desku se považuje celá konstukce (jak spodní rovná část tak hořejší část - i když to na obrázku vypadá oddělené...)

Omlouvám se za chybu - v poslední rovnici posledním členu má být místo 2k2 - 2k3. Pružiny k3 jsou paralelně zapojeny , jejich tuhosti se tedy sčítají.

3)  Délky pružin i tíhová síla tam musí také být (bohužel)..tohle mě bude hodně trápit až budu model přepisovat do state-space reprezentace..

Offline

 

#10 01. 06. 2014 01:26 — Editoval karlB (01. 06. 2014 01:27)

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

ještě ke členu $\frac{b_{1}}{m_{desky}} (x_{2}-x_{1})$ ... nemělo by tam být $\frac{b_{1}}{m_{desky}} (x_{3}-x_{1})$ jako že se jedná o tření mezi deskou (poloha x3) a tělesem m1 ?

samozřejmě s derivací..

Offline

 

#11 01. 06. 2014 20:36

FliegenderZirkus
Příspěvky: 544
Škola: RWTH Aachen
Reputace:   25 
 

Re: Odvození dif. rovnice z modelu

Zdravím a posílám obrázek, ke kterým rovnicím jsem dospěl já a jak. Je potřeba si ale uvědomit, že předpoklad tlumení úměrného rychlosti nebude v případě suchého tření splněn ani přibližně! Musely by tam skutečně být kapalinové tlumiče.
//forum.matweb.cz/upload3/img/2014-06/47655_Soustava_tri_hmot.jpg

Offline

 

#12 02. 06. 2014 11:11

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

Zdravím Vás,

mockrát Vám děkuji! Moc mi to pomohlo :-)

Rovnice jsem vyzkoušel namodelovat a funguje to jak má.

Ještě malýá poddotaz: Kdybych chtěl systém linearizovat, jak by vypadaly vstupy pro linearizaci? Myslím, že by vstupem mhl být moment sil a výstupem polohy x1, x2, x3 ale nejsem si tím jistý..

Offline

 

#13 02. 06. 2014 12:09

FliegenderZirkus
Příspěvky: 544
Škola: RWTH Aachen
Reputace:   25 
 

Re: Odvození dif. rovnice z modelu

Já bych řekl, že ve smyslu definice na Wikipedii je ta soustava diferenciálních rovnic lineární. Co se potom myslí linearizací?

Zase mě napadá, že jsme měli tření modelovat jako coulombovské, protože pak by rovnice skutečně byly nelineární. Směr třecí síly se pak totiž skokově mění podle znaménka rychlosti.

Taky mi není jasný pojem moment sil, když máme jednorozměrný model (hmotné body se pohybují po přímce).

Offline

 

#14 03. 06. 2014 17:27 — Editoval karlB (04. 06. 2014 12:06)

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

Ano, máte pravdu lineární je. Linearizaci nechte být, týká se modelování v matlabu a už jsem na to přišel. Vstupem je velikost síly a výstupem jsou polohy x1,x2,x3. Linearizace se dělá kvůli řiditelnosti systému. Více o linearizaci je možné nalézt zde

Pro zmíněnou linearizaci potřebuji zjistit rovnovážný stav tohoto modelu.
Výpočet rovnovážného se provede z odvozených rovnic (derivace jsou v rovnovážném stavu nulové)

Jen si nejsem jistý jak..

znamená to nahrazení všech derivací v rovnicích nulou. Tedy z první rovnice zbyde:

$0=M_{1}g-k_{1}(x_{1}-x_{2})$

a takto vyjádřím i ostatní, dostanu soustavu 3 rovnic o třech neznámých a výsledné vyjádření x1, x2, x3 je onen hledaný rovnovážný stav?

Offline

 

#15 07. 06. 2014 09:31 — Editoval FliegenderZirkus (07. 06. 2014 09:32)

FliegenderZirkus
Příspěvky: 544
Škola: RWTH Aachen
Reputace:   25 
 

Re: Odvození dif. rovnice z modelu

Neměl jsem čas, tak pardon, že odpovídám až teď. Ano ve statice jsou rychlosti i zrychlení nulová, zbydou jen síly v pružinách a tíhy. S napsanou rovnicí souhlasím.

No a my chceme tu soustavu regulovat tak, aby zůstala v rovnovážné poloze? A kde bude nějaký akční prvek? Nejlepší by bylo napsat kompletní zadání.

Offline

 

#16 07. 06. 2014 10:33

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

Tady je celé zadání (je takové strohé, musel jsem se doptávat)

Model „Mechanický systém s pružinami“

    1) Pro daný systém sestavte
      -  model pomocí SimMechanics
      -  model pomocí diferenciálních rovnic (odvoďte Newton-Eulerovou metodou)
    2) Porovnejte výsledky obou metod pro různé počáteční podmínky, model v Simulinku realizujte pomocí bloku   State-Space
    3) Určete takové počáteční podmínky, aby byl systém v rovnovážném stavu
    4) Vykreslete frekvenční charakteristiku sytému, jestliže vstupem je síla působící na vrcholek konstrukce ve svislém směru a výstupem je poloha x1 (závaží M1)
    5) Určete alespoň 1 konfiguraci senzorů polohy a/nebo rychlosti takovou, že systém bude
        -  pozorovatelný
        -  nepozorovatelný

m =10 kg, M_1 = M_2 = 2 kg, k_1 = 30; k_2 = 50; k_3 = 40; b_1 = b_2 = 0,7. Klidové délky pružin k_1; k_2; k_3 jsou a_1 = a_2 = a_3  = 5


Ta rovnovážná poloha je důležitá pro linearizaci, která se v matlabu provádí v rovnovážném pracovním bodě, stavu.
Linearizací získám matice A,B,C,D které popisují stavový prostor a následně určitým blokem v Simulinku (State-Space) je možné původně namodelovaný model (SimMechanics) nahradit právě tímto linearizovaným - State-Space.
Nevím jestli máte nějaké zkušenosti s modelováním v Matlabu - SimMechanics a Simulink, pokud by Vás to zajímalo, mohu model zaslat.

Offline

 

#17 07. 06. 2014 12:23

FliegenderZirkus
Příspěvky: 544
Škola: RWTH Aachen
Reputace:   25 
 

Re: Odvození dif. rovnice z modelu

Automatické řízení jsem měl celkem dávno, ale něco málo si pamatuju. Jestli máš konkrétní dotaz, tak zkusím odpovědět.

Offline

 

#18 09. 06. 2014 10:37

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

Ještě bych jeden konkrétní dotaz měl:

Týká se posleního bodu - řiditelnosti.

Dejme tomu, že mám k dispozici matice popisující stavový prostor tj. A, B, C, D..

Matici pozorovatelnosti (O) získám násobením matic A a C.


Má-li matice O plnou hodnost, je systém řiditelný..Teď už jen zjistit pro které vektory (kombinace senzorů polohy) ano a pro které ne.

Vím, že by pro konfiguraci senzorů polohy (vektoru x) mělo platit:

$O(A,C) x(t_{0})=y(t_{0})$

X bude zřejmě vektor, který hledáme ale nevím co dosadit za $y(t_{0})$

Offline

 

#19 09. 06. 2014 12:11

FliegenderZirkus
Příspěvky: 544
Škola: RWTH Aachen
Reputace:   25 
 

Re: Odvození dif. rovnice z modelu

karlB napsal(a):

Ještě bych jeden konkrétní dotaz měl:

Týká se posleního bodu - řiditelnosti.
(...)
Matici pozorovatelnosti (O) získám násobením matic A a C.
(...)
Má-li matice O plnou hodnost, je systém řiditelný

Tak nevím, řešíme řiditelnost nebo pozorovatelnost? :D

Každopádně si nejsem jistý, jak se má postupovat. Jestliže pozorovatelnost závisí pouze na maticích stavového popisu (resp. na hodnosti nějakého jejich součinu), tak musí poloha senzorů mít vliv na tyto matice, ne? Tomu úplně neorzumím. Třeba se ozve některý z kolegů.

Offline

 

#20 09. 06. 2014 17:48

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

pardon, jedná jen o pozorovatelnost:-)

Vím jen, že matice se kterými se pracuje, by měli zůstat stejné. Zatím jsem se k ničemu, co by mi v tomto směru nějak pomohlo nedopátral.

Offline

 

#21 10. 06. 2014 10:46

karlB
Příspěvky: 28
Škola: VSE
Pozice: student
Reputace:   
 

Re: Odvození dif. rovnice z modelu

Tak už vím jak to je.


Matice/vektor C se volí tak aby po n-tém vynásobení s maticí A byla hodnost výsledné matice pozorovatelnosti nižší než plná hodnost. Teď je otázkou jak např. v matlabu na takový vektor přijít..

Offline

 

Zápatí

Powered by PunBB
© Copyright 2002–2005 Rickard Andersson