Do Projektu iAutomatyka dołączyli:

https://iautomatyka.pl/wp-content/uploads/2015/01/wirtualne-laboratorium-grafika.jpg

Kurs regulacji PID: Wpływ Nieliniowości Na Regulację – cz. 34/34

autor: MAREK GONET.

Artykuł z serii: Kurs regulacji PID - Wirtualne Laboratorium


Rozdz. 34.1 Wstęp
Najbardziej wpływową nieliniowością na przebieg regulacji jest nasycenie we wzmacniaczu mocy (inna nazwa – człon wykonawczy) za regulatorem. Dlatego tytuł rozdziału mógłby być bardziej jednoznaczny – „Wpływ nasycenia członu wykonawczego na jakość regulacji”.
Zaraz będziesz wiedział o co chodzi. Nawet jeden element nieliniowy „psuje” cały schemat . Nie możesz więc bezkrytycznie stosować poznanej dotąd teorii, w której wszystkie elementy były liniowe.
Przypomnienie o transmitancjach.
Gdy używasz transmitancji to znaczy, że siedzisz w obiektach liniowych. W obiektach tych sygnał wyjściowy y(t) w stanie ustalonym jest proporcjonalny do sygnału wejściowego y(t)

Rys. 34-1
Charakterystyka obiektu liniowego w stanie ustalonym.
A w stanie nieustalonym (inaczej przejściowym) obiektami liniowymi rządzą równania różniczkowe liniowe, w których wszystkie współczynniki przy y(t), y'(t), y”(t),… x(t), x'(t),… są stałe. np. takie:

Rys. 34-2
Równanie różniczkowe obiektu liniowego

Rys. 34-3
Przykłady równań różniczkowych, które nie są liniowe
Czerwone to pluskwy przez które równanie różniczkowe nie jest liniowe, czytaj – trudne do rozwiązania.
Wróćmy do równania różniczkowego liniowego z Rys. 34-2, któremu odpowiada transmitancja G(s).

Rys. 34-4
Transmitancja G(s) którą opisuje równanie różniczkowe z Rys. 34-2
Dla ćwiczenia skojarz parametry 1, 2 i 3 licznika i 1, 2 mianownika powyższej transmitancji G(s) z równaniem różniczkowym liniowym z Rys. 34-2. Jak masz kłopoty to spójrz na Rys. 19-8 w rozdz. 19.
Do tej pory cały czas zajmowaliśmy się obiektami liniowymi, nie licząc rozdziału o sterowaniu dwupołożeniowym . I tylko dla nich pojęcie transmitancji G(s) ma sens. Transmitancje G(s) odpowiadają równaniom różniczkowym liniowym, a te zostały już dawno rozgryzione przez matematyków. Dlatego teoria zachowania się obiektów dynamicznych liniowych jest bardzo elegancka. Z transmitancji G(s) można względnie łatwo przewidzieć jak się będzie zachowywał układ po zamknięciu pętlą sprzężenia zwrotnego. Czy będzie stabilny? Jaki będzie stan ustalony? Jak dobrać parametry regulatora PID? Jeżeli na w schemacie blokowym pojawi się chociaż 1 element nieliniowy, to (prawie na 100%) cały obiekt stanie się nieliniowy. Czyli taki, który nie da już się opisać równaniem różniczkowym liniowym.

Rozdz. 34.2 Elementy nieliniowe
Rozdz. 34.2.1 Wstęp
Przykłady elementów nieliniowych.

Rys. 34-5
Typową cechą każdego obiektu rzeczywistego (a nie idealnego!), jest to że prędzej czy później dochodzi do stanu nasycenia. Np. dla wzmacniacza nasyceniami są napięcia zasilania +/-15V. Widać to na powyższych charakterystykach. No, może oprócz Rys. 34-5d.
Typowe charakterystyki nieliniowe
Rys. 34-5a – „coraz szybciej rosnąca z nasyceniem” – Nazwa dziwna ale trafna.
Rys. 34-5b – liniowa z nasyceniem
Rys. 34-5c – „coraz wolniej rosnąca z nasyceniem”.
Rys. 34-5d – histereza. Zauważ, że wartość y(t) zależy od tego czy x(t) rośnie czy maleje. Tu akurat nie widać +/- nasycenia. Pojawiłyby się przy większych amplitudach sygnału wejściowego x(t).
Rys. 34-5e – liniowa ze strefą martwą i nasyceniem
W następnych 2 podpunktach na wejścia elementu:
-liniowego z nasyceniem (b)
-liniowego ze strefą martwą i nasyceniem (e)
podamy czarną sinusoidę x(t) i zbadamy czerwoną odpowiedź y(t).

Rozdz. 34.2.2 Badanie elementu liniowego z nasyceniem
Wywołaj PID/20_regulacja_i_nasycenia/01_nasycenie.zcos

Rys. 34-6
Element liniowy z nasyceniem i wymuszeniem sinusoidalnym x(t)
Wciśnij start

Rys. 34-7
Nasycenie występuje przy x(t)=+/-0.5.
Przed osiągnięciem nasycenia nachylenie wynosi 1 dlatego x(t)=y(t),
dla x(t)>0.5 y(t)=+0.5
dla x(t)<0.5 y(t)=-0.5
Gdyby nasycenie było +/-1 to cały czas byłoby x(t)=y(t) i widziałbyś czerwoną sinusoidę y(t).
Gdyby nasycenie było +/-0 to cały czas byłoby y(t)=0 i widziałbyś płaskie czerwone y(t)=0 i czarną sinusoidę x(t)

