+ All Categories
Home > Documents > Transformatorul,Sinteza1,optimizare,2011

Transformatorul,Sinteza1,optimizare,2011

Date post: 10-Jul-2015
Category:
Upload: ana-nistor
View: 74 times
Download: 0 times
Share this document with a friend
18
 Transformatoare trifazate cu conversie radiala (infasurari cilindrice concentrice) Partea II Sinteza(Proiectarea) 1 Breviar 1.1 Solutia constructiva -miez magnetic trifazat,nesimetric (3 coloane-c,2juguri-j),dispuse in acelasi plan -infasurari cilindrice,dispuse concentric (fig.1) JT JT JT JT JT JT IT IT IT IT IT IT a1 a1 a2 a2 Lj H F Hb D1 2 D a12 a22 ho2 Fig.1 In cazul infasurarilor concentrice,(fig.1),infasurarea de joasa tensiune(JT,avand w 1 spire,parcurse de curentul de faza I f1 )este dispusa interior,langa miez,iar infasurarea de inalta tensiune(IT),cu w 2 spire,parcurse de curentul de faza I f2, la exterior. Principalele dimensiuni indicate in fig 1 sunt urmatoarele: -diametrul mediu (D 1,2 ) si inaltimea medie a celor doua infasurari(H  b ) -grosimile radiale ale infasurarilor:a 1 (J.T.,in interior,langa miez) si a 2 (I.T,in exterior) -distantele de izolatie radiale: a01 (miez –JT),a12 (JT-IT), a22 (IT-IT) si axiale:h10 (JT-miez),h20 (miez-IT) -dimensiunile ferestrei: H (inaltimea),F (largimea) -dimensiunile coloanei: D -diametrul circumscris treptelor, -dimensiunile jugului: L  j -lungimea jugului 1.2 Vectorul Poynting si dimensiunile principale ale transformatoarelor Transformatorul electric este un convertor electromagnetic static,in care se realizeaza o conversie de tip electric-electric(E-E). Puterea electromagnetica este preluata de la sursa de catre infasurarea primara si transmisa infasurarii secundare, prin spatiul (cilindric)dintre cele doua infasurari- sediul campului magnetic de scapari, cu  suprafata A  P -aria Poynting,(m 2 )- cu densitatea de suprafata a puterii electromagnetice σ   H  E S  s  p × = (VA/m 2 ) –vectorul Poynting,orientat radial (fig.2). (1)  Dimensiunile principale ale transformatorului cu conversie radiala sunt dimensiunile cilindrului de transfer: diametrul 1 2  D si inaltimea b  H ,respectiv diametrul si inaltimea medie a infasurarilor, iar  B  H  D / 12 π β = -o marime adimensionala,  factorul de geometrie (arhitectura) a transformatorului.
Transcript
Page 1: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 1/18

Transformatoare trifazate cu conversie radiala(infasurari cilindrice concentrice)

Partea II Sinteza(Proiectarea)1 Breviar

1.1 Solutia constructiva-miez magnetic trifazat,nesimetric (3 coloane-c,2juguri-j),dispuse in acelasi plan-infasurari cilindrice,dispuse concentric (fig.1)

JT JT JT JT JT JTIT IT IT IT IT IT

a1 a1a2 a2

Lj

H

F

Hb

D12

D

a12 a22 ho2

Fig.1

In cazul infasurarilor concentrice,(fig.1),infasurarea de joasa tensiune(JT,avand w1 spire,parcurse decurentul de faza If1)este dispusa interior,langa miez,iar infasurarea de inalta tensiune(IT),cu w2

spire,parcurse de curentul de faza If2, la exterior.

Principalele dimensiuni indicate in fig 1 sunt urmatoarele:-diametrul mediu (D1,2) si inaltimea medie a celor doua infasurari(H b)-grosimile radiale ale infasurarilor:a1 (J.T.,in interior,langa miez) si a2 (I.T,in exterior)-distantele de izolatie radiale: a01 (miez –JT),a12 (JT-IT), a22 (IT-IT) si axiale:h10 (JT-miez),h20 (miez-IT)-dimensiunile ferestrei: H (inaltimea),F (largimea)-dimensiunile coloanei: D -diametrul circumscris treptelor,-dimensiunile jugului: L j-lungimea jugului

1.2 Vectorul Poynting si dimensiunile principale ale transformatoarelor

Transformatorul electric este un convertor electromagnetic static,in care se realizeaza o conversie de tip

electric-electric(E-E). Puterea electromagnetica este preluata de la sursa de catre infasurarea primara sitransmisa infasurarii secundare, prin spatiul (cilindric)dintre cele doua infasurari- sediul campuluimagnetic de scapari, cu suprafata A P -aria Poynting,(m2)- cu densitatea de suprafata a puterii

electromagnetice σ H E S s p ×= (VA/m2) –vectorul Poynting,orientat radial (fig.2). (1)

Dimensiunile principale ale transformatorului cu conversie radiala sunt dimensiunile cilindrului detransfer: diametrul 12 D si inaltimea b H ,respectiv diametrul si inaltimea medie a infasurarilor, iar

