81 3.4. Circuite de ordinul doi 3.4.1. Introducere Circuitele ... - derivat
81 3.4. Circuite de ordinul doi 3.4.1. Introducere Circuitele ... - derivat
81 3.4. Circuite de ordinul doi 3.4.1. Introducere Circuitele ... - derivat
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>3.4.</strong> <strong>Circuite</strong> <strong>de</strong> <strong>ordinul</strong> <strong>doi</strong><br />
<strong>3.4.</strong>1. <strong>Introducere</strong><br />
<strong>Circuite</strong>le care contin doua elemente dinamice ( doua con<strong>de</strong>nsatoare, doua bobine sau un<br />
con<strong>de</strong>nsator si o bobina) se numesc circuite <strong>de</strong> <strong>ordinul</strong> <strong>doi</strong>. Un circuit liniar <strong>de</strong> <strong>ordinul</strong> <strong>doi</strong> contine<br />
rezistoare liniare, elemente dinamice liniare, surse comandate liniar si surse in<strong>de</strong>pen<strong>de</strong>nte.<br />
x<br />
x =<br />
x<br />
⎡ 1 ⎤<br />
⎣⎢<br />
2 ⎦⎥<br />
Un astfel <strong>de</strong> circuit poate fi caracterizat prin ecuatia <strong>de</strong> stare x = Ax + u() t cu<br />
vectorul var iabilelor <strong>de</strong> stare ( u<br />
C<br />
pentru con<strong>de</strong>nsatoare si i<br />
L<br />
pentru bobine)<br />
u () t<br />
ut () = vectorul marimilor <strong>de</strong> int rare<br />
u () t<br />
⎡ 1 ⎤<br />
−<br />
⎣⎢<br />
2 ⎦⎥<br />
<strong>81</strong><br />
a a<br />
A = matricea <strong>de</strong> stare<br />
a a<br />
⎡ 11 12 ⎤<br />
−<br />
⎣⎢<br />
21 22 ⎦⎥<br />
Ecuatiei <strong>de</strong> stare x = Ax + u() t ii corespund ecuatii diferentiale scalare <strong>de</strong> <strong>ordinul</strong> <strong>doi</strong> pentru x1si x 2 :<br />
.. .<br />
x = T x − ∆x<br />
+ u ( )<br />
0<br />
1 1 1 a t cuconditiaa<br />
≠<br />
12<br />
.. .<br />
x = T x − ∆x<br />
+ u ( t)<br />
cuconditiaa<br />
≠ 0<br />
2 2 2 b<br />
21<br />
un<strong>de</strong>: T=a 11 +a 22 se numeste urma matricei A, ∆=a 11 a 22- a 12 a 21 este <strong>de</strong>terminantul matricei A si<br />
.<br />
ua () t =− a<br />
22<br />
u<br />
1<br />
() t + a<br />
12<br />
u<br />
2<br />
() t + u () t<br />
1<br />
.<br />
u () t = a u () t − a u () t + u () t<br />
b 22 1 11 2 2<br />
Daca a 12 =a 21 =0 ecuatia <strong>de</strong> stare se reduce la doua ecuatii diferentiale <strong>de</strong> <strong>ordinul</strong> intai. Pentru a<br />
<strong>de</strong>duce ecuatiile <strong>de</strong> <strong>ordinul</strong> II se consi<strong>de</strong>ra a12≠0, a21≠0 si se <strong>de</strong>riveaza in raport cu timpul ecuatia<br />
<strong>de</strong> stare<br />
.. . . .<br />
x 1 = a x + a x + u1(<br />
t)<br />
11 1 12 2<br />
Se inlocuiesc x .<br />
si x. cu expresiile date <strong>de</strong> ecuatia <strong>de</strong> stare si rezulta:<br />
1 2<br />
..<br />
.<br />
x a ( a x a x u ( t)) a ( a x a x u ( t)) u () t<br />
1 11 11 1 12 2 1 12 21 1 22 2 2 1<br />
..<br />
.<br />
x1( a<br />
11<br />
a a ) x ( a a a a ) x a u ( t) a u ( t) u () t<br />
dar a x x a x u () t<br />
2<br />
12 21 1 11 12 12 22 2 11 1 12 2 1 1<br />
= + + + + + +<br />
= + + + + + +<br />
12 2<br />
= −<br />
1 11 1<br />
−<br />
1<br />
..<br />
x ( a a a ) x ( a a )( x a x u t a u t a u t u t<br />
x a a x a a a a x a u t a u t u t<br />
.<br />
.<br />
()) () () ()<br />
..<br />
( ) .<br />
1 11<br />
.<br />
( ) () () ()<br />
2<br />
= +<br />
12 21 1<br />
+<br />
11<br />
+<br />
22<br />
1−<br />
11 1<br />
−<br />
1<br />
+<br />
11 1<br />
+<br />
12 2<br />
+<br />
1<br />
=<br />
1 11<br />
+<br />
22<br />
−<br />
1 11 22<br />
−<br />
12 21 1<br />
−<br />
22 1<br />
+<br />
12 2<br />
+<br />
1
In mod similar se obtine si ecuatia <strong>de</strong> <strong>ordinul</strong> <strong>doi</strong> pentru x 2 .<br />
<strong>3.4.</strong>2. Scrierea ecuatiilor <strong>de</strong> stare ale circuitelor liniare<br />
Orice circuit liniar invariant in timp <strong>de</strong> <strong>ordinul</strong> <strong>doi</strong> poate fi consi<strong>de</strong>rat ca un cuadripol<br />
diport rezistiv liniar N (care contine rezistoare liniare si surse in<strong>de</strong>pen<strong>de</strong>nte) cu elementele<br />
dinamice conectate la porti. Se urmareste scrierea unui sistem <strong>de</strong> doua ecuatii diferentiale <strong>de</strong><br />
<strong>ordinul</strong> intai avand ca necunoscute variabilele <strong>de</strong> stare (tensiunile con<strong>de</strong>nsatoarelor si curentii prin<br />
bobine).<br />
Circuitul cu doua con<strong>de</strong>nsatoare<br />
Pentru fiecare con<strong>de</strong>nsator se poate scrie<br />
. i . i<br />
u =<br />
1<br />
u<br />
2<br />
1 C 2<br />
1<br />
C<br />
2<br />
−<br />
=−<br />
Daca N are o solutie unica pentru orice u 1 , u 2<br />
atunci N are o reprezentare controlata in tensiune (vezi paragraful 2.4.3.3.):<br />
⎡i<br />
g g u is t<br />
1⎤<br />
=<br />
11 12 1 1<br />
i g g u is t<br />
⎣⎢<br />
2⎦⎥<br />
21 22 2 2<br />
⎡ ⎤<br />
⎣⎢<br />
⎦⎥ ⎡ ⎤<br />
+<br />
⎣⎢<br />
⎦⎥ ⎡ () ⎤<br />
⎢⎢<br />
⎥⎥ ()<br />
⎣ ⎦<br />
Se exprima i1 si i2 in functie <strong>de</strong> u 1,<br />
C1 si , u 2 , C2 si se obtine:<br />
g g<br />
is () t<br />
11 12<br />
1<br />
⎡<br />
.<br />
⎤<br />
u C C u C<br />
⎢⎢<br />
1⎥⎥<br />
1 1 1 1<br />
. g g u i<br />
u<br />
s t<br />
2 21 22<br />
()<br />
⎣ ⎦<br />
2 2<br />
C<br />
2<br />
C<br />
2<br />
C<br />
2<br />
=<br />
⎡ − − ⎤ ⎡ − ⎤<br />
⎢⎢⎢⎢⎢<br />
⎥⎥⎥⎥⎥<br />
⎢⎢⎢⎢⎢<br />
⎥⎥⎥⎥⎥<br />
⎡ ⎤<br />
+<br />
− −<br />
−<br />
⎣⎢<br />
⎦⎥<br />
⎣<br />
⎦ ⎣ ⎦<br />
ceeace reprezinta ecuatiile <strong>de</strong> stare pentru acest circuit.<br />
Circuitul cu doua bobine<br />
82
di<br />
di<br />
Pentru fiecare bobina se scrie: u =− L<br />
1<br />
, u =−L<br />
2<br />
1 1 dt 2 2 dt<br />
Daca N are o solutie unica pentru orice i1 , i2 , atunci N are o reprezentare controlata in curent<br />
(vezi<br />
paragraful 2.4.3.3.) si ecuatia <strong>de</strong> stare este:<br />
⎡ −r r<br />
us t<br />
11<br />
−<br />
12 ⎤ ⎡ − () ⎤<br />
1<br />
⎡ . ⎤ ⎢<br />
⎥ ⎢ ⎥<br />
⎢i<br />
L L i L<br />
1 ⎥ 1 1 1 1<br />
.<br />
=<br />
⎢<br />
⎥⎡<br />
⎤ ⎢ ⎥<br />
⎢ ⎥ ⎢−r<br />
r i u<br />
i 21<br />
− ⎥⎢<br />
⎥ + ⎢<br />
22<br />
− s () t ⎥<br />
2<br />
⎣⎢<br />
2 ⎦⎥<br />
⎢<br />
⎥⎣<br />
⎦ ⎢ 2 ⎥<br />
⎣⎢<br />
L<br />
2<br />
L<br />
2 ⎦⎥<br />
⎢<br />
⎣<br />
L<br />
⎥<br />
2 ⎦<br />
Circuitul cu o bobina si un con<strong>de</strong>nsator<br />
Tensiunea la bornele bobinei este u<br />
2<br />
=− Li<br />
2<br />
.<br />
si curentul prin con<strong>de</strong>nsator este i<br />
1<br />
=− Cu<br />
1<br />
.<br />
.<br />
Daca N are o solutie unica pentru orice u1 si i2 atunci N are o reprezentare hibrida corespunzatoare<br />
(vezi capitolul 2) si ecuatia <strong>de</strong> stare este:<br />
⎡ −h h<br />
is () t<br />
. 11<br />
−<br />
12 ⎤ ⎡ − ⎤<br />
1<br />
⎡<br />
u<br />
⎤ ⎢<br />
⎥ ⎢ ⎥<br />
⎢ C C u<br />
1⎥<br />
1 1 1<br />
C<br />
1<br />
⎢.<br />
⎥<br />
=<br />
⎢<br />
⎥⎡<br />
⎤ ⎢ ⎥<br />
⎢−h<br />
h i u<br />
i 21<br />
− ⎥⎢<br />
⎥ + ⎢<br />
22<br />
− s () t ⎥<br />
2<br />
⎣⎢<br />
2 ⎦⎥<br />
⎢<br />
⎥⎣<br />
⎦ ⎢ 2 ⎥<br />
⎣⎢<br />
L<br />
2<br />
L<br />
2 ⎦⎥<br />
⎢<br />
⎣<br />
L<br />
⎥<br />
2 ⎦<br />
<strong>3.4.</strong>3. Raspunsul unui circuit liniar la excitatie nula<br />
In cazul excitatiei nule u(t)=0 ecuatia <strong>de</strong> stare <strong>de</strong>vine<br />
.<br />
x = Ax sau<br />
⎡<br />
.<br />
x<br />
⎤<br />
⎢ ⎥ a a x<br />
⎢ . ⎥<br />
=<br />
a a x<br />
x<br />
⎣⎢<br />
⎦⎥<br />
⎡ ⎤<br />
⎢ ⎥<br />
⎣ ⎦<br />
⎡<br />
1 11 12 1<br />
⎤<br />
⎢ ⎥<br />
21 22 ⎣ 2<br />
⎦<br />
2<br />
83
Primul pas pentru rezolvarea ecuatiei <strong>de</strong> mai sus este <strong>de</strong>terminarea valorilor proprii λ1 si λ2 ale<br />
matricei A. λ 1 si λ 2 sunt solutiile ecuatiei.<br />
⎡a<br />
<strong>de</strong>t<br />
11<br />
− λ a<br />
12 ⎤<br />
= λ<br />
2<br />
− ( a a ) ( a a a a )<br />
2<br />
T<br />
a<br />
21<br />
a<br />
11<br />
+<br />
22<br />
λ +<br />
11 22<br />
−<br />
12 21<br />
= λ − λ + ∆ = 0<br />
⎣⎢<br />
22<br />
− λ<br />
⎦⎥<br />
s1,<br />
2<br />
T<br />
= ±<br />
2<br />
2<br />
T<br />
− ∆<br />
4<br />
1<br />
s1 si s2 sunt frecventele naturale ale circuitului (vezi capitolul 7). Daca ∆≠ T<br />
2<br />
atunci s1≠s2 4<br />
sunt sau numere reale sau numere complex conjugate. Pasul urmator in rezolvarea ecuatiei <strong>de</strong> stare<br />
este <strong>de</strong>terminarea unor vectori proprii η1 si η2 un<strong>de</strong><br />
η<br />
η<br />
1<br />
=<br />
11<br />
η<br />
12<br />
⎡ ⎤<br />
⎣⎢<br />
⎦⎥<br />
si<br />
η<br />
η<br />
2<br />
=<br />
21<br />
η<br />
22<br />
⎡ ⎤<br />
⎣⎢<br />
⎦⎥<br />
Prin <strong>de</strong>finitie η1 si η2 sunt <strong>doi</strong> vectori nenuli care satisfac relatiile: Aη 1<br />
= s<br />
1<br />
η<br />
1<br />
siAη 2<br />
= s<br />
2<br />
η<br />
2<br />
.<br />
Daca se cunosc valorile proprii s1≠s2 si vectorii asociati lor η1 si η2 solutia ecuatiei <strong>de</strong> stare poate<br />
fi scrisa xt k e st s t<br />
() = (<br />
1<br />
) + ( k e )<br />
1<br />
1 2 2<br />
η η<br />
2<br />
un<strong>de</strong> k1 si k2 sunt constante arbitrare <strong>de</strong>terminate <strong>de</strong><br />
conditiile initiale x(0). Intr-a<strong>de</strong>var<br />
xt ke st s k e s t s k e st A k e s t .<br />
( ) = (<br />
1<br />
) + ( ) = ( ) + ( ) A = Ax.<br />
1<br />
11 2 2<br />
2 2 1 1<br />
1 2 2<br />
η η η η<br />
2<br />
<strong>3.4.</strong>4. Comportarea calitativa a unui circuit liniar cu excitatie nula<br />
Solutia x(t) are doua componente: x1(t) si x2(t). Evolutia circuitului plecand <strong>de</strong> la o stare<br />
initiala (x1(0), x2(0)) poate fi reprezentata printr-o curba in planul <strong>de</strong> coordonate x1, x2 (planul<br />
fazelor). Aceasta curba se numeste traiectorie si se obtine eliminand timpul din expresiile lui x1(t)<br />
si x2(t).<br />
84
Evolutia circuitului corespunzatoare mai multor stari initiale poate fi reprezentata printr-o<br />
multime <strong>de</strong> traiectorii in planul fazelor. Aceste traiectorii formeaza un portret <strong>de</strong> faza.<br />
Starea <strong>de</strong> echilibru este o stare initiala care ramane nemodificata in cursul evolutiei<br />
circuitului. Aceasta stare corespun<strong>de</strong> unui punct <strong>de</strong> echilibru x1Q, x2Q din planul fazelor astfel<br />
incat daca x1(0)=x1Q si x2(0)=x2Q atunci x1(t)=x1Q si x2(t)=x2Q. In consecinta, in punctul <strong>de</strong><br />
echilibru avem x = 0 si x = 0.<br />
1 2<br />
Punctele <strong>de</strong> echilibru se pot <strong>de</strong>termina rezolvand sistemul <strong>de</strong> ecuatii Ax=0 adica<br />
⎡a11<br />
a12⎤⎡x1Q⎤<br />
⎢<br />
0<br />
⎣a21<br />
a<br />
⎥⎢<br />
22 ⎦⎣x<br />
⎥ =<br />
2Q<br />
⎦<br />
Daca ∆≠0 acest sistem admite numai solutia banala si originea este singurul punct <strong>de</strong> echilibru<br />
adica x1Q = x2Q =0. Daca ∆=0 atunci avem o infinitate <strong>de</strong> puncte <strong>de</strong> echilibru care satisfac ecuatia<br />
a11 x1Q +a12 x2Q =0.<br />
In continuare vom arata ca doua traiectorii nu se pot intersecta intre ele intr-un punct care<br />
nu este punct <strong>de</strong> echilibru. Fie un circuit neliniar <strong>de</strong> <strong>ordinul</strong> <strong>doi</strong> avand ecuatiile <strong>de</strong> stare<br />
xC = f ( x , x ) si xC = f ( x , x ) . Panta tangentei la traiectorie<br />
1 1 1 2 2 2 1 2<br />
dx<br />
dx<br />
2<br />
1<br />
poate fi calculata ca<br />
x<br />
x<br />
2<br />
1<br />
f2( x1, x2)<br />
= ceea ce reprezinta o valoare unica intr-un punct dat.<br />
f ( x , x )<br />
1 1 2<br />
Intr-un punct <strong>de</strong> intersectie a doua traiectorii ar exista, in mod evi<strong>de</strong>nt, doua pante, <strong>de</strong>ci<br />
traiectoriile nu se pot intersecta. De la aceasta regula fac exceptie punctele <strong>de</strong> echilibru in care<br />
f1( x1, x2) = 0 si f2( x1, x2)<br />
= 0 <strong>de</strong>ci dx2<br />
dx<br />
0<br />
= (ne<strong>de</strong>terminat) si <strong>de</strong>ci pot exista mai multe pante.<br />
0<br />
1<br />
Asa cum se va ve<strong>de</strong>a in continuare portretul <strong>de</strong> faza, in care este reprezentata evolutia circuitului<br />
pornind din orice stare initiala, constituie o imagine sintetica a comportarii calitative a circuitului.<br />
s t s t<br />
Fie xt ()= ke<br />
1<br />
+ k e 1<br />
1 2 2 η η<br />
2<br />
solutia ecuatiei <strong>de</strong> stare Cx = Ax cu valorile proprii s1 si s2<br />
si vectorii proprii η1 si η2. Valorile proprii <strong>de</strong>termina comportarea calitativa a circuitului. In<br />
continuare se discuta toate cazurile posibile pentru s1 si s2.<br />
85
Cazul 1. Matricea A are valori proprii reale si distincte respectiv T<br />
4<br />
T T<br />
2<br />
T T<br />
2<br />
s<br />
1<br />
= + − ∆ si s<br />
2 4 2<br />
= − −∆<br />
. Exista trei posibilitati:<br />
2 4<br />
86<br />
2<br />
> ∆ si ∆ ≠ 0 si<br />
s t s t<br />
a) s2 0 Cu un rationament asemanator rezulta ca daca t →+∞, k<br />
1<br />
e 1<br />
dominant si traiectoriile tind catre ∞ si sunt paralele cu η 1 si daca t k e<br />
η<br />
1<br />
<strong>de</strong>vine<br />
s t<br />
→−∞, 2<br />
2<br />
η<br />
2<br />
<strong>de</strong>vine<br />
dominant si traiectoriile pleaca din origine tangente la η 2 . In acest caz se spune <strong>de</strong>spre origine ca<br />
este un nod instabil.
s t<br />
c) s2
a) pentru α=0 solutia ecuatiei <strong>de</strong> stare <strong>de</strong>vine<br />
r<br />
xt () = k( rcos( + t) − isin( + t) cu r=<br />
si i<br />
r<br />
⎡ ⎤<br />
⎢ ⎥ =<br />
⎣ ⎦<br />
⎡<br />
η 1 η<br />
2 η ϕ ω η ϕ ω η<br />
η<br />
η<br />
⎢<br />
2 ⎣η<br />
si portretul <strong>de</strong> faza este o elipsa care are pe η<br />
i<br />
si η r drept directii conjugate. In acest caz<br />
starea <strong>de</strong> echilibru x=0 se numeste centru si corespun<strong>de</strong> raspunsului fara pier<strong>de</strong>ri. Cele doua<br />
variabile <strong>de</strong> stare vor fi: x () t = A cos( ω t + θ ) , x () t = A cos( ω t + θ ) , un<strong>de</strong> A , A , θ , θ<br />
sunt constante.<br />
1 1 1<br />
88<br />
2 2 2<br />
i1<br />
i2<br />
⎤<br />
⎥<br />
⎦<br />
1 2 1 2<br />
b)daca α 0 x(t) are aceeasi expresie ca la punctul b iar portretul <strong>de</strong> faza contine spirale<br />
logaritmice care tind spre infinit cand t →∞.<br />
In acest caz, starea <strong>de</strong> echilibru x=0 este un focar instabil.
<strong>3.4.</strong>5. Ecuatiile <strong>de</strong> stare ale circuitelor neliniare<br />
Ecuatiile <strong>de</strong> stare pentru un circuit neliniar <strong>de</strong> <strong>ordinul</strong> <strong>doi</strong> scrise sub forma<br />
⎧ x1<br />
= f1(<br />
x1,<br />
x2<br />
, t)<br />
x = f ( x,<br />
t)<br />
sau⎨<br />
se numesc ecuatii <strong>de</strong> stare in forma normala.<br />
⎩x2<br />
= f 2 ( x1,<br />
x2<br />
, t)<br />
Scrierea ecuatiilor <strong>de</strong> stare in forma normala rezulta pentru orice circuit neliniar corect mo<strong>de</strong>lat.<br />
Multe meto<strong>de</strong> numerice <strong>de</strong> rezolvare a ecuatiilor diferentiale neliniare sunt formulate pentru forma<br />
normala a ecuatiilor.<br />
In cazul in care circuitul contine numai elemente invariante in timp si surse in<strong>de</strong>pen<strong>de</strong>nte<br />
<strong>de</strong> curent continuu, variabila timp nu mai apare explicit in ecuatii si acestea sunt <strong>de</strong> forma<br />
⎧ x1<br />
= f1(<br />
x1,<br />
x2<br />
)<br />
x = f ( x)<br />
sau⎨<br />
si se numesc ecuatii <strong>de</strong> stare autonome, iar circuitul se numeste<br />
⎩x2<br />
= f 2 ( x1,<br />
x2<br />
)<br />
circuit autonom. Un circuit in care parametrii surselor in<strong>de</strong>pen<strong>de</strong>nte sunt functii neconstante <strong>de</strong><br />
timp se numeste circuit neautonom.<br />
neliniare:<br />
Pentru scrierea ecuatiilor <strong>de</strong> stare se iau in consi<strong>de</strong>rare cele trei configuratii posibile:<br />
Alegerea variabilelor <strong>de</strong> stare se face in functie <strong>de</strong> tipurile con<strong>de</strong>nsatoarelor si bobinelor<br />
- pentru con<strong>de</strong>nsatoare controlate in tensiune, variabila <strong>de</strong> stare este tensiunea;<br />
- pentru con<strong>de</strong>nsatoare controlate in sarcina, variabila <strong>de</strong> stare este sarcina;<br />
- pentru bobinele controlate in curent, variabila <strong>de</strong> stare este curentul;<br />
- pentru bobinele controlate in flux, fluxul magnetic este variabila <strong>de</strong> stare.<br />
In continuare se proce<strong>de</strong>aza la fel ca in cazul circuitelor liniare. Daca presupunem ca diportul<br />
rezistiv N are o solutie unica pentru orice valori ale parametrilor surselor in<strong>de</strong>pen<strong>de</strong>nte <strong>de</strong> la porti<br />
alese corespunzator, atunci exista urmatoarele reprezentari ale lui N:<br />
- reprezentarea controlata in tensiune (circuitul cu doua con<strong>de</strong>nsatoare)<br />
i1 = i1( u1, u2, t)<br />
i = i ( u , u , t)<br />
2 2 1 2<br />
89
- reprezentarea controlata in curent (circuitul cu doua bobine)<br />
u1 = u1( i1, i2, t)<br />
u = u ( i , i , t)<br />
2 2 1 2<br />
- o reprezentare hibrida (circuitul cu con<strong>de</strong>nsator si bobina)<br />
i1 = i1( u1, i2, t)<br />
u = u ( u , i , t)<br />
2 2 1 2<br />
care se utilizeaza la scrierea ecuatiilor <strong>de</strong> stare. Iata cateva exemple:<br />
i) Fie un circuit autonom cu doua con<strong>de</strong>nsatoare controlate in tensiune avand ecuatiile<br />
constitutive q = q< ( u ) si q = q< ( u ). Exprimam pe i si i in functie <strong>de</strong> variabilele <strong>de</strong> stare<br />
u si u<br />
1 2 :<br />
i<br />
i<br />
1<br />
2<br />
1 2 1 2 2 2<br />
90<br />
1 2<br />
dq1<br />
dq1<br />
dq1( u1)<br />
=− =− uC 1 =− C1( u1)C u1 cu C1( u1)<br />
=<br />
dt du du<br />
1<br />
dq2<br />
dq2<br />
dq2( u2)<br />
=− =− uC 2 =− C2( u2)C u2 cu C2( u2)<br />
=<br />
dt du du<br />
2<br />
1<br />
si ecuatiile <strong>de</strong> stare sunt: C<br />
< 1<br />
u<br />
( , ), C<br />
<<br />
1<br />
=− i u u u<br />
i ( u , u )<br />
C ( u ) 1 1 2 2<br />
=−<br />
1 1<br />
C<br />
2<br />
( u<br />
2<br />
) 2 1 2<br />
ii) Fie un circuit autonom cu doua con<strong>de</strong>nsatoare controlate in sarcina avand ecuatiile<br />
constitutive u = u< ( q ) si u = u< ( q ). Variabilele <strong>de</strong> stare fiind q1 si q2 rezulta<br />
1 1 1 2 2 2<br />
i<br />
1<br />
=− q<br />
1<br />
, i<br />
2<br />
=−q<br />
2<br />
si ecuatiile <strong>de</strong> stare sunt:<br />
qC =− i< ( u ( q ), u ( q )) sau qC=− ~<br />
i ( q , q ) si qC=− ~<br />
i ( q , q ).<br />
1 1 1 1 2 2 1 1 1 2 2 2 1 2<br />
iii) Fie un circuit autonom cu o bobina controlata in curent <strong>de</strong> ecuatie constitutiva<br />
φ = φ<<br />
( i ) si o bobina controlata in flux <strong>de</strong> ecuatie constitutiva i = i ( φ ). Exprimam pe<br />
1 1 1<br />
u si u<br />
1 2<br />
1<br />
2<br />
2 2 2<br />
in functie <strong>de</strong> variabilele <strong>de</strong> stare i1 si φ 2 tinand seama ca tensiunea asociata dupa<br />
regula <strong>de</strong> la receptoare cu curentul din bobina este u<br />
C<br />
( ) < ( ,<<br />
1 1<br />
i<br />
( )) ~<br />
1 =− u1 i1 i2<br />
φ 2 =− u1( i1,<br />
φ 2)<br />
L1 i1 L1( i1) C φ =− u< ( i , i ( φ )) =−u~<br />
( i , φ )<br />
2 2 1 2 2 2 1 2<br />
Celelalte cazuri se trateaza similar.<br />
d<br />
=<br />
dt<br />
φ .<br />
<strong>3.4.</strong>6. Comportarea calitativa a unui circuit neliniar in jurul unui punct <strong>de</strong> echilibru<br />
Starea <strong>de</strong> echilibru este data <strong>de</strong> solutiile ecuatiilor f ( x , x ) = 0 si f ( x , x ) = 0.<br />
1 1 2 2 1 2<br />
Pentru circuite autonome neliniare curbele corespunzatoare celor doua ecuatii se pot intersecta in
mai multe puncte (puncte <strong>de</strong> echilibru) si <strong>de</strong>ci exista mai multe stari <strong>de</strong> echilibru. In continuare se<br />
studiaza com-<br />
portarea unui circuit neliniar autonom in jurul unui punct <strong>de</strong> echilibru Qx (<br />
1Q , x ).<br />
2Q<br />
Conceptele <strong>de</strong> traiectorie si portret <strong>de</strong> faza introduse la studiul circuitelor liniare <strong>de</strong> <strong>ordinul</strong><br />
<strong>doi</strong> pot fi utilizate si pentru circuitele neliniare. Traiectoria este curba care se obtine in planul<br />
x1 − x2<br />
eliminand pe t din expresiile x1() t si x2(). t Acesata curba se poate vizualiza aplicand<br />
semnalele x1() t si x2() t pe placile <strong>de</strong> <strong>de</strong>flexie pe verticala si pe orizontala ale unui osciloscop.<br />
Portretul <strong>de</strong> faza ofera informatii <strong>de</strong>spre comportarea calitativa a circuitului; <strong>de</strong> exemplu, o<br />
traiectorie inchisa inseamna ca circuitul oscileaza si o spirala care converge spre punctul <strong>de</strong><br />
echilibru inseamna ca oscilatiile sunt amortizate.<br />
Ecuatiile <strong>de</strong> stare ale circuitului sunt: x = f ( x) sau<br />
⎧x1<br />
= f1( x1, x2)<br />
⎨<br />
⎩x2<br />
= f ( x1, x2)<br />
Se <strong>de</strong>zvolta f1 si f2<br />
in serie Taylor in jurul punctului <strong>de</strong> echilibru Qx ( Q, x Q)<br />
∂ f1<br />
x1 = f1( x1Q, x2Q)<br />
+<br />
∂ x1<br />
∂ f1<br />
Q( x1 − x1Q)<br />
+<br />
∂ x2<br />
Q( x2 − x2Q)<br />
+ ...<br />
∂ f 2<br />
x2 = f2( x1Q, x2Q)<br />
+<br />
∂ x<br />
∂ f 2<br />
Q( x1 − x1Q)<br />
+<br />
∂ x<br />
Q( x2 − x2Q)<br />
+ ...<br />
1<br />
2<br />
91<br />
1 2 .<br />
Deoarece Q este punctul <strong>de</strong> echilibru f ( x , x ) = 0 si f ( x , x ) = 0.<br />
Notam<br />
x<br />
1<br />
= x<br />
1<br />
− x<br />
1Q , x<br />
2<br />
= x<br />
2<br />
− x<br />
2Q<br />
si<br />
a<br />
a<br />
11<br />
21<br />
∂ f<br />
=<br />
∂ x<br />
∂ f<br />
=<br />
∂ x<br />
1 1Q 2Q 2 1Q 2Q<br />
1<br />
1<br />
2<br />
1<br />
Q<br />
Q<br />
si<br />
si<br />
a<br />
a<br />
12<br />
22<br />
∂ f<br />
=<br />
∂ x<br />
1<br />
2<br />
∂ f<br />
=<br />
∂ x<br />
Daca se pot neglija termenii <strong>de</strong> ordin superior din <strong>de</strong>zvoltarea Taylor ecuatiile <strong>de</strong> stare se pot scrie:<br />
x = a x + a x<br />
x = a x + a x<br />
1 11 1 12 2<br />
2 21 1 22 2<br />
Astfel s-a aproximat ecuatia <strong>de</strong> stare neliniara cu o ecuatie liniara in jurul punctului <strong>de</strong> echilibru.<br />
2<br />
2<br />
Q<br />
Q
In jurul punctului Q, portretul <strong>de</strong> faza al circuitului neliniar autonom este similar cu<br />
portretul <strong>de</strong> faza asociat acestei ecuatii liniare. Deci comportarea calitativa a circuitului neliniar in<br />
jurul unui punct <strong>de</strong> echilibru se poate studia cu ajutorul circuitului liniar asociat ecuatiei <strong>de</strong> stare <strong>de</strong><br />
mai sus. Daca, <strong>de</strong> exemplu, pentru circuitul liniar originea este nod stabil sau instabil acelasi lucru<br />
se poate spune si <strong>de</strong>spre punctul <strong>de</strong> echilibru al circuitului neliniar.Aceeasi afirmatie se poate face<br />
si pentru focare. Atunci cand originea este centru in ecuatiile neliniare intervin termenii <strong>de</strong> ordin<br />
superior care nu mai pot fi neglijati si nu se mai poate face interpretarea comportarii calitative a<br />
circuitului neliniar in functie <strong>de</strong> mo<strong>de</strong>lul liniarizat.<br />
<strong>3.4.</strong>7. Oscilatii in circuitele neliniare<br />
Orice raspuns al unui circuit autonom neliniar este <strong>de</strong>terminat atat <strong>de</strong> starea initiala cat si <strong>de</strong><br />
excitatii (sursele in<strong>de</strong>pen<strong>de</strong>nte <strong>de</strong> curent continuu). Daca intr-un astfel <strong>de</strong> circuit exista cel putin un<br />
raspuns (tensiune sau curent) care este functie periodica <strong>de</strong> timp spunem ca circuitul oscileaza (este<br />
un oscilator).<br />
Un oscilator simplu utilizat in multe aplicatii este circuitul RLC serie in care elementele<br />
dinamice sunt liniare cu parametri pozitivi (L>0, C>0) si rezistorul este neliniar. Starea initiala<br />
C 2 L 2<br />
(uC(0) si I2(0)) <strong>de</strong>termina o energie totala W(<br />
0)<br />
= uC<br />
( 0)<br />
+ iL<br />
( 0) ≥ 0 acumulata in elementele<br />
2 2<br />
dinamice la t=0. Daca rezistorul neliniar este strict pasiv (u ⋅I>0) acesta va absorbi o putere<br />
pozitiva care se va transforma ireversibil in caldura. Energia disipata <strong>de</strong> rezistor provine din W(0)<br />
<strong>de</strong>ci daca rezistorul este strict pasiv lim Wt ( ) = 0 si uC() t ⎯⎯⎯ → 0, i t<br />
t→∞<br />
t L () ⎯⎯⎯ → 0 <strong>de</strong>ci nu<br />
→∞ t→∞<br />
poate exista nici un raspuns periodic. Rezulta ca rezistorul trebuie sa fie activ.<br />
Se poate arata ca daca un rezistor controlat in curent cu ecuatia constitutiva u= u() i<br />
satisface conditiile<br />
i) u
atunci circuitul RLC serie oscileaza. De exemplu caracteristica<br />
se poate obtine cu un circuit cu un amplificator operational (vezi paragraful 2), iar caracteristica<br />
3<br />
i<br />
u = − i are o alura similara.<br />
3<br />
Ecuatiile <strong>de</strong> stare ale circuitului RLC serie sunt<br />
i<br />
C L<br />
uC<br />
=− = f1( uC, iL)<br />
C<br />
0,<br />
∆=<br />
a11a22− a21a12 = . Deoarece T>0 rezulta ca λ 12 , sunt sau<br />
L<br />
LC<br />
numere reale pozitive sau numere complexe cu partea reala pozitiva. Daca λ 1 > λ 2 > 0 atunci<br />
originea este nod instabil si traiectoriile pornesc toate din origine<br />
2
in<strong>de</strong>partandu-se <strong>de</strong> zona din jurul acesteia. Daca λ1, λ 2 ∈C originea este focar instabil si<br />
traiectoriile se in<strong>de</strong>parteaza <strong>de</strong> zona din jurul originii cand t →∞. Pe masura ce punctul <strong>de</strong><br />
functionare( <strong>de</strong> coordinate uC(t) si iL(t)) se in<strong>de</strong>parteaza <strong>de</strong> origine valorile lui iL cresc si punctul <strong>de</strong><br />
functionare pe caracteristica ui 0 ) . Acestui rezistor ii corespund valori proprii cu<br />
partea reala negativa. Ca urmare intr-o zona care nu este in jurul originii traiectoriile nu tind spre<br />
infinit. Deoarece circuitul are un singur punct <strong>de</strong> echilibru in origine traiectoriile nu pot sa<br />
convearga <strong>de</strong>cat spre acest punct. Rezulta ca, <strong>de</strong>oarece traiectoriile nu se pot intersecta intre ele<br />
<strong>de</strong>cat in punctul <strong>de</strong> echilibru, trebuie sa existe o curba limita spre care tind aceste traiectorii.<br />
Aceasta curba este inchisa astfel incat parcurgand-o se obtin forme <strong>de</strong> unda periodice pentru<br />
iLsiuC . Acest rationament este doar o justificare din consi<strong>de</strong>rente fizice fara a fi o <strong>de</strong>monstratie a<br />
existentei si unicitatii ciclului limita. Aparitia oscilatiilor este ilustrata in continuare prin cateva<br />
exemple.<br />
Exemplul 1 - Oscilatorul liniar este circuitul RLC serie in care rezistorul are R = 0.<br />
iL<br />
uC<br />
2 1<br />
Daca ui () = 0 ecuatiile <strong>de</strong> stare sunt uC<br />
=− , iL<br />
= si cu notatiile ω 0 = solutia are<br />
C C<br />
LC<br />
forma iL() t = Asin( ω 0t − θ ), uC() t = ω 0LAcos( ω 0 t −θ)<br />
un<strong>de</strong> A si θ <strong>de</strong>pind <strong>de</strong> conditiile<br />
initiale. Eliminand timpul, din expresiile iL() t si uC() t rezulta ecuatia traiectoriei<br />
2<br />
L<br />
i<br />
2<br />
uC<br />
+ = A<br />
2<br />
L ω<br />
0 2<br />
2<br />
care este o elipsa.<br />
Portretul <strong>de</strong> faza contine o multime <strong>de</strong> elipse. In toate punctele unei elipse energia totala este<br />
C 2 L 2 1 2<br />
aceeasi Wt () = uC () t + iL () t = LA <strong>de</strong>ci oscilatia consta intr-un transfer al energiei<br />
2 2 2<br />
acumulate intre con<strong>de</strong>nsator si bobina. Aceste traiectorii <strong>de</strong> nivel energetic constant se numesc<br />
orbite.<br />
Observatii<br />
94
i) oscilatorul liniar este un mo<strong>de</strong>l i<strong>de</strong>alizat care nu tine seama nici <strong>de</strong> pier<strong>de</strong>rile din<br />
con<strong>de</strong>nsatoarele si bobinele reale si nici <strong>de</strong> rezistentele firelor <strong>de</strong> legatura<br />
ii) spre <strong>de</strong>osebire <strong>de</strong> ciclul limita, in orice vecinatate a unei orbite exista traiectorii inchise (orbitele<br />
foarte apropiate).<br />
Exemplul 2 Oscilatorul Van <strong>de</strong>r Pol este un ciruit RLC serie in care ui ()<br />
3 3<br />
i<br />
( / )<br />
stare sunt u ,<br />
C i L uC − iL −iL<br />
C =− L =<br />
. Facem schimbarea <strong>de</strong> variabila τ =<br />
L<br />
du<br />
dt<br />
du 1 di di 1<br />
= ⋅ , = ⋅ si notand ε =<br />
dτ LC dt dτ LC<br />
C<br />
L<br />
C C L L<br />
du i di ⎡ ⎛ 3<br />
i ⎞ ⎤<br />
C L L<br />
L<br />
=− , = ε ⎢uC<br />
−⎜ −iL<br />
dτ<br />
ε dτ<br />
⎜<br />
⎟ ⎥ .<br />
⎣⎢<br />
⎝ 3 ⎠ ⎦⎥<br />
95<br />
3<br />
i<br />
= −i.<br />
Ecuatiile <strong>de</strong><br />
3<br />
t<br />
LC<br />
si rezulta:<br />
ecuatiile <strong>de</strong> stare <strong>de</strong>vin:<br />
Aceste ecuatii nu au solutie analitica. Pentru anumite valori ale lui ε se pot <strong>de</strong>termina solutii<br />
analitice aproximative.<br />
Daca ε → 0 se pot face urmatoarele aproximatii:<br />
in ecuatia <strong>de</strong> <strong>ordinul</strong> 2 pentru iL<br />
2<br />
L<br />
2<br />
d i<br />
dτ<br />
2 diL<br />
=−iL −ε ( iL<br />
−1)<br />
dτ<br />
di ⎡ ⎛ 3<br />
i ⎞⎤<br />
L 2 L<br />
se poate consi<strong>de</strong>ra ca ε = ε ⎢uC<br />
−⎜−iL dτ<br />
⎜<br />
⎟<br />
⎟⎥<br />
≈ 0 si rezulta<br />
⎣⎢<br />
⎝ 3 ⎠⎦⎥<br />
2<br />
L<br />
2<br />
d i<br />
dτ<br />
=− i .<br />
Aceasta ecuatie are o solutie <strong>de</strong> forma iL = Acos(<br />
τ −θ<br />
) din care rezulta<br />
u<br />
1 A<br />
=− i dt =− −<br />
C ∫ sin( τ θ ) <strong>de</strong>ci traiectoria este o elipsa.<br />
ε<br />
C L<br />
In solutia ecuatiei simplificate A <strong>de</strong>pin<strong>de</strong> <strong>de</strong> conditiile initiale. Se poate arata (<strong>de</strong> exemplu prin<br />
integrare numerica) ca pentru ε < 01 , solutia ecuatiilor <strong>de</strong> stare nesimplificate are un ciclu limita<br />
eliptic corespunzator valorii A=2. Pentru valori 03 , < ε < 07 , ciclul limita este o elipsa <strong>de</strong>formata.<br />
L
Daca ε →∞ rezulta<br />
duC<br />
≈ 0<br />
dτ<br />
<strong>de</strong>ci uC = ct.<br />
In acest caz, daca<br />
uC −u( iL)<br />
≠ 0<br />
diL<br />
→ ∞.<br />
Deci daca uC−u( iL)<br />
≠ 0 traiectoriile vor fi paralele cu axa iL. Sensul<br />
dτ<br />
parcursului dinamic pe aceste traiectorii este:<br />
diL<br />
daca uC > u( iL),<br />
→ +∞ si iL creste<br />
dτ<br />
diL<br />
daca uC < u( iL),<br />
→ −∞<br />
dτ<br />
si iL sca<strong>de</strong>.<br />
Pentru u u i valoarea diL<br />
C ≈ ( L)<br />
→ ∞⋅0<br />
nu se poate <strong>de</strong>termina.<br />
dτ<br />
1 diL<br />
Se poate insa consi<strong>de</strong>ra ca uC − u( iL)<br />
= → 0<br />
ε dτ<br />
Deci pentru ε →∞ portretul <strong>de</strong> faza este:<br />
96<br />
<strong>de</strong>ci si caracteristica ui ( L ) este traiectorie.<br />
Acest portret <strong>de</strong> faza justifica regula <strong>de</strong> salt <strong>de</strong> la oscilatorul <strong>de</strong> relaxare realizat cu un<br />
circuit <strong>de</strong> <strong>ordinul</strong> I cu con<strong>de</strong>nsator liniar si rezistor cu caracteristica ui () (vezi paragraful 3.3.6.).<br />
Intr-a<strong>de</strong>var daca L → 0 in oscilatorul Van <strong>de</strong>r Pol rezulta ε →∞ iar circuitul <strong>de</strong>vine <strong>de</strong> <strong>ordinul</strong> I.<br />
Saltul lui I se face conform portretului <strong>de</strong> faza al circuitului <strong>de</strong> <strong>ordinul</strong> II. Timpul <strong>de</strong> salt este<br />
practic nul <strong>de</strong>oarece variatia in timp a lui iL este foarte rapida di ⎛<br />
⎜<br />
⎝ dt<br />
L →∞<br />
⎞<br />
⎟. Mo<strong>de</strong>lul <strong>de</strong> <strong>ordinul</strong> I<br />
⎠<br />
care are puncte <strong>de</strong> impas este un mo<strong>de</strong>l incorect. Mo<strong>de</strong>lul corect este cel <strong>de</strong> <strong>ordinul</strong> II in care<br />
intervine si inductivitatea L. In realitate aceasta inductivitate exista fiind un element parazit <strong>de</strong><br />
circuit asociat firelor <strong>de</strong> legatura dintre con<strong>de</strong>nsator si rezistorul neliniar.