PRORAČUN NAPONA KRATKOG SPOJA TRANSFORMATORA METODOM KONAČNIH ELEMENATA FINITE ELEMENT CALCULATION OF SHORT-CIRCUIT VOLTAGE

Transcription

1 Leonardo Štrac, dipl. ing. Končar energetski transformatori d.o.o. PRORAČUN NAPONA KRATKOG SPOJA TRANSFORMATORA METODOM KONAČNIH ELEMENATA SAŽETAK Sa sveprisutnom težnjom za skraćivanjem vremena projektiranja i izrade transformatora i smanjivanja njegovih dimenzija, sve je veća potreba za modeliranjem i simuliranjem transformatora pomoću računala. U ovom članku opisana je metodologija proračuna napona kratkog spoja energetskih transformatora te metoda konačnih elemenata. Za proračun je korišten programski paket Ansoft Maxwell. Uspoređen je D i 3D model s rezultatima mjerenja za transformator 00MVA. Na osnovu usporedbe dana je preporuka za modeliranje. Ključne riječi: energetski transformator, napon kratkog spoja, metoda konačnih elemenata FINITE ELEMENT CALCULATION OF SHORT-CIRCUIT VOLTAGE ABSTRACT In production of power transformers, there is an increased tendency of cutting down projecting and manufacturing time as well as shrinking transformers dimensions. This process is increasing the need for computer modeling and simulation of transformers. Article describes methods for short-circuit voltage calculation and the finite element method. Calculations and finite method analysis are made with software package Ansoft Maxwell. D and 3D models of 00MVA transformer are compared with the experimental results and based on those conclusions it is given a recommendation for modeling. Key words: power transformer, short-circuit voltage, finite element method 1. UVOD Na tržištu energetskih transformatora konkurencija je sve veća. Zahtjev za smanjenjem vremena izrade transformatora stalno je prisutan. Težnja za profitom, konkurencija i standardi stalno smanjuju dimenzije transformatora, no pri tome ostaje zahtjev za kvalitetom i ugovorenim tehničkim karakteristikama transformatora. S obzirom na karakter proizvodnje energetskih transformatora komadna proizvodnja po narudžbi i cijenu proizvoda red veličine , nemoguće je napraviti prototip koji bi služio za ispitivanje, provjeru projektantskih jednadžbi i rješenja. Iz tog je razloga nužno modeliranje transformatora pomoću računala i simuliranje pojava. Napon kratkog spoja važna je tehnička karakteristika transformatora. U pravilu, kupac ga točno definira s obzirom na karakteristike svoje energetske mreže ili pak transformatora koji će raditi u 1

2 paralelnom radu. Ugovorom su često propisane tolerancije i novčani penali u slučaju odstupanja od ugovorenih vrijednosti. Sve je to samo razlog više za provjeru projektnih rješenja na računalnom modelu. Končar energetski transformatori od 001. godine posjeduje simulacijski paket Ansoft Maxwell kojim se mogu analizirati magnetostatski, elektrostatski, tranzijentni i termički problemi, utjecaj vrtložnih struja, izmjeničnog elektromagnetskog polja i dr. u -D i 3-D modelima.. NAPON KRATKOG SPOJA Napon kratkog spoja je onaj napon koji će uz kratko spojeni sekundarni namot kroz primarni namot potjerati nazivnu struju. Najčešće se izražava kao postotak nazivnog napona: I Z S Z u = 100 = 100, (1) n k n k k% gdje su Z k impedancija kratkog spoja, a S n, I n, i U n nazivne vrijednosti snage, struje i napona. Nadalje, djelatna i jalova komponenta napona kratkog spoja računaju se prema jednadžbama: u 100 I R n k r% = () u 100 I X n k σ% =, (3) gdje su R k i X k radna i jalova komponenta impedancije kratkog spoja, I n nazivna struja. U pokusu kratkog spoja u kojem se u k% određuje, narinuti napon je puno manji od nazivnog, pa su i gubici u jezgri puno manji, te se kod ovog pokusa zanemaruju. Amperzavoji primarnog i sekundarnog namota su približno jednaki i suprotni te je gotovo sav tok u transformatoru rasipni. Napon kratkog spoja i rasipni induktivitet vezani su jednadžbom: L = σ u u 100 ω S n k% r%, (4) snaga. gdje su U k% i U r% napon kratkog spoja i njegova radna komponenta, a U n i S n nazivni napon i Rasipni induktivitet može se izračunati iz energije rasipnog magnetskog toka W σ što je osobito pogodno za računanje metodom konačnih elemenata: 1 W L I W L = σ σ (6) I σ = σ (5) 1 Wσ = Bσ Hσ dv, (7) V gdje je L σ rasipni induktivitet, B σ i H σ indukcija i jakost rasipnog magnetskog polja, a I struja namota.