B H D /

12π β = -o marime adimensionala, factorul de geometrie (arhitectura) a transformatorului.

Page 2: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 2/18

w1

Sw 1 Sw 2

w2

i1 i2

SpEs

Hs

Infasurare

Primara(receptor )

Infasurare

Secundara(generator )

a1 a2

Hb

a01

D

ApHs

(A/m)

a1 a2

Hs=w1*I1/Hb=w2*I2/Hb

a12

D12

a12

Fig.2

În relaţia (1), Es (V/m) reprezintă intensitatea câmpului electric indus prin transformare (un vector orientat in lungul spirei in sensul curentului la infasurarea generatoare si in sens invers la infasurarea

receptoare), Hσ (A/m) intensitatea câmpului magnetic de scăpări (un vector axial pentru configuratiadin fig.2.a, cu distributia trapezoidala din fig.2b), iar B P H D A ⋅=12

π (m2) -aria laterală a cilindrului descăpări pe care o vom defini ca suprafata,de transfer, Poynting.Vectorul Poynting P

S rezulta astfel radial, fiind indreptat spre infasurarea secundara, regim degenerator a transformatorului,transformatoarele cu infasurari concentrice putand fi denumite catransformatoare cu conversie de tip radial.

În regim sinusoidal , se defineşte puterea electromagnetică aparentă SΣca semiprodusul valorilor devârf s E ˆ şi σ

H ˆ respectiv ca produsul valorilor efective ale acestora (extinzând formal aceasta noţiuneşi în ceea ce priveşte câmpul magnetic ), astfel încât se poate scrie urmatoarea relatie putere-dimensiuni:

P P BS e AS H D H E S S ⋅=⋅⋅⋅=≅Σ 12π

σ (VA) (2)

Intensitatea câmpului electric indus Es(V/m)), se obţine aplicând legea inducţiei electromagnetice uneispire (medii) de diametru D12, traversată de câmpul magnetic principal, respectiv din tensiunea de spira

12/ DU E SpS

π ≅ (V/m), unde 2/ Fe Fe sp B AU ⋅⋅=ω (V/sp) (3)unde BFe reprezintă valoarea, de vârf, a inducţiei magnetice în miezul magnetic (coloana), iar AFe ariasectiunii de fier corespunzatoare.Intesitatea câmpului magnetic de scăpări 2ˆ

σ σ H H = (sau S H -fig.2b)se determină din legea

circuitului magnetic

1221

11 Ai A A H

I w H

B

′=≈=⋅

≅σ (A/m), unde 21 A A

, reprezintă pânza de curent (cu 12 ii≈′

).Înlocuind aceste rezultate în relaţia (2) se obţine:

C s p B

B

s p

BS e S I U I U w H D H

I w

D

U H D H E S ≈⋅=⋅⋅=⋅⋅

⋅=⋅⋅⋅=⋅ 11111 2

11

1 2

1 2π

π

π σ

(VA) (4)

respectiv puterea aparentă (pe o coloană),transferata intre cele doua infasurari ale transformatorului.

Page 3: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 3/18

Impedanta Poynting P Z se defineste ca raportul dintre S E si σ H obtinandu-se expresia :

2

1

1

w

Z

H

E Z nS

P ⋅

==

β σ

(Ω) unde f n

f n

n I

U Z

,1

,1

1 = (Ω) si B H

D12π β = (5)

1.3 Relaţia pierderilor Joule (PJ ) si a fluxului termic (Qt )

