UNIVERSITATEA POLITEHNICA BUCURETI FACULTATEA DE ENERGETIC CATEDRA SISTEME ELECTROENERGETICE

METODA PERTURBAIILOR SINGULARE N SIMULAREA UNUI GRUP ENERGETIC

COORDONATOR TIINIFIC: Conf. Dr. Ing. Clin Soare ABSOLVENT: Ruu Matei
BUCURETI

2011

CUPRINS

1. INTRODUCERE 2. METODA PERTURBATIEI SINGULARE UTILIZARE. 2.1 Introducere 2.2 Prezentarea metodei perturbaiei singulare pentru sisteme cu doua scale de timp 2.3 Separarea iterativ a scalelor de timp 3. SEPARABILITATEA N SCARA DE TIMP. EXEMPLIFICRI 3.1 Introducere 3.2 Regulatorul automat de tensiune 3.2.1 Necesitatea i avantajele introducerii regulatorului automat de tensiune 3.2.2 Separabilitatea mrimilor de stare utilizate n descrierea matematic a regulatorului automat de tensiune 3.3 Transformatorul electric 3.3.1 Elementele constructive ale transformatorului monofazat 3.3.2 Separarea scalelor de timp aplicat transformatorului monofazat 4. MODELAREA DINAMICII UNUI GRUP ENERGETIC CUPLAT LA SISTEMUL ENERGETIC 4.1 Noiuni generale despre maina sincron 4.2 Ecuaiile regimului tranzitoriu al generatorului sincron 4.2.1 Ecuaiile generatorului sincron n mrimi instantanee reale 4.2.2 Transformarea Park. Ecuaiile funcionale ale generatorului sincron n sistemul de referin rotativ (dk, qk, o)

