Wykład KCh3.
3. Równanie ruchu układu napędowego- część III
3.1. Człony M0 oraz Me występujące w równaniu ruchu układu napędowego
3.2. Równanie (3.1) - statyczna charakterystyka silnika
3.3. Równanie (3.1) - moment oporowy zależny od kąta obrotu silnika
3.4. Stabilny i niestabilny punkt pracy układu napędowego
3.5. Uwagi wstępne dotyczące układów o więcej niż jeden stopniach swobody
3.1. Człony M0 oraz Me występujące w równaniu ruchu układu napędowego
Z poprzedniego wykładu wiadomo, że symbol M0 występujący w równaniu ruchu
(3.1)
oznacza moment oporowy maszyny roboczej. Pod tą nazwą kryje się zarówno moment użyteczny (niezbędny do wykonania procesu technologicznego) jak i suma momentów niepożądanych wynikających z tarcia w maszynie roboczej oraz w silniku. Typowe charakterystyki maszyn roboczych pokazano w sposób jakościowy na rys.3.1. Linia 1
odpowiada maszynom, w których moment oporowy jest stały (np pompom o stałej wydajności, maszynom wyciągowym, walcarkom). Linia 2 charakteryzuje wentylatory, pompy odśrodkowe i inne maszyny, w których opory ruchu są wprost proporcjonalne do kwadratu prędkości obrotowej. Linia 3 odpowiada - przykładowo - urządzeniom odwijającym, w których jest zachowana stała siła naciągu oraz stała prędkość liniowa taśmy lub drutu; maszyny takie są stosowane między innymi w przemyśle papierniczym.
Wyróżnia się bierne i czynne momenty (siły) oporowe. Momenty (siły) oporowe bierne są związane z wykonywaniem procesu technologicznego i pokonywaniem tarcia. Mówiąc inaczej: praca momentu (siły) oporowego biernego jest pracą użyteczną lub pracą tarcia. Moment (siła) oporowy bierny ma zawsze zwrot przeciwny do zwrotu prędkości. Obrazują to na rys.3.2 linie A1 i A2.
Czynny moment (siła) oporowy jest niezależny od zwrotu prędkości. Obrazuje to linia B na rys.3.2. W przypadku czynnego momentu (siły) oporowego należy wyróżnić dwie sytuacje, które łatwo wyjaśnić na przykładzie windy. Jeżeli winda porusza się do góry, to moment oporowy przeciwdziała ruchowi. Jeżeli natomiast winda zjeżdża w dół, to moment oporowy staje się momentem czynnym; „pomaga” temu ruchowi, a silnik może oddawać energię do źródła. Innym przykładem jest siła spowodowana odkształceniem sprężystym elementu maszyny. W obydwu tych przypadkach mamy do czynienia z energią potencjalną. W przypadku windy - z energią potencjalną związaną z polem grawitacyjnym, w przypadku elementu maszyny - z energią potencjalną odkształceń sprężystych. Energia ta jest najpierw akumulowana, a następnie może być oddana.
Dyskusję na temat członu Me występującego w równaniu ruchu układu napędowego należy rozpocząć od wyróżnienia statycznych i dynamicznych charakterystyk silnika elektrycznego. Charakterystyka statyczna silnika to taka, którą uzyskano na hamowni zmieniając bardzo wolno obciążenie silnika. „Bardzo wolno” oznacza, że zmiana obciążenia nie powoduje ani zauważalnych momentów dynamicznych (przyspieszeń), ani procesów przejściowych w silniku. Charakterystykę statyczną silnika indukcyjnego pokazano schematycznie na rys.3.3. W przybliżeniu jest ona
opisana tzw. wzorem Klossa
(3.2)
gdzie s oznacza tak zwany poślizg, równy
(3.3)
- prędkość kątowa,
- prędkość kątowa synchroniczna.
Symbole ωk , sk , Mk oznaczają krytyczna prędkość kątową, krytyczny poślizg i krytyczny (maksymalny) moment silnika. Poślizg krytyczny można obliczyć ze wzoru
gdzie Mn oznacza moment znamionowy. Dla poślizgów znacznie mniejszych od poślizgu krytycznego zachodzi zależność sk/s >>s/sk , dzięki czemu wzór Klossa można napisać w postaci
. (3.4)
Podstawiając (3.3) do (3.2) a następnie do (3.4) otrzymuje się
(3.2a)
oraz
(3.4a)
Charakterystyka statyczna silnika innego niż indukcyjny ma oczywiście odmienną postać, jednak jest to zawsze zależność algebraiczna wiążąca moment Me z prędkością kątową silnika.
3.2. Równanie (3.1) - statyczna charakterystyka silnika
Rozpatrzmy najpierw najprostszy przypadek równania ruchu. Odpowiada on założeniu, że obydwa momenty, Me i M0, mają wartości stałe. Zastosujmy to założenie do oszacowania czasu rozruchu silnika indukcyjnego. Przyjmijmy, że moment M0 jest znany i mniejszy od momentu rozruchowego Mr (rys.3.3). W czasie rozruchu prędkość kątowa wzrasta od
do
, gdzie
oznacza prędkość ruchu ustalonego, czyli sytuację, w której Me = M0. Korzystając z charakterystyki silnika można powiedzieć, że rozruch rozpoczyna się od punktu (0, Mr) i kończy w punkcie (
, M0). Symbol Mr oznacza moment rozruchowy silnika. Przyjmijmy, i jest to istota uproszczenia, że w czasie rozruchu moment silnika jest stały i wynosi
.
Równanie (3.1) przyjmuje postać
Zauważmy, że równanie to umieją (ściśle: powinni umieć) rozwiązać uczniowie szkół średnich, gdyż opisuje ono ruch obrotowy jednostajnie przyspieszony. My rozwiążemy to równanie nieco „mądrzej”, czyli, jak niżej
gdzie C jest stałą, która należy obliczyć z warunku początkowego: dla t = 0 jest ω = 0. Uzyskuje się
Rozruch, którego czas wynosi tr , kończy się gdy prędkość osiągnie wartość
. Stąd wynika czas rozruchu
(3.5)
Zadanie
O jaki kąt obróci się w czasie rozruchu wał silnika?
Uzyskanemu rozwiązaniu można zarzucić, że w rzeczywistości moment silnika nie jest stały, a więc przyspieszenie też nie jest stałe. Rozwiązanie to ma jednak ogromną zaletę: jest proste. W technice bardzo często przedkłada się prostotę rozwiązania nad jego dokładność.
Zastosujmy dokładniejszy opis analityczny silnika indukcyjnego, czyli równanie (3.4a), do rozwiązania problemu zdefiniowanego poniżej. Niech w układzie napędowym, który znajduje się w stanie ustalonym zmieni się, w chwili t = 0, moment oporowy z M01 na M02. Z równania (3.4a) wynika, że momentowi M01 odpowiada prędkość kątowa
gdzie
Równanie ruchu ma więc postać
czyli
Oznaczając
otrzymujemy równanie
którego rozwiązaniem jest
Korzystając z warunków początkowych (dla t = 0 jest ω = ω1) otrzymuje się
czyli
Zauważmy, że gdy
to
, co pozwala uzyskać rozwiązanie końcowe w postaci
(3.6)
Zadanie
Silnik indukcyjny ma następujące parametry: Moc znamionowa 17 kW, obroty znamionowe 1460 obr/min, moment rozruchowy 128 Nm, moment krytyczny 220 Nm. Silnik jest początkowo obciążony momentem 50 Nm. W pewnej chwili obciążenie wzrasta gwałtownie do wartości 125 Nm. Narysować wykres obrazujący zależność prędkości kątowej od czasu aż do chwili, w której osiągnięta prędkość będzie wynosić 99,5% prędkości ω2. Jakie wnioski wynikają z rozwiązania?
Zależność (3.5) dotycząca rozruchu została wyprowadzona przy założeniu, że moment silnika jest stały, a zależność (3.6) dotycząca skokowej zmiany momentu oporowego - przy założeniu, że moment silnika jest liniową funkcją jego poślizgu. Można byłoby zamiast tych założeń próbować wprowadzić do równania ruchu równanie Klossa, czyli wzór (3.2). Jednak zauważmy, że zgodnie z równaniem Klossa moment silnika jest nieliniową funkcja poślizgu, a więc nieliniową funkcją prędkości kątowej. Jest to poważne utrudnienie, gdyż równanie ruchu przestaje być równaniem liniowym. Większość równań nieliniowych nie ma rozwiązań, które można by było - tak jak (3.5) i (3.6) zapisać w postaci wzoru analitycznego. Nieliniowe równanie różniczkowe można rozwiązać, ale tylko numerycznie, a rozwiązanie numeryczne dotyczy konkretnego przypadku. Z tego powodu - jeżeli to tylko nie jest niezbędne - należy unikać założeń, które prowadzą do nieliniowych równań ruchu. Zalecenie to dotyczy nie tylko napędu elektrycznego.
Zadanie
Przyjmijmy, że ktoś „uparł się” wprowadzić wzór Klossa w postaci (3.2) lub (3.2a) do równania ruchu i oceńmy to postępowanie w przypadku analizy skokowej zmiany momentu oporowego.
Zadanie
Ocenić celowość wprowadzenia wzoru Klossa do analizy rozruchu silnika.
3.3 Równanie (3.1) - moment oporowy zależny od kąta obrotu silnika
W niektórych maszynach moment oporowy zależy od kata obrotu wału wejściowego. Są to najczęściej maszyny, których element roboczy jest napędzany mechanizmem
dźwigniowym, przykładowo, mechanizmem korbowym (rys.3.4). Typowymi przykładami są sprężarka tłokowa, prasa korbowa, itp. Łatwo zauważyć, że moment M1 niezbędny do pokonania siły F działającej na tłok (w prasie - na suwak) zależy od kąta α określającego położenie kątowe korby, czyli
Równanie ruchu należy w tym przypadku napisać
wprowadzając po lewej stronie
zamiast
. Jeżeli funkcja M1(α) ma postać jak na rys.3.5, to silnik przez cały czas pracuje w warunkach przejściowych. Złagodzenie warunków pracy silnika osiąga się stosując pomiędzy silnikiem i napędzaną maszyną koło zamachowe. Zwiększa się w ten sposób moment bezwładności J, a więc zmniejsza ujemne przyspieszenie (opóźnienie) kątowe
. W przypadku zmian momentu M0 jak na rys.3.5, powyższe równanie jest równaniem nieliniowym. Można je zlinearyzować zastępując, w dość żmudny sposób, krzywe odcinkami prostych.
3.4. Stabilny i niestabilny punkt pracy układu napędowego
Stabilny układ (niekoniecznie napędowy) to taki, w którym małe odchylenie od stanu, w którym znajduje się układ, jest samoczynnie likwidowane i układ wraca do stanu wyjściowego. Bardzo często pojęcie stabilności wyjaśnia się rozpatrując zachowanie kulki na powierzchni wklęsłej (stabilne położenie kulki) oraz na powierzchni wypukłej (położenie niestabilne). Stabilność układu napędowego można natomiast wyjaśnić posługując się wykresami momentów (silnika Me i oporowego M0) w funkcji prędkości kątowej. Rozpatrzmy
dwie różne sytuacje pokazane na rys.3.6a i 3.6b i przeanalizujmy odchylenia układu od punktu równowagi, czyli od punktu S (rys.3.6a) przecięcia charakterystyk Me(ω) i M0(ω). Załóżmy że moment oporowy wzrósł o ΔM0 . Zgodnie z charakterystyką M0(ω), wzrostowi temu odpowiada wzrost prędkości kątowej. Powinien mu towarzyszyć wzrost momentu silnika Me. Tak jednak nie jest, gdyż wraz ze wzrostem prędkości kątowej moment silnika zmniejsza się. Oznacza to, że w punkcie P układ napędowy nie może pracować i powróci do stabilnego punktu pracy S. Inaczej jest w sytuacji pokazanej na rys.3.6b. Jeżeli moment oporowy na chwilę by się zmniejszył, to układ powinien przejść do punktu R. Praca układu w tym punkcie nie jest jednak możliwa, ponieważ zwiększonej prędkości kątowej towarzyszy większy moment silnika. Nadwyżka momentu silnika w porównaniu z momentem oporowym spowoduje wzrost obrotów i dalsze zmniejszenie momentu oporowego. Oznacza to, że punkt N jest niestabilnym punktem pracy silnika.
Z powyższego rozumowania wynika, że aby punkt pracy układu napędowego był stabilny, musi być spełniony warunek
.
Powyższy warunek pozwala stwierdzić, że w przypadku silnika indukcyjnego (rys.3.7) punkt S1 jest punktem pracy stabilnej. Podobnie stabilnym punktem pracy jest punkt S2 jednak w punkcie tym długotrwała praca silnika nie jest możliwa z uwagi na duży prąd płynący w uzwojeniach silnika, który spowodowałby niedopuszczalny wzrost temperatury tych uzwojeń.
Przeanalizujmy jeszcze jedną sytuację zakładając, że na silnik, pracujący stabilnie w punkcie S1, zaczyna działać zwiększony moment oporowy; pokazuje to w sposób umowny strzałka skierowana od punktu S1 w górę. Na rys.3.7 widać, że silnik odpowiadając spadkiem prędkości i zwiększeniem momentu przechodzi do nowego punktu stabilnej pracy. Dzieje się tak aż do osiągnięcia momentu maksymalnego (krytycznego). Jeżeli moment oporowy nadal wzrasta, to silnik po prostu się zatrzymuje. W analizie tej pominęliśmy nagrzewanie się silnika. W rzeczywistości silnik nie może długo pracować obciążony momentem większym od momentu znamionowego (nominalnego).
3.5. Uwagi wstępne dotyczące układów o więcej niż jeden stopniach swobody
Rozpocznijmy od definicji: „Stopień swobody jest to liczba niezależnych zmiennych opisujących jednoznacznie stan modelu układu fizycznego.” Z definicji wynika, że wszystkie analizowane dotychczas modele układów napędowych miały jeden stopień swobody. Oscylator harmoniczny, którego model pokazano na rys.1.1 (Wykład KCh1), ma też jeden stopień swobody.
W przytoczonej powyżej definicji warte omówienia jest bardzo często używane w nauce (i nie tylko) pojęcie „model”. Zacznijmy od pytania - żartu: Co przedstawia rys.3.8?
Odpowiedź „kota” jest błędna. Odpowiedzią prawidłową jest: „Na rys.3.8 przedstawiono model graficzny kota”. Żart ten można uznać za wprowadzenie do rys. 3.9, który przedstawia klasyfikację modeli. Wszystkie rozpatrywane dotychczas modele należały do kategorii „symboliczny” i pod-kategorii „matematyczny”, zarówno modele statyczne jak i dynamiczne. Definicja modelu jest następująca: „A model is a simplified representation of a system (or process, or theory) intended enhance our ability to understand, predict, and possibly control the behavior of the system”. Zdanie to zaczerpnięto z ksiażki F.Neelamkavil'a “Computer simulation and modeling” celowo bez tłumaczenia.
Ważne jest podkreślenie, że „najlepszy” model nie istnieje. Model powinien być najprostszy z możliwych a jednocześnie wystarczający do opisania zjawisk, które dzięki modelowi mają być zrozumiane lub przewidziane. Jest to wskazówka bardzo ogólna, ale trudno podać lepszą. Ogólność powyższej wskazówki jest właśnie powodem, że w budowie modeli, szczególnie matematycznych, gra doświadczenie. Do tworzenia modeli o prostej strukturze „zmusza” niepełna znajomość wartości liczbowych wielkości występujących w równaniach opisujących model. Istotnie. Stosunkowo łatwo można zmierzyć masę, trudniej
moment bezwładności i sprężystość, jeszcze trudniej cechy tłumiące elementu układu podobnie jak wartości liczbowe wielkości występujących w schematach zastępczych silników elektrycznych.
Pamiętając o powyższych uwagach, wróćmy do układów o więcej niż jeden stopniach swobody rozpatrując układ pokazany na rys.3.10. W skład układu wchodzą dwie masy: m1 i m2. Są one połączone sprężyną k2: sprężyna k1 łączy masę m1 z nieruchomą ścianą. Zachowanie
układu będzie znane jeżeli w każdej chwili znane będą położenia mas m1 i m2. Położenia te określone są dwoma współrzędnymi x1 i x2, co oznacza, ze układ ma dwa stopnie swobody. Na rys.3.10 narysowane zostały siły działające na obydwie masy. Chcąc uzyskać równania ruchu wystarczy zsumować „narysowane' siły i sumę przyrównać do zera.
Układ o dwóch stopniach swobody jest opisany układem dwóch równań różniczkowych i wymaga znajomości położeń i prędkości obydwu mas w chwili uznanej za początkową, czyli
Podobnie, dwa stopnie swobody ma układ pokazany na rys.3.11 składający się z dwóch obracających się elementów (1 i 2) połączonych elastycznym wałem. Niech elementy 1 i 2 cechują się momentami bezwładności J1 i J2. Wał łączący elementy ma sztywność kątowa k12 (sztywność skrętna mierzona jest w Nm/rad). Odkształcenie sprężyste wału powoduje, że na elementy 1 i 2 działają momenty skręcające k12(α1-α2). równe co do wartości bezwzględnej, lecz przeciwnie skierowane.
Na element 1 działa także moment MS, zaś na element 2 - moment oporowy M02. Stan układu wymaga znajomości dla dowolnej chwili t zarówno kąta obrotu α1 elementu 1 jak i kąta obrotu α2 elementu 2.
Zadanie
Napisać równania ruchu dla układu pokazanego na rys.3.11. Podać, w ogólnej postaci, warunki początkowe.
Trzy masy połączone elementami sprężystymi lub tłumiącymi maja trzy stopnie swobody, cztery masy - cztery, itd. Bryła (np. pręt, płyta, układ prętów, itp.) wykonana z rzeczywistego, sprężystego materiału ma nieskończenie wiele stopni swobody.
9
ω
A1
Rys.3.2
B
Me
Mk
A2
M0
ω
Rys.3.1
Rys.3.11
J1d2α1/dt2
MS
k12(α1-α2)
J2d2ω2/dt2
2
J2
0
3
2
1
0
ω
M0
1
J1
ω1
ω2
k12
M02
ωS
0
ω<0
ω>0
0
s
1
Rys.3.3.
ωk
sk
Mr
ω
Me
α
Rys.3.5.
6π
4π
2π
0
M1 (=-M0)
α
Rys.3.4
M1
F
ω
M0
ω
M0
Me
Rys.3.6
(a)
(b)
ΔM0
ΔM0
S
N
P
R
Rys.3.8
ω
Me
0
ωS
M02
Model
S2
M01
S1
Rys.3.7
Fizyczny
Symboliczny
Mentalny
Składa się z elementów, które można dotknąć.
Zależności matematyczne, symboliczne lub logiczne, mapy, wykresy, rysunki, itp.
Istnieje tylko w umyśle; niemożliwy do przekazania.
Matematyczny
Niematematyczny
Statyczny
Dynamiczny
Równania algebraiczne
Równania różniczkowe
Rys.3.9
F
m2 (d2x2/dt2)
k2(x2-x1)
m1 (d2x1/dt2)
k1x1
F2
m2
m1
k1
x1
+x1, + dx1/dt, +d2x1/dt2 +x2, +dx2/dt, +d2x2/dt2
k2
x2
k2(x2-x1)
Rys.3.10