Daca se considera infasurarea (concentrica) ca o panza de curent de grosime ( ,in m), inaltime HB sidiametru mediu D12, si desitatea (medie,pentru cele doua infasurari) de curent J ,sunt valabile relatiile; )/()/()()/(

2 m A Am A H mm A J ==∆⋅σ

, (6)

unde B H I w A /⋅= , in A/m, este densitatea lineara de curent a panzei de curent .

Pierderile Joule din infasurare,respectiv din panza de curent,se pot scrie astfel sub forma :)()( 2

1212

22 A J m A H D A J H D J V J P P B BCu J ⋅⋅⋅=⋅⋅⋅⋅=∆⋅⋅⋅⋅=⋅⋅= ρ π ρ π ρ ρ (W) (7)

Fluxul termic unitar , pe suprafata de transfer(Poynting) a puterii electromagnetice(AP ),in acelasi timpsi suprafata de referinta pentru transmisia caldurii prin convectie, se poate astfel determina cuurmatoarea relatie:

A J H D P Q B J t ⋅⋅=⋅= ρ π )/(12 (W/m2) (8)

respectiv, fluxul termic unitar este proportional cu sarcina termica a transformatorului ( A J ⋅ ).

Din expresia de mai sus, se poate obtine relatia de calcul a pierderilor Joule(nominale) , in u.r. :

r E

E

E

J

H E

A J

S

Q

H DS

H DQ

S

P p

S

J

S S P

t

B P

Bt

e

J J ==

⋅=

⋅⋅==

⋅⋅

⋅⋅==

ρ ρ

π

π

σ 12

12 (u.r) (9)

astfel încât pierderile Joule nominale, exprimate în u.r.,respectiv rezistenta infasuraii,in u.r.,suntnumeric egale cu raportul dintre intensităţile câmpului electrocinetic (EJ) respectiv indus (Es).Relatia (9) este fundamentala, stabilind cea mai simpla si concisa exprimare in u.r. a pierderilor Joulenominale–respectiv a rezistentei infasurarii (in unitati relative r =pJ ), ca raportul dintre intensitateacampului electric de conductie ( J E J ⋅= ρ ) si intensitatea campului electric indus( S E ).Pe de alta parte, avand in vedere relatia pierderilor de scurtcircuit ( J J J kn p p p p 221 ≅+= ) relatia (9) semai poate scrie sub forma :

kn

t

J

t

S S P p

Q

p

Q A E H E S

21

≅=== σ (VA/m2) (10)

respectiv, vectorul Poynting este numeric egal cu raportul dintre fluxul termic unitar (Q t, W/m2) şi pierderile Joule (pJ1, u.r.) din înfăşurare, stabilindu-se astfel o relaţie directă între solicitarea termică atransformatorului (proporţională cu fluxul termic unitar Qt), pierderile de scurtcircuit (pkn) şi vectorul

Poynting (S p).

1.4 Relaţia tensiunii de scurtcircuit (uk )

Tensiunea de scurtcircuit nominala (in u.r.),se determina in functie de componentele activa si reactiva:2

,

2

,, r k ak nk uuu += (u.r.),unde componenta activa nk k ak pr r r u ,

,

21, =+== (u.r.)

(11)

Page 4: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 4/18

Componenta reactiva a tensiunii (nominala) de scurtcircuit ( r k u

, )este, (in u.r.),numeric egala cu

reactanta de scurtcircuit,respectiv cu reactanta totala de scapari: ,

2,1,,1, σ σ σ x x x xuk r k

+=== (u.r),unde:

nk k Z L x ,1,1,1 /⋅=ω (u.r.), cu σ σ σ Λ⋅=+=

2

1

,

,2,1,1 w L L L k (12)Permeanta totala a canalului de scapari se calculeaza cu relatia generala Bo H A /

σ σ µ ⋅≅Λ (H),

unde σ Α (m2)-aria echivalenta strabatuta de campul magnetic de scapari, este o suprafata inelara cu

diametrul mediu 2,1 D si grosimea echivalenta 3/)( 212,1 aaaa ++=σ (fig.2), σ σ π a D A ⋅= 12 ,

respectiv:

S

o

P

o

n

o

B

o

n

r k E

H a

Z

a

Z

aw

H

Da

Z

w xu σ

σ σ σ σ

σ µ ω µ ω β µ ω π µ ω

⋅⋅⋅=

⋅⋅

=

⋅⋅⋅⋅

=

⋅⋅

==

,1

2

112

,1

2

1, (u.r.) (13)

1.5. Relaţia factorului (mediu) de umplere a înfăşurărilor (k u )

Se definesc factorii de umplere(subunitari) ai celor doua infasurari(joasa tensiune,inalta tensiune) 1,uk

si 2,uk ca raportul dintre aria neta de material conductor(cupru sau aluminiu),ocupata,in sectiune, deconductoarele celor doua infasurari( 11 wS w ⋅ ,si 22 wS w ⋅ ),si suprafata de gabarit a acestora, B H a ×1 ,si

B H a ×2 (fig.1),respectiv:

B

W

u

B

W

u H a

S wk si

H a

S wk

2

22

2

1

11

1 == (<1,si 2,uk < 1,uk ,avand in vedere ponderea izolatiei in inf.I.T ) (14)

In functie de factorii de umplere individuali, 1,uk si 2,uk se defineste factorul mediu de umplere al

celor doua infasurari:

B

W W u

H aa

S wS wk

)( 21

2211

+

+= =

)(

2

)(

2

)(

22

211211

11

21

1,1

21

21

aa J

H

H aa J

I w

H aa

S w

k k

k k

B B

w

uu

uu

+⋅

⋅≅

⋅+⋅

⋅⋅=

⋅+

⋅≅

+

σ

unde s-a avut in vedere ca,pentru o aceeasi densitate de curent, J J J =≅ 21 ,avem 11 wS w ⋅ ≅ 22 wS w ⋅

Din relatia de definitie (14) se poate explicita suma ( 21aa + ):

J k

H aa

u ⋅

⋅=+

σ 221 =

σ µ ω

ρ ρ

a pk

u

Z pk knu

kr

P knu ⋅⋅⋅⋅

⋅⋅

=

⋅⋅

0

44(m) , cu

3

2112

aaaa

++=σ (15)

unde s-a inlocuit ρ ⋅

⋅=

2

S kn E p J (conform relatiei 9), si

σ µ ω a

u Z kr

P ⋅⋅

=

0

.(Ω),

Din ecuatia (15) rezulta valoarea largimii echivalente a canalului de scapari( σ a )

22

1212(

2

1 Raaa ++⋅=

σ (m), undeukn

r k

k p

u R

⋅⋅⋅

⋅⋅=

ω µ

ρ

0

,2

3

16(m2) (16)

Rezulta astfel expresii de calcul a dimensiunilor radiale ale infasurarilor:

21

2

211)(

uu

u

k k

k aaa

+⋅+= , respectiv

21

1

212 )(uu

u

k k

k aaa

+⋅+= (m) ,unde )(3 1221 aaaa −=+ σ (17)

relatii analitice importante pentru grosimea infasurarilor (in mod uzual estimate cu o relatieempirica),evaluata astfel in concordanta cu principalele date de proiectare ale

Page 5: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 5/18

transformatorului:materialul infasurarilor si temperatura de functionare( ρ ),frecventa( f ),parametride scurtcircuit( nk

p, , nk u

, ),si factorul mediu de umplere a infasurarilor( 1,uk , 2,uk ).

1.6 Relatia tensiunii de spira

Pornind de la expresia de calcul a componentei reactive a tensiunii de scurtcircuit,(13),

n

k kr Z

aw xu

1

2

10 β ω µ σ ⋅⋅⋅== ,unde c f n f n f nn S U I U Z // 2

,1,1,1,1 == , si sp f n U wU ⋅≅ 1,1 (18)

se obtine o relatie alternativa pentru tensiunea de spira,in functie de puterea pe coloana (S c) si factorulde geometrie :

r k oc sp uaS U ,. /)( β µ ω σ ⋅⋅⋅⋅= (V) , unde Bm H D /π β = (19)

respectiv sectiunea de fier a miezului(coloanei):2,

22

Fer k

oc

Fe

sp

Fe

Bu

aS

B

U A

⋅⋅

⋅⋅⋅⋅=

⋅=

ω

β µ

ω

σ (m2) (20)

2. Algoritmi de dimensionare2.1 Algoritmul tensiunii de spiraRelatia de baza este relatia tensiunii de spira (19),care permite generarea analitica a secventelor decalcul pentru dimensiuni geometrice,solicitari electromagnetice,consumuri de materiale,etc.,in functiede factorul de geometrie ,si principalele date de proiectare(putere,frecventa,parametri de scc.,etc).Astfel, pentru factorul de geometrie variabil ( 5.....5.0∈β ),se calculeaza urmatoarele secvente:

a) Tensiunea de spira )/)( ,. r k oc sp uaS U β µ ω σ ⋅⋅⋅⋅= , unde ),,,,(,,

ω ρ σ ur k nk

k u p f a = , (21)

b) Sectiunea de fier (utila) a miezului magnetic: )/(2 Fe FeBUsp A ⋅⋅= ω ,(se adopta . Fe B )

(22)Observatie. Valoarea inductiei magnetice in miez,BFe, se poate evalua in functie de performantele

dorite la mersul in gol, respectiv po,n si io,n (in u.r.),avand in vedere urmatoarele relatii:)()/(.).( Fe Fe Feo B pVAkg cr u p

⋅≅ si )()/(.).( Fe Fe Feo BqVAkg cr uq∗

⋅≅ ,respectiv ∗∗≅ Fe Feoo q pq p //

(23)unde ∗

Fe p (w/kg) si ∗

Feq (Var/kg) sunt pierderile specifice de putere activa,respectiv reactiva;Avand in vedere expresia curentului de mers in gol:

2

,

2

.0).( r oao iir ui += , unde oao pi =, (u.r.) si or o qi =, (u.r.)

(24)valoarea corespunzatoare a inductiei magnetice(BFe) poate rezulta din urmatoarea egalitate:

)()()()()( ,

22

, Fe P Fe Feno Fe Fe Fe Feno FeQ B f B pi Bq B p p B f =⋅=+⋅=∗∗∗

(25)

c) Diametrul coloanei(circulara): )/()4( ,cu Fe k A D ⋅⋅= π , unde g Feiz cu k k k ⋅= ,, ,(26)

- k iz,Fe ,factorul de izolatie a tolelor (0.9....0.95)

Page 6: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 6/18

- c Fe g A Ak /= ,unde 4/2

D Ac π = factorul geometric de umplere a coloanei,in functie de numarulde trepte ( 4/π .....1)

d) Dimensiunile principale (diametrul mediu D1,2 si inaltimea medie a infasurarilor H b):- dD D D +=2,1 (m), unde 12101 22 aaadD ++= , (27)

