A4 Modele układów dynamicznych


Modele układów dynamicznych

Przystępując do budowy modelu matematycznego danego procesu (układu), czyli do napisania równań matematycznych opisujących ten proces, należy przede wszystkim sprecyzować wszystkie założenia upraszczające, określić za­kresy zmian poszczególnych parametrów i wielkości oraz dokonać wyboru zmien­nych opisujących ten proces. Należy podkreślić, że wybór tych zmiennych nie jest jednoznaczny, istnieje wiele równoważnych sposobów wyboru zmiennych opisujących rozpatrywany proces. Następnie korzystając z odpowiednich praw fizyki (takich jak podstawowe prawa mechaniki, prawa obwodów elektrycznych, równania bilansu masy i energii), chemii, biologii, ekonomii itd., piszemy rów­nania wiążące zmiany w czasie procesów, zachodzące pod wpływem zewnętrz­nych oddziaływań lub nagromadzonej w układzie energii (niezerowych warun­ków początkowych lub brzegowych). W przypadku układów złożonych z ele­mentów o różnej naturze fizycznej, np. układów elektromechanicznych, przy wypisywaniu równań dynamiki wygodnie jest skorzystać z zasady Hamiltona, zwanej również zasadą ekstremum działania. Sposoby układania równań dy­namiki opisujących układy o parametrach skupionych i rozłożonych objaśnimy na prostych przykładach tych układów o różnej naturze fizycznej występują­cych w technice (elektrycznej, mechanicznej, elektromechanicznej itp.), biologii i ekonomii.

Układy elektryczne

W układach (obwodach) elektrycznych nieliniowych za zmienne stanu wygodnie jest przyjąć ładunki na kondensatorach i strumienie skojarzone cewek, a w li­niowych — napięcia na kondensatorach i prądy w cewkach.

0x01 graphic

Rysunek 1. Obwód elektryczny: schemat elektryczny i schemat blokowy

Metodą bezpośrednią układania równań dynamiki układów elektrycznych jest metoda oparta na pierwszym i drugim prawie Kirchhoffa. Metoda ta polega na napisaniu na podstawie praw Kirchhoffa odpowiednich równań i przekształ­ceniu ich na postać równań stanu. Istotę tej metody objaśnimy na przykładzie obwodu elektrycznego, którego schemat podano na rysunku 1. Za zmienne stanu przyjmujemy napięcie UC na kondensatorze i prąd Ul w cewce, a za odpowiedź ~ napięcie U na oporniku o rezystancji R3. Wymuszeniami w tym obwodzie są napięcia źródłowe e1 i e2. Na podstawie pierwszego i drugiego prawa Kirchhoffa dla obwodu tego możemy napisać równania

0x01 graphic

= 0

= 0

= 0

(1)

(2)

(3)

Biorąc pod uwagę, że

0x01 graphic

(4)

oraz korzystając z równania (1), rugujemy z równań (2) i (3) prądy iC oraz i3 Otrzymamy wówczas

0x01 graphic

= 0

(5)

= 0

(6)

Z równania (6) otrzymujemy

0x01 graphic
(7)

Podstawiając zależność (7) do równania (5) otrzymujemy

0x01 graphic
(8)

Zapisując równania (7), (8) w postaci jednego równania macierzowego otrzymujemy

0x01 graphic
(9)

Z porównania równania (9) z równaniem 0x01 graphic
wynika, że w tym przypadku 0x01 graphic
, macierze A i B są równe

0x01 graphic

0x01 graphic

Aby wyznaczyć równanie wyjścia y=Cx+Du oraz elementy macierzy C i D dla tego obwodu, na podstawie drugiego prawa Kirchohoffa, piszemy równanie

0x01 graphic
(10)

Podstawiając zależność (7) do równania (10), otrzymamy następujące równanie wyjścia

0x01 graphic

W tym przypadku macierze C i D są równe

0x01 graphic

W wielu przypadkach metodą wygodniejszą i szybszą od omówionej wyżej me­tody bezpośredniej, jest metoda włączania idealnych źródeł napięcia i prądu.

Układy mechaniczne

