Analityczne i graficzne kryteria stabilności liniowych UAR


0x01 graphic

Akademia Górniczo - Hutnicza

W Krakowie

Wykonał:

Dariusz

Krakowski

KATEDRA AUTOMATYKI , NAPĘDU I URZĄDZEŃ PRZEMYSŁOWYCH AGH

Wydział:

EAIiE

Rok akad.:

1998 / 99

Rok studiów:

II

Kierunek:

Elektrotechnika

Grupa:

3

Temat ćwiczenia:

Analityczne i graficzne kryteria stabilności liniowych UAR.

Data wykonania:

Data zaliczenia:

Ocena:

A. Cel ćwiczenia:

Zapoznanie się istniejącymi kryteriami stabilności układów liniowych oraz ich praktycznym wykorzystaniem przy badaniu stabilności obiektów.

B. Przebieg ćwiczenia:

  1. Kryterium Routha - Hurwitza.

  1. opis transmitancyjny - funkcja roots(), routh(), pzmap():

Równanie charakterystyczne badanego obiektu:

s4 + s3 + s2 + s +1 = 0

pierwiastki równania charakterystycznego:

0.3090 + 0.9511i -0.8090 + 0.5878i

0.3090 - 0.9511i -0.8090 - 0.5878i

0x01 graphic
0x01 graphic

odpowiedź układu bieguny układu

na skok jednostkowy na płaszczyźnie zespolonej

0x01 graphic
0x01 graphic

Wszystkie współczynniki są jednakowego znaku, więc pierwszy warunek jest spełniony.

Następnie sprawdzamy znaki poszczególnych podwyznaczników:

0x01 graphic
0x01 graphic

Jak widać wyznacznik całej macierzy wyszedł ujemny, więc z tego wynika wg kryterium Hurwiztza, że układ jest niestabilny.

0x01 graphic
0x01 graphic

0x01 graphic

0x01 graphic

Jak łatwo zauważyć nastąpiła podwójna zmiana znaku w pierwszym wierszu macierzy a więc układ wg kryterium Rutha powinien posiadać 2 pierwiastki charakterystyczne w dodatniej półpłaszczyźnie.

Układ z takim równaniem charakterystycznym jest układem niestabilnym. Przemawia za tym kryterium Hurwitza jak i Routha.

  1. wyznaczenie „K” dla którego układ ten jest układem stabilnym:

0x01 graphic

Równanie charakterystyczne układu:

s3 + 2s2 + 4s + K = 0 dla K∈(0,20)

Macierz wzmocnień i biegunów dla K zmieniającego się od 0 do 20:

wzmocnienie K: pierwiastki równania charakterystycznego:

0 0 -1.0000 - 1.7321i -1.0000 + 1.7321i

1.0000 -0.8576 - 1.6661i -0.8576 + 1.6661i -0.2848

2.0000 -0.6806 - 1.6332i -0.6806 + 1.6332i -0.6389

3.0000 -0.5000 - 1.6583i -0.5000 + 1.6583i -1.0000

4.0000 -0.3522 - 1.7214i -0.3522 + 1.7214i -1.2956

5.0000 -0.2370 - 1.7946i -0.2370 + 1.7946i -1.5260

6.0000 -0.1443 - 1.8669i -0.1443 + 1.8669i -1.7113

7.0000 -1.8664 -0.0668 - 1.9355i -0.0668 + 1.9355i

8.0000 -2.0000 0.0000 - 2.0000i 0.0000 + 2.0000i

9.0000 -2.1179 0.0589 - 2.0606i 0.0589 + 2.0606i

10.0000 -2.2236 0.1118 - 2.1177i 0.1118 + 2.1177i

11.0000 -2.3198 0.1599 - 2.1717i 0.1599 + 2.1717i

12.0000 -2.4082 0.2041 - 2.2229i 0.2041 + 2.2229i

13.0000 -2.4902 0.2451 - 2.2717i 0.2451 + 2.2717i

14.0000 -2.5667 0.2833 - 2.3182i 0.2833 + 2.3182i

15.0000 -2.6386 0.3193 - 2.3628i 0.3193 + 2.3628i

16.0000 -2.7064 0.3532 - 2.4056i 0.3532 + 2.4056i

17.0000 -2.7707 0.3854 - 2.4468i 0.3854 + 2.4468i