Rozdz. 34.2.3 Badanie elementu liniowego ze strefą martwą i nasyceniem
Wywołaj PID/20_regulacja_i_nasycenia/2_strefa_martwa

Rys. 34-8
Element liniowy ze strefą martwą i nasyceniem przy wymuszeniu sinusoidalnym x(t)
Wciśnij start

Rys. 34-9
Strefa martwa występuje występuje dla -0.5….+0.5.
Gdyby strefa martwa była zerowa (inaczej- nie ma jej), to cały czas byłoby x(t)=y(t) i widziałbyś czerwoną sinusoidę y(t).
Gdyby strefa martwa była w zakresie -1…+<1, to cały czas byłoby y(t)=0 i widziałbyś płaskie czerwone y(t)=0 i czarną sinusoidę x(t)

Rozdz. 34.3 Jak działa idealny układ automatyki?
Czyli układ bez nieliniowości. Był to najbardziej typowy eksperyment eksperyment, ale powtórka nie zaszkodzi.

Rys. 34-10
Typowy układ regulacji z regulatorem PID i obiektem dwuinercyjnym.
Może to być np. regulator temperatury z wyjściem na grzałkę 1 kW. Grzałka zanurzona jest w niedużym zbiorniku z cieczą. Załóżmy, że ciecz jest idealna tzn. nie zamarza ani nie paruje. Spodziewamy się, że w stanie ustalonym będzie zerowy uchyb e(t), czyli y(t)=x(t)=+1.
Wciśnij start

Rys. 34-11
Niech +1 to +100°C.
I słusznie się spodziewamy bo jest składowa całkująca I. Temperatura wyjściowa już po 1.5 sek doszła do +100°C! Potem nastąpiło niewielkie przeregulowanie rzędu +10 °C i po 10 sek wystąpił stan ustalony y(t)=x(t)=+100°C. Za to sygnał sterujący sPID(t) ma duże przeregulowania. Moc „dodatnia”- grzanie jest chwilami większa od +2, moc „ujemna”-chłodzenie mniejsza od -1. Dlatego temperatura y(t) bardzo szybko, już po 1.5 sek doszła do stanu ustalonego, chociaż obiekt ma 2 inercje 3 sek i 5 sek!
To jest właśnie dobra robota ujemnego sprzężenia zwrotnego, ogólniej – automatyki. Gdyby to był układ otwarty, czytaj bez sprzężenia zwrotnego, czytaj bez automatyki, to trwałoby to pół minuty a nie tylko 1.5 sek. Ty też gotując zupę stosujesz ujemne sprzężenie zwrotne przez siebie. Podkręcasz gaz na full a potem zmniejszasz, gdy temperatura osiągnie punkt wrzenia. Trochę się rozgadałem, wróćmy więc do przebiegu.
Wyraźnie widać sygnał sterujący sPID(t) obcięty przez oscyloskop na poziomie -1 i +2 czyli -100°C…+200°C
Aby zobaczyć jak wielkie są te przeregulowania sygnału sterującego sPID(t), zmieńmy zakres oscyloskopu z -1…+2 na -20…+160. Czyli na -2000°C…+16000°C! Przypominam, że w naszym kursie istnieje temperatura -2000°C, a nawet niższa.
Wywołaj PID/20_regulacja_i_nasycenia/04_PID_ideal_pelen_widok.zcos
Otrzymasz ten sam schemat co poprzednio, tylko inne są teraz nastawy oscyloskopu.
Wciśnij start

Rys. 34-12
Co prawda sygnały x(t), y(t) wielkości ok. +100°C są teraz ledwo widoczne, za to widzisz całe sPID(t). Na początku sygnał sterujący sPID(t) waha się w granicach -5…+160. Czyli przez krótki moment moc dodatnia (grzejna) „chce” podgrzać ciecz do +16000°C a za chwilę moc ujemna „chce” schłodzić do -500°C!
W stanie ustalonym grzałka sterująca ciągnie 1kW, za to w czasie piku aż 160 kW! Żaden konstruktor nie da przecież takiej wielgachnej grzały! Najwyżej 1.5 kW. Także rzadko stosuje się moc ujemną-chłodzenie.
Dlatego ograniczeniem będzie węższy zakres 0…1.5 zamiast -5…+160 . Nawet jeszcze węższy. Przecież przy wyłączonej mocy temperatura spada nie do 0°C lecz do temperatury otoczenia np. +20°C.
Tym właśnie różni się teoria od praktyki*. Systemy liniowe które do tej pory rozpatrywaliśmy, są eleganckie, łatwe (względnie, względnie!) i dają odpowiedź na wiele ogólnych pytań. Stabilność, Hurwitz itd. Natomiast systemy rzeczywiste – nieliniowe będą odbiegać od liniowych.
Jak daleko odbiegać? Na ile ich liniowe przybliżenie jest wiarygodne? O tym będzie mowa za chwilę.
* Kiedyś w zacnej firmie automatycznej widziałem plakat z jej dewizą:
Praktyk – Wszystko działa, ale nie wie dlaczego.
Teoretyk– Nic nie działa, ale wie dlaczego.
My łączymy teorię z praktyką. Nic nie działa i nie wiemy dlaczego.
Od razu ich polubiłem.

Rozdz. 34.4 Jak działa rzeczywisty układ automatyki?
Rozdz. 34.4.1 Wstęp
W życiu nie ma ideałów. Transmitancja G(s) jest tylko przybliżeniem systemu rzeczywistego. Nasuwa się pytanie. Jak daleko odbiega ono od rzeczywistości. Pierwszym przybliżeniem systemu rzeczywistego jest coś takiego.

