Aps - Chap9 (Power Flow)
Aps - Chap9 (Power Flow)
steady-state
balanced single-phase network
network may contain hundreds of nodes and
branches with impedance X specified in per unit on
MVA base
Nodal solution
1
1
=
Z ij rij + jxij
see Fig.6.1
Bus admittance
equations
the admittance is
based on bus-tobus: see Fig.6.2
if no connection
between bus-tobus, leave as zero
node voltage
equation is in the
form
I bus = YbusVbus
Node-voltage matrix
Ibus=YbusVbus
Y11 Y12 Y1i Y1n
V1
I1
Node-voltage matrix
Yii = yij
ji
j =0
Yij = Y ji = yij
when the bus currents are known, bus voltages are unknown,
bus voltage can be solved as
1
Vbus = Ybus
I bus
Z bus = Ybus
Node-voltage matrix
Ybus
0
0 j12.5
j12.5
j12.5
0
0
Gauss-Seidal
Newton-Raphson
Quasi-Newton
Gause-Seidal method
GAUSE-SEIDAL METHOD
Nature of Gause-Seidal method
GAUSE-SEIDAL METHOD
GAUSE-SEIDAL METHOD
f1 ( x1 , x2 , , xn ) = c1
x1 = c1 + g1 ( x1 , x2 , , xn )
f 2 ( x1 , x2 , , xn ) = c2
x2 = c2 + g 2 ( x1 , x2 , , xn )
f n ( x1 , x2 , , xn ) = cn
xn = cn + g n ( x1 , x2 , , xn )
j =0
j =1
I i = Vi yij yijV j
ji
Pi jQi
Ii =
Vi *
ji
ji
Pi sch jQisch
(k )
+
y
V
ij k
Vi *( k )
=
yij
n
n
= Re Vi *( k ) Vi ( k ) yij yijV j( k )
j =0
j =1
ji
n
n
*( k ) ( k )
(k )
V
yij yijV j
= Im Vi
i
j =0
j =1
j i
ji
Vi
( k +1)
Vi *( k )
+ yijV j( k )
yij
ji
n
n
= Im Vi *( k ) Vi ( k ) yij yijV j( k )
j =0
j =1
j i
ji
Pi sch jQi( k )
+ yijVk( k )
*( k )
Vi
=
yij
ji
solve Vi
stopping criteria
Vi
solve Pi
n
n
*( k ) ( k )
(k )
V
yij yijV j
= Re Vi
i
j =0
j =1
j i
ji
n
n
*( k ) ( k )
(k )
V
yij yijV j
= Im Vi
i
j =0
j =1
ji
Pi ( k +1)
solve Qi
Qi( k +1)
( 0)
(x )
( 0) 2
+ = c
f (x
( 0)
df
+ x ) = f ( x ) + x ( 0 ) = c
dx
( 0)
( 0)
c
(
1
)
(
0
)
the new approximation of x x = x +
( 0)
df
dx
( 0)
c ( 0 )
df
dx
( 0)
(k )
c ( k )
(k )
df
dx
x ( k +1) = x ( k ) + x ( k )
for more information, see Ex.(6.4)
Newtons method converges faster than Gauss-Seidal, the root
may converge to a root different from the expected one or
diverge if the starting value is not close enough to the root
( 0)
( 0)
f1
f1
f1
( 0)
( 0)
( 0)
( 0)
( 0)
xn ( 0 ) = c1
x2 + +
f1 ( x + x ) = f1 ( x ) + x1 +
x1
x2
xn
( 0)
( 0)
( 0)
( 0)
( 0)
( 0)
f 2
f 2
f 2
( 0)
( 0)
( 0)
( 0)
( 0)
xn ( 0 ) = c2
x2 + +
f 2 ( x + x ) = f 2 ( x ) + x1 +
x1
x2
xn
f
f
f
( 0)
( 0)
( 0)
f n ( x ( 0 ) + x ( 0 ) ) = f n ( x ( 0 ) ) + n x1 + n x2 + + n xn = cn
x1
x2
xn
f 2
( 0)
c
f
2 2 =
x1
( 0)
cn f n
( 0)
f n
x
1
f1
x2
( 0)
f 2
x2
( 0)
f n
x2
( 0)
f1
xn ( 0 )
x1
( 0)
f 2 ( 0 )
x2
xn
( 0)
xn
( 0)
f n
xn
( 0)
C ( k )
x1( k )
(k )
x
= 2
(k )
xn
c1 f1( k )
(k )
c
f
= 2 2
(k )
cn f n
J (k )
f ( k )
1
x1
(k )
f 2
= x1
(k )
n
x
1
f1
x2
(k )
f 2
x2
(k )
f n
(k )
(k )
f1
xn
(k )
f 2
xn
(k )
f n
xn
S ji = V j I *ji
j =1
V j Yij cos( ij i + j )
j =1
J 4 V
P J 1
Q = J
3
V
j i
Pi
= Vi V j Yij sin ( ij i + j )
j
Energy Conversion Lab
ji
J 4 V
P J 1
Q = J
3
Pi
= 2 Vi Yii cos ii + V j Yij cos(ij i + j )
Vi
j i
Pi
= Vi Yij cos( ij i + j )
Vj
ji
V
j i
V j Yij cos( ij i + j )
Qi
= Vi V j Yij cos( ij i + j )
j
Energy Conversion Lab
ji
J1
J
3
J 2
J 4 V
Qi
= 2 Vi Yii sin ii V j Yij sin ( ij i + j )
Vi
j i
Qi
= Vi Yij sin ( ij i + j )
Vj
ji
j=
n1
V j Yij cos( ij i + j )
Pi ( k ) = Pi sch Pi ( k ) ,
V
j =1
V j Yij cos( ij i + j )
Pi ( k ) = Pi sch Pi ( k )
Energy Conversion Lab
J 2
J 4 V
J1
J
3
Qi( k )
P
Q =
J1
J
3
J 2
J 4 V
P
Q =
J1
0
J 4 V
Q
= B '' V
Vi
P
,
V
V = [B"]
Q
V