29.06.2013 Views

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

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!