Cap.2. Aplicatie: Analiza circuitelor electrice liniare (c.c. si …Tipuri de circuite / probleme...

1/45

IntroducereAnaliza n c.c.Analiza n c.a.

Cap.2. Aplicatie: Analiza circuitelor electriceliniare (c.c. si c.a.)

Prof.dr.ing. Gabriela Ciuprina

Universitatea "Politehnica" Bucuresti, Facultatea de Inginerie Electrica

Suport didactic pentru disciplina Metode numerice,Facultatea de Inginerie Electrica, 2017-2018

Gabriela Ciuprina Analiza circuitelor electrice liniare (c.c, c.a)

2/45


Cuprins

1 IntroducereModelareSimulare

2 Analiza circuitelor rezistive liniare n c.c.Formularea problemeiMetoda nodala clasica

3 Analiza circuitelor liniare n c.a.Formularea problemeiSimilitudinea cu c.c.Caracteristici de frecventa


Notes

Notes

3/45


ModelareSimulare

Circuitele electrice sunt modele ale realitatii

Circuitele electrice

modele ale realitatii;

contin elemente ideale, obtinute prin idealizareaelementelor reale;

reprezinta o multime de elemente ideale conectate ntreele pe la borne (terminale).


4/45


ModelareSimulare

Circuitele electrice sunt alcatuite din elemente ideale

Elementele ideale de circuit electric

sunt caracterizate de marimi electrice definite la borne(curenti, tensiuni sau potentiale);

se definesc functional, printr-o relatie caracteristica(constitutiva) ntre marimile definite la borne.

Modelarea nu este obiectul teoriei circuitelor, ea presupuneanaliza cmpului electromagnetic.

C = 0As


Notes

Notes

5/45


ModelareSimulare

Exemple de elemente ideale

Cele mai frecvent folosite:

liniare dipolare: R, L, C, conductorul si izolatorul perfect;

parametrice: K (comutatorul);

neliniare rezistive : SIT, SIC, DP;

liniare multipolare: SICU, SUCI, SUCU, SICI, AOP, M;

neliniare multipolare: AOPn.


6/45


ModelareSimulare

Exemple de elemente ideale

u

u

u

u

i

i

i

i

+


Notes

Notes

7/45


ModelareSimulare

Modelarea componentelor din circuitele reale


8/45


ModelareSimulare

Determinarea raspunsului sub actiunea unei excitatii

Simulare = simulare numerica (cu ajutorul calculatorului)

Simularea

determinarea marimilor de interes (tensiuni, curenti) dincircuit;

determinarea raspunsului sub actiunea unui semnal deexcitatie cunoscut.


Notes

Notes

9/45


ModelareSimulare

Determinarea raspunsului sub actiunea unei excitatii

O simulare facuta cu succes presupune

buna formulare a circuitului (solutia sa existe si sa fieunica); este echivalenta cu buna formulare a problemeimatematice asociate;

conceperea sau alegerea unui algoritm numeric robustpentru rezolvare.


10/45


ModelareSimulare

Algoritmul de rezolvare

Algoritmul potrivit pentru rezolvare depinde de

caracteristicile elementelor de circuit (liniare/neliniare,rezistive/reactive);

tipul marimilor din circuit (constante - c.c., sinusoidale -c.a., periodice, oarecare).


Notes

Notes

11/45


ModelareSimulare

Tipuri de circuite / probleme matematice

Tip de circuit

1 Circuite rezistiveliniare/neliniare n c.c.)

2 Circuite liniare n regimsinusoidal (c.a.);

3 Circuite liniare/neliniare nregim tranzitoriu;

4 Circuite liniare/neliniare nregim periodic;

5 Oscilatoare (frecvente derezonanta.)

Problema matematica1 Sisteme de ec. algebrice

liniare/neliniare, n IR;2 Sisteme de ec. algebrice

liniare, n complex.3 Sisteme ODE, lin./nelin. cu

conditii initiale.4 Superpozitie de c.a./ODE