W układach mechanicznych za zmienne stanu wygodnie jest przyjąć współrzęd­ne określające położenia poszczególnych ciał oraz prędkości tych ciał. Metodą bezpośrednią układania równań dynamiki układów mechanicznych jest metoda oparta na prawach Newtona (lub na zasadzie d'Alamberta).

0x01 graphic

Rysunek 2. Układ mechaniczny złożony z dwóch ciał o masach m1 oraz m2 i dwóch sprężyn

Metoda ta polega na napisaniu na podstawie praw Newtona odpowiednich równań i przekształceniu tych równań w postać równań zmiennych stanu. Istotę tej metody obrazuje prosty przykład układu mechanicznego, złożo­nego z dwóch ciał o masach odpowiednio m1 i m2 i dwóch sprężyn o współ­czynnikach sprężystości k1 i k2 (rys. 2). Zakładamy, że na ciało o masie m2 działa siła zewnętrzna f(t), a opory tarcia ciał są proporcjonalne do prędkości, przy czym r1 i r2 są współczynnikami tarcia odpowiednio ciała o masie m1 i m2. Wymuszeniem w tym przypadku jest siła zewnętrzna f(t), a za zmienne stanu przyjmujemy:

x1 - współrzędną określającą położenie ciała o masie m1

x2 - prędkość ciała o masie m1 (x2=0x01 graphic
);

x3 - współrzędną określającą położenie ciała o masie m2;

x4 - prędkość ciała o masie m2 (x4 = 0x01 graphic
);

Niech wektorem odpowiedzi y będzie

0x01 graphic

Biorąc pod uwagę siły działające na poszczególne ciała, na podstawie drugiego prawa Newtona (lub zasady d'Alamberta) można zapisać równaniami:

0x01 graphic

= 0

= f(t)

(11)

(12)

Uwzględniając, że 0x01 graphic
=x2 i 0x01 graphic
=x4, oraz przekształcając równania (11), (12) w równania zmiennych stanu, otrzymamy

0x01 graphic

W tym przypadku macierze A i B są równe

0x01 graphic

Równanie wyjścia ma postać

0x01 graphic

a macierze C i D są równe

0x01 graphic

Układy elektromechaniczne

Układając równania dynamiki układów złożonych, wygodnie jest skorzystać z zasady Hamiltona. Zgodnie z tą zasadą układ holonomiczny i zachowawczy (układ, w którym nie występuje rozpraszanie energii oraz energia nie jest dostarczana do układu) jest określony funkcją skalarną L, zwaną funkcją Lagrange'a, współrzędnych uogólnionych q1, q2, …, qn, prędkości uogólnionych 0x01 graphic
i czasu t

0x01 graphic

Funkcja ta jest różnicą energii kinetycznej T i energii potencjalnej U układu

L= T - U

W myśl zasady Hamiltona ruchu układu (przejście z jednego położenia w drugie) przebiega tak, aby funkcjonał

0x01 graphic
(13)

przyjmował wartość ekstremalną, zwykle minimalną. Funkcjonał (13) osiąga ekstremum wtedy, gdy są spełnione równania Lagrange'a w postaci

0x01 graphic
(14)

lub

0x01 graphic

przy czym q jest wektorem o składowych q1,q2,…,qn, a T oznacza transpozycje wektora.

W przypadku ogólnym, gdy do układu jest dostarczona energia i występuje rozpraszanie energii, układ jest określany dwiema funkcjami skalarnymi: funkcją Lagrange'a i funkcją Rayleigh F, będącą różnicą energii rozproszonej w układzie Fr i energii Fd dostarczonej do układu.

F = Fr - Fd

W tym przypadku równania Lagrange'a przyjmują postać

0x01 graphic

Sposób korzystania z równań Lagrange'a przy układaniu równań dynamiki układu obrazuje przykład przetwornika elektromechanicznego, przed­stawionego schematycznie na rysunku 3. Na rdzeń ruchomy o masie m i przekroju A, osadzony w rdzeniu nieruchomym, działa siła przyciągania elektromagne­tycznego, zależna od natężenia prądu z, oraz siła sprężyny o współczynniku sprężystości k. Przyjmujemy następujące założenia upraszczające:

Obliczamy energię kinetyczną i potencjalną dla obwodu elektrycznego i me­chanicznego. Energia kinetyczna obwodu elektrycznego

0x01 graphic

przy czym indukcyjność uzwojenia przetwornika jest określana zależnością

0x01 graphic

gdzie:

D - stała zależna od przekroju rdzenia, liczby zwojów oraz przenikalności magnetycznej powietrza;

d - grubość przekładki antymagnetycznej,

x = x(t) -odległość rdzenia ruchomego od rdzenia nieruchomego;

i = 0x01 graphic
- natężenie prądu w uzwojeniu przetwornika;

q - ładunek elektryczny.

0x01 graphic

Rysunek 3. Schemat przetwor­nika elektromechanicznego

Obwód elektryczny nie zawiera elementów pojemnościowych, wobec tego energia potencjalna równa się zeru

Ue= 0

Energia kinetyczna obwodu mechanicznego

0x01 graphic

a energia potencjalna

0x01 graphic

przy czym x0 jest odległością rdzenia ruchomego od nieruchomego w położeniu naturalnym sprężyny (siła sprężyny równa zeru). Funkcja Lagrange'a w tym przypadku jest równa

0x01 graphic

Aby z kolei wyznaczyć funkcje Rayleigh, obliczamy różnice energii rozproszonej i dostarczonej do układu dla obwodu elektrycznego i mechanicznego. Dla obwodu elektrycznego

0x01 graphic

przy czym R jest rezystancją uzwojenia, a u=u(t) napięciem doprowadzonym do uzwojenia. Z kolei dla obwodu mechanicznego mamy

0x01 graphic

przy czym Rm jest oporem mechanicznym. Funkcja Rayleigh w tym przypadku równa się więc

0x01 graphic

Za składowe wektora q przyjmujemy ładunek elektryczny q oraz odległość rdzenia ruchomego x

0x01 graphic

Na podstawie (14) dla i=1 oraz i=2 otrzymujemy

0x01 graphic
(15)

0x01 graphic
(16)

Przyjmując za zmienne stanu x1=q; x2=0x01 graphic
=i, x3=x, x4=0x01 graphic
, równania (15), (16) możemy napisać w postaci

0x01 graphic

W tym przypadku równania stanu układu są równaniami nieliniowymi. Przyjmując za składowe wektora odpowiedzi y natężenie prądu w uzwojeniu przetwornika i oraz odległość x rdzenia ruchomego od nieruchomego, otrzymamy równanie wyjścia w postaci

0x01 graphic

Powyższe przykłady podane dla układów elektromechanicznych, można uwzględnić również jako układy elektropneumatyczne, elektrohydrauliczne, elektroakustyczne itp.

Charakterystyki czasowe członów dynamicznych

Teoria automatycznego sterowania klasyfikuje układy sterowania pod względem właściwości dynamicznych, opisanych równaniami różniczkowymi. Te same równania mogą opisywać układy o różnej strukturze fizycznej (np. układ masa - sprężyna z tłumieniem drgań i obwód elektryczny RLC) - mówi się wówczas o analogii pomiędzy tymi układami. Analogie pozwalają na budowę i badanie modeli układów zamiast samych układów.

Układy opisane liniowymi równaniami różniczkowymi o stałych współczynnikach nazywają się układami liniowymi stacjonarnymi. Jeżeli współczynniki te zmieniają się w czasie, lecz nie są zależne od wielkości wejściowych ani wyjściowych układ nazywany jest niestacjonarnym. Układy opisane równaniami nieliniowymi noszą nazwę układów nieliniowych.

Dla układu liniowego, stacjonarnego i jednowymiarowego, tj. o jednym wejściu i jednym wyjściu zależność pomiędzy sygnałem wyjściowym y(t) (odpowiedzią układu) i wejściowym x(t) (wymuszeniem) określona jest ogólnym równaniem:

0x01 graphic
(17)

gdzie:

0x01 graphic
- stałe współczynniki, zależne od struktury i od wartości parametrów układu;

0x01 graphic
- stałe współczynniki, zależne od źródła sygnału wejściowego oraz od wartości parametrów układu i jego struktury.

Rząd n najwyższej pochodnej sygnału wyjściowego występującej w równaniu nazywamy rzędem układu.

Poddając obie strony równania różniczkowego (17) przekształceniu Laplace'a dla zerowych warunków początkowych dostaniemy:

0x01 graphic
(1.2)

gdzie 0x01 graphic
(1.3)

0x01 graphic
(1.4)

Stosunek transformaty Laplace'a sygnału wyjściowego Y(s) układu do transformaty Laplace'a sygnału wejściowego X(s), przy zerowych warunkach początkowych nazywamy transmitancją operatorową układu.

0x01 graphic
(1.5)

Transmitancja operatorowa jest funkcją wymierną zmiennej zespolonej s. Ma ona postać ilorazu dwóch wielomianów stopnia m oraz n, przy czym dla układów realizowalnych fizycznie zawsze stopień wielomianu licznika m jest niższy lub co najwyżej równy stopniowi wielomianu mianownika n. Transmitancja operatorowa układu nie zależy od transformat wielkości wejściowej i wyjściowej. Dla danego układu jest ona wielkością stałą, zależną jedynie od natury fizycznej układu, a więc od równania różniczkowego i parametrów układu (współczynniki wielomianów N(s) i M(s) są przeważnie prostymi funkcjami parametrów - pojemności, indukcyjności, rezystancji, masy itp.). Można zatem powiedzieć, że transmitancja operatorowa określa właściwości dynamiczne układu. Znając transmitancję układu można wyznaczyć odpowiedź układu na dowolny sygnał wejściowy x(t):

0x01 graphic
(1.6)

gdzie L-1 - operator odwrotnego przekształcenia Laplace'a.

Najważniejsze charakterystyki czasowe

Charakterystyką czasową układu nazywamy przebieg w czasie odpowiedzi układu na określony sygnał wejściowy, podany na wejście układu będącego w stanie równowagi.

Stosowanie tych samych sygnałów wejściowych do badania różnych układów pozwala na porównanie właściwości dynamicznych tych układów. Do opisywania i porównywania własności dynamicznych układów oprócz charakterystyk czasowych stosuje się także charakterystyki częstotliwościowe, będące tematem innego ćwiczenia.

W zależności od rodzaju zastosowanego sygnału wejściowego wśród charakterystyk czasowych można rozróżnić następujące:

Charakterystyka skokowa jest to odpowiedź y(t)=h(t) układu, na którego wejście doprowadzony został sygnał skokowy x(t) (rys.1.1.a) opisany równaniem:

0x01 graphic
(1.7)

gdzie funkcja skoku jednostkowego: 0x01 graphic
(1.8)

Transformata wymuszenia skokowego ma postać 0x01 graphic
(1.9)

więc odpowiedź skokowa członu 0x01 graphic
(1.10)

Charakterystyka impulsowa układu jest to odpowiedź y(t)=k(t) układu, na którego wejście doprowadzony został sygnał w postaci impulsu Diraca x(t)=δ(t) (impuls o jednostkowej energii, nieskończonej amplitudzie i nieskończenie krótkim czasie trwania):

0x01 graphic
(1.11)

Ponieważ 0x01 graphic
(1.12)

więc odpowiedź impulsowa członu: 0x01 graphic
(1.13)

Z zależności (1.13) wynika, że charakterystyka impulsowa układu, zwana także funkcją wagi, jest odwrotną transformatą Laplace'a transmitancji układu. Impuls Diraca przedstawia rys.1.1.b.

Charakterystyka liniowo-czasowa jest to odpowiedź y(t)=v(t) układu, na którego wejście doprowadzony został sygnał x(t) liniowo zależny od czasu (rys.1.1.c):

0x01 graphic
(1.14)

Ponieważ 0x01 graphic
(1.15)

więc charakterystyka liniowo-czasowa członu:

0x01 graphic
(1.16)

0x01 graphic
Rys. 1.1 Sygnał x(t) podawany na wejście układu w celu uzyskania charakterystyki:

a) skokowej b) impulsowej c) liniowo-czasowej

