L V V L U S y S: 2.1.3 Performanţele Sistemelor Automate Liniare Şi Continue
L V V L U S y S: 2.1.3 Performanţele Sistemelor Automate Liniare Şi Continue
L V V L U S y S: 2.1.3 Performanţele Sistemelor Automate Liniare Şi Continue
y (t ) = u (t − τ) , τ ≥ 0. (2.133)
Din ecuaţia (2.133) rezultă că răspunsul unui element cu timp mort la un semnal
treaptă unitară are aspectul din figura 2.25,b.
Exemple de elemente cu timp mort:
• Instalaţiile de transport de masă şi energie la care mărimea de ieşire rezultă numai
după un anumit timp τ, reprezintă elemente cu timp mort (v. fig. 2.26,a,b).
l
v S0 (u ) S1 ( y )
v
u1 (u ) u2 ( y )
Cap de Cap de l
înregistrare ştergere
(a) (b)
Fig.2.26 Exemple de elemente cu timp mort.
• Eroarea staţionară
Acest indice de calitate caracterizează precizia de funcţionare a sistemelor automate
în regim staţionar.
Considerând un sistem automat liniar, cu parametrii constanţi, se poate defini eroa-
rea (sau abaterea) staţionară ε st ca fiind:
ε st = u − y st , (2.134)
unde y st reprezintă valoarea staţionară (sau stabilizată) a mărimii de la ieşire:
y st = lim y (t ) , (2.135)
t →∞
48
SCOPUL ŞI METODELE ANALIZEI
ε st = 0 , (2.136)
iar pentru sistemele de tipul celui din figura 2.27,b se impune limitarea erorii staţionare la
o valoare impusă ε st imp (valoarea maximă admisibilă din considerente legate de buna
desfăşurare a procesului tehnologic automatizat); deci se impune performanţa:
ε st ≤ ε st imp , (2.137)
t
0
(b)
Fig.2.27 Răspunsul indicial al unui sistem de ordinul doi: ε st = 0 (a); ε st ≠ 0 (b).
• Suprareglajul
49
ANALIZA SISTEMELOR AUTOMATE LINIARE ŞI CONTINUE
σ = y1max − y st . (2.138)
y1max − y st
σ 00 = 100 . (2.139)
y st
σ[ 0 0 ] ≤ σ imp [ 0 0 ] , (2.140)
• Gradul de amortizare
σ3 σ
δ = 1− = 1− 3 , (2.141)
σ1 σ
50
SCOPUL ŞI METODELE ANALIZEI
Cu cât δ este mai mare (deci mai apropiat de 1) cu atât amortizarea oscilaţiilor din
timpul regimului tranzitoriu este mai rapidă, deci procesul tranzitoriu al sistemului auto-
mat este de mai bună calitate.
după cum s-a considerat şi în figura 2.27,a, deci o gamă de ± 5 0 0 în jurul valorii staţio-
nare y st .
În scopul asigurării rapidităţii necesare desfăşurării procesului de reglare automată,
pentru timpul de răspuns t r se impune o performanţă de forma:
t r ≤ t r imp , (2.145)
unde t r imp este valoarea maximă admisibilă impusă de specificul instalaţiei tehnologice
automatizate şi al procesului tehnologic respectiv.
u
y
y u
D
1
0,95 B
y st = u
0,05 A t
0 C tc D′
ts
tm
Fig.2.28 Răspunsul indicial al unui sistem de ordinul doi cu punerea în evidenţă a timpului de creştere ( t c ),
timpului primului maxim ( tm ) şi a timpului primei atingeri a valorii staţionare ( t s ).
– Al doilea procedeu: se trasează tangenta la curba y în punctul curbei care are or-
donata 0,05 y st sau în punctul de inflexiune al curbei, stabilindu-se punctele de intersecţie
ale tangentei cu axa absciselor şi cu orizontala dusă prin ordonata y = y st ; pentru t c se
adoptă o valoare egală cu proiecţia pe axa timpului a segmentului de tangentă cuprins
între cele două puncte de intersecţie.
Definirea timpului primului maxim t m şi a timpului primei atingeri a valorii staţio-
nare t s , rezultă din figura 2.28 reprezentată pentru cazul când ε st = 0 şi deci y st = u .
• Eroarea staţionară
ε st = 0 , (2.146)
ε st ≠ 0 . (2.147)
52
SCOPUL ŞI METODELE ANALIZEI
• Abaterea maximă
0,95 yst
u, y
p yst = u
1
t
0 tr
u (a)
y ν1 = ν
p
y
u, y u ν3 ε st ≠ 0
ν2
y st
p yst
1
t
0
(b)
Fig.2.29 Răspunsul unui sistem de ordinul doi la o variaţie în treaptă a perturbaţiei: ε st = 0 (a); ε st ≠ 0 (b).
unde ν imp este valoarea impusă. Pentru exprimarea în procente a condiţiei (2.148), aba-
terea maximă se raportează la y st .
• Gradul de amortizare
În cazul variaţiei unei perturbaţii p (v. fig. 2.29,a şi 2.29,b) gradul de amortizare ρ
se defineşte ca diferenţa dintre unitate şi raportul amplitudinilor a două semioscilaţii suc-
cesive de acelaşi semn ale răspunsului y la variaţia respectivă a perturbaţiei, considerate
în raport cu valoarea staţionară y st :
53
ANALIZA SISTEMELOR AUTOMATE LINIARE ŞI CONTINUE
ν3 ν
ρ = 1− = 1− 3 , (2.149)
ν1 ν
ρ ≥ ρ imp , (2.150)
Pentru regimul sinusoidal, respectiv regimul staţionar corespunzând aplicării unor sem-
nale sinusoidale de frecvenţă variabilă la intrarea sistemului, se impune ca performanţă
lărgimea de bandă, caracterizând proprietăţile de filtru ale sistemului, respectiv com-
portarea acestuia în raport cu perturbaţiile de înaltă frecvenţă.
Performanţa lărgimii de bandă va fi expusă în cadrul §2.4.2, după prezentarea ca-
racteristicilor de frecvenţă.
Pentru analiza sistemelor automate este folosită o varietate relativ largă de metode.
O primă metodă constă în stabilirea ecuaţiei diferenţiale care leagă mărimea de ie-
şire y şi derivatele ei (în raport cu timpul) de mărimea de intrare u şi derivatele ei (în ra-
port cu timpul) şi în rezolvarea acestei ecuaţii, soluţia reprezentând răspunsul y al siste-
mului la un semnal de intrare u dat. Având variaţia în timp a mărimii de ieşire y , toate
performanţele (staţionare şi tranzitorii) pot fi determinate, verificând astfel dacă sunt în-
deplinite condiţiile impuse acestor performanţe.
În mod analog pot fi stabilite ecuaţii în care intervin mărimea de ieşire, derivatele
ei, o perturbaţie şi derivatele ei, rezultând răspunsul y , la variaţia unei perturbaţii, răs-
puns care permite verificarea performanţelor răspunsului la perturbaţie.
Metoda devine greu utilizabilă pentru sistemele industriale complexe deoarece ne-
cesită rezolvarea unor ecuaţii diferenţiale de grad superior. În ultimul timp, însă calcu-
latoarele electronice îşi aduc un aport preţios în acest domeniu.
O a doua metodă, larg folosită, utilizează funcţiile de transfer, deci transformata
Laplace directă şi inversă, pentru analiza sistemelor automate; această metodă permite de
asemenea să se obţină răspunsul sistemului la o variaţie a mărimii de intrare sau la o
variaţie a unei perturbaţii, ceea ce permite verificarea îndeplinirii performanţelor impuse.
Utilizarea funcţiilor de transfer permite să se stabilească anumite corespondenţe
între domeniul variabilei complexe şi domeniul timpului, corespondenţe utile atât pentru
analiza, cât şi pentru sinteza sistemelor automate.
54
METODA INTEGRĂRII ECUAŢIILOR DIFERENŢIALE ÎN ANALIZA SISTEMELOR AUTOMATE
O serie de elemente (sisteme) sunt descrise de ecuaţii diferenţiale liniare de ordinul întâi.
În mod convenţional, aceste elemente (sisteme) se numesc de ordinul întâi, ordinul ecu-
aţiei diferenţiale determinând, în acest mod, ordinul sistemului automat.
Forma generală a ecuaţiei diferenţiale liniare pentru un element (sistem) de ordinul
întâi este:
dy
a0 + a1 y = b1u , (2.151)
dt
55
ANALIZA SISTEMELOR AUTOMATE LINIARE ŞI CONTINUE
dy
T + y = Ku , (2.152)
dt
a0
T= , (2.153)
a1
b1
K= . (2.154)
a1
y (t ) = y f (t ) + y l (t ) , (2.155)
y l (t ) = 0 , (2.156)
astfel că:
y (t ) = y f (t ) = y st = ct. , (2.157)
dy f
T + y f = Ku , (2.158)
dt
În conformitate cu relaţia (2.157), primul termen este nul şi ecuaţia (2.158) devine:
y f (t ) = Ku (t ) = K1+ (t ) . (2.159)
Tp + 1 = 0 , (2.160)
1
p=− . (2.161)
T
Pentru componenta liberă va rezulta o expresie de forma:
t
−
pt
yl (t ) = Ce = Ce T , (2.162)
y (t ) = 0 pentru t = 0 ,
C = −K . (2.164)
scurtă. e
T
0,63K
57
ANALIZA SISTEMELOR AUTOMATE LINIARE ŞI CONTINUE
Elementele (sistemele) a căror modele matematice sunt ecuaţii diferenţiale de ordinul doi,
de forma generală:
d2 y dy
a0 + a1 + a 2 y = b2 u , (2.166)
dt 2 dt
d2 y dy
T2 2
+ 2ξT + y = Ku , (2.167)
dt dt
sau:
d2 y dy
2
+ 2ξωn + ω n 2 y = Kω n 2 u , (2.168)
dt dt
a0 1 a1 b
T= = , ξ= , K= 2 . (2.169)
a2 ωn 2 a 2 a0 a2
J 1 D
T= = , ξ= , K = 1, (2.170)
K T ωn 2 KT J
Dacă se ţine seama că în noul regim staţionar, după anularea componentei libere, va
rezulta:
y f (t ) = y st (t ) = ct , (2.172)
u (t ) = ct. = 1 (2.173)
şi variaţia în timp a mărimii y (t ) trebuie să aibă, în noul regim staţionar, acelaşi aspect ca
şi variaţia în timp a mărimii u (t ) .
Considerând deci că regimul tranzitoriu s-a încheiat, adică:
y l (t ) = 0 (2.174)
şi:
y (t ) = y st (t ) = y f (t ) = ct. , (2.175)
ecuaţia (2.168) (în care în continuare se va lua K = 1 ) pentru noul regim staţionar se poa-
te scrie:
d2 y f dy f
2
+ 2ξωn + ωn 2 y f = ωn 2 u . (2.176)
dt dt
Având în vedere relaţia (2.172), primii doi termeni sunt nuli şi ecuaţia (2.176)
devine:
ωn 2 y f = ωn 2 u , (1.177)
y f (t ) = y st (t ) = u (t ) = 1 . (2.178)
ε st = u − y st = 0 . (2.179)
59
ANALIZA SISTEMELOR AUTOMATE LINIARE ŞI CONTINUE
p 2 + 2ξω n p + ω n 2 = 0 , (2.180)
În practică toţi factorii care intervin în expresiile (2.169) sunt pozitivi, de aceea fac-
torul de amortizare ξ nu poate lua valori negative, ci numai valorile:
ξ ≥ 0. (2.182)
p1, 2 = ± jω n , (2.183)
pentru:
0 < ξ < 1, (2.184)
p1, 2 = −ω n , (2.186)
y (t ) = 0 pentru t = 0 . (2.189)
dy
= 0 pentru t = 0 . (2.190)
dt
60
METODA INTEGRĂRII ECUAŢIILOR DIFERENŢIALE ÎN ANALIZA SISTEMELOR AUTOMATE
u
y
2
y b
y1max y
u σ a
1
y c
d yst = u = 1
y
t
0 tm
Fig.2.31 Răspunsul indicial al unui sistem de ordinul doi pentru diverse valori ale factorului de amortizare ξ.
1 + C1 + C 2 = 0 , (2.191)
respectiv:
C1 + C 2 = −1 . (2.192)
În mod analog aplicând a doua condiţie iniţială (2.190) la expresia (2.188) rezultă:
p1C1 + p 2 C 2 = 0 . (2.193)
−1 1 1 −1
0 p2 − p2 p 0 p1
C1 = = , C2 = 1 = . (2.194, 2.195)
1 1 p 2 − p1 1 1 p 2 − p1
p1 p2 p1 p2
În ipoteza (2.184), care este cea mai frecventă în practică, se poate verifica uşor, în-
locuind valorile p1 şi p 2 din (2.181) în relaţiile (2.194) şi (2.195), că C1 şi C 2 sunt mă-
rimi complex conjugate.
Conform cu (2.181), considerând:
p1 = −ξω n + jω n 1 − ξ 2 , (2.196)
iar:
p 2 = −ξω n − jω n 1 − ξ 2 (2.197)
61
ANALIZA SISTEMELOR AUTOMATE LINIARE ŞI CONTINUE
p1 = p 2 = ξ 2 ωn 2 + ω n 2 (1 − ξ 2 ) = ωn , (2.199)
ωn 1 − ξ 2
1 p1 Re sin ϕ = = 1 − ξ2 , (2.201)
p 2 − ωn ξω n ωn
0
respectiv:
p2 − ωn 1 − ξ 2 1 − ξ2
tgϕ = . (2.202)
0 ξ
0 1 p2
Fig.2.32 Reprezentarea rădăcinilor ecuaţiei
Înlocuind relaţiile (2.200) şi (2.201) în
caracteristice în planul complex (2.198), expresia răspunsului sistemului de-
pentru diverse valori ale lui ξ . vine:
e − ξωnt
y (t ) = 1 − (cos ϕ sin ωn 1 − ξ 2 t + sin ϕ cos ωn 1 − ξ 2 t ) =
1 − ξ2
(2.203)
e −ξωnt
= 1− sin(ωn 1 − ξ 2 t + ϕ),
2
1− ξ
e −ξωnt 1 − ξ2
y (t ) = 1 − sin(ωn 1 − ξ 2 t + arctg ). (2.204)
1− ξ 2 ξ
62
METODA INTEGRĂRII ECUAŢIILOR DIFERENŢIALE ÎN ANALIZA SISTEMELOR AUTOMATE
Factorul de amortizare ξ are o influenţă hotărâtoare asupra regimului tranzitoriu din punct
de vedere al amortizării (ceea ce a determinat şi denumirea acestui factor). În acest regim
prezintă interes, după cum s-a mai menţionat, doi indici de calitate: suprareglajul şi dura-
ta regimului tranzitoriu (în cazul în care cauza apariţiei regimului tranzitoriu a fost varia-
ţia semnalului de la intrare).
Factorul de amortizare ξ exercită deci o influenţă determinantă asupra supraregla-
jului, iar în cazul sistemului considerat (sistem automat de ordinul doi descris de ecuaţia
(2.167) sau (2.168)) rezultă că suprareglajul (v. fig. 2.31, curba a) este determinat numai
de factorul ξ .
Pentru ξ = 0 , răspunsul sistemului de ordinul doi la un semnal de intrare treaptă
unitară (2.204), are forma:
π
y (t ) = 1 − sin(ωn t + ) = 1 − cos ωn t , (2.205)
2
pentru ϕ = π 2 (v. fig. 2.32)
Regimul tranzitoriu în acest caz este denumit regim oscilatoriu neamortizat sau
întreţinut şi este reprezentat prin curba b din figura 2.31. Întrucât mărimea de ieşire y (t )
nu ajunge la o valoare staţionară constantă sistemul nu este stabil (de fapt, sistemul se
află la limita de stabilitate – v. subcap. 2.4).
Datorită faptului că răspunsul sistemului neamortizat are pulsaţia ωn , acest factor a
primit denumirea menţionată anterior de pulsaţie naturală a sistemului neamortizat.
Pentru 0 < ξ < 1 , răspunsul sistemului are expresia analitică (2.204) şi este reprezen-
tat în figura 2.31, curba a. Rezultă deci în acest caz un regim oscilatoriu subamortizat.
Pentru ξ = 1 , din expresia (2.204) rezultă o nedeterminare; prin înlăturarea acesteia
se obţine răspunsul (v. fig. 2.31, curba c):
y (t ) = 1 − e − ωnt (1 + ω n t ) , (2.206)
deci se obţine un regim tranzitoriu aperiodic, denumit regim aperiodic critic amortizat.
Dacă se consideră ωn = ct. , în cazul diferitelor ipoteze făcute mai sus asupra valori-
lor parametrului ξ în gama 0 ≤ ξ ≤ 1 , se constată că rădăcinile ecuaţiei caracteristice se
deplasează pe un semicerc de rază ωn (v. fig. 2.32).
În cazul ξ > 1 , din expresia (2.181) rezultă că rădăcinile ecuaţiei caracteristice sunt
reale şi diferite, iar componenta corespunzătoare regimului tranzitoriu va fi formată din
două exponenţiale de semne diferite. Răspunsul sistemului, prezentat grafic în figura
2.31, curba d, are expresia analitică:
e − ξωnt
y (t ) = 1 − sh (ωn ξ 2 − 1 t − ϕ) . (2.207)
ξ2 −1
63
ANALIZA SISTEMELOR AUTOMATE LINIARE ŞI CONTINUE
Regimul tranzitoriu este tot aperiodic, dar de durată mai lungă decât în cazul ξ = 1 ;
acest regim este denumit regim aperiodic supraamortizat.
Revenind la cazul (2.184), se poate determina dependenţa dintre suprareglajul σ şi
factorul de amortizare ξ . În acest scop, se anulează derivata dy dt , determinându-se toa-
te momentele de maxim şi de minim ale răspunsului y, printre care şi tm , timpul la care
apare suprareglajul σ .
Întrucât:
y (t m ) = y st + σ = 1 + σ = y 1max ,
rezultă:
σ = y (t m ) − 1. (2.208)
dy
=
dt
1
=− [ −ξωn e −ξωnt sin(ωn 1 − ξ 2 t + ϕ) + ωn 1 − ξ 2 e −ξωnt cos(ωn 1 − ξ 2 t + ϕ)] = 0,
1 − ξ2
2 1 − ξ2
tg (ωn 1 − ξ t + ϕ) = = tgϕ , (2.209)
ξ
ceea ce conduce la relaţia din care se obţin toate momentele de maxim şi de minim ale
răspunsului y:
ω n 1 − ξ 2 t = Kπ , cu K = 0,1,2,… (2.210)
π
tm = , (2.211)
2
ωn 1 − ξ
1 2π π
tm = = . (2.212)
2 ω 1 − ξ2 ω 1 − ξ2
n n
64
METODA INTEGRĂRII ECUAŢIILOR DIFERENŢIALE ÎN ANALIZA SISTEMELOR AUTOMATE
65
ANALIZA SISTEMELOR AUTOMATE LINIARE ŞI CONTINUE
Înlocuind relaţiile (2. 214) şi (2, 218) în expresia (2.141) se obţine gradul de amorti-
zare pentru sistemele de ordinul II:
3πξ
− 2 πξ
1−ξ 2 −
σ3 e 1−ξ 2
δ = 1− = 1− πξ
= 1− e . (2. 219)
σ −
1−ξ 2
e
Se constată că la fel ca şi suprareglajul σ , gradul de amortizare δ depinde numai de
ξ şi nu depinde de ωn .
Factorul de amortizare ξ exercită o anumită influenţă şi asupra duratei regimului
tranzitoriu t r , însă ţinând seama că în practică factorul ξ are o gamă limitată de variaţie
( 0,5 < ξ < 0,9 ) determinată de faptul că suprareglajele admisibile sunt cele corespun-
zătoare acestei game, influenţa factorului ξ asupra duratei t r este redusă.
În paragraful 2.1.3 s-a arătat că regimul tranzitoriu se poate considera încheiat atunci
când mărimea de ieşire intră într-o gamă de ± 5% din valoarea sa staţionară şi ulterior nu
o mai părăseşte. Aceasta înseamnă, în cazul sistemului automat analizat (sistem automat
de ordinul doi descris de ecuaţia diferenţială (2.168)), că regimul tranzitoriu sa va încheia
când componenta corespunzătoare regimului tranzitoriu yl (t ) = ytr (t ) a mărimii de ieşire
va intra în banda ± 0,05 y st şi nu o va mai părăsi.
Din expresiile (2.171) şi (2.203) rezultă:
e − ξωnt
y tr (t ) = − sin(ω n 1 − ξ 2 t + ϕ) (2.220)
2
1− ξ
şi deci, având în vedere şi relaţia (2.178), condiţia încheierii regimului tranzitoriu devine:
e − ξωnt
ytr (t ) = sin(ω n 1 − ξ 2 t + ϕ) ≤ 0,05 y st = 0,05 . (2.221)
1 − ξ2
Ţinând seama de faptul că valorile absolute ale sinusului sunt limitate la unitate,
condiţia (2.221) poate fi înlocuită prin alta acoperitoare:
e − ξωnt
≤ 0,05 , (2.222)
1 − ξ2
66
METODA INTEGRĂRII ECUAŢIILOR DIFERENŢIALE ÎN ANALIZA SISTEMELOR AUTOMATE
e − ξωntr
= 0,05 . (2.223)
2
1− ξ
− ξω n t r = ln(0,05 1 − ξ 2 ) , (2.224)
respectiv:
ln(0,05 1 − ξ 2 )
tr = . (2.225)
− ξω n
e − ξωnt1
y trt1 (t ) =− sin(ω n 1 − ξ 2 t1 + ϕ) , (2.226)
1 − ξ2
e − ξωnt2
y trt2 (t ) = − sin(ω n 1 − ξ 2 t 2 + ϕ) . (2.227)
1 − ξ2
Pentru simplificare se consideră valorile maxime ale funcţiei sinus, în acest caz se
obţine:
y trt1 (t )
= e ξωn (t2 −t1 ) , (2.228)
y trt2 (t )
de unde:
67
ANALIZA SISTEMELOR AUTOMATE LINIARE ŞI CONTINUE
1 y trt1 (t )
t 2 − t1 = ln . (2.229)
ξω n y trt2 (t )
1 1 3 1 1 4
tr = ln ≅ , sau t r = ln ≅ . (2.230)
ξω n 0,05 ξω n ξωn 0,02 ξωn
π
tr = tm = . (2.231)
ωn 1 − ξ 2
Practic realizarea valorii lui ξ din relaţiile (2.232) este costisitoare; din acest motiv
se admite că regimul tranzitoriu să se sfârşească după un număr dat de oscilaţii, de exem-
plu, t r = (1 ÷ 5)T (T reprezintă perioada oscilaţiei amortizate):
B
Im 2π
T= ,
A ω 1 − ξ2
n
° unde ω n 1 − ξ 2 reprezintă pulsaţia oscilaţiei a-
45
ϕ = 72 Re ξ1 = 0,44 , ϕ1 = 64° ,
0
respectiv:
Relaţia (2.230) este acoperitoare deoarece, pentru diferite valori ξ , conduce la va-
lori mai mari pentru t r , decât cele din tabelul 2.1, durata reală a regimului tranzitoriu
fiind astfel mai mică decât cea obţinută cu relaţia (2.230); aceasta înseamnă că dacă du-
rata t r , rezultată din relaţia aproximativă (2.230) va satisface performanţa (2.145), cu atât
mai mult durata reală t r (mai mică decât cea obţinută prin calcul) va satisface performan-
ţa menţionată.
De altfel, după cum s-a specificat cu ocazia relaţiei (2.222) şi datele din tabelul 2.1
sunt acoperitoare.
Transformarea directă Laplace permite ca studiul unei funcţii f (t ) să fie trecut din do-
meniul real al timpului (domeniul variabilei reale) în domeniul variabilei complexe s =
= σ + jω .
Prin aplicarea transformatelor Laplace, în condiţii iniţiale nule, ecuaţiile diferenţiale
care descriu funcţionarea elementelor sistemelor automate devin ecuaţii algebrice, cu po-
linoame de variabilă complexă s .
Fie f (t ) o funcţie măsurabilă reală sau complexă, definită pe dreapta reală şi s un
parametru real sau complex.
Definiţia 1. Dacă:
+∞
− st
∫ f (t )e dt , (2.234)
−∞
Se notează:
∞
F ( s ) = ‹[ f (t )] = ∫ f (t )e − st dt . (2.236)
0
69