curso

Ing. Moiss Ventosilla Zevallos 1

ANALISIS DE FLUJO DE CARGA
METODO DE NEWTON RAPSON

Semana 6, Clase 6

2
 Teorema de Taylor
f ( x) c
f ( x 0 x 0 ) c
f ( x0 ) 2 f ( x0 ) 3 f ( x0 )
f ( x ) f ( x0 ) ( x x0 ) ( x x0 )
2
( x x0 )3 .... c
x x 2
x 3

f ( x0 )
f ( x ) f ( x0 ) ( x x0 ) c
x
f ( x0 )
f ( x0 ) c x
x
1
f ( x0 )
x f ( x0 ) c f ( x0 )(J )1
x
Determinacin de x

x( k 1) x( k ) x(k )

1
x ( k 1) x ( k ) J ( k ) f ( x0 )( k )

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 3

Elemplo
Calcular la raz de la ecuacin f(x)=x2 64 = 0, 0.005
Inicio del proceso
x(0) 5
f (0) (5) 52 64 39
Iteracin 1
f ( x (0) )
J
(0)
2x (0) 10
x
1
x (J (0) )1 f (0) (5) ( 39) 3.9
(0)

10
x(1) x(0) x 0 5 3.9 8.9
f (1) (8.9) 8.92 64 15.21
Iteracin 2
f ( x (1) )
J
(1)
2x (1) 17.8
x
1
x (J (1) )1 f (1) (8.9)
(1)
(15.21) 0.8545
17.8
x(2) x(1) x(1) 8.0455
f (2) (8.0455) 8.04552 64 0.7301
Iteracin 3
f (3) (8.0001) 8.00012 64 0.002
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 4
 f(x)

f(0)

<0.001 f(1)
f(2)
X(2) X
(1)
X(1) X(0) X(0) x

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 5

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 6
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 7
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 8
 Ley de Kirchof
2
q
1

n Ip=Ip1+Ip2++Ipq++Ipn

Vp V1 Vp V2 Vp Vq Vp Vn
Ip ... ...
zp1 zp 2 zpq zpn

V1 V2 V 1 1 1 1 Vn
Ip ... q ... ... ... p
V ...
zp1 zp 2 zpq z z z z zpn
p1 p 2 pq pn
En trminos de admitancia
n
Ip y p1V1 y p 2V2 ... y pqVp ... y pnVn
q 1

En trminos de elementos de una matriz

Ip Yp1V1 Yp2V2 ... YppVp ... YpnVn

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 9
Concepto: Potencia compleja

Potencia compleja constante

entregada a la carga.

I
I
V Carga P&Q
constantes.
S V I

S P jQ VI cos jVI sen Q = P tan

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 10

Concepto: Modelado
Lneas de transmisin
Circuito Pi
Transformadores
Impedancia serie ubicacin de taps
Circuito Pi
Generadores P,V
Potencia activa constante
Potencia reactiva por determinar y dentro de rango
Voltaje constante
Cargas P,Q
Potencia compleja constante

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 11

Formulacin de las ecuaciones de flujo
Considerando todas las barras tenemos la ecuacin matricial

I 1 Y11 Y12 .. Y1 p .. Y1n 1 Y1n V1

I Y Y22 Y2 p Y2 n 1 Y2 n V2
2 21
. .
Y
Ypn V p
p p1
I Yp 2 Ypp Ypn 1
.

n 1 Yn 11 Yn 12
I Yn 1 p Yn 1n 1 Yn 1n Vn 1
I Yn1 Ynn Vn
n Yn 2 Ynp Ynn 1

En forma compacta

[I] = [Y][V]

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 12

 Newton Rapson Cartesiano, barras P-Q, [Y]
n
Pp Qp V p
*
Y
q 1
V
pq q

haciendo
Ypq = Gpq + jBpq
Ip = cp + jdp
Vp = ep + j fp
n n
Ip cp jd p YpqVq (Gpq jBpq )(eq jfq )
q 1 q 1

n
c p epG fpBpp (Gpq eq Bpq fq )
q 1
qp
n
d p epBpp fpBpq (Gpq eq Bpq fq )
q 1
qp
n
Pp ep (Gppep Bppfp ) fp (Bppep Gppfp ) (ep (Gpq eq Bpqfq ) fp (Bpqeq Gpqfq ))
q 1
q np
Qp fp (Gppep Bppfp ) ep (Bppep Gppfp ) (fp (Gpq eq Bpq fq ) ep (Bpq eq Gpq fq ))
q 1
qp
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 13
Newton Rapson Cartesiano, barras P-Q, [Y]
Pp Pp (e1,e2 ... en ,f1,f2 ... fn )

Qp Qp (e1,e2 ... en ,f1,f2 ... fn )

(0) (0) (0)
Pp Pp Pp
Pp Pp (e , e
(0) (0)
... e , f , f
(0) (0) (0)
... f (0)
) e (0)
e (0)
... en(0)
e1 e2 en
1 2 n 1 2 n 1 2

(0) (0) (0)

Pp Pp Pp
f (0)
f (0)
... fn(0)
f1 f2 fn
1 2

(0) (0) (0)

