+még1 bonyolultabb
+még1 bonyolultabb
+még1 bonyolultabb
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