cu conditii de periodicitate.5 Calcul de valori proprii

(analiza modala).


12/45


ModelareSimulare

Scopul acestui curs

ntelegerea:

modului n care se dezvolta instrumentele software pentruanaliza circuitelor electrice;

importantei bunei formulari a problemei (circuitului) cetrebuie rezolvata;

modului n care se genereaza automat sistemele derezolvat;

faptului ca fundamentul simularii numerice a circuitelorelectrice l constituie disciplina Metode numerice Algoritmi.


Notes

Notes

13/45


Formularea problemeiMetoda nodala clasica

Problema fundamentala

Contin: rezistoare (R), surse ideale de tensiune (SIT) si curent(SIC), surse comandate liniar (SUCU, SUCI, SICU, SUCI).

Problema fundamentala a analizei acestor circuite

Se dau:

topologia circuitului (schema/tabel de descriere(netlist)/matrice de incidenta sau apartenenta);

valorile parametrilor (rezistentele, valorile surselor).

Se cer:

curentii si tensiunile din fiecare latura;

puteri.


14/45



Conditii de buna formulare

TeoremeTopologice:

Pentru ca circuitul sa fie bine formulat este necesar saexiste un arbore normal;

Daca circuitul nu are surse comandate si toate rezistoarelesunt strict pozitive, atunci este necesar si suficient saexiste un arbore normal.

Algebrice:

Pentru ca circuitul sa fie bine formulat este necesar sisuficient ca matricea sistemului de ecuatii algebrice liniare,asamblat printr-o metoda sistematica sa fie nesingulara.

Q1: Ce este un arbore normal?


Notes

Notes

15/45



Metode de rezolvare sistematice

metoda ecuatiilor Kirchhoff :(

metoda potentialelor nodurilor :) (daca nu sunt sursecomandate matricea coeficientilor este simetrica sidiagonal dominanta)

metoda curentilor ciclici :| (daca nu sunt surse comandatematricea este simetrica, necesita definirea unui sistem debucle independente convenabil ales)

= metoda potentialelor nodurilor ("tehnica nodala")


16/45



Tratarea SRT

Laturi standard: Rk ikek

uk

(nik ) (nfk )

Formularea problemei

Se dau:

topologia: N, L, (nik ,nfk , k = 1, . . . ,L);

toate rezistentele Rk , k = 1, . . . ,L, presupuse nenule,

toate t.e.m. ek , k = 1, . . . ,L

Se cer:

uk k = 1, . . . ,L

ik k = 1, . . . ,L

puterea consumata si puterea generata n circuit.


Notes

Notes

17/45



Ecuatii

Rk ik

ek

uk

(nik ) (nfk )

Kirchhoff clasic:

k(n)

Aik = 0, n = 1, . . . ,N 1, (1)

k[b]

Auk = 0, b = 1, . . . ,L N + 1, (2)

uk = Rk ik ek , k = 1, . . . ,L, (3)

2L ecuatii cu 2L necunoscuteGabriela Ciuprina Analiza circuitelor electrice liniare (c.c, c.a)

18/45



Necunoscute

Rk ik

ek

uk

(nik ) (nfk )

Schimbare de variabila - necunoscutele sunt:vk , k = 1, . . . ,N, vN = 0 (prin conventie)Kirchhoff II:

k[b]

Auk = 0, b = 1, . . . ,L N + 1, (4)

uk = vnik vnfk , k = 1, . . . ,L. (5)


Notes

Notes

19/45



Notatii

u = [ u1 u2 . . . uL ]T IRL1

i = [ i1 i2 . . . iL ]T IRL1

v = [ v1 v2 . . . vN1 ]T IRN11

e = [ e1 e2 . . . eL ]T IRL1

R = diag([ R1 R2 . . . RL ]) IRLL

(6)

Kirchhoff I:Ai = 0, (7)

A = (aij)i=1,N1;j=1,L este matricea incidentelor laturi-noduri -matrice topologica, (N 1) L

aij =