- β β π //12 wb L D H == ,unde Lw este lungimea medie a spirelor (de joasa si inalta tensiune)

e) Solicitari electromagnetice,si termice:-vectorul Poynting )/(

2,1 bc P H DS S ⋅= π (VA/m2), )/( 2,1 DU E

sp s π = (V/m), s P E S A H /==σ ,A/m- densitatea (medie) de curent in cele doua infasurari : )2/(

,ρ ⋅⋅= snk

E p J ,(A/mm2), (28)- fluxul termic pe suprafata de transfer(Poynting): A J Q pt ⋅⋅= ρ , ,

f) Dimensiunile ferestrei(FxH):-largimea ferestrei: 222,11, )(2 aaaa F

o+++= (m) ,unde 21 aaa += (29)

-inaltimea ferestrei: 2,2ob h H H += (m),

g) Consumuri de materiale active(cupru,fier)-consumul de cupru(aluminiu) in infasurari: 2

,/ J P d V d M nk CuCuCuCu ⋅⋅=⋅= ρ (kg)

-consumul specific de cupru: )/(/2 J d pS M c CuknnCuCu ⋅⋅== ρ (kg/VA) (30)

-consumul de fier -in cele 3 coloane: Fe Fec Fe

d H AM ⋅⋅=3, (kg),-in cele doua juguri: Fe Fe j Fe

d F D AM ⋅+⋅= )23(2, (kg),unde s-a considerat Fec Fe j Fe A A A == ,,

-consumul specific de fier: n j Fec Fe Fe S M M c /)( ,, += (kg/VA)

h) Costul specific al materialelor active