4.2.3 Cazul particular al raportrii la un referenial legat de rotor ( 4.3 Ecuaiile funcionale ale generatorului sincron n mrimi relative

4.4 Modelul matematic al generatorului sincron prin evidenierea tensiunilor electromotoare

5. ANALIZA REGIMULI DINAMIC AL UNUI GRUP ENERGETIC PRIN METODA PERTURBAIILOR SINGULARE 5.1 Introducere 5.2 Schema grupului electric cuplat la sistemul energetic. Parametrii tehnici ai grupului, precum i ai liniei de cuplare 5.3 Reprezentarea sistemic a grupului energetic analizat 5.4 Aplicarea procedurii iterative de separare 5.5 Elemente necesare simulrii dinamice a grupului energetic prin metoda perturbaiei singulare

6. PROGRAM NUMERIC DE SIMULARE DINAMICA A UNUI GRUP ENERGETIC PRIN METODA PERTURBAIILOR SINGULARE

7. VALIDAREA METODEI PRIN SIMULARE NUMERIC

8. CONCLUZII

BIBLIOGRAFIE

1.

INTRODUCERE

Creterea exploziv a consumului de energie electric si termic n absolut toate sectoarele economiei, precum i concentrarea acestui consum, fie n unitai mari consumatoare de energie, fie in zone cu densitai mari din punct de vedere al consumului energetic, au impus conditii din ce in ce mai pretentioase in realizarea si exploatarea ansamblului de instalaii i mijloace tehnice din sistemele electroenergetice actuale. Dezvoltarea interconexiunilor i modificarea structurilor marilor sisteme energetice interconectate, creterea puterilor unitare ale agregatelor generatoare din centrale electrice, ca i a valorilor nominale ale tensiunilor, liniilor de transport i distribuie a energiei electrice, au determinat condiii tot mai severe de funcionare a instalaiilor de reglare i de automatizare n vederea creterii stabilitii sistemului i pentru prevenirea avariilor. Totodata au fost concepute, realizate i impementate n instalaii noi mijloace de automatizare, conducere i prelucrare automat a datelor, mult perfecionate n raport cu echipamentele clasice din exploatare. n activitatea de proiectare a noilor echipamente s-au impus introducerea de noi metode i procedee de analiz. Un loc important n cadrul acestora l ocupa studiul regimurilor tranzitorii determinate de apariia unor perturbaii n funcionarea normala a sistemului, cu urmari nefaste asupra nivelelor de frecvena i tensiune i asupra siguranei n funcionare. Descrierea matematic a dinamicii fenomenelor tranzitorii presupune determinarea unor modele matematice ct mai simple. Obiectivul proiectului este acela de a preyenta analiza regimului dinamic al unui grup energetic cuplat la sistemul energetic prin metoda perturbaiei singular pentru separarea pe scale de timp. Sistemul testat este constituit dintr-un generator sincron ce debiteaz pe bare de putere infinit, prevzut cu un regulator automat de tensiune cu rolul de a mentine tensiunea constant la generator i implicit la consumatori, independent sau slab dependent la perturbaii. Pentru cazul studiat, reglarea tensiunii pe barele generatorului este realizat prin vibraia tensiunii electromotoare, prin variaia excitaiei generatorului sincron reglat.
5

n cadrul metodei de analiz este prezentat un algoritm iterativ de separare a modurilor lente i rapide astfel nct sunt nlturate inconsecvenele aproximrii clasice de stare cvasistaionar, iar precizia modelelor de ordin inferior este mbuntait sistematic la fiecare pas iterativ. Pentru separarea scalelor de timp se introduce un parametru specific metodei, numit parametru de scalare, care are o valoare numeric mic, prin intermediul cruia se separ subsistemele cu dinamici mult diferite. Acest parametru rezult ca o perturbaie inerent a modelului i anume parametrul este nul conducnd la erori. Dar realizarea calculelor pe baza acestei ipoteze ne conduce la obinerea unor modele separate de ordin inferior. Cu ajutorul metodei iterative se corectez i se mbuntete precizia acestor modele, iar n plus noile variabile de stare introduse au aceeai semnificaie fizic ca variabilele originale. Metoda poate fi utilizat att pentru sistemele bine amortizate, ct i pentru sistemele cu amortizare uoar a oscilaiilor de frecven nalt, dar interpretrile procedurii de reducere sunt diferite n cadrul celor dou tipuri de sisteme. ntr-unul dintre capitole, metoda va fi expus prin exemplele a dou sisteme elementare, componente ale unui grup energetic: regulatorul automat de tensiune i transformatorul electric. Se va constata c nu ntotdeauna variabilele iniial alese se preteaz la o separare direct pe scale de timp. Din acest motiv, una dintre problemele majore o constituie identificarea i realegerea variabilelor lente i rapide. Totalitatea strilor unui model sistemic pe mai multe scale de timp se spune c este separabil pe stri daca prtile rapide ale unora dintre aceste stri sunt mai mici n comparaie cu prile lor lente i cu prile rapide ale altor stri. Analiza sistemului dinamic al unui generator sincron ce debiteaz pe bare de putere infinit, produs de scurtcircuitul trifazat care conduce la pierderea unei linii, pornete de la modele matematice ale subsistemelor componente i anume modelul matematic scris n mrimi relative, avnd drept referin axele rotorului, care evideniaz tensiunile electromotoare de inducie, mpreun cu ecuaiile de stare ale regulatorului dinamic de tensiune, scrise pe baza schemei bloc a acestuia, n care se obine un model de ordinul 7. Asupra modelului liniarizat de ordinul 7 se va aplica procedura de separare a variabilelor i se vor introduce scalele de timp. Rezultatele simulrii post-defectului pe baza subsistemelor lente i rapide indic ntr-un mod evident faptul c modelul cvasistaionar conduce la erori substaniale n variaia mrimilor de stare liniarizate faa de modelul exact, n timp ce utiliznd metoda iterativ cu un singur nivel de corecie aplicat subsistemelor rapid i lent, eroarea este insesizabil.
6

Plecnd de la rezultatele obinute, metoda analizei pe mai multe scale de timp poate fi extins pentru studiul regimului dinamic al unui sistem energetic complex, dar aceasta implic o procedur algebric mult mai sistematic pentru separarea strilor, precum i luarea n considerare a interaciunilor ce intervin ntre diferitele subsisteme componente ale ntregului sistem. De asemenea, sunt n curs de dezvoltare mbunatairi n ceea ce privete calculele prezentate n algoritmii folosii n procedura de identificare i validare.

2. METODA PERTURBAIEI SINGULARE. UTILIZARE

2.1 INTRODUCERE

Modelarea sistemelor de scar implic studierea interaciunilor care se produc ntre fenomene de dinamici foarte diferite. Astfel, n descrierea matematic a mainii sincrone, ecuaiile difereniale care corespund prii mecanice sunt mult mai lente dect cele care corespund prii electrice. Astfel de modele sunt de ordin ridicat, greu de rezolvat i se urmrete reducerea ordinului i separarea scrilor de timp prin metoda perturbaiei singulare, utiliznd agregarea, analiza modal i tehnici similare. Modelarea acestor sisteme se bazeaz pe ipoteza c, variabilele lente rmn constante n timpul tranziiilor i atunci cnd ncep s varieze, tranziiile rapide au ajuns la un nivel cvasistaionar. Pornind de la aceast ipotez, vom mpari variabilele de stare n mrimi lente, notate cu x i marimi rapide, notate cu z, modelul complet de stare fiind:

Se are n vedere numai variaia mrimilor z, considernd mrimile x constante. La momente mari de timp se studiaz numai mrimile x, ecuaiile difereniale pentru z reducndu-se la ecuaii algebrice sau transcendente, ce se scriu n forma z = 0. Modelul cvasistaionar se scrie:

Conform relatiei (2.4) zs este o mrime variabil n timp ceea ce contrazice aproximarea staionar pentru care zs este constant, aa cum implic = 0. Mai mult,

conditia iniial pentru z a trebuit s fie ignorat n relaia (2.4) deoarece nu mai putea fi ndeplinit. Pentru a nltura inconsistenele modelului clasic cvasistaionar, se va utiliza o metod iterativ de separare a variabilelor n mrimi lente i rapide, metod care mbuntete sistematic i precizia modelelor de ordin inferior. Aceast metod modific ipoteza de stare cvasistaionar ntr-o scar multitimp, potrivit pentru sistemele perturbate singular, aplicnd modificarea ipotezei cvasistaionare la fiecare pas al iteraiei cu modelul obinut la pasul precedent. Astfel, precizia modelelor pentru prile lente si rapide este mbuntit la fiecare pas. Iteraiile sunt nlanuite, dar mai uor de interpretat dect modelele standard de cretere asimptotic. Utilizarea succesiv a modificrii ipotezei cvasistaionare se poate face fr reinere n teoria perturbaiilor singulare.

2.2 PREZENTAREA METODEI PERTURBAIEI SINGULARE PENTRU SISTEME CU DOU SCRI DE TIMP

Pentru ilustrarea metodei vom considera un sistem cu dou scri de timp: scara de timp t corespunztoare fenomenelor lente i scara de timp ce descrie fenomenele rapide. care are o valoare variabila de

Raportul celor dou scri de timp este n general un parametru pozitiv

mic i care va fi parametrul principal al analizei asimptotice. Utiliznd timp se poate scrie n funcie de variabila t ca fiind:

a crui valoare iniial t din scara de timp t. Parametrul

= 0 este aleas astfel nct s corespund unei valori particulare

va fi cu att mai mic cu ct diferena dintre cele dou scri de timp este mai mic, va fi mai mare pentru un interval

este mai mare. Pe de alta parte, cu ct (t t) dat. La limit, cnd interval infinit in . Atunci cnd

, chiar un interval mai mic in t este extins la un este suficient de mare, fenomenele rapide ating starea

staionar. Aceasta nu contrazice ipoteza conform creia ( t-t ) este suficient de mic pentru a considera variabilele lente ca fiind constante. n acest fel cu ipoteza cvasistaionar, dar fr contradiciile acesteia. Daca se cunoate c dinamica mrimilor z este a lui x, atunci z este cu aproximativ redimensionat astfel: G= i ca atare f si g vor fi de acelai ordin de mrime. Modelul (2.1), (2.2) devine:
9

este echivalent

ori mai rapid n timp dect cea

mai mare n timp dect x si G poate fi

(2.6)

Raionamentul calitativ de mai sus este bazat pe nite estimri empirice ale derivatelor i . Cnd aceast informaie nu este disponibil, se verific parametrii

fizici, cum ar fi constantele de timp ce nmagazineaz energie, pentru a determina care dintre mrimi sunt lente i rapide. Nu toate variabilele de stare pot fi separate n acest fel, dar acolo unde separarea este posibil, un model (2.7), (2.8) poate fi obinut exprimnd constantele mici de timp i inversul coeficienilor de amplificare mari ca multiplu de un singur parametru mic . La limit, cnd , modelul (2.7), (2.8) scris n scara de timp t definete starea

cvasistaionar xs(t), zs(t) ca fiind:

(2.10)

Acest model este identic cu modelul (2.3), (2.4) dar originea i sensul su sunt diferite. Principala diferen const c partea stng, ci prin a relatiei (2.10) este ndeplinit acum nu prin , ci prin

10

Pentru a obine componentele lui x si z, se rescrie (2.7), (2.8) n scara de timp rapid

i se verific din nou la limit

. Rezult c

= 0, ceea ce nseamn ca x este singurele variaii rapide sunt =0

constant n scara de timp rapid. C atare, atunci cnd

cele ale lui z de la starea cvasistaionar zs. Notm aceste variaii cu zf =z-zs i facem n (2.11) i (2.12) obinnd:

Valoarea momentan fixat t a fost aleas ca fiind t0 i prin urmare constantele modelului sunt t0, x0, zs0 = zs(t0) care descriu fenomenele rapide din jurul lui t0. Utiliznd relaiile (2.9), (2.10), ca model lent i (2.13) ca model rapid, se aproximeaz mrimile x si z prin:

(2.15)

11

n care

este exprimat in scara de timp t. Perturbaiile singulare se preteaya mmult

mai bine la aceste reyultate dect la orice alt model simplificator. Parametrul cu care se opereaz, nentlnit la celelelalte modele este parametrul de scalare . Att timp ct modelele complete (2.7), (2.8), (2.11), (2.12) sunt exacte, separarea modelelor de ordin mai mic (2.9), (2.10), (2.13) conduce la erori, deoarece presupune n loc de inerent a modelului. Aproximarea (2.14), (2.15) poate fi mbuntait construind extinderea asimptotic n Este important ca fiecare termen al cresterii s fie calculat la , meninnd . Acest parametru este numit singular i rezult ca o perturbaie

avantajul de a avea modele separate de ordin inferior. Pentru aceasta, se va utiliza o metod iterativ de separare a scrilor de timp pentru un sistem liniar, metod care se va ilustra pe cazul concret al unui grup energetic, la care se poate aplica o procedur de corecie pentru sistemele neliniare.

2.3 SEPARAREA ITERATIV A SCALELOR DE TIMP

Ca un caz special al sistemului (2.7), (2.8), se consider un sistem de forma:

i se presupune c matricea D este inversabil. Starea cvasistaionar este definit de (2.17): , de unde rezult conform relaiei

adic:
12

Valorile reale x, z vor diferi de xs, zs n principal prin partea lor rapid. Pentru a gsi componenta rapid a mrimii z introducem variabila ca o diferen ntre z i zs:

Deci: de unde rezult c:

nlocuind pe z i z din relaiile (2.19) , (2.20) n relaiile (2.16) , (2.17) se obine:

Aceasta reprezint modelul de tipul (2.7), (2.8) n care

joac rolul lui z.

Principala diferen ntre cele dou modele o constituie prezena slab a lui x n ecuaiile lui deoarece matricea notat cu C1 este 0( ), adic o funcie care dup mprirea cu

rmne mrginit n vecintatea valorii Starea cvasistaionar a lui obtinut din relaia (2.22):

este 0( ), ceea ce arat c

este predominant rapid.

13

n mod asemntor se introduce variabila cauzat de presupunerea cvasistaionar se scrie:

reprezentnd eroarea

innd cont de relaia (2.23) aceast eroare

i n continuare se va nlocui variabila (2.21), (2.22).

n funcie de variabilele

i x n modelul

Repetnd aceast procedur dupa k pai, se obine:

care va conduce la sistemul:

ale crui matrice sunt definite de:

Se observ c n aceste relaii,

a fost redus, de aceast data la

, intrarea x

la fiecare model rapid succesiv este mai slab. Pentru a reconstrui z din i x, se observ din (2.25) c:

Diagrama bloc de reprezentare a iteraiilor este prezentat n figura 2.1. Viteza de integrare n bucl rapid este mare datorit factorului de amplificare mare, Intrarea ,

de la mrimea x la fiecare model rapid succesiv este mai slab. La limit, cnd k
14

modelul devine cascad rapid lent n care

conine toate modurile lente iar

toate modurile rapide. Utiliznd relaiile (2.28), (2.29), (2.30) se poate determina matricile subsistemelor lente i rapide astfel:

se obin din termenii matricilor subsistemului original fr transformri modale

ru condiionate. Un alt avantaj al acestei metode este acela c sensul fizic al variabilelor de stare iniiale se pstreaz. Astfel, din (2.31) observm c definirea lui x rmne aceeai, n timp ce noua variabil de stare are acelai sens ca i variabila z.

x C

Ck

Dk

15

Fig 2.1 Diagramele bloc de reprezentare ale iteraiilor

Se observ c dup k iteraii, subsistemul lent (2.26) al modelului are nc intrarea rapid complet. Se urmrete s se decupleze total cele dou subsisteme i pentru din relaia (2.27) n funcie de

aceasta se procedeaz n felul urmtor: se exprim i x:

care duce la schimbarea de variabil:

ca parte lenta a lui x. nlocuind pe x din relaia (2.33) n relaia (2.32) subistemul lent va deveni:

Prin urmare, intrarea rapid

a fost redus, deoarece i ,

este

. , n relaia

nlocuind mrimea x n funcie de (2.27), se obine:

16

Se definete unde rezult:

ca parte lent a lui

i se continua procedura n mod iterativ, de

unde:

Micorarea intrrii rapide a fost realizat, pentru c dup fiecare iteraie, reduce cu un grad al lui i tinde ctre 0 cnd j .

se

Cu acestea, subsistemele lente i rapide ale subsistemului rezultant sunt doar slab cuplate ntre ele deoarece este , iar este :

Se observ c

reprezint aproximaii ale matricilor i i pe baza relaiei (2.36) se obine:

Se reconstituie x din

Astfel variabila lenta

este partea dominant a lui x, partea rapid fiind are acelai sens fizic ca i x.

importan deosebit o constituie faptul ca

17

Aceste consideraii rmn valabile i dac se utilizeaz D. Aceasta nseamn c

, n loc de C i

, care este foarte important n analiza asimptotic de

convergen i de validare, nu trebuie s fie prezent explicit n iteraii. Obiectivul de a reduce un sistem la prile lente i rapide cuplate ntre ele este atins. n sistemul transformat (2.40), (2.41) termenii de cuplare i au valori foarte

mici i prin urmare se pot neglija, astfel nct n locul sistemului original complet (2.40), (2.41) se vor utiliza subsistemele decuplate de ordin mai mic:

cu valorile iniiale

obinute pe baza valorilor x0, z0 conform relaiei (2.31) i . Simularea lui se poate realiza n

obinute din relaia (2.42) pe baza valorilor x0 i scara de timp rapid . Eroarea i este

. Utiliznd mrimile

n locul variabilelor

se obin aproximrile corespunztoare pentru x(t), z(t) prin evaluarea lui din relaia (2.31). Eroarea este este ,

din relaia (2.42) i a lui unde i=min(j,k), n timp ce eroarea

unde l=min(j,k). Obinndu-

se acest model, se studiaz n continuare o nou simplificare i anume pstrarea a numai uneia din cele dou relaii prezente n sistemul (2.43). n general este necesar s se calculeze pentru matrici , i sumele din relaiile (2.31) i (2.42).

Un algoritm alternativ al acestei metode a fost prezentat de Kokotovic (1975), care pornete de la ecuaiile (2.31) i (2.42) ce sugereaz schimbarea de variabil nlocuind variabila z prin n relaia (2.17) se obine:

18

Rezult:

Unde cu M s-a notat matricea Conform relaiei (2.16) avem:

Pentru a decuple complet starea x de starea nct M=0. Vom avea:

n relaia (2.44), se allege matricea L astfel

Ceea ce indic faptul c L poate fi determinat iterativ: unde Sistemul (2.16), (2.17) dupa k iteraii are forma (2.26), (2.27) n care matricile sunt definite astfel:

Se observ c dac k Analog se face substituia

atunci cnd

este suficient de mic.

sugerat n relaia (2.42), n (2.27):

19

Rezult:

Unde cu N s-a notat matricea:

n continuare se exprim mrimea x n funcie de

care se

nlocuiete n relaia (2.27), i care va avea urmtoarea form:

Pentru a decupla

de , n relaia (2.46), se face N=0 i rearanjnd expresia se obine:

Relaia sugereaz c se poate determina H iterativ, ca fiind: unde Sistemul (2.26), (2.27) dup j-H iteraii are forma (2.40), (2.41) n care matricile au expresiile:

Determinarea iterativ a matricilor L i H evit inversarea repetat a matricilor i care sunt necesare relaiilor (2.25) i (2.26). Aceste relaii arat c matricile L

i H se obin dintr-o succesiune de aproximaii cvasistaionare. Pentru calcului matricilor , se vor utiliza n continuare matricile L i H determinate iterativ.

3. SEPARABILITATE N SCARA DE TIMP. EXEMPLIFICRI

20

3.1 INTRODUCERE

n capitolul anterior s-a prezentat o metod iterativ de separare a scrilor de timp bazndu-se pe faptul c variabilele de stare utilizate n descrierea modelului matematic sunt separabile, adic prin examinarea parametrilor fizici se poate determina care dintre mrimile de stare sunt rapide i care sunt lente. Se va observa c exist situaii n care setul de variabile de stare alese nu pot fi separate n acest fel. n cazul n care variabilele de stare alese iniial nu sunt separabile, se redefinesc variabilele de stare astfel nct modelul s poat fi scris sub forma (2.7), (2.8). n continuare se va prezenta problema separrii pe dou sisteme elementare componente ale unui grup energetic, i anume: un regulator de tensiune la care variabilele de stare ale modelului sunt separabile i un transformator ale crui mrimi de stare alese iniial nu sunt separabile, dar la care printr-o transformare adecvat se ajunge la un model cu variabilele separabile.

3.2 REGULATORUL AUTOMAT DE TENSIUNE

3.2.1

NECESITATEA I AVANTAJELE AUTOMAT DE TENSIUNE

INTRODUCERII

REGULATORULUI

Pentru a se asigura o repartiie de putere convenabil consumatorilor, concomitent cu o producie raional a energiei electrice este absolut necesar ca tensiunea la consumatori, precum i nodurile generatoare ale sistemului energetic, s fie ct mai constant posibil. n funcionarea normal a unui sistem energetic se produc permanent variaii ale sarcinilor active i reactive, care determin variaii de tensiune, att la barele consumatorilor de energie electric, ct i la nodurile de distribuie sau generatoare. Perturbaiile din reelele electromagnetice sunt de asemenea nsoite de importante variaii ale tensiunii. La scurtcircuitele din sistemul energetic se produce scderea nivelelor de tensiune pe barele centralelor i staiilor electrice de transformare sau interconexiune i se creeaz pericolul pierderii sincronismului n funcionare i stabilitii sistemului energetic. Deconectarea de ctre protecia prin relee a unor linii radiale de nalt tensiune
21

sau deconectarea unor consumatori importani sunt nsoite de importante variaii ale tensiunii la bornele generatoarelor sau ale consumatorilor. Obiectivul regulatorului automat de tensiune este de a menine constant sau cel puin ntre limitele prestabilite a tensiunii n nodurile generatoare, independent sau slab dependent de perturbaii. Meninerea constant a tensiunii aduce mari avantaje n exploatarea sistemului energetic i in funcionarea consumatorilor de energie electric, n particular. Astfel, regulatorul automat de tensiune prin reducerea creterii unghiului electric , mrete

stabilitatea static a sistemului energetic; crete totodat valoarea puterii maxime ce se poate debita n regim static de funcionare. De asemenea, regulatorul automat de tensiune mrete sensibil stabilitatea dinamic a funcionrii agregatelor din sistemul energetic i a sistemului n ansamblu. Datorit regulatorului automat de tensiune se evit producerea fenomenului de avalan de tensiune, caracterizat prin scderea continu i aproape total a tensiunii pe bare, ca urmare a unui defect sau regim anormal, nerezolvat n timp, ct i a ncercrilor nereuite de autopornire a motoarelor asincrone, alimentate la barele respective. n impedanele motoarelor care refuz s autoporneasc se produc cderi de tensiune suplimentare. Se tie c la variaia cuplului motoarelor asincrone prezint o

discontinuitate, obinndu-se puncte stabile de funcionare i astfel, autopornirea motorului asincron este asigurat. Reglarea automat a teniunii mai prezint i alte avantaje, cum ar fi: mrirea sensibilitii proteciei prin relee, n special n regim minim de scurtcircuit, precum i creterea indicilor de calitate ai energiei electrice produse i a fiabilitii instalaiilor i a sistemului energetic n ansamblu. n continuare, se prezint schema de reactane echivalent, pentru o singur faz, n cazul unui generator sincron G, care debiteaz pe barele B (figura 3.1).

xd

xs

22

Ed

Figura 3.1 SCHEMA DE REACTANE ECHIVALENT N CAZUL GENERATORULUI SINCRON G n aceast schem se utilizeaz urmtoarele notaii: Xd - reactana sincron longitudinal a generatorului innd seama de saturaie Xs - reactana elementelor de legtur dintre generator i bare Ed - tensiunea electromotoare eficace pe faz ( n spatele reactanei Xd ) V - tensiunea eficace pe faz pe barele B - unghiul electric intern ntre V i Ed Ig valoarea eficace a curentului debitat de generatorul G Pg , Qg - puterea activ, respectiv reactiv debitate de generatorul G Pc , Qc - puterea activ, respective reactiv consumat - factorul de putere la generator i consumator Diagrama fazorial a caderilor de tensiune corespunztoare schemei din figura de mai sus este reprezentat in figura 3.2: y

XIg Ed

23

Ig

Figura 3.2 DIAGRAMA FAZORIAL A CDERILOR DE TENSIUNE Considernd proieciile fazoriale pe axele Ox, Oz, se pot scrie ecuaiile:

unde:

Se observ c:

prin

s-a notat componenta reactiv a curentului generat.

Exprimnd:

se va obine:

Pentru meninerea constant a tensiunii pe barele generatorului sau a staiei electrice, exist dou posibiliti:
-

variia tensiunii electromotoare Ed ( Ed i generatorului sincron reglat;

variabile), prin variaia excitaiei

24

variaia reactanei de legtur X ( X si

variabile), prin variaia raportului de

transformare al transformatoarelor din staiile electrice, prevzute cu reglaj sub sarcin Cazul = constant nu se ia n considerare, deoarece complic mult problema reglrii

automate a tensiunii; reglarea se face ns unilateral, acionnd asupra lui Ed sau X i rezultnd corespunztor.

3.2.2

SEPARABILITATEA MRIMILOR DE STARE UTILIZATE N DESCRIEREA MATEMATIC A REGULATORULUI AUTOMAT DE TENSIUNE

Se consider un regulator automat de tensiune RAT care pentru reglarea tensiunii, sacioneaz asupra excitaiei generatorului sinscron. Schema standard a unui astfel de regulator este reprezentat n figura urmtoare:

AMPLIFICATOR

EXCITATIE

Vref

V Efd

VR

COMPENSATOR CU REACTIE

25

Rf

Figura 3.3 REGULATORUL AUTOMAT DE TENSIUNE

n figura de mai sus s-au folosit urmtoarele notaii: V tensiunea la bornele generatorului sincron Vref tensiunea de referi Efd tensiunea de excitaie KA factorul de amplificare a regulatorului de tensiune KE factorul de amplificare al excitaiei KF factorul de amplificare al compensatorului cu reacie Rf mrimea de stare a compensatorului cu reacie TA constanta de timp a regulatorului de tensiune TE constantade timp a excitaiei TF constanta de timp a compensatorului cu reacie VR tensiunea la ieirea regulatorului de tensiune SE ( Efd ) saturaia excitaiei

26

Se va utiliza saturaia excitaiei o relaie ce ine cont de neliniariti:

Valorile numerice ale constantelor sunt:

Utiliznd funciile de transfer ale componentelor, schema cu blocuri a regulatorului automat este:

Vref

Efd

27

Rf

Se calculeaz funciile de transfer ale buclelor: pentru amplificator:

pentru excitaie (nu se ia in considerare saturaia SE (Efd) ):

pentru bucla compensatorului:

28

Vref

VR +

Efd

Rf

Funciile de transfer ale buclelor se vor scrie astfel:

Din relaiile (3.4), (3.5), (3.6) rezult c:

29

de unde se obin ecuaiile de stare sub forma:

n care Pentru a separa mrimile de stare se va identifica buclele rapide i lente din figura 3.3 prin examinarea constantelor regulatorului automat de tensiune. Deoarece TA este mult mai mic dect TE (TA << TE) i cu mult mai mic dect TF bucla de amplificare este rapid, iar variabila de stare VR va fi n principal rapid. Din analiza parametrilor fizici nu rezult natura mrimii de stare are o parte rapid, aceasta nu va trece prin calea mare), ci prin calea paralel cu . Dac