0 daca nodul i nu apartine laturii j ;+1 daca nodul i este nod initial pentru latura j ;1 daca nodul i este nod final pentru latura j .


20/45



Ecuatii scrise compact

Kirchhoff I (KCL):Ai = 0, (8)

Kirchhoff II (KVL):u = AT v, (9)

Joubert (relatii constitutive):

u = Ri e. (10)

Daca R este inversabila (Rk 6= 0, k = 1,L)

i = R1(u + e). (11)

AR1AT v = AR1e. (12)

Gnv = jn. (13)


Notes

Notes

21/45



Sistem de ecuatii

Gnv = jn. (14)

Gn conductante nodale; jn injectii de curent n noduri.

Gn = AR1AT IR(N1)(N1) (15)

Gnii =

k(i)

1Rk

, Gnij =

k(i);k(j)

1Rk

pentru i 6= j .

jn = AR1e IR(N1)1 (16)

jnk =

m(k)

A em

Rm


22/45



Proprietatile matricei Gn

Gn: simetrica, diagonal dominanta si pozitiv definita dacarezistentele sunt pozitiveA IRnn este pozitiv definita daca ea este simetrica si daca xT Ax > 0 pentru orice vector real, nenul x IRn1.

R1 = diag([ 1/R1 1/R2 . . . 1/RL ]). (17)

Simetria:

GTn =(

AR1AT)T

=(

AT)T (

R1)T

(A)T = AR1AT = Gn

Pozitiv definire: Fie x vector coloana arbitrar, nenul.

xT Gnx = xT AR1AT x = yT R1y =

L

k=1

y2kRk

> 0,

unde y = AT x are componentele yk , k = 1, . . . ,L.Gabriela Ciuprina Analiza circuitelor electrice liniare (c.c, c.a)

Notes

Notes

23/45



Etapele algoritmului

etapa de preprocesare n care se descrie problema si seasambleaza sistemul de ecuatii de rezolvat;

etapa de rezolvare n care se apeleaza o procedurapropriu-zisa de rezolvare a sistemului de ecuatii rezultat("solver");

etapa de postprocesare n care se calculeaza alte marimide interes.


24/45



Structuri de date

Rk ik

ek

uk

(nik ) (nfk )

; declaratii date - varianta Antreg N ; numar de nodurintreg L ; numar de laturitablou ntreg ni[L] ; noduri initiale ale laturilortablou ntreg nf[L] ; noduri finale ale laturilortablou real R[L] ; rezistentetablou real e[L] ; tensiuni electromotoare

n vederea obtinerii unui algoritm simplu, vom presupune ca:sensul de referinta al curentului unei laturi este identic cucel al t.e.m de pe latura;toate laturile sunt orientate cf. regulii de la receptoare.


Notes

Notes

25/45



Structuri de date

Rk ik

ek

uk

(nik ) (nfk )

Se recomanda agregarea datelor:

; declaratii date - varianta Bnregistrare circuit

ntreg N ; numar de nodurintreg L ; numar de laturitablou ntreg ni[L] ; noduri initiale ale laturilortablou ntreg nf[L] ; noduri finale ale laturilortablou real R[L] ; rezistentetablou real e[L] ; tensiuni electromotoare


26/45



Matrice rare

Gn si jn sunt foarte rare.

Exemplu:daca pp. 4 laturi care concura la un nod, atunci densitateamatriceid = 5n/n2 = 5/n, (pentru n 1000 d = 0.5 %).

Pentru simplitate:

; declaratii variabile utiletablou real Gn[N, N] ; stocata rartablou real jn[N] ; stocat rartablou real v [N] ; vectorul potentialelor


Notes

Notes

27/45



Citire date

Rk ik

ek

uk

(nik ) (nfk )

functie citire_date_B (); declaratii...citeste circuit.N, circuit.Lpentru k = 1,circuit.L

citeste circuit.nik , circuit.nfkciteste circuit.Rk , circuit.ek

ntoarce circuit


28/45



