03.01.2013 Views

+még1 bonyolultabb

+még1 bonyolultabb

+még1 bonyolultabb

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

4. gyakorlat<br />

Zárlatszámítás szimmetrikus<br />

összetevők módszerével<br />

A gyakorlat tartalma:<br />

− viszonylagos egységek használata<br />

− 1FN (2F zárlatszámítás)<br />

Feltételezések:<br />

− 50 Hz, kvázi-stacionárius állapot<br />

Hálózat adatai, a modell paraméterei saját feszültségszinten<br />

Adott az alábbi hálózat, amely az egyszerűség kedvéért megegyezik a 3. gyakorlatban már látottal:<br />

A<br />

d<br />

B C D E<br />

H Yg y<br />

D y<br />

Mögöttes<br />

hálózati<br />

táppont<br />

Un = 120 kV<br />

Sz = 1200 MVA<br />

Földelt csp.<br />

T1<br />

transzformátor<br />

UN1/UK1 =<br />

132/22 kV<br />

ε1 = 10%<br />

Sn1 = 25 MVA<br />

Ygy6+d<br />

20 kV<br />

szabadvezeték<br />

Lv1 =10 km<br />

r1 = 0,3 Ω/km<br />

x1 = 0,3 Ω/km<br />

T2 transzformátor<br />

UN2/UK2 = 20/0,4 kV<br />

ε2 = 5%<br />

Sn2 = 160 kVA<br />

Dyg5<br />

0,4 kV szabadvezeték<br />

Lv2 =1 km<br />

r2 = 0,4 Ω/km<br />

I. körzet II. körzet III. körzet<br />

Modell elemei, saját feszültségszinten (transzformátoroké a kisebb feszültségen):<br />

Mögöttes<br />

hálózati<br />

táppont<br />

T1<br />

transzformátor<br />

20 kV<br />

szabadvezeték<br />

T2 transzformátor<br />

0,4 kV szabadvezeték<br />

UG , jXH jXTRK1 Rv1 + jXv1 jXTRK2 Rv2<br />

UG = Un / √3 = 120 kV / √3 = 69,3 kV<br />

jXH = j Un 2 / Sz = j (120 kV) 2 /1200 MVA = j 12 Ω<br />

jXTRK1 = j UK1 2 / Sn1∙ε1 /100 = j (22 kV) 2 /25 MVA * 0,1 = j 1,94 Ω<br />

Rv1 + jXv1 = Lv1 * (r1 + jx1) = 10 km * (0,3 + j 0,3) Ω/km = 3 + j 3 Ω<br />

jXTRK2 = j UK2 2 / Sn2∙ε2 /100 = j (0,4 kV) 2 /0,160 MVA * 0,05 = j 0,05 Ω<br />

Rv2 = Lv2 * r2 = 1 km * 0,4 Ω/km = 0,4 Ω<br />

Legyen minden hálózati elem negatív sorrendű impedanciája azonos a pozitív sorrendűvel.<br />

Legyen továbbá a távvezetékek zérus sorrendű impedanciája a pozitív sorrendűnek a duplája, és<br />

legyen a háromtekercselésű transzformátor zérus sorrendű impedanciája a pozitív sorrendűnek a<br />

másfélszerese. A többi elem zérus sorrendű impedanciája egyezzen meg a pozitív sorrendűvel.<br />

1


A hálózat sorrendi modellje viszonylagos egységekben<br />

Hálózat pozitív, negatív és zérus sorrendű modelljének elemei, viszonylagos egységekben:<br />

+<br />

-<br />

0<br />

Válasszuk a következő alapmennyiségeket: Ua (kif) = 0,4 kV (vonali) , Sa = 160 kVA<br />

Innen a többi alapegység (ld. a 3. gyakorlat anyagában):<br />

Ua (köf) = Ua (kif) * UN2/UK2 = 0,4 kV * 20/0,4 = 20 kV<br />

Ua (nf)<br />

= Ua (köf) * UN1/UK1 = 20 kV * 132/22 = 120 kV<br />

= [Ua (kif) ] 2 / Sa = 1 Ω<br />

Za (köf) = [Ua (köf) ] 2 / Sa = 2,5 kΩ<br />

Za (nf)<br />

= [Ua (nf) ] 2 / Sa = 90 kΩ<br />

Za (kif)<br />

Ia (kif)<br />

Ia (köf)<br />