Charakterystyki czasowe członów podstawowych

Członem układu automatyki nazywamy urządzenie lub układ o wyodrębnionym wejściu i wyjściu będący częścią składową tego układu. Okazuje się, że istnieje ograniczona ilość liniowych członów podstawowych, a wszystkie inne układy liniowe można przedstawić jako ich połączenie; schemat układu przedstawiający te połączenia nazywa się schematem strukturalnym (blokowym).

Poniżej podano transmitancje oraz charakterystyki skokowe wszystkich członów podstawowych, oraz charakterystyki impulsowe i liniowo-czasowe dla niektórych członów.

Człon bezinercyjny (proporcjonalny) P

Transmitancja członu ma postać 0x01 graphic
(1.17)

gdzie k - współczynnik wzmocnienia, określony jako stosunek odpowiedzi do wymuszenia.

W członie bezinercyjnym w każdej chwili czasu sygnał wyjściowy jest proporcjonalny do sygnału wejściowego.

Odpowiednie charakterystyki czasowe dane są wzorami:

- skokowa 0x01 graphic
, 0x01 graphic
(1.18)

- impulsowa 0x01 graphic
, 0x01 graphic
(1.19)

- liniowo-czasowa 0x01 graphic
, 0x01 graphic
(1.20)

0x01 graphic