-costul specific al infasurarilor CuCuCu ck λ ⋅= , (u.m./VA) , (31)-costul specific al miezului magnetic Fe Fe Fe ck λ ⋅= , (u.m./VA),

unde Cuλ si Feλ (u.m./kg), reprezinta costul specific al materialului conductor(cupru saualuminiu),respectiv al miezului magnetic(tole)

i) Parametri de mers in gol (po ,io)- pierderile de putere activa, in miez: )( Fe Fe Fe Fe B pc p ∗

⋅= = aoo i p ,= (u.r.) (32)

- pierderile de putere reactiva, in miez: )( Fe Fe Fe Fe Bqcq ∗

⋅= r oi

,= (u.r.)

- curentul de mers in gol: 2,

2, r oaoo iii += (u.r.)

j) Criterii de optimizare

-criteriul parametrilor de mers in gol: nonoo po p p pd ., /)( −= < poε , nonoooi iiid ,.,, /)( −= < poε (33)

-criterii economice: FeCum ccc += ≅ min; FeCum k k k += ≅ min;2.2 Algoritmul PoyntingUn algoritm de dimensionare alternativ se poate construi avand in vedere relatia(10), dintre vectorulPoynting(SP), fluxul termic (Qt) pe suprafata de transfer(AP) a puterii electromagnetice-proportional cu sarcina termica A*J, si pierderile Joule nominale din infasurare( 2/,,2,1 nk n jn j

p p p =≅ ).

Page 7: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 7/18

Astfel,daca se impun solicitarea termica Qt(W/m2), si pierderile de scurtcircuit pk,n (u.r.),algoritmul Poynting corespunde urmatoarelor secvente specifice,in functie de factorul de geometrie :a) vectorul Poynting si elementele sale(independente de !): nk t P

pQS ,/2= (VA/m2), P P s Z S E ⋅= (V/m), s P E S A H /==σ (A/m), r k o P ua Z ,/

σ ωµ = (Ω) (32)

b) dimensiunile principale: P c P S S A /= (m2), β / P b A H = (m), π β /2,1 b H D ⋅= (m) si dD D D −= 12 (m) (33)

c) solicitarile electrice si magnetice,in infasurari si miez: ρ 2/

, snk E p J ⋅= (A/m2), )/( 2

,12 Dk D E B f c s Fe ⋅⋅= (T),unde 24/,, cu f c k k ⋅=ω (34)In continuare se determina dimensiunile ferestrei,consumul de materiale active, parametri de mers ingol,conform secventelor f,g,h,i din algoritmul precedent.Observatie: Utilizarea algoritmului Poynting conduce la solicitari electrice(J) constante,independentede factorul de geometrie ( ),respectiv consum de cupru constant(rel.28),astfel incat criteriile deoptimizare privesc exclusiv miezul magnetic(consumul de fier,pierderile in fier),unde inductiamagnetica este variabila(rel.34).

2.3 Date de proiectare orientative

a) distante de izolatie,in functie de tensiunea nominala a infasurarii (IT)

Un (kV) Uinc (kV) a01 (cm) a12 (cm) a22 (cm) h10=h20 (cm)≤ 1 5 0,5 0,9 1 1,86 25 1,5 1,2 1 1,810 35 1,8 1,5 1,4 2,515 45 2 1,7 1,7 320 55 2,3 1,8 2 4,535 85 3 2,7 3 7

b) factorul geometric de umplere a coloanei(k g),in functie de numarul de trepte

Sn(kVA) 5÷ 15 15÷ 45 45÷ 5000 5000÷ 20000nt 4 5 6 7kg 0,886 0,908 0,92 0,93

c)valori recomandate pentru inductia magnetica-table laminate la rece Sn (kVA) 10-50 75-250 315-500 630-1000 >1000BFe,c (T) 1.21-1.55 1.48-1.65 1.55-1.675 1.56-1.69 1.57-1.7

d) proprietati fizice materiale conductoare(cupru,aluminiu)

Propr. Fiz. ρ (Ω m)la 20oC

α(1/oC)

d(kg/m3)

σ(N/m2)

c(J/kg0C)