Rys. 34-13
W następnych punktach spróbuję uzasadnić, że nieliniowość wprowadzana jest główne przez Wzmacniacz Mocy. Dlatego kolejno rozpatrzymy pod tym względem Wzmacniacz Mocy, Regulator PID i Obiekt

Rozdz. 34.4.2 Wzmacniacz Mocy z Nasyceniem
Schemat różni się od Rys. 34-10 tylko Wzmacniaczem Mocy w środku schematu. Do tej pory zakładaliśmy, że sygnał z regulatora sPID(t) zawsze ma wystarczającą moc do wysterowania obiektu. Wiadomo, że tak nie jest. Każdy wzmacniacz ma swoje nasycenia .
Nawet elektrownia w Bełchatowie nie da więcej niż 4 000 MW. To też jest gigantyczny wzmacniacz mocy. Wejściem jest pokrętło (nie ważne, że komputerowe) mocy w Centralnej Dyspozytorni, a wyjściem moc oddana do sieci.
Wzmacniacz Mocy na schemacie to nic innego jak element liniowy z nasyceniem z Rys. 34-5b. Nie musi występować w wersji elektrycznej. Może to być np. zawór z siłownikiem. Sygnał z regulatora steruje stopniem otwarcia zaworu, pośrednio więc przepływem np. dopływem pary do turbiny. Przepływ F nie będzie dowolny, tylko w granicach 0…Fmax Czyli znowu mamy jakiś tam Wzmacniacz Mocy.
Elektryczny wzmacniacz mocy lub zawór lub coś tam jeszcze, noszą wspólną nazwę- Człon Wykonawczy.
No dobrze, mamy Człon Wykonawczy, który jako element nieliniowy „psuje” cały schemat. „Psuje” czyli nie można obliczyć transmitancji całego układu zamkniętego. A co z PID-em i obiektem?

Rozdz. 34.4.3 Regulator PID
To już prawie idealny układ liniowy o transmitancji

Rys. 34-14
Na schemacie z Rys.34-13 nastawy regulatora są następujące:
Kp=10
Ti =7 sek
Td=2 sek (różniczkowanie, jak w każdym porządnym regulatorze jest rzeczywiste)
Jeszcze komentarz do „prawie idealny układ liniowy”.
Współczesne regulatory PID realizowane są przeważnie w technice mikroprocesorowej.
Czyli odpowiedź regulatora nie jest ciągła lecz kwantowana przez zegar. Jednak przy częstym kwantowaniu odpowiedź regulatora możemy traktować jako ciągłą. Na wyjściu regulatora jest przetwornik cyfrowo/analogowy. Czyli na wzmacniacz mocy wchodzi sygnał ciągły.
Drugim założeniem jest to, że wzmacniacz mocy wcześniej wejdzie w nasycenie niż PID. Dlatego z czystym sumieniem możemy traktować PID jako transmitancję Gpid(s)

Rozdz. 34.4.4 Obiekt G(s)
Obiektem może być piec, kolumna rektyfikacyjna albo rakieta. Każdy z nich ma oczywiście inne własności dynamiczne. Na Rys. 34-13 przykładowym obiektem jest człon dwuinercyjny. Czyli milcząco założyliśmy, że jest to człon liniowy o transmitancji G(s).
Dla obiektu są słabsze argumenty dotyczące liniowości niż dla PID-u. Charakterystyka statyczna większości typowych obiektów jest „coraz wolniej rosnąca z nasyceniem” –>Rys. 34-5c.
Co z tego wynika? To że mianownik, czyli nachylenie charakterystyki równe jest jedności dla zerowego punktu pracy. a potem się zmniejsza. Tak samo stałe czasowe mogą się zmieniać w pewnych granicach. Tak że parametry „prawdziwej” transmitancji mogą zależeć od punktu pracy i prawdziwa wartość jest jakby rozmyta w granicach.

Rys. 34-15
Rys. 34-15a-na podstawie obserwacji obiektu, lub na podstawie modelu matematycznego obiektu zakładamy, że jest to obiekt dwuinercyjny o wzmocnieniu k=1, stałych czasowych T1=10 sek i T2=5 sek.
Rys. 34-15b-Prawdziwe parametry mogą być,w zależności od punku pracy i błędu pomiaru, w zakresach k=0.95…1.05, T1=9.5…10.5 sek i T2 4.5…5.5 sek.
Wiele obiektów, zwłaszcza wieloinercyjnych, da się przybliżyć jako inercję z opóźnieniem.
Pierwsze pytanie – po co? Lepszy jest model dokładny niż przybliżony. Tak jak lepiej być zdrowym i bogatym niż biednym i chorym.
Po pierwsze, łatwiej wyznaczyć jej parametry K, T i To obserwując odpowiedź na skok jednostkowy.
Po drugie dla tak wyznaczonej transmitancji, są już gotowe parametry nastaw regulatora PID. Wtedy odpowiedź układu sterowanego przez PID jest optymalna pod jakimś tam kryterium. Np. najkrótszy czas regulacji przy możliwie małych przeregulowaniach.
O inercji z opóźnieniem To i wzmocnieniu K była już mowa dużo wcześniej. Dla przypomnienia pokażę odpowiedź y(t) na skok jednostkowy x(t).

Rys. 34-16
Inercja z opóźnieniem

