Ingineria Reglarii Automate - Note de Curs
Ingineria Reglarii Automate - Note de Curs
Ingineria Reglarii Automate - Note de Curs
NOTE DE CURS
REV.01
1
CUPRINS
Regulatoare convenţionale
2
INTRODUCERE - elemente recapitulative
1.NOŢIUNEA DE SISTEM DINAMIC
Proces şi sistem
Procesul (fizic) reprezintă tranziţia unui sistem dintr-o stare în
alta, starea sistemului fiind descrisă prin parametrii acestuia. In
cadrul unui proces tehnologic putem avea transfer transfer masic
şi/sau transfer de energie, notat generic cu Wi la intrarea în
proces şi cu We la ieşirea din proces.
Wi We
PROCES
fig.1
Funcţionarea normală a procesului (neglijând pierderile ) este dată
de relaţia în regim staţionar
Wi - We = 0 (1.1)
Ecuaţia (1.1) caracterizează regimul de echilibru / regimul
staţionar.
În cazul existenţei, în cadrul procesului, a unor fenomene de
acumulare/ dezacumulare vom avea:
Wi - We ≠ 0 (1.2)
Aceste fenomene sunt descrise de mărimile de stare, astfel vom
avea:
ẋ=W i −W e (1.3)
Acumulare are loc atâta timp cât ẋ≠0
t
x t =x 0∫ W i−W e d (1.4)
0
3
x t =−
W
a
W
a
x 0 e at (1.7)
W
Termenul constant − a descrie regimul staţionar iar cel variabil
descrie regimul tranzitoriu care se amortizeză asimtotic (a<0).
Pentru regim staţionar axW i −W e =0 (1.8)
Cu cât „nivelul de acumulare” este mai scazut cu atât consumul este
mai mic.
Fără autoechilibrare
x(t)
Regim
W
− staţionar
a
Cu autoechilibrare
tg α'=ΔW+ax(0)
tg α=ΔW
t
fig. 2
Mărimile de intrare pot fi grupate în două categorii:
comenzi - notate u(t) – care acţionează în mod voit asupra
sistemului; valorile comenzilor sunt cunoscute de către
proiectant, t-timp.
perturbaţii - notate v(t) – care acţionează în mod nedorit
asupra sistemului; de multe ori perturbaţiile au un caracter
aleator. v(t)
u(t) y(t)
x(t) stare
fig.3
Considerăm : Wi = bu şi We=-ev (pp că mărimea perturbatoare este de
sens opus consumului), x mărimea de stare, z mărimea de calitate
(ieşirea optimă/prescrisă).
ẋ=axbuev
y =cx (1.9)
z=dx
4
Noţiunea de sistem dinamic
5
Un sistem dinamic liniar continuu şi invariant în timp poate fi
descris printr-un model matematic de tip ecuaţie diferenţială
liniară de ordinul n, cu coeficienţi constanţi şi reali:
a n y ( n ) (t ) a n 1 y ( n1) (t ) a0 y (t ) bm u ( m ) (t ) bm1u ( m1) (t ) b0 u (t ) , (2.1)
unde:
u(t) reprezintă mărimea de comandă aplicată,
y(t) desemnează mărimea de ieşire,
iar t reprezintă variabila independentă timp.
6
lui Q(s) se numesc zerourile sistemului.
7
DEFINIREA UNUI SISTEM DE REGLARE AUTOMATĂ (SRA)
Perturbaţii
Mărimea Mărimea
de intrare Comandă de ieşire
Dispozitiv
Proces
comandă
8
Sistem în circuit închis sau Sistem de Reglare Automată (cu
reacţie) (SRA):
v
Element de r ε u Elemente m
Regulator Procesul y*
impunere + de execuţie
(R) condus (P)
(EI) -- (EE)
y
Traductoare (Tr)
SRA - reprezentare funcţională
9
referinţei r. Aducerea valorii lui y la valoarea impusă prin r
se asigură de către regulatorul R, care elaborează o comandă
în tensiune u, aplicată elementului de execuţie EE. Elementul
de execuţie, prin mărimea m, acţionează asupra procesului
modificând mărimea reglată y* în sensul dorit. Valoarea lui y*
este măsurată cu ajutorul traductorului Tr şi convertită în
semnalul electric y , care prin calea de reacţie negativă este
adusă la intrarea sistemului. Diferenţa dintre referinţa r şi
măsura yr este abaterea/eroarea ε, mărime care se aplică
regulatorului R. Dacă ε ≠ 0 (ε >0 sau ε<0) regulatorul
elaborează o comandă u pe baza legii de reglare (algoritmului
de reglare) implementate în structura regulatorului, până când
y atinge valoarea impusă de referinţa r. Dacă ε = 0 abaterea
este nulă (r = y) atunci şi comanda regulatorului este de
asemenea nulă.
Diferenţele dintre referinţă şi mărimea reglată apar ca
urmare a modificării valorii mărimii de referinţă sau prin
acţiunea perturbaţiilor.
Ţinând seama de sensul fluxului informaţional,
regulatorul R, elementul de execuţie E şi procesul tehnologic
P, formează calea directă a sistemului automat, iar
traductorul calea de reacţie.
Sistemul automat îşi îndeplineşte cerinţele impuse dacă
abaterea sau eroarea de reglare , definită ca fiind diferenţa
dintre referinţa r şi ieşirea măsurată y tinde către zero în
regim staţionar, indiferent de acţiunea perturbaţiilor
externe care acţionează asupra procesului condus. Această
condiţie se exprimă astfel:
lim ε(t )=lim [ r ( t)− y(t) ]
t→∞ t→∞
v
Partea fixă
+ ε u y
r R EE P Tr
--
10
Structura unui SRA cu reacţie unitara
Clasificarea SRA
Există mai multe criterii de clasificare a SRA:
A. După modul de evoluţie în timp a mărimii de referinţă;
B. După modul de reprezentare a informaţiei în regulator;
C. După numărul mărimilor reglate;
D. După distribuţia spaţială a mărimii reglate
11
proceselor lente cu timp mort. Prezenţa unui număr mare de
constante de timp în funcţia de transfer a procesului face
dificilă utilizarea unor algoritmi de reglare tipizaţi.
Pentru procesele tehnologice la care se pot evidenţia
mărimi intermediare măsurabile, iar funcţia de transfer a
procesului poate fi scrisă ca un produs de funcţii de transfer
care nu conţin mai mult de două constate de timp, se recomandă
reglarea în cascadă. Admiţând că procesul condus este
decompozabil în subprocese interconectate cauzal, cu variabile
intermediare accesibile măsurării şi cu dinamică
descrescătoare de la intrare spre ieşire, se poate alcătui o
structură de reglare în cascadă folosind un număr de
regulatoare egal cu numărul variabilelor măsurate din proces.
În cazul reglării în cascadă are loc o reglare simultană
a mai multor mărimi din cadrul procesului, ceea ce determină o
reducere însemnată a duratei procesului tranzitoriu, mai ales
dacă mărimile intermediare sunt alese astfel încât să răspundă
mai repede decât mărimea de ieşire la perturbaţiile care
acţionează asupra procesului.
V2(s) V1(s)
U2(s) Y2(s)
W(s) E1(s) U1(s) E2(s)
+ Y1(s)
HR1(s) HR2(s) HE(s) HP2(s) HP1(s)
_ _
Yr1(s) Yr2(s)
HT2(s)
HT1(s)
12
care la rândul ei determină cauzal evoluţia variabilei de
ieşire y1 din proces. Regulatorul automat secundar, cu funcţia
de transfer HR2(s), este destinat reglării variabilei y2 şi
compensării acţiunii perturbaţiei v2, iar regulatorul
principal cu funcţia de transfer HR1(s) are rolul de a asigura
realizarea funcţiei de reglare în raport cu referinţa w. De
asemenea furnizează referinţa pentru regulatorul secundar.
Cele două regulatoare din cadrul acestei structuri
funcţionează în regim de urmărire (master–slave), adică
mărimea de comandă a regulatorului principal este mărimea de
referinţă pentru regulatorul secundar. Se recomandă ca bucla
interioară să aibă o viteză de răspuns mai mare decât bucla
principală.
Acordarea regulatoarelor în cadrul structurilor de
reglare în cascadă se iniţiază cu bucla interioară,
deconectând regulatorul principal şi se continuă cu
regulatorul principal, folosind procedurile de acordare
specifice .
13
STRUCTURA DE REGLARE DUPĂ PERTURBAŢIE ŞI
STRUCTURA DE REGLARE COMBINATĂ
uV
Rv
v
-
w u m y
+ R E P
+
-
yr
Structura de reglare combinată.
14
intrarea părţii fixate a sistemului să rezulte compensarea
influenţei exercitate de perturbaţia v asupra mărimii de
ieşire y. În acest scop, întrucât perturbaţia este considerată
cu semnul plus, este necesar ca mărimea uV sumată algebric în
cel de-al doilea sumator cu mărimea de ieşire u a
regulatorului R (care realizează reglarea după eroare) să fie
aplicată sumatorului cu semnul minus.
Reglarea după perturbaţie este mult folosită în
conducerea instalaţiilor de încălzire a clădirilor, unde
aportul de căldură asigurat de corpurile de încălzire trebuie
să ţină seama, printre altele în primul rând, de temperatura
exterioară, care reprezintă perturbaţia principală. Dacă
reglarea s-ar efectua numai în funcţie de temperatura
interioară din clădire, atunci când temperatura exterioară
variază, evident se modifică şi temperatura interioară după un
timp dat de constantele de timp ale clădirii (care au valori
importante de ordinul orelor) şi în mai mică măsură de timpul
mort de ordinul minutelor necesar agentului termic pentru
circulaţia prin instalaţia de încălzire. Eroarea staţionară a
sistemului automat (diferenţa dintre temperatura interioară
dorită în clădire şi temperatura interioară reglată) va fi cea
mai mare parte a timpului diferită de zero şi este posibil să
atingă valori importante care crează disconfort pentru
ocupanţii clădirii. De asemenea, performanţa cea mai
importantă a sistemului automat, stabilitatea, poate fi mult
diminuată. Reglarea după temperatura exterioară (perturbaţia)
este o reglare cu anticipaţie, care modifică cantitatea de
căldură produsă de instalaţia de încălzire funcţie de valoarea
temperaturii exterioare, fără să se aştepte modificările
nedorite ale temperaturii interioare.
Reglarea după perturbaţie se utilizează de asemenea la
sistemele de reglare automată a temperaturii apei din piscine.
Problema stabilităţii
15
a comportării sale, se poate obţine prin determinarea
variaţiei în timp a mărimii de ieşire y ca urmare a unei
variaţii a mărimii de referinţă w, sau ca urmare a variaţiei
unei perturbaţii v. Cunoscând astfel răspunsul sistemului la
variaţiile mărimilor de intrare, pot fi măsuraţi indicii de
calitate ai regimului staţionar şi tranzitoriu, respectiv pot
fi determinate performaţele staţionare şi tranzitorii,
stabilindu-se în consecinţă dacă sistemul analizat poate fi
sau nu poate fi utilizat în cazul concret dat.
O condiţie necesară, dar nu suficientă, pentru ca un
sistem automat să poată fi utilizat în practică este
stabilitatea sistemului, respectiv proprietatea sistemului de
a restabili, prin acţiunea sa, un nou regim staţionar, atunci
când datorită variaţiei mărimilor de intrare a fost scos
dintr-un regim staţionar anterior. Astfel, regimul tranzitoriu
al unui sistem automat stabil are o durată limitată.
Un sistem instabil nu este utilizabil, deoarece nu poate
îndeplini scopul pentru care este creat, acela de a realiza pe
cale automată o anumită lege de dependenţă între mărimea de
ieşire şi cea de intrare. La sistemele instabile, mărimea de
ieşire are variaţii necontrolate, deoarece un nou regim
staţionar nu mai este restabilit după ieşirea dintr-un regim
staţionar anterior.
Determinarea performanţelor unui sistem are deci sens
numai în cazul când sistemul este stabil şi, de aceea, analiza
comportării sistemelor cuprinde o verificare prealabilă a
stabilităţii.
16
CRITERIILE DE PERFORMANŢĂ IMPUSE UNUI SRA
v1
v2
+ ε u y
r HR(s) Σ HF(s) Σ
--
v1 + y
HF(s)
--
HR(s)
Ieşirea sistemului în raport cu perturbaţia v1
y
v2 +
--
HR(s) HF(s)
17
Performanţele generale unui sistem sunt definite pentru regimul
tranzitoriu şi cel staţionar:
Criteriile generale de performanta ale unui sistem sunt determinate
prin
analiza raspunsului in timp:
► răspuns tranzitoriu:
suprareglaj
timp tranzitoriu tt (durata regimului tranzitoriu)
factor de amortizare
timp de creştere tc
timp de întârziere ti
si prin analiza raspunsului in frecvenţă:
► performanţele în domeniul frecventelor :
stabilitatea sistemului
precizia în regim staţionar sau eroarea staţionara
st
banda de frecventa
marginea de fază M
marginea de amplitudine Mc
pulsaţia de rezonanta R
valoarea de vârf a modulului Mv
y
yref +
18
Ex: răspunsul oscilant este îmbunătăţit daca:
¥ ¥
yref -
-
+
+ +
O
t
19
Performanţele sistemelor automate determinate pe
baza răspunsului indicial
k
H (s)
1 sT
y 1 (t ) £-1 , k1(t ) k e T
k 1 e T
pentru t 0,
s s 1
T
t
y (t ) k 1 e
1 T
Răspunsul indicial de mai sus are două componente:
20
componenta permanentă;
componenta tranzitorie.
Componenta permanentă se determină folosind teorema
valorii limită finale:
k
y p (t ) lim sY ( s ) lim s k pentru t 0
s 0 s 0 s (1 sT )
Componenta tranzitorie a răspunsului indicial este:
t
y t (t ) k e T
d
t t
dy 1 (t ) 1 k
t 0 k (1 e T
t 0 0 k e T t 0
dt dt T T
k
unde este unghiul cu abscisa al tangentei de pantă maximă.
T
w
,y st
1 w
k
T t
ttr
21
tranzitoriu şi eroarea staţionară.
Timpul tranzitoriu reprezintă timpul necesar pentru ca
răspunsul tranzitoriu al sistemului să intre în banda (10,05)
yst fără a o mai părăsi ulterior.
Tangenta în origine intersectează răspunsul staţionar
după un timp egal cu constanta de timp a sistemului.
y st y p k
t
t
0,95k k 1 e T
tt
t t T ln 0,05 3T
e T
1 0,95
tt
tt
ln e T
ln 0,05 ln 0,05
T
Eroarea staţionară a sistemului se defineşte ca
diferenţa dintre mărimea de intrare (referinţa) şi valoarea
staţionară a mărimii de ieşire a sistemului:
Pentru intrare treaptă unitară:
εst =r (t)− y st (t )=1−k
ω 2n
H (s)=
s2 + 2 ξ ωn s+ ω2n
U(s) Y(s)
p 2 n p 0
2 2
n
p1, 2 n j n 1 2
22
h(t) 0<ξ<1
ξ=0
t
ξ>1
Răspunsul indicial al unui sistem oscilant de ordinul II pentru
diferite valori ale lui ξ
h(t) ε st
1.05 σ
1
0.95
0.05
tc
tt
23
Suprareglajul sau abaterea dinamica maxima σ este diferenţa
intre valoarea maxima a ieşirii şi valoarea de regim staţionar:
σ= y max − y st
y max − y st
. Se poate defini procentual ca reprezentând σ= 100
y st
πξ
−
Pentru sisteme de ordinul II σ II =e √ 1− ξ =>
2
σ= f ( ξ)
4
t t≈
ξ ωn
24
Modelele matematice ale elementelor tip ale unui SRA sunt:
mort
Element cu timp e−τ s
mort şi dy H (s)=K 0
T + y (t)=K 0 u (t−τ) Ts+ 1
întârziere de dt
ordinul I
Element cu timp 2
d y dy K 0 e−τ s
mort şi T 1T 2 2
+ (T 1+ T 2) + ω2n y=K 0 u (t−τ) H (s )=
dt dt (T 1 s +1)(T 2 s+1)
întârziere de
ordinul II
Element de du H (s)=Ts+ 1
y ( t)=T + u (t)
anticipaţie de dt
ordinul I
Element de d u du
2 2
H (s)=s + 2 ξ ω n s+ ω n
2
anticipaţie de y (t)= 2 + 2 ξ ωn + ω 2n u
dt dt
ordinul II
Regulator
proporţional-
intergral-
u(t)=K R ε+
( 1
Ti ∫ ε(t )dt+ T d
dε
dt ) H ( s)=K R 1+
( 1
Ti s
+Tds
)
derivativ
Element 2
d y dy K0
integrator cu T 2
+ =K 0 u (t ) H ( s)=
dt dt s( 1+ T i s)
întârziere de
ordinul I
Circuit de T 1 s+ 1
corecţie de
anticipaţie /
T1
dy
dt
du
+ y (t)= K 0 T 2 + u (t)
dt ( ) H ( s)=K 0
T 2 s+ 1
întârzâiere
25
ANALIZA SISTEMELOR DE REGLARE AUTOMATĂ
PE BAZA METODELOR DE FRECVENŢĂ
26
partea reală σ =0, deci H(s)=H(σ+jω)=H(jω).
Un semnal de intrare, periodic sau nu, poate fi considerat
ca fiind o suprapunere de semnale sinusoidale rezultate din
analiza sa armonică. Răspunsul unui sistem liniar la un astfel
de semnal poate fi obţinut prin suprapunerea răspunsurilor la
semnalele sinusoidale care îl compun.
Considerăm intrarea W(s) şi ieşirea Y(s) pentru un
sistem.automat
w y
SA
w=Wcos(ωt+φi) (4.1)
y=Ycos(ωt+φe) (3.2)
Aplicând relaţiile cunoscute de la trigonometrie se obţin
w=W(cosωt cosφi – sinωt sinφi)
(3.3)
V(ω)=Im H(jω)
Reprezentarea defazajului
27
Y w
W y
I t
e
W
W ( s) ( s cos i sin i ) (3.5)
s 2
2
Y
Y ( s) ( s cos e sin e ) (3.6)
s 2
2
Y j cos e sin e
H ( j ) (3.8)
W j cos i sin i
Y cos e j sin e
H ( j ) (3.9)
W cos i j sin i
28
Ţinând seama de formulele lui Euler
Y e j e Y j ( e i )
H ( j ) e (3.10)
W e j i W
Y ( ) j [ e ( ) i ]
H ( j ) e H ( )e j ( ) (3.12)
W
H ( j ) H ( )e j ( ) H Re ( ) jH Im ( ) (3.13)
H ( ) [ H Re ( )] 2 [ H Im ( )] 2 (3.14)
H Im ( )
( ) arctg (3.15)
H Re ( )
29
adică defazarea între semnalul de la ieşire şi cel de la
intrare.
30
Im Im
K
α=3 K
θ=3
Re
θ=0 Re
α=1
ω→0 ω→∞
31
ca un sistem SLIT sa fie stabil în stare închisa, este ca
punctul de coordonate (-1, j0) sa fie înconjurat în sens
trigonometric de caracteristica amplitudine-faza a sistemului
deschis (trasata pentru ω variind de la -∞ la +∞), de un număr
de ori egal cu numărul polilor situaţi în semiplanul drept al
funcţiei de transfer a sistemului închis.
Criteriul simplificat este o particularizare a celui
general şi porneşte de la ipoteza sistemului SLIT stabil în
stare deschisa ( numărul polilor din semiplanul drept al
planului complex (σ, jω) este 0).
32
stânga punctul (-1, j0).
Pulsaţia de antifază ωπ .
Pulsaţia de antifază ωπ , reprezintă cea mai mică pulsaţie
pentru care caracteristica complexă
Hd ( jω) taie semiaxa reală negativă, Se obţine din,
ωπ = min {ω|ϕd (ω) = −π}
Banda de pulsaţie ωb .
Banda de pulsaţie, marcată prin valoarea maximă ωb , se
defineşte ca fiind intervalul de pulsaţie [0,ωb ], unde ωb este
cea mai mică valoare pentru care A(ω)<
√ (2) A(0) ∀ω> ωc
2
34
Cu cat Mc şi Mj sunt mai mari cu atât gradul de stabilitate
al sistemului automat este mai mare. Cu cat locul de transfer
este mai aproape de origine, lăsând mult în stânga punctul
(-1, j0), cu atât sistemul este mai stabil.
AdB = 20log H ( jw )
argH j
Caracteristici Bode
La ambele caracteristici, pulsaţia în abscisă se consideră în
scară logaritmică.
Fie H(j) răspunsul la frecvenţă al unui sistem, exprimat ca
produsul funcţiilor de transfer H1 ( jw ) şi H 2 ( jw ) . Caracteristicile
Bode pentru H ( jw ) = H1 ( jw ) × H 2 ( jw ) se deduc pe baza celor aferente
funcţiilor H1 ( jw ) şi H 2 ( jw ) , astfel:
20lg|H(jω)|=20lg|H1(jω)+20|H2(jω)||
arg H ( jw ) = arg H1 ( jw ) + arg H 2 ( jw )
sau
AdB (w ) = A1dB (w ) + A2dB (w ) , j (w ) = j1 (w ) + j 2 (w ) (111)
35
însumare, din cele aferente funcţiilor de transfer H1 ( jw ) şi
H 2 ( jw ) .
Drept cosecinţă, pentru construcţia diagramelor Bode se
factorizează atât numitorul cât şi numărătorul în elemente de
ordinul unu sau doi.
La numărător, pentru caracterstica amplitudine -
pulsaţie: Fiecare element de ordinul 1 introduce o pantă de
+20 dB, iar fiecare elent de ordinul 2 introduce o pantă de
+40 dB.
La numitor, pentru caracterstica amplitudine - pulsaţie:
Fiecare element de ordinul 1 introduce o pantă de -20 dB, iar
fiecare elent de ordinul 2 introduce o pantă de -40 dB.
La numărător, pentru caracterstica fază - pulsaţie:
Fiecare element de ordinul 1 introduce o pantă de +90º, iar
fiecare elent de ordinul 2 introduce o pantă de +180º.
La numitor, pentru caracterstica fază - pulsaţie: Fiecare
element de ordinul 1 introduce o pantă de -90º, iar fiecare
elent de ordinul 2 introduce o pantă de -180º.
36
REGULATOARE CONVENŢIONALE
H R (s )=K R
Ε(s) U(s)
u( t)=K R ε(t) ;ε(t)=1 →u (t )=K R
u(t)
KR
t
Pentru un sistem format dintr-un regulator şi o parte
fixată reprezentată printr-un element de întârzâiere de
ordinul I avem:
1
1 s
s
v
y
ε u KF
r
+
HR(s) Σ H F (s)=
--
T F s+ 1
y 37
Pentru bucla închisă, în raport cu referinţa r de tip
treaptă unitară avem:
H d (s ) KRK F 1 K0
H 0 ( s)= = =
H d (s )+ 1 K R K F + 1 TF T 0 s+ 1
s+ 1
K R K F+ 1
Sistemul H0 rămâne unul de ordinul I cu următoarele
performanţe:
1 K R KF TF
εst = iar K 0 = ; T 0=
K R K F+ 1 K R K F+ 1 K R K F+ 1
KF
T F s+ 1 KF KF 1
H 0v= = ; Y v ( s)=
KR T F s + 1+ K F K R T F s+ 1+ K F K R s
1+ K F
T F s+ 1
Y st v =lim [s Y v (s)]=lim s
¿
s→0 s→0 [ KF
]
1
T F s+ 1+ K F K R s
=
KF
K F K R+ 1
Creşterea factorului de amplificare KR determina o
reducere a erorii staţionare (deci o creştere a preciziei) şi
o reducere a constantei de timp a sistemului (viteza mai buna
de răspuns).
P
KR ↑ εst ↓ T0 ↓ K0 ↑ yst_v ↓
Cu cat KR creste, cu atât efectul perturbaţiei este scăzut
(răspunsul ieşirii în regim staţionar scade): Un regulator P se
poate alege atunci când procesul conţine cel puţin un element
integrator, astfel se asigura eroare staţionară nulă deci o
buna comportare a sistemului în regim staţionar. Pentru
procese cu mai multe constante de timp, alegerea unui
regulator P, poate duce la instabilitatea sistemului.
38
Regulatorul integrator (I)
t
Legea de reglare: u( t)= 1 ∫ ε( t)dt ,Ti - constanta de
Ti 0
integrare
1 1
H R ( s)= sau H R ( s)= dacă se consideră şi
Ti s T i s( τ s+ 1)
întârzâierea proprie a regulatorului.
u(t) ideal
KR
real
t
Răspunsul indicial
KR 1 KR t
Y (s )= = 2
→ y ( t)=
Tis s T is Ti
Abaterea răspunsului real indicial al unui regulator I, în
raport cu răspunsul ideal, este cu atât mai mare cu cat constanta
de timp proprie regulatorului este mai mare
Răspunsul unui regulator integrator este o rampa cu panta
1/Ti.
Regulatorul integrator are un caracter de memorie deoarece la
o comanda u(t) nenula poate fi trimisa spre proces chiar daca
intrarea în regulator este nula.
Prezenta polului în origine în funcţia de transfer a
regulatorului asigura o buna comportare în regim staţionar a SRA la
intrarea treapta unitara însa gradul de stabilitate al sistemului
poate să scadă.
Regulatoarele integratoare se folosesc în combinaţie cu cele
de tip P.
39
Regulatorul proporţional-integrator (PI)
( )
t
Legea de reglare: u(t)=K R ε(t)+ 1 ∫ ε(t)dt ,Ti - constanta de
Ti 0
(
H R ( s)=K R 1+
1
Tis ) sau H R ( s)=
(
K R 1+
1
T is ) dacă se consideră şi
τs+ 1
întârzâierea proprie a regulatorului (τ - constanta de timp
proprie a regulatorului).
Efectul integrator determina asigurarea preciziei
răspunsului (eroare staţionară zero) iar efectul proporţional
duce la creşterea vitezei de răspuns a SRA.
u(t)
ideal
2KR
real
KR
ε
1
Ti t
40
KR
Y (s )=
s
1
+ 2
s Ti
→ y (t)=K R 1+
( ) t
Ti
u
KF
( )
+ ε 1 y
r H R (s )=K R 1+
Tis Σ H F ( s)=
T F s+ 1
--
ω2n =
KR K F
TiT F
; ξ=
1
2 √ T i T F (1+ K R K F )
KR K F TF
41
Răspunsul sistemului este compus din răspunsul unui
sistem oscilant de ordin II la care se adăuga mărimi
proporţionale cu derivata acestuia, cu coeficient de
proporţionalitate egal cu Ti. Ca urmare a prezentei acţiunii
derivative, durata procesului tranzitoriu se reduce. Poziţia
zeroului influenţează puternic suprareglajul şi durata reg.
tranzitoriu.
Pe măsura apropierii zeroului de 0 suprareglajul creste
iar durata regimului tranzitoriu scade.
42
Regulatorul proporţional-derivativ (PD)
u
+ ε KF y
r H R (s )=K R ( 1+ T d s ) Σ H F ( s)=
T F s+ 1
--
K R ( 1+ T d s )
H R ( s)=K R ( 1+ T d s ) sau H R ( s)= dacă se consideră şi
τ s+ 1
întârzâierea proprie a regulatorului (τ - constanta de timp
proprie a regulatorului).
Răspunsul indicial al regulatorului PD
u(t)
ideal
real
KR
ε
1
43
Analiza unui SRA cu REG de tip PD şi parte fixată de
tip „elem. de întârziere de ordinul I”:
KF
H d (s )=K R (1+ T d s)
În raport cu referinţa: T F s+ 1=K R K F Astfel,
Se alegeT d =T F
constanta derivativă compensează constanta dominantă a întârzâierii
părţii fixate.
KRK F
H 0 ( s)=
1+ K R K F
1
εst = T 0 ( s)=0
KRKF
u
+ ε KF y
r H R (s )=K R ( 1+ T d s ) Σ H F (s)=
(T 1 s+ 1)(T 2 s+ 1)
--
y
KRK F
H 0 ( s)= 2
T 1 T 2 s + ( T 1+ T 2 ) s+ K R K F + 1
K R KF 1
K 0 (s)= ; ε st = Acestea sunt identice cu cele pentru
1+ K R K F 1+ K R K F
regulator P.
44
Răspunsul indicial al regulatorului PID
u(t) ideal
real
D P
I
KR
ε
1
t
1
1 s
s v
ε u
KF y
r
+
--
H R (s )=K R
( 1+
1
Tis
+ Td s
) Σ H F (s)=
(T 1 s+ 1)(T 2 s+ 1)
46
K R (Θ1 s+ 1)(Θ2 s+ 1)
H R ( s)=
K Θ s ( τ s+ 1)
K
K τ K RK F
H 0 ( s)= = , unde s−a notat K =
s( τ s + 1)+ K 1 K Θ
s 2+ τ s+ τ
47
ACORDAREA REGULATOARELOR PENTRU PROCESE RAPIDE
+ ε u y
r HR(s) HF1(s) Σ HF2(s)
--
Y R (s ) Y V (s)
⏞
Y ( s )=H 0 (s )R (s )+ ⏞
H 0v (s)V ( s )
48
{ }
∣H 0 ( j ω)∣=1
{ Y V ( s )=H 0V ( s)V (s)=0 }{
Y R ( s )=H 0 ( s) R(s)=R(s) → H 0 (s)=1 →
H 0V ( s )=0 }
∣H 0V ( j ω)∣=0 ptr.orice ω .
arg ( H 0 (s ))=0
KF
H IF ( s )= n
(1)
(T Σ s+ 1) ∏ (T i s+ 1)
i=1
sau
KF
H IIF ( s)= n (2)
s (T Σ s+ 1) ∏ (T i s+ 1)
i=1
nde:
KF - reprezintă coeficientul de transfer al parţii fixate
Ti – reprezintă constantele dominante (mari) de timp
(constante de timp principale):
Ti < 10
T∑ - reprezintă suma constantelor de timp parazite care sunt
mult mai mici decât cele predominante:
1
TS = min {Ti }
10
49
intr-unul de ordinul II în care se elimina contribuţia în
răspuns a constantelor de timp dominante. In ceea ce priveşte
anularea efectului perturbaţiilor, criteriul modulului propune
forme analitice de legi de reglare care plasează un element
integrator înaintea punctului de aplicaţie al acestora.
Prezenta elementului integrator înaintea punctului de
aplicaţie a perturbaţiilor conduce la anularea efectului
acestora daca ele sunt de tip treapta unitara.
1
Y V = H 0V ( s) →Y st =lim
s
1
s→0 s ( )
H 0V =lim H 0V=0
s→0
∏ (θk s+ 1)
{ }
m=n
I
H (s )=
R
1
unde θ k =T k ptr. varianta (1) a părţii
θs
θ=2K F T Σ
fixate,
şi
m
∏ (θ k+ 1)
{ }
m=n
II
H (s)=
R
1
unde θk =T k ptr. varianta (2) a părţii fixate,
θs
θ=2K F T Σ
Pentru cazul I:
I 1
H d (s )= iar funcţia de tranfer ptr. bucla închisă va
2TΣ s (T Σ s + 1)
fi:
1
I 1 2T2Σ ω 2n
H 0 ( s)= = =
1+ 2T2Σ s 2+ 2TΣ S s 2+ 1 s+ 1 s 2+ 2 ξ ωn s+ ω2n
TΣ 2T2Σ
50
{ }
1
ωn=
√ 2T Σ
ξ= √ =0,7
2
2
II s
Pentru cazul II: H d ( s)= 2 iar funcţia de tranfer ptr.
2TΣ s (T Σ s+ 1)
bucla închisă va fi:
1
II 1 2T2Σ ω 2n
H 0 (s)= 2 2
= = 2 2
1+ 2TΣ s + 2T Σ s s 2+ 1 s + 1 s + 2 ξ ωn s+ ωn
TΣ 2T 2Σ
{ }
1
ωn=
√ 2T Σ
√ 2
ξ= =0,7
2
4
tt ≈ ≈ 8TΣ
ξ ωn
51
PROIECTAREA SRA PRIN METODA ALOCĂRII POLI-ZEROURI
K 0 ω2n
H 0 ( s)=
s2 +2 ξ ωn s+ω2n
52
transpunerea performanţelor pentru un sistem de ordinul II cu
un pol suplimentar.
K 0 p 3 ω2n
H 0 ( s)= 2
(s + 2 ξ ωn s+ ω 2n)⋅( s+ p 3)
53
sistemului), (banda de frecventa sau lărgimea de banda ce
indica comportarea sistemului fata de perturbaţiile de
frecventa înaltă)
p1
Sistem cu doi poli:
ωn
φ
p2
A. Eroarea staţionara nulă la intrarea treapta unitara :
εst =0
e st = lim e ( t ) = lim éë s × e ( s ) ùû = lim éë s × ( R ( s ) - Y ( s ) ) ùû = lim éë s × (1 - H 0 ( s ) ) × R ( s ) ùû
t ®¥ s ®0 s ®0 s ®0
é 1ù
e st = lim ê s × (1 - H 0 ( s ) ) × ú = lim éë(1 - H 0 ( s ) ) ùû = 1 - H 0 ( s ) = 0
s ®0
ë s û s®0
H 0 (0)=1 pentru reacţie unitara.
H0 ( s) =
C
e st = 0 Þ H 0 ( s ) =
p1 p2 p1 p2
= 2
( s + p1 )( s + p2 ) (1 + p1 )(1 + p2 ) s + ( p1 + p2 ) s + p1 p2
Dacă dorim să punem in evidență elementele caracteristice ale
unui sistem oscilant de ordinul II:
K 0 ω2n
p 1,2=−ξ ωn √ 1−ξ
2
H 0 ( s)= 2 2
unde
s + 2 ξ ωn s+ω n
σ 100%
πξ
−
= e √ 1−ξ ⋅100
2
σ II
0
1 ξ 54
( ) ( ) ( )
Restricția privind timpul de tranziție va plasa polii in
stânga dreptei dată de ξ ωn
4 4
Ex: t t=10 sec → ≤10 → ξ ωn≥ → ξ ω n≥0.4
ξω n 10
wn2 2
H 0 ( jw B ) = =
2 →
(w - wB2 ) + 4xwn2wB2
2 2
n
√
ω B = ω n 1−2 ξ2 + √ 2−4 ξ 2 +4 ξ 4 ≤ ωB impus
56
( ( )) ( ( () ( ) ))
( ( ( )) ( )) ( ( ))
SISTEME CU DOI POLI ȘI CU UN ZERO
p1
ωn
Cp1 p2 ( s + z )
Sistemele au forma: H 03 =
( s + p1 )( s + p2 ) z φ
p2
1
H 03 (0)=1 →C = pentru reacţie unitara.
z
p1 p 2 ω2n
( s+ z) ( s+ z)
z , de asemenea z
H 03 ( s)= H 03 ( s)= 2 2
(s+ p 1)(s+ p 2) s +2 ξωn s+ωn
ω 2n
(s+ z )
H 03 ( s)= 2
z
2
= 2
ω2n
2
+
1
2
ω2n
2
s +2 ξ ωn s+ωn s + 2ξ ωn s+ω n z s +2 ξ ω n s+ωn
ceea ce corespunde
(
s
)
1 1 dy II
H 03 ( s)=Y ( s)+ s Y ( s)→ y 3 (t )= y II (t )+
z z dt
Oscilaţiile procesului tranzitoriu cresc ceea ce conduce
implicit la scăderea duratei regimului tranzitoriu. Aceasta
duce la modificarea la modificarea lărgimii de bandă.
Contracţia în domeniul timpului corespunde dilatării în
domeniul frecventelor şi invers.
ωn
B. Condiţia de suprareglaj: H 03 ( s)= f (λ) , λ=
∈[0,2]
z
Introducerea unui zero este justificata daca λ Î [0,2] . Daca
zeroul este plasat spre +¥ , efectul introducerii lui este
neglijabil.
dy 3
Pentru a determina suprareglajul se pune condiţia: =0
dt
58
−ξ π−( γ−ϕ)
1 1 1 2ξ 1
εv = + − = Rezulta ca pentru un sistem de ordinul II,
p 1 p2 z ω n z
prin introducerea unui zero suplimentar eroarea la viteza se
reduce comparativ cu cea a sistemului necompensat
C p1 p 2
Sistemele au forma: H 04 (s)=
(s+ p 1)(s+ p2 )(s+ p 3)
ω 2n p 3 p1 p 2 p3
H 04 (s)= 2 2
sau H 04 (s)=
(s +2 ξ ωn s+ω n)(s+ p3) (s+ p 1)(s+ p2 )(s+ p 3)
Raspunsul indicial:
p3
Y 4 ( s)=Y 2 (s)+ C → y 4 (t)= y II (t )+e− p t C
3
s + p3
In acest caz suprareglajul, timpul de tranzitie si eroarea la
viteza cresc. Datorita cresterii timpului de tranzitie se
reduce largimea de banda.
Polii suplimentari trebuie introduşi pentru realizabilitatea
fizica a sistemului astfel încât: ∣p 3∣≥(5 ⋯10) ξ ωn
59
ω2n p3 p1 p 2 p 3
(s+ z) (s+ z )
z sau z
H 05 ( s)= 2 2
H 05 ( s)=
(s +2 ξ ω n s+ωn )(s+ p 3) (s+ p 1)(s+ p 2)( s+ p3 )
1 1 1 1
εv = + + −
p 1 p2 p 3 z
60