λ(W/moC)

Material

Cupru 1,75⋅ 10-8 3,9⋅ 10-3 8900 (2,1÷ 4,5)108 388 375

Page 8: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 8/18

Aluminiu3,1⋅ 10-8 3,7⋅ 10-3 2700 (0,7÷ 2,3)108 885 205

e)solicitari electromagnetice si termice,recomandate(infasurari din cupru,racire in ulei,cl.iz.A)

Sn kVA) 1~5 5~10 10~50 50~100 100~1000 1600~10000 >10000J(A/mm2) 2.35 2.5 2.65 2.8 2.9 3.1 3.7A(A/cm) 165 190 225 260 325 500 625Qt(W/m2) 820 1000 1250 1500 2000 3300 4900

Obs: Fluxul termic Qt s-a calculat cu relatia: A J QCut ⋅⋅=

75, ρ ,pentru Ω⋅=

−8

75,101.2

Cu ρ m,

f) factori de umplere ai infasurarilor

Fig.33. Program de calcul (date de proiectare, algoritmul Usp, solutia optima)

3.1 Program de calcul(matlab)

function TrafoSintezaOptimizareSAE2011

%S1 Sinteza optimala

1.bobinaj cilindricJOAS Ă TENSIUNE

2.bobinaj spiralat

200 400 600 800 1000 1200 16001400

0,2

0,4

0,6

0,8 (1) (2)525 V400 V 400 V525 V

3,6 kV

(Sn)

k u

kVA

200 400 600 800

NALT ATENSIUNE

1000 1200 16001400

0,2

0,4

0,6

0,8

(S n)

k u

kVA

0,5

3kV

6kV

15kV

30kV

Page 9: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 9/18

%S1.1 Date nominale

Sn=1000000

fn=50

U1n=525 %stea

U2n=10000 %triunghi

Pkn=10500 % W

ukn=0.05 %u.r.

Pon=1700 %W

ion=0.014 %u.r.

%date nominale prelucrate

pkn=Pkn/Sn

uka=pkn

ukr=sqrt(ukn^2-pkn^2)

pon=Pon/Sn

ioa=pon

ior=sqrt(ion^2-ioa^2)

%conditii de compatibilitate

Cpu=pkn/ukn %<1

Cpi=pon/ion %<1

%date material infasurari

ro1=3.7e-8%rezistivitate aluminiu,la 75 grade%ro1=2.1e-8;%rezistivitate cupru,la 75 grade

dcu=2700 %densitatea aluminiului

%dcu=8900 %densitate cupru

xcu=2%cost unitar,aluminiu (u.m./kg)

%xcu=3;%cost unitar,cupru (u.m./kg)

%S1.2.Date constructive(trafo trifazat cu inf.cilindr.conc)

%1.2.1 date constructive miez (trifazat,3 coloane)

Sc=Sn/3;%puterea electromagnetica,pe coloana

kizfe=0.95%factor de izolatie tole(oxizi ceramici

kg=0.92 %factor geometric de umplere a coloanei(nt=6)

kuc=kg*kizfe;%factor global de umplere a coloanei%proprietati material tole- tabla Armco M6 0.35 mm

dfe=7600;%densitate

xfe=1;%cost unitar(u.m/kg)

%caracteristici pierderi de putere activa pfet(b)-W/kg,si %reactiva,qfet(b)-VAr/kg

Bfe=[1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8];

pfet=[0.5 0.6 0.7 0.82 0.95 1.1 1.32 1.6 2];% W/kg,(vezi pg 95,carte)

qfet=[1.46 2.1 3.1 4.5 6.7 10.1 15.8 25 48]; %VA/kg inclusiv intrefieruri

%aproximare analitica curbe pierderi tabla laminata la rece

a1=.5 ;n1=2.3;%propunere pentru p10=0.5W/kg

ppfe=a1.*(Bfe.^n1);% expr analitica p10=0.44

b1=1.46; c1=0.0002; n2=2;n3=21; %propunere pentru q10=1.46VAR/kg

b1=1.46; c1=0.04; n2=2.2;n3=12; %varianta

qqfe=b1.*(Bfe.^n2)+c1.*(Bfe.^n3);%expr analitica

plot(Bfe,ppfe,'-',Bfe,pfet,'.')%verificare aproximare analiticaplot(Bfe,qqfe,'-',Bfe,qfet,'.')%verificare aproximare analitica

%S1.2.2 Date constructive infasurari(cilindrice concentrice,stratificat))

dcu=2700;%Aluminiu

%dcu=8900; %cupru

xcu=2;%cost unitar,aluminiu (u.m./kg)

%xcu=3;%cost unitar,cupru (u.m./kg)

ro1=3.7e-8;%rezistivitate aluminiu,la 75 grade

%ro1=2.1e-8;%rezistivitate cupru,la 75 grade

Page 10: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 10/18

ro2=ro1;

ku1=.75;%factor de umplere infasurare j.t.

ku2=.5;%factor de umplere infasurare i.t

%S1.2.3.Schema de izolatie(trafo in ulei,cl.iz.A)