Rozdz. 34.4.5 Wnioski
W nawiązaniu do Rys. 34-13 stwierdziliśmy że:
Regulator PID jest obiektem liniowym o transmitancji Gpid(s)–>Rys. 34-14
– Z pewną dozą wyrozumiałości, obiekt G(s) też jest obiektem liniowym. Tym bliżej liniowości im mniejszy sygnał wejściowy.
Często obiekty przybliżane są inercją z opóźnieniem.
Człon Wykonawczy jest elementem liniowym z nasyceniem, czyli jest elementem nieliniowym. Głównie on „psuje” cały układ pozbawiając go liniowości. Czyli nie da się policzyć transmitancji układu zamkniętego Gz(s) z Rys. 34-13.
– Cała do tej pory wyłożona teoria–>badanie stabilności, dobór nastaw regulatora itd…, dotyczy układów liniowych.
Biorąc pod uwagę to, że w rzeczywistych układach regulacji zawsze występuje Człon Wykonawczy możemy dojść do wniosku cała teoria sterowania to tylko piękna błyskotka matematyczna. Nie ma w niej korzyści materialnych, tylko duchowe. Na szczęście nie jest tak do końca.

Rozdz. 34.5 Porównanie układu idealnego i rzeczywistego
Rozdz. 34.5.1 Wstęp
Będziemy badać odpowiedź na skoki jednostkowe x(t)=+1, oraz zakłócenia z(t)= +0.4 i z(t)=-0.4 przy nasyceniach członu wykonawczego 0…+1.5 i 0…+5 Na zdrowy rozum spodziewamy się, że im później wystąpi nasycenie, tym bardziej idealny jest system. Także zmniejszenie amplitudy skoku x(t) powoduje ten sam efekt. A co z zakłóceniem z(t)? Okaże się, że reakcje na małe zakłócenie zakłócenie są podobne dla układów rzeczywistych i idealnych. Tym bardziej im mniejsze zakłócenia.
W następnych doświadczeniach regulator PID o nastawach Kp=10, Ti=7 sek, Td=1.5 sek będzie sterował obiektem dwuinercyjnym o parametrach K=1, T1=10 sek i T2=5 sek. Powyższe nastawy regulatora dobraliśmy wcześniej metodą prób i błędów jako optymalne dla tego obiektu.

Rozdz. 34.5.2 Nasycenie 0…+1.5, x(t)=+1, z(t)=+0.4 (grzanie)
Mamy 2 takie same zbiorniki z identyczną cieczą obok siebie. Parametry regulatorów też są identyczne. Różnica jest tylko w mocy grzałek.
Pierwszy system-idealny ma moc nieograniczoną. Czyli nie występują żadne nasycenia. Może on dać nawet moc grzania +1000 000kW i moc chłodzenia -1000 000kW.
Drugi system-rzeczywisty ma już prawdziwy człon wykonawczy tzn. grzałkę z mocą ograniczoną +1.5 kW. Dodatkowo sygnał ujemny z regulatora nie chłodzi tylko wyłącza grzałkę. Temperatura spada wtedy do otoczenia np. do 0°C. Czyli zakres regulatora to 0…+1.5 kW.
Zakłóceniem jest włączenie grzałki o mocy +0.4 kW. Wartość zadana x(t) i zakłócenie z(t)+0.4 działają równolegle na 2 systemy.
Na pewno Idealny będzie lepszy od Rzeczywistego. Ale o ile lepszy? Czy to będzie 100:0 czy tylko 5:2?
Wywołaj PID/20_regulacja_i_nasycenia/5_PID_nas0…1.5_skok+1_zak+0.4.zcos

Rys. 34-17
Na górze jest system idealny, czyli z członem wykonawczym bez nasyceń.
Na dole jest system rzeczywisty z członem wykonawczym o nasyceniach 0…+1.5.
Jak to interpretować? Maksymalna wartość regulatora odpowiada mocy grzałki +1.5 kW, która w stanie ustalonym da temperaturę cieczy +150°C. Dolne nasycenie 0 kW to wyłączona grzałka i temperatura spada wtedy do temperatury otoczenia, np. do 0ºC.
Wciśnij „start”

Rys. 34-18
Do 60 sek mamy tylko samo x(t)=1, w 60 sek wystąpi zakłócenie +0.4-grzanie.

SYSTEM IDEALNY-człon wykonawczy bez nasycenia
Sygnał wyjściowy yi(t) i sygnał sterujący sPIDi(t) i w miarę szybko doprowadził do stanu yi(t)=1 w którym uchyb ustalony jest zerowy.
Charakterystyczne są bardzo duże wahania sygnału sterującego sPID(t) na początku skoku x(t)–>Rys. 34-12. W 60 sekundzie mamy zakłócenie z(t)=+0.4–> włączenie dodatkowej grzałki o mocy 0.4 kW. Regulator PID prawidłowo ocenił sytuację i zmniejszył moc też o 0.4 kW. Był oczywiście kilkunastosekundowy stan przejściowy, ale wszystko skończyło się tak jak trzeba.

