Matriz Admitancia
Matriz Admitancia
Matriz Admitancia
• Conceptos Básicos.
Los 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 Y11 Y12 Y13 Y14 V1 I 1
incomparablemente menor que el necesario para reflejar Y
los cambios en la estructura topológica de la red en la
21 Y22 Y23 Y24 V2 I 2
matriz impedancia de nodos.
Además, para minimizar el número de elementos de
Y31 Y32 Y33
Y34 V3 I3
relleno existen técnicas de reasignación de números a los
nodos. Y41 Y42 Y43 Y44 V4 I 4
Matriz Admitancia
• Factorización Triangular.
El método de eliminación de Gauss para resolver un conjunto linealmente independiente de
ecuaciones lineales se puede visualizar como una descomposición de la matriz de
coeficientes [A] en tres factores, a saber: una matriz de permutación [P] asociada con la
estrategia de pivote, una matriz triangular inferior unitaria [L] que contiene esencialmente
los multiplicadores usados y almacenados en las posiciones que se vuelven cero durante el
proceso y la matriz triangular superior final [U].
Matemáticamente:
[A] =[ P ] [L] [U] L11
L L22
L
21
i 1
L31 L32 L33
u ij a pi j lir u rj i≤j
r 1 L41 L42 L43 L44
y teniendo en cuenta que lii = 1 1 U 12 U 13 U 14
j 1 1 U U
a pi j lir urj U 23 24
lij r 1
i≤j 1 U 34
U ij
1
Matriz Admitancia
• Formación de la Matriz Admitancia de Nodos.
i pq y pq y m v pq
i y y rs v rs
rr m
Puesto que inicialmente se considera que se trata de dos elementos sin acoplamiento mutuo entre
ellos de admitancias propias ypq e yrs únicamente queda pendiente el efecto de los términos y m en
dichas ecuaciones, las cuales se deben reajustar como se describe a continuación:
i p ... y m V r V s ... 0
i q ... y m Vr V s ... 0 Y pr Y pr y m
i r ... y m V p V q ... 0 Y ps Y ps y m
i ... y m V p Vq ... 0
s
Yqr Yqr y m
Como regla general, cada admitancia mutua Yqs Yqs y m
ym se debe sumar a los elementos de la matriz
admitancia ubicados en la fila y columna que Yrp Yrp y m
corresponden a los nodos de envío, (p, r) y (r, Yrq Yrq y m
p), y a los de recibo, (q, s) y (s, q), y se resta
en las localizaciones definidas por el nodo de Y ps Ysp y m
envío de un elemento y el de recibo del otro, Ysq Ysq y m
(p, s), (s, p), (q, r) y (r, q).
Matriz Admitancia
• Modificación de la Matriz Admitancia.
• Eliminación de nodos (Reducción de Kron)
En ciertas ocasiones es necesario solo conocer los voltajes de ciertas barras del sistema, por ejemplo
si se tratan de subestaciones de importancia. Se puede aplicar la eliminación gaussiana para reducir
las ecuaciones de Ybarra del sistema a un conjunto que contenga sólo las barras deseadas.
En el proceso de eliminación gaussiana, se quita secuencialmente del sistema original de N
ecuaciones con N incógnitas, una variable de voltaje de barra en cada etapa.
Por ejemplo, se pueden escribir las ecuaciones de admitancias de nodo cuando I1 = O en un
sistema de cuatro barras, en la forma
Los sistemas de potencia a gran escala tienen solamente un pequeño número de líneas transmisión
que se conectan a cada subestación de gran capacidad. La relación del número ramas al de nodos,
en la gráfica de la red de tales sistemas, es 1.5 aproximadamente y Ybarra correspondiente tiene,
principalmente, elementos que valen cero. En efecto, si hay 750 ramas en una red de 500 nodos
(sin incluir el nodo de referencia), el número total de elementos que no son cero es de
(500+2x750) = 2 000, porque cada nodo tiene un elemento diagonal asociado y cada rama da
origen a dos elementos colocados simétricamente fuera de la diagonal.
Desde el punto de vista de la velocidad computacional, la exactitud y el almacenamiento de
información, es deseable procesar sólo los elementos que no son cero en Ybarra y evitar rellenar
nuevos elementos diferentes cero en el transcurso de la eliminación gaussiana y de la factorización
triangular.