PRz AiS W6
UKAADY REGULACJI I DYNAMIKA
Transmitancja układu ze sprzężeniem zwrotnym. Serwomechanizm ze sterowaniem
prądowym. Układ elektro pneumatyczny. Dwa układy regulacji. Transformaty Laplace a.
Przebieg aperiodyczny. Przebieg aperiodyczny krytyczny. Przebieg oscylacyjny. Ekstrema
odpowiedzi oscylacyjnej. Wpływ dodatkowego zera i bieguna. Urządzenia.
TRANSMITANCJA UKAADU ZE SPRZŻENIEM ZWROTNYM
1. Połączenia szeregowe i równoległe
2. Sprzężenie zwrotne
Y(s)
Szukane
W(s)
Y RGoE RGoW RGoHY,
ponieważ E W HY
Y(s) RGo
W(s) 1 RGo H
Transmitancja układu z ujemnym sprzężeniem zwrotnym jest równa stosunkowi
transmitancji toru głównego przez jeden plus iloczyn transmitancji toru głównego i
transmitancji sprzężenia.
Pytanie. W jaki sposób zakłócenia wpływają na wyjście?
Yz
1
H(s)
Yz Go(Z U) GoZ GoU GoZ GoHRYz , ponieważ U RHYz .
Yz (s) Go
Yz (1 RGoH) GoZ
Z(s) 1 RGoH
Tym razem w torze głównym znajduje się a w torze sprzężenia RH.
Go,
SERWOMECHANIZM ZE STEROWANIEM PRDOWYM
1. Konstrukcja i połączenia
Serwomechanizm jest układem stabilizacji położenia silnika (elektrycznego,
hydraulicznego, pneumatycznego). Serwomechanizm z silnikiem elektrycznym
wykonuje ruch o zadany kąt lub liczbę obrotów.
Zastosowanie: obrabiarki (3D), suwnice, podajniki, roboty (5D, 6D), stacje radarowe
itp.
Mechanika serwomechanizmu
Realizacja nadajnika położenia:
enkoder (tarcza kodowa +
licznik fotoelektryczny)
potencjometr wieloobrotowy
rezolwer (maszyna elektryczna,
rodzaj selsyna)
Schemat elektryczny serwomechanizmu analogowego
2
Problem nastawienia sprzężenia zwrotnego
Jakie wartości r, należy ustawić, aby przebiegi regulacyjne były aperiodyczne
krytyczne o stałej czasowej T = 0.2s ?
Typowe przebiegi wielkości regulowanych:
Przebieg aperiodyczny pierwiastki rzeczywiste jednokrotne (różne)
Przebieg aperiodyczny krytyczny pierwiastki rzeczywiste wielokrotne
Przebieg oscylacyjny pierwiastki zespolone w mianowniku transmitancji
2. Elementy układu regulacji
Sterownik prÄ…du
U U 0
e
i
1
e ir 0
1
r
i e
r
ir
Silnik prądu stałego
Dane: MemN = 20 Nm, iN = 10 A,
J = 0.0005 kg m2 (moment bezwładności wirnika i przekładni)
Równanie momentów
d
J Mem Mo, Mo 0
dt
20 Nm
Mem = ks i, ks 2
10 A
Mo
-
Mem
i
1
ks Js 1
s
Silnik sterowany prÄ…dowo schemat blokowy
3
Przekładnia redukcyjna
= 0.01
Nadajnik położenia (potencjometr)
maksymalny kÄ…t obrotu, np. 5 obrotów zakres wejÅ›cia, = 5 · 2 rd
M M
Up = 12 V napięcie na potencjometrze zakres wyjścia
Proporcja:
y
U
, stÄ…d p
Up M y
y
U
M p
M
y
U
p
M
Współczynnik przeliczeniowy prostego elementu układu regulacji
oblicza siÄ™ jako zakres WY przez zakres WE .
zakres wyjścia przez zakres wejścia .
Tachogenerator
2
nN 12000 obr/min, UTGN 10V , 12000 rd/s
N
60
UTG
UTGN
3
kTG 7.95 10 V/rd/s
kTG
N
zakres przez zakres
Potencjometr tachogeneratora
UTG
UTG
² d" 1
Sumator połączenie przewodów elektrycznych
e = w y UTG błąd regulacji, w wielkość zadana
y wyjście układu (zmienna procesowa)
Rolę sumatora w układzie zrealizowanym cyfrowo pełni procesor sterownika.
4
3. Transmitancja serwomechanizmu
Schemat blokowy
I
E
max
Sprzężenie
tachometryczne
Sprzężenie pozycyjne
Schemat blokowy serwomechanizmu sterowanego prÄ…dowo
Simulink
Gain 1/r, Gain2 1/0.0005, Gain3 0.01, Gain4 12/(5*2*pi), Gain5 ß*0.00795
Upraszczanie schematu
Y (s)
Szukane dla Mo 0 - oddziaływanie sygnałów W i Mo można rozpatrywać
W (s)
niezależnie ze względu na zasadę superpozycji obowiązującą w układzie liniowym.
1) Połączenia szeregowe
max
2) Zwinięcie sprzężenia wewnętrznego 3) Standardowy układ ze sprzężeniem
zwrotnym
W
5
Tor główny dla sprzężenia pozycyjnego
ksÄ… U
1
p
AB Jr Åšmax Én2
s2
1 ks ² kTG s2 2¾Éns
1 AC
1
s Jr
Zbiorcze oznaczenia
ks U
ks kTG
p
2
, 2
n n
Jr
Jr
max
Transmitancja
2
n
2
Y(s) s2 2 s
n n
2 2
W(s)
s2 2 s
n n n
1 1
s2 2 s
n
2
współczynnik tłumienia
n
Y(s) W(s)
2 częstotliwość naturalna
n
s2 2 s
n n
4. Projektowanie dobór nastaw
Wymagania: a) przebiegi aperiodyczne krytyczne
b) stała czasowa 0.2
= 0 pierwiastek podwójny (przebiegi aperiodyczne krytyczne)
< 0 pierwiastki zespolone (przebiegi oscylacyjne)
> 0 pierwiastki rzeczywiste różne (przebiegi aperiodyczne zwykłe)
Serwomechanizmy nastawia siÄ™ na przebiegi aperiodyczne krytyczne.
= 0:
2 2 2 2 2
4 4 4 ( 1) 0
1
n n n
Tzam = 0.2 stała czasowa układu zamkniętego
2
2
1
n
n
= 1 :
2
s2 2 s (s )2 ( 1 s 1)2
n n n
n
1
Tzam 0.2,
5
n
n
Wniosek. Z warunków = 1, n = 5 należy wyznaczyć dwie niewiadome r, .
6
Obliczenia
ks U
2 0.01 12
p
2
25 r 0.61
n
Jr 25 0.0005 10
M
Rezystor 0.61 &! w torze sprzężenia sterownika prądu.
ks kTG
2 5 0.0005 0.61
2 2 1 5
0.192
n
3
Jr 2 7.95 10
Potencjometr należy ustawić na 19.2 %.
Simulink (schemat jak poprzednio)
Zmiany wielkości zadanej
Pulse Generator:
Amplitude 1
Period 5
Pulse Width (% of period) 50
Zakłócenie skokowe
Step
Final value 1
7
5. Transmitancja względem zakłócenia
Yz (s)
Szukane dla W=0
Z(s)
1) Schemat wyjściowy
Z
Yz
2) Przesunięcie Z na początek 3) Zwinięcie pętli wewnętrznej
Z
AB Én2
1 ks
Z
1 AC s2 2¾Éns
A
s Jr
Yz
Yz
Pętla ze sprzężeniem zwrotnym jest taka sama jak poprzednio. Transmitancja względem
zakłócenia ma więc postać:
Yz (s) r 1
Z(s) ks 1 2
s 1
n
1
Dane: r 0.61, ks 2, 5 ( 1)
n
0.2
8
6. Matlab
omn=5
dwaksiomn=10
ks=2
r=0.61
l=omn^2
m=[1 dwaksiomn 0]
l=[0 0 l]
t=0:0.02:2;
y=step(l,l+m,t);
plot(t,y), grid
yz=step(-(r/ks)*l,l+m,t);
plot(t,y,t,yz), grid
Sprzężenie zwrotne w prostym serwomechanizmie zmniejsza wpływ zakłócenia, ale
go całkiem nie eliminuje (wada). Mówi się, że pozostaje błąd ustalony (uchyb).
Wyjaśnienie l+m w step():
2
l
n
m s2 2 s
n
l
Licznik układu zamkniętego z jednostkowym
l
Gzam m
sprzężeniem zwrotnym jest taki sam jak licznik
l
l m
1 1
układu otwartego, a mianownik jest sumą
m
licznika i mianownika y = step(l, l+m, t).
Sumowany licznik i mianownik muszą mieć w Matlabie taką samą długość.
UKAAD ELEKTRO PNEUMATYCZNY
1. Realizacja techniczna
Zastosowanie układów elektro pneumatycznych
Regulacja poziomu w strefie zagrożonej wybuchem (petrochemia), strefie
zawilgoconej (cukrownictwo) i innych aplikacjach, gdzie nie jest wskazane lub wręcz
zabronione stosowanie aparatury elektrycznej.
9
0.02 ... 0.1 MPa
do 1.5 km 4 ... 20 mA
STEROWNIA
Y
do 1.5 km 4 ... 20 mA
6 mm
q
A
PPP
0.14 MPa
Elektro pneumatyczny układ regulacji poziomu
s s
Elementy składowe
ZR zawór regulacyjny z grzybkiem stożkowym
SP siłownik pneumatyczny
PPP pneumatyczny przetwornik poziomu
WZDP wzmacniacz dysza przesłona (odpowiednik wzmacniacza operacyjnego)
PE przetwornik pneumatyczno elektryczny
EP przetwornik elektro pneumatyczny
Trendy rozwojowe
Inteligentne przetworniki i urządzenia wykonawcze z własnymi procesorami do
prostych obliczeń, diagnostyki i komunikacji HART (1100/2200 Hz)
Zastępowanie sygnału analogowego 0/4...20 mA przez sygnał cyfrowy, np.
31.25 kbit/s (standardy protokołów: Profibus PA, Foundation Fieldbus, WFIP)
Przegląd urządzeń automatyki przemysłowej
www.katalogautomatyki.pl
2. Elementy układu na schemacie blokowym
Zakres WY - przez - zakres WE
3
Przetwornik poziomu PPP: 10 m 0.02...0.1 MPa kPPP 0.08 8 10 MPa/m
10
Przetwornik PE: 0.08 MPa 4...20 mA kPE 16 2 102 mA/MPa
0.08
2
Przetwornik EP: 16 mA 0.08 MPa kEP 0.5 10 MPa/mA
Siłownik SP: 0.08 MPa 3 cm kSP 3 37.8 cm/MPa
0.08
10
Zawór ZR: 3 cm 108 m3/h kZR 108 36 m3/h cm
3
Dane zbiornika (poprzedni wykład): 5 m2, h 10 m, q 108 m3/h
k kz S(s)
H(s) Qwe(s)
Ts 1 Ts 1
k 0.194 m m3/h , T 0.939 h, kz 0.8037 m/cm2 (8037 m/m2)
Schemat blokowy uwzględniający
urzÄ…dzenia
kZ
k1=1.309
E
H
1
kEP kSP kZR k
PID
Ts 1
[mA]
kPE
kPPP
k2=1.6
2
k1 kEP kSP kZR k 0.5 10 37.8 36 0.194 1.309 m/mA
3
k2 kPE kPPP 2 102 8 10 1.6 mA/m
Simulink
Simulink Extras Additional Linear PID Controller
Gain 1.309, Gain2 0.8037
Parametry zob. dalej
11
3. Regulator PID
Transmitancja
k
1 Td s 1 1 1
p
k 1 k 1 Td s k k T s k ki kd s
p p
Td D>1 p {ð 1ð 3ð
2ðd
Tis Tis s2ð p p s
3ð
s 1 1ðTi D
P
D I
kp wzmocnienie, kp, ki, kd wzmocnienia P, I, D
Ti czas całkowania, Simulink
Td czas różniczkowania, PID Controller (idealny)
D dzielnik Td (5& 10) PID Controller (with Approximate Derivative)
Regulator PI
Regulator PID przeważnie występuje jako PI. Otrzymuje się go poprzez ustawienie Td = 0.
1 kP Tis 1
kp(1 ) kp kp , Ti czas całkowania
{ð
Tis Tis Tis
{ð
P
I
Do sterowania obiektami inercyjnymi I rzędu, także z opóznieniem,
wystarczajÄ… regulatory PI.
4. Projektowanie
Schemat blokowy
Z="S
[cm2]
kZ
PI
k2
W 1 Y
kp 1 k1
Ts 1
[mA] Tis
[mA]
[mA]
kp, Ti
Jak dobrać ?
Najprostszym sposobem doboru Ti jest eliminacja stałej czasowej.
P
Y(s)
Transmitancja (Z=0)
W (s)
Ti T
Eliminacja:
Tis 1 k1k2
12
kP
Tis Ts 1
kPk1k2
k1 k2 1.309 1.6 2.1, K
T
K
K 1
Gzam s
K 1
s K
1 s 1
s K
1
Dane. Niech Tzam 0.1, co spowoduje, że układ zamknięty będzie prawie 10 razy
K
szybszy niż układ otwarty (T=0.939 h; zwykle jednak poprzestaje się na kilkukrotnym,
a nie 10 krotnym zwiększeniu szybkości).
1
Zatem 0.1 K 10 wzmocnienie układu otwartego
K
kp k1 k2
T K 0.939 10
Wzmocnienie kP regulatora PI K kP 4.5
T k1 k2 2 1
Nastawami regulatora PI są kp 4.5, Ti 0.939 (eliminacja stałej czasowej)
Simulink
Skok wielkości zadanej, a potem skok zakłócenia
13
Y(s)
Transmitancja zakłóceniowa (W = 0)
S(s)
Y(s) kz Ts 1
1
S(s) kP k1 Ts 1
s 1
K
Matlab
k1=1.309
kz=0.8037
T=0.939
kp=4.5
K=10
t=0:0.025:2.5;
l=[0 K]
m=[1 0]
y=step(l, l+m, t);
plot(t,y), grid
Odpowiedz zakłóceniowa:
lz=-kz*T*conv([1 0], l)
mz=kp*k1*conv([T 1], l+m)
yz=step(lz, mz, t);
plot(t, y, t, yz), grid
Regulator stopniowo eliminuje wpływ zakłócenia, ale trwa to wyraznie dłużej niż
czas regulacji tr 4 Tzam 0.4 . Dominującą stałą czasową jest bowiem T , a nie
1/ K.
14
DWA UKAADY REGULACJI
1. Serwomechanizm ze sterowaniem prÄ…dowym
Schemat uproszczony
y
w k
s2
k
p
k , Td
J
1 s
sprzężenie pozycyjne i tachometryczne
Transmitancja układu
k
2
Y(s) k
s2 n
2
(1 s) k
s2 2 s
W(s) s2 ks k
n n
1
s2
gdzie:
- częstotliwość drgań naturalnych (niegasnących), 1
2
n
k, k
n
- współczynnik tłumienia, 2
2. Układ regulacji poziomu
Rozdzielony regulator PI
y
I w
O
I O
k
1 k
p
y
w
s Ti
Ts 1
P
kp
P
kp kp
1 1 1
PI: kp (1 ) kp kp ki , gdzie ki
Tis s Ti s Ti
Transmitancja układu
O ki k
I
Y(s) IO
1 OP s Ts 1
O k ki k
W(s) 1 PO IO
1 I 1 kp
1 OP Ts 1 s Ts 1
15
kik
2
kik
n
T
2
1 kpk
Ts2 s(1 kpk) kik kik s2 2 s
n n
s2 s
T T
Gdyby cały regulator PI umieścić w torze głównym, to nie otrzymalibyśmy
standardowej transmitancji 2-go rzędu (nastąpiłaby zmiana w liczniku). Struktura
rozdzielona jest preferowana przez praktyków dla uniknięcia przeregulowania (zob.
dalej). Przeregulowania można też uniknąć filtrując wielkość zadaną.
Wniosek. Proste serwomechanizmy i układy automatyzacji procesów są opisane
standardowymi transmitancjami 2-go rzędu.
1. Pierwiastki mianownika bieguny
Transformata Laplace a odpowiedzi skokowej
2
1 1
n 2 2 2 2 2
W(s) , Y(s) , 4 4 4 ( 1)
n n n
2
s s2 2 s s
n n
Zależnie od współczynnika tłumienia pierwiastki mianownika są:
> 1 rzeczywiste różne
= 1 jednakowe (pierwiastek podwójny)
< 1 zespolone
Kształt odpowiedzi zależy przede wszystkim od współczynnika tłumienia , a czas od
częstotliwości .
n
TRANSFORMATY LAPLACE2 A
1. Wzory ogólne
t
F(s)
f ( )d
o
s
16
2. Transformaty elementarne
3. Funkcje trygonometryczne
4. Stała czasowa
17
PRZEBIEG APERIODYCZNY
1. ¾ > 1 " > 0
2
Dwa pierwiastki rzeczywiste różne s1, s2 , takie że s1 s2 n
Rozkład na ułamki proste:
2
1 R0 R1 R2
n
Y(s)
(s s1)(s s2) s s s s1 s s2
Residua R0, R1, R2 oblicza się metodą przysłaniania.
2. Pierwiastki jednokrotne
Przypadek ogólny
tzw. metoda przysłaniania
Kropki po prawej stronie rozwinięcia F(s) reprezentują ułamki odpowiadające
pierwiastkom pozostałej części mianownika.
3. Obliczenia
p1 0, p2 s1 , p3 s2
2 2 2 2
n n n n
R0 |s 0 1
2
(s s1)(s s2) ( s1) ( s2) s1s2 n
Oznaczmy: s1 1 , s2 1
T1 T2
1
2
s1s2 s2 T2 T1
R1 n |s s1
s(s s2) s1 (s1 s2) s1 s2 1 1 T1 T2
T1 T2
Podobnie
1
2 2
s1 T1 T2
n
R2 n |
s s2 s2 (s2 s1) s2 s1
1 1
s(s s1) T1 T2
T2 T1
18
Odpowiedz
t t
T1 T1 T2 T2
y(t) 1 e e
T1 T2 T1 T2
Jest to tzw. przebieg aperiodyczny zwykły (o dwóch stałych czasowych).
PRZEBIEG APERIODYCZNY KRYTYCZNY
1. ¾ = 1 " = 0
Jeden pierwiastek podwójny s1,2 n
2
1 R0 R1 R2
n
Y(s)
(s )2 s s s (s )2
n n n
2. Pierwiastki wielokrotne
Przypadek ogólny
Rezidua Ri oblicza się również metodą przesłaniania, ale zaczynając od końca , tj.:
,
&
Ponieważ , więc
W szczególności dla m=3 mamy
,
19
1
gdzie R3 1( p), R2 2( p), R1 3( p),
2
3. Obliczenia
2
1
n
p0 0, p , Y(s)
n
(s )2 s
n
2 2
2 2
1 d
n n n
R0 |s 1, R2 n |s , R1 ( ) |s |s 1
n
s2 n
(s )2 0 s 1! ds s
n n
n
1 1
n
Y(s)
s s
(s )2
n
n
Na podstawie drugiej tabeli transformat mamy
t t
n n
y(t) 1 e te
n
1
Niech T . Wtedy
n
t
t
T
y(t) 1 (1 )e
T
Jest to przebieg aperiodyczny krytyczny. W nastawianiu serwomechanizmów chodzi o
uzyskanie właśnie takiego przebiegu.
Zapisując Y(s) za pomocą stałej czasowej T jako
1 1
Y(s)
(Ts 1)2 s
powyższy wzór otrzymuje się wprost z czwartej tabeli transformat Laplace a.
PRZEBIEG OSCYLACYJNY
1. ¾ < 1 " < 0
Dwa pierwiastki zespolone
2
s1,2 j 1 j ,
n n
2
gdzie: , 1 częstotliwość drgań tłumionych
n n
Ponieważ (s s1)(s s2 ) (s )2 2 , więc
2
1
n
Y(s)
(s )2 2 s
2. Pierwiastki zespolone
Przypadek ogólny
20
gdzie
Ponieważ , więc
Przysłaniając obydwa czynniki i podstawiając s j otrzymuje się
,
zatem
Na podstawie trzeciej tabeli transformat
Dlatego S i C sÄ… nazywane amplitudami sinusa i cosinusa .
3. Obliczenia
Zapis wyjściowy
R0 C(s ) S
Y(s)
s (s )2 2
2 2 2 2
n n n n
R0 |s 1
2
(s )2 2 0 2 2 ( n )2 ( n 1 2 )2 n
Amplitudy S, C
S C(s )
t t
Se sin t Ce cos t
,
(s )2 2 (s )2 2
Według wzoru ogólnego
2 2
1 1 1
n n
S jC (s )2 2 Y(s) ( )
s j
s j
s j
2
j 1
n n
( j ) j j j S jC
2 2
2
2
1 1
n
Zatem
S , C 1
2
1
Odpowiedz skokowa
t
y(t) 1 e (S sin t C cos t)
Przekształcenia trygonometryczne:
sin( ) sin cos cos sin , gdzie t
21
S C
S sin t C cos t S2 C2 (sin t cos t ) S2 C2 sin( t ) ,
2
S22ð4ð2 S C2
C
1ð4ð 3ð 1ð4ð 3ð
2ð4ð
cos sin
gdzie
2
1
sin C 1
2
arctg arctg arctg , S C2
cos S 2
1
1
t
n
y(t ) 1 e sin( t )
2
1
Jest to przebieg oscylacyjny.
4. Matlab
omn=1
y
l=omn^2
% >1
ksi=2
m=[1 2*ksi*omn omn^2]
t=0:0.1:20;
y=step(l,m,t);
% =1
ksi=1
m=[1 2*ksi*omn omn^2]
y1=step(l,m,t);
% <1
ksi=0.5
t
m=[1 2*ksi*omn omn^2]
y2=step(l,m,t);
plot(t,y,t,y1,t,y2),grid
Uwaga. Na lekko oscylacyjne przebiegi nastawia się układy automatyzacji
procesów technologicznych (ciśnienie, poziom, przepływ) ze względu na dobre
tłumienie zakłóceń.
EKSTREMA ODPOWIEDZI OSCYLACYJNEJ
1. Wyznaczenie ekstremum
W celu wyznaczenia ekstremów odpowiedzi oscylacyjnej II rzędu można zastosować jedną z
dwóch metod:
dy
określenie wyrażenia na pochodną w sposób tradycyjny i przyrównanie jej do zera
dt
(metoda często pracochłonna),
skorzystanie z właściwości transformaty Laplace a, tzn. ze wzoru na transformatę pochodnej:
22
1
2
2
n
dy
1 1 t
n
L s Y(s) L e sin t 0
dt (s )2 2
ti , 2 , 3 ..., ti i , ti i , ti i ,
i 1, 2, 3...
2
1
n
t1
2
1
n
2
1
t1 t1 t1 t1
A1 y(t1) 1 e S sin t1 C cos t1 e S sin C cos C e e e ,
ponieważ , , t1
C 1
n
2
1
n
2. Parametry odpowiedzi oscylacyjnej
czas narastania tn (10 do 90%)
przeregulowanie p% A1 100% ,
A3
czas regulacji tr ,
stopień tłumienia d
A1
Podstawowe znaczenie majÄ… p% i tr.
Przykład
1
n
0.3
1
T 3.333
n
tr 4T 13.332
23
3. Przeregulowanie p%
p%
ln
2
100
1
p% A1 100% e 100%
p%
2 2
ln
100
100
Matlab p%( )
p% 90
80
70
ksi=0.01:0.01:0.99;
60
p=exp(-pi*ksi./sqrt(1-ksi.*ksi))*100;
50
plot(ksi,p),grid
40
30
20
Dzielenie i mnożenie z kropką
10
(wektory).
0
0 0.2 0.4 0.6 0.8 1
Wzór przybliżony p% (1 ) 100%, 0.6
0.6
100
p%
90
80
Matlab
70
60
ksi=0.01:0.01:0.6;
50
p=exp(-pi*ksi./sqrt(1-ksi.*ksi))*100; 40
30
paproks=(1-ksi/0.6)*100;
20
plot(ksi,p,ksi,paproks),grid
10
0
0 0.1 0.2 0.3 0.4 0.5 0.6
4. Czas regulacji tr
t
t
T
e
4
T
tr
0.0 1.0 2%
n
1.0 0.3679
2.0 0.1353
3.0 0.0495
5%
3 4.6
tr , tr
4.0 0.0183
2%
5% 1%
n n
4.6 0.0100 1%
5. Serwomechanizm i regulator PI
Dla serwomechanizmu prądowego (zob. wcześniej) znalezć k i , gdy dane są Tzam= 0.1, p%= 0.
p% = 0 oznacza 1 przebiegi aperiodyczne krytyczne.
4 4
tr 4 Tzam 0.4 10
n
n n
24
2
k
n 2
Gzam , 2 k , k k 100, 0.2
n n
2
s2 k s k s2 2 s
n n
Dobrać nastawy regulatora PI w układzie regulacji poziomu dla danych p% = 20%, tr = 0.5.
1
ki k
s Ts 1
kp
ki k
2
T n
k , ki ?
Gzam , (k = 2.1, T = 0.939 h)
p
2
1 k k
ki k
s2 2 s
n n
s2 p s
T T
p%
ln
100
4
0.456, tr 0.5 17.54
n
p% 20%
n
2 2
ln
100
1 k k
ki k
p
2
17.542 ki 137.6 , 2 k 6.7
n n p
T T
kp 6.7
Ti 0.0487 h 2.92 min 175 s
ki 137.6
WPAYW DODATKOWEGO ZERA I BIEGUNA
1. Wpływ zera
Układ regulacji poziomu ze standardowym regulatorem PI (nierozdzielonym) ma postać jak
niżej
k
1 k 2.1, 17.54
n
kp 1
Ti 0.0487, 0.456
Ts 1
Tis
k 6.7
p
1
s
Ti
1
kp 1 kp
PI:
Tis s
kkp
s 1/Ti k
Tis 1
k
p
k k s 1/Ti
TTi
p
s Ts 1
Gzam
s 1/Ti k
Ts2 s kkp s 1/Ti s2 1 kkp s kkp
1 k
p
s Ts 1
T TTi
25
kkp 2 1 kkp
, 2
n
TTi n T
17.542 0.0487s 1
dodatkowe zero w liczniku
Gzam
s2 2 0.456 17.54 s 17.542
Matlab
y
omn=17.54
ksi=0.456
l=omn^2
m=[1 2*ksi*omn omn^2] pierwotnie
t=0:0.01:2;
y=step(l,m,t);
Ti=0.0487
% dodatkowe zero
l=conv(l,[Ti 1])
t
y1=step(l,m,t);
plot(t,y,t,y1), grid
Jak widać przeregulowanie wzrosło (niekorzystnie).
Uwaga. Przez rozdzielenie PI na P+I, a PID na PI+D (ewentualnie I+PD) można uniknąć
wprowadzania zer do licznika transmitancji układu zamkniętego, czyli nadmiernych
przeregulowań.
2. Wpływ bieguna
2
1
n
Gzam
2
Ts 1 s2 2 s
n n
Matlab
y
omn=17.54
ksi=0.456
T=0.125
l=omn^2
m=[1 2*ksi*omn omn^2]
t=0:0.01:2;
y=step(l,m,t);
% dodatkowy biegun
m=conv(m,[T 1])
y2=step(l,m,t);
plot(t,y,t,y2),grid
t
Uwaga. Pojawienie siÄ™ dodatkowego bieguna powoduje eliminacjÄ™, a przynajmniej
zmniejszenie przeregulowania (korzystnie). Zatem przeregulowanie można
wyeliminować przez dodanie filtru wielkości zadanej przed wprowadzeniem jej na
układ ze sprzężeniem zwrotnym.
26
URZDZENIA
www.katalogautomatyki.pl
1. Serwomechanizm
Enkodery obrotowe
Silniki napędowe Sterownik napędu
Robot
27
2. Układ regulacji
Radarowe przetworniki Regulatory PID
poziomu
Siłownik pneumatyczny Siłowniki elektryczne
28
Wyszukiwarka
Podobne podstrony:
W6 Układy regulacji i dynamika AiS 2013W6 Dwa układy regulacjiW1 Układy kombinacyjne AiS 201309 Układy regulacjiUklady regulacji DUN2 Dyskretne układy regulacji, rozdział 3 i 4 Funkcje dyskretne Równania różnicoweid497W8 Linie pierwiastkowe Evansa AiS 2013Uklady regulacji napedow? DUNUklady regulacji zRegulamin Sztafety Studenckiej 201309 Uklady regulacjiC w6 zmienne dynamiczne wskazniki funkcji2013 nr 26 Gaz łupkowy – nowe regulacjezadania zestaw 5 dynamika uklady nieinercjalneUklady Dynamiczne Zad ser II p14 Statyczne i dynamiczne wlasciwosci regulatorowwięcej podobnych podstron