SYSTEM RZECZYWISTY człon wykonawczy z nasyceniem
Poprzednio sygnał sterujący sPIDi(t) mógł hulać w granicach -nieskończoność …+nieskończoność , obecnie człon wykonawczy daje tylko skromne sPIDr(t) w granicach 0…+1.5. Dlatego zamiast szpili 160 jest zielony płaszczak o amplitudzie +1.5 od początku skoku do około 18 sek.
Dobre i to. W tym czasie yr(t) dąży do +1.5 zamiast do 160. Dlatego yr(t) wolniej narasta. W dodatku jest duże przeregulowanie do ok. +1.3, bo hamowanie od członu różniczkującego D jest słabsze niz w idealnym. Zauważ, że potem zielony sygnał sterujący sPIDr(t) już nie wchodzi w nasycenia i dlatego yr(t) dochodzi do 1 w stanie ustalonym (uchyb zerowy), jak w każdym porządnym systemie liniowym czyli idealnym.
W 60 sekundzie wystąpiło zakłócenie z(t)=+0.4. Przebiegi yr(t) i yi(t) idealnie się teraz pokrywają, ponieważ nie zadziałało żadne nasycenie. System rzeczywisty myśli teraz, że jest idealnym.
Mimo tego, że przeregulowanie i czas regulacji są większe systemie rzeczywistym to nie są 100 razy większe, jakby wynikało z maksymalnych sygnałów sPIDi(t)=+160 i sPIDr(t)=+1.5.
Wniosek
Nieliniowość, w tym głównie nasycenie wzmacniacza mocy (członu wykonawczego), pogarsza oczywiście jakość regulacji, ale nie tak bardzo!

Uwaga której możesz nie czytać, bo już było coś na ten temat
Co robią Kp, Ti, Td?

Kp=10-składowa proporcjonalna
Szybko sprowadza sygnał wyjściowy yi(t) do wartości 0.91

Rys. 34-19
Czyli składowa Kp nie zapewnia zerowego uchybu regulacji. A na wykresie widać przecież, że w stanie ustalonym po kilkunastu sekundach jest yi(t)=1, czyli jest zerowy uchyb regulacji. Która składowa PID to zrobiła? Oczywiście, całkująca Ti.

Ti=7sek-składowa całkująca
Sprowadza uchyb do zera. Pytanie. Skoro składowa całkująca sprowadza uchyb do zera, to po kiego licha jest Kp? Odpowiedź jest prosta.
Kp szybko sprowadza yi(t) do wartości 0.91. Składowa I też od samego początku sprowadza uchyb do zera, ale robi to powoli.
W skrócie i niezbyt precyzyjnie – Kp szybko robi yi(t)=0.91, a za resztę, czyli za doprowadzenie yi(t) do 1, odpowiada właśnie Ti.

Td=1.5 sek-składowa różniczkująca
Jeszcze bardziej niż Kp przyspiesza dojście do yi(t)=1. Na samym początku skoku Kp daje kopa w sygnale sPID(t) do wartości 10 zaś za resztę czyli za 150 odpowiada składowa różniczkująca, czyli Td. Przecież toto różniczkuje skok x(t). I tak całe szczęście, że jest to różniczkowanie rzeczywiste a nie idealne. Wtedy mielibyśmy nieskończenie wielką szpilę, a nie nędzne sPID(t)=160.
Wróćmy do tego kopa od skoku x(t). Dzięki niemu yi(t) bardzo szybko rośnie. Tak jakby chciał dojść nie 1 ale do 160. Za chwilę x(t) jest płaskie–> sPID szybko spada. Ale za to działa różniczkowanie od narastającego yi(t). I to działa w kierunku przeciwnym niż początkowego kopa 160! Przez chwilę daje nawet wartości ujemne sPID(t). A co to oznacza? Że składowa Td działa teraz hamująco czyli stabilizująco. Zapobiega w ten sposób przeregulowaniom a nawet niestabilnościom

Rozdz. 34.5.3 Nasycenie 0…+1.5, x(t)=+1, z(t)=-0.4 (chłodzenie)
Schemat różni się tylko zakłóceniem. Poprzednio było z(t)=+0.4-grzanie, teraz jest z(t)=-0.4-chłodzenie
Wywołaj PID/20_regulacja_i_nasycenia/06_PID_nas0…1.5_skok+1_zak-0.4
Otrzymasz schemat różniący się od poprzedniego tylko znakiem zakłócenia.
Wciśnij „start”

Rys. 34-20
Z oczywistych powodów przebiegi do 60 sek są identyczne z Rys. 34-18. Na zakłócenie z(t)=-0.4-chłodzenie, regulator zareagował prawidłowo zwiększając moc grzania o +0.4. Sygnał sterujący nie wszedł w nasycenie, regulator „myśli”, że system jest liniowy. Dlatego sygnały sPIDi(t) i sPIDr(t) pokrywają się.

Rozdz. 34.5.4 Nasycenie 0…+5, x(t)=+1, z(t)=+0.4 (grzanie)
Ten i następny schemat różnią się od poprzednich tylko nasyceniami. Teraz będzie większe 0…+5 zamiast 0….+1.5. Większe tzn. system jest bardziej „idealny”, czyli spodziewamy się lepszej regulacji. Mamy teraz człon wykonawczy o większej mocy–>czytaj droższy.
Wywołaj PID/20_regulacja_i_nasycenia/07_PID_nas0…+5_skok+1_zak+0.4
Otrzymasz schemat, który nie różni się wizualne od poprzednich.
Wciśnij „start”

Rys. 34-21
Na Rys. 34-18 widziałeś wpływ nasycenia +1.5 na sygnał sterujący rzeczywisty sPIDr(t)– (za wzmacniaczem mocy). Tu jest większe górne nasycenie – +5 i dlatego nie widzisz tego zielonego p łaszczaka sPIDr(t). W każdym bądź razie większe maksymalne sPIDr(t) spowodowało szybszą regulację i mniejsze przeregulowanie. Można powiedzieć że system zbliżył się do idealnego. Reakcja na zakłócenie nie poprawiła się, bo i tak była taka, jak dla systemu idealnego.