18.0000 -2.8320 0.4160 - 2.4866i 0.4160 + 2.4866i

19.0000 -2.8904 0.4452 - 2.5249i 0.4452 + 2.5249i

20.0000 -2.9463 0.4732 - 2.5621i 0.4732 + 2.5621i

Aby badany przez nas układ należał do układów stabilnych musi on mieć ujemne części rzeczywiste pierwiastków charakterystycznych. Jak z tego wynika i widać z powyższych danych układ jest stabilny dla 1 K 7, dla K=8 układ jest na granicy stabilności a dla K>8 układ się rozbiega i staje się niestabilny.

Położenie biegunów układu dla K (0,20):

0x01 graphic

s3 + 2s2 + 4s + K = 0 dla 0x01 graphic

wg kryterium Hurwitza:

- wszystkie współczynniki muszą istnieć i być tego samego znaku stąd K musi być dodatnie,

- wyznacznik główny macierzy i podwyznaczniki muszą być także dodatnie,

0x01 graphic
0x01 graphic

więc jak widać z powyższych przeliczeń K musi należeć do przedziału (0,7), aby układ był stabilny.

wg kryterium Routha:

0x01 graphic
0x01 graphic

Jak widać z powyższych obliczeń, aby układ był stabilny to wzmocnienie K musi należeć do przedziału < 0,8 > .

Stąd nasuwa się następujący wniosek, że kryterium Routha inaczej traktuje granice stabilności bowiem przyjmuje za układ stabilny układ na granicy stabilności.

Położenie biegunów układu stabilnego czyli dla K∈(0,7)

0x01 graphic

c) opis w przestrzeni stanu - funkcje Matlaba poly() oraz eig()

0x01 graphic

Równanie charakterystyczne powyższej macierzy stanu ma następującą postać (funkcja poly):

s3 - 6s2 - 7s - 52 = 0

pierwiastki równania charakterystycznego:

-0.8821 + 2.4330i

-0.8821 - 2.4330i

7.7642

bieguny układu odpowiedź układu na płaszczyźnie zespolonej na skok jednostkowy

0x01 graphic
0x01 graphic

Jak także widać z tej metody układ jest niestabilny. Zarówno widać to z położenia pierwiastków równania charakterystycznego jak i z odpowiedzi na skok jednostkowy (układ się nam rozbiega do nieskończoności).

  1. Kryterium Nyquista.

Określenie kiedy (dla jakiego K) poniższy układ (opisany poniższą transmitancją) jest stabilny, niestabilny oraz marginalnie stabilny:

0x01 graphic

Jeżeli układ otwarty jest stabilny, a jego charakterystyka amplitudowo-fazowa nie obejmuje punktu (-1, j0) przy zmianie pulsacji ω od minus nieskończoności do plus nieskończoności, to układ zamknięty także będzie stabilny. W przeciwnym przypadku układ zamknięty będzie układem niestabilnym. Jeśli charakterystyka amplitudowo-fazowa układu otwartego przechodzi przez punkt (-1, j0), to układ zamknięty jest na granicy stabilności.

W naszym przypadku układ otwarty znajduje się na granicy stabilności, więc można by przyjąć, że jest prawie stabilny.

Wybrane charakterystyki:

0x01 graphic
0x01 graphic

0x01 graphic
0x01 graphic

Na podstawie powyższych charakterystyk można wywnioskować, że układ jest stabilny dla K∈(1,7), dla K=8 układ znajduje się na granicy stabilności, natomiast dla K > 8 układ staje się układem niestabilnym. Stabilność układu można ocenić także po rozmieszczeniu pierwiastków charakterystycznych. Jak widać układ otwarty jest na granicy stabilności, natomiast układ zamknięty jest stabilny dla K < 8 , a niestabilny dla K > 8 (obejmuje zasięgiem punkt (-1+j0) ).

0x01 graphic
0x01 graphic

0x01 graphic
0x01 graphic

Inną metodą stabilność układu można określić przy pomocy charakterystyk Bodego. Jeśli układ otwarty jest stabilny, a jego charakterystyka amplitudowa logarytmiczna przecina oś 0 dB przy pulsacji mniejszej od pulsacji przecięcia osi -π przez charakterystykę fazową, to układ zamknięty jest stabilny. W przeciwnym przypadku układ zamknięty jest niestabilny.

dla K=5 dla K=8

0x01 graphic
0x01 graphic

dla K=13

0x01 graphic

Na podstawie charakterystyk Bodego można stwierdzić, że badany układ jest stabilny dla K<8, niestabilny dla K>8 a natomiast dla K=8 układ jest na granicy stabilności.

  1. zbadanie stabilności układu:

0x01 graphic

W powyższym układzie: T1 = 1 [s] T2 = 10 [s] k = 1

Odpowiedź na skok jednostkowy Charakterystyka amplitudowa i fazowa

0x01 graphic
0x01 graphic

Rozmieszczenie biegunów i zer badanego układu:

0x01 graphic

Układ znajduje się na graniacy stabilności dla Ti = 0,91 [s]. Jak widać na wykresie odpowiedzi układu na skok jednostkowy na wyjściu układu pojawiają się stałe oscylacje. Także można to odczytać z położenia pierwiastków równania charakterystycznego na płaszczyźnie zmiennej zespolonej (widzimy, że układ posiada pierwiastki na osi urojonej). Przy zmniejszeniu Ti pojniżej 0,91 [s] układ staje się niestabilny, a jeżeli Ti będzie większe niż 0,91 [s] to ten układ będzie stabilny.

  1. Stabilizacja układu niestabilnego.

0x01 graphic

Transmitancja obiektu G(s):

0x01 graphic

Jest to układ niestabilny, gdyz posiada 2 sprzężone pierwiastki zespolone w prawej półpłaszczyźnie. Do jego kompensacji użyto kompensatora o transmitancji:

0x01 graphic

Na podstawie poniższego wykresu można ocenić czy układ z takim kompensatorem będzie kiedykolwiek stabilny.

0x01 graphic

Jak widać z wykresu nie ma możliwości wprowadzenia takiego wzmocnienia kompensatora, aby układ stał się stabilny. Więc możemy łatwo wywnioskować, że układ z takim kompensatorem nigdy nie będzie stabilny, ponieważ nigdy wszystkie pierwiastki równania charakterystycznego nie będą leżeć w lewej półpłaszczyźnie.

Wnioski:

Bardzo ważnym problemem w układach automatycznej regulacji jest stabilność badanego układu. Jest to podstawowe wymaganie stawiane układowi. Aby zbadać stabilność układu wcześniej musieliśmy jego pierwiastki równania charakterystycznego. Lecz gdy układ opisuje równanie wyższego rzędu wyliczenie pierwiastków może stać się prawie że niemożliwe. I tutaj nam z pomocą przychodzą poznane na ćwiczeniu kryteria stabilności. W ćwiczeniu tym poznaliśmy i porównaliśmy najczęściej stosowane kryteria stabilności takie jak kryterium Routha-Hurwitza czy Nyquista. Dzięki nim nie musimy znać wartości pierwiastków równania charakterystycznego, tylko wystarczy nam poznać ich położenie na wykresie zmiennej zespolonej. Kryteria określaja nam czy badany przez nas układ jest stabilny oraz co trzeba zrobić, aby stał się on stabilny. Przekonaliśmy się również, że poprzez właściwy dobór odpowiedniego kompensatora możliwe staje się ustabilizowanie układu niestabilnego.

- 2 -



Wyszukiwarka

Podobne podstrony:
Analityczne i graficzne kryteria stabilności liniowych UAR
LAB08 Analityczne i graficzne kryteria stabilno ci liniowych UAR
Wykład 6 Stabilność liniowych układów automatyki (2013)
3 Opracowanie wyników pomiaru metodą analityczno graficzną (Langa)
07 Stabilność liniowych stacjonarnych układów sterowania
Elektornika teoria Stabilizatory liniowe
Definicja i kryteria stabilności
obiegi skladow analityczna i graficzna-moje skoncz, Szkoła, Semestr 5, Technologia pasażerskich prze
L5 Badanie stabilności liniowego układu 3 rzędu z opóźnieniem Wpływ wartości opóźnienia na stabilno
W8 Stabilizatory Liniowe
Stabilność liniowych U
Chorowski,podstawy automatyki, Hurwitza kryterium stabilności
Chorowski,podstawy automatyki, kryteria stabilności
Wykład 6 Stabilność liniowych układów automatyki (2013)
WYKLAD 06 STABILIACJA LINIOWA
07 Stabilność liniowych stacjonarnych układów sterowania
Stabilizatory liniowe

więcej podobnych podstron