Asamblarea sistemului de ecuatii

Orientata pe laturi:

Rk ik

ek

uk

(nik ) (nfk )

nik nfk

nik +1/Rk 1/Rk

nfk 1/Rk +1/Rk

nik ek/Rk

nfk +ek/Rk

Contributia unei laturi k la matricea conductantelor nodale (stnga) si la vectorul injectiilor de curent (dreapta).


Notes

Notes

29/45



Preprocesare

Rk ik

ek

uk

(nik ) (nfk )

procedura nodalRE_v1 (circuit, Gn, t); asambleaza sistemul de ecuatii pentru un circuit; cu laturi de tip R,E folosind tehnica nodala; parametri de intrare:; circuit - structura de date ce descrie circuitul; parametri de iesire:; Gn - matricea conductantelor nodale si; jn - vectorul injectiilor de curent; declaratii....L = circuit.L ; pentru simplificarea scrierii algoritmuluiN = circuit.Nni = circuit.ninf = circuit.nfR = circuit.Re = circuit.e


30/45



Preprocesare

Rk ik

ek

uk

(nik ) (nfk )

procedura nodalRE_v1 (circuit, Gn, jn)....Gn = 0jn = 0; asambleaza sistempentru k = 1, L ; parcurge laturi

i = nik ; nodul initial al laturii kj = nfk ; nodul final al laturii kGnii = Gnii + 1/RkGnjj = Gnjj + 1/RkGnij = Gnij 1/RkGnji = Gnji 1/Rkjni = jni ek/Rkjnj = jnj + ek/Rk

retur


Notes

Notes

31/45



Preprocesare

Observatii:

am folosit pseudocod simplificat pentru a scrie anulareacomponentelorAtentie! varianta

pentru i = 1,Npentru j = 1,N

Gnij = 0

scrisa pentru "instructiunea" Gn = 0 va umple completmatricea Gn.

pentru a evita repetarea unor calcule, se pot memoravalorile 1/Rk si ek/Rk .


32/45



Preprocesare - varianta a II-a

Rk ik

ek

uk

(nik ) (nfk )

procedura nodalRE_v2 (circuit, Gn, jn)....; anuleaza componentele:A = 0 ; matricei incidente laturi noduriG = 0 ; matricei diagonale R1

; asambleaza sistempentru k = 1, L ; parcurge laturi

i = nik ; nodul initial al laturii kj = nfk ; nodul final al laturii kAik = 1Ajk = +1Gkk = 1/Rk

Gn = A G AT ; apel proceduri speciale pentru matrice rarejn = A G eretur


Notes

Notes

33/45



Rezolvare

Sistemul asamblat are dimensiunea N N, nodul dereferinta nefiind tratat special.

Sistemul de rezolvat trebuie sa aiba dimensiunea N 1.Dupa rezolvare trebuie adaugata o componenta n plusvectorului potentialelor: vN = 0.

Exemplu:

Gauss (N 1,G,t ,v)vN = 0

Q2: Cum implementati aceasta idee n Matlab/Octave ?


34/45



Rezolvare

Metode posibile de rezolvare:

directe (Gauss, factorizare) - nu introduc erori detrunchiere, dar matricele se umple n cursul algoritmului;

iterative (Jacobi, Gauss-Seidel, SOR) - matricele sipastreaza gradul de raritate, dar apar erori de trunchiere sieventuale probleme de convergenta;

semiiterative (gradienti conjugati, GMRES, etc) -avantajoase daca matricea sistemului este simetrica sipozitiv definita (daca nu exista surse comandate).


Notes

Notes

35/45



Postprocesare

Rk ik

ek

uk

(nik ) (nfk )

procedura postprocesare_circuitRE (circuit, v )...Pc = 0 ; puterea consumataPg = 0 ; puterea generatapentru k = 1, L ; parcurge laturi

u = vnik vnfk

; tensiunea laturiic = (u + ek )/Rk ; curentul prin laturascrie "Latura" k "are tensiunea" u "si curentul" cPc = Pc + Rk c