Rozdz. 34.5.5 Nasycenie 0…+5, x(t)=+1, z(t)=-0.4 (chłodzenie)
Schemat różni się od poprzedniego tylko zakłóceniem z(t)=-0.4 (chłodzeniem).
Wywołaj PID/20_regulacja_i_nasycenia/08_PID_nas0…+5_skok+1_zak-0.4
Otrzymasz schemat, który nie różni się wizualne od poprzednich.
Wciśnij „start”

Rys. 34-22
Różnica jest tylko w reakcji na zakłócenie z(t)=-0.4-chłodzenie.

Rozdz. 34.5.6 Wnioski
W porównaniu do systemu Idealnego system Rzeczywisty wolniej i z większymi przeregulowaniami dochodzi do stanu ustalonego. Jest to oczywiście pewną wadą, ale da się z tym żyć. Zwłaszcza, że jest zapewniony zerowy uchyb ustalony i tłumienie zakłóceń jest prawie identyczne jak w idealnym.
Weź też pod uwagę to, że system automatyki głównie tłumi zakłócenia z(t), zaś zmiany wartości zadanej x(t) są znacznie rzadsze.
Chociaż system idealny daje na początku skoku ponad 100 razy większego kopa , niż rzeczywisty (160 do 1.5), to reakcja nie jest 100 razy gorsza-czytaj wolniejsza. Dlatego (z pewną wyrozumiałością) możemy traktować systemy rzeczywiste tak jak liniowe. Ale tylko wtedy, gdy x(t) jest między nasyceniami członu wykonawczego–> w naszym przypadku x(t)=0…+1.5.
Ten warunek nie jest trudny do spełnienia. A to znaczy, że dla systemów rzeczywistych możemy korzystać z całej bogatej wiedzy dla systemów liniowych.
Wracając zaś do piłki nożnej, to relacje między systemem Idealnym a Rzeczywistym, to bardziej Niemcy-Polska niż Niemcy-Gibraltar.

Rozdz. 34.6 Jak się zachowuje system, gdy wartość zadana x(t) jest większa od nasycenia?
Rozdz. 34.6.1 Pierwsze podejście-tylko x(t), yi(t) i yr(t)
Maksymalna moc za regulatorem to +0.95 kW, która jest w stanie podgrzać ciecz do temp. +95°C. Ale chcę, żeby regulator podgrzał ciecz do +105 °C. Wiadomo, że jest to fizycznie niemożliwe! Tak jakby kierownik żądał od pracownika wykonania zadania powyżej jego kompetencji.
Wywołaj PID/20_regulacja_i_nasycenia/09_PID_nas0…+0.95_skoki_-1…2_simple

Rys. 34-23
Tym razem nie będzie to pojedynczy skok jednostkowy wartości zadanej x(t), lecz sygnał z generatora schodków. Regulator PID rzeczywisty ma na wyjściu wzmacniacza mocy ograniczenia 0…+0.95. Najpierw x(t)=+0.3 (poniżej nasycenia) +0.95, potem x(t)=1.05 (powyżej nasycenia!) a na koniec spadnie do x(t)=+0.5 (znowu poniżej nasycenia).
Spróbuj przewidzieć przebieg yr(t). Reakcja yi(t) systemu idealnego jest oczywista.
Wciśnij „start”

Rys. 34-24
okres 1. x(t)=+0.3
Wartość zadana x(t) jest mniejsza od nasycenia=+0.95. Zadziała więc sprzężenie zwrotne. Sygnały yi(t) oraz yr(t) dojdą więc do wartości x(t)=yi(t)=yr(t)=+0.3, chociaż regulator rzeczywisty (czyli z nasyceniem) zrobi to wolniej.

okres 2. x(t)=+1.05
Wartość zadana x(t) każe osiągnąć regulatorowi rzeczywistemu stan yr(t)=+1.05, czyli temperaturę na wyjściu+105°C.
Ale jak biedaczysko PID ma to zrobić, skoro człon wykonawczy ma nasycenie +0.95? Nie da się! Po prostu yr(t) będzie dążyło do +0.95. Nie zadziała ujemne sprzężenie zwrotne! Tak, jakby to był układ otwarty.
Widać to zresztą na przebiegu, gdzie dopiero pod koniec okresu 2 sygnał yr(t) osiągnął stan ustalony +0.95. Regulator idealny nie ma żadnych ograniczeń. Tu ujemne sprzężenie zwrotne szybko doprowadziło do stanu yi(t)=x(t)=1.05

okres 3. x(t)=+0.5
Wartość zadana x(t) spadła do +0.5. Błyskawicznie zareagował PID idealny dając yi(t)=x(t)=+0.5. Nie ma się czemu dziwić. PID rzeczywisty też doprowadził do stanu yr(t)=x(t)=+0.5. Zrobił to oczywiście wolniej. Ale skąd wziął się dodatkowo czas martwy, którego nie było w okresie 1?
A może by tak zbadać sygnały sterujące za regulatorem idealnym i rzeczywistym?

Rozdz. 34.6.2 Drugie podejście-dokładniejsze
Do poprzedniego schematu dodamy obserwację sygnałów sterujących z sPIDi(t) i sPIDr(t). Może pomogą odpowiedzieć na pytanie skąd się wziął czas martwy na Rys. 34-24?
Wywołaj PID/20_regulacja_i_nasycenia/10_PID_nas0…+0.95_skoki_-1…2.zcos