Qp Qp Qp
Qp Qp (e , e
(0) (0)
... e , f , f
(0) (0) (0)
... f(0)
) e (0)
e (0)
... en(0)
e1 e2 en
1 2 n 1 2 n 1 2

(0) (0) (0)

Qp Qp Qp
f(0)
f(0)
... fn(0)
f1 f2 fn
1 2

PP(0) Pp Pp (e1(0),e2(0) ... en(0) ,f1(0) ,f2(0) ... fn(0) )

QP(0) Qp Qp (e1(0), e2(0) ... en(0) ,f1(0) ,f2(0) ... fn(0) )

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 14

NR Cartesiano, barras P-Q, [Y]
P1 P1 P1 P1
e ... ...
en 1 f1 fn 1
1
P1 . ... . . ... . e1
.
Pn 1 ...
Pn 1 Pn 1
...
Pn 1 .
Pn 1 e1 en 1 f1 fn 1 en 1

Q1 Q1 Q1 Q1 Q1 f1
... ...
. e1 en 1 f1 fn 1 .

Qn 1 . ... . . ... . fn 1
Q Qn 1 Qn 1 Qn 1
n 1 ... ...
e1 en 1 f1 fn 1

P J1 J2 e H N e
Q J
J4 f J L f
3

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 15

NR Cartesiano, barras P-Q, [Y]
Formacin de la matriz H Formacin de la matriz N

PP PP pq
H pq epGpq fpBpq pq
Npq epBpq fpGpq
eq fq

PP PP
Hpp epGpq fpBpq cp Npp epBpq fpGpq d p
ep fp
Formacin de la matriz J Formacin de la matriz L

QP QP
J pq fpGpq epBpq pq Lpq fpBpq epGpq pq
fq
eq
QP QP
J pp fpGpp epBpp d p Lpp epGpp fpBpp cp
ep fp
1
ek HN Pk
k k
f JL Q
ek+1=epk+epk Ppk+1 = Pp(esp) Ppk+1
fk+1=fpk+fpk p=1,..,n-1 Qpk+1 = Qp(esp) Qpk+1
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 16
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 17
 1 2

3
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 18
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 19
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 20
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 21
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 22
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 23
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 24
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 25
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 26
CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 27
 Ley de Kirchof
2
q
1

n Ip=Ip1+Ip2++Ipq++Ipn

Vp V1 Vp V2 Vp Vq Vp Vn
Ip ... ...
zp1 zp 2 zpq zpn

V1 V2 V 1 1 1 1 Vn
Ip ... q ... ... ... p
V ...
zp1 zp 2 zpq z z z z zpn
p1 p 2 pq pn
En trminos de admitancia
n
Ip y p1V1 y p 2V2 ... y pqVp ... y pnVn
q 1

En trminos de elementos de una matriz

Ip Yp1V1 Yp2V2 ... YppVp ... YpnVn

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 28
Formulacin de las ecuaciones de flujo
Considerando todas las barras tenemos la ecuacin matricial

I 1 Y11 Y12 .. Y1 p .. Y1n 1 Y1n V1

I Y Y22 Y2 p Y2 n 1 Y2 n V2
2 21
. .
Y
Ypn V p
p p1
I Yp 2 Ypp Ypn 1
.

n 1 Yn 11 Yn 12
I Yn 1 p Yn 1n 1 Yn 1n Vn 1
I Yn1 Ynn Vn
n Yn 2 Ynp Ynn 1

En forma compacta

[I] = [Y][V]

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 29

Newton Rapson Cartesiano, barras P-Q, [Y]
n
Pp Qp V p
*
Y
q 1
V
pq q
n n
Ip cp jd p YpqVq (Gpq jBpq )(eq jfq )
q 1 q 1

haciendo n
c p epG fpBpp (Gpq eq Bpq fq )
Ypq = Gpq + jBpq q 1
qp
Ip = cp + jdp
n
Vp = ep + j fp d p epBpp fpBpq (Gpq eq Bpq fq )
q 1
qp

n
Pp ep (Gppep Bppfp ) fp (Bppep Gppfp ) (ep (Gpq eq Bpqfq ) fp (Bpqeq Gpqfq ))
q 1
q np
Qp fp (Gppep Bppfp ) ep (Bppep Gppfp ) (fp (Gpq eq Bpq fq ) ep (Bpq eq Gpq fq ))
q 1
qp

PP(0) Pp Pp (e1(0),e2(0) ... en(0) ,f1(0) ,f2(0) ... fn(0) ) 1

ek HN Pk
k k
QP(0) Qp Qp (e1(0), e2(0) ... en(0) ,f1(0) ,f2(0) ... fn(0) ) f JL Q

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 30

Newton Rapson Polar

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 31

Newton Rapson Polar

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 32

Newton Rapson Polar
Determinar el Jacobiano para la red

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 33

Newton Rapson Polar-Cartesiano

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 34

Newton Rapson Polar-Cartesiano

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 35

 Elementos fuera de la diagonal

Newton Rapson Polar-Cartesiano

MATRICES N y L
Elementos de la diagonal

Elementos fuera de la diagonal

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 36

Metodo Desacoplado (FDM)

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 37

Metodo Desacoplado (FDM)
MATRICES H y L

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 38

Metodo Desacoplado (FDM)

CURSO: ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA I 39

FIN CLASE 6

40