Scientia Et Technica

ISSN: 0122-1701
[email protected]
Universidad Tecnológica de Pereira
Colombia

GARCÉS RUIZ, ALEJANDRO; GALVIS MANSO, JUAN CARLOS; GALLEGO R., RAMÓN ALFONSO
FLUJO DE CARGA ARMÓNICO PARA SISTEMAS RADIALES
Scientia Et Technica, vol. X, núm. 25, agosto, 2004, pp. 65-70
Universidad Tecnológica de Pereira
Pereira, Colombia

Disponible en: http://www.redalyc.org/articulo.oa?id=84911685013

Cómo citar el artículo

Número completo
Sistema de Información Científica
Más información del artículo Red de Revistas Científicas de América Latina, el Caribe, España y Portugal
Página de la revista en redalyc.org Proyecto académico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
Scientia et Technica Año X, No 25, Agosto. UTP. ISSN 0122-1701 65

FLUJO DE CARGA ARMÓNICO PARA SISTEMAS RADIALES

RESUMEN ALEJANDRO GARCÉS RUIZ
En el estudio de sistemas con cargas no-lineales el flujo de carga debe ser Estudiante X semestre
modificado para incluir la respuesta del sistema a distintas frecuencias armónicas. Ingeniería Eléctrica U.T.P.
En este documento se plantean las principales metodologías de solución al flujo [email protected]
de carga armónico y se presenta un método enfocado hacia los sistemas radiales.
Para esto se utiliza un algoritmo de barrido iterativo. El procedimiento se JUAN CARLOS GALVIS MANSO
presenta para sistemas balanceados pero puede ser expandido fácilmente a Estudiante X semestre
sistemas desbalanceados. Ingeniería eléctrica U.T.P.
[email protected]
PALABRAS CLAVES: Armónicos, flujo de carga armónico, barrido iterativo,
penetración armónica, cargas no-lineales. RAMÓN ALFONSO GALLEGO R
Ingeniero Electricista, Ph.D.
ABSTRACT Profesor Titular
The conventional power flow has been reformulated to permit the inclusion of Facultad de Ingeniería Eléctrica U.T.P.
non-linear loads. This paper presents a review of the different single phase [email protected]
harmonic load flow formulations and shows a methodology for radial distribution
systems. The proposed methodology uses a back / forward sweep method Grupo de Investigación en Planeamiento
applied to harmonic load flow. de Sistemas Eléctricos
Universidad Tecnológica de Pereira
KEYWORDS: Harmonics, harmonic power flow, iterative sweeper, harmonic
penetration, non-linear loads.

transmisión (excluyendo el caso de co-generadores) y su

1. INTRODUCCIÓN modelo consiste en un equivalente de thèvenin en el cual
el voltaje de Thévenin se considera de forma senoidal
El flujo de carga armónico, al igual que el flujo de carga pura y la impedancia de Thèvenin involucra el efecto de
convencional, busca establecer el estado del sistema a la distorsión en el nodo generador.
partir de los parámetros de los elementos lineales que lo
conforman, de información sobre demanda en los nodos, V0
de la potencia generada por las máquinas síncronas, de la
topología del sistema y de las características de los
Nodo
elementos no lineales que dan origen a voltajes y Zth Fuente(0)
corrientes de frecuencias múltiplos de la fundamental en
el sistema. Los elementos lineales y no lineales deben slack
ser modelados considerando la variación que éstos Vth
pueden presentar con la frecuencia.

El documento presenta en su primera parte algunos Figura 1. Equivalente de red del sistema de distribución
modelos para los elementos del sistema, posteriormente
se exponen las principales metodologías que han sido La corriente distorsionada que es inyectada por las cargas
utilizadas para plantear y resolver el problema de flujo de no lineales hace parte de la corriente total que circula por
carga armónico. Finalmente, se propone una solución Zth distorsionando V0. Zth es una impedancia RL serie y
para un sistema de prueba utilizando la técnica de barrido su valor para cada frecuencia armónica esta dado por la
de corrientes, aprovechando la característica radial propia expresión:
de los sistemas de distribución.
(h)
Z th = Rth + j ⋅ hX th (1)
2. MODELAMIENTO
donde h es el orden de la frecuencia armónica
Los elementos del sistema (generadores, considerada; Rth y Xth son parámetros a frecuencia
transformadores, líneas, condensadores, reactores, cargas fundamental y de secuencia positiva.
lineales y no lineales) deben ser modelados considerando
su dependencia con la frecuencia. Para las cargas lineales se considera un equivalente RL
paralelo, donde los valores de R y L son calculados con
En el caso de las redes de distribución los generadores se base a factores asociados al tipo de carga [1] (comercial,
representan por equivalentes de red del sistema de residencial, industrial). Si se tienen datos de potencia se
Fecha de recepción: 31 Mayo de 2004
Fecha de aceptación: 23 Juliode 2004
66 Scientia et Technica Año X, No 25, Agosto 2004. UTP