a01=.018;a12=.015;a22=.014;h01=.025;h02=h01; %U2n=10kV

% calcul grosime(radiala) infasurari

miuo=pi*4e-7;

omega=2*pi*fn;

ku=2*ku1*ku2/(ku1+ku2);%factor mediu de umplere al celor 2 infasurari

RR=16*ro1*ukr/(3*miuo*2*pi*fn*pkn*ku);

as=(a12+sqrt(a12^2+RR))/2;%as=a12+(a1+a2)/3,grosime echiv.canal

scapari,inf.concentrice

a=3*(as-a12); %a=a1+a2 ,grosimea infasurarilor

%S2 Algoritmi de optimizare

%S.2.1.Algoritmul Usp (Afe)

%estimare inductie in fier Bfen

Bfe=[1:.05:1.8];

qqfe=b1.*Bfe.^n2+c1.*Bfe.^n3;

ppfe=a1.*Bfe.^n1;fp=ion.*ppfe;

fq=pon.*sqrt(qqfe.^2+ppfe.^2);

plot(Bfe,fp,'.',Bfe,fq,'-');

1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.80

0.01

0.02

0.03

0.04

0.05

0.06

0.07

0.08

0.09

%Fig M-1

Bfen=1.55;%pt fp=fq

%factor de geometrie variabil (beta=pi*D12/Hb=0.5....3.5)

for n=1:40

beta(n)=0.5+n*0.1 ;%se modifica beta=0.5-4.5

% tens de spira, aria sect de fier,diametrul coloanei

Usp(n)=sqrt(miuo*omega*Sc*as*beta(n)/(ukr));

Afe(n)=sqrt(2)*Usp(n)/(omega*Bfen);

D(n)=sqrt(4.*Afe(n)/(pi*kuc));

%dimensiuni principale si solicitari electromagnetice

dD=a+a12+2*a01;

D12(n)=D(n)+dD;%diametrul mediu al celor doua infasurari

Hb(n)=(pi*D12(n))/beta(n);%inaltimea medie a celor doua infasurari

Page 11: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 11/18

Ap(n)=(pi*D12(n))*Hb(n);%suprafata(aria)de transfer (Poynting)

Sp(n)=Sc/Ap(n);%densitatea de putere (vectorul Poynting,VA/m2)

Es(n)=Usp(n)/(pi*D12(n));%intensitatea campului electric (de transformare)V/m

Hs(n)=Sp(n)/Es(n);%intensitatea campului magnetic(de scapari),A/m

Zp(n)=Es(n)/Hs(n);%impedanta Poynting

J(n)=(pkn*Es(n))/(2*ro1);%densitatea medie de curent,

Qt1(n)=(ro1*J(n))*Hs(n);%fluxul termic unitar,pe suprafata Poynting

plot(beta,D12,'+',beta,Hb,'.',beta,D,'*')

0.5 1 1.5 2 2.5 3 3.5 4 4.50

0.2

0.4

0.6

0.8

1

1.2

1.4

1.6

1.8

% Fig.M2

%dimensiuni fereastra(F*H)

F(n)=(2*(a+a01+a12)+a22)*n/n; %largimea ferestrei, constanta

H(n)=Hb(n)+h01+h02; %inaltimea ferestrei

Lj(n)=2*F(n)+3*D(n); %lungimea jugului

%consumuri si costuri specifice fier/cupru

Mfec(n)=3*dfe*Afe(n)*H(n);%kg

Mfej(n)=2*dfe*Lj(n)*Afe(n);%Afej=Afec

cfec(n)=Mfec(n)/Sn;%kg/VA

cfej(n)=Mfej(n)/Sn;

cfe(n)=cfec(n)+cfej(n);%consum specific de fier(tole) kg/VA

kfe(n)=xfe*cfe(n);%costul specific al fierului, u.m/VA

ccu(n)=dcu*pkn/(ro1*J(n)^2);%consum specific de cupru(aluminiu),kg/VA

kcu(n)=xcu*ccu(n);%cost specific cupru(aluminiu)

cm(n)=cfe(n)+ccu(n);%consum specific de materiale active(fier,cupru)

km(n)=kfe(n)+kcu(n);%cost specific de materiale active

plot(beta,cm,'-',beta,km,'+')

Page 12: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 12/18

0.5 1 1.5 2 2.5 3 3.5 4 4.51.6

1.7

1.8

1.9

2

2.1

2.2

2.3

2.4

2.5x10

-3

%Fig.M3

% pierderi in fier si putere de magnetizare

ppfen=a1*Bfen^n1;%pierderi specifice,putere activa,in miez,W/kg

qqfen=b1*Bfen^n2+c1*Bfen^n3;%pierderi specifice,putere reactiva,VAr/kg

pfe(n)=cfe(n)*ppfen;%pierderi de putere activa in fier,in u.r.,(Po/Sn))qfe(n)=cfe(n)*qqfen;%pierderi de putere reactiva(magnetizare),u.r.(Qo/Sn))

io(n)=sqrt(pfe(n)^2+qfe(n)^2);%curent de magnetizare in u.r.(I1o/I1n)

