Date post: | 10-Jul-2015 |
Category: |
Documents |
Upload: | ana-nistor |
View: | 74 times |
Download: | 0 times |
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.
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.
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)
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
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)
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 =≅ ).
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
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
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
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
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,'+')
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'
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
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 =
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
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
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 =
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