Ia (nf)<br />

UG (ve) jXH (ve) jXTRK1 (ve) Rv1 (ve) + jXv1 (ve) jXTRK2 (ve) Rv2 (ve)<br />

~ UA1 UB1 UC1 UD1 UE1<br />

= Sa / [√3 Ua (kif) ] = 230,9 A<br />

= Sa / [√3 Ua (köf) ] = 4,620 A<br />

= Sa / [√3 Ua (nf) ] = 0,770 A<br />

Alapmennyiség \<br />

feszültségszint<br />

jXH (ve) jXTRK1 (ve) Rv1 (ve) + jXv1 (ve) jXTRK2 (ve) Rv2 (ve)<br />

NF Köf kif<br />

Sa Sa = 160 kVA Sa = 160 kVA Sa = 160 kVA<br />

Ua 120 kV 20 kV Ua (kif) = 0,4 kV<br />

Za 90 kΩ 2,5 kΩ 1 Ω<br />

Ia 0,770 A 4,620 A 230,9 A<br />

Ezek segítségével a modell paraméterei:<br />

UA2 UB2 UC2 UD2 UE2<br />

jXH (ve) j 1,5 XTRK1 (ve) 2 (Rv1 (ve) + jXv1 (ve) ) jXTRK2 (ve) 2 Rv2 (ve)<br />

X X<br />

UA0 UB0 UC0 UD0 UE0<br />

UG (ve) = UG / (Ua (nf) / √3) = 69,3 kV / (120/ √3)kV = 1 („/√3”, mert a vonali alapot számoltuk ki eddig)<br />

jXH (ve) = jXH / Za (nf) = = j 0,133∙10 -3<br />

jXTRK1 (ve) = jXTRK1 / Za (kof) = = j 0,774∙10 -3<br />

Rv1 (ve) + jXv1 (ve) = (Rv1 + jXv1) / Za (kof) = (1,2 + j 1,2)∙10 -3<br />

jXTRK2 (ve) = jXTRK2 / Za (kif) = j 0,05 Ω / 1 Ω = j 0,05<br />

Rv2 (ve) = Rv2 / Za (kif) = 0,4 Ω / 1 Ω = 0,4<br />

2


Különböző zárlatok leképezése a sorrendi modellben<br />

Ih0,1,2: a hibahelyi áram szimmetrikus összetevői<br />

1FN(a) zárlat az „E” kisfeszültségű végponton<br />

+<br />

-<br />

0<br />

1FN(a) zárlat a „C” csomóponton<br />

+<br />

-<br />

0<br />

~ UA1 UB1 UC1 UD1 UE1<br />

UA2 UB2 UC2 UD2 UE2<br />

X X<br />

UA0 UB0 UC0 UD0 UE0<br />

Látható, hogy ebben az esetben a zárlati áram nulla lesz.<br />

Ih1<br />

~ UA1 UB1 UC1 UD1 UE1<br />

Ih2<br />

UA2 UB2 UC2 UD2 UE2<br />

Ih0<br />

X X<br />

UA0 UB0 UC0 UD0 UE0<br />

Ih1<br />

Ih2<br />

Ih0<br />

3


2F (b,c) zárlat az „E” csomóponton<br />

+<br />

-<br />

0<br />

Látható, hogy a hibahelyi áram zérus sorrendű összetevője nulla: Ih0 = 0<br />

2F (b,c) zárlat a „C” csomóponton<br />

+<br />

-<br />

0<br />

~ UA1 UB1 UC1 UD1 UE1<br />

UA2 UB2 UC2 UD2 UE2<br />

X X<br />

UA0 UB0 UC0 UD0 UE0<br />

Látható, hogy a hibahelyi áram zérus sorrendű összetevője nulla: Ih0 = 0<br />

Ih1<br />

~ UA1 UB1 UC1 UD1 UE1<br />

Ih2<br />

UA2 UB2 UC2 UD2 UE2<br />

X X<br />

Ih0<br />

UA0 UB0 UC0 UD0 UE0<br />

Ih1<br />

Ih2<br />

Ih0<br />

4


2FN (b,c) zárlat az „E” csomóponton<br />

+<br />

-<br />

0<br />

Látható, hogy a hibahelyi áram zérus sorrendű összetevője NEM nulla: Ih0 ≠ 0<br />

2FN (b,c) zárlat a „C” csomóponton<br />