2 ; adauga contributia laturii la PcPg = Pg + ek c ; adauga contributia laturii la Pg

scrie Pc, Pgretur

Q3: Cum implementati postprocesarea n Matlab/Octavefolosind operatii cu matrice?


36/45


Formularea problemeiSimilitudinea cu c.c.Caracteristici de frecventa


Contin:

rezistoare liniare (R);

bobine liniare (L);

bobine liniare cuplate (M);

condensatoare liniare (C);

surse ideale de tensiune (SIT);

surse ideale de curent (SIC);

surse comandate liniar (SUCU, SUCI, SICU, SUCI).

SIT sau SIC au variatii de forma:

y(t) = Y

2 sin(t + ). (18)

unde are aceeasi valoare pentru toate marimile.


Notes

Notes

37/45




Problema fundamentala a analizei circuitelor de c.a.

Se dau:

topologia circuitului (schema/tabel de descriere(netlist)/matrice de incidenta sau apartenenta);

valorile parametrilor (rezistentele, bobinele, cuplajele,condensatoarele, valorile surselor: frecventa, valorileefective, fazele initiale).

Se cer:

curentii si tensiunile din fiecare latura (valori efective, fazeinitiale);

puteri (active, reactive, aparente, defazaje).


38/45



Similitudinea cu c.c.

Metoda de analiza se bazeaza pe reprezentarea ncomplex.

y(t) = Y

2 sin(t + ) Y = Y ej. (19)

Ideea: ecuatiile similare:

Circuitul de c.c. Circuitul de c.a.

TK1(A)

k(n)ik = 0

(A)k(n)

Ik = 0

TK2(A)

k[b]uk = 0

(A)

k[b]Uk = 0

SRT uk = Rk ik ek Uk = Z k Ik EkSRC ik = Gk uk + jk Ik = Y k Uk + JkSUCI ek = rkm im Ek = zkm ImSICU jk = gkmum Jk = ykmUmSUCU ek = kmum Ek = kmUmSICI jk = km im Jk = km Im


Notes

Notes

39/45



Reprezentarea n complex a elementelor ideale

Rezistor (R) Bobina (L) Condensator (C)Impedanta complexa Z R jL 1/(jC)Admitanta complexa: Y 1/R 1/(jL) jC

Defazajul: 0 /2 /2Impedanta: Z R L 1/(C)Admitanta: Y 1/R 1/(L) C

Rezistenta de c.a.: R R 0 0Reactanta: X 0 L 1/(C)

Conductanta de c.a.: G 1/R 0 0Susceptanta: B 0 1/(L) C


40/45



Algoritm

Similar cu cel din c.c.:

n loc de rezistente se lucreaza cu impedante complexe;

parametrii surselor sunt tot valori constante, dar complexe,obtinute din reprezentarea n complex a variatiilor care sedau.

Diferente fata de algoritmul din c.c.:

n etapa de preprocesare: citirea datelor de descriere sireprezentarea lor n complex;

n etapa de asamblare, apar n plus bobinele cuplate, carecontribuie la sistem cu urmatoarele stampile:


Notes

Notes

41/45



Algoritm

Cuplaje

nij nfj nik nfk

Am

[

+1 1 0 00 0 +1 1

]

j k

Bm

nij

nfj

nik

nfk

+1 01 0

0 +10 1

j k

Zm

j

k

[

jLjj jLjkjLkj jLkk

]

em Nu contribuie

jn

[

ijik

]


42/45



Caracteristici de frecventa

n multe aplicatii practice intereseaza reprezentareacaracteristicilor de frecventa: comportarea semnalelor de iesirepentru un interval al frecventelor semnalelor.Variante de implementare:

1 Se lucreaza simbolic, cu parametrul si se obtin expresiisimbolice ale marimilor de iesire care apoi se evalueazanumeric;

2 Se lucreaza numeric, pentru frecvente din intervalul deinteres se rezolva mai multe probleme de c.a.


Notes

Notes

43/45