Rys.1.2. Charakterystyki czasowe członu proporcjonalnego

a) skokowa b) impulsowa c) liniowo-czasowa


Przykładem realizacji członu proporcjonalnego jest rezystancyjny dzielnik napięcia z rys.1.3.

Transmitancja czwórnika:

0x01 graphic
(1.21)

0x01 graphic

Rys.1.3. Czwórnik proporcjonalny


Człon inercyjny pierwszego rzędu

0x01 graphic
Transmitancja członu: 0x01 graphic
(1.22)

gdzie T - stała czasowa

Odpowiedź czasowa członu na skutek pewnej bezwładności (inercji) charakteryzuje się występowaniem stanu przejściowego, po zaniknięciu którego sygnał wyjściowy staje się proporcjonalny do sygnału wejściowego (ze współczynnikiem proporcjonalności k).

Dla odpowiedzi skokowej członu mamy:

0x01 graphic
, 0x01 graphic
(1.23)

Stała czasowa T charakteryzuje prędkość zmian przebiegu przejściowego. Jest to czas, po upływie którego odpowiedź skokowa osiąga wartość (1-1/e)⋅k⋅a=0.632⋅k⋅a. Interpretację geometryczną stałej czasowej przedstawia rys.1.4

0x01 graphic

Rys.1.4. Charakterystyka skokowa członu inercyjnego I-go rzędu