+<br />

-<br />

0<br />

~ UA1 UB1 UC1 UD1 UE1<br />

UA2 UB2 UC2 UD2 UE2<br />

X X<br />

UA0 UB0 UC0 UD0 UE0<br />

Látható, hogy a hibahelyi áram zérus sorrendű összetevője nulla: Ih0 = 0<br />

Ih1<br />

~ UA1 UB1 UC1 UD1 UE1<br />

Ih2<br />

UA2 UB2 UC2 UD2 UE2<br />

X X<br />

Ih0<br />

UA0 UB0 UC0 UD0 UE0<br />

Ih1<br />

Ih2<br />

Ih0<br />

5


Zárlatszámítás: 1FN(a) zárlat az „E” kisfeszültségű végponton<br />

Számítsuk ki az „A”, „C” és „D” pontokon a fázisfeszültségeket, valamint az itt folyó áramokat (a<br />

tápponttól a végpont felé felvett pozitív iránnyal).<br />

A „ve” jelölést a továbbiakban az egyszerűség kedvéért elhagyjuk.<br />

�� ��,��,�� � �� ∑ � �<br />

�<br />

2���� � ������ � ��� � ���� � ������ � ���� � 2��� � ������ 1<br />

�<br />

� 0,618 � �0,0595 � 0,621 � � 5,5°<br />

1,602 � �0,154<br />

UA1 = UG1 – jXH Ih1 = 1 – j 0,133∙10 -3 0,621 � � 5,5° = 1<br />

UA2 = UG2 – jXH Ih2 = 0 – j 0,133∙10 -3 0,621 � � 5,5° = 0<br />

UA0 = UG2 – jXH 0 = 0 – 0 = 0<br />

UC1 = UG1 – (jXH + jXTRK1 + Rv1 + jXv1) Ih1 = 0,999 � � 0,07°<br />

UC2 = 0 – (jXH + jXTRK1 + Rv1 + jXv1) Ih2 = 0,0015 � � 125°<br />

UC0 = = 0<br />

UD1 = UG1 – (jXH + jXTRK1 + Rv1 + jXv1 + jXTRK2) Ih1 = 0,997 � � 1,8°<br />

UD2 = 0 – (jXH + jXTRK1 + Rv1 + jXv1 + jXTRK2) Ih2 = 0,032 � � 96,8°<br />

UD0 = 0 – jXTRK2 Ih0 = 0,031 � � 95,5°<br />

Mint az sejthető volt, a NF és Köf oldali feszültségekre semmilyen hatással nincs a kif végponti zárlat,<br />

és a kif tápponti feszültséget is alig befolyásolja.<br />

A zárlati (hibahelyi) áram: Ih = IEa = Ih0 + Ih1 + Ih2 = 3Ih0,h1,h2 = 1,86 � � �, �°<br />

[ Ennek dimenzionális értéke Ih,dim = Ih Ia (kif) = 230,9 � � �1,86 � � 5,5°) = 430,4 � � � 5,5° ]<br />

Ha a hálózat különböző pontjain a fázisfeszültségekre vagy a fázisáramokra is kíváncsiak vagyunk,<br />

akkor figyelembe kell vennünk a transzformátorok kapcsolási csoportját, forgatását.<br />

A k kapcsolási óraszám definíciója, ha U N+ és U K+ a transzformátor nagyobb ill. kisebb feszültségű<br />

oldalán mérhető pozitív sorrendű feszültség viszonylagos egységben (áramokra hasonlóan):<br />

� �� � � �� � �����°<br />

Negatív sorrendű feszültségekre, viszonylagos egységben (és áramokra hasonlóan):<br />

� �� � � �� � �����°<br />

Zérus sorrendű feszültségekre viszonylagos egységben (és áramokra hasonlóan), ha k = {2,6,10}:<br />

� �� � � �� � �����°<br />

Jelen esetben referenciának a kif oldalt választjuk, mert ott lépett fel a zárlat, ezért ott mindig<br />

teljesülnie kell a sorrendi modellek összekötése során már teljesített<br />

UE1+UE2+UE0 = 0 és Ih1 = Ih2 = Ih0<br />

egyenleteknek. (Ha ott elforgatnánk egymáshoz képest a pozitív és a negatív sorrendű<br />

mennyiségeket, akkor a fenti egyenlőségek nem teljesülnének.)<br />