Lectura obligatorie pentru aceasta saptamna

Cap.5 din[1] Gabriela Ciuprina, Mihai Rebican, Daniel Ioan - Metode numerice in ingineria electrica - Indrumar de

laborator pentru studentii facultatii de Inginerie electrica, Editura Printech, 2013, disponibil la

http://mn.lmn.pub.ro/indrumar/IndrumarMN_Printech2013.pdf

Nota: daca folositi Matlab, nu aveti voie la acest curs sa folositi operatii cu vectori si matrice dect pentru validarearezultatelor, nu pentru implementarea procedurilor.De exemplu, ntr-o prima varianta puteti folosi mldivide (backslash) pentru rezolvarea sistemului de ecuatii asamblat,pentru a verifica programul, dar n final nlocuiti-o cu una din procedurile de rezolvare pe care le-ati implementat voi.


44/45



Simulatoare de circuit

Free and Open Source

NgSpice (are si varianta online), GnuCap, CircuitLogix,LTSpice, MultiSim, TopSpice, MacSpice, Xyce (opensource, SPICE-compatible, high-performance analogcircuit simulator)

Licensed/Paid Circuit simulation software

Spectre (Cadence), PSpice, MultiSim, SiMetrix, TINA

Vedeti sihttp://www.circuitstoday.com/circuit-design-and-simulation-softwares

https://en.wikipedia.org/wiki/Electronic_circuit_simulation


Notes

Notes

http://mn.lmn.pub.ro/indrumar/IndrumarMN_Printech2013.pdfhttp://www.circuitstoday.com/circuit-design-and-simulation-softwareshttps://en.wikipedia.org/wiki/Electronic_circuit_simulation

45/45



Tema pentru bonus

1 Scrieti un program pentru analiza circuitelor de curent alternativ pentru circuite care contin rezistoare,bobine necuplate, condensatoare si surse independente de tensiune.

2 Alegeti pentru testarea codului un exemplu simplu (de exemplu, dar nu obligatoriu, un filtru pasiv adica faraA.O., din lista http://sim.okawa-denshi.jp/en/Fkeisan.htm. Structura de date pentru circuitul de test ales va fiinstantiata ntr-o functie (nu se vor cere date de la tastatura).

3 Verificati solutia comparnd-o cu o solutie de referinta care poate fi: analitica sau obtinuta cu un instrumentde tipul calculator online http://sim.okawa-denshi.jp/en/CRlowkeisan.htm

4 Verificati solutia comparnd-o cu un simulator de circuit de tipul spice - va recomandam: ngspice varianta online disponibila la http://www.ngspice.com/ sau

LTSpice http://www.linear.com/designtools/software

Scrieti un raport care sa rezolve punctele de mai sus. Este obligatoriu ca raportul sa aiba: o pagina de titlu, uncuprins generat automat, o lista de referinte. Dati o structura coerenta raportului.Fisierele care rezolva tema se vor organiza ntr-un folder numit NumePrenume_grupa. n acest folder vor existaurmatoarele subfoldere care vor contine fisiere relevante: raport, surse, spice. Folderul se arhiveaza (zip) si sencarca pe moodle.

Termenul de predare a acestei temei va fi anuntat pe moodle.


Notes

Notes

http://sim.okawa-denshi.jp/en/Fkeisan.htmhttp://sim.okawa-denshi.jp/en/CRlowkeisan.htmhttp://www.ngspice.com/http://www.linear.com/designtools/software

IntroducereModelareSimulare

Analiza circuitelor rezistive liniare n c.c.Formularea problemeiMetoda nodala clasica

Analiza circuitelor liniare n c.a.Formularea problemeiSimilitudinea cu c.c.Caracteristici de frecventa

Cap.2. Aplicatie: Analiza circuitelor electrice liniare (c.c. si …Tipuri de circuite / probleme...

Documents

Transcript of Cap.2. Aplicatie: Analiza circuitelor electrice liniare (c.c. si …Tipuri de circuite / probleme...