Pozostałe charakterystyki czasowe:

- impulsowa 0x01 graphic
, 0x01 graphic
(1.24)

- liniowo-czasowa 0x01 graphic
, 0x01 graphic
(1.25)

0x01 graphic

Rys.1.5. Charakterystyki czasowe członu inercyjnego I-go rzędu

a) impulsowa b) liniowo-czasowa


Transmitancja czwórnika

0x01 graphic
(1.26)

(k=1, T=RC)

0x01 graphic

Rys.1.6. Czwórnik RC inercyjny.


Człon całkujący idealny I

Transmitancja członu: 0x01 graphic
(1.27)

gdzie Ti - czas całkowania

W członie całkującym idealnym sygnał wyjściowy jest proporcjonalny do całki sygnału wejściowego.

Odpowiedź skokowa ma postać: 0x01 graphic
, 0x01 graphic
(1.28)

Jeżeli na wejściu członu całkującego idealnego pojawi się sygnał stały to sygnał wyjściowy będzie narastał w funkcji czasu liniowo. Współczynnik k reprezentuje stosunek pochodnej względem czasu (prędkości) odpowiedzi do wartości wymuszenia, stąd też nazywany jest wzmocnieniem prędkościowym.

Pozostałe charakterystyki czasowe:

- impulsowa 0x01 graphic
, 0x01 graphic
(1.29)

- liniowo-czasowa 0x01 graphic
, 0x01 graphic
(1.30)

0x01 graphic

Rys.1.7. Charakterystyki czasowe członu całkującego idealnego

a) skokowa b) impulsowa

Człon całkujący z inercją

Ściśle rzecz biorąc nie jest to człon podstawowy, gdyż można go zrealizować jako szeregowe połączenie członów całkującego idealnego i inercyjnego. Ze względu na praktyczne znaczenie jest on jednak tutaj przedstawiony. Transmitancja członu:

0x01 graphic
(1.31)

Charakterystyki czasowe:

- skokowa 0x01 graphic
, 0x01 graphic
(1.32)

- impulsowa 0x01 graphic
, 0x01 graphic
(1.33)

- liniowo-czasowa 0x01 graphic
,0x01 graphic
(1.34)

0x01 graphic

Rys.1.8. Charakterystyki czasowe członu całkującego rzeczywistego

a) skokowa b) impulsowa

Człon różniczkujący idealny D