3 3. METODA KONAČNIH ELEMENATA 3.1. Općenito Rješavanje elektromagnetskih problema metodom konačnih elemenata uključuje nekoliko osnovnih koraka. Prvi korak je crtanje modela koji dimenzijama odgovara fizičkom modelu. Model može biti prikazan u D kartezijevom sustavu, u rotacijsko simetričnom odnosno cilindričnom sustavu ili u 3D kartezijevom sustavu. Drugi korak je definiranje materijala za svaki dio modela. Materijali se definiraju isključivo kroz svoje električne i magnetske karakteristike (vodljivost, permeabilnost i sl.) jer njihova mehanička i druga svojstva nemaju utjecaja na raspodjelu elektromagnetskog polja u prostoru. Sljedeći korak je definiranje izvora električnih i magnetskih polja: električni potencijali, naboj, struje ili pak eksplicitno navedeni izvori polja. Sada je moguće pristupiti rješavanju problema. Da bi se prikazalo električno ili magnetsko polje u raznim područjima modela, model je podijeljen na velik broj manjih elemenata trokuta u D modelu, odnosno tetraedara u 3D modelu. Polje u svakom elementu predstavljeno je zasebnim polinomom. Skup svih trokuta ili tetraedara predstavlja mrežu konačnih elemenata. 3.. Ansoft Maxwell Ansoft Maxwell je program za rješavanje elektromagnetskog polja metodom konačnih elemenata. U programu je moguće modelirati D i 3D modele te rješavati probleme s vremenski promjenjivim magnetskim poljem te elektrostatske i magnetostatske probleme. Maxwell koristi iterativnu metodu stvaranja mreže usitnjavajući mrežu na mjestima s najvećom greškom. Tako se izbjegava zagušenje kapaciteta računala, a opet osigurava ciljana točnost rješenja. Točnost izračunatog polja ovisi o veličini pojedinih elemenata, a točnost cjelokupnog rješenja ovisi o točnosti polja u svakom pojedinom elementu. Što je veličina elementa manja, to je polje ujednačenije unutar njegovih granica. Međutim, generiranja rješenja polja uključuje invertiranje matrice sa približno onoliko elemenata koliko ima vrhova tetraedara (trokuta) što za mreže s velikim brojem elemenata zahtijeva znatnu količinu memorije i procesne snage Magnetostatski problem Rješavanje magnetostatskog problema odvija se u dva koraka. U prvom se računa gustoća struje J u vodičima prema zadanim izvorima istosmjerne struje. Ako je gustoća struje J eksplicitno zadana, ovaj korak se preskače. Na objektima koji su definirani kao savršeni vodiči računaju se površinske struje jer polje u takve objekte ne prodire. S obzirom da u statičnom polju nema gibanja naboja, vrijedi jednadžba (8). Kada se u (8) uvrsti (9), dobije se jednadžba (10) koja se rješava u prvom koraku rješavanja magnetostatskog problema. ρ div J = = 0 (8) t J= σ E = σ grad ϕ (9) div ( σ grad ϕ ) = 0, (10) gdje je J gustoća struje, ρ gustoća naboja, σ specifična vodljivost, a ϕ električni potencijal. U drugom koraku računa se statično magnetsko polje H s gustoćom struje kao izvorom polja prema jednadžbama: 3