Rys. 34-25
Wciśnij „start”

Rys. 34-26
Sygnały systemu idealnego są oczywiste. Silniejsze i krótsze sPIDi(t) spowodowało jak zwykle szybką odpowiedź yi(t). Wyraźny jest za to w okresie 2 wpływ nasycenia wzmacniacza mocy na obcięcie sPIDr(t) regulatora rzeczywistego. Tu nie działa już sprzężenie zwrotne, i dlatego odpowiedź jest jak dla układu otwartego. Tylko skąd wziął się ten czas martwy w okresie 3?
To zaobserwujmy jeszcze sPIDir(t) z przed wzmacniacza mocy.

Rozdz. 34.6.3 Trzecie podejście-jeszcze dokładniejsze
Wywołaj PID/20_regulacja_i_nasycenia/11_PID_nas0…+0.95_skoki_-1…2_more.zcos

Rys. 34-27
Wciśnij „start”

Rys. 34-28
System Idealny
Przebieg żółtego sygnału sterującego sPIDi(t) jest oczywisty i nie wymaga komentarza.
System Rzeczywisty
Przypominam, że obserwujemy zielony sygnał sPIDir(t) z przed regulatora mocy. Literka i w nazwie sPIDir(t) sugeruje, że sygnał jest za częścią idealną rzeczywistego regulatora PID
okres 1 x(t)=+0.3
Do czasu tz=6.5 sek zielony sygnał sPIDir(t) jest większy od nasycenia +0.95. Objawi się to obcięciem sygnału za wzmacniaczem mocy niebieskiego sPIDr(t). Wtedy yr(t) dąży do +0.95 jak w układzie otwartym. Sprzężenie zwrotne pojawi po tz=6.5 sek, gdy sPIDi(t) spadnie poniżej +0.95. Potem sprzężenie zwrotne szybko sprowadzi sygnał wyjściowy yr(t) do wartości yr(t)=x(t)=0.3. Szybko ale oczywiście wolniej niż yi(t) w systemie Idealnym.

okres 2 x(t)=+0.95
Sygnał sPIDir(t) skoczył do góry. A co się z nim dzieje dalej, to nie widać bo za mały zakres oscyloskopu. Tak czy owak sygnał za wzmacniaczem mocy jest obcięty na poziomie sPIDr(t)=+0.95 i yr(t) dąży do +0.95 (inaczej +95°C ) tak jak w układzie otwartym (bez sprzężenia zwrotnego).

okres-3 x(t)=+0.5
Chyba coś świta z czasem martwym. Widać, że w tym czasie yir(t) jest większe od +0.95. Wtedy yr(t) jest w przybliżeniu obcięte na poziomie też +0.95. Gdy yir(t) stanie się mniejsze od +0.95 to sprzężenie zwrotne doprowadzi yr(t) do +0.5
Czy przyczyną jest zielony sygnał sterujący sPIDir(t)? A jak on wygląda , gdy nie jest obcięty przez oscyloskop?

Rozdz. 34.6.4 To samo tylko większy zakres oscyloskopu
Wywołaj PID/20_regulacja_i_nasycenia/12_PID_nas0…+0.95_skoki_-10…140_more.zcos
Otrzymamy wizualnie ten sam schemat, ale z innym zakresem oscyloskopu -100…+140, poprzednio -1+2. Co do systemu idealnego, to wszystko jest jasne. Dlatego interesuje nas tylko zielony sygnał sterujący sPIDir(t) systemu rzeczywistego.
Wciśnij „start”

Rys. 34-29
Ze względu na zakres oscyloskopu, sygnały x(t), sPIDr(t), yi(t), yr(t) są teraz tyci, ledwo widoczne.
Za to widzisz całe sygnały sterujące sPIDi(t) i sPIDir(t).
Żółty sygnał sterujący sPIDi(t) jest przesłonięty przez zielony sygnał sterujący sPIDir(t) i dlatego widzisz tylko zielony sygnał sterujący sPIDir(t).
okres 1 x(t)=+0.3
Opadająca zielona szpila yir(t) do czasu tz=6.5 sek. Nie widziałeś tego na Rys. 34-28.

okres 2 x(t)=+0.95
Cały czas nie działa sprzężenie zwrotne bo sPIDir(t)>0.95. Dlatego mamy klasyczną odpowiedź PID-u na dodatni skok jednostkowy w układzie otwartym!.

okres 3 x(t)=+0.5
Wartość zadana x(t) spadła z +1.05 na +0.5. Czyli mamy znów klasyczną odpowiedź PID-u na ujemny skok jednostkowy x(t). Wyraźnie widać teraz, że czas martwy to dochodzenie składowej całkującej regulatora PID do takiej wartości w której sPIDir(t) stanie się mniejsze od nasycenia +0.95. Potem sprzężenie zwrotne sprowadzi yr(t) do wartości yr(t)=x(t)=0.5–> zerowy uchyb regulacji.

Rozdz. 34.6.5 Wnioski
1. Jeżeli regulator rzeczywisty „każe” obiektowi wejść w stan który jest możliwy do spełnienia* przez wzmacniacz mocy, to regulacja działa w sposób podobny do systemu idealnego. Odpowiedź yr(t) jest oczywiście wolniejsza (gorsza) niż dla idealnego, ale w końcu yr(t)=x(t) czyli zapewniony będzie zerowy uchyb regulacji.
Uwaga* „możliwy do spełnienia”–>wzmacniacz mocy nie wejdzie w nasycenia tak jak w okresie 1 na Rys. 34-28 czyli gdy x(t)=0…+0.95
2. Jeżeli regulator rzeczywisty „każe” obiektowi wejść w stan który nie jest możliwy do spełnienia, czyli gdy x(t)>0.95 to sygnał wyjściowy yr(t) wolno dojdzie do stanu nasycenia tak jak w układzie otwartym–> okres 2 na Rys. 34-28.
3. Jeżeli regulator wychodzi ze stanu nasycenia, tak jak w okresie 3 na Rys. 34-28 to wystąpi wtedy czas martwy. Przyczyną tego jest składowa całkująca regulatora–>okres 3 na Rys.34-29.

Jak walczyć z czasem martwym?
Czas martwy To wprowadza do systemu opóźnienie które jest oczywiście zjawiskiem niekorzystnym. Bo czy chciałbyś prowadzić samochód, w którym kierownicą skręcisz w lewo, a koła zareagują dopiero po To=1 sek?
Obecnie regulatory wykonane są w technice mikroprocesorowej. Dlatego gdy regulator stwierdzi, że pomimo wzrastającego od składowej całkującej sygnału sterującego, sygnał wyjściowy nie wzrasta , to wyłącza on składową całkującą. Czyli w stanie nasycenia z regulatora PID robi się regulator PD.
Jest to przykład regulatora o zmiennej strukturze. To przełączenie może wystąpić np. od wejścia wzmacniacza mocy w stan nasycenia, a to jest łatwe wykrycia przez regulator mikroprocesorowy.


Więcej z serii: Kurs regulacji PID - Wirtualne Laboratorium


Utworzono: / Kategoria:
  • Autor: MAREK GONET
  • Wydział Elektryczny - specjalność Automatyka Przemysłowa skończyłem na Politechnice Warszawskiej w latach 1966...72.  Potem zajmowałem się projektowaniem automatyki i różnymi rzeczami. Finał kariery zawodowej to sygnalizacja pożaru, sieci strukturalne i inne wynalazki... Nie to jednak tygrysy lubią najbardziej. W 2008 roku trafiła mi się wreszcie dobra fucha - emerytura. Mam teraz dużo wolnego czasu i stąd ten kurs. Nie ma w nim dużo teorii, za to wszystko staram się wytłumaczyć tak, jak to czuję. Często przeginam w uproszczeniach np. "Automatyka to głównie dochodzenie sygnału wyjściowego y(t) do wartości zadanej x(t) i tłumienie zakłóceń z(t)" Z drugiej strony jestem pewien, że kilkaset doświadczeń z układami regulacji automatycznej na pewno nie zaszkodzi! Może nawet poczujesz radochę strojąc regulator PID. Emocje nie są oczywiście takie same jak na prawdziwym obiekcie. Gdy zawory stukają, siłowniki syczą i alarmy wyją. W dodatku z nieznanych powodów tryska kwasem siarkowym prosto na personel. W naszym laboratorium nie podlegamy już takim emocjom. Za to jest bezpiecznie i kierunek właściwy. Zachęcam do stawiania pytań. Nie wiem czy na wszystkie odpowiem. Tu liczę na pomoc. Zwłaszcza Tych którym zawory stukają, siłowniki syczą i  alarmy wyją...
  • Profil Autora
  • http://Kowieńska%209

Reklama

Newsletter

Zapisz się i jako pierwszy otrzymuj nowości!



PRZECZYTAJ RÓWNIEŻ



NAJNOWSZE PUBLIKACJE OD UŻYTKOWNIKÓW I FIRM

Reklama



POLECANE FIRMY I PRODUKTY
  • ÖLFLEX® CLASSIC 110 – elastyczny przewód sterowniczy do różnych zastosowań, w płaszczu z PVC, aprobata VDE, odporność na oleje, 300/500 V, również do YSLY lub YY CPR: informacje pod adresem www.lapppolska.pl Certyfikat zgodności VDE z...
  • PR200 to uniwersalne i łatwe w obsłudze urządzenie zaprojektowane w plastikowej obudowie do montażu na szynie DIN jako alternatywa dla PLC. Przekaźnik jest dostępny w kilku wersjach dla napięcia stałego i przemiennego. Jest wyposażony w cyf...
  • Bezpieczny, prosty w obsłudze i energooszczędny Seria serwowzmacniaczy Mitsubishi Electric MELSERVO MR-J4 wraz z kompatybilnymi modułami pozycjonującymi oraz zaawansowanymi kontrolerami motion, umożliwia konstruktorom maszyn i urządzeń oraz...
  • Seria EX-Z Czujniki z serii EX-Z to jedne z najmniejszych urządzeń tego typu na świecie. Najcieńszy model posiada grubość jedynie 3 mm co zostało osiągnięte przez zastosowanie nowych półprzewodników i dzięki temu wyeliminowanie przewodów. B...
  • Poniższy poradnik jest zbiorem schematów połączeń elektrycznych. W poradniku zapoznamy się z podstawami wprowadzenia do systemów przekaźnikowych, sekwencji przełączeń przekaźników, porównania systemów przekaźnikowych z systemami tradycyjnym...
  • Selektor napędów Panasonic umożliwia przeglądanie napędów z serii MINAS, wyszukiwanie ich w prosty sposób, a nawet porównywanie ze sobą. Dzięki wyszukiwaniu po słowach kluczowych i przy użyciu funkcji filtrowania, potrzeba zaledwie sekund a...