Transmitancja układu ma postać: 0x01 graphic
(1.35)

gdzie Td - czas różniczkowania

W członie różniczkującym idealnym sygnał wyjściowy jest proporcjonalny do pochodnej sygnału wejściowego względem czasu.

Ponieważ stopień licznika transmitancji jest wyższy od stopnia mianownika człon ten jest niemożliwy do zrealizowania w praktyce i może być modelowany jedynie w przybliżeniu.

Charakterystyki czasowe: skokowa i liniowo-czasowa są postaci:

0x01 graphic
, 0x01 graphic
(1.36)

0x01 graphic
, 0x01 graphic
(1.37)

0x01 graphic

Rys.1.9. Charakterystyki czasowe członu różniczkującego idealnego

a) skokowa b) liniowo-czasowa

2.2.6 Człon różniczkujący z inercją

Człon różniczkujący rzeczywisty jest układem złożonym z szeregowo połączonych członów: inercyjnego i różniczkującego idealnego. Ma on duże znaczenie praktyczne, gdyż każdy fizycznie realizowalny człon różniczkujący posiada pewną inercję.

Transmitancja członu: 0x01 graphic
(1.38)

Charakterystyki czasowe:

- skokowa 0x01 graphic
, 0x01 graphic
(1.39)

- liniowo-czasowa 0x01 graphic
, 0x01 graphic
(1.40)

0x01 graphic

Rys.1.10. Charakterystyki czasowe członu różniczkującego rzeczywistego

a) skokowa b) liniowo-czasowa


Przykładem realizacji członu różniczkującego rzeczywistego jest czwórnik z rys.1.11.

Transmitancja czwórnika:

0x01 graphic
(1.41)

(k=RC, T=RC)

0x01 graphic

Rys.1.11. Czwórnik RC różniczkujący rzeczywisty


Człon opóźniający (opóźnienie transportowe)

Transmitancja: 0x01 graphic
(1.42)

Sygnał na wyjściu członu opóźniającego pojawia się nie w chwili doprowadzenia sygnału wejściowego, lecz po upływie czasu oznaczonego przez τ0.

Charakterystyki czasowe:

- skokowa 0x01 graphic
0x01 graphic
(1.43)

- impulsowa 0x01 graphic
0x01 graphic
(1.44)

- liniowo-czasowa 0x01 graphic
0x01 graphic
(1.45)

0x01 graphic

Rys.1.12. Charakterystyki czasowe członu opóźniającego

  1. skokowa b) impulsowa

Człon oscylacyjny drugiego rzędu

Transmitancja członu jest postaci: 0x01 graphic
(1.46)

gdzie ωn - pulsacja drgań nietłumionych

ζ - współczynnik tłumienia

Odpowiedź skokowa członu: 0x01 graphic
(1.47)

0x01 graphic
(1.48)

gdzie wielkość 0x01 graphic
jest tzw. pulsacją drgań tłumionych.

Przebieg czasowy odpowiedzi skokowej członu jest przebiegiem oscylacyjnym o pulsacji ωd.

O charakterze oscylacji decyduje współczynnik tłumienia drgań ζ (w zależności od wartości ζ wykładnik potęgi funkcji wykładniczej we wzorze (1.48) jest ujemny, dodatni lub równy zeru).

Możemy wyróżnić trzy przypadki:

Dla ζ2>1 człon przestaje być oscylacyjnym i staje się członem inercyjnym drugiego rzędu (szeregowe połączenie dwóch członów inercyjnych pierwszego rzędu).

0x01 graphic

Rys.1.13.Charakterystyka skokowa członu oscylacyjnego drugiego rzędu dla różnych wartości

współczynnika tłumienia ζ.

Odpowiedź impulsowa członu:

0x01 graphic
(1.49)

0x01 graphic
(1.50)

Człony korekcyjne pierwszego rzędu

Ogólna postać transmitancji członów korekcyjnych jest następująca:

0x01 graphic
(1.51)

gdzie k - współczynnik wzmocnienia

T1, T2 - stałe czasowe

Zależnie od tego, która stała czasowa jest większa, człon korekcyjny przyspiesza lub opóźnia fazę w układzie korygowanym.

Człon opóźniający fazę