4 rot H = J (11) div B = 0 (1) B = µµ H, (13) 0 r gdje je H jakost magnetskog polja, J gustoća struje, B magnetska indukcija, a µ 0 i µ r permeabilnost vakuuma i relativna permeabilnost. Energija magnetskog polja računa se u linearnim materijalima prema (14). W L 1 = B H dv (14) V gdje je H jakost magnetskog polja, a B magnetska indukcija. 4. D MODEL TRANSFORMATORA ZA PRORAČUN NAPONA KRATKOG SPOJA Modelirani autotransformator snage je 00MVA i ima pet namota: visokonaponski (VN) niskonaponski (SN), tercijar koji nije opteretiv već služi za stabilizaciju (NN) te dva regulacijska (FRA i FRB) čiji se zavoji dodaju ili oduzimaju visokonaponskom namotu radi promjene napona visokonaponske mreže. Simulirano je šest radnih položaja: protjecane su kombinacije namota SN-VN, NN-VN, NN-VN±FRA/B i SN-VN±FRA/B. Zavoji regulacijskih namota dodaju se zavojima visokonaponskog kada je visokonaponska mreža na višem naponu od nazivnog plus položaj, a oduzimaju kada je visokonaponska mreža na naponu nižem od nazivnog minus položaj regulacije. Rješenja elektromagnetskih polja D modela dobivena su u osam koraka adaptivnog usitnjavanja mreže i imaju približno trokuta. linija 5 linija 4 linija 3 z linija r linija 1 FRB FRA NN SN VN JEZGRA Slika 1: crtež D modela u cilindričnom koordinatnom sustavu s ucrtanim pomoćnim linijama 1 5 4

5 Tablica I: protjecanja namota za 00MVA (protjecanje NN za 4MVA) efektivno protjecanje [Az] kombinacija NN SN VN FRA FRB VN(-)/NN VN(+)/NN VN(0)/NN VN(-)/SN VN(+)/SN VN(0)/SN Nakon dobivanja konvergentnih rješenja, za svaku kombinaciju protjecanja izračunata je ukupna energija u prostoru problema prema jednadžbi 15. Za rješavanje takvih problema koristi se Maxwellov field calculator u postprocesoru. W uk = Wdrdz, (15) r z gdje je W energija u prostoru. Kad je poznata ukupna energija u prostoru, napon kratkog spoja lako se izračuna prema jednadžbi: ω Wuk π f Wuk 00 π Wuk uk% = = =, (16) I U I U S f f f f f gdje su S f, U f i I f fazna snaga, napon i struja, W uk ukupna energija u prostoru, f frekvencija, a ω kružna frekvencija. Tablica II: snaga, energija u prostoru i napon kratkog spoja za različite kombinacije protjecanja S [MVA] W uk [J] u k [%] VN(-)/NN 4 375,89,95 VN(+)/NN 4 468,68 3,68 VN(0)/NN 4 410,6 3, VN(-)/SN ,75 10,34 VN(+)/SN ,1 15,65 VN(0)/SN ,47 1,16 5

6 Slika : z komponenta indukcije B z [T] za kombinacija protjecanja SN-VN Slika 3: magnetski tok Φ[Wb] za kombinacija protjecanja SN-VN 6

