440 likes | 1.14k Views
AUTOMATYKA i ROBOTYKA (wykład 4). Wykładowca : dr inż. Iwona Oprzędkiewicz Nazwa wydzia łu: WIMiR Nazwa katedry: Katedra Automatyzacji Procesó w AGH. Podstawowe człony dynamiczne obiekt bezinercyjny. Przykład fizyczny. Schemat równoważni:. x(t). y(t). a. b. i(t). R. y(t). u(t). C.
E N D
AUTOMATYKAiROBOTYKA(wykład 4) Wykładowca : dr inż. Iwona Oprzędkiewicz Nazwa wydziału:WIMiRNazwa katedry:Katedra Automatyzacji Procesów AGH
Podstawowe człony dynamiczneobiekt bezinercyjny Przykład fizyczny. Schemat równoważni: x(t) y(t) a b
i(t) R y(t) u(t) C Podstawowe człony dynamiczneobiekt inercyjny I rzędu Przykład fizyczny. Schemat dwójnika RC: Zakładamy, że sygnałem sterującym jest napięcie zasilające u(t), a sygnałem wyjściowym – spadek napięcia na kondensatorze y(t) Po przekształceniu w dziedzinie zmiennej zespolonej otrzymujemy:
Podstawowe człony dynamiczneobiekt inercyjny II rzędu Przykład fizyczny. Schemat procesu mieszania w zbiornikach: Roztwór o natężeniu objętościowym i stężeniu przechodzi przez dwa zbiorniki – mieszalniki o objętościach c1 oraz c2.
Podstawowe człony dynamiczneobiekt inercyjny II rzędu Jeżeli przyjmiemy całkowite wymieszanie, to dla stężeń 1 oraz 2 w poszczególnych zbiornikach możemy sformułować następujące równania bilansowe: Przyjmujemy, że sygnałem wyjściowym jest stężenie w drugim zbiorniku 1. Sygnałem wejściowym stężenie zadane . Po przekształceniach i transformacji otrzymanego równania otrzymamy:
k y(t) 1 2 3 4 u(t)=1(t) czas Podstawowe człony dynamiczneobiekt inercyjny n-tego rzędu Gdzie: k – współczynnik wzmocnienia T1 … Tn – stałe czasowe. Charakterystyki czasowe
Q(ω) ω=0 k P(ω) 1 2 3 4 Podstawowe człony dynamiczneobiekt inercyjny n-tego rzędu Charakterystyki amplitudowo-fazowe :
y(t) Ak T T/k t Podstawowe człony dynamiczneobiekt różniczkujący rzeczywisty Transmitancja obiektu: gdzie: T – czas różniczkowania, k – współczynnik wzmocnienia Charakterystyka czasowa:
w Q( ) w P( ) kT /T d Podstawowe człony dynamiczneobiekt różniczkujący rzeczywisty
20logM(ω) 20log(Td/T) +20dB/dekadę ω=1/T Φ(ω) ω /2 /4 Podstawowe człony dynamiczneobiekt różniczkujący rzeczywisty Charakterystyki częstotliwościowe logarytmiczne modułu i fazy.
Podstawowe człony dynamiczneobiekt różniczkujący rzeczywisty Przykład fizyczny. Schemat dwójnika RC: i(t) C R y(t) u(t)
y(t) y(t) u(t)=1(t) 1 Ti czas Podstawowe człony dynamiczneobiekt całkujący idealny Transmitancja obiektu: gdzie: Ti – czas całkowania. Charakterystyka czasowa
Q(ω) P(ω) ω=0 Podstawowe człony dynamiczneobiekt całkujący idealny Charakterystyka amplitudowo-fazowa :
20logM(ω) -20dB/dekadę Φ(ω) ω -/2 Podstawowe człony dynamiczneobiekt całkujący idealny Charakterystyki częstotliwościowe logarytmiczne modułu i fazy.
Podstawowe człony dynamiczne Logarytm modułu jest najczęściej mierzony w decybelach [dB], przy czym 1 dB jest równy 20 log M(). Charakterystyka częstotliwościowa logarytmiczna jest więc linią prostą o ujemnym nachyleniu. Nachylenie tej charakterystyki możemy łatwo obliczyć: załóżmy na początku, że rozważamy dwie wartości pulsacji, powiązane z sobą następująco: 1 = , 2 = 101. Wtedy –(20log2T-20log1T) = -20 ( log101T - log1T ) = -20( 1 + log1T - log1T ) = -20dB. Zmianę częstotliwości w stosunku 1 : 10 nazywamy dekadą. Stąd mówimy, że nachylenie charakterystyki wynosi –20 dB/dekadę.
Podstawowe człony dynamiczneobiekt całkujący idealny Przykład fizyczny. Jako przykład fizyczny obiektu całkującego rozważmy zbiornik o stałym polu przekroju równym S, z wymuszonym dopływem i odpływem. Załóżmy, że natężenie dopływu jest równe Fd. Oznaczmy gęstość cieczy w zbiorniku przez , a poziom cieczy przez h. Wtedy na podstawie bilansu masy możemy zapisać równanie stanu tego systemu: S Fd h
Podstawowe człony dynamiczneobiekt całkujący z inercją I rzędu Transmitancja obiektu: gdzie: T – stała czasowa, kv - współczynnik wzmocnienia prędkościowego . Charakterystyka czasowa
Podstawowe człony dynamiczneobiekt całkujący z inercją I rzędu Charakterystyka amplitudowo-fazowa :
Podstawowe człony dynamiczneobiekt całkujący z inercją I rzędu 20logM(ω) -20dB/dekadę -40dB/dekadę Φ(ω) ω=1/T ω -/2 - Charakterystyki częstotliwościowe logarytmiczne modułu i fazy.
Podstawowe człony dynamiczneobiekt całkujący z inercją I rzędu Przykład fizyczny. Przykładem obiektu całkującego z inercją jest silnik prądu stałego przy założeniu, że zbiornik energii pola magnetycznego ( indukcyjność uzwojeń ) jest pomijalnie mały w porównaniu ze zbiornikiem energii kinetycznej ruchu obrotowego ( wirujące masy ). Wtedy, przyjmując że prąd wzbudzenia jest stały i żadne opory ruchu nie występują otrzymujemy następujący schemat: Schemat silnika prądu stałego:
Podstawowe człony dynamiczneobiekt całkujący z inercją I rzędu Przyjmujemy, że sygnałem wejściowym jest napięcie zasilania u(t), a sygnałem wyjściowy – kąt obrotu wału (t) Zakładamy, że: Równanie silnika można zapisać w następującej postaci: Po przeprowadzeniu transformacji Laplace’a otrzymujemy: Definiując transmitancję jako: Otrzymujemy:
Podstawowe człony dynamiczneobiekt oscylacyjny Transmitancja obiektu: gdzie: k – współczynnik wzmocnienia, T0 – okres drgań własnych, - współczynnik tłumienia . Warunek wystąpienia oscylacji: <1
Podstawowe człony dynamiczneobiekt oscylacyjny y(t) 1 u(t)=1(t) Charakterystyka czasowa czas
Q(ω) P(ω) Podstawowe człony dynamiczneobiekt oscylacyjny Charakterystyka amplitudowo-fazowa : ω=0 k
20logM(ω) 20log(k) -20dB/dekadę Φ(ω) ω ω=1/T0 -/2 - Podstawowe człony dynamiczneobiekt oscylacyjny
Podstawowe człony dynamiczne -obiekt oscylacyjnyPrzykład – siłownik pneumatyczny membranowy Sygnałem wejściowym siłownika jest ciśnienie Pz podawane na membranę wejściową. Siła wywierana przez ciśnienie jest wprost proporcjonalna do ciśnienia oraz powierzchni membrany. Sygnałem wyjściowym jest przesunięcie trzpienia x. A - powierzchnia membrany, m - masa części ruchomych ( membrana i trzpień ), k - stałą sprężystości sprężyny podpierającej, R - współczynnik oporów ruchu części ruchomych.
Podstawowe człony dynamiczne -obiekt oscylacyjnyPrzykład – siłownik pneumatyczny membranowy • Transmitancję operatorową rozważanego układu wyznaczymy na podstawie bilansu sił występujących w nim: • Oznaczmy siłę pochodzącą od ciśnienia wejściowego przez Fp. Fp(t) = Apz(t) • Siła sprężystości sprężyny jest proporcjonalna do przesunięcia trzpienia Fs(t)=kx(t) • Siła oporu części ruchomych występuje tylko podczas ruchu i w rozważanym przypadku można uznać, że jest ona proporcjonalna do prędkości: FR(t)=Rv(t) • jest siła bezwładności. Jest ona opisana powszechnie znanym wzorem: Fb(t)=ma(t)
Podstawowe człony dynamiczne -obiekt oscylacyjnyPrzykład – siłownik pneumatyczny membranowy • Transmitancję operatorową rozważanego układu wyznaczymy na podstawie bilansu sił występujących w nim: • Oznaczmy siłę pochodzącą od ciśnienia wejściowego przez Fp. Fp(t) = Apz(t) • Siła sprężystości sprężyny jest proporcjonalna do przesunięcia trzpienia Fs(t)=kx(t) • Siła oporu części ruchomych występuje tylko podczas ruchu i w rozważanym przypadku można uznać, że jest ona proporcjonalna do prędkości: FR(t)=Rv(t) • jest siła bezwładności. Jest ona opisana powszechnie znanym wzorem: Fb(t)=ma(t)
Podstawowe człony dynamiczne -obiekt oscylacyjnyPrzykład – siłownik pneumatyczny membranowy Bilans sił można zapisać następująco: Fp = Fs+FR+Fb Po uwzględnieniu wcześniejszych zależności otrzymujemy: Apz(t) = kx(t) + Rv(t) + ma(t) Wiedząc, że: Otrzymujemy: Transformata Laplace’a powyższego równania, przy założeniu zerowych warunków początkowych na x oraz będzie mieć następującą postać: APz(s) = kX(s) + RsX(s) +ms2X(s) Jeżeli teraz przypomnimy, że wyjściem układu jest sygnał x, a wejściem – sygnał pz, to widzimy, że transmitancja operatorowa układu będzie mieć postać:
Podstawowe człony dynamiczneobiekt opóźniający Transmitancja obiektu: gdzie: - opóźnienie (czas martwy) obiektu. Charakterystyka czasowa: y(t) y(t) u(t)=1(t) 1 czas
Podstawowe człony dynamiczneobiekt opóźniający Charakterystyka amplitudowo-fazowa : Q(ω) 1 1 1 P(ω) 1
Podstawowe człony dynamiczneobiekt opóźniający Charakterystyki częstotliwościowe logarytmiczne modułu i fazy. 20logM(ω) M() = 1 L()=20logM()=0 Φ(ω) ω Φ(ω)=- ω
Podstawowe człony dynamiczneobiekt opóźniający Przykład fizyczny. Z elementami opóźniającymi najczęściej spotykamy się podczas opisu wszelkiego rodzaju procesów transportu, np. z użyciem przenośników taśmowych. Rozważamy układ pokazany na rys. 2.28. Materiał sypki na przenośnik jest podawany w punkcie a, a do zbiornika podawany jest w punkcie b, odległym od a o długość l. Taśmociąg jako element opóźniający.
Podstawowe człony dynamiczneobiekt opóźniający W rozważanym układzie możemy zauważyć, że jeśli prędkość przesuwu taśmy taśmociągu jest stała i równa v, wielkością wejściową w układzie jest masa materiału podawana na wejście w punkcie a, a wyjściem układu jest masa podawana do zbiornika w punkcie b to opóźnienie wnoszone przez ten element jest równe: = l/v . Jeżeli oznaczymy masę substancji podawaną w punkcie a przez ma, a masę podawaną do zbiornika w punkcie b przez mb, to zależność pomiędzy tymi masami jako funkcja czasu może tu być zapisana w uproszczeniu ( przy założeniu braku strat po drodze ) następująco: mb(t) = ma(t-)
Podstawowe człony dynamiczneobiekt I rzędu z opóźnieniem y(t) y(t) k u(t)=1(t) czas T Transmitancja obiektu: gdzie: • - opóźnienie (czas martwy) obiektu, k – wzmocnienie obiektu, T – stała czasowa obiektu. Charakterystyka czasowa
Podstawowe człony dynamiczneobiekt I rzędu z opóźnieniem Q(ω) k P(ω) Charakterystyka amplitudowo-fazowa :
Podstawowe człony dynamiczneobiekt I rzędu z opóźnieniem 20logM(ω) 20log(k) -20dB/dekadę ω=1/T Φ(ω) ω Charakterystyki częstotliwościowe logarytmiczne modułu i fazy.
Modele zastępcze obiektów dynamicznych k y(t) y(t) u(t)=1(t) czas Załóżmy, że mamy eksperymentalnie wyznaczoną odpowiedź skokowa nieznanego obiektu wysokiego rzędu.
Modele zastępcze obiektów dynamicznych Budowa poprawnie działającego układu sterowania nie wymaga znajomości dokładnego modelu obiektu. W wielu sytuacjach wystarczy model przybliżony, mający postać np. transmitancji zastępczej z opóźnieniem. Model zastępczy Kupfmullera I rzędu:
Identyfikacja parametrów modelu: k y(t) u(t)=1(t) czas Metoda graficzna: y(t) ym(t) T
Model zastępczy Kupfmullera II rzędu: Inne modele zastępcze obiektów dynamicznych
Model zastępczy Strejca bez opóźnienia: Inne modele zastępcze obiektów dynamicznych
Model zastępczy Strejca z opóźnieniem: Inne modele zastępcze obiektów dynamicznych