Tehát a kif feszültségszinten az �� � feszültséget és az �� áramot nem forgatjuk el.<br />

A NF és Köf feszültségszinteken az alábbiak szerint kell eljárni:<br />

UC1,forg = UC1 e +j150° = 0,999 �149,9°<br />

UC2,forg = UC2 e - j150° = 0,0015 � � 85°<br />

UC0,forg = UC0 e + j180° = 0<br />

� �<br />

6


IC1,forg = Ih1 e +j150° = 0,621 �144,5°<br />

IC2,forg = Ih2 e - j150° = 0,621 � � 155,5°<br />

UA1,forg = UA1 e +j150° e +j180° = 1 � � 30°<br />

UA2,forg = UA2 e - j150° e - j180° = 0<br />

IA1,forg = Ih1 e +j150° e +j180° = 0,621 � � 35,5°<br />

IA2,forg = Ih2 e - j150° e - j180° = 0,621 � 24,5°<br />

Innen a viszonylagos egységben számított fázismennyiségek kiszámíthatók:<br />

�<br />

�<br />

�<br />

�<br />

��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� 1 1 1<br />

� � �1<br />

�� �<br />

1 � �� � �<br />

1 1 1<br />

� � �1<br />

�� �<br />

1 � �� � �<br />

1 1 1<br />

� � �1<br />

�� �<br />

1 � �� � �<br />

1 1 1<br />

� � �1<br />

�� �<br />

1 � �� � �<br />

��� ���,���� ���,���� 0<br />

���,���� ���,���� 0<br />

���,���� ���,���� 0<br />

���,���� ���,���� 0,9998 � � 149,85°<br />

� � � 0,9976 � � 29,9°<br />

1 � � 90°<br />

�<br />

1,0759 � ���, �°<br />

� � � 1,0759 � � �, �° �<br />

0<br />

1 � � 30°<br />

� � �1<br />

� � 150° �<br />

1 � 90°<br />

1,0759 � � 5,5°<br />

� � �1,0759<br />

� 174,5° �<br />

0<br />

A viszonylagos egységben értett fázismennyiségekből az adott helyre (csomópontra) vonatkozó<br />

alapmennyiségekkel történő szorzással megkaphatók a dimenzionális értékek.<br />

Pl. UCa,dim = UCa Ua (köf) /√3 = 0,9998 � � 149,85° *20kV/√3 = 11,544 �� � � 149,85°<br />