7 5. 3D MODEL TRANSFORMATORA ZA PRORAČUN NAPONA KRATKOG SPOJA 3D model transformatora razlikuje se od D rotacijsko-simetričnog modela u izgledu jarma, dok su namoti i stup jezgre jednaki. Model je prepolovljen tri puta po tri ravnine simetrije tako da je proračun napravljen s osminom cijelog modela. Crtež nacrta, bokocrta i nebocrta prikazana je na slici 8. Protjecanja namota dana su u tablici I. Rješenja elektromagnetskog polja 3D modela dobivena su u šesnaest koraka adaptivnog usitnjavanja mreže, a imaju približno tetraedara. Slika 4: crtež 3D modela s ucrtanim pomoćnim linijama Za svaku kombinaciju protjecanja izračunata je u postprocesoru ukupna energija prostora (jednadžba 17) te je napon kratkog spoja izračunat prema jednadžbi 17. Izračunate vrijednosti prikazane su u tablici III. W uk = Wdxdydz, (17) dx dy dz gdje je W energija u prostoru. Tablica III: snaga, energija u prostoru i napon kratkog spoja za različite kombinacije protjecanja S [MVA] W uk [J] u k [%] VN(-)/NN 4 366,00,87 VN(+)/NN 4 449,4 3,53 VN(0)/NN 4 397,16 3,1 VN(-)/SN ,13 10,01 VN(+)/SN ,80 14,88 VN(0)/SN ,90 11,7 7

8 Slika 5: z komponenta indukcije B[T] za kombinaciju protjecanja SN-VN0 Slika 6: indukcija B[T] za kombinaciju protjecanja SN-VN0 8

9 6. USPOREDBA REZULTATA D I 3D MODELA S MJERENJIMA Tablica IV pokazuje da 3D model daje točnije rezultate nego D. Srednje odstupanje simuliranih od mjerenih vrijednosti kod D modela je: p =4,33%, dok je kod 3D modela p =1,43%. Ti podaci govore da je 3D proračun pouzdaniji i ima manje rasipanje rezultata. Još je važnije pogledati apsolutna odstupanja koja kod D proračuna mogu prelaziti i 1%. Tablica IV: usporedba rezultata simulacije D i 3D modela s mjerenjima S [MVA] uk [%] - D uk [%] - 3D uk [%] - mjereno (mj-d)/mj [%] (mj-3d)/mj [%] VN(-)/NN 4,95,87,9-1,03 1,71 VN(+)/NN 4 3,68 3,53 3,54-3,95 0,8 VN(0)/NN 4 3, 3,1 3,14 -,55 0,64 VN(-)/SN 00 10,34 10,01 9,70-6,60-3,0 VN(+)/SN 00 15,65 14,88 14,53-7,71 -,41 VN(0)/SN 00 1,16 11,7 11,68-4,11-0,34 7. ZAKLJUČAK Usporedbom rezultata simulacije i mjerenja pokazuju da proračun metodom konačnih elemenata može dati dovoljno točne rezultate za potrebe projektiranja transformatora. Za primjene s manjom zahtijevanom točnošću svakako je bolje odlučiti se za D model. D model je bitno jednostavniji i brži u svim koracima izrade projekta, a rezultati se nalaze u okviru ±10% odstupanja od mjerenih vrijednosti. Pored toga, zahtjev za količinom memorije i procesnom snagom bitno je manji u odnosu na 3D modeliranje. Kada je potrebna veća točnost proračuna, svakako treba odabrati 3D model koji puno bolje opisuje fizikalni izgled modeliranog objekta, te je stoga i odstupanje od mjerenih vrijednosti gotovo tri puta manje nego kod D modela. Srednje postotno odstupanje p =1,43% i najveće apsolutno odstupanje 0,35% prihvatljive su vrijednosti, te se 3D model može sa sigurnošću prihvatiti kao projektantski alat. 8. LITERATURA [1] T. Kelemen: Transformator, Tehnička enciklopedija svezak 13, Zagreb [] Ansoft Maxwell 9: Online help, June 00. [3] Z. Haznadar, Ž. Štih: Elektromagnetizam, Zagreb [4] H. Wang: Finite Element Analysis of Internal Winding Faults in Distribution Transformers, IEEE transactions on power delivery, vol. 16, no. 3, July