PROGRAMA NACIONAL DE FORMACIÓN
EN ELECTRICIDAD
(PNFE)
UNIDAD CURRICULAR: SISTEMAS ELÉCTRICOS DE POTENCIA
T4–ELEC–1M
MODELOS DE CÁLCULOS PARA FLUJO DE POTENCIA
DOCENTE RESPONSABLE: PARTICIPANTES:
Ing. Pedro Rosal TSU. Luis Leal
C.I.: V-28.031.365
TSU. Jesús Ambrosino
C.I.: V-27.940.201
Ciudad Bolívar, 17 de Abril del 2023
1. EXPLIQUE Y REALICE LA DEFINICIÓN DEL CONCEPTO DE FLUJO DE
POTENCIA
El flujo de potencia es un análisis numérico del flujo de energía eléctrica en un sistema
interconectado, donde se considera que ese sistema está en operación normal de estado
estable. Los estudios de flujo de potencia o flujo de carga son importantes para planificar la
expansión futura de los sistemas de potencia, así como para determinar la mejor operación
de los sistemas existentes. La principal información que se obtiene del estudio de flujo de
potencia es la magnitud y el ángulo de fase del voltaje en cada barra, y la potencia real y
reactiva que fluye en cada línea. Una vez que se conocen los voltajes de todos los nodos, se
pueden calcular los flujos de línea y las pérdidas. El punto de partida para resolver
problemas de flujo de potencia es identificar las variables conocidas y desconocidas en el
sistema. Con base en estas variables, los nodos se clasifican en tres tipos:
- Barras de referencia: también conocida como barra Slack, en esta barra se conoce la
magnitud de voltaje y el ángulo de fase, y los valores de potencia activa y reactiva neta se
desconocen.
- Barras P-Q: también conocidas como barras de carga, en este tipo de barra se conocen los
valores de potencia activa y reactiva neta, pero los valores de voltaje y ángulos son
desconocidos.
- Barras P-V: a estas barras también se las conoce como barras de generación en las cuales
los valores de potencia activa neta y magnitud de voltaje son conocidos, mientras los
valores de potencia reactiva y el ángulo no son especificados.
2. REPRESENTACIÓN DE LOS ELEMENTOS DE SEP EN LOS ESTUDIOS DE
FLUJO DE CARGA (FC)
3. CORRIENTES QUE LLEGAN A UN NODO EN FUNCIÓN DE LAS
POTENCIAS
Suponga un nodo k que posee n barras conectadas y donde cada una de ellas posee una
corriente orientada saliendo de la barra k, y una corriente única entrando a la barra k que
n
puede ser escrita por I K =V K Y KK + ∑ Y kj V j
j=l
j ≠k
La potencia que es inyectada a la barra k viene dada por:
¿
S K =V K I K
S¿ K =V ¿ K I K
Donde
[ ]
n
S =V V k Y kk + ∑ Y kj V j =Pk − j Qk
¿ ¿
k k
j=1
J≠k
{[ ]}
n
Pk =ℜ V ¿k V k Y kk + ∑ Y kj V j
j=1
J ≠k
{[ ]}
n
Qk =−ℑ V k V k Y kk +∑ Y kj V j
¿
j=1
J≠k
Se puede ver fácilmente que el voltaje de la barra k, puede ser escrito como:
[ ]
n
1 P k − jQk
Vk= ¿ +∑ Y kj V j
Y kk Vk j=1
J≠k
4. CORRIENTES QUE LLEGAN A LOS NODOS EN FUNCIÓN DE LOS
VOLTAJES EN LOS NODOS Y DE LAS ADMITANCIAS DE LAS LÍNEAS
Como se sabe de la teoría de redes eléctricas, el desempeño de una red puede
determinarse mediante un conjunto de R voltajes linealmente independientes o N corrientes
linealmente independientes, donde R y N son el rango y la nulidad asociados con dicha red.
Para que se cumpla la condición de independencia lineal los R voltajes deben ser los
correspondientes a los arcos de enlace, es decir, los arcos que no forman parte dl árbol.
Una forma de modelar la red es a través de un conjunto de n voltajes, que bien pueden
escogerse como los voltajes entre cada uno de los nodos que no son el de referencia, y la
referencia misma, ya que también estos forman un arbolo. Otra forma es por medio de un
conjunto de n-1 corrientes de enlace que bien podrían ser las corrientes que pasan por los
elementos de la estructura radial. En lo sucesivo, a estos voltajes uy corrientes se les
llamaran voltajes nodales y corrientes de rama.
Los voltajes que provocan las corrientes de rama se pueden expresar como una
combinación n lineal de los voltajes nodales, ya que estos últimos forman un conjunto de
voltajes linealmente independientes. Así pues, la diferencia de potencial de cualquier rama
se expresa como:
V П k =V П −V k
k k
k ∈ { 2, 3 , … , n }
Donde es evidente queV П k es el voltaje de la rama que conecta al nodo k con su predecesor
k
Πk.
La ecuación anterior expresada bien puede utilizarse en forma recursiva para el cálculo de
voltajes nodales, como en el método de barrido de Shirmohammadi. La extensión para toda
la red es:
V П k =C V k
k
Donde:
VΠkk: Vector de voltajes de rama de orden (n − 1) × 1.
C: Matriz de conexión con posibles entradas {0, 1, −1} y de orden (n − 1) × n.
VK: Vector de voltajes nodales de orden n × 1.
Por otra parte, es posible expresar las corrientes nodales como una combinación lineal de
las corrientes de rama, ya que estas últimas forman un conjunto de corrientes linealmente
independientes.
pk
I k =−J k + ∑ J σ k
x
x=1
k ∈ {1 , 2 , … , n }
Donde p es el número de sucesores del nodo k.
La ecuación bien puede utilizarse en forma recursiva para el cálculo de corrientes nodales,
como en el método de barrido de Shirmohammadi. Ahora, para toda la red se tiene:
I k =B J k
Donde:
IK: Vector de corrientes nodales de orden n × 1
B: Matriz de conexión, con posibles entradas {0, 1, −1} y de orden n × (n − 1).
Jk: Vector de corrientes de Rama de orden (n − 1) × 1.
Ahora, de la ley de Ohm, se tiene que las corrientes de rama son
Jk = yk V П k
k
k ∈ { 2, 3 , … , n }
Donde yk es la admitancia de la rama que conecta al nodo k y su predecesor Πk, y en forma
matricial para todas las ramas
J k =Y k V П k
k
Ahora, combinando las tres fórmulas anteriormente expresadas se obtiene una ecuación que
expresa las corrientes nodales como función de los voltajes nodales:
I k =B Y k C V k
De lo anterior, como Ik y VK son los vectores de corriente y voltaje nodal respectivamente,
es de esperar que
BY k C=Y
Haciendo una sustitución en la ecuación, se tiene
I k =Y V k
O simplemente,
I =Y V
Eliminando el sub índice k por simplicidad, donde Y es la matriz de admitancia nodal o
Ybus, como se le conoce normalmente en la literatura, en la cual cada elemento de la
diagonal es la suma de las admitancias incidentes en cada nodo, y los elementos fuera de la
diagonal; el negativo de la admitancia entre nodos físicamente conectados.
Donde y•k es la admitancia del nodo k al nodo de referencia, comúnmente llamada
admitancia en derivación, siendo normalmente de naturaleza capacitiva, ya que se debe, en
la mayoría de los casos, a la capacitancia de la l nea o a equipos de compensación reactiva
o a ambos. Es común que en una red de distribución se desprecie y•k, si es debida solo a la
capacitancia de la línea, ya que su valor es peque o comparado con las admitancias de rama.
Lo anterior provoca que el rango de la matriz de admitancia nodal disminuya en una
unidad, es decir, Rank (Y) = n − 1, lo cual es evidente al ver que la suma de los elementos
de cualquier la o columna es cero.
Finalmente, se puede representar en forma explícita como:
O desacoplando
5. FLUJO DE POTENCIA ENTRE DOS NODOS
Para un enlace (i-k) entre dos nodos cualesquiera de un SEP que tiene una impedancia de
enlace Zik, y unos voltajes de nodo Vi∠δi y Vk∠δk
V i−V k
La potencia del nodo “i” es: Si=Vi∗I iK =V i =Pi − jQ i pu
Z ik
V k −V i
La potencia del nodo “k” es: Sk=Vk∗I iK =V k =P k − j Q k pu
Z ik
Las pérdidas de potencia activa y reactiva se calculan sumando las expresiones anteriores:
∆ Sik =Si + S k =∆ Pik +∆ Q ik
6. EXPLICACIÓN DETALLADA DE LA FORMULACIÓN COMPLETA DE LA
MATRIZ DE ADMITANCIA.
Conceptos Básicos.
Algoritmos convencionales para la sistematización de los estudios de corto circuito en
sistemas eléctricos de potencia trifásicos exigen la formación de la matriz impedancia de
nodos, la cual se caracteriza por la alta densidad de elementos diferentes de cero. En
contraste, el número de éstos y el esfuerzo computacional necesario para obtener su
inversa, la matriz admitancia de nodos, es mucho menor. Puesto que la matriz admitancia
de nodos [Y] es diagonalmente dominante y definida positiva se puede despreciar el error
acumulado por redondeo de las operaciones en punto flotante. De otra parte, el esfuerzo
computacional requerido para modificar la matriz admitancia de nodos es
incomparablemente menor que el necesario para reflejar los cambios en la estructura
topológica de la red en la matriz impedancia de nodos. Además, para minimizar el número
de elementos de relleno existen técnicas de reasignación de números a los nodos.
Obtención de la Matriz Admitancia de Barras
Se multiplica la ecuación I ramabx 1= y ramabxb U ramabx 1por AT: AT
AT Irama = AT yrama Uramay
Teniendo en cuenta las expresiones (Ibarra4x1 = AT4x8 Irama8x1) y (Urama8x1 = A8x4
Ubarra4x1) resulta:
Ibarra = (AT yrama A) Ubarra
Ibarra4x1 = Ybarra4x4 Ubarra4x1
Dónde: Ybarra4x1 = AT4x8 yrama8x8 A8x4 y quedaría de la siguiente manera:
La Matriz Admitancia de Barra en un SEP de n barras
La ecuación matricial compleja encontrada en Ibarra 4x1 = Ybarra4x4 Ubarra4x1 se puede
generalizar para un SEP de n barras:
Donde n es el número de barras del SEP, sin considerar el nodo tierra (nodo 0). En forma
compacta se puede escribir: Ibarranx1 = YbarranxnUbarranx1
Es decir, el vector de corrientes de barra (corrientes inyectadas netas en la barra) es igual al
producto de la matriz admitancia de barra y el vector de tensiones de barra (tensiones de
fase), donde:
Ybarranxn = ATnxb yramabxb Abxn
Resultando que el elemento ii de la matriz admitancia de barra bale: Y ii =∑
KεΩ
Y iK
Ω = conjunto de los nodos vecinos del nodo i (incluido el nodo tierra). El elemento fuera
de la diagonal:Y ij =Y ji =−Y ij
La matriz admitancia de barra Ybarra es cuadrada de orden nxn, es simétrica respecto de
la diagonal principal y es esparza para los SEP reales (aproximadamente un 90-95 % de sus
elementos son iguales a cero). Si todos los elementos fueran diferentes de cero, esto
significaría que cada nodo del SEP está conectado físicamente con todos los otros nodos.
Para almacenar los datos diferentes de cero en la memoria del computador, se hace uso de
técnicas llamadas de almacenamiento compacto. La ecuación matricial:
Se deduce que la corriente inyectada en el nodo genérico i bale:
n
Ii = Yi1 U1 + Yi2 U2 +..........+ Yii Ui +.........+ Yij Uj +......... + Yin Un = ∑ Y iK U k
K−I
Y de aquí resulta:
a) Yii = Ii/Ui con Uk = 0 para todo k = 1, 2,…, n y k≠i
El elemento de admitancia de barra en la posición ii (diagonal principal), Y ii, es igual a
la corriente inyectada en el nodo i, Ii, cuando en el mismo se tiene aplicada una tensión de 1
pu, y el resto de los nodos está cortocircuitado a tierra (Uk = 0).
b) Yij = Ii/Uj con Uk = 0 para todo k = 1, 2,…, i,…, n y k≠j
El elemento de admitancia de barra en la posición ij (fuera de la diagonal), Yij, es igual
a la corriente inyectada en el nodo i, Ii, cuando se tiene aplicada una tensión de 1pu en el
nodo j y el resto de los nodos está cortocircuitado a tierra (Uk = 0).
A partir del SEP de 3 nodos que se muestra a continuación, determinar la matriz de
admitancia.
Los datos correspondientes a las líneas son los siguientes:
LINEA RESISTENCIA REACTANCIA IMPEDANCIA
L12 0.030 0.3 0.030 +j0.3
L23 0.060 0.1998 0.060 + j0.1998
Como primer paso se debe de realizar la conversión de las impedancias a admitancias por
medio de la inversa. La siguiente Tabla resume los valores obtenidos.
LINEA IMPEDANCIA ADMITANCIA
L12 0.030 +j0.3 0.33 - j3.3
L23 0.060 + j0.1998 1.38 - j4.59
Por lo que ya se cuenta con los elementos para formar la matriz de admitancia Y BUS ∈
C3x3. A partir del diagrama del sistema eléctrico de potencia se forma la matriz de
admitancia, tomando en cuenta las admitancias que existen entre nodos, donde cada fila
pertenecerá a la relación de un nodo con respecto a todos los demás, por ejemplo, la hilera
que posee un color de fuente amarillo indicará los valores de admitancia del nodo 1, el L12
sería el nodo 1 con 1 ya que para este serían todas las admitancias que lo toquen, y
solamente existe el contacto a través de la línea 1-2, por lo tanto será el valor positivo de la
admitancia de la línea 1-2. Ahora, el nodo 1 respecto al nodo 2 será el valor negativo de la
admitancia que tenga en común -L12. Y debido a que el nodo 1 no toca al nodo 3, o no
tienen nada en común, el valor es 0. Y así sucesivamente con todos los demás nodos y sus
relaciones con las otras barras.
[ ]
L12 −L12 0
3 ×3
Y BUS ∈C = −L12 L12 +L23 −L23
0 −L23 L23
Y al reemplazar los valores de la Tabla que contiene los datos de admitancia, la matriz Y BUS
resultante es:
[ ]
0.33− j3.3 −0.33+ j3.3 0
3 ×3
Y BUS ∈C = −0.33+ j3.3 1.71− j7.89 −1.38+ j 4.59
0 −1.38+ j 4.59 1.38− j 4.59
Se obtiene la matriz de admitancia simétrica.
Se muestran los valores de la conductancia eléctrica G y la susceptancia B.
A partir del SEP de 4 nodos que se muestra a continuación, determinar la matriz de
admitancia.
Los datos correspondientes a las líneas son los siguientes:
LINEA RESISTENCIA REACTANCIA IMPEDANCIA
L12 0.01008 0.05040 0.01008 +j0.05040
L23 0.00744 0.03720 0.00744 + j0.03720
L24 0.00744 0.03720 0.00744 + j0.03720
L34 0.01272 0.06360 0.01272 + j0.06360
Como primer paso se debe de realizar la conversión de las impedancias a admitancias por
medio de la inversa. La siguiente Tabla resume los valores obtenidos.
LINEA IMPEDANCIA ADMITANCIA
L12 0.01008 +j0.05040 3.815629 - j19.078144
L23 0.00744 + j0.03720 5.169561 - j25.847809
L24 0.00744 + j0.03720 5.169561 - j25.847809
L34 0.01272 + j0.06360 3.023705 - j15.118528
Por lo que ya se cuenta con los elementos para formar la matriz de admitancia YBUS ∈ C4x4
A partir del diagrama del sistema eléctrico de potencia se forma la matriz de admitancia
[ ]
L 12+ L13 −L 12 −L 13 0
Y BUS ∈C 4×4
= −L12 L12+ L 24 0 −L 24
−L13 0 L13+ L 34 −L 34
0 −L 24 −L 34 L24 + L 34
Y al reemplazar los valores de la Tabla que contiene los datos de admitancia, la matriz
YBUS resultante es:
[ ]
8.985190− j 44.835953 −3.815629+ j 19.078144 −5.169561+ j 25.847809 0
Y= −3.815629+ j 19.078144 8.985190− j 44.835953 0 −5.169561+ j 25.847809
−5.169561+ j25.847809 0 8.193267− j 40.863838 −3.023705+ j 15.118528
0 −5.169561+ j25.847809 −3.023705+ j15.118528 8.193267− j 40.863838
Cada columna reflejará el orden que se va siguiendo de los nodos con respecto a todos
los demás, por ejemplo, en una matriz Y BUS ∈ C3x3, la primera fila corresponde al nodo 1,
significa que el primer valor pertenece a la relación que existe entre el nodo 1 con 1, el
segundo valor al nodo 1 con 2, y el tercer nodo 1 con 3, razón por lo cual las siguientes dos
filas deberán seguir el mismo orden, formando sus columnas de acuerdo a la misma
secuencia.
[ ]
1 con 1 1 con 2 1 con 3
3 ×3
Y BUS ∈C = 2 con 1 2 con 2 2 con 3
3 con 1 3 con 2 3 con 3
● Número de nodo por fila.
● Orden por columna.
7. RESOLUCIÓN MATEMÁTICA DE UN PROBLEMA DE UN SISTEMA DE
POTENCIA PARA CALCULAR EL FC POR LOS DOS MÉTODOS (GAUSS
SEYDEL Y NEWTON RAPHSON)
Gauss-Seidel
A partir del SEP encuentre la solución de flujos de potencia por medio del método Gauss-
Seidel.
Los datos son los siguientes:
NODO VOLTAJE P. P. G. Q.
GENERAD CONSUMO GENERADA CONSUMO
A
1 1∠0 - 0 - 0
2 v2 0 0.6 0 0.3
PASO 1: A partir de la información que se proporciona obtener la matriz de admitancia
YBUS.
Y 12 =Z12−1=0− j 4
−1
Y 1=Z 1 =0− jO .116
Por lo que la matriz YBUS ∈ C2x2 en forma rectangular es igual a
Y BUS =G+ jB= |0−0+j 4.166
j4
0+ j 4
0− j 4 |
Donde al separar los términos se obtiene
G∈R
2 ×2
=|00 00|
Y
jB ∈C
2 ×2
|− j+4.166
=
j4
j4
−j4 |
Separar los términos facilita la formulación de las ecuaciones de flujo de potencia.
NODO TIPO VARIABLES VARIABLES
CONOCIDAS DESCONOCIDAS
1 Generación (slack) P, v1 Q, θ1
2 Carga P, Q v2, θ2
Las variables que se plantean corresponden a las del nodo de carga, v2 y θ2.
Se plantea el método de Gauss-Seidel.
n
Pi− jQi=v
¿
i ∑ Y ij v j
j=1
Donde i representa el número de nodo y n es la cantidad de nodos que existen en el SEP.
ITERACION 1: Para realizar la iteración por medio de Gauss-Seidel se expande la
ecuación anterior tomando en cuenta el nodo de las variables y la cantidad de nodos. Para i
= 2 y n = 2 se tiene
(0)
Pi− jQ2=v 2 ¿
Donde al despejar para el voltaje v2 se tiene la expresión
v 2=
[
1 P2− jQ2
Y 22 v2
(0)
−Y 12 v 1
]
Siendo v2 el resultado en magnitud y ángulo.
Al considerar inicio de perfil plano y los datos del sistema se tiene
v 2(0)=1 ∠0 °
De la tabla de datos se obtiene
P2=PG −PC =0−0.6=−0.6
Q2=QG −QC =0−0.3=−0.3
v1 =1∠0 °
De la matriz de admitancias se obtienen los siguientes valores:
Y 12 = j 4
Y 22=− j 4
Estos datos se reemplazan en la ecuación que se obtuvo anteriormente.
Por lo que el valor del voltaje v2 es
v 2=
− j4 [
1 −0.6+ jO .3
1∠ 0°
− j 4 (1 ∠ 0 °)
]
Y el resultado es
v 2=0.9370 ∠−9.21°
Y se obtiene la información de magnitud y ángulo del nodo correspondiente.
ITERACIÓN 2: Se repite el procedimiento reemplazando el valor del perfil plano por el
resultado anterior, por lo que la ecuación ahora es
v 2=
1
[ −0.6+ jO .3
− j 4 0.9370∠−9.21 °
− j 4 (1 ∠0 ° )
]
Y el resultado es
v 2=0.9619 ∠−10.22 °
Y se obtiene la información de magnitud y ángulo del nodo correspondiente.
ITERACIONES POSTERIORES: Se repite el procedimiento hasta obtener una variación
pequeña en los resultados. Aquí se muestran solamente las primeras 4 iteraciones del
proceso de Gauss-Seidel.
ITERACIÓN 3:
v 2=
1
[ −0.6+ jO .3
− j 4 0.9619∠−10.22 ° ]
− j 4 ( 1 ∠ 0° ) =0.9655 ∠−9.97 °
ITERACIÓN 4:
v 2=
1
[ −0.6 + jO .3
− j 4 0.9655∠−9.97 ° ]
− j 4 ( 1∠ 0 ° ) =0.9648 ∠−9.93°
Después de 4 iteraciones el resultado es el siguiente
v 2=0.9648 ∠−9.93°
Y se han encontrado el valor de magnitud y ángulo del nodo de carga.
Newton-Raphson.
A partir del SEP encuentre la solución de flujos de potencia por medio del método Newton-
Raphson.
Los datos son los siguientes:
NODO VOLTAJE P. P. G. Q.
GENERAD CONSUMO GENERADA CONSUMO
A
1 1∠0 - 0 - 0
2 v2 0 0.6 0 0.3
PASO 1: A partir de la información que se proporciona obtener la matriz de admitancia
YBUS.
−1
Y 12 =Z12 =0− j 4
−1
Y 1=Z 1 =0− jO .116
Por lo que la matriz YBUS ∈ C2x2 en forma rectangular es igual a
|
Y BUS =G+ jB= 0+ j 4.166 0− j 4
0− j 4 0+ j 4 |
Donde al separar los términos se obtiene
G ∈ R2 ×2 = 0 0
0 0 | |
Y
jB ∈C
2 ×2
= | j−4.166
j4
−j4
j4|
Separar los términos facilita la formulación de las ecuaciones de flujo de potencia en forma
rectangular.
NODO TIPO VARIABLES VARIABLES
CONOCIDAS DESCONOCIDAS
1 Generación (slack) P, v1 Q, θ1
2 Carga P, Q v2, θ2
Las variables que se plantean corresponden a las del nodo de carga, v2 y θ2. Debido a esto,
se deben plantear dos ecuaciones para la solución de las dos variables, que se formulan a
partir de las expresiones de flujos de potencia en forma rectangular.
A partir de las ecuaciones en forma rectangular
n
Pi=v i ∑ v j ( Gij cos θ ij + Bij sinθ ij )
j=1
n
Qi=v i ∑ v j ( Gij sin θij −Bij cos θij )
j=1
Donde i representa el número de nodo y n es la cantidad de nodos que existen en el SEP.
De acuerdo al Angulo del nodo θ2 se plantea
2
(P ¿ ¿ G 2−P D 2)=v 2 ∑ v j ( G2 j cos θ 2 j +B 2 j sin θ 2 j ) ¿
j=1
Y al expandir se obtiene
( P ¿ ¿ G 2−P D 2)=v 2 [v 1 ( G21 cos θ21 + B21 sin θ21 ) + v 2 ( G22 cos θ22+ B22 sin θ22) ]¿
Donde se sabe que el ángulo θ 22 = 0 y que la matriz G es nula
(P ¿ ¿ G 2−P D 2)=v 2 ¿¿
Y al reducir términos se obtiene
(P ¿ ¿ G 2−P D 2)=v 1 v 2 B21 sinθ 2 ¿
Remplazando datos de la Tabla y matriz YBUS se obtiene:
(0−0.6)=( 1 ) v 2 ( 4) sinθ 2
Simplificando se tiene
4 v 2 sinθ 2+ 0.6=0
De acuerdo con el voltaje del nodo v2 se plantea
2
(Q ¿ ¿ G 2−Q D 2)=v 2 ∑ v j ( G2 j cos θ 2 j +B 2 j sinθ 2 j ) ¿
j=1
Y al expandir se obtiene
(Q ¿ ¿ G 2−Q D 2)=v 2 [v 1 ( G 21 sin θ21−B21 cos θ 21) + v 2 ( G 22 sin θ22−B22 cos θ 22) ]¿
Donde se sabe que el ángulo θ 22 = 0 y que la matriz G es nula
(Q ¿ ¿ G 2−Q D 2)=v 2 ¿ ¿
Y al reducir términos se obtiene
2
(Q ¿ ¿ G 2−Q D 2)=−v 1 v 2 B21 cos θ2 +4 v 2 ¿
Remplazando datos de la Tabla y matriz YBUS se obtiene:
2
(0−0.3)=4 v 2−(1)v 2(4)cos θ2
Simplificando se tiene
2
4 v 2−4 v 2 cos θ 2+0.3=0
Las ecuaciones para determinar las variables v2 y θ 2 son
g1 (θ2 , v 2 )=4 v 2 sinθ 2+ 0.6=0
g2 (θ2 , v 2 )=4 v 22−4 v 2 cos θ2 +0.3=0
En este punto comienza la formulación del método de Newton-Raphson debido a que ya se
cuenta con las ecuaciones no lineales.
La solución se obtiene a partir de la expresión (8) obtenida en la formulación teórica
| || |
(0) (0)
(0) ∆ x1 ∆ g1
J = (0)
= (0)
∆ x2 ∆ g2
(0 )
Donde ∆ x 1,2 representa el valor del ángulo θ 2 en radianes y el voltaje v2 respectivamente.
PASO 2: Se establecen los valores iniciales de las variables desconocidas
(0)
θ2 =0 rad
(0)
V 2 =1.0 pu
Se debe plantear una tolerancia, que en este caso es ∈=10−3.
PASO 3: Se debe reemplazar en la ecuación (8) de la explicación anterior.
| || |
(0 ) (0 )
(0) ∆ θ2 ∆ g1
J = (0 )
= (0 )
∆ v2 ∆ g2
Que representa la formulación para determinar los incrementos.
Iteración 1: Se determina la matriz Jacobiana J 0 R2 ×2 a partir de las derivadas parciales de
g1 y g2, que se expresa como
| ||
0
∂ g1 ∂ g1
|
θ =0 rad , v =1.0
0 ∂ θ1 ∂ v2 4 v cos θ2 4 sinθ 2 2 2
J= = 2
∂ g2 ∂ g2 4 v 2 sin θ2 8 v 2−4 cos θ2
∂ θ2 ∂ v2
Y al sustituir los valores iniciales establecidos, se obtiene
(0)
J = |44cos 0
sin 0
4 sin 0
=
8−4 cos 0 0
4
| | 04|
El valor numérico de la matriz Jacobiana.
Por otro lado, la función de error se expresa por medio de
[ ][ ][ ]
(0) (0 ) ( 0) ( 0) (0 )
∆ g1 0−g 1(θ 2 , v 2 ,u) b 1−h1 (θ2 , v 2 , u)
(0)
= (0 ) ( 0) = ( 0) (0 )
∆ g2 0−g 2(θ 2 , v 2 ,u) b 2−h2 (θ2 , v 2 , u)
Donde al realizar la evaluación con los valores de inicio se obtiene el vector en forma
numérica
[ ]|
∆ g(0)
1
(0)
∆ g2
=
−0.6−4 sin 0
−0.3−4−4 cos 0
=
−0.6
−0.3 || |
Por lo tanto, la expresión (8) del método se representa por medio
| || || |
( 0)
4 0 ∆ θ2 −0.6
= =
0 4 ∆ v2
( 0)
−0.3
Y al realizar el despeje para el vector de los incrementos
| ||
1
| |
0
|
(0 )
∆ θ2 4 −0.6
=
∆v
(0 )
1 −0.3
2 0
4
Se obtiene la solución numérica.
PASO 4: Al resolver para los incrementos se obtiene
(0)
∆ θ 2 =0.15 rad
∆ v(0)
2 =−0.075
PASO 5: Ahora se obtiene el valor de la siguiente iteración
(1) (0) (0)
θ2 =θ2 +∆ θ 2 =0+ (−0.15 ) =−0.15 rad
(1) (0) (0 )
v 2 =v 2 + ∆ v 2 =1+ (−0.075 ) =0.925
PASO 6: Debido a que no cumple aun con el valor de se debe realizar una segunda
iteración, por lo que los valores encontrados deben de sustituirse nuevamente en la
ecuación (8), que expresa el método iterativo.
Iteración 2: Se sustituyen los nuevos valores en J
(1 )
J =
[ 4(0.925) cos−0.15 4 sin−0.15
=
3.658 −0.598
4( 0.925)sin−0.15 8(0.925)−4 cos −0.15 −0.553 3.445 ]| |
Y en el vector de error
[ ]| || |
(1 )
∆ g1 −0.6−4 (0.925)sin −0.15 −0.0471
= =
∆g
(1 )
2
−0.3−4(0.925) 2−4( 0.925) cos−0.15 −0.0640
EL vector de incrementos para la segunda iteración es
|∆ v | | | |−0.0471
−0.0640|
(1 ) −1
∆ θ2 3.658 −0.598
=
(1 )
2
−0.553 3.445
Y al resolver se obtiene
∆ θ(1)
2 =−0.0163rad
∆ v(1)
2 =−0.075
Donde se observa que el incremento es menor al de la iteración anterior.
Por último en la segunda iteración, se suma el incremento al valor anterior
(2) (1) (1)
θ2 =θ2 +∆ θ2 =−0.15+ (−0.0163 )=−0.1663rad
v(2) (1) (1)
2 =v 2 + ∆ v 2 =0.925+ (−0.0212 ) =0.9038
Y se obtienen los valores de inicio para la siguiente iteración.
Iteración 3: Como en el paso anterior se sustituye el nuevo valor en la matriz Jacobiana
(2 )
J =
[ 4( 0.9038)cos−0.1663 4 sin−0.1663
=
3.565 −0.662
4(0.9038) sin−0.1663 8(0.9038)−4 cos −0.1663 −0.598 3.286 ]| |
Y en la función de error
[ ]| || |
(1 ) −3
∆ g1 −0.6−4 (0.9038) sin−0.1663 −1.55 × 10
= =
(1 )
∆ g2 −0.3−4(0.9038) 2−4(0.9038) cos−0.1663 −2.093× 10−3
Se despeja para el vector de incrementos
| ||
∆ θ(22 )
|| |
−1 −3
3.565 −0.662 −1.55 × 10
=
∆ v2 (2 )
−0.598 3.286 −2.093× 10
−3
Y se obtiene la solución
(2) −3
∆ θ 2 =−0.5724 ×10 rad
∆ v(2)
2 =−0.7411× 10
−3
Debido a que los incrementos están dentro del rango de ∈ se puede detener en este punto
las iteraciones.
Los siguientes valores para las variables son
θ2 =θ 2 +∆ θ 2 =−0.0163+ (−0.5724 ×10 )=−0.1668 rad
(3) (2) (2) −3
2 =v 2 + ∆ v 2 =0.9038+ (−0.7411 ×10 )=0.9031
v(3) (2) (2) −3
Y debido a que el valor de error está dentro del rango, serían los valores que satisfacen
hasta este punto los requerimientos del ejercicio.