Jeżeli w wyrażeniu na transmitancję członu korekcyjnego dwie stałe czasowe zastąpi się jedną i współczynnikiem α równym stosunkowi T1/T2, to transmitancja członu opóźniającego przyjmie następującą postać:

0x01 graphic
(1.52)

Przebieg odpowiedzi skokowej członu opisany jest wyrażeniem

0x01 graphic
(1.53)

i przedstawiony jest na rys.1.14.a.

0x01 graphic

Rys.1.14. Człon korekcyjny opóźniający fazę

  1. charakterystyka skokowa b) realizacja fizyczna

Przykładem członu korekcyjnego opóźniającego fazę jest czwórnik RC z rys.1.14.b.

Dla przedstawionego układu parametry T i α określone są wzorami

0x01 graphic
, 0x01 graphic
(1.54)

Człon przyspieszający fazę (forsujący)

W przypadku, gdy w ogólnym wyrażeniu na transmitancję członu korekcyjnego stała czasowa T2 jest większa od stałej czasowej T1 tzn. α<1 człon korekcyjny przyspiesza fazę, a jego transmitancję określa się następująco:

0x01 graphic
(1.55)

Przebieg odpowiedzi skokowej członu wyrażony jest równaniem:

0x01 graphic
(1.56)

Przykładową realizację członu przyspieszającego fazę przedstawia rys.1.15.b.

Parametry czwórnika oblicza się następująco

0x01 graphic
, 0x01 graphic
(1.57)

0x01 graphic

Rys. 1.15. Człon korekcyjny przyspieszający fazę

a) charakterystyka skokowa b) realizacja fizyczna

Człon korekcyjny drugiego rzędu opóźniająco - przyspieszający fazę

W pewnych przypadkach zachodzi potrzeba stosowania korekcji zapewniającej przy niższych częstotliwościach opóźnienie, a przy wyższych przyspieszenie fazy. Można wtedy zastosować korektor, którego działanie jest analogiczne do szeregowego połączenia członu opóźniającego i członu przyspieszającego fazę. Transmitancja takiego członu jest następująca

0x01 graphic
(1.58)

Odpowiedź skokowa członu

0x01 graphic
gdzie 0x01 graphic
(1.59)

Przykładem opisanego członu jest czwórnik RC z rys.1.16.b.

Transmitancja członu ma postać:

0x01 graphic
(1.60)

Po wprowadzeniu następujących oznaczeń :

0x01 graphic
, 0x01 graphic
(1.61)

wyrażenie na transmitancję przyjmie postać ogólną (1.58).

(Wartość liczbowa stałej α jest w przybliżeniu równa 0x01 graphic
).

0x01 graphic
Rys.1.16. Człon korekcyjny opóźniająco - przyspieszający fazę

a) charakterystyka skokowa b) realizacja fizyczna

Literatura

  1. J. Mazurek: Podstawy automatyki, PWN, Warszawa 1996.

  2. T. Kaczorek: Podstawy teorii sterowania WNT 2005;

  3. W. Pełczewski: Teoria sterowania, WNT, Warszawa, 1980.

  4. S.Węgrzyn: Podstawy automatyki, PWN, Warszawa 1972.

Laboratorium Teorii Sterowania

- 4 -

Ćwiczenie 1 (CS) - Charakterystyki czasowe członów dynamicznych

Laboratorium Teorii Sterowania

4

Ćwiczenie 1 - Charakterystyki czasowe członów dynamicznych - 13 -

- 6 - Ćwiczenie 1 (CS) - Charakterystyki czasowe członów dynamicznych

- 14 -

Ćwiczenie 1 (CS) - Charakterystyki czasowe członów dynamicznych



Wyszukiwarka

Podobne podstrony:
03 ScilabControl, 2 ROK, 3ci SEMESTR, Modele ukladow dynamicznych, materialy na lab i cw
galar,modele układów dynamicznych, równanie 1 rzędu
uklady1, 2 ROK, 3ci SEMESTR, Modele ukladow dynamicznych, wyklad
Analiza Algorytmów Genetycznych jako Ukladow Dynamicznych 08 Kotowski PhD p72
23 Analogie układów dynamicznych o różnej strukturze fizyczne
Modelowanie układów dynamicznych na elektronicznej maszynie analogowej, STUDIA - Kierunek Transport,
Charakterystyki czasowe układów dynamicznych

więcej podobnych podstron