po(n)=pfe(n);%pierderi de mers in gol,in u.r.(Pon/Sn)

%f)criterii de optimizare

dpo(n)=(po(n)-pon)/pon;%eroare relativa,fata de valoarea nominala,pon,impusa

dio(n)=(io(n)-ion)/ion;%eroare relativa,fata de valoarea nominala,ion,impusa

plot(beta,dpo,'-',beta,dio,'.')

plot(beta,cm,'-',beta,km,'.') %fig.M3

0.5 1 1.5 2 2.5 3 3.5 4 4.5-0.1

0

0.1

0.2

0.3

0.4

0.5

0.6

% Fig.M4

grid

end 'S.3 Solutia optima'

% Sol.opt. Alg.Usp(solutia optima,dpo,dio < epsipo,epsiip,sau cm sau km =min)'

%epsipo =0.15, epsiio=0,2,erori admisibile,fata de parametri de mers in

%gol(pon,ion)impusi

betax=1.7 %cm,km=min

nx=round((betax-0.5)/0.1)%nr ordine ciclul beta

'solicitari electromagnetice si termice'

Page 13: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 13/18

Bfe=Bfen %(T)

Spx=Sp(nx) %(VA/m2)

Esx=Es(nx) %(V/m)

Hsx=Hs(nx) %(A/m)

Zpx=Es(nx)/Hs(nx) %ohm

Jx=J(nx) %A/m2

Qtx=Qt1(nx) %(W/m2

'dimensiuni principale'

Apx=Ap(nx) %m2

D12x=D12(nx)%m

Hbx=Hb(nx) %m

epsiAp=(pi*D12x*Hbx-Apx)/Apx %u.r.

Dx=D(nx) %m

Afex=Afe(nx)%m2

Uspx=Usp(nx) %V

'dimensiuni fereastra'

Fx=F(nx) %m

Hx=H(nx) %m

Ljx=Lj(nx)%m

'consumuri specifice materiale active'

ccux=ccu(nx) *1000 %kg/kVAcfex=cfe(nx)*1000 %kg/kVA

'parametri de gol'

pox=po(nx) %u.r.

epsipox=(pox-pon)/pon %u.r.

iox=io(nx) %in u.r.

epsiiox=(iox-ion)/ion %u.r.

'parametri de scc'

pkx=ccu(nx)*ro1*Jx^2/dcu

epsipkx=(pkx-pkn)/pkn

ukrx=omega*miuo*as/Zpx

ukax=pkx

ukx=sqrt(ukax^2+ukrx^2)

epsiukx=(ukx-ukn)/ukn

3.2 Rezultate numerice (Aplicatii)

%observatii:

1.se va nota semnificatia marimilor,unitatea de masura,interpretari,etc

2.se vor specifica si interpreta marimile din graficele anexate(fig.M1,M2,M3,M4)

3.se va comenta alegerea solutiei optime (fig.M3,M4),respectiv factorul beta optim

4.se va desena si cota,pentru solutia optima, schita reprezentativa din

Fig.1(miez magnetic,infasurari,dimensiuni de gabarit)

DATE NOMINALE

Sn =

1000000

fn =

50

Page 14: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 14/18

U1n =

525

U2n =

10000

Pkn =

10500

ukn =

0.0500

Pon =

1700

ion =

0.0140

pkn =

0.0105

uka =

0.0105

ukr =

0.0489

pon =

0.0017

ioa =

0.0017

ior =

Page 15: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 15/18

0.0139

Cpu =

0.2100

Cpi =

0.1214

ro1 =

3.7000e-008

dcu =

2700

xcu =

2

kizfe =

0.9500

kg =

0.9200

SOLUTIA OPTIMA

betax =

1.7000

nx =

12

1.solicitari electromagnetice si termice

Bfe =

1.5500

Page 16: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 16/18

Spx =

4.3540e+005

Esx =

11.7896

Hsx =

3.6931e+004

Zpx =

3.1923e-004

Jx =

1.6728e+006

Qtx =

2.2858e+003

2.dimensiuni principale

Apx =

0.7656

D12x =

0.3631

Hbx =

0.6711

epsiAp =

0

Dx =

0.2385

Page 17: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 17/18

Afex =

0.0391

Uspx =

13.4499

3.dimensiuni fereastra

Fx =

0.2272

Hx =

0.7211

Ljx =

1.1700

4.consumuri specifice materiale active

ccux =

0.2738

cfex =

1.3369

5. parametri de gol

pox =

0.0018

epsipox =

0.0774

iox =

Page 18: Transformatorul,Sinteza1,optimizare,2011

5/11/2018 Transformatorul,Sinteza1,optimizare,2011 - slidepdf.com

http://slidepdf.com/reader/full/transformatorulsinteza1optimizare2011 18/18

0.0155

epsiiox =

0.1079

6.parametri de scc

pkx =

0.0105

epsipkx =

-1.6521e-016

ukrx =

0.0489

ukax =

0.0105

ukx =

0.0500

epsiukx =

2.7756e-016


Recommended