( deoarece constanta de timp TF este i se

. De aceea nu se tine cont de calea cu

analizeaz sistemul rmas. n continuare, se va utiliza factorul KE =constant drept echivalent liniarizat al factorului de amplificare din excitatory. Schema sistemului este urmtoarea:

30

Vref

+ +

Efd

Funcia de transfer a buclei este:

n regim staionar avem:

Se observ c factorul de amplificare al buclei mare, iar

este mare deoarece

este

este mic. Rezult c Efd , ca semnal ntr-o bucl cu factor de

amplificare mare va avea o partid rapid. Mrimea rmas Rf va fi mrime lent.

Urmtorul pas este verificarea dac aceast alegere a variabilelor de stare poate fi adaptat relaiilor (2.7), (2.8). Pentru c fenomenele rapide sunt cauzate de valorile mici ale parametrilor , considerm astfel nct:

Trebuie alese scrile de timp astfel nct , adic Ca atare, se va utiliza

Se constat c

dac KA

ca o variabil de stare rapid. Vom avea:

Cu valorile numerice menionate, nlocuind n relaiile (3.10), (3.11), (3.12), modelul regulatorului de tensiune n forma (2.7), (2.8) este:

Starea cvasistaionar se obine pentru anterioare se determin modelul lent:

. Drept urmare, din relaiile

32

Pentru determinarea modelului rapid, se rescriu ecuaiile n scala de timp rapid i se examineaz la limit, cnd conform procedurii de reducere (2.11) (2.15)

prezentate n capitolul anterior. Se obine astfel modelul rapid sub forma:

3.3 TRANSFORMATORUL ELECTRIC 3.3.1. ELEMENTELE CONSTRUCTIVE ALE TRANSFORMATORULUI MONOFAZAT

Transformatorul monofazat este un aparat de current alternative care transform o putere electric alternativ, putere primar, de anumii parametri, n alta putere electric alternativ, putere secundar, parametrii fiind schimbai i frecvena rmnnd aceeai. Circuitele electrice ale aparatului ntre care are loc transferul de putere are n general numr diferit de spire, ceiu doi factori ai puterii, tensiunea i curentul sufer prin transformare schimbri inverse: daca tensiunea crete, curentul se micoreaz i invers. Funcionarea transformatorului se bazeaz pe legea induciei electromagnetice. n esen un transformator const dintr-un fier pe care sunt fixate dou nfurri izolate una de alta i fa de miezul de fier. nfurarea care primete energie de la surs este nfurarea primar, iar cealalt care cedeaz energie unei reele, este nfurare secundar. n studiul transformatorului, toate mrimile care se refer la nfurarea primar se numesc primare i se noteaz cu indicele 1, iar cele care se refer la nfurarea secundar se numesc secundare i se noteaz cu indicele 2. Principalele elemente constructive ale transformatorului sunt: sistemul magnetic, nfurrile, prizele de reglaj al tensiunii, izolatoarele de trecere i cuva de ulei. Sistemul magnetic al transformatorului reprezint asnamblul de poriuni feromagnetice pe care se nchid circuitele sale magnetice. nfaurrile transformatorului se execut din conductoare de cupru sau aluminiu electrolitic, a crui form i seciune depind de curentul nominal al transformatorului.

33

Numrul de spire al nfurrii primare i secundare se fixeaz n funcie de tensiunea nominal a transformatorului i de raportul de transformare. Cnd transformatoarele sunt instalate n puncte diferite ale reelei primare de distribuie, se poate ntmpla ca tensiunile alimentare n acele puncte, din cauza cderilor inegale de tensiune pe reea s fie diferite; nseamn c vor fi diferite i tensiunile sub care transformatoarele alimenteaz reeaua secuindar. Pentru a remedia acest neajuns, nfurrile de nalt tensiune ale transformatoarelor de putere sunt prevzute cu dou sau mai multe prize suplimentare de reglaj pe lng cea principal, care permit s se aplice nfurrii primare o tensiune a crei valoare s fie aproximativ diferit de cea

nominal. Izolatoarele de trecere asigur izolaia necesar buloanelor de born ale transformatorului, care fac legtura ntre nfurri i barele de nalt tensiune i de joas tensiune a instalaiei. Cuva de ulei asgur transformatorului o rcire bun dar i o izolaie perfect. Regimul nominal al transformatorului este regimul definit prin ansamblul valorilor mrimilor electrice sau de alta natur, nscrise pe plcua indicatoare a transformatorului i care caracterizeaz funcionarea n condiiile prescrise. Cu alte cuvinte, regimul nominal de funcionare este regimul de sarcin pentru care a fost proiectat i construit transformatorul. Pentru transformatoarele de putere n ulei, funcionarea n regim nominal este fixat de urmtoarele date nominale: puterea, tensiunile, raportul de transformare, cureni, tensiunea de scurtcircuit i frecvena, n condiiile de rcire ale transformatorului. La transformatoarele cu mai multe prize, regimul nominal este acela care corespunde prizei nominale (priz nominal).

3.3.2 SEPARAREA SCALELOR DE TIMP APLICAT TRANSFORMATORULUI MONOFAZAT Se consider schema monofazat a unui transformator reprezentat astfel:

34

R1 N1 V ii N2

R2

i2

Figura 3.4 TRANSFORMATORUL MONOFAZAT

Pentru obinerea ecuaiilor generale de funcionare ale transformatorului se folosesc inductanele totale, presupuse a fi constante, ale cror dou circuite ale transformatorului se gsesc ntr-n mediu omogen, de permeabilitate constant ( ). Se noteaz cu L1 inductana totala primar, cu L2 indcutana total secundar, cu M indcutana mutual a celor dou nfurri i cu l1 , l2 inductanele de scpri ale celor dou nfurri ale transformatorului. Se consider c raportul dintre inductanele de dispersie l1 , l2 i inductanele proprii L1 , L2 este egal cu acelai parametru mic :

Fluxurile de legtur

se vor exprima astfel:

n care i1 , i2 reprezin curenii din nfurarea primar i secundar. innd cont de relaia (3.20), ecuaiile (3.21) se scriu:
35

Se cunoate c inductanele de dispersie sunt:

ceea ce conduce la:

nlocuind aceste relaii n ecuaiile fluxurilor (3.22) se obine:

Matricial, sistemul se scrie:

de unde se exprim curenii

36

adic:

Notm cu:

Ecuaiile (3.28) devin:

Aplicnd celor dou circuite teorema a doua a lui Kirchhoff, se obin relaiile:

nlocuind relaiile pentru cureni, modelul transformatorului devine:

Pe baza notaiei, avem:

37

de unde rezult:

adic:

din care:

Prin urmare, ecuaiile (3.32) devin:

n continuare, se introduc urmtoarele notaii:

care nlocuite n relaiile (3.34) conduc la:

38

Se obine c parametrul mic ambele mrimi de stare avem:

nmulete ambele derivate i deci sunt rapide. Astfel, de exemplu, dac lum

n care ambele derivate vor tinde la infinit Analiznd modelul (3.35), scris sub forma:

atunci cnd

se observ c matricea sistemului:

care are determinantul:

Devine o matrice singular lent ascuns.

, ceea ce indic posibilitatea existenei unui fenomen

Deoarece ambele variabile de stare sunt rapide, se va forma starea lent eliminndu-se fenomenele rapide din
39

Acest lucru este echivalent cu a defini

starea x ca o combinaie liniar de cu parametrul .

astfel nct derivata lui x s nu fie multiplicat

n cele ce urmeaz se va determina o transformare liniar x adecvat, iar pentru simplificare se va aproxima:

adic:

nmulind ultima ecuaie cu h i scznd-o din prima, se obine:

de unde rezult c:

Realiznd schimbarea de variabil

se determin c:

Am gsit transformarea de variabil

astfel nct noua variabil x este o

variabil lent, deoarece derivata sa nu este multiplicat cu . Se observ c variabila care se va nota z, z= se pstreaz ca variabil rapid.

40

Prin urmare, modelul iniial (3.35) devine:

Noul model este un model de tipul (2.7), (2.8), al crui model lent se determin conform procedurii prezentate n capitolul anterior, adic la limit cnd

de unde:

Prin urmare avem:

Modelul lent este:

a crui semnificaie fizic este aceea de fluxul de legtur al tranformatorului ideal:

iar modelul rapid va fi:

41

representnd fluxul pierderilor de dispersie. Astfel s-a demonstrat c printr-o alegere adecvat a variabilelor de stare, modelul matematic poate fi adus la forma (2.7), (2.8), adic la un model cu variabile de stare lente i rapide.

4. MODELAREA DINAMICII UNUI GRUP ENERGETIC CUPLAT LA SISTEMUL ENERGETIC

4.1 NOIUNI GENERALE DESPRE MAINA SINCRON

Maina sincron este una dintre cele mai rspndite maini electrice ce funcioneaz n curentul alternativ la o vitez constant, numit vitez de sincronism. ntre acestea, frecvena f a relaiei i numrul p de perechi de poli ai mainii, exist relaia: Dup destinaia sa, maina sincron poart numele de generator sincron sau alternator, motor sincron i compensator sincron. Maina sincron are un domeniu foarte larg de utilizare ca geenrator sincron de curent alternativ trifazat de frecven industrial, deoarece n prezent energia electric se produce aproape n exclusiviate sub form trifazat. Maina sincron se utilizeaz ca motor sincron atunci cnd se cere o turaie constant, motorul sincron prezentnd n
42

exploatare urmtoarele avantaje:

ridicat, siguran n funcionare i economie. n

plus, maina sincron poate fi utilizat ca generator de putere reactiv atunci cnd funcioneaz n regim de motor supraexcitat. Acest lucru conduce la ameliorarea factorului de putere al reelei la care este cuplat. Dup forma rotorului, mainile sincrone se mpart n dou categorii : cu poli apareni (ntrefierul neuniform) i cu poli necai (ntrefierul uniform). Mainile de turaie lent sau mijlocie ( <1500 rot/min) se execut cu poli apareni, iar cele de mare turaie (3000 rot/min) precum i cele de 1500 rot/min din unitile mari, cu poli necai. Generatoarele sincrone moderne sunt acionate de regul de turbine cu abur sau cu gaze, motoare Diesel i de turbine hidraulice. Turbinele termice sunt maini de turaie mare; generatoarele sincrone antrenate de ele se execut cu poli necai i poart numele de turbogeneratoare. Motoarele Diesel sunt maini cu turaie mic i mijlocie, iar generatoarele antrenate de acestea se numesc generatoare pentru grupuri Diesel electrice i se execut cu un numr relativ mare de poli apareni. Turbinele hidraulice fac parte, n general din categoria mainilor cu turaie redus, iar generatoarele respective, numite hidrogeneratoare se construiesc cu un numr relativ mare de poli apareni i cu arbore vertical, la fel ca i turbina care le antreneaz. Statorul mainii sincrone cu poli apareni const dintr-un miez de fier de form cilindric, care constituie partea sa activ i poart nfurarea indus dintr-o carcas metalic, care fixeaz maina de fundaie i susine miezul din scuturi, care mpreun cu carcasa i miezul creeaz sistemul de camere i canale de ventilaie utilizat n rcirea mainii. Rotorul cu poli apareni este construit dintr-un butuc de font sau oel, montat pe arbore i din polii inductori, fixai pe butuc prin buloane. Polii sunt formai din electromagnei excitai n curent continuu. Bobinele polilor sunt legate n serie astfel nct s se magnetizeze succesiv unul nord i altul sud. Oricare ar fi tipul constructiv al rotorului, pentru a alimenta cu curent continuu nfurarea acestuia n timp ce rotorul se nvrtete, capetele nfurrii rotorice sunt legate la dou inele colectoare din bronz sau din font, fixate pe axul mainii i izolate ntre ele i fa de ax. Pe inele alunec dou perii care aduc curentul de la un mic generator de curent continuu, numit excitator, montat de obicei pe acelai arbore cu rotorul. Atunci cnd generatorul este conectat pe partea statoric la o reea trifazat simetric, nfurrile statorice vor prezenta un sistem trifazat simetric i echilibrat de cureni, succesiunea acestora fiind identic cu cea a tensiunilor electromotoare induse de cmpul nvrtitor de excitaie.

43

n funcionarea normal a generatoarelor sincrone pot surveni o serie de defecte, ditnre care cele mai frecvente sunt scurtcicuitele. Exist dou regimuri de scurtcircuit: regimuri staionare i regimuri tranzitorii. Se poate pune un generator sincron n regim staionar de scurtcircuit dac i se reduce treptat pn la zero impedana circuitului exterior, meninndu-i-se constant turaia i curentul de excitaie. Dac scderea impedanei circuitului exterior are loc brusc, aa cum se ntmpl de obicei n exploatare, scurtcircuitul se numete scurtcircuit brusc sau instantaneu i n funcionarea alternatorului apare regimul tranzitoriu de scurtcircuit. Prin urmare, regimul tranzitoriu apare numai n cazul scurtcircuitelor brute i precede regimul permanent, stabilizat de scurtcircuit. 4.2 ECUAIILE REGIMULUI TRANZITORIU AL GENERATORULUI SICRON

Pentru studiul regimurilor tranzitorii ale mainii sincrone se utilizeaz pe scar larg ecuaiile Park-Blondel, care rezult din ecuaiile n mrimi instantanee reale ale diferitelor circuite statorice i rotorice dup o transformare adecvat. n scopul deducerii acestor ecuaii, se va utiliza un generator sincron cu poli apareni care posed o colivie de amortizare plasat n piesele polare cu inele de scurtcircuitare complete, astfel nct aceasta colivie acioneaz att dupa axa longitudinal, ct i dup cea transversal. Se presupune c sunt variabile ipotezele teoriei fizice i anume: - circuitul magnetic al mainii are o caracteristic liniar de magnetizare, se neglijeaz saturaia i fenomenul de histerezis, iar intesitatea cmpului magnetic n prile de fier ale circuitului magnetic se consider nul; - nu se iau n considerare dect armonicele fundamentale ale undelor densitilor liniare ale soluiilor de excitaie i de reacie; - se admite o simetrie perfect constructiv, magnetic i electric; - viteza la care este antrenat generatorul se consider constant.

4.2.1

ECUAIILE GENERATORULUI SINCRON N MRIMI INSTANTANEE REALE

Se consider schema principal a unui generator sincron cu poli apareni din figura 4.1, care evideniaz mrimi absolut necesare descrierii matematice. Posed trei nfurri statorice de faz, notate a, b, c, decalate uniform la periferia interioar a statorului. Rotorul are nfurarea de excitaie f, axa unui pol nordic rotoric fiind denumit axa longitudinal a rotorului. Tot pe rotor se va lua n considerare dou nfurri n scurtciruit de amortizare. Prima, notat cu g acioneaz dup axa
44

longitudinal ca i nfurarea de excitaie. A doua nfurare de amortizare caracterizat prin indicele h acioneaz dup o ax de simetrie interpolar situat imediat nainte, n sensul micrii rotorului, n raport cu axa longitudinal aleas. Decalajul electric dintre acele axe este de 90 . Fie o vitez unghiular instantanee a rotorului, considerate n sensul invers micrii acelor de ceasornic. Unghiul fcut la un moment oarecare t de axa longitudinal a rotorului cu axa de simetrie a unei bobine a fazei a, considerate ca referin spaial a fost notat cu . Desigur, i

45

Sistemul rotator legat de rotor (d,q,o) se poate raporta n general la sistemul rotitor de vitez unghiular i care este notat cu ( . Se poate obine astfel

generalizarea caracterizrii si a unui ansamblu de generatoare. Cazul din figur corespunde cazului n care exist o singur pereche de poli (p=1) i pulsaia tensiunii electice produse este egal cu viteza unghiular mecanic cazul general cnd generatorul are p perechi de poli vom avea: unde cu s-a notat viteza unghiular mecanic a rotorului ( care este fixat de turbina pulsaia tensiunii electrice. n

de antrenare), iar cu

n scrierea ecuaiilor funcionale se va adopta convenia de sensuri pozitive ale mrimilor variabile de la circuitul generator pentru nfurrile statorice i convenia de la circuitul receptor pentru nfurrile rotorice. Ecuaiile lui Ohm pentru rotor i stator sunt:

46

unde:

reprezint matricile corespunztoare rezistenelor, tensiunilor i curenilor ce caracterizeaz rotorul, respectiv statorul sincron. De obicei, rezistenele r ale statorului au valori mici i de aceea acestea n general se neglijeaz. Fluxurile magnetice prezente n ecuaiile (4.2), (4.3) se exprim astfel:

n care inductanele proprii i mutuale sunt reprezentate prin matricile:

47

Unde:

Datorit structurii circuitului magnetic i dispunerii spaiale a nfurrilor, elementele prezentate n relaiile (4.9) primesc scrierea:

unde:

i n general:

respectiv:

48

Se poate observa c partea electric a generatorului sincron se descrie prin ecuaii difereniale, dar n care coeficienii sunt funcii de timp deoarece:

i de exemplu pentru: = constant rezult c:

Ca urmare, sistemul dynamic ce i corespunde este un sistem variant acest neajuns se nltur apelnd la o schimbare adecvat, cea mai cunoscut fiind cea propus de Park. Pentru a putea descrie matematic i partea mecanic, se menioneaz c deducerea cuplului electromagnetic rezistent se face cu relaia:

49

4.2.2 TRANSFORMAREA PARK. ECUA IILE FUNCIONALE N SISTEMUL DE REFERIN ROTATIV (

Integrarea ecuaiilor funcionale difereniale n mrimi relative instantanee ale generatorului sincron cu poli apareni, n diferite situaii concrete, este dificil din cauza inductivitilor proprii ale circuitelor statorice i a celor mutuale dintre circuitele statorice i rotorice care sunt variabile n funcie de timp. Datorit acestui caracter variat se prefer exprimarea n coordonate rotative n care dispare acest neajuns; transformarea respectiv se numete transformarea Park i n cele ce urmeaz ea este facut fa de ssitemul rotativ cu viteza unghiular Aceasta

transformare se aplic tuturor mrimilor statorice i rotorice i este exprimat sub forma a dou matrici: - pentru stator:

pentru rotor:

Prin aceasta schimbare de variabil, mrimile statorice de faz, notate cu indice a,b,c vor fi nlocuite cu alte mrimi avnd indicii Sd, Sq, So. Uzual se notez:

50

Pentru a determina invervesele matricilor transpusele acestora:

se calculeaz determinanii i

rezult:
51

Din ecuaiile (4.2), (4.3) cu transformrile (4.18), (4.19) se determin:

din care se obine:

Calculnd, se obin urmtoarele rezultate:

52

Se noteaz:

unde:

n urma acestor notaii, ecuaiile (4.25) se scriu:

Aceste ecuaii reprezint o prim form a ecuaiilor funcionale, pe partea electric a generatorului sincron. Acestora li se asociaz i ecuaiile (4.7), (4.8), trecute n coordinate rotative ( . Astfel, ecuaiile fluxurilor magnetice (4.20) devin:

n care se vor folosi notaiile:

53

(4.32) cu:

De asemenea:

cu:

n aceste condiii relaiile dintre fluxuri i cureni, scrise n coordonatele (

sunt:

54

Cuplul electromagnetic n coordonatele ( i (4.17), transformrile (4.18) i (4.19).

se obine plecnd de la relaiile (4.16)

Pentru evaluarea relaiei (4.38) se determin:

55

de unde:

n relaia (4.43) s-a notat:

unde: reprezint inductana sincron longitudinal, respectiv transversal i care se pot determina astfel:

Continund:

56

ntruct

sunt mrimi scalare i tinnd cont de relaia (4.9) rezult:

Deoarece Rezult c: ntr-o scriere compact cuplul electromagnetic devine :

Evaluarea parametrului m este necesar n cazul n caracterizrii complete a dinamicii generatorului sincron, care pe partea mecanic va fi descris de:

unde

este cuplul dezvoltat de turbinele de antrenare. Cnd se analizeaz numai

reglajul de tensiune (putere reactiv), deoarece dinamica indus de (4.51) este mult mai lent, se poate considera c =const., adic =const.

4.2.3

CAZUL PARTICULAR AL RAPORTRII LA UN REFEREN IAL LEGAT DE ROTOR

57

Se consider: = (constant) (4.52)

i se analizeaz elementele definite anterior: Din relaia (4.19) rezult:

Iar din (4.23):

n acest caz pe cele dou axe ale rotorului vor fi dispuse urmtoarele nfurri: f- nfurarea de excitaie g - nfurarea de amortizare d - raportarea circuitului statoric la rotor pe axa q: h nfurarea de amortizare q raportarea circuitului statoric la rotor Transformarea (4.53) are ca efect:

58

De asemenea, din relaia (4.27) rezult:

Ecuaiile corespunztoare statorului, n conformitate cu relaiile (4.29) se scriu:

iar cele din rotor:

Din relaiile (4.58) se deduc ecuaiile echivalente:

Ultimele ecuaii din (4.59) se egaleaz cu 0 pentru c amortizrile g i h la generatorul sincron sunt n scurtcircuit. Ecuaiile de legtur ntre fluxuri i cureni se deduc dup particularizarea ecuaiei (4.52) aplicat relaiilor (4.31), (4.32), (4.34) i (4.35). n urma acestei particularizri, expresiile inductanelor proprii i mutuale devin:
59

Conform relaiilor anterioare i innd seama c:

Rezult pentru stator:

60

iar pentru rotor:

Relaiile (4.57) i (4.59) mpreun cu (4.65) i (4.67) reprezint ecuaiile funcioanle ale generatorului sincron, pe partea electric a acestuia, n cazul n care intereseaz numai reglajul de tensiune. Aceasta se realizez cu ajutorul regulatorului automat de tensiune care menine constant tensiunea n nodurile generatoare ntr-un mod independent sau slab dependent de perturbaii. Pentru eventuala extindere a modelului i pe partea electromagnetic se evalueaz i cuplul electromagnetic rezistent. Plecnd de la relaia (4.50), n ipoteza (4.52) se obine:

Dup efectuarea calculelor n relaia (4.68) se obine:

iar dac se ine cont de (4.65) rezult:

61

Relaie fundamental pentru evaluarea cuplului electromagnetic rezistent, valabil i n cazul regimului sincron menionat prin relaia (4.52). O alta form de exprimare a cuplului m care s implice numai curenii se poate deduce din (4.68) dac vectorul prezent i se aplic transformarea:

a crui invers este:

n relaia (4.68) se noteaz:

Discutnd despre un singur generator sincron, condiia (4.52) poate fi interperetat ca reprezentnd raportarea la un sistem rotitor legat rigid de rotorul mainii; n acest caz, relaiile funcionale se simplific mult i coeficineii nu mia depind de timp. n cazul mai multor generatoare sincrone cuplate sau chiar n cazul unui singur generator fa de un sistem de putere infinit, alegerea sistemului de referin comport dou variante:
62

raportarea fiecrei maini la sistemul propriu legat de rotorul propriu al mainii; raportarea tuturor la sistemul rotativ de vitez unghiular sincron (al sistemului de putere infinit).

Prima manier are avantajul c simplific ecuaiile funcionale ale fiecrei maini i de aceea este utilizat frecvent. Trebuie specificat c toate consideraiile fcute pn acum sunt valabile i pentru cazul turbogeneratoarelor sincrone, la care, datorit absenei polilor apareni, se ndeplinesc condiiile:

Urmrind ecuaiile (4.57), (4.59), (4.56) se constat c componentele homopolare ale statorului nu intervin n celelalte ecuaii, n cazul raportrii la sistemul rotativ propriu. De aceea, modelul funcional se scrie compact astfel:

Ultima ecuaie din (4.75), datorit datelor constructive ale statorului, se poate scrie:

63

n ecuaiile funcionale (4.75), ce se rescriu imediat sub form de sistem dinamic liniar, mrimea de stare este format din fluxuri sau cureni, mrimea de execuie este tensiunea de excitaie , iar sunt mrimi intermediare ce se vor elimina cu

ajutorul ecuaiilor ce corespund liniei de cuplare a generatorului sincron la reea.

4.3 ECUAIILE FUNCIONALE ALE GENERATORULUI SINCRON N MRIMI RELATIVE

Mrimile ce caracterizeaz un sistem (tensiune, curent, putere activ sau reactiv etc) pot fi calculate mult mai uor dac sunt reprezentate n uniti relative. Valoarea unitii relative, pentru oricare dintre mrimi, se calculeaz ca fiind raportul dintre mrimea respectiv i o mrime de baz de aceeai natur fizic, arbitrar aleas. n alegerea condiiilor de baz se impune ca operaiile de calcul s fie ct mai simple iar ordinul de mrime al unitilor relative s permit utilizarea comod a acestora. Ecuaiile generatorului sincron sunt scrise n general utiliznd mrimi relative. Aceasta se face cu ajutorul mrimilor ce corespund regimurilor nominale n felul urmtor:
1. Toate tensiunile din stator se impart (relativizeaz) prin tensiunea maxima de

regim nominal de pe faza generatorului

, eventual tensiunea de la

transformatorul de cuplare la linia de nalt tensiune:

O tensiune relativizat (oricare dintre

etc) este:

2. Toi curenii din stator se relativizeaz prin

care este curentul nominal

maxim pe faz.

64

Relativizarea unui curent de stator ( relaia:

etc) se face cu

3. Puterea se relativizeaz prin:

4. Rezistenele i reactanele se relativizeaz prin:

care se pot determina astfel:

5. Timpul se relativizeaz astfel:

i atunci timpul relativizat este:

De asemenea:

65

i pulsaia relativ este:

Se observ c:

6. Cuplul se relativizeaz prin:

iar mrimea relativizat va fi de forma:

7. Pentru fluxuri i inductiviti mrimea cu care se face relativizarea se deduce astfel:

de unde rezult:

de unde fluxul relativ va fi:

66

Corespunztor rezult:

Scriind:

se determin

i deci:

Deoarece reactana relativ este:

n cazul n care

se boserv c:

care este o relaie cunoscut i anume, numai n sincronism, reactanele i inductanele conduc la aceeai relativ 8. Momentul de inerie conduce la:

Iar mrimea relativ va fi:

Mrimile de baz ale circuitului de excitaie sunt: relaia:

n care exist

n mod asemntor se aleg mrimile de baz ale circuitelor de amortizare, ntre aceste mrimi existnd relaiile:
67

Cu acestea, mrimile raportate se vor defini astfel:

Expresiile fluxurilor se simplific dac n relaiile (4.105) se impun condiiile:

68

care sunt variabile numai dac:

Prin urmare, valorile impedanelor legate de mrimile

nu sunt independente, ci trebuie

Modelul funcional al generatorului sincron n mrimi relative se poate scrie astfel:

Dac se aleg curenii de baz permit simplificarea relaiilor (4.112), se obine:

ntr-o form particular care s

69

Pe partea electromagnetic:

Prin urmare, modelul generatorului sincron n mrimi relative se obine nlocuind relaiile simplificatoare (4.113) n ecuaiile (4.112), dar nu se vor mai specifica mrimile relative prin asterisc i nu se vor mai nota cu indicii S mrimile statorice raportate la cele 2 axe d, q:

n care:

unde:
-

reactana de dispersie a unei faze a indusului reactana corespunztoare fluxului longitudinal reactana corespunztoare fluxului transversal

70

reactana sincron longitudinal reactana sincron transversal

Pe partea electromagnetic:

4.4

MODELUL MATEMATIC AL GENERATORULUI EVIDENIEREA TENSIUNILOR ELECTROMOTOARE

SINCRON

PRIN

Ipotezele care stau la baza scrierii ecuaiilor funcionale ale generatorului sicnron sunt: considerm rezistena statorului ca fiind nul: r = 0 generatorul sincron are amortizare numai pe axa q

Din demonstraiile anterioare, modelul generatorului sicnron este:

71

innd cont de ipotezele anterioare, modelul se poate rescrie sub forma urmtoare:

nmulind prima ecuaie din (4.19) cu cont de relaia (4.120), se obine:

, cea de-a doua cu

i innd

care se rescriu:

Notnd:

72

rezult:

Pentru circuitele statorice ale generatorului sincron, regimul normal de funcionare este regimul permanent sinusoidal, simetric i de succesiune direct. Regimul care se realizeaz n mod frecvent n exploatarea mainii sincrone, este regimul tranzitoriu cvasisinusoidal al circuitelor simetrice. Regimul permanent sinusoidal n circuitele a, b, c implic regimul staionar al circuitelor d, q ( derivatele i curenii prin amortizri sunt egale cu 0 ). Se obin astfel:

Regimul tranzitoriu cvasisinusoidal n circuitele a, b, c implic anumite aproximaii n ecuaiile difereniale ale circuitelor d, q. Aceste aproximaii sunt: se neglijeaz tensiunile electromotoare induse prin transformatoare fa de

tensiunile electromotoare induse prin rotaie. n acest fel ecuaiile difereniale ale circuitelor statorice devin ecuaii algebrice, la fel ca n regimul permanent sinusoidal, ns ecuaiile difereniale ale circuitelor rotorice rmn neschimbate. Ecuaiile statorice n regim tranzitoriu se scriu:

n care:

73

reprezint reactanele tranzitorii reprezint tensiunile electromotoare tranzitorii Din ecuaiile (4.134), (4.135), (4.136), (4.137) rezult:

innd cont c:

se obine:

Din ecuaia (4.132), utiliznd ecuaia (4.143) i din (4.133) utiliznd (4.122) se obine:

74

Din aceste ultime patru ecuaii, se obine expresia modelului generatorului sincron cu evidenierea tensiunilor electromotoare induse. Aceasta este:

5. ANALIZA REGIMULUI DINAMIC AL UNUI GRUP ENERGETIC PRIN METODA PERTURBAIEI SINGULARE

75

5.1

INTRODUCERE ntr-un capitol anterior s-a prezentat metoda perturbaiei singulare, introdus

pentru a elimina inconsistenele aproximarii clasice cvasistaionare i s-a propus o procedura de separare a scalelor de timp ce folosete o schem iterativ i care reprezint o cale de mbunatire a metodelor reduse. Se va analiza modelul de ordin 7 al unui generator ce debiteaz pe bare de putere infinit. Acest model surprinde dinamicile fluxului de legturi, a prtii mecanice a sistemului de excitaie i prin scrierea acestui model se va utiliza modelul folosit n descrierea regulatorului automat de tensiune, prezentat in capitolul 3, precum i modelul generatorului sincron din capitolul anterior. Acest exemplu va ilustra cteva caracteristici ale procedurii de separare i va introduce scalele de timp, iar prin aplicarea schemei iterative se valideaz alegerea variabilelor lente sau rapide i se mbuntete precizia subsistemelor lente i rapide.

5.2

SCHEMA GRUPULUI ENERGETIC CUPLAT LA SISTEMUL

ENERGETIC. PARAMETRII TEHNICI I SISTEMICI AI GRUPULUI PRECUM I AI LINIEI DE CUPLARE Se analizeaza cazul unui singur generator sincron cuplat la retea a crui schem este prezentat n figura 5.1. Generatorul este prevzut cu un regulator automat de tensiune, a crei schem a fost prezentat in capitolul 3 (figura 3.3.) i al crui model este descris prin ecuaiile (3.10), (3.11), (3.12). n plus, fa de funcia ce descrie saturaia excitaiei, se consider o exponentiala pentru a ine seama de neliniariti:

Se consider c la un moment dat apare un scurtcircuit trifazat ce dureaz cinci cicluri pe linia a, n apropierea barei 2, caz n care protecia declaneaz ntreruptoarele liniei a, care va fi astfel scoas din funciune.

76

FIG 5.1 SCHEMA GENERATORULUI SINCRON CE DEBITEAZ PE BAZA DE PUTERE INFINIT

Mrimile prezentate n aceast schem sunt: - tensiunea de excitaie a generatorului - tensiunea de referin a regulatorului automat de tensiune - tensiunea pe bara 3 - unghiul intern al mainii - tensiunea electromotoare tranazitorie longitudinal - tensiunea electromotoare tranazitorie transversal - reactana tranazitorie reactana liniei de transport n afara acestor mrimi n descrierea matematic a sistemului mai apar i alte notate cu:

77

D - amortizarea mainii H - constanta de inerie a mainii - constanta de timp longitudinal a circuitului deschis - constanta de timp transversal a circuitului deschis - reactana sincron longitudinal - reanctana sincron transversal - reactanele liniilor a i b Valorile nominale ale constantelor sunt prezentate n tabelul 5.1 TABEL 5.1. DATELE MAINII SINCRONE

La acestea se adaug constantele regulatorului de tensiune grupate n tabelul 5.2 TABELUL 5.2. CONSTANTELE REGULATORULUI DE TENSIUNE

5.3 REPREZENTAREA SISTEMIC A GRUPULUI ENERGETIC ANALIZAT

78

Modelul matematic al sistemului va include ecuaiile electromecanice, ecuaiile fluxlui de legtur i ecuaiile regulatorlui de tensiune. Tensiunile induse n axele d i q de ctre fluxul de legtur sunt:

n care prin Y s-a notat admitana liniei de transmisie.

Ecuaiile electromecanice ale sistemului sunt:

Ecuaiile regulatorului de tensiune au fost prezente n capitolul 3, acestea fiind:

De reinut faptul c tensiunea generatorului V rezult din diagrama fazorial i este definit prin relaia:

unde:

Ecuaiile (5.1), (5.2), (5.3), (5.4), (5.5), (5.6) i (5.7) constituie modelul matematic de ordinul 7 al mainii ce debiteaz pe bare de putere infinit, variabilele de stare ale acestui model fiind:
79

Se observ c modelul obinut este un model neliniar, neliniaritile fiind introduse de prezena funciilor sin i cos, precum i de saturaia excitaiei. Se calculeaz admitana liniei de transmisie Y:

de unde se obine:

Ca s determinm puterea

se va face un calcul de regim staionar, iar pentru aceasta

se va ine cont de configuraia reelei din figura 5.1, precum i de datele acesteia:

unde:

rezult:

de unde se determin puterea activ i reactiv a sistemului:

80

Pentru modelul determinat (5.1), (5.2), (5.3), (5.4), (5.5), (5.6), (5.7), se nlocuiesc valorile constantelor prezentate n tabelul 5.1, 5.2 i puterea activ calculat anterior i se obine:

n continuare se liniarizeaz modelul:

81

Se calculeaz valorile de regim staionar ale lui Din diagrama fazorial se determin c:

i .

Se cunoate:

de unde rezult c nlocuind, se obine:

Se cunoate c:

unde:

82

Deci:

adic:

nlocuind valorile de regim staionar, modelul matematic liniarizat (5.16), (5.17), (5.18), (5.19), (5.20), (5.21) i (5.22) se sscrie:

Mai departe, pentru separarea variabilelor de stare, se va ine seama de consideraiile fcute anterior la modelarea regulatorului dinamic de tensiune, care ilustreaz faptul c variabilele cel mai lent. Analiznd ecuaiile fluxului de legtur i valorile constantelor de timp, se observ c axa transversal are o constant de timp axa longitudinal variabil lent. Mrimile de stare rmase , respectiv decide dac aceast ipotez este corect. Prin urmare, se consider , ca fiind lente, rmnnd 5 variabile rapide. , se consider rapide, iar ulterior se va . Astfel, se alege , mult mai mic dect drept variabil rapid, iar drept i constituie subsistemul rapid, iar subsistemul

Reordonnd ecuaiile modelului cu separarea de variabile propuse:

mrimi lente mrimi rapide

Modelul liniarizat se rescrie n matricea:

Valorile

proprii

ale

matricii

sistemului

se

determin

ca

fiind:

, ceea ce confirm existena a dou mrimi de stare lente i rapide i 5 mrimi de stare rapide.

5.4 APLICAREA PROCEDURII ITERATIVE DE SEPARARE PE STRI

Schema iterativ descris n capitolul 3 se bazeaz pe modelul matematic scris n forma strilor separabile:

n care x este predominant lent, iar z surprinde fenomenele trnazitorii rapide. Parametrul de scalare present explicit n cadrul metodei, va fi n continuare nglobat n matricile C i D ale sistemului. Pentru cazul de fa vectorii x i z sunt:

84

Matricile A, B, C, D se obin din partiionarea matricei sistemului, care se va nota cu M, astfel nct:

ntregul sistem se descompune ntr-un subsistem lent:

i un subsistem rapid:

unde k i j indic numrul iteraiilor operate asupra variabilelor lente i rapide. Se reia procedura iterativ prezentat n capitolul respectiv, ns acum fr explicitarea parametrului . Se calculeaz nti pornind de la:

85

Se obine:

Se calculeaz

din:

Se obine:

Iteraiile (5.37) i (5.38) sunt folosite pentru a corecta variabilele rapide timp ce iteraiile (5.39), (5.40) se utilizeaz pentru corecia variabilelor lente

, n

Alegerea numrului de pai j i k depinde de cerinele preciziei variabilelor lente i rapide. Astfel, dac este necesar o exactitate mai bun a variabilelor lente, n timp ce micile inexactiti ale variabilelor rapide pot fi tolerate, atunci este necesar un numr mnai mare de iteraii pentru calculul matricei H i invers. Unul din avantajele metodei l constituie faptul variabilele a variabilelor originale x, n timp ce originale z. n plus, matricile de calcul rein natura fizic

sunt de aceeai natur fizic cu variabilele i sunt funcie de matricile originale A, B,

C, D, ceea ce permite n continuare un calcul simplu.

86

5.5 ELEMENTE NECESARE SIMULRII DINAMICE A GRUPULUI ENERGETIC PRIN METODATA PERTURBAIEI SINGULARE

Simularea dinamic a grupului energetic va pune n eviden faptul c modelul corelat numai dup o singur iteraie aproximeaz mult mai bine modelul exact, dect modelul aproximrii de stare cvasistaionar. Cu ajutorul programului de simulare dinamic, se vor trasa variaiile n timp ale mrimilor de stare, n conformitate cu cele tri modele studiate, ceea ce se reduce la reprezentarea grafic a soluiilor sistemului de ecuaii:

n care vectorul y i matricea N se particularizeaz n funcie de modelul studiat. Se cunoaste c soluia sistemului (5.41) este de forma:

iar pentru trasarea graficului se va aplica un algoritm de discretizare al crui pas h se alege astfel nct:

n care Rezult c:

reprezint valorile proprii ale matricei N.

Ca atare, soluia discretizat va fi:

Se introduce notaia: i obinem:

87

n cazul modelului exact, vectorul y i matricea N vor fi:

Modelul aproximrii de stare scasistaionar se obine din particularizarea matricii M a sistemului n conformitate cu variabilele lente i rapide:

Din care pentru modelul lent se reine numai influena modurilor lente, iar pentru modelul rapid numai cea a modurilor rapide, adic:

88

n acest caz se va aplica algoritmul de discretizare prezentat anterior att sistemului (5.44) ct i ssitemului (5.45), determinndu-se astfel variaia n timp a mrimilor de stare rapide i lente ale modelului cvasistaionar, cu precizarea c pasul de discretizare k n cazul modelului lent, se alge astfel nct

n care

reprezint valorile proprii ale matricii A.

Modelul iterativ, prezentat n paragraful (5.4) i exprimat prin relaiile (5.34), (5.35), (5.36) va fi ilustrat n cazul unei singure iteraii k = 1, j = 1, procedura de discretizare fiind aplicat matricilor i prezentate n cazul modelului respectiv.

89

6. PROGRAM NUMERIC DE SIMULARE DINAMIC A ENERGETIC PRIN METODA PERTURBAIEI SINGULARE

UNUI

GRUP

clear all clc disp( ' SIMULAREA DINAMICA A UNUI GRUP ENERGETIC') disp(' Matricea sistemului') M=[-0.58 0 0 -0.269 0 0.2 0 0 -1 0 0 0 1 0 0 0 -5 2.12 0 0 0 0 0 0 0 314 0 0 -0.141 0 0.141 -0.2 -0.28 0 0 0 0 0 0 0 0.838 2 -173 66.7 -116 40.9 0 -66.7 -16.7] disp('Valorile proprii ale matricii M:') Lambda_m=eig(M) %Matricile pentru modelul lent si rapid in regimul cvasistationar: As=M(1:2,1:2); Ds=M(3:7,3:7); %Variatiile initiale ale marimilor initializate: X0=[-0.01;0.02;0;0.12;0;0;2]; X01=X0(1:2); X0r=X0(3:7); %Pasul de discretizare: h= 0.08; %discretizarea matricii M si a celor de regim cvasistationar: Md=expm(M*h); k=0.32; Asd=expm(As*k); Dsd=expm(Ds*h); Xd(1:7,1)=eye(7)*X0; for j=1:1:7; for i=2:1:101; Mdi=Md^i; Xd(j,i)=Mdi(j,1:7)*X0; end; end; Xlsd(1:2,1)=eye(2)*X01; for j =1:1:2; for i=2:1:101; Asdi=Asd^i; Xlsd(j,i)=Asdi(j,1:2)*X01; end;

90

end; Xrsd(1:5,1)=eye(5)*X0r; for j =1:1:5; for i=2:1:101; Dsdi=Dsd^i; Xrsd(j,i)=Dsdi(j,1:5)*X0r; end; end; %Partitionarea matricii M, a sistemului liniar: A=M(1:2,1:2); B=M(1:2,3:7); C=M(3:7,1:2) ; D=M(3:7,3:7); %Aplicarea Metodei PERTURBATIEI SINGULARE % Iteratia in cazul partii lente: L1=inv(D)*C; A1=A-B*L1; C1=C-D*L1+L1*A1; D1=D+L1*B; %Iteratia in cazul partii rapide H11=B*inv(D1); A11=A1-H11*C1; B11=B-H11*D1+A11*H11; D11=D1+C1*H11; %Discretizarea noilor matrici ale modelului: A11d=expm(A11*h); D11d=expm(D11*h); X11(1:2,1)=eye(2)*X01; Z11(1:5,1)=eye(5)*X0r; for i=2:1:101 A11di=A11d^i; D11di=D11d^i; zita1=A11di*X01; ita1=D11di*X0r; X11(1:2,i)=zita1+H11*ita1; Z11(1:5,i)=ita1-L1*X11(1:2,i); end; axis([0 100 -0.1 0.15]) i=0:1:100 ; plot(i,Xd(1,i+1),i,Xlsd(1,i+1),i,X11(1,i+1));grid; xlabel ( 'timp/0.08[sec]' ) ; ylabel('[pu]'); title('Variatia lui e_q'); h=figure; axis([0 100 -0.02 0.06])

91

i=0:1:100; plot(i,Xd(2,i+1),i,Xlsd(2,i+1),i,X11(2,i+1));grid; xlabel('timp/O.08 [sec]'); ylabel('[pu]'); title('Variatia lui R__f '); h=figure; axis([0 100 -0.15 0.05]) i=0:1:100; plot(i,Xd(3,i+1),i,Xrsd(1,i+1),i,Z11(1,i+1));grid; xlabel('timp/0.08 [sec]'); ylabel('[pu]'); title('Variatia lui e_d'); h= figure; axis([0 100 -0.1 0.15]) i=0:1:100; plot(i,Xd(4,i+1),i,Xrsd(2,i+1),i,Z11(2,i+1));grid; xlabel('timp/0.08 [sec]'); ylabel('[pu]'); title('Variatia lui delta'); h=figure; axis([0 100 -0.05 0.4]) i=0:1:100; plot(i,Xd(5,i+1),i,Xrsd(3,i+1),i,Z11(3,i+1));grid; xlabel('timp/0.08 [sec]'); ylabel('[pu]'); title('Variatia lui omega'); h=figure; axis([0 100 -0.15 0.3]) i=0:1:100; plot(i,Xd(6,i+1),i,Xrsd(4,i+1),i,Z11(4,i+1));grid; xlabel('timp/0.08 [sec]'); ylabel('[pu]'); title ('Variatia lui E__fd') ; h=figure; axis([0 100 -0.6 2]) i=0:1:100; plot(i,Xd(7,i+1),i,Xrsd(5,i+1),i,Z11(5,i+1));grid; xlabel('timp/0.08 [sec]'); ylabel('[pu]'); title('Variatia lui V_R');

SIMULAREA DINAMIC A UNUI GRUP ENERGETIC Matricea sistemului:

92

M=

-0.5800

0 -0.2690 0 0 2.1200 0 314.0000 0

0.2000 0 0 0 0 2.0000

0 -1.0000 0 0 -0.1410 0 0

1.0000 0 0 0

0 -5.0000 0 0 0

0.1410 -0.2000 -0.2800 0 0 0 0.8380

-173.0000 66.7000 -116.0000 40.9000

0 -66.7000 -16.7000

Valorile proprii ale matricii M:

Lambda_m =

-8.1992 + 7.8150i -8.1992 - 7.8150i -0.8348 + 7.6004i -0.8348 - 7.6004i -0.3332 + 0.5847i -0.3332 - 0.5847i -3.9875

93

7. VALIDAREA METODEI PRIN SIMULARE NUMERIC Variaia lui

Variaia lui Rf

94

Variaia lui

Variaia lui delta

Variaia lui omega

Variaia lui Efd

96

Variaia lui VR

97

8. CONCLUZII

n lucrarea de fa se prezint metoda perturbaiei singulare, al crui avantaj este c elimin inconsistenele aproximrii clasice de stare cvasistaionar, modelul redus obinut fiind de ordin inferior. Se prezint o metod iterativ de separare a scalelor de timp, cu ajutorul creia se mbuntete sistematic precizia modelelor reduse. Se arat, prin exemplul concludent al modelului unui transformator monofazat, c n cazul n care variabilele alese de stare nu sunt separabile, atunci printr-o redefiniere adecvat a mrimilor de stare se poate obine un model cu variabile separabile. Modelarea i separarea scalelor de timp este ilustrat n cazul unui grup energetic, rezultatele simulrii post-defectului aprut la acesta fiind foarte suggestive. Astfel se constat c modelul classic de stare cvasistaionar introduce erori foarte mari n liniarizarea mrimilor lente i , n timp ce dup o singur iteraie a subsistemelor

rapid i lent ( j=k=1 ), eroarea este mult mai mica. Rspunsul mrimilor rapide pentru ambele modele rapde: cel cvasistaionar i cel corectat dup o iteraie, se apropie foarte mult de rspunsul pentru modelul exact, dar cel corectat l aproximeaz mai bine. Avantajele metodei perturbaiei singulare se evideniaz mai pregnant n cazul modelului
98

lent dect n cazul modelului rapid, dar trebuie subliniat faptul c i n acest caz eroarea raspunsului pentru modelul corectat, conform metodei perturbaiei singulare, fa de rspunsul exact este mai mic dect pentru rspunsul modelului necorectat cvasistaionar.

BIBLIOGRAFIE

1. Clin S., Dumitrache I. Regulatoare automate, Ed. Didactic i Pedagogic, Bucureti 1985; 2. Corlateanu V. Curs de maini electrice. Vol 2, Maini sincrone, Ed. Institutului Politehnic Iai; 3. Gheorghiu I., Fransua A. Tratatat de maini electrice, Ed. Academiei, Bucureti 1972; 4. Ghinea M., Fireteanu V. Matlab. Calcul Numeric grafic aplicaii , Ed. Teora, Bucureti 1995: 5. Kokotovic P.V. Singular perturbation and iterative separation of time scale; 6. Kokotovic P.V. Application of singular perturbation techniques to control problems; 7. Lazu C. Maini electrice, Ed. Didactic i Pedagogic, Bucureti 1966
99

8. Radulet R., Opaschi M. Proiectarea hidrogeneratoarelor i a motoarelor sincrone, Ed Tehnic, Bucureti 1980

100