ICa,dim = ICa Ia (köf = 1,0759 � 174,5° *4,62A = 4,97 � � 174,5°<br />

Érdemes megfigyelni, hogy ICc = 0. Ennek fizikai magyarázatát a T2 transzformátor kapcsolási<br />

csoportjának segítségével érthetjük meg. Először lássuk be, hogy az alábbi transzformátor Dyg5-ös,<br />

kövessük nyomon a fázisfeszültségek alakulását:<br />

A<br />

B<br />

C<br />

C c<br />

A<br />

Ezt követően hozzunk létre 1FN(a) földzárlatot a kisebb feszültségű oldalon:<br />

B<br />

b<br />

b<br />

a<br />

a<br />

c<br />

7


A<br />

B<br />

C c<br />

Látható, hogy a csillag oldalon csak egy tekercsben folyhat áram, és transzformátor csak akkor lehet<br />

mágnesesen kiegyenlített, hogyha a delta oldalon ennek megfelelő ellengerjesztés alakul ki, vagyis az<br />

A és a B fázisokban folyik csak áram. A delta oldali áramok olyanok, mintha a transzformátor delta<br />

oldali kapcsain egy 2F(A,B) zárlat következett volna be.<br />

(Megjegyzés: Számításaink során az ICa és az Ih = IEa áramok ellentétes irányúaknak adódtak. Ennek<br />

oka az, hogy a számításhoz az egyszerűség kedvéért minden gyűjtősínen/csomópontban „balról<br />

jobbra” vettük fel az áramok pozitív irányát, ld. Ih1, Ih2 és Ih0 a korábbi ábrákon. Azonban a<br />

transzformátor „természetes” referenciairányai szimmetrikusak, pl.:<br />

TR<br />

a<br />

b<br />

Ennek oka pedig a transzformátor reciprocitása és<br />

szimmetriája, vagyis az áttételtől eltekintve egyformán<br />

kell viselkednie mindkét oldalról nézve.)<br />

8


Zárlatszámítás: 2FN(b,c) zárlat a „C” csomóponton<br />

�� �� � �� ∑ � �<br />

�� 2���� � ������ � ��� � ����� �<br />

1<br />

� 206,2 � � 60,3°<br />

�2,4 � �4,215�10�� (Itt még mindig mindent viszonylagos egységben számolunk!!)<br />

Ih2 = -Ih1 , Ih0 = 0<br />

UA1 = UG1 – jXH Ih1 = 0,9762 � � 0,79°<br />

UA2 = 0 – jXH Ih2 = 0,0275 � 29,7°<br />

UA0 = = 0<br />

UC1 = UG1 – (jXH + jXTRK1 + Rv1 + jXv1) Ih1 = 0,5 � 0°<br />

UC2 = 0 – (jXH + jXTRK1 + Rv1 + jXv1) Ih2 = UC1<br />

UC0<br />

UD1 = UC1<br />

UD2 = UC2<br />

UD0<br />

= UC1 (ha nem 2FN, hanem csak 2F zárlat lenne, akkor itt 0 állna!)<br />

= 0<br />

Jelen esetben referenciának a Köf oldalt választjuk, mert ott lépett fel a zárlat, ezért ott mindig<br />

teljesülnie kell a sorrendi modellek összekötése során már teljesített<br />

(UC1=UC2=UC0) és (Ih1 + Ih2 + Ih0 = 0)<br />

egyenleteknek. (Ha ott elforgatnánk egymáshoz képest a pozitív és a negatív sorrendű<br />

mennyiségeket, akkor a fenti egyenlőségek nem teljesülnének.)<br />

UA1,forg = UA1 e +j180° = 0,9762 � 179,2°<br />

UA2,forg = UA2 e - j180° = 0,0275 � � 150,3°<br />

IA1,forg = Ih1 e +j180° = 206,2 � 119,6°<br />

IA2,forg = Ih2 e - j180° = 206,2 � � 60,3°<br />

UD1,forg = UD1 e -j150° = 0,5 � � 150°<br />

UD1,forg = UD1 e -j150° = 0,5 � 150°<br />

�<br />

�<br />

�<br />

�<br />

�<br />

��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� ��� 1 1 1<br />

� � �1<br />

�� �<br />

1 � �� �, �<br />

� � � � � 0 �<br />

1 1 1<br />

0<br />

� � �1<br />

�� 1 �<br />

�<br />

�� 0<br />

0<br />

1<br />

� � �1<br />

1<br />

�<br />

� �����<br />

� ����<br />

� � ���, �° �<br />

��� ��� � ��, �°<br />

1<br />

� �<br />

1 � �� 0<br />

�1<br />

� ����,����<br />

� � � 0.9768 � 57,6° �<br />

���,���� 0.9523 � � 60,0°<br />

1<br />

� � �1<br />

1<br />

�<br />

1<br />

� �<br />

1 � �� 0<br />

� ����,����<br />

� � �<br />

0<br />

��� � ��, �° �<br />

��� � � ���, �°<br />

1 1 1<br />

� � �1<br />

�� �<br />

1 � �� � �<br />

� ��,����<br />

0<br />

���,���� ���,���� � � �<br />

��, ���<br />

0<br />

�, ���<br />

�<br />

9


Matlab kód<br />

%% --------------- Alapadatok ---------------<br />

UnH = 120e3; SzH = 1200e6;<br />

UN1 = 132e3; UK1 = 22e3; eps1 = 10; Sn1 = 25e6;<br />

Rv1 = 3; Xv1 = 3;<br />

UN2 = 20e3; UK2 = 0.4e3; eps2 = 5; Sn2 = 160e3;<br />

Rv2 = 0.4;<br />

UG = UnH/sqrt(3);<br />

XH = UnH^2/SzH;<br />

XTR1K1 = UK1^2/Sn1*eps1/100;<br />

XTR2K2 = UK2^2/Sn2*eps2/100;<br />

a = exp(1j*120*pi/180); a2 = exp(-1j*120*pi/180);<br />

T = [1 1 1; 1 a2 a; 1 a a2 ]; % szimm.ot. ==> fázis<br />

%% --------------- Viszonylagos egységek ---------------<br />

% --------------- ALAPOK ---------------<br />

Sa = 160e3;<br />

Uakif = 0.4e3;<br />

Uakof = Uakif * UN2/UK2<br />

Uanf = Uakof * UN1/UK1<br />

Zakif = Uakif^2 / Sa<br />

Zakof = Uakof^2 / Sa<br />

Zanf = Uanf^2 / Sa<br />

Iakif = Sa / (sqrt(3) * Uakif)<br />

Iakof = Sa / (sqrt(3) * Uakof)<br />

Ianf = Sa / (sqrt(3) * Uanf)<br />

% --------------- MODELL paraméterek ---------------<br />

UGve = UG / (Uanf / sqrt(3))<br />

XHve = XH / Zanf<br />

XTR1K1ve = XTR1K1 / Zakof<br />

Zv1ve = (Rv1 + 1j*Xv1) / Zakof<br />

XTR2K2ve = XTR2K2 / Zakif<br />

Rv2ve = Rv2 / Zakif<br />

%% ------------------ 1FN(a) : E -------------------<br />

I = UGve/(2*(1j*XHve + 1j*XTR1K1ve + Zv1ve + 1j*XTR2K2ve + Rv2ve) +<br />

1j*XTR2K2ve + 2*Rv2ve), [abs(I), 180/pi*angle(I)]<br />

UA1 = UGve - 1j*XHve * I, [abs(UA1), 180/pi*angle(UA1)]<br />

UA2 = 0 - 1j*XHve * I, [abs(UA2), 180/pi*angle(UA2)]<br />

UA0 = 0 - 1j*XHve * 0, [abs(UA0), 180/pi*angle(UA0)]<br />

UC1 = UGve - (1j*XHve + 1j*XTR1K1ve + Zv1ve) * I, [abs(UC1), 180/pi*angle(UC1)]<br />

UC2 = 0 - (1j*XHve + 1j*XTR1K1ve + Zv1ve) * I, [abs(UC2), 180/pi*angle(UC2)]<br />

UC0 = 0<br />

UD1 = UGve - (1j*XHve + 1j*XTR1K1ve + Zv1ve + 1j*XTR2K2ve) * I<br />

[abs(UD1), 180/pi*angle(UD1)]<br />

UD2 = 0 - (1j*XHve + 1j*XTR1K1ve + Zv1ve + 1j*XTR2K2ve) * I<br />

[abs(UD2), 180/pi*angle(UD2)]<br />

UD0 = 0 - (1j*XTR2K2ve) * I, [abs(UD0), 180/pi*angle(UD0)]<br />

10


UC1forg = UC1 * exp(1j*150*pi/180), [abs(UC1forg), 180/pi*angle(UC1forg)]<br />

UC2forg = UC2 * exp(-1j*150*pi/180), [abs(UC2forg), 180/pi*angle(UC2forg)]<br />

IC1forg = I * exp(1j*150*pi/180), [abs(IC1forg), 180/pi*angle(IC1forg)]<br />

IC2forg = I * exp(-1j*150*pi/180), [abs(IC2forg), 180/pi*angle(IC2forg)]<br />

UA1forg = UA1 * exp(1j*150*pi/180) * exp(1j*180*pi/180)<br />

[abs(UA1forg), 180/pi*angle(UA1forg)]<br />

UA2forg = UA2 * exp(-1j*150*pi/180)* exp(-1j*180*pi/180)<br />

[abs(UA2forg), 180/pi*angle(UA2forg)]<br />

IA1forg = I * exp(1j*150*pi/180) * exp(1j*180*pi/180)<br />

[abs(IA1forg), 180/pi*angle(IA1forg)]<br />

IA2forg = I * exp(-1j*150*pi/180)* exp(-1j*180*pi/180)<br />

[abs(IA2forg), 180/pi*angle(IA2forg)]<br />

vmi = T * [UC0; UC1forg; UC2forg];<br />

UCa = vmi(1); UCb = vmi(2); UCc = vmi(3);<br />

[abs(UCa), 180/pi*angle(UCa)], [abs(UCb), 180/pi*angle(UCb)], [abs(UCc), 180/pi*angle(UCc)]<br />

vmi = T * [0; IC1forg; IC2forg];<br />

ICa = vmi(1); ICb = vmi(2); ICc = vmi(3);<br />

[abs(ICa), 180/pi*angle(ICa)], [abs(ICb), 180/pi*angle(ICb)], [abs(ICc), 180/pi*angle(ICc)]<br />

vmi = T * [UA0; UA1forg; UA2forg];<br />

UAa = vmi(1); UAb = vmi(2); UAc = vmi(3);<br />

[abs(UAa), 180/pi*angle(UAa)], [abs(UAb), 180/pi*angle(UAb)], [abs(UAc), 180/pi*angle(UAc)]<br />

vmi = T * [0; IA1forg; IA2forg];<br />

IAa = vmi(1); IAb = vmi(2); IAc = vmi(3);<br />

[abs(IAa), 180/pi*angle(IAa)], [abs(IAb), 180/pi*angle(IAb)], [abs(IAc), 180/pi*angle(IAc)]<br />

UCadim = UCa * Uakof/sqrt(3); [abs(UCadim), 180/pi*angle(UCadim)]<br />

ICadim = ICa * Iakof; [abs(ICadim), 180/pi*angle(ICadim)]<br />

%% ------------------ 2F(b,c) : C ------------------<br />

Ih1 = UGve/(2*(1j*XHve + 1j*XTR1K1ve + Zv1ve ) ), [abs(Ih1), 180/pi*angle(Ih1)]<br />

Ih2 = -Ih1; Ih0 = 0;<br />

UA1 = UGve - 1j*XHve*Ih1, [abs(UA1), 180/pi*angle(UA1)]<br />

UA2 = 0 - 1j*XHve*Ih2, [abs(UA2), 180/pi*angle(UA2)]<br />

UA0 = 0;<br />

UC1 = UGve - (1j*XHve + 1j*XTR1K1ve + Zv1ve) * Ih1, [abs(UC1), 180/pi*angle(UC1)]<br />

UC2 = 0 - (1j*XHve + 1j*XTR1K1ve + Zv1ve) * Ih2, [abs(UC2), 180/pi*angle(UC2)]<br />

UC0 = UC1;<br />

UD1 = UC1;<br />

UD2 = UC2;<br />

UD0 = 0;<br />

UD1forg = UD1 * exp(-1j*150*pi/180), [abs(UD1forg), 180/pi*angle(UD1forg)]<br />

UD2forg = UD2 * exp(+1j*150*pi/180), [abs(UD2forg), 180/pi*angle(UD2forg)]<br />

UA1forg = UA1 * exp(1j*180*pi/180), [abs(UA1forg), 180/pi*angle(UA1forg)]<br />

UA2forg = UA2 * exp(-1j*180*pi/180), [abs(UA2forg), 180/pi*angle(UA2forg)]<br />

IA1forg = Ih1 * exp(1j*180*pi/180), [abs(IA1forg), 180/pi*angle(IA1forg)]<br />

IA2forg = Ih2 * exp(-1j*180*pi/180), [abs(IA2forg), 180/pi*angle(IA2forg)]<br />

vmi = T * [UC0; UC1; UC2];<br />

UCa = vmi(1); UCb = vmi(2); UCc = vmi(3);<br />

[abs(UCa), 180/pi*angle(UCa)], [abs(UCb), 180/pi*angle(UCb)], [abs(UCc),<br />

180/pi*angle(UCc)]<br />

11


vmi = T * [0; Ih1; Ih2];<br />

ICa = vmi(1); ICb = vmi(2); ICc = vmi(3);<br />

[abs(ICa), 180/pi*angle(ICa)], [abs(ICb), 180/pi*angle(ICb)], [abs(ICc), 180/pi*angle(ICc)]<br />

vmi = T * [UA0; UA1forg; UA2forg];<br />

UAa = vmi(1); UAb = vmi(2); UAc = vmi(3);<br />

[abs(UAa), 180/pi*angle(UAa)], [abs(UAb), 180/pi*angle(UAb)], [abs(UAc), 180/pi*angle(UAc)]<br />

vmi = T * [0; IA1forg; IA2forg];<br />

IAa = vmi(1); IAb = vmi(2); IAc = vmi(3);<br />

[abs(IAa), 180/pi*angle(IAa)], [abs(IAb), 180/pi*angle(IAb)], [abs(IAc), 180/pi*angle(IAc)]<br />

vmi = T * [UD0; UD1forg; UD2forg];<br />

UDa = vmi(1); UDb = vmi(2); UDc = vmi(3);<br />

[abs(UDa), 180/pi*angle(UDa)], [abs(UDb), 180/pi*angle(UDb)], [abs(UDc), 180/pi*angle(UDc)]<br />

12

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!