Charakterystyki Czasowe Czlonow Dynamicznych
Charakterystyki Czasowe Czlonow Dynamicznych
Charakterystyki Czasowe Czlonow Dynamicznych
Wydział Elektryczny
Zespól Automatyki (ZTMAiPC)
Ćwiczenie 1
CS
Charakterystyki czasowe członów dynamicznych
1. Cel ćwiczenia
• Zapoznanie się z charakterystykami czasowymi podstawowych członów dynamicznych.
• Przeprowadzenie pomiarów charakterystyk czasowych członów w postaci układów
elektrycznych RC.
• Nabycie umiejętności określania parametrów transmitancji członów na podstawie
zmierzonych charakterystyk czasowych.
2. Podstawy teoretyczne
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:
d n y(t ) dy (t ) d m x(t ) dx(t )
an ⋅ n
+ ! a1 ⋅ + a0 ⋅ y (t ) = bm ⋅ m
+ !b1 ⋅ + b0 ⋅ x (t ) (1.1)
dt dt dt dt
gdzie:
ai , i = 0! n - stałe współczynniki, zależne od struktury i od wartości parametrów układu;
b j , j = 0! m - 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 (1.1) przekształceniu Laplace’a dla zerowych
warunków początkowych dostaniemy:
M ( s) ⋅ Y ( s) = N ( s) ⋅ X ( s) (1.2)
gdzie M ( s ) = a n s n + ! a1 s + a 0 (1.3)
N ( s ) = bm s m + !b1 s + b0 (1.4)
0 dla t < 0
gdzie funkcja skoku jednostkowego: 1(t ) = (1.8)
1 dla t ≥ 0
a
Transformata wymuszenia skokowego ma postać L[x(t )] = X ( s ) = (1.9)
s
a
więc odpowiedź skokowa członu h(t ) = L−1 [H ( s )] = L−1 [G ( s ) X ( s )] = L−1 G ( s ) (1.10)
s
2.1.2 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):
0 dla t ≠ 0
∞
x(t ) = δ(t ) = , ∫ δ(t )dt = 1 (1.11)
∞ dla t = 0 −∞
Ponieważ X ( s ) = L[δ(t )] = 1 (1.12)
więc odpowiedź impulsowa członu: k (t ) = L−1 [K ( s )] = L−1 [G ( s ) X ( s )] = L−1 [G ( s )] (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.
0 t < 0
x(t ) = (1.14)
b ⋅ t t ≥ 0
b
Ponieważ X ( s ) = L[x(t )] = (1.15)
s2
więc charakterystyka liniowo-czasowa członu:
b
v(t ) = L−1 [V ( s )] = L−1 [G ( s ) X ( s )] = L−1 G ( s ) 2 (1.16)
s
arctan (b)
t t t
Rys. 1.1 Sygnał x(t) podawany na wejście układu w celu uzyskania charakterystyki:
a) skokowej b) impulsowej c) liniowo-czasowej
W członie bezinercyjnym w każdej chwili czasu sygnał wyjściowy jest proporcjonalny do sygnału
wejściowego.
k⋅a
arctan (k⋅b)
t t t
R1
k
Transmitancja członu: G( s) = (1.22)
Ts + 1
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:
t
k a −
H (s) = , h(t ) = k ⋅ a ⋅ (1 − e T
) ⋅ 1(t ) (1.23)
(Ts + 1) s
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
h(t) h(t)
k⋅a
k⋅a
0.95⋅k⋅a T1
T2
0.632⋅k⋅a T3
T1 < T2 < T3
t t
T
3T
Rys.1.4. Charakterystyka skokowa członu inercyjnego I-go rzędu
Ćwiczenie 1 (CS) – Charakterystyki czasowe członów dynamicznych
-4-
Laboratorium Teorii Sterowania
0.368 k
T arctan (k⋅b)
t t
T T
Rys.1.5. Charakterystyki czasowe członu inercyjnego I-go rzędu
a) impulsowa b) liniowo-czasowa
Transmitancja czwórnika
x(t) C y(t)
1
G( s) = (1.26)
1 + RCs
(k=1, T=RC)
Rys.1.6. Czwórnik RC inercyjny.
k 1
Transmitancja członu: G( s) = = (1.27)
s Ti s
gdzie Ti – czas całkowania
W członie całkującym idealnym sygnał wyjściowy jest proporcjonalny do całki sygnału wejściowego.
k ⋅a
Odpowiedź skokowa ma postać: H ( s) = , h(t ) = k ⋅ a ⋅ t ⋅1(t ) (1.28)
s2
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:
k
- impulsowa K (s) = , k (t ) = k ⋅ 1(t ) (1.29)
s
b k ⋅b 2
- liniowo-czasowa V ( s) = k 3 , v( t ) = ⋅ t ⋅ 1(t ) (1.30)
s 2
a) h(t) b) k(t)
arctan (k⋅a)
t t
Rys.1.7. Charakterystyki czasowe członu całkującego idealnego
a) skokowa b) impulsowa
Charakterystyki czasowe:
k ⋅a −
t
- skokowa H (s) = , h ( t ) = k ⋅ a ⋅ t − T ⋅ (1 − e T
) ⋅ 1(t ) (1.32)
s (Ts + 1)
2
t
k −
- impulsowa K (s) = , k (t ) = k ⋅ (1 − e T
) ⋅ 1(t ) (1.33)
s(Ts + 1)
k ⋅b t t2 −
t
- liniowo-czasowa V ( s ) = 3 , v(t ) = k ⋅ a ⋅ T ⋅ 1 − + 2 − e ⋅ 1(t )
2 T
(1.34)
s (Ts + 1) T 2T
a) h(t) b) k(t)
k⋅a⋅T
k
0.632⋅k
arctan ( k⋅a)
T t T t
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:
H ( s ) = ka , h(t ) = k ⋅ a ⋅ δ(t ) (1.36)
b
V ( s) = k , v(t ) = k ⋅ b ⋅ 1(t ) (1.37)
s
a) h(t) b) v(t)
k⋅b
t t
Rys.1.9. Charakterystyki czasowe członu różniczkującego idealnego
a) skokowa b) liniowo-czasowa
a) h(t) b) v(t)
k⋅a
T k⋅b
0.632⋅k⋅b
0.368 ⋅k⋅a
T
t t
T T
Rys.1.10. Charakterystyki czasowe członu różniczkującego rzeczywistego
a) skokowa b) liniowo-czasowa
Transmitancja: G ( s ) = k ⋅ e − τ0 s (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:
k ⋅ a − τ0 s
- skokowa H (s) = e h(t ) = k ⋅ a ⋅ 1(t − τ0 ) (1.43)
s
- impulsowa K ( s ) = k ⋅ e − τ0 s k (t ) = k ⋅ δ(t − τ 0 ) (1.44)
k ⋅b
- liniowo-czasowa V ( s ) = 2 e − τ0 s v(t ) = k ⋅ b ⋅ (t − τ 0 ) ⋅ 1(t − τ 0 ) (1.45)
s
a) h(t) b) k(t)
k⋅a
t t
τ0 τ0
kω 2n
Transmitancja członu jest postaci: G (s) = (1.46)
s 2 + 2ζω n s + ω 2n
gdzie ωn - pulsacja drgań nietłumionych
ζ - współczynnik tłumienia
a ⋅ k ⋅ ω 2n
Odpowiedź skokowa członu: H ( s) = (1.47)
( s 2 + 2ζω n s + ω 2n ) s
e − ζωnt 1− ζ2
h(t ) = k ⋅ a ⋅ 1 −
sin ω d ⋅ t + arctan ⋅ 1(t ) (1.48)
1 − ζ 2 ζ
Ćwiczenie 1 (CS) – Charakterystyki czasowe członów dynamicznych
-8-
Laboratorium Teorii Sterowania
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 0<ζ<1 i amplituda oscylacji maleje tzw. drgania tłumione (rys.1.13.a,b),
- dla ζ=0 występują oscylacje o stałej amplitudzie (rys.1.13.c),
- dla -1<ζ<0 amplituda oscylacji rośnie do nieskończoności (rys.1.13.d).
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).
a) b)
h(t) h(t) ζ>0
e − ζωnt hmax ζ1
k ⋅ a ⋅ (1 + ) ζ2
hmax 1 − ζ2
ζ3
k⋅a
k⋅a
t t
π \ ωd
c) d)
h(t) ζ=0 h(t) ζ<0
e− ζω nt
2⋅k⋅a k ⋅ a ⋅ (1 + )
1 − ζ2
k⋅a k⋅a
t
− ζω nt
e
k ⋅ a ⋅ (1 − )
t 1 − ζ2
π \ ωn
e − ζωn t 1− ζ2
k (t ) = k
sin ω d ⋅ t + arctan ⋅ 1(t ) (1.50)
1− ζ2 ζ
1 + T1 s
G( s) = k (1.51)
1 + T2 s
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.
a
R2
a x(t) y(t)
α C2
t
αT
Rys.1.14. Człon korekcyjny opóźniający fazę
a) 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
R1 + R2
T = R2 C2 , α= (1.54)
R2
b) 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:
1 + Ts
G( s) = α , α <1 (1.55)
1 + αTs
Przebieg odpowiedzi skokowej członu wyrażony jest równaniem:
−
t
− αtT
h(t ) = a ⋅ α ⋅ 1 − e αT + e ⋅ 1(t )
(1.56)
Przykładową realizację członu przyspieszającego fazę przedstawia rys.1.15.b.
Parametry czwórnika oblicza się następująco
R2
T = R1C1 , α= (1.57)
R1 + R2
a) b)
h(t)
R1
x(t) y(t)
a⋅α C1 R2
t
α⋅T
Rys. 1.15. Człon korekcyjny przyspieszający fazę
a) charakterystyka skokowa b) realizacja fizyczna
a) h(t) b) R1
a
R2
T1T2 T C1 y(t)
t0 = ln α 2 1 x(t)
T
h(t0) αT1 − 2 T2
α C2
t0
t
Literatura
1. J. Mazurek: Podstawy automatyki, PWN, Warszawa 1996.
2. W. Pełczewski: Teoria sterowania, WNT, Warszawa, 1980.
3. S.Węgrzyn: Podstawy automatyki, PWN, Warszawa 1972.
3. Przebieg ćwiczenia
Ćwiczenie polega na badaniu charakterystyk skokowych członów dynamicznych zrealizowanych w
formie czwórników RC. Charakterystyki poszczególnych członów zdejmowane są w układzie
pomiarowym pokazanym na rys.1.17.
C1
Generator fali
prostokątnej R1
R2
kanał 0
Uwe
Uwe Uwy kanał 1
C2
Punkty 3-7 powtórzyć dla kolejnych czwórników RC. Czwórniki należy łączyć w następującej
kolejności:
• proporcjonalny
• inercyjny I-go rzędu
• różniczkujący rzeczywisty
• opóźniający fazę
• przyspieszający fazę
• opóźniająco-przyspieszający fazę
Schematy poszczególnych członów pokazano na rys.1.18.
a) b) c)
C1
R1 R1
R2 C2 R2
d) e) f)
C1 C1
R1
R2 R2
R1 R2 R1
C2 C2
4. Opracowanie sprawozdania
Sprawozdanie powinno zawierać dla każdego badanego czwórnika:
1. Schemat połączeń i wartości elementów RC,
2. Teoretycznie wyprowadzenie transmitancji operatorowej G(s) członu. Wyrażenie parametrów
transmitancji: T, k (ewentualnie α) jako funkcji odpowiednich elementów RC i wyznaczenie ich
wartości liczbowej (patrz przykład obliczeniowy) - dla każdego z członów parametry transmitancji
oblicza się dla dwóch wybranych wcześniej wartości elementów RC,
3. Graficzne wyznaczenie z wykresów parametrów transmitancji i porównanie ich z wartościami
uzyskanymi w punkcie poprzednim,
4. Omówienie zaobserwowanego wpływu zmian wartości poszczególnych elementów R i C członu
na kształt charakterystyki skokowej (powiązać parametry RC ze stałą czasową T,
współczynnikiem wzmocnienia k, ewentualnie współczynnikiem α).
Sprawdzić zgodność obserwacji z wyliczeniami teoretycznymi. Omówić możliwe przyczyny
ewentualnych rozbieżności wyników doświadczalnych i teoretycznych oraz przedyskutować granicę
zaufania do wyników identyfikacji eksperymentalnej.
Przykład obliczeniowy
Człon inercyjny I-go rzędu
1) Schemat połączeń
R
1
X(s)=Uwe(s) Y(s)=Uwy(s)
Cs
a ka
0.638ka
uwe(t) uwy(t)