realizan las conversiones necesarias teniendo en cuenta la para un sistema más extenso de ecuaciones
tensión a frecuencia fundamental. correspondiente al FCA. Esta propuesta requiere un
modelo analítico de las cargas no-lineales que pueda ser
Para transformadores y líneas de distribución el modelo introducido en el jacobiano.
serie RL es adecuado aunque puede considerarse el
efecto capacitivo de las líneas debido a que este varía con El conjunto de ecuaciones correspondiente es el
la frecuencia. siguiente:

Los compensadores, en los sistemas de distribución están En las cargas lineales se conoce normalmente la potencia
básicamente constituidos por bancos de condensadores, a frecuencia fundamental, por lo tanto se plantean las
su modelo es una impedancia constante obtenida a partir ecuaciones de potencia nodal:
( ) ⋅ ∑Y
del valor de la capacitancia o de la potencia reactiva a N
(1) *
Pc − j ⋅ Qc − Vc ⋅ Vk =0
(1) (1) (1) (1)
frecuencia fundamental y el voltaje de barra a frecuencia ck (2)
k =1
fundamental.
El superíndice indica el orden armónico (1 significa
frecuencia fundamental).
El modelamiento de las cargas no lineales depende del
tipo de carga que se considere aunque la regla general es
En las cargas no-lineales se tiene el modelo de los
utilizar una fuente de corriente conformada por las
elementos y el balance de corrientes:
componentes armónicas que la carga inyecta al sistema. N
Los modelos de estas cargas son importantes para la
f (1) (VL(1) ,VL( h ) , β ) − ∑ YLk ⋅Vk(1) = 0 (3)
(1)
convergencia y la exactitud del método. Se recomienda
k =1
usar modelos que consideren la interacción armónica.
En los generadores se tiene las restricciones de potencia y
tensión:
3. FLUJO DE CARGA ARMÓNICO (FCA)
 
PG(1) − Re (VG(1) ) ⋅ ∑ YGk ⋅ Vk(1)  = 0 (4)
N
* (1)
Para una red radial con N nodos, el numero mínimo de
incógnitas que determinan el estado del sistema es
 k =1 
2 ⋅ N PQ + 1 ⋅ N PV donde NPQ es el numero de nodos de VG = constante
(1)
(5)
carga del sistema y NPV es el numero de nodos de En el nodo slack se tiene el valor de tensión (magnitud y
generación. ángulo) a frecuencia fundamental y los valores de
potencia se pueden obtener una vez determinadas las
Cuando una red posee cargas no-lineales el problema se tensiones nodales del sistema. Para frecuencias
hace más complejo, ya que además de las incógnitas a armónicas se requiere hacer cumplir el balance de
frecuencia fundamental aparecen las tensiones y los corrientes:

( )
N
1
⋅ Vs + ∑ Ysk ⋅ Vk = 0
ángulos a frecuencias armónicas de cada nodo del (h) ( h) (h)
sistema para cada una de las frecuencias armónicas (h) Req + h ⋅ j ⋅ X eq k =1
(6)
como se muestra en la tabla 1.
En las cargas lineales se tiene el balance de corrientes y
VARIABLES VARIABLES NUMERO DE la impedancia es calculada con los datos de potencia a
ELEMENTOS frecuencia fundamental:
CONOCIDAS DESCONOCIDAS INCÓGNITAS

( )
Slack V(1) y θ(1) V(h) y θ(h) 2•h 1 N
⋅ Vc + ∑ Yck ⋅ Vk = 0
( h) (h) (h)
Generadores V(1)yP(1) Q(1) y θ(1)
2•NG + 2•NG•h
(NG) generados V(h) y θ(h) Rc + h ⋅ j ⋅ X c k =1
(7)
Cargas (NC) P(1) y Q(1) V y θ
(1) (1)
2•NC + 2•NC•h Rc y Xc son calculados así:
(Vc )
V(h) y θ(h)
Cargas no- Parámetros (β) V(1) y θ(1) (1) 2
2•NL + 2•NL•h
lineales (NL) V(h) y θ(h) Rc + j ⋅ X c = (8)
Total 2 • [ h + ( NG + NC + NL ) • ( 1+h )] Pc(1) − j ⋅ Qc(1)
Tabla 1. Variables del flujo de carga armónico
En las cargas no lineales se tiene el balance de corrientes
de acuerdo con el modelo para cada nivel armónico h:
Existen diversas técnicas para solucionar este problema, N
f ( h ) (VL(1) ,VL( h ) , β ) − ∑ YLk
(h)
la mayoría de estas utilizan algoritmos basados en ⋅ Vk( h ) = 0 (9)
Newton-Raphson tradicional. k =1
3.1 Flujo De Carga Completo El sistema completo de ecuaciones F(X) = 0 puede ser
resuelto para las h frecuencias armónicas consideradas.
El primer enfoque a este problema fue propuesto por Xia
[2] y es una extensión del algoritmo de Newton-Raphson
Scientia et Technica Año X, No 25, Agosto 2004. UTP 67

Este enfoque tiene ventajas en cuanto a convergencia, sin Para agilizar el proceso se puede utilizar una reducción
embargo, el tamaño del jacobiano y la necesidad de del sistema de tal manera que solo considere las cargas
modelos analíticos completos para las cargas no lineales no-lineales: La YBUS depende exclusivamente de las
limitan su uso. Por esta razón existen metodologías variaciones de frecuencia fundamental (a diferencia de
alternas para resolver el problema. las corrientes que dependen tanto de la frecuencia
fundamental como de los armónicos) esto hace que entre
3.2 Penetración Armónica una iteración y otra del sub-problema Fh la matriz YBUS
no cambie. Por esta razón se puede reducir el sistema a
Un enfoque distinto es suponer la división del problema un sistema equivalente que contenga únicamente los
F(X) = 0 en dos sub-problemas (F1,Fh). Estos se nodos con cargas no-lineales, como se muestra en la
resuelven de manera separada bajo las siguientes ecuación (11).
aproximaciones:
• El modelo de las cargas no-lineales depende de las
[YBUS '] = [YNN ] − [YNL ]⋅ [YLL ]−1 ⋅ [YLN ] (11)
tensiones de frecuencia fundamental. Esta metodología puede ser complementada realizando
• La potencia de las cargas lineales esta asociada a la iteraciones adicionales de flujo de carga a frecuencia
frecuencia fundamental. fundamental para considerar la interacción con la
Estas suposiciones implican que el FC puede ser aplicado frecuencia fundamental. El algoritmo completo se
para encontrar los parámetros a frecuencia fundamental. muestra en la figura 1.
Una vez finalizado este paso se tienen definidos los Flujo de Carga
Frecuencia Fundamental
modelos de las cargas lineales y no lineales.

Para encontrar las tensiones armónicas se debe construir

la matriz YBUS para cada frecuencia y resolver la YBUS
ecuación (10):
(h)
I ( h ) = YBUS ⋅ V (h) (10)
Modelos
No-Lineales

En la expresión anterior las inyecciones de corriente son

conocidas porque dependen de la carga. En las cargas
I(h)=YBUS(h)V(h)
no-lineales la corriente es obtenida del modelo mientras
que en las cargas lineales la corriente es cero debido a No
que la impedancia equivalente de estas cargas, (ecuación Tolerancia ?
8), es considerada en la matriz YBUS.
Si
Esta metodología tiene claras ventajas por su simplicidad
pero no considera la interacción armónica, es decir, la Si Cambios
Requiere otro
Modelos a frecuencia
relación entre las corrientes armónicas y los voltajes Flujo fundamental
Fundamental ? ?
armónicos en los modelos de las cargas no-lineales.
Una ventaja adicional es que no se requieren las No

derivadas del modelo para las cargas no-lineales. Esto

Fin
posibilita el modelamiento mediante algoritmos o
funciones no analíticas.
Figura 2. Diagrama de flujo del algoritmo de penetración
armónica iterativa.
3.3 Penetración Armónica Iterativa
4. MÉTODO DE BARRIDO ITERATIVO
Esta metodología tiene la misma estructura básica de la
penetración armónica, y adicionalmente tiene pasos que
Hasta ahora, se han mostrado tres métodos para el flujo
llevan en cuenta la interacción armónica, para esto se
de carga armónico. Entre estos los más destacados son
hace un subproceso iterativo en el cual se actualizan los
los que plantean dos subproblemas: Flujo de carga a
modelos de las cargas no-lineales de acuerdo al cambio
frecuencia fundamental y penetración armónica.
de las tensiones armónicas.
Finalizado el cálculo de penetración armónica, los
En el flujo de carga a frecuencia fundamental pueden ser
voltajes armónicos cambian respecto a sus valores
aplicadas las metodologías tradicionales, sin embargo, en
iniciales y esto produce un cambio en los modelos de las
el caso de redes radiales de distribución, existen
cargas no-lineales haciendo que las corrientes también se
algoritmos más rápidos especialmente diseñados para
modifiquen. El cambio en las corrientes modifica
este tipo de redes, uno de estos métodos es el barrido
nuevamente las tensiones armónicas. Este proceso
iterativo de corrientes.
iterativo se repite hasta que los modelos no cambien de
una iteración a otra.
68 Scientia et Technica Año X, No 25, Agosto 2004. UTP

Este método aplica de manera separada la primera y la la suma de los elementos de cada columna es igual a
segunda ley de Kirchhoff: la primera ley se usa en un cero, esto se expresa en la ecuación (17).
barrido hacia arriba en donde se encuentra los flujos de
corriente por cada línea. La segunda ley se usa en el
[A ]+ [A ]⋅ e
0
T T
N −1 N −1 =0 (17)
barrido hacia abajo encontrando las tensiones en los donde eN-1 es un vector columna cuyos elementos son
nodos. iguales a uno.

En general para cualquier red lineal, las leyes Esta característica permite que la ecuación (17) pueda ser
expresada de la siguiente forma:
[AN −1 ] ⋅ [YB ] ⋅ [AT N −1 ]⋅ (VN −1 − V0 ⋅ e N −1 ) = I N −1
fundamentales de circuitos son expresadas de la siguiente
forma:
Primera ley de Kirchhoff: (18)
I N = [A]⋅ I B La matriz YBUS está definida por AN-1YBATN-1. La
(12)
Segunda ley de Kirchhoff: Voltajes de los elementos en ecuación (18) puede ser separada en dos para resolver el
función de los voltajes de nodo: flujo de carga:
V B = [A] ⋅ V N
T
(13) [AN −1 ] ⋅ I B = I N −1 (19)
Relación voltaje corriente:
I B = [YB ]⋅VB
[YB ] ⋅ [A T
N −1 ]⋅ (V N −1 − V0 ⋅ e N −1 ) = I B (20)
(14)
En donde: La matriz AN-1 es un matriz triangular superior, si el
• YB: Matriz de admitancias de rama. sistema es ordenado. En este caso resolver IB es
• VB: Vector de tensiones por cada rama. equivalente a realizar un barrido hacia arriba mientras
• IB: Vector de corrientes por cada rama. que resolver VN-1 es realizar un barrido hacia abajo.
• VN: Vector de tensiones nodales.
• IN: Vector de inyección de corrientes nodales. El ordenamiento nodal es importante pues permite
• A: Matriz de incidencia de nodos. sistematizar el proceso al calcular el flujo de corriente
por una línea a partir de los flujos en las líneas inferiores
Usualmente, los elementos paralelo son ignorados en el y de la inyección de corriente en el nodo de recibo,
modelamiento de las líneas de distribución, no obstante, (barrido hacia arriba). Igualmente, permite calcular el
cualquier elemento paralelo puede ser reemplazado por voltaje en nodos inferiores a partir del voltaje en un nodo
una inyección de corriente nodal haciendo uso del superior (barrido hacia abajo).
teorema de sustitución. A(1)

Usando las ecuaciones (12), (13) y (14) se obtiene:

[A] ⋅ [YB ] ⋅ [A]T ⋅ VN = I N (15)
B(2) C(3) D(4)
La tensión en el nodo fuente (V0) es conocida, por lo
tanto, si la red se ordena de tal forma que este nodo sea el
primero, se puede modificar la ecuación (15) para
encontrar las tensiones en los nodos siguientes:

 V 
[AN −1 ] ⋅ [YB ] ⋅ [AT 0 ]
AT N −1 ⋅  0  = I N −1
V N −1 
(16) E(5) F(6) G(7) H(8)

Figura 3: esquema de una red ordenada.

La matriz de incidencia de nodos (A) es de tamaño
{NxB} en donde B es el número de ramas del sistema. 4.1 Barrido Hacia Arriba
En el caso de sistemas radiales B = N-1, Por tal motivo,
la matriz reducida (AN-1) es una matriz cuadrada. Las líneas se modelan inicialmente a través de su
La matriz aumentada (A) se construye de tal forma que equivalente serie. Esto permite calcular la corriente por
una rama orientada aportará un valor de 1 a la posición las líneas superiores a partir de las líneas inferiores y la
de la matriz correspondiente al nodo de envío y un valor carga correspondiente.
de –1 al nodo de recibo. Cada columna tendrá solo dos Si se asume un valor inicial para los voltajes de 1 pu se
elementos distintos de cero con valor de ± 1 , por lo tanto pueden calcular las corrientes en todos los nodos del
sistema, según el caso:
Scientia et Technica Año X, No 25, Agosto 2004. UTP 69

Para las cargas lineales

V0 = VTH − I TH ⋅ ZTH (24)
∗
 P + jQ 
I = jθ  (21) Este valor V0 servirá de base para el barrido de voltajes
 Ve  que se explica a continuación.
Para los nodos con impedancia constante (filtros): 4.3 Barrido Hacia Abajo
V
I= Conocida la tensión en el nodo fuente se calcula la
Z (22) tensión en los nodos inferiores. Para esto se hace uso del
ordenamiento nodal.
V = V − I ⋅ (R + j ⋅ X )
Para las cargas no lineales se utiliza la función (25)
correspondiente al modelo de fuente de corriente K K −1 K
para la frecuencia fundamental.
En la ecuación (25) se puede observar que para conocer
Después de haber sido calculadas las corrientes en los el valor de la tensión en un nodo inferior (K) se requiere
nodos se procede a determinar los flujos de corriente por la tensión en el nodo superior, la corriente en la rama de
las líneas, así: alimentación correspondiente y los parámetros del
modelo serie de dicha línea. La corriente por la línea es
I K = I CK + ∑ I i , K una información conocida en el barrido de corrientes
i (23) hacia arriba.

en donde: Si se compara con el método de Newton-Raphson este

IK : Corriente por la línea K método puede converger en un menor número de
ICK : Corriente por los elementos conectados al nodo K. iteraciones. La aplicación del barrido de corrientes es
Ii,K: Corriente por la línea i que esta conectada al nodo K. más rápida porque no requiere de la inversión de matrices
o factorización LDU (la cual esta implícita en la fase de
La figura 4 ilustra el significado de cada una de estas ordenamiento nodal). Esto es especialmente importante
corrientes. en el flujo de carga armónico pues se requiere ejecutar
K múltiples FC para encontrar la solución del problema.
K-1

I1,K Así mismo, los FC que se ejecutan después de la primera

IK iteración no se inicializan en 1pu ya que se conocen los
resultados del flujo anterior. Los FC darán convergencia
R+jX en una o dos iteraciones a partir de la segunda iteración
global.
ICK
Nodo Nodo
I2,K 5. INICIALIZACIÓN DEL FLUJO ARMÓNICO
Envío recibo
DE CARGA
Carga
CK
La selección de las condiciones iniciales en el caso del
Figura 4. Corrientes en un sistema de distribución. FCA es mucho más compleja que las del flujo
convencional. No obstante, en [3] se hacen algunas
Si se desea aplicar el modelo completo de las líneas recomendaciones para seleccionar estos valores:
(modelo PI) se adiciona el efecto capacitivo a través • La amplitud y fase de las tensiones a frecuencia
admitancias en paralelo a cada lado de la línea. fundamental se seleccionan de la manera habitual en
1 pu a 0º.
4.2 Actualización Del Nodo Slack • La magnitud de los voltajes armónicos en los
elementos no-lineales y a las frecuencias de estudio
Cuando las corrientes en las N-1 líneas han sido se asumen en un valor cercano a 0.01 pu y 0º.
actualizadas se debe corregir el valor del voltaje en la La segunda recomendación no debe ser tomada de
subestación (Nodo fuente). La figura 1 muestra el manera estricta pues depende del tipo de cargas no-
modelo utilizado. lineales que se consideren, además, las magnitudes de las
tensiones tienden a disminuir a medida que la frecuencia
La corriente equivalente en el nodo fuente Ith es la suma aumenta.
de las corrientes en todos los elementos del sistema, su
valor actualizado se tiene una vez finalizado el barrido No todas las tensiones armónicas deben ser inicializadas,
hacia arriba. Para conocer el nuevo valor de la tensión ya que los elementos no-lineales solo producen
V0 se utiliza la ecuación (24): contaminación armónica a determinadas frecuencias.
70 Scientia et Technica Año X, No 25, Agosto 2004. UTP

En el caso de aplicar el método de Newton Raphson Tabla 2. Resultados en pu del flujo de carga armónico.
completo, la inicialización de los voltajes armónicos es
fundamental pues de no hacerlo se pueden presentar Potencia base: 1 MW
singularidades indeseables en la matriz jacobiana. Voltaje base: 13.8 kV en IND1
Estos resultados fueron comprobados usando un software
Si se pretende aplicar penetración armónica, la de simulación en el dominio del tiempo (Simulink).
inicialización de los voltajes armónicos no tiene ninguna Tiempo de cálculo del flujo de carga armónico: 0.7 s.
repercusión en el resultado final del proceso (a diferencia Tiempo de cálculo de Simulink: 109 s.
de la penetración armónica iterativa). Esto se debe a que
el método no considera la interacción armónica. 7. CONCLUSIONES

Una fase de inicialización es propuesta en [4], esta utiliza La metodología de barrido iterativo puede ser aplicada al
el concepto de penetración armónica tradicional flujo de carga armónico acelerando el tiempo de cálculo
empleando voltajes planos (1 pu) a frecuencia y produciendo resultados comparables con los obtenidos
fundamental para obtener un aproximado de las con metodologías más exactas como las que utilizan el
corrientes armónicas (usando el modelo porcentual de las dominio del tiempo.
cargas no lineales). De esta forma se tienen valores
iniciales de cada una de las tensiones armónicas del La metodología de barrido hacia atrás permite obviar las
sistema. Este tipo de inicialización será adoptada en el expresiones analíticas para las cargas no-lineales al no
presente trabajo, con el método de penetración armónica tener que construir un jacobiano.
iterativa.
El modelamiento de los elementos no lineales
6. EJEMPLO DE APLICACIÓN considerando la interacción armónica permite utilizar el
método de penetración armónica iterativa.
La figura 5 muestra un sistema radial con una carga no-
lineal (convertidor de 6 pulsos). El convertidor de 6 8. BIBLIOGRAFÍA
pulsos es modelado con la misma filosofía expuesta en
[5]. Este modelamiento permite considerar la interacción [1] XU, Wilson. Harmonic modeling and simulation
armónica necesaria para utilizar penetración armónica task force reference documents. Disponible en internet:
iterativa. http://www.ee.ualberta.ca/pwrsys/IEEE/download.html

Convertidor
[2] XIA, Daozhi and HEYDT G.T. Harmonic power
de 6 pulsos flow studies part I – Formulation and solution. En:
IEEE Transactions on power apparatus and systems. Vol
PAS 101 Nº6 (Jun 1982).
Figura 5. Ejemplo de una red radial.
[3] ________. Harmonic power flow studies part II
Los datos del sistema son: Implementation and practical application. En: IEEE
• Línea: Z = 0.02 + j 0.06 [Ω] transactions on power apparatus and systems. VOL PAS
• Transformador: Z = (0.0098+j0.0762) [pu] 101 Nº6 (Jun 1982).
20[MVA].
• Nodo 1: V = 69 kV Barraje infinito. [4] SAINZ SAPERA, Luis. Estudio de la formulación y
• Nodo 2 (IND1): V = 13.8 kV resolución del problema del flujo de armónico de cargas.
Carga: 10 MW fp = 0.85 ↓ Barcelona. 1995. Tesis Doctoral. Universidad
• Nodo 3(IND2): V = 13.8 kV politécnica de Cataluña. Departamento de ingeniería
Carga: 3MW, fp: 0.9↓ eléctrica.
Convertidor de 6 pulsos sin filtro capacitivo.
P = 5 MW [5] BOIX ARAGONÈS, Oriol. Estudio y modelización
Los resultados del flujo de carga armónico usando en régimen permanente de cargas no lineales para el
penetración armónica iterativa se muestran en la tabla 2 análisis armónico de redes eléctricas. Barcelona. 1996.
276p. Tesis doctoral. Universidad politécnica de
Cataluña. Departamento de ingeniería eléctrica.
VARIABLE FRECUENCIA ARMÓNICA
1 5 7 11 13 17
V(IND1) 0.96 0.02 0.01 0.01 0 0
V(IND2) 0.95 0.02 0.01 0.02 0.01 0.01
I(Convertidor) 5.22 1.11 0.61 0.36 0.37 0.11