Artykuł z serii: Kurs regulacji PID - Wirtualne Laboratorium
Rozdział 31.1 Wstęp
W rozdziałach 26..30 omawialiśmy zasadę działania regulatorów P, PD, I, PI i PID. Sądzę, że dobierając ręcznie metodą prób i błędów nastawy Kp, Ti i Td dobrze zrozumiałeś ich role. Przy okazji poznałeś pierwszą metodę doboru nastaw. Metoda prób i błędów, używaną także w innych dziedzinach dalekich od automatyki.
Przykład
Jesteś bardzo ostrożny jadąc pierwszy raz na nartach. Skręcasz dużymi zakosami, długo hamujesz, przewracasz się. Podobny jesteś do obiektu z regulatorem PID z bardzo ostrożnymi nastawami. Małe Kp i Td, całkowania prawie nie ma czyli bardzo duże Ti. Ucząc się nabierasz doświadczenia, ruchy są coraz bardziej precyzyjne. A co to oznacza w języku automatyki? To,że poprawiasz nastawy swojego prywatnego regulatora Kp, Ti, Td zakodowanego w zwojach mózgowych. W końcu śmigasz slalomem między tyczkami. Twój regulator jest już prawie optymalny. Podejrzewam, że Radwańska, Wiliams czy inna Woźniacki, najlepsze tenisistki świata, mają nastawy swoich regulatorów bardzo zbliżone do siebie. Ale tylko mistrzyni ma trochę lepsze. Kończę z tą naiwną i mocno uproszczoną filozofią.
Na czym to stanęło? Że pierwszym podejściem do strojenia regulatorów (czyli dobierania Kp, Ti i Td) jest metoda prób i błędów–> inaczej metoda ręczna. Przejdźmy teraz do metod poważnych.
Przedstawię trochę staromodne metody wymyślone jeszcze w czasie II wojny światowej przez niejakiego Zieglera i Nicholsa:
– Metoda odpowiedzi skokowej
– Metoda cyklu granicznego
Każda jest przykładem innej szkoły. Pierwsza wymaga opisu matematycznego, czytaj transmitancji Go(s) obiektu. Druga tylko zbadania obiektu, a ściślej zmierzenia okresu drgań przy pewnym wzmocnieniu tzw. krytycznym. Jest więc mniej wybredna co do dokładnej znajomości obiektu.
Rozdział 31.2 Metoda odpowiedzi skokowej
Rozdział 31.2.1 Wstęp
Wymaga znajomości modelu matematycznego danego obiektu.
Model możemy uzyskać poprzez:
– dokładną analizę matematyczną poprzez znajomość procesów fizyko-chemicznych. Jest to oczywiście zadanie dla ambitnych.
– nie wnikamy w żadne termodynamiki, tylko badamy odpowiedź na skok jednostkowy. Zadanie mniej ambitne i łatwiejsze i dlatego ją omówimy.
Wystąpią w niej 4 etapy.
1-Badanie odpowiedzi obiektu na skok w układzie otwartym i wyznaczanie transmitancji przybliżonej Gp(s) jako inercji T i opóźnienia To
2-Wyznaczenie z odpowiednich tabel nastaw optymalnych Kp, Ti i Td
3-Sprawdzenie odpowiedzi na skok układu zamkniętego
4-Próby znalezienia lepszych nastaw metodą ręczną
Rozdział 31.2.2 Etap 1–Badanie odpowiedzi obiektu na skok w układzie otwartym
Wywołaj PID/17_dobór_nastaw_PID/01_obiekt_4T.zcos
Rys. 31-1
Jest to typowy człon wieloinercyjny, tu czteroinercyjny
Na prawdziwym obiekcie nie znasz oczywiście tej transmitancji. Domyślasz się tylko, że jest to człon wieloinercyjny. Jakie jest to „wielo”, jakie są stałe czasowe? Tego nie wiesz. Dodam tylko, że dużo procesów, zwłaszcza w przemyśle chemicznym, ma taki charakter .
Wciśnij „start”
Rys. 31-2
Skup się na czarnym skoku x(t) i czerwonej odpowiedzi y(t). Widzisz charakterystyczne kolanko, inaczej punkt przegięcia. Przez ten punkt została poprowadzona styczna i z tej konstrukcji zostały wyznaczone parametry K, T i To transmitancji zastępczej Gp(s). Temat ten został dokładnie omówiony w rozdz. 14 Rys. 14-9.
Rys. 31-3
Wyznaczone parametry obiektu K T i To
Rozdział 31.2.3 Etap 2-Wyznaczenie z odpowiednich tabel nastaw optymalnych Kp, Ti i Td
Wg Pana Zieglera jest kilka optymalnych kombinacji Kp, Ti, Td dla tej samej odpowiedzi skokowej. Jak to? Przecież optymalny zestaw Kp, Ti i Td. Zwycięzca, może być tylko jeden. Święta prawda, ale jedno optimum jest dla najkrótszego czasu bez przeregulowań. Drugie to też najkrótszy czas, ale np. z przeregulowaniem 20%, a trzecie to kryterium minimum całki z kwadratu uchybu. Już sama nazwa trzeciego budzi niechęć. Dlatego nastroimy regulator wg kryterium nr 2 a ono rzecze tako:
Dla regulatora PID i obiektu z Rys. 31-3 następujące nastawy regulatora PID zapewnią najkrótszy czas regulacji przy przeregulowaniu nie większym niż 20%.
Rys. 31-4
Wzory i obliczone nastawy dla obiektu o parametrach z Rys. 31-3.
Są pewne ograniczenia. Mianowicie opóźnienie To musi być w zakresie 0.15T…0.6T. Czyli np. czysta inercja albo czyste opóźnienie odpada. Warunek jest oczywiście spełniony bo 2.1 sek<4 sek<8.4 sek. Sprawdźmy więc jak wygląda odpowiedź na skok z w/w nastawami.
Rozdział 31.2.4 Etap 3- Sprawdzenie odpowiedzi na skok układu zamkniętego
Wywołaj PID/17_dobór_nastaw_PID/02_Kp4.2_I8_D1.6.zcos
Rys. 31-5
Układ zamknięty z nastawami optymalnymi z Rys. 27-4.
Wciśnij „start”
Rys. 31-6
Co by tu powiedzieć? Po pierwsze przeregulowanie to 50% a nie 20%. Po drugie czas regulacji to ponad 1 min. Porządny regulator powinien dać odpowiedź szybszą niż w układ otwarty (bez regulatora) na Rys. 31-2.
To ja to sobie tak tłumaczę. Nastawy wyznaczone przez kryterium Zieglera-Nicholsa to tylko pierwsze przybliżenie. To tak jakby ktoś zabrał nas helikopterem w okolice Giewontu i powiedział. Teraz szukaj szczytu-optimum sam. Co prawda nie wysadził nas przy samym krzyżu na Giewoncie, ale znacznie ułatwił robotę. Zawęził się obszar poszukiwań. Dalej będziemy więc szukać nastaw optymalnych metodą prób i błędów z punktu startowego Kp=4.2 Ti=8 sek Td=1.6 sek. Będzie mniejsza liczba kombinacji Kp, Ti, Td do zbadania, niż gdybyśmy zaczynali z przypadkowego (czytaj „dalszego”) punktu startowego.
Rozdział 31.2.5 Etap 4-Próby znalezienia lepszych nastaw metodą ręczną
Wykonałem kilka prób z różnymi nastawami. W końcu znalazłem takie nastawy regulatora
Wywołaj PID/17_dobór_nastaw_PID/03_Kp3_I7_D3.zcos
Rys. 31-7
Kp=3 Ti=8 sek Td=1.6 sek
Rys. 31-8
Jest poprawa? Jest. Nawet przeregulowanie jest w okolicach 20%. Ale dlaczego pierwszy strzał (Rys. 31-6) nie był trafiony? Próbowałem nieporadnie bronić Zieglera. Że to tylko pierwsze przybliżenie do optymalnych nastaw, że… itd. Nie wiem.
Rozdział 31.3 – Metoda cyklu granicznego
Rozdział 31.3.1 Wstęp
W rozdziale 31. 2 wyznaczyliśmy optymalne nastawy Kp, Ti i Td na podstawie modelu matematycznego opisującego obiekt. Nie ważne, że był to model uproszczony-transmitancja zastępcza Gp(s) typu inercja z opóźnieniem z Rys. 31-3. Nie ważne, że powstał na podstawie doświadczenia, a nie rozważań teoretycznych. Ważne jest to, że korzystaliśmy z modelu obiektu. To była pierwsza szkoła strojenia regulatorów.
Metoda cyklu granicznego jest przykładem drugiej szkoły. Tu nie jest potrzebny model model. Wystarczy tylko zrobić „coś” obiektowi i obserwować jego zachowanie. Tym „cosiem” może być np. stopniowe zwiększanie wzmocnienia w układzie zamkniętym do momentu, w którym stanie się niestabilny.Innymi słowy – stanie się generatorem.
Rozdział 31.3.2 Próba nr 1
Wywołaj PID/17_dobór_nastaw_PID/04_Kp3_wzbudzenie.zcos
Rys. 31-9
Czym różni się schemat od Rys. 31-7?
1-Poprzez zera zostały wyłączone składowe – całkująca i różniczkująca. Czyli z regulatora PID zrobił się regulator P
2-Wymuszenie x(t) w 3 sekundzie jest teraz impulsem Diraca (właściwie pseudoimpulsem)
Przypominam, że układy niestabilne mogą się „nie ruszać”. Tak jak zaostrzony ołówek postawiony pionowo na stole. Może tak sobie stać do końca świata, chociaż jest to układ ewidentnie niestabilny. Dlatego żeby zaobserwować niestabilność, dajemy na wejście impuls Diraca.
Wciśnij „start”
Rys. 31-10
Pukneliśmy obiekt w nos młotkiem Diraca. Zabujało i układ wrócił ze stanu początkowego 0 do stanu końcowego (równowagi) 0.
Czyli dla K=3 system jest stabilny.
Zwiększmy wzmocnienie na Kp=5.
Rozdział 31.3.3 Próba nr 3
Wywołaj PID/17_dobór_nastaw_PID/05_Kp5_wzbudzenie.zcos
Otrzymasz taki sam schemat, tylko Kp=5
Wciśnij „start”
Rys. 31-11
Dłużej się buja. Wniosek-Zbliżyliśmy się do niestabilności ale układ nadal jest stabilny.
Zwiększmy wzmocnienie na Kp=6.27
Rozdział 31.3.4 Próba nr 4
Wywołaj PID/17_dobór_nastaw_PID/06_Kp6.27_wzbudzenie_graniczne.zcos
Otrzymasz taki sam schemat, tylko Kp=6.27
Wciśnij „start”
Rys. 31-12
O!!! buja się bez końca ze stałą amplitudą. Wniosek-jesteśmy na granicy stabilności. Co za przypadek, że wybraliśmy akurat Kp=6.27!
Okres oscylacji Tosc=16.3 sek.
Z tego doświadczenia zapamiętajmy 2 wartości:
– wzmocnienie Kp=6.27 które nazwiemy wzmocnieniem krytycznym Kkr=6.27
– okres oscylacji Tosc=16.3
Obydwa parametry będą potrzebne do wyznaczenia nastaw optymalnych Kp Ti i Td w p. 31.3.8
Dalsze zwiększanie wzmocnienia nie jest nam już potrzebne. Ale zróbmy to dla zwykłej ciekawości.
Zwiększmy wzmocnienie na Kp=6.5
Rozdział 31.3.5 Próba nr 5
Wywołaj PID/17_dobór_nastaw_PID/07_Kp6.5_wzbudzenie.zcos
Otrzymasz taki sam schemat, tylko Kp=6.5
Wciśnij „start”
Rys. 31-13
Tego się chyba spodziewaliśmy. Amplituda powoli rośnie do nieskończoności. Gdyby to był układ rzeczywisty, a nie liniowy idealny, to amplituda rosłaby do +/- nasycenia
Rozdział 31.3.6 Wyznaczenie parametrów optymalnych
Rys. 31-14
Parametry optymalne Kp, Ti i Td obliczamy wg. bardzo prostych wzorów–>Rys. 31-14a. Powinny zapewnić najkrótszy czas regulacji przy przeregulowaniu nie większym niż 30%.
Gdy do tych wzorów podstawimy Kkr=6.27 i Tosc=16.3 sek z p.31.3.5 Próba nr 4 to otrzymamy nastawy optymalne –> Rys. 31-14b. Proponuję sprawdzić kalkulatorem.
Rozdział 31.3.7 Odpowiedź na skok przy nastawach optymalnych
Wywołaj PID/17_dobór_nastaw_PID/8_Kp3.76_I8.15_2.04.zcos
Rys. 31-15
Nastawy są z Rys. 31-14b czyli są optymalne wg metody cyklu granicznego
Wciśnij „start”
Rys. 31-16
Nie rzuca na kolana, ale w porównaniu do poprzedniej Metody odpowiedzi skokowej z p.31.2 odpowiedź jest bliższa optymalnej. Oczywiście nie znaczy to, że metoda jest lepsza. Może przy innym obiekcie jest inaczej? Nie zajmę tu stanowiska.
Możemy próbować poprawić parametry metodą ręczną, tak żeby otrzymać odpowiedź nie gorszą niż z Rys. 31-8.
