Aportaciones Al Estudio Del Flujo Armónico de Cargas
Aportaciones Al Estudio Del Flujo Armónico de Cargas
Aportaciones Al Estudio Del Flujo Armónico de Cargas
TESIS DOCTORAL
Aportaciones al estudio del flujo armnico de cargas
TESIS DOCTORAL
A mi familia
ndice General
1 Introduccin
1.1
1.2
Objetivos de la tesis................................................................................................. 3
1.2.1 Calidad en el suministro elctrico ............................................................... 4
1.2.2 Modelizacin de dispositivos no lineales .................................................... 5
1.2.3 Formulacin del flujo armnico de cargas .................................................. 7
1.3
Estructura de la tesis................................................................................................ 8
11
2.2
49
125
4.1
Introduccin......................................................................................................... 126
4.2
4.3
4.4
169
5.1
5.2
173
191
2.1
2.2
2.3
iv
ndice de Figuras
2.1. Aplicacin de una funcin ventana a una seal ......................................................... 16
2.2. Muestreo de seales para su procesamiento .............................................................. 17
2.3. Representacin del STDFT de una seal g(t) ............................................................ 18
2.4. Intensidad consumida por un rectificador AC/DC de seis pulsos ideal ..................... 21
2.5. Representacin del clculo del STHD ....................................................................... 25
2.6. Variacin tpica del STHD en funcin de T para una seal aperidica................... 26
2.7. Alimentacin de un horno elctrico de arco .............................................................. 34
2.8. STHD de un horno elctrico de arco en funcin del ancho de ventana T ............... 35
2.9. Prdidas de potencia en funcin del STHD ............................................................... 36
2.10. Ilustracin del aliasing ............................................................................................... 38
2.11. a) Funciones de orden-0, g1(t), y orden-1, g2(t); b) Funciones cannica-0, g3(t),
y cannica-1, g4(t). ..................................................................................................... 40
2.12. Curvas de error de la funcin de orden-0 g1(t) = etU(t)........................................... 40
2.13. Intensidad consumida por una fuente de alimentacin monofsica, iM ..................... 41
2.14. Contenido armnico de iM.......................................................................................... 41
2.15. Curvas de error para la intensidad iM . ........................................................................ 42
2.16. Intensidad consumida por una fuente de alimentacin trifsica, iT. .......................... 43
2.17. Contenido armnico de iT. ......................................................................................... 43
2.18. Curvas de error para la intensidad iT.......................................................................... 44
2.19. Intensidad consumida por una lmpara de descarga, iL. ............................................ 45
2.20. Contenido armnico de iL. ......................................................................................... 45
2.21. Curvas de error para la intensidad iL.......................................................................... 45
2.22. Intensidad consumida por un rectificador AC/DC de seis pulsos ideal, iR6............... 46
2.23. Intensidad consumida por un rectificador AC/DC de doce pulsos ideal, iR12............ 46
2.24. Contenido armnico de iR6......................................................................................... 47
2.25. Contenido armnico de iR12. ...................................................................................... 47
2.26. Curva de error para las intensidades iR6 e iR12............................................................ 47
3.1. Esquema del rectificador monofsico con filtro capacitivo........................................51
3.2
viii
ndice de Tablas
2.1. Lmites de la distorsin de la tensin segn la normativa IEEE................................ 22
2.2. STDFT de una seal de clase A ................................................................................. 31
2.3. STDFT de una seal de clase I................................................................................... 32
2.4. Clculo del K-factor y del valor eficaz usando funciones ventana............................ 33
2.5. K-factor y valor eficaz para diferentes valores de ancho de ventana......................... 34
2.6. STHD del horno elctrico de arco para intervalos de tiempo consecutivos .............. 35
3.1. Lmites armnicos para equipos de Clase D. ............................................................. 94
3.2. Lmites armnicos para equipos de Clase A ............................................................ 121
3.3. Lmites armnicos para equipos de corriente superior a 16 A (Estadio 1). ............. 122
3.4. Lmites armnicos para equipos de corriente superior a 16 A (Estadio 2) .............. 122
4.1. Datos e incgnitas del flujo de cargas convencional ............................................... 127
4.2. Datos e incgnitas del flujo armnico de cargas...................................................... 129
4.3. Datos e incgnitas de la penetracin armnica ........................................................ 131
4.4. Datos e incgnitas de la penetracin armnica iterativa.......................................... 135
ix
4.5. Datos e incgnitas del flujo de cargas armnico simplificado ................................ 140
4.6. Datos e incgnitas del CHLFF. ................................................................................ 144
4.7. Datos e incgnitas del CHLFH................................................................................. 146
4.8. Datos e incgnitas del ejemplo para la penetracin armnica ................................. 151
4.9. Datos e incgnitas del ejemplo para la penetracin armnica iterativa................... 153
4.10. Datos e incgnitas del ejemplo para el flujo de cargas armnico simplificado....... 155
4.11. Datos e incgnitas del ejemplo para el flujo armnico de cargas completo
( componente fundametnal) ......................................................................................159
4.12. Datos e incgnitas del ejemplo para el flujo armnico de completo
(componentes fundamental y armnicas) .................................................................160
4.13. Datos e incgnitas del ejemplo numrico estudiado ................................................ 164
Captulo 1
Introduccin
Introduccin
Introduccin
Introduccin
Introduccin
Captulo 2
Calidad en el suministro elctrico
Las aportaciones de este captulo se centran, por una parte, en el estudio y la
caracterizacin de un nuevo ndice de calidad para seales aperidicas, el cual llamaremos
Distorsin Armnica de Ventana Mvil (STHD, Short-Term Harmonic Distortion) [13,
14], y por otra, en el estudio de la cuantificacin del error que se comete en la
determinacin de los ndices de calidad para seales peridicas, al aplicar la transformada
discreta de Fourier para el clculo de las tensiones armnicas del sistema [15, 16].
Tal como se menciona en la introduccin de la presente tesis, la problemtica que
implica el aumento de la distorsin de las tensiones y corrientes en los circuitos elctricos,
debido al incremento del nmero de dispositivos no lineales, ha derivado, entre otras cosas,
en la definicin de ndices y normativas cuyo objetivo es cuantificar y limitar la distorsin
de las ondas elctricas a niveles que la red y los receptores puedan tolerar. La aplicacin de
las normativas anteriores supone el anlisis previo del espectro frecuencial de las seales
11
del sistema, tensiones y corrientes, con el fin de determinar su distorsin. En este sentido,
si la seal es peridica, la transformada de Fourier (FT, Fourier Transform) para seales
continuas, la transformada discreta de Fourier (DFT, Discrete Fourier Transform) para
seales discretas o la transformada rpida de Fourier (FFT, Fast Fourier Transform) para
seales discretas con un nmero de muestras potencia de dos, seran las herramientas a
utilizar [7, 9, 70]. En el caso de seales aperidicas, existen herramientas como la
transformada de Fourier de ventana mvil (STFT, Short-Time Fourier Transform) para
seales continuas, la transformada discreta de Fourier de ventana mvil (STDFT, ShortTime Discrete Fourier Transform) o el anlisis mediante la transformada wavelet que
permiten analizar el espectro frecuencial de la seal [20, 23, 24].
Histricamente, para el caso de seales peridicas han existido numerosos ndices de
calidad [1, 6, 17] y las principales recomendaciones existentes para seales peridicas se
encuentran en las normas de la Comisin Electrotcnica Internacional (CEI), las normas de
la CENELEC (Comit Europeo de Normalizacin Electrotcnica), las recomendaciones
fijadas por el IEEE (Institute of Electrical and Electronics Engineers) y las
recomendaciones del CIGRE (Congress Internationale des Grand Rseaux lectriques
Haute Tension) [1]. La mayora de ndices de calidad y, en consecuencia, de normativas,
estn basadas en el desarrollo de Fourier de la seal, por lo que dejan de ser tiles en el
momento de trabajar con seales aperidicas [19]. As, para el caso de seales aperidicas
no existen tantos estudios realizados y el nmero de ndices de calidad es reducido [6].
Desgraciadamente, existen muchas cargas que no consumen corrientes peridicas. Entre
ellas se encuentran los hornos de arco elctricos, los soldadores elctricos, laminadoras de
acero, y ciertos variadores de velocidad para esquemas de control especializados. Existen
otros casos en los cuales la corriente de carga es no peridica: la carga puede variar, puede
estar bajo controles que causen una variacin en ella, o la tensin que produce puede estar
sufriendo variaciones no peridicas. Debido a que pocos procesos son realmente
estacionarios, las variaciones en las corrientes de carga y la aparicin de aperiocidades son
sucesos corrientes. Por este motivo, es deseable extender la definicin de ndices de
calidad existentes para contemplar dicha aperiodicidad.
12
Frente a todo lo anterior en este captulo se realizan dos estudios. Por un lado, el estudio
de un nuevo ndice de calidad para seales aperidicas, y por otro, el anlisis del error
cometido al utilizar la transformada de Fourier en el clculo de la distorsin armnica de
las seales peridicas en un sistema elctrico de potencia.
1
2
G ( ) e j t d
( )
(2.1)
( )
(2.2)
1
Gk e jk1t
T k =
(2.3)
13
(2.4)
1
2
G ( ) e
j n
(2.5)
ge
n =
j n
(2.6)
1
N
N 1
Gk e
jk
2
Ts n
T
( n = 0 N 1)
(2.7)
k =0
Gk = g n e
n=0
14
jk
2
Ts n
T
( k = 0 N 1)
(2.8)
Gk = DFT { g n } = g n e
jk
2
Ts n
T
( k = 0 N 1)
(2.9)
n =0
Se observa que los trminos Gk as obtenidos corresponden a una frecuencia fk = k/T = kf1,
donde f1 = 1/T es la resolucin en el espectro frecuencial, o frecuencia mnima, del clculo
realizado.
El clculo de la DFT comporta el tratamiento de seales discretas de duracin bastante
larga y en tales casos es muy importante emplear procedimientos computacionales
eficientes. En este sentido si el nmero de muestras de la seal discreta tiene un nmero de
muestras potencia de dos se emplea la transformada rpida de Fourier, FFT, [7, 70].
En el caso de que la seal de tiempo continuo, g(t), presente un contenido armnico
variable con el tiempo, se recurre a la transformada de Fourier de ventana mvil, STFT, si
se trabaja en tiempo continuo, o a la transformada discreta de Fourier de ventana mvil,
STDFT, si se trabaja en tiempo discreto a partir del muestreo de la seal g(t).
15
g(t)
SEAL
t
1.0
w(t t)
VENTANA
g(t)w(t t)
FRACCIN
DE SEAL
( )
(2.10)
m-sima se podra expresar como gn,m = g((m1)T + nTs)w(nTs) con n = 0, ..., N-1,
m = 0, 1, ..., siendo Ts el tiempo de muestreo y T = NTs la duracin de la ventana.
n=0
n=N1
....
Gk , m = STDFT { g n ,m } = g n ,m e
jk
2
Ts n
T
( k = 0 N 1)
(2.11)
n=0
Al igual que suceda en la expresin (2.9) los trminos Gk,m corresponden a una frecuencia
fk = k/T = kfmin, donde fmin = 1/T sera la resolucin en el espectro frecuencial, o
frecuencia mnima, del clculo realizado.
Tal como se ha comentado al caracterizar la STFT, la funcin ventana utilizada es la
funcin ventana rectangular de amplitud unitaria, es decir,
1 n = 0, , N 1
wn = w ( nTs ) =
0 caso contrario
(2.12)
V1
Vm
Ventana m-sima
(N-1) fmin
.....
.....
Valor Gk,m
fmin
k fmin
Frecuencia k
.....
0
g(t)
N muestras
.....
.....
mT
As,
- Para una ventana m se obtendran Gk,m valores, con k = 0 N-1, y la interseccin
marcada en el dibujo captara la componente Gk,m que corresponde a la frecuencia
fk = kfmin y a un instante de tiempo tm = mt.
- Los resultados de la DFT, es decir los valores Gk,m, sern distintos para cada
ventana.
18
f min T
1
4
(2.13)
19
seal por lo que no son aplicables en todos los casos y, en muchas ocasiones, estos ndices
pueden ocultar informacin en lugar de mostrarla.
Para tensiones y corrientes peridicas, se han desarrollado un amplio conjunto de
ndices de calidad de potencia, los ms comunes estn definidos para seales que se
pueden descomponer segn series de Fourier, y su definicin se expresa a partir de los
trminos fundamental y armnicos de la seal. En referencia a los ndices de seales
peridicas, ya han sido histricamente estudiados, [1, 6, 17], sobre la base del desarrollo de
Fourier de la seal, y slo se har un breve resumen.
En referencia a los ndices para seales no peridicas, para las cuales no existe el
desarrollo de Fourier, aparecen en la bibliografa estudios al respecto [6, 19, 21, 25]
aunque el conjunto de indicadores no est tan definido ni desarrollado como en el caso
peridico.
2.1.2.1 Seales peridicas
Tal y como se ha visto, toda onda distorsionada que sea peridica, de periodo T, puede
descomponerse en suma de ondas sinusoidales de frecuencia mltiplo de la frecuencia
fundamental f1 = 1/T segn (2.3) a partir de (2.4) si la seal es continua o segn (2.7) a
partir de (2.8) si la seal es discreta. Se denomina componente fundamental a la
componente de frecuencia f1, y componentes armnicas (o armnicos) a las de frecuencias
que son mltiplo de f1. As, los armnicos es la forma matemtica de describir, en rgimen
permanente, la distorsin de las ondas de tensin y corriente que permite abordar el estudio
de los sistemas de potencia en presencia de cargas no lineales. Por ejemplo, para la forma
de onda mostrada en la Figura 2.4, su expresin en funcin de sus componentes
fundamental y armnicas es,
i (t ) =
donde 1 = 2f1.
20
= 1; k = 6n + 1, n = 0,1, 2,
2 3
I D k cos ( k1t ) k
k = k
k = 1; k = 6n 1, n = 1, 2,
(2.14)
i(t)
ID
2/3
/3
4/3
2
5/3
Figura 2.4. Intensidad consumida por un rectificador AC/DC de seis pulsos ideal
THD =
HDh =
(V )
h
h=2
V1
Vh
V1
(2.15)
(2.16)
La CEI llama a la expresin (2.15) factor de distorsin (DF, Distortion Factor). Ambos
ndices presentan la ventaja de que son fciles y rpidos de calcular. Sin embargo, en su
definicin se puede perder informacin del espectro total de la seal ya que no describe de
forma unvoca la seal.
Un valor comnmente utilizado como frontera entre alta y baja distorsin es el
correspondiente a un THD igual a 5%. Sin embargo, este valor es demasiado alto para ser
21
Tensin
Distorsin de tensin
individual (%)
Distorsin de la tensin
total THD (%)
Hasta 69 kV
3.0
5.0
69.001 kV a 161 kV
1.5
2.5
Ms de 161 kV
1.0
1.5
Adems del THD existen otros ndices, como el factor de influencia telefnica (TIF,
Telephone Influence Factor), el cual es una variacin del THD. En la definicin de este
ndice, los armnicos tienen un peso en funcin de la respuesta del odo humano y de la
manera en que afectan a sistemas de comunicacin adyacentes (tensiones inducidas), y en
el denominador aparece el valor eficaz de la tensin,
TIF =
(w V )
k 2
k =1
(V )
(2.17)
k 2
k =1
(w V )
k =1
k 2
(2.18)
donde wk son los pesos relativos de los armnicos y son utilizados a fin de obtener una
indicacin aproximada de la interferencia causada por cada armnico.
22
2.2.1 Caracterizacin
La primera idea es utilizar la transformada STDFT de un modo similar a como se utiliza
la DFT para seales peridicas al definir el THD, expresin (2.15). De este modo, el nuevo
ndice contendr informacin de las diferentes componentes frecuenciales presentes en la
seal, no slo de los componentes armnicos, y podr ser usado tanto para caracterizar
seales peridicas como aperidicas.
Dada una seal, se define el ndice Distorsin Armnica de Ventana Mvil (STHD,
Short Term Harmonic Distortion), para cada una de las m ventanas para las que se ha
calculado la STDFT segn (2.11), como
N 1
STDFT
STHDm =
'
k =1
k k1
2
k ,m
STDFTk1 ,m
(2.19)
donde,
- N es el nmero de frecuencias para las que se ha calculado la STDFT.
- k1 = int(f1/fmin) es la posicin de la frecuencia fundamental de la seal, f1, en el
espectro de frecuencias calculado.
24
V1
Vm
(N-1)fmin
.....
.....
kfmin
Gk,m STDFTk,m
.....
.....
fmin
Frecuencia k
.....
.....
0
STHD0
.....
STHDm
.....
1/(fmin)lim
STHD
STHD se incrementa
con el ancho de la
ventana
Ancho de ventana, T
Figura 2.6. Variacin tpica del STHD en funcin de T para una seal aperidica
26
El teorema de Parseval establece la relacin entre la energa de una seal g(t) y el valor de la magnitud de
su transformada de Fourier.
27
DF
1 + THDI
(2.20)
donde THDI hace referencia a la distorsin armnica total de la corriente consumida i(t).
Para el caso aperidico, esta relacin no est definida, pero de forma similar a la obtencin
de la expresin (2.19), se puede obtener una expresin aproximada igual a la anterior en
funcin del STHD,
PFm
DF
1 + STHDIi ,m
(2.21)
2 f1 dt RMS
K factor =
i 2
(2.22)
RMS
donde [g]RMS indica el valor eficaz (RMS, Root Mean Square) de la seal g(t).
Dada una intensidad i(t) peridica pero no sinusoidal, la expresin anterior resulta,
28
K factor =
( hI )
h =1
(I )
h =1
(2.23)
TIFm =
( w STDFT )
k =1
(2.24)
STDFT
k =1
V Tm =
k ,m
2
k ,m
( w STDFT )
k =1
k ,m
(2.25)
29
(2.26)
30
Para esta seal se calcul la STFDT, que se muestra en la Tabla 2.2, siguiendo el patrn
de la Figura 2.5. Las filas de la Tabla 2.2 corresponden a diferentes frecuencias espaciadas
con una frecuencia fmin = 60 Hz. Las columnas corresponden a ventanas espaciadas un
intervalo de tiempo T = 16.67 ms. El valor de STHD es evaluado en cada ventana y
mostrado en la ltima fila de la Tabla. Las filas correspondientes a las frecuencias
comprendidas entre 180 y 980 Hz, y entre 1140 y 1440 Hz, se han eliminado por motivos
de claridad.
Se observa como los resultados se aproximan bastante al valor terico que se debera
haber obtenido, es decir 0.2, pues en este caso los ndices STHD y THD deberan coincidir.
f (Hz)
V1
V2
V3
V4
V5
V6
1500
0.226
0.227
0.226
0.226
0.227
0.226
1080
0.095
0.076
0.095
0.095
0.076
0.095
1020
16.881
16.888
16.881
16.888
16.881
16.888
960
0.487
0.491
0.487
0.487
0.491
0.487
120
0.222
0.224
0.222
0.222
0.224
0.222
60
83.700
83.701
83.700
83.700
83.701
83.700
0.400
0.401
0.400
0.400
0.401
0.400
STHDm
0.2027
0.2028
0.2027
0.2027
0.2028
0.2027
(2.27)
31
f (Hz)
V1
V2
V3
V4
V5
V6
1500
0.467
0.234
0.356
0.406
0.471
0.095
840
3.230
2.891
3.047
3.113
3.240
2.787
780
15.06
14.71
14.87
14.94
15.07
14.611
720
5.051
5.410
5.252
5.184
5.038
5.501
120
0.224
0.619
0.459
0.487
0.096
0.705
60
83.312
84.121
83.958
83.038
83.530
84.208
0.023
0.827
0.665
0.250
0.239
0.912
STHDm
0.201
0.197
0.198
0.201
0.200
0.196
32
(2.28)
En las Tablas 2.4 y 2.5 se muestra el clculo del K-factor y del valor eficaz de las
seales anteriores utilizando funciones ventana rectangulares. En la primera tabla, se
utiliz una ventana de ancho igual a 1/60 segundos y se compar el resultado con el
obtenido considerando una duracin de la muestra de un segundo. Los resultados no varan
apreciablemente entre ventanas, pero su valor resulta inferior al correspondiente a un
segundo.
De la Tabla 2.5 se desprende que a medida que aumenta el ancho de ventana
considerado, los valores del K-factor y del valor eficaz van tendiendo a los valores
correspondientes a los considerados en el caso de una duracin de un segundo. Esto es una
consecuencia del hecho de que una ventana de duracin mayor se traduce en una
resolucin en el dominio de la frecuencia ms fina, es decir, el espectro frecuencial de la
seal calculado es ms exacto.
Caso
Seal g1(t)
Seal g2(t)
K-factor
[g1]RMS
K-factor
[g2]RMS
t=1s
4.4400
0.7745
6.842
0.7211
V1
4.0412
0.7738
6.661
0.7201
V2
4.0416
0.7738
6.953
0.7209
V3
4.0420
0.7738
6.602
0.7204
V4
4.0424
0.7738
6.643
0.7202
V5
4.0428
0.7738
6.712
0.7202
V6
4.0433
0.7738
6.937
0.7201
Tabla 2.4. Clculo del K-factor y del valor eficaz usando funciones ventana
33
Seal g1(t)
Ancho de
ventana
(s)
Seal g2(t)
K-factor
RMS
K-factor
RMS
1/60
4.04
0.7738
6.661
0.7201
2/60
4.26
0.7738
6.815
0.7205
3/60
4.33
0.7738
6.798
0.7205
4/60
4.34
0.7744
6.780
0.7209
5/60
4.36
0.7743
6.814
0.7208
6/60
4.38
0.7742
6.840
0.7208
Tabla 2.5. K-factor y valor eficaz para diferentes valores de ancho de ventana (seales g1 y g2 )
100 MVA
ABIERTO
ALIMENT. 138 kV
PT
V
34.5 kV
34.5 kV
CT
I
50 MVA AC
HORNO DE ARCO
SVC
34
0,4
0,35
0,3
STHD
0,25
0,2
0,15
0,1
0,05
0
0
0,5
1,5
Figura 2.8. STHD de un horno elctrico de arco en funcin del ancho de ventana T
STHDm
V1
V2
V3
V4
0.415
0.344
0.389
0.280
Tabla 2.6. STHD del horno elctrico de arco para intervalos de tiempo consecutivos
35
La aplicacin del STHD para evaluar las prdidas de potencia activa en el sistema es
posible para casos de seales aperidicas o cargas altamente variables. Esto se ilustra a
partir del caso del horno de arco elctrico discutido previamente. A partir de las medidas
de corrientes realizadas en el sistema y los datos del mismo, se pueden calcular las
prdidas de potencia instantneas, y a partir de stas, la potencia activa perdida, como su
valor medio. En la Figura 2.9 se representan las prdidas en el transformador CT en
funcin del STHD para diferentes puntos de operacin. Se puede observar que el STHD
puede ser utilizado como un indicador de las prdidas de potencia debido a la seal en el
sentido que cuanto mayor es el ndice, mayores son las prdidas. Por tanto, es posible que
la ntima relacin entre las prdidas y el STHD conviertan el ndice en una herramienta
adecuada en aplicaciones de los estndares de calidad de suministro elctrico.
Prdidas(kW )
60
50
40
0.25
0.30
0.35
S TH D
36
(2.29)
G ( n )
n =
= k s
(2.30)
37
g(t)
GM
FT
t
|G(/Ts)|
gn
GM/Ts
DFT
n
2/(NTs)
/Ts
2/Ts
/Ts
Ts Gk G ( k s )
G ( k s )
(2.31)
38
valor
de
g(t)
en
una
discontinuidad
es
0.5[g(+)+g()];
funcin,
g1(t) = etU(t)
es
una
funcin
de
orden-0,
la
segunda,
( Z ) 1
E ( Z ,1) =
sin 2 ( Z )
2
(2.32)
39
(2.33)
Z constante
g1(t) = etU(t)
g2(t) = cos(t)Rect(t)
(a)
g3(t)
g4(t)
(b)
Figura 2.11. a) Funciones de orden-0, g1(t), y orden-1, g2(t); b) Funciones cannica-0, g3(t), y
cannica-1, g4(t).
ENg(h)
101
100 2
3
-1 4
10
1. N = 16
2. N = 32
3. N = 64
4. N = 128
5. N = 256
6. Ideal, E(Z,0)
5
10-2
10
-3
0.1
0.2
0.3
Z = h/N
0.4
0.5
iM (A)
10
5
0
10
15
-5
20 t (ms)
-10
-15
-20
Ih (A)
4.0405
4.0258
4.0391
3.7881
3.7717
3.7870
3.3198
3.3041
3.3192
3.5
2.7010
2.6944
2.7012
2.5
2.5
Exacto
N = 64
N = 256
Ih (A)
4.5
2.0144
2.0270
2.0155
Exacto
N = 64
N = 256
1.3422
1.3815
1.3461
1.5
0.7622
0.8186
0.7645
1
0.5
0
11
13
0.3190
0.3772
0.3215
19
41
ENf(h)
101
1
100
2
10-1
1. N = 64
2. N = 256
3. N = 640
4. N = 1500
5. Ideal, E(Z,1)
3
10-2
4
10-3 5
0
0.04
0.08
0.12
Z = h/N
0.16
0.2
42
iT (A)
5
0
15
10
20 t (ms)
-5
-10
-15
Figura 2.16. Intensidad consumida por una fuente de alimentacin trifsica, iT.
4.5
4.1995
4.0763
4.1912
Ih (A)
3.5
Exacto
N = 64
N = 256
3.4008
3.3082
3.3955
2.7228
2.7627
2.7284
3
2.5
0.6882
0.8548
0.6983
1.2750
1.3222
1.2825
1.5
1
Exacto
N = 64
N = 256
0.8
0.6
0.2618
0.3592
0.2717
0.2173
0.2281
0.2231
0.4
Ih (A)
0.2
0
11
13
0.1761
0.1823
0.1842
19
43
1
0
10
2
-1
1. N = 64
2. N = 256
3. N = 640
4. N = 1500
5. Ideal, E(Z,1)
10
3
10-2
10-3 5
0
0.04
0.08
0.12
Z = h/N
0.16
0.2
44
0.2
0
10
20 t (ms)
15
0.2
0.4
0.6
0.8
0.6409
0.6408
0.6409
9
8
Exacto
N = 64
N = 256
0.6
Ih (mA)
0.4
0.3
0.0595
0.0599
0.0595
0.2
0.1
0
0.0109
0.0114
0.0110
3.2
3.6
3.2
4
3
2.4
2.9
2.4
2
1
0
11
13
19
ENg(h)
0.0214
0.0219
0.0215
Exacto
N = 64
N = 256
4.4
4.9
4.5
0.5
Ih (A)
6.6
7.1
6.6
100
10-1
1. N = 32
2. N = 64
3. N = 256
4. Ideal, E(Z,1)
-2
10
3
-3
10
0.04
0.08
0.12
Z = h/N
0.16
0.2
funcin E(Z,1) ajusta el error cometido prcticamente para cualquier valor de N, y se puede
utilizar dicha funcin para estimar y acotar el error cometido al muestrear la seal.
2.3.3.4 Rectificadores AC/DC de 6 y 12 pulsos
Procediendo del mismo modo que en los casos anteriores, la forma de onda
correspondiente a la corriente consumida por los rectificadores de seis y doce pulsos
ideales, iR6 e iR12, se muestran en las Figuras 2.22 y 2.23 respectivamente, siendo ambas
funciones cannicas-0. Las Figuras 2.24 y 2.25 muestran su contenido armnico de forma
exacta y para dos muestreos, N = 64 y N = 256, para los primeros ocho armnicos.
Tomando como referencia la expresin del error en funciones cannicas-0, (2.32), se
calcul el error cometido para distinto nmero de muestras comprobndose que para ambas
intensidades (iR6 e iR12) la funcin E(Z,0) ajustaba el error cometido de forma exacta para
200
150
100
iR6 (A)
50
0
10
15
20 t (ms)
-50
-100
-150
-200
Figura 2.22. Intensidad consumida por un rectificador AC/DC de seis pulsos ideal, iR6.
400
300
iR12 (A)
200
100
0 0
-100
10
15
20
t (ms)
-200
-300
-400
Figura 2.23. Intensidad consumida por un rectificador AC/DC de doce pulsos ideal, iR12.
46
cualquier valor de N, Figura 2.26. Por tanto se puede utilizar dicha funcin para estimar y
acotar el error cometido al muestrear la seal.
240
200
20
194.06
192.05
193.59
Exacto
N = 64
N = 256
15
14.92
10.98
14.32
Ih (A)
Ih (A)
160
120
80
38.81
39.69
39.21
40
0
27.72
24.69
27.18
7
Exacto
N = 64
N = 256
11.41
9.67
11.67
10.21
5.47
9.54
10
8.43
5.14
8.62
17.64
17.27
17.97
11
17
13
19
23
20
388.13
387.68
388.10
Exacto
N = 64
N = 256
Ih (A)
300
250
200
150
35.28
30.42
34.97
100
50
0
11
29.85
24.16
29.49
16.87
7.38
16.23
13
23
12
Exacto
N = 64
N = 256
11.08
1.71
10.11
10.49
2.71
9.46
8.25
6.13
6.96
4
0
35
25
37
47
101
ENg(h)
Ih (A)
350
16
15.52
5.37
14.82
100
10-1
10-2
10-3
0.04
0.12
0.08
Z = h/N
0.16
0.2
Captulo 3
Rectificadores no controlados con filtro
capacitivo
Las aportaciones de este captulo se centran en el estudio y modelizacin de las fuentes
de alimentacin monofsicas [28-30,32] y trifsicas con filtro capacitivo [31].
Los estudios centrados en la problemtica de los armnicos en los sistemas elctricos de
potencia han ido recibiendo cada vez ms atencin debido al incremento de cargas no
lineales conectadas a los sistemas de distribucin [1, 6]. Entre estas cargas se encuentran
las que utilizan rectificadores monofsicos
id (t)
ZT
ZS
+
i(t)
u(t)
P
C
uC (t) carga
_
_
PCC
Figura 3.1. Esquema del rectificador monofsico con filtro capacitivo
(PCC: punto de acoplamiento comn)
51
52
- La carga del rectificador, que consume una corriente id(t), y se puede modelizar
mediante una fuente de corriente constante. La potencia activa consumida por la
carga es igual a P.
As, en base a lo anterior en este apartado primero se realiza una completa
caracterizacin del comportamiento del dispositivo [28, 30]. A partir de un elevado nmero
de simulaciones realizadas con un programa desarrollado a tal efecto, se presenta el
comportamiento del rectificador de forma grfica tal y como se desarrolla en [26] para el
rectificador trifsico. Dichos grficos permiten analizar de forma clara la influencia de los
diferentes parmetros del dispositivo en su comportamiento.
Posteriormente se desarrolla un modelo analtico muy sencillo que determina el
comportamiento del rectificador monofsico cuando su punto de operacin se define a
partir de su corriente de continua consumida o de su potencia consumida [29, 32]. Para
ambos casos se calcula el ancho y el valor de pico del pulso de corriente, adems de los
armnicos de la corriente de lnea. El modelo asume como hiptesis un valor infinito para
el condensador de filtrado (tensin de continua constante), tensin de alimentacin
perfectamente sinusoidal, y tiene en cuenta el valor de la impedancia del lado de alterna.
En el caso de que la potencia consumida por el convertidor sea un dato conocido, las
expresiones a las que se llega no son lineales. Para este caso, se ha implementado un
algoritmo que analiza el comportamiento del convertidor a partir de ese dato de partida.
El comportamiento del rectificador tambin se caracteriza de manera grfica para los
dos supuestos anteriores, corriente de continua conocida o consumo de potencia conocida,
a partir de las ecuaciones obtenidas anteriormente y mediante la definicin de variables
normalizadas.
La bondad del modelo presentado se ha validado mediante su comparacin grfica con
un modelo exacto que tiene en cuenta un valor no infinito para el condensador de filtrado
(tensin de continua con rizado) (ver apndice).
53
En ltimo lugar, el modelo propuesto permite estudiar de una manera fcil y cmoda los
lmites de contaminacin armnica para rectificadores monofsicos, tomando como
referencia los lmites dados por el estndar IEC 1000-3-2 para equipos de Clase D, [74].
Para estos dispositivos la norma limita la corriente armnica mxima admisible por vatio
consumido, por lo que se deber calcular el ratio anterior a partir de los resultados
obtenidos con el modelo, armnicos de corriente de lnea, a fin de aplicar dicha norma.
SCV U 2 / X
U2
=
=
P
P
X P
(3.1)
XC
1
=
X
2 f LC
(3.2)
54
xC = 0
xC = 2.5
uC (t)
RSC = 25
i(t)
u(t)
RSC = 100
t
RSC = 1000
Figura 3.2 Formas de onda de la corriente de lnea del lado de alterna i(t)
y de la tensin del lado de continua uC(t)
Se puede observar en la Figura 3.2 la fuerte influencia que tiene el parmetro RSC en la
forma de la corriente de lnea en comparacin con el parmetro xC. Por este motivo, RSC es
uno de los parmetros en funcin del cual se presentan todos los resultados obtenidos.
Por otro lado, el parmetro xC influye en el rizado que presenta la tensin de continua.
De este modo, xC deviene el segundo parmetro en funcin del cual se presentan los
resultados.
Al final del presente apartado, se considera y estudia la influencia que la resistencia de
entrada presenta en el comportamiento del dispositivo. Para ello se utiliza la relacin
resistencia-reactancia del lado de alterna, r, como parmetro de estudio. La definicin del
parmetro r es la siguiente:
r=
R
X
(3.3)
55
As se realizar el anlisis del dispositivo para los siguientes rangos de los parmetros
anteriores:
- Relacin de cortocircuito: 5 RSC 1000.
- Relacin de reactancia: xC = 0, 1.5, 2.5, 3.5, y 4.5.
- Relacin
resistencia-reactancia:
En
primer
lugar
se
estudia
el
caso
evolucin de t, tal como se vio en la Figura 3.2, es tal que a medida que se incrementa el
valor de RSC, el ancho del pulso disminuye, mientras que la influencia de xC se va haciendo
cada vez menor.
Se observa de la figura que los anchos habituales de la corriente consumida por estos
dispositivos, t 3 ms 54, se alcanzan con un valor de RSC aproximadamente igual a 70.
As valores inferiores de esa relacin llevan en realidad a situaciones poco frecuentes.
Tanto en esta Figura 3.3 como en las sucesivas, los tramos de curva representados con
lnea discontinua corresponden a casos en los que el valor del rizado de la tensin de
continua supera el 40%.
140
120
g
t (grados)
100
80
60
40
20
10
10
10
RSC = U /(X P)
57
135
C = UC/U (%)
130
125
(a)
120
115
h
110
105
10
10
10
RSC = U /(X P)
50
45
40
up = Up/U (%)
35
30
25
(b)
20
15
g
10
5
0
10
10
10
RSC = U /(X P)
Figura 3.4. Caracterizacin de la tensin de continua, con xC = 0 (n), 1.5(g), 2.5 (),
3.5 (t), 4.5 (u): a) valor medio, c, b) valor pco a pico (rizado), up
58
59
118
80
114
75
112
g1 (%)
i1 (%)
85
116
110
70
65
108
g
106
104
55
102
100
60
10
50
10
10
10
10
RSC = U /(X P)
1
0.98
10
RSC = U /(X P)
80
75
0.96
70
0.92
65
(%)
cos 1 (%)
0.94
0.9
60
0.88
t
u
55
0.86
50
0.84
0.82
10
10
10
45
10
RSC = U 2/(XPd)
10
2
RSC = U /(X Pd )
10
90
90
h=3
70
h=9
30
20
h=5
10
10
10
h = 11
20
10
40
30
50
10
10
10
RSC = U /(XPd)
RSC = U /(XPd)
80
70
60
50
ih (%)
ih (%)
ih (%)
40
10
60
50
70
60
80
80
40
30
h=7
10
0
20
h = 13
10
10
10
RSC = U /(XPd)
Figura 3.6. Corrientes armnicas de lnea ih = Ih/I1, con xC = 0 (n), 1.5 (g), 2.5 (), 3.5 (t),
4.5 (u): a) secuencia homopolar, h = 3, 9; b) secuencia negativa, h = 5, 11;
c) secuencia positiva h = 7, 13
- La clara dominancia del tercer armnico (dominio del armnico de menor orden)
y nicamente para relaciones de cortocircuito elevadas (reactancia de alterna o
potencias consumidas reducidas) los armnicos quinto y sptimo alcanzan
magnitudes comparables a la del tercer armnico.
- Los rangos habituales de emisin de armnicos son 50% i3 95%,
40% i5 85%, 10% i7 75%, y 5% i9 60%.
- Valores elevados de los elementos de filtrado, inductancia y capacidad, reducen la
cantidad de armnicos bajos emitidos.
- En concreto, valores elevados de la reactancia del lado de alterna (relacin de
reactancia pequea) reducen dicha cantidad de armnicos [33]. As, una forma de
61
1/ 2
(3.4)
160
THD (%)
140
120
100
80
g
60
40
10
10
10
RSC = U 2/(X P)
20
40
60
100
80
Re (%)
210
u
240
t
1000
5
500 330
250
25
50
h 100
RSC
Im (%)
300
270
63
En la Figura 3.8 (h = 3), se puede observar que el valor normalizado para el tercer
armnico, i3, vara entre un 50% y un 95%, y el ngulo de fase, 3, vara entre 270 y 360.
En general, a medida que la relacin de cortocircuito aumenta, el ngulo de fase tambin lo
hace, pero se ha comprobado que para valores elevados del condensador de filtro, hasta un
valor de RSC igual a 500, el ngulo de fase decrece. En la Figura 3.9 (h = 5 y 7), se puede
comprobar la misma tendencia para las variables i5, i7, 5, y 7: 40% i5 85%,
10% i7 75%, y 180 5 360, 270 7 330.
22.5
67.5
45
90
Re (%)
210
u
u
25
t
t
50
240
330
100
1000
250
500
RSC
Im (%)
300
270
25
20
40
210
60
80
0
Re (%)
u
u
h
h
50
100
240
t
t
330
250
500
RSC
1000
Im (%)
300
270
Finalmente, como se ilustra en las Figuras 3.8 y 3.9, la dispersin en los ngulos de fase
de las corrientes armnicas es debida principalmente al nivel de potencia consumida y el
valor de la reactancia de entrada segn lo expuesto en [33]. Tal como se comentar
posteriormente, la constatacin de esta evolucin es importante en los estudios sobre la
cancelacin de armnicos entre dispositivos no lineales [36].
3.1.2.1.6 Lmites para las emisiones de corriente armnica
La norma establecida por la Comisin Electrotcnica Internacional IEC 1000-3-2 [74]
prescribe los lmites de contaminacin armnica para equipos entre los que se engloban las
fuentes de alimentacin monofsicas. Dentro de esta norma, las fuentes de alimentacin
con una potencia de entrada inferior a 600 W se consideran como equipos de Clase D, y el
lmite de su emisin armnica se define a partir de la corriente mxima admisible por vatio
consumido (mA/W), es decir, la relacin entre las Ih y P. Las condiciones en la norma son
aplicables, adems, para equipos diseados para ser conectados a sistemas de 220/380 V,
230/400 V y 240/415 V operando a 50 60 Hz.
Un examen de la normativa indica que la limitacin en el tercer armnico es la ms
restrictiva, por lo que a partir de las Figuras 3.5 (i1 = I1/Ili) y 3.6 (i3 = I3/I1) y para U igual a
3.5
3.4 mA/W
I3 / P (mA/W)
2.5
1.5
10
10
2
RSC = U /(XP)
Figura 3.10. Emisin de corriente armnica (h = 3) por vatio consumido (lnea continua: U = 220V,
lnea discontinua: U = 240 V), con xC = 0 (n), 1.5(g), 2.5 (), 3.5 (t), 4.5 (u)
65
Zu
Ztr
+
Zu
Pu
Pu
I hN
,
N I h1
(3.5)
donde Ih1 es la corriente armnica consumida cuando slo hay una carga conectada (N = 1)
e IhN es la corriente armnica consumida cuando hay N cargas idnticas a la anterior
conectadas en paralelo. Los grficos de la Figura 3.12 han sido calculados para el caso
mostrado en la Figura 3.11, con N igual a 5, 10, 25 y 75 cargas no lineales idnticas
(potencia unitaria Pu e impedancia de la lnea Zu iguales). Para ajustar el estudio del
circuito propuesto al del circuito mostrado previamente en la Figura 3.1, la impedancia Z
67
se calcula como NZtr + Zu, y se considera una nica carga de valor igual a la potencia total
consumida por las N cargas.
Si no existiese el fenmeno de la atenuacin, IhN debera coincidir con NIh1 (las
intensidades Ih1 tienen la misma fase debido a que las cargas son idnticas), pero esto no es
as, ya que un incremento en el nivel de potencia consumida por la carga afecta a la
corriente consumida por el convertidor.
Los grficos muestran los factores de atenuacin correspondientes a las corrientes
armnicas predominantes (AF3, AF5 y AF7) para un rango usual de valores de la relacin de
cortocircuito (50 < RSC < 500), dos relaciones de resistencia-reactancia de entrada (r = 0 y
r = 2) y tensin de continua sin rizado (xC = 0).
De los resultados obtenidos se observa que:
- La atenuacin es ms elevada para los armnicos de mayor orden.
- El efecto de la atenuacin tiende a incrementarse con N, pero para un nmero de
cargas superior a 75 este incremento no es significativo.
- A medida que la relacin de cortocircuito se incrementa de modo que se tiende a
un bus stiff (por ejemplo, sistema de impedancia cero), el factor de atenuacin
tiende a la unidad.
- Cuanto mayor es el valor de la relacin resistencia-reactancia de alterna, r = R /X,
mayor es el efecto de atenuacin, pero su influencia es menor que la de la relacin
de cortocircuito.
- Los valores usuales de la relacin de reactancias proporcionan factores de
atenuacin ligeramente mayores que los correspondientes al valor ideal xC igual a
0. Estos valores se encuentran ms prximos a medida que la relacin de
cortocircuito aumenta.
- El caso presentado en [35] se analiza para N = 5 (RSC = 115.20, xC = 2.62, r = 1.60
y AF3 = 0.91, AF5 = 0.72, AF7 = 0.48) y N = 10 (RSC = 57.60, xC = 1.85, r = 1.60 y
68
90
90
80
80
70
70
AFh (%)
100
AFh (%)
100
60
50
60
50
40
40
30
10
50
30
Punto A
20
100
RSC = U 2/(XP)
10
50
500
90
90
80
80
70
70
AFh (%)
100
AFh (%)
100
60
50
30
30
20
N = 25
100
RSC = U 2/(XP)
100
500
RSC = U 2/(XP)
500
50
40
10
50
N = 10
60
40
20
Punto B
20
N=5
10
50
N = 75
100
RSC = U 2/(XP)
500
Figura 3.12. Factores de atenuacin para las corrientes armnicas en funcin de RSC y r (r = 0, lnea
continua, y r = 2, lnea discontinua), AF3 (), AF5 ( ) and AF7 ( )
AF3 = 0.84, AF5 = 0.54, AF7 = 0.23). Los puntos A y B indicados en la Figura 3.12
muestran los resultados de ambos casos, respectivamente.
El efecto de la cancelacin de corrientes armnicas se produce al conectar a un mismo
sistema cargas con diferentes niveles de consumo de potencia, o diferentes caractersticas
de lneas de entrada (diferente valor de reactancia, diferente valor de relacin de
resistencia-reactancia), debido a la diversidad en los ngulos de fase (Figuras 3.8 y 3.9)
[36].
En la presente tesis slo se ha analizado y presentado resultados para el efecto de la
atenuacin, con lo que los grficos adjuntos seran vlidos solamente en el caso de
considerar circuitos como el mostrado en la Figura 3.11 con N cargas no lineales idnticas.
69
uC
IDC
En la realizacin del estudio se han considerado dos casos para definir el modelo,
dependiendo de cmo se considere la impedancia del lado de alterna. Estos casos son los
siguientes:
- Caso resistivo: la impedancia se considera una resistencia, Z = R.
- Caso inductivo: la impedancia se considera una inductancia, Z = jX = j L.
70
De acuerdo con la Figura 3.14, el anlisis realizado se basa en las dos hiptesis
siguientes:
- La tensin de alimentacin, u, es perfectamente sinusoidal, siendo U su valor
eficaz.
- El valor del condensador de filtrado, C, es lo suficientemente elevado como para
considerar que la tensin en bornes del condensador es constante y, segn la Figura
3.14, de valor
U C = 2 U cos
(3.6)
UC
i, caso inductivo
u
ImL
I mR
i, caso resistivo
= t
tR
tL
T1
T2
T3
2
sT3 = 2 U 1 + ,
2
sT1 = sT2 = 2 U
(3.7)
uT2 ( ) = 2 U ( 2 )
2
uT3 ( ) = 2 U 1 + ( 2 ) .
2
uT1 ( ) = 2 U
(3.8)
72
2R
U
( 2 )
2R
(3.9)
L 0
2 2L
1
U
4 2 2 2 )
iT2 ( ) = iT1 ( ) +
uT2 ( ) d =
(
L
2 2L
iT1 ( ) =
iT3 ( ) = iT2 ( 2 ) +
(3.10)
1
U 2
2
uT3 ( ) d =
1 + ( 2 )
2
2
L
2L
iT1 ( ) d + iT2 ( ) d +
2 +
iT3 ( ) d = I DC
(3.11)
2 R U
U R
R
d +
( 2 R ) d = I DC
R
2R
2R
(3.12)
R = 1.6439 3
I DC R
U
(3.13)
73
I DC R
U
, R =
t R
2
(3.14)
2 L U
U L
L
4 L 2 L2 2 ) d +
2 d +
(
L
2 2 L
2 2 L
2 L + L U
2
L 2
+
L 1 + 2 ( 2 L ) d = I DC
2 L
2 L
(3.15)
L = 5.4667 4
I DC L
U
(3.16)
I DC L
U
, L =
t L
+2
(3.17)
En resumen, las expresiones resultantes que determinan el ancho de pulso son, para los
casos resistivo e inductivo, respectivamente:
t R = 0.0105
t L = 0.0494
I DC R
t R
, R =
2
U
I DC L
t L
, L =
U
+2
(3.18)
74
I mR =
U R2
U L3
; I mL =
.
2R
2 L
(3.19)
ImR
ImL
Im
(a)
2L L
(b)
b/2
b/2
(c)
b
b
2
2
(3.20)
donde el valor mximo se considera igual a los valores mximos de la forma de onda
original, es decir, ImR = Im o ImL = Im, segn sea el caso resistivo o inductivo,
respectivamente.
El parmetro b es el ancho normalizado del pulso de corriente. La manera de calcular el
valor de b es considerando que las reas bajo las diferentes formas de onda son de idntico
valor. As, para el caso resistivo, se impone que las reas bajo las formas de onda
75
correspondientes a las Figuras 3.15.a y 3.15.c son iguales, resultando bR = 0.5R. De forma
anloga para el caso inductivo, se obtiene que bL = 0.7798L.
As, la expresin propuesta en [38] para el clculo de los armnicos de la corriente de
lnea viene dada por
Ih =
Im
Kh
2
cos hb
4b
2
; Kh =
1 ( hb )2
( h = 1,3,5, 7,)
(3.21)
(3.22)
PUNTO DE
TRABAJO
U, R, L, IDC
CARACTERIZACIN
DE LA CORRIENTE
t R = 0.0105
I DC R
, t L = 0.0494
U
2
I mR
U tR
=
, I mL =
2R 2
I DC L
U
tL
2 L 2 +
U
I m = Min { I mR , I mL }
tL > tR
No
t R
b = bR = 0.5
Ih =
t L
b = bL = 0.7798
2+
cos hb
2
2
2 1 ( hb )
I m 4b
CLCULO DE LAS
CORRIENTES
ARMNICAS
77
PR
U CU
; t L = 0.0494
PL
U CU
(3.23)
PR
t R
, R =
2
2
U
PL
t L
, L =
2
U
+2
(3.24)
78
PUNTO DE
TRABAJO
DATOS
(U, R, L, P)
(0)
I DC
=
t R = 0.0105
I DC R
U
DATOS
(U, R, L, IDC)
P
2 2U
, t L = 0.0494
I DC L
U
tL > tR
No
I DC =
No
P
t R
2U cos
I DC =
P
t L
2U cos
2+
CARACTERIZACIN DE LA CORRIENTE
79
L
,
R
(3.25)
; Z=
L
1 + Q2
Q
(3.26)
I DC Z
U 1 + Q2
; t L = 0.0494
I DC ZQ
U 1 + Q2
(3.27)
I DC
I SC 1 + Q 2
; t L = 0.0494
I DC Q
I SC 1 + Q 2
(3.28)
Finalmente, si tomamos iDC como la relacin entre la corriente continua consumida por
la carga alimentada y la corriente de cortocircuito del dispositivo,
iDC =
I DC
I
= DC
I SC (U / Z )
(3.29)
El desarrollo del modelo exacto que caracteriza el comportamiento de la fuente de alimentacin monofsica
se encuentra recogido en el apndice (apartado A.1).
80
iDC
1+ Q
; t L = 0.0117
iDC Q
1 + Q2
(3.30)
Utilizando los mismos parmetros Q e iDC en las expresiones del valor mximo de la
corriente de lnea (3.19), se obtienen las siguientes relaciones:
I mR
2
1+ Q2
= 1.923 3 iDC
I SC
i3 1 + Q 2
I mL
= 1.534 4 DC
I SC
Q
(3.31)
iDC
1+ Q
; bL = 1.009 4
iDC Q
1 + Q2
(3.32)
81
CU
I DC
(3.33)
82
t (ms)
Q = L/R = 2.0
Q = L/R = 0.1
(a)
3
2
1
0
10
0.35
Imax/ISC
0.3
0.25
0.2
(b)
Q = L/R = 0.1
Q = L/R = 2.0
0.15
0.1
0.05
0
0
10
0.9
0.8
0.7
Ih/I1
0.6
iDC1 = 500
0.5
iDC1
0.4
0.3
(c)
= 100
iDC1 = 10
0.2
iDC1 = 30
0.1
0
1
11
Armnico (h)
I mL
3
1.534 4 iDC
, bL 1.009 4 iDC
I SC
(3.34)
de forma que no dependen del valor de Q. Por este motivo, en la Figura 3.18c, las curvas
obtenidas no dependen de Q, ya que se han obtenido para valores de Q mayores que 2.0.
Esto se ha realizado as porque es habitual considerar nicamente la reactancia debido al
elevado valor del ratio Q que presentan estos sistemas.
3.1.4.2 Caracterizacin por la potencia consumida
Para validar el modelo propuesto caracterizado por la potencia consumida en el lado de
continua, de manera anloga al modelo de la corriente de continua, se utiliza la variable
normalizada Q definida en el apartado anterior, y una segunda variable anloga a iDC que
reflejar el dato relativo a la potencia consumida por la carga conectada a la fuente. Esta
84
P
U I SC
(3.35)
= 0.00936
t L
t L 4 cos
= 0.0107
2.839
p
3
1 + Q2
Q
(3.36)
1 + Q2
1 + Q2
I mR
2
3
= 1.5267 p
I SC
cos 2 R
1 + Q2
I mL
3
4
= 1.1829 p
I SC
Q cos3 L
(3.37)
b
bR 3 cos R = 0.3671
0.25
p
3
bL
bL 4 cos
= 0.9253
0.7798
1 + Q2
4
(3.38)
1 + Q2
85
86
t (ms)
Q = L/R = 2.0
Q = L/R = 0.1
(a)
A
3
2
1
0
10
12
p (%) = 100P/UISC
0.45
0.4
0.35
Imax/ISC
0.3
0.25
(b)
0.2
Q = L/R = 0.1
Q = L/R = 2.0
0.15
0.1
0.05
0
10
12
p (%) = 100P/UISC
1
0.9
0.8
Ih/I1
0.7
p = 0.1%
0.6
p 2%
0.5
(c)
p = 0.5%
0.4
0.3
0.2
p = 12%
p = 1.5%
0.1
0
1
11
Armnico (h)
La Figura 3.20a muestra el grfico del ancho del pulso de corriente obtenido mediante
el modelo propuesto para los casos resistivo e inductivo, tR y tL, segn (3.18), y
mediante el mtodo exacto considerando BC = 25, t.
La Figura 3.20b muestra el grfico del valor mximo de la corriente de entrada obtenida
mediante el modelo propuesto para los casos resistivo e inductivo, ImR e ImL, segn (3.19), y
mediante el mtodo exacto considerando BC = 25, Im.
As, las Figuras 3.20a y 3.20b permiten analizar la bondad del mtodo propuesto y
cmo se adapta al caso de una impedancia con componente resistiva e inductiva. A partir
de la observacin de los resultados obtenidos se desprende que el error numrico cometido
es aceptable siguiendo el criterio expuesto en (3.22) en el clculo del ancho, del valor
mximo del pulso de corriente, y el parmetro b.
88
6.0
iDC = 10 %
tL
5.5
t (ms)
5.0
4.5
tR
(a)
4.0
tL
3.5
3.0
iDC = 1 %
2.5
tR
2.0
1.5
0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
1.6
1.8
2.0
Q
0.55
0.50
iDC = 10 %
ImR
0.45
Im (A)
0.40
0.35
0.30
ImL
0.25
0.20
ImR
iDC = 1 %
0.15
Im
0.10
0.05
0
(b)
Im
ImL
0.2
0.4
0.6
0.8
1.0
1.2
1.4
1.6
1.8
2.0
En los grficos, el anterior criterio se ha indicado marcando con trazo continuo los
intervalos de curva, tR y tL en la Figura 3.20a, y ImR e ImL en la Figura 3.20b, que mejor
aproximan las curvas t e Im, respectivamente. Tal como se coment al presentar la Figura
3.18, para valores usuales de la impedancia, Q > 0.5, se observa que las expresiones
correspondientes al caso inductivo ajustan el resultado exacto de una manera aceptable.
89
90
90
u
t
i3 (%)
80
70
g
60
50
40
30
10
10
10
RSC = U /(XPd)
91
92
93
Corriente armnica
mxima admisible por
vatio (mA/W), (Ih/P)lim
3.4
1.9
1.0
0.5
11
0.35
13 h 39
3.85/h
Ih / P U
( I h / P )lim 230
(3.39)
94
X P
U2
(3.40)
obtenindose las Figuras 3.22a (para los armnicos tercero, quinto y sptimo) y 3.22b (para
los armnicos noveno, undcimo y decimotercero). En dichos grficos, por ejemplo, el
valor ihp = 1 corresponde al lmite armnico indicado por la normativa [74] considerando
una tensin de alimentacin U = 230 V.
Los armnicos ms crticos (h = 3, 5 y 7) se estudian con ms detalle en el entorno de
U = 230 V (220 240 V) en la Figura 3.23, lo que supone trabajar con un rango de valores
de x(%) ms restringido. Por un lado, se dispone de los resultados correspondientes al
empleo del modelo propuesto (condensador de filtrado ideal, BC = ); por otro lado, se
obtienen los resultados correspondientes al empleo del modelo exacto, en el que se
considera la existencia de un condensador no ideal. Como en los casos anteriores, se
considera un valor de la relacin de susceptancias igual a BC = 25. A partir de la Figura
3.23, se deduce que para un valor dado de ihp, el modelo propuesto proporciona un valor de
la reactancia, x(%), ms pequeo que el correspondera al obtenido con el modelo exacto.
Este aspecto debe ser tenido en cuenta a la hora de aplicar el modelo propuesto. Por
ejemplo, para x igual a 2.25%, el valor real de i3p es igual a 1.05, mientras que el valor
obtenido para i3p con el modelo propuesto es igual a 1.00, por lo que realmente se est
excediendo en un 5% el lmite armnico indicado en la normativa.
Por otro lado, el estudio anterior se puede realizar tambin numricamente mediante los
algoritmos expuestos en el apartado 3.1.3.
3.1.6.1 Ejemplo 1: Cumplimiento de la normativa
Conocida las condiciones de operacin y carga del sistema, U y P, as como la
reactancia del mismo, x(%), se desea averiguar su situacin desde el punto de vista de la
emisin armnica que produce.
El procedimiento a seguir es el mismo que se explic en apartados anteriores, puesto
que se trata de caracterizar la corriente consumida por el dispositivo y, posteriormente,
comparar los resultados obtenidos con los valores lmite que dicta la normativa. As,
95
ihp
2.0
(a)
1.5
h=3
1.0
h=5
0.5
0
1.5
ihp
(b)
3
2
h=9
h = 11
1
0
0.1 0.5 1.0
1.5 2.0
h = 13
2.5 3.0 3.5 4.0
4.5 5.0
x = 2.75%
ihp
1.2
1.0
h=3
h=7
0.8
0.6
1.0
h=5
Modelo propuesto
(BC = )
Modelo exacto (BC = 25)
1.4
1.8
2.2
2.6
3.0
(3.41)
k +1)
( k +1)
x1
1 (k )
k
k +1
k
x1 + x2( ) ; x2( ) = x2( )
2
(k )
( k +1)
= x1 ; x2
1 k
k
= x1( ) + x2( )
2
) 0
1 k
k
si p x1( ) + x2( )
2
1 k
k
si p x1( ) + x2( )
2
>0
(3.42)
97
x U 2 2.252302
=
= 11.90
100 P
100100
(3.43)
98
id(t)
Z
ZS
+
u(t)
_
ZS
i(t)
ZT
ZT
C
ZS
uC(t) carga
ZT
_
99
100
Por otro lado, el modelo propuesto permite estudiar de una manera fcil y cmoda los
lmites de contaminacin armnica para rectificadores trifsicos, tomando como referencia
los lmites dados por el estndar IEC 1000-3-2 para equipos de Clase A si el rectificador
consume una corriente inferior a 16 A [74], y por el estndar IEC 1000-3-4 si la corriente
consumida es superior a 16 A [75].
En estos dispositivos, a diferencia de los rectificadores monofsicos, la norma limita las
corrientes armnicas mximas admisibles por lo que los resultados obtenidos por el
modelo desarrollado sern directamente comparables con las normas.
P
+
uC
IDC
101
(3.44)
UC
i, caso inductivo
u
ImL
I mR
i, caso resistivo
= t
tR
tL
T1
T2
T3
2
sT3 = 2 U 1 + ,
2
sT1 = sT2 = 2 U
(3.45)
uT2 ( ) = 2 U ( 2 )
2
uT3 ( ) = 2 U 1 + ( 2 ) .
2
uT1 ( ) = 2 U
(3.46)
103
valor uC, segn [37], la aproximacin realizada permitir obtener un modelo del
rectificador menos exacto pero extremadamente ms sencillo.
A partir de las expresiones de la tensin u, la corriente de lnea del lado de alterna, i,
se puede calcular a partir de las conocidas relaciones entre la tensin en bornes de una
resistencia o una inductancia y la corriente que circula por ella. Para el caso resistivo se
obtiene, en cada intervalo:
U
22 R
U
( 2 )
22 R
uT 1
=
2R
u
iT2 ( ) = T 2 =
2R
iT3 ( ) = 0
iT1 ( ) =
(3.47)
2 L
uT1 ( ) d =
iT2 ( ) = iT1 ( ) +
U
2
2 2 ( 2 L )
U
u ( ) d =
( 4 2
2 L
2 2 ( 2 L )
iT3 ( ) = iT2 ( 2 ) +
T2
u ( ) d =
2 L
2
T3
2 )
(3.48)
2
U
2
1 + 2 ( 2 )
2 ( 2 L )
I DC R
t R
, R =
2
U
I DC L
t L
, L =
U
+2
(3.49)
U R2
U L3
; I mL =
2 ( 2R )
2 ( 2L )
(3.50)
ImR
ImL
(a)
Im
2L L
(b)
b/2
b/2
(c)
onda del pulso de la corriente de lnea para el caso resistivo, Figura 3.27a, e inductivo,
Figura 3.27b, caracterizada a partir de los parmetros R, ImR, o L, ImL, y que ya han sido
calculados en las expresiones 3.49 y 3.50, se aproximar a la forma de onda casi-sinusoidal
[38], Figura 3.27c.
La senoide, Figura 3.27c, viene dada por la expresin
i ( ) = I m cos
b
b
b
2
2
(3.51)
donde el valor mximo se considera igual a los valores mximos de la forma de onda
original, es decir, ImR = Im o ImL = Im, segn sea el caso resistivo o inductivo,
respectivamente.
El parmetro b es el ancho normalizado del pulso de corriente. La manera de calcular el
valor de b es considerando que las reas bajo las diferentes formas de onda son de idntico
valor. As, para el caso resistivo, se impone que las reas bajo las formas de onda
correspondientes a las Figuras 3.27a y 3.27c son iguales, resultando bR = 0.5R. De forma
anloga para el caso inductivo, se obtiene bL = 0.7798L.
As, la expresin propuesta en [38] para el clculo de los armnicos de la corriente de
lnea, suponindola constituida por un nico pulso, vena dada por la expresin (3.21). En
el caso de la fuente de alimentacin trifsica la corriente de lnea est constituida por dos
pulsos iguales decalados /3 por lo que el desarrollo armnico sera la suma de los
desarrollos presentados en la expresin (3.21) con el desfase adecuado, es decir,
ITh = Ih + Ih1h/3, lo que lleva a ITh = 3 Ih, h = 1, 5, 7, 11, ... Por ello, la expresin
para el clculo de los armnicos de la corriente de lnea de la fuente trifsica viene dado
por
I Th = 3
106
Im
Kh
2
cos hb
4b
2
; Kh =
1 ( hb )2
h = 1,5, 7,11
(3.52)
(3.53)
Una vez definido el criterio a seguir en la resolucin del caso general, se expone el
procedimiento analtico a seguir en la resolucin de un problema de este tipo (Figura 3.28):
a partir de los datos del dispositivo (U, R, L, IDC), en primer lugar se caracteriza la corriente
consumida por la fuente de alimentacin suponiendo que el problema a resolver es
puramente resistivo y puramente inductivo, (3.49) y (3.50), respectivamente; a
continuacin se aplica el criterio (3.53), caracterizando definitivamente la corriente de
entrada; finalmente, se calculan las componentes armnicas de dicha corriente, (3.52).
3.2.2.2 Caracterizacin por la potencia consumida
Aunque los fabricantes dan el contenido armnico a la potencia nominal del
rectificador, dichos dispositivos suelen trabajar al 35% 60% de su potencia nominal [27]
por lo que es interesante determinar dicho contenido armnico a otras potencias.
En el caso de definir como dato de la fuente de alimentacin la potencia que debe
suministrar a la carga de continua conectada, el ancho del pulso de corriente se puede
obtener a partir de la expresin (3.49) reemplazando su consumo de corriente en el lado de
107
PUNTO DE
TRABAJO
U, R, L, IDC
CARACTERIZACIN
DE LA CORRIENTE
tR = 9.143103
I mR =
I DC R
, tL = 4.150102
U
2
I DC L
U
tL
U
tR
, I mL =
2 ( 2R ) 2
2 ( 2L ) 2 +
U
I m = Min { I mR , I mL }
tL > tR
No
t R
b = bR = 0.5
t L
b = bL = 0.7798
2+
CLCULO DE LAS
CORRIENTES
ARMNICAS
I Th
cos hb
I m 4b
2
= 3
2
2 1 ( hb )
PR
t R
3
= 8.14510 3 2 , R =
U
2
t L
PL
t L
2
t L 4 cos
= 3.80610 4 2 , L =
U
+2
+2
108
(3.54)
PUNTO DE
TRABAJO
DATOS
(U, R, L, P)
(0)
=
I DC
t R = 9.143103
I DC R
U
DATOS
(U, R, L, IDC)
P
3 2U
, t L = 4.150102
I DC L
U
tL > tR
No
I DC =
No
P
t R
2U cos
I DC =
P
t L
2U cos
2+
CARACTERIZACIN DE LA CORRIENTE
109
L
R
(3.55)
; Z=
L
1 + Q2
Q
(3.56)
El desarrollo del modelo exacto que caracteriza el comportamiento de la fuente de alimentacin trifsica se
encuentra recogido en el apndice (apartado A.2).
110
t R = 9.143103 3
I DC Z
; t L = 4.150102 4
U 1 + Q2
I DC ZQ
U 1 + Q2
(3.57)
I DC
I SC 3 1 + Q 2
; t L = 4.150102 4
I DC Q
I SC 3 1 + Q 2
(3.58)
Finalmente, si tomamos iDC como la relacin entre la corriente continua consumida por
la carga alimentada y la corriente de cortocircuito de la red,
iDC =
I DC
I
= DC
U
I SC
3Z
(3.59)
iDC
1+ Q
; t L = 8.593103 4
iDC Q
1 + Q2
(3.60)
Utilizando los mismos parmetros Q e iDC en las expresiones del valor mximo de la
corriente de lnea (3.50), se obtienen las siguientes relaciones:
I mR
2
1+ Q2
= 0.876 3 iDC
I SC
3
iDC
1 + Q2
I mL
4
= 0.655
I SC
Q
(3.61)
iDC
1+ Q
; bL = 0.7414
iDC Q
1 + Q2
(3.62)
armnicas de la corriente se calculan a partir de los valores del ancho y valor mximo del
pulso de la corriente, lo que se traduce tambin en una dependencia final de las variables
anteriores. As, las variables iDC y Q sern las que se utilicen para presentar la validacin
grfica del modelo propuesto.
El modelo propuesto se valid a partir de su comparacin con el modelo exacto del
dispositivo, para el rango de valores de la variable iDC comprendido entre el 0% y el 1.4%
y para dos valores de la variable Q, 0.1 y 2.0.
El lmite superior de la variable iDC viene fijado por la condicin de conduccin
discontinua; as, dicha condicin lleva a que el ancho del pulso de la corriente no pueda
exceder un sexto del periodo, es decir, t < T/6 = 3.33 ms a 50 Hz, por ser el momento en
el que los dos pulsos de intensidad se solaparan [26, 27]. Imponiendo esta condicin en la
expresin 3.49 para el caso inductivo se obtiene un resultado de iDC = 1.68%.
Los lmites de la variable Q se han adoptado por similitud con el estudio del rectificador
monofsico aunque es conocido que los rectificadores trifsicos suelen presentar un valor
elevado, Q 5 [27], y por tanto el estudio para Q reducidos no tendra tanta utilidad como
en el caso monofsico.
Respecto al modelo exacto utilizado, se considera la presencia de un condensador
caracaterizado con el parmetro,
BC =
CU
I DC
(3.63)
La Figura 3.30 recoge los resultados obtenidos para los dos casos comentados, exacto y
aproximado, para las variables que caracterizan la corriente consumida por la fuente de
alimentacin: su ancho del pulso, su valor mximo (expresado en base a la corriente de
cortocircuito) y sus componentes armnicos (expresados en base a la corriente
fundamental consumida). Finalmente, de los dos estudios exactos que se realizaron, 2.5% y
10% de rizado, se muestra en los grficos el ms extremo (10% de rizado, es decir, BC
igual a 5).
De la observacin de los grficos, se desprenden conclusiones similares a las
presentadas para el rectificador monofsico, es decir,
- Los grficos anteriores, y en general su determinacin a partir de las variables Q e
iDC, se revelan como una forma sencilla, cmoda y muy til de caracterizar una
fuente de alimentacin trifsica a partir de sus parmetros, o como herramienta de
diseo, en el sentido de poder determinar parmetros de la carga que permitan
cumplir unos requisitos mnimos en relacin con la contaminacin armnica que
el dispositivo va a producir en el sistema.
- El error del modelo aproximado propuesto en este captulo es aceptable en un
amplio rango de valores de impedancia de entrada. Solamente el error cometido
en el clculo del valor mximo de la corriente, para pequeos valores de Q
(impedancias puramente resistivas) y elevados valores de iDC (elevado consumo
del rectificador) se podra considerar como no aceptable.
- Cuanto mayor es el valor de iDC, mayor es el ancho del pulso de corriente y el
valor mximo de la misma. Analticamente, esto se refleja en las expresiones
(3.50) y (3.51), respectivamente.
- Tal como se ha comentado, los rectificadores para potencias elevadas suelen
presentar un valor de Q 5 [27], por lo que las corrientes armnicas, Figura
3.30c, se han obtenido considerando un valor de Q superior a 2.0. Segn se ve en
esa figura y como era de esperar, cuanto ms estrecho es el pulso de corriente,
mayor es su contenido armnico.
113
2.8
2.6
2.4
Q = L/R = 2.0
t (ms)
2.2
(a)
Q = L/R = 0.1
2
1.8
1.6
1.4
1.2
0.2
0.4
0.6
0.8
1.2
1.4
0.03
Imax/ISC
0.025
0.02
Q = L/R = 0.1
(b)
Q = L/R = 2.0
0.015
0.01
0.005
0.2
0.4
0.6
0.8
1.2
1.4
iDC-1 = 325
0.9
0.8
iDC-1 = 60
0.7
iDC-1 = 500
ITh/I1
0.6
(c)
0.5
iDC-1 = 200
0.4
0.3
0.2
0.1
0
1
11
Armnico (h)
Figura 3.30. Caracterizacin de la corriente de lnea del lado de alterna a partir de la corriente
consumida: a) ancho del pulso de corriente en funcin de iDC para dos valores de Q; b) valor
mximo de corriente en funcin de iDC para dos valores de Q; c) armnicos de corriente para
diferentes valores de iDC y Q > 2.0
114
I mL
3
0.655 4 iDC
, bL 0.741 4 iDC
I SC
(3.64)
de forma que no dependen del valor de Q. Por este motivo, en la Figura 3.30c, las curvas
obtenidas no dependen de Q, ya que se han obtenido para valores de Q mayores que 2.0.
Esto se ha realizado as porque es habitual considerar nicamente la reactancia de entrada
debido al elevado valor del ratio Q que presentan estos sistemas [27].
115
P
U I SC
(3.65)
3
= 6.78210
p
3
t L
3
t L 4 cos
= 7.88010
2.839
1 + Q2
Q
(3.66)
1 + Q2
1 + Q2
I mL
= 0.505 4 p 3
I SC
Q cos3 L
(3.67)
116
p
3
1 + Q2
p
Q
1 + Q2
(3.68)
117
2.8
2.6
Q = L/R = 2.0
2.4
t (ms)
2.2
(a)
A
Q = L/R = 0.1
2
1.8
1.6
1.4
1.2
0.2
0.4
0.6
0.8
1.2
1.4
1.6
p (%) = 100P/(UISC)
0.035
0.03
Imax/ISC
0.025
0.02
(b)
Q = L/R = 0.1
Q = L/R = 2.0
0.015
0.01
0.005
0.2
0.4
0.6
0.8
1.2
1.4
1.6
p (%) = 100P/(UISC)
1
0.9
p (%) = 0.2
B
0.8
0.7
p (%) = 1.5
ITh/I1
0.6
(c)
0.5
p (%) = 0.6
0.4
0.3
0.2
0.1
0
1
p (%) = 0.85
11
Armnico (h)
118
119
En este punto, el proceso se repite a partir del paso 2 hasta que se alcanza su
convergencia, esto es, la diferencia entre los dos ltimos valores de IDC calculados difiere
en una cantidad menor a una tolerancia fijada. En el presente ejemplo, despus de tres
iteraciones, la solucin a la que se llega es IDC = 9.081 A y t = 2.092 ms. Una vez
obtenido el ancho del pulso, se pueden obtener el valor mximo de la corriente (3.50) y sus
componentes armnicos (3.52).
La caracterizacin de la corriente tambin se podra realizar grficamente a partir la
Figura 3.31a. Para ello se calcula el valor de la variable normalizada p,
p = 10030.15000/4002 = 0.540%, y se accede al grfico correspondiente al valor de
Q = 2, resultando el punto A marcado en la Figura, el cual indica la solucin para el ancho
de pulso correspondiente. De forma anloga, de los grficos mostrados en las Figuras
3.31b y 3.31c se obtienen el resto de valores que caracterizan la corriente que consume el
sistema.
3.2.4.2 Ejemplo 2: Diseo del sistema: clculo de la impedancia del lado de alterna
En este segundo ejemplo, el problema que se resuelve consiste en calcular la
impedancia que debera tener el sistema para obtener un determinado contenido de
armnicos. Los datos del sistema son los siguientes: U = 415 V (f = 50 Hz) y P = 1 kW. Se
debe calcular la impedancia Z necesaria (tomando como hiptesis que Q = 2) para limitar
el quinto armnico segn IT5/IT1 = 0.8.
En este caso, la resolucin se realiza grficamente tal y como se detalla a continuacin:
1. El valor de p correspondiente a IT5/IT1 = 0.8 se obtiene a partir de la Figura 3.31c.
Como se muestra en ella, el punto B corresponde al caso estudiado, obtenindose
un valor de p aproximadamente igual al 0.6%.
2. El valor de la impedancia del lado de alterna se obtiene a partir de la definicin
del ratio de potencia p, Z = 0.64152/(10031000) = 0.569 .
120
1.14
0.77
11
0.33
13
0.21
17 h 37
0.1515/h
En dicha tabla slo se han considerado los armnicos impares h = 6n 1 con n = 1, 2, ...
que son los emitidos por el rectificador en condiciones equilibradas.
Los equipos trifsicos, entre los que estaran las fuentes de alimentacin trifsicas, de
tensin nominal mxima igual a 600 V, y cuyo consumo sea superior a 16 A son tratados
en la Norma Internacional CEI 1000-3-4 [75]. En la norma los equipos se clasifican en tres
estadios:
- Estadio 1: Equipos con conexiones usuales cuyo ratio SSC/SEQ 33. Segn la
norma, los armnicos de la corriente de entrada no debern exceder las relaciones
Ih/I1 dadas en la Tabla 3.3. Estos lmites suelen imponer serias restricciones de
funcionamiento para muchos equipos.
- Estadio 2: Equipos que no cumplan los lmites del estadio anterior cuyo ratio
SSC/SEQ > 33. As, para estos equipos se aceptan valores de Ih/I1 mayores
121
ih (%) = 100Ih/I1
10.7
7.2
11
3.1
13
17
1.2
19
1.1
Tabla 3.3. Lmites armnicos para equipos de corriente superior a 16 A (Estadio 1).
Mnima relacin
SSC/SEQ
ih (%) = 100Ih/I1
i5
i7
i11
i13
66
14
11
10
120
16
12
11
175
20
14
12
250
30
18
13
350
40
25
15
10
450
50
35
20
15
600
60
40
25
18
122
123
de
la
corriente
consumida
por
el
dispositivo,
de
la
corriente
consumida
por
el
dispositivo,
124
Captulo 4
Formulacin del flujo armnico de
cargas
La aportacin de este captulo supone una completa recopilacin, revisin y ordenacin
de los mltiples enfoques de la formulacin del flujo armnico de cargas, en condiciones
equilibradas y desequilibradas, que han ido apareciendo histrica y de forma relativamente
desordenada en la bibliografa. El estudio pretende clasificar las numerosas formas de
tratar el problema armnico en el dominio de la frecuencia y es complementado con el
balance, en condiciones equilibradas, de los datos, incgnitas y ecuaciones planteadas en
cada formulacin existente [46, 47]. Asimismo, para finalizar el captulo se analizan, a
partir de dos ejemplos, las dos hiptesis ms importantes en el enfoque del problema
armnico,
125
4.1 Introduccin
En los estudios de los sistemas elctricos de potencia (EPS, Electric Power Systems) es
necesaria, para un estado de carga dado, la determinacin de las tensiones fundamentales
de la red en rgimen permanente. Esta determinacin, en general se realiza mediante la
resolucin del flujo de cargas convencional.
En el transcurso de los ltimos aos, ha existido un creciente inters en la obtencin, en
rgimen permanente, de las tensiones de red a las frecuencias fundamental y armnicas
debido al incremento del nmero de dispositivos no lineales en los sistemas elctricos de
potencia, lo que ha conducido a la modificacin del flujo de cargas convencional para
contemplar adems del clculo de las tensiones fundamentales el de las armnicas y el
comportamiento de los dispositivos no lineales presentes en la red. Este nuevo enfoque del
problema se denomina genricamente flujo armnico de cargas.
Por otro lado, en un sistema elctrico de potencia, las tensiones trifsicas en los nudos
de la red se espera que sean simtricas de forma que el estudio habitual de la distorsin
armnica del sistema se realiza mediante el flujo armnico de cargas en condiciones
equilibradas. Sin embargo, ciertas asimetras originadas por la no transposicin de las
lneas de transmisin, cargas trifsicas desiguales, cargas monofsicas elevadas, etc.
pueden producir cierto grado de distorsin en la simetra de las tensiones de lnea. Esto,
unido al hecho de que, en sistemas desequilibrados, la generacin y propagacin de
armnicos es un poco ms compleja, ya que, por ejemplo, se pueden producir armnicos
126
NUDOS
TIPO
DATOS
INCGNITAS
N INCGNITAS
SLACK
Vi 1
----
2ag
PV
Pi, Ui
Vi 1
2(g - 1)
g+1ac
PQ
Pi, Qi
Vi 1
2(c - g)
Tabla 4.1. Datos e incgnitas del flujo de cargas convencional (Slack: Nudo de referencia,
PV: Nudos de generadores PV, PQ: Nudos de cargas lineales PQ)
127
( i = 2,
... g )
(4.1)
n 1 1
S i = V Y ij V j
j =1
1
i
( i = g + 1,
... c )
(4.2)
128
a partir de los datos que fijan el funcionamiento de las cargas no lineales y de unas
condiciones de demanda y carga dadas.
Si suponemos un sistema con n nudos, donde hay cargas no lineales, NLD, y
consideramos hm armnicos, la Tabla 4.2 ilustrara las incgnitas del problema a resolver,
donde h = 3, 5, ... y k = 1, 3, 5, ... son los rdenes de los armnicos, y Bi j y i j con j = 1, ...
r, son los datos y parmetros de los dispositivos no lineales, respectivamente.
NUDOS
TIPO
DATOS
INCGNITAS
N INCGNITAS
SLACK
V1 1
Vi h
2 hm
2ag
PV
Pi, Ui
Vi k
2hm
g +1 a c
PQ
Pi, Qi
V i k, Y i 1
2(hm+1)(cg)+2(cg)
c +1 a n
NLD
2(hm+1)(nc)+r(nc)
Tabla 4.2. Datos e incgnitas del flujo armnico de cargas (Slack: Nudo de referencia, PV: Nudos
de generadores PV, PQ: Nudos de cargas lineales PQ, NLD: Nudos de cargas no lineales)
Las tcnicas para analizar el problema armnico se clasifican en dos grandes tipos:
mtodos en el dominio del tiempo y mtodos en el dominio de la frecuencia. Tambin
existen mtodos hbridos que trabajan en ambos dominios aprovechando las ventajas
propias de cada uno de ellos [53, 67, 72]. En estas metodologas, la parte lineal del sistema
se representa en el dominio de la frecuencia y los componentes no lineales en el dominio
del tiempo.
Los mtodos en el dominio del tiempo, como el EMTP [4], se basan en la resolucin,
por integracin numrica, de las ecuaciones diferenciales que representan el
comportamiento del sistema. Estos mtodos son muy tiles en el tratamiento de los
dispositivos no lineales pero el esfuerzo de clculo para alcanzar la solucin de rgimen
permanente puede ser considerable y, adems, no permiten definir el problema en trminos
de consumo de potencia.
129
Los mtodos en el dominio de la frecuencia son una extensin natural del flujo de
cargas convencional, donde se incluye el tratamiento de los dispositivos no lineales y de
las tensiones armnicas. Estos mtodos, denominados genricamente flujo armnico de
cargas (HLF, Harmonic Load Flow) son los ms extensamente utilizados en la resolucin
de problemas con armnicos [1, 2, 6].
En principio, y tal como se realiza en el CLF, los mtodos en el dominio de la
frecuencia requieren el planteamiento de un sistema no lineal de tantas ecuaciones como
incgnitas presenta el problema, Tabla 4.2. Esto lleva a tener que programar un sistema de
ecuaciones de difcil resolucin, debido al gran nmero de incgnitas a tratar, que puede
comportar ciertos inconvenientes como son unos largos tiempos de ejecucin en el clculo
de la solucin, problemas de convergencia numrica y problemas de capacidad del
computador. Es por esta razn que aparecen en la bibliografa mltiples formas de enfocar
el problema armnico en el dominio de la frecuencia en una bsqueda de un compromiso
entre la simplificacin en su tratamiento y la fiabilidad de los resultados finales.
Una vez se ha caracterizado el estado del sistema, incluidos los dispositivos no lineales
presentes en el mismo, las tensiones armnicas pueden ser calculadas mediante el mtodo
de los nudos (VN, Voltage Nodes method) a partir de la matriz de admitancias que
caracteriza el conjunto red-cargas pasivas, YBmh, y de la inyeccin de las intensidades
armnicas de los dispositivos no lineales. El proceso queda reflejado en la Figura 4.1.
CLFm
VN
(YBmh )
Los datos e incgnitas del mtodo de Penetracin Armnica para un sistema de n nudos
VN
CLFm
NUDOS
TIPO
DATOS
INCGNITAS
N INCGNITAS
SLACK
V1 1
----
2ag
PV
Pi, Ui
Vi 1
2(g-1)
g+1 a c
PQ
Pi, Qi
Vi 1
2(c-g)
c+1 a n
NLD
(2+r)(n-c)
SLACK
X1
V1 h
2hm
2ag
PV
Xi
Vi h
2 hm(g-1)
g+1 a c
PQ
Yi h
Vi h
2 hm(c-g)
c+1 a n
NLD
Ii h
Vi h
2 hm(n-c)
I i = d i V i , i1 ,... ir
1
( i = c + 1,... n )
(4.3)
Bi j = nl j V i , i1 ,... ir
1
( i = c + 1,...n ;
j = 1,...r )
(4.4)
Resuelto el sistema, las corrientes armnicas inyectadas en la red por los dispositivos no
lineales, Iih, con h = 3,5,..., se pueden obtener a partir de funciones del tipo (4.3).
Finalmente, la consideracin de la potencia nicamente debida a la componente
fundamental, permite calcular las admitancias de las cargas PQ directamente a partir de la
solucin del flujo de cargas convencional modificado:
1
i
Y =
Si
(V )
1 2
; Yi = f i Yi,h
( i = g + 1,... c ;
h = 3,5,...)
(4.5)
donde fi es una funcin que depende del modelo de la admitancia armnica de la carga [ 45,
51, 71].
132
Es importante volver a remarcar que las expresiones (4.3, 4.4 y 4.5) slo dependen de la
tensin fundamental de los nudos debido a la hiptesis de que las tensiones armnicas no
influyen en el comportamiento de las cargas no lineales ni en la potencia consumida.
El mtodo VN se basa en la resolucin del sistema lineal correspondiente a la expresin
matricial IBh = YBmhUBh para las frecuencias armnicas, h = 3, 5,... Dicho mtodo se
completa incorporando las corrientes inyectadas por los dispositivos no lineales al vector
de corrientes de nudo IBh , y las admitancias armnicas de las cargas y los generadores a la
matriz de admitancias de red YBh. Todo ello calculado a partir de los resultados del CLFm y
gracias a las hiptesis realizadas y ya comentadas.
Esta formulacin permite abordar el problema armnico de una manera tan sencilla
como la del flujo de cargas convencional pero el hecho de no considerar la sensibilidad de
los dispositivos no lineales frente a las tensiones armnicas se traduce en una
sobreestimacin de la contaminacin armnica en la red [34, 42, 47] y slo puede ser
aplicado cuando la distorsin sea reducida.
Por esta razn, la necesidad de resolver el problema de una manera ms precisa ha
llevado a considerar la interaccin entre la red y los dispositivos no lineales y a desarrollar
un nuevo conjunto de formulaciones. Estas metodologas todava conservan la hiptesis de
que la potencia solamente se define a partir de la onda fundamental de la tensin y
corriente.
Actualmente, la mayor parte de las formulaciones del flujo armnico de cargas se basan
en ambas consideraciones, es decir,
- Interaccin armnica: las expresiones (4.3 y 4.4) que caracterizan el
comportamiento de los dispositivos no lineales y las corrientes fundamental y
armnicas inyectadas k = 1, 3, 5 ..., se deben modificar para que reflejen su
dependencia de las tensiones armnicas:
133
Bi j = nl j V i , V i , i1 ,... ir
k
Ii = d
( V ,V
1
i
h
i
, i1 ,... ir
( i = c + 1,...n ;
j = 1,...r ; h = 3,5,...)
( i = c + 1,... n ;
h = 3,5,...)
(4.6)
(4.7)
- Potencia fundamental: las admitancias de las cargas pueden ser calculadas a partir
de la expresin (4.5). Esto permite mantener la idea original del mtodo HP, es
decir, desacoplar la resolucin del flujo de cargas convencional y la determinacin
de las tensiones armnicas.
CLFm
IHA
(YBRh )
VN
(YBmh )
134
VN
IHA
CLFm
NUDOS
TIPO
DATOS
INCGNITAS
N INCOGNITAS
SLACK
V1 1
----
2ag
PV
Pi, Ui
Vi 1
2(g-1)
g+1 a c
PQ
Pi, Qi
Vi 1
2(c-g)
c+1 a n
NLD
(2+r)(n-c)
SLACK
X1
----
2ag
PV
Xi
----
g+1 a c
PQ
Yi h
----
c+1 a n
NLD
(2hm+r)(n-c)
SLACK
X1
V1 h
2hm
2ag
PV
Xi
Vi h
2hm(g-1)
g+1 a c
PQ
Yi h
Vi h
2hm(c-g)
c+1 a n
NLD
Ii h
----
135
As, una vez obtenida la solucin del CLFm, los generadores de los nudos PV se
sustituyen por sus reactancias, y las cargas correspondientes a los nudos PQ por sus
impedancias, (4.5), del mismo modo que en la penetracin armnica. La caracterizacin
pasiva del comportamiento armnico de los generadores PV y de las cargas PQ se puede
realizar gracias a la hiptesis de potencia fundamental, incorporndose sus admitancias
armnicas a la matriz de admitancias de la red YBh, para formar la matriz YBmh.
Posteriormente, el sistema se reduce a los nudos no lineales, Figura 4.3, mediante la
eliminacin Gaussiana de la matriz de admitancias modificada correspondiente a la parte
lineal de la red, obtenindose la matriz YBRh. Este procedimiento permite reducir el estudio
a los nudos donde se encuentran conectados los dispositivos no lineales, i = c+1, ... n.
g
2
.....
g+1
hX g
Yg+1h
.....
hX2
1
Red
YBh
hX 1
c
Y ch
c+1
.....
n
Red pasiva
Y Bmh
c+1
.....
Y BRh
Ic+1h = d h(Vc+1h)
n
Inh = d h(Vnh)
136
(I )
h
i
( ) ,( )
h
= d Vi, Vi
(( V )
h
i
(4.8)
Las corrientes calculadas permiten el clculo de las tensiones armnicas de los nudos no
lineales para la siguiente iteracin:
V
....
Vh
n
h
c +1
+1
= (Y
h 1
BR
....
Ih
n
h
c +1
d h V ch+1
h 1
= ( YBR
) ....
h
h
d Vn
( )
(4.9)
137
138
CLF
HA
(Ethk, Zthk )
no
(Vi1)CLF = (Vi1)HA ?
s
VN
(YBmh )
Figura 4.4. Estructura del mtodo del flujo de cargas armnico simplificado
139
VN
HA
CLF
NUDOS
TIPO
DATOS
INCGNITAS
N INCGNITAS
SLACK
V1 1
----
2ag
PV
Pi, Ui
Vi 1
2(g-1)
g+1 a c
PQ
Pi, Qi
Vi 1
2(c-g)
c+1 a n
NLD
Ii 1
Vi 1
2(n-c)
SLACK
X1 , V1 1
----
2ag
PV
Xi , Vi 1
----
g+1 a c
PQ
Yi h
----
c+1 a n
NLD
(2(hm+1)+r)(n-c)
SLACK
X1
V1 h
2hm
2ag
PV
Xi
Vi h
2hm(g-1)
g+1 a c
PQ
Yi h
Vi h
2hm(c-g)
c+1 a n
NLD
Ii h
----
Ethk = ; Zthk =
k
k
k
Z th n c , 1 Z th n c , n c
E th n c
(4.10)
140
Red pasiva
Y Bm1
g
.....
Vg
g+1
Yg+11
Red
Y B1
c+1
c
Y c1
.....
.....
V2 1
Eth1
Zth1
V1 1
c+1
.....
g
2
g+1
.....
hXg
Yg+1h
.....
hX 2
hX 1
Red
Y Bh
c
Y ch
c+1
Ethh = 0
Zthh
c+1
.....
.....
n
Bn j = nl j(Vnk, n1, ... nr)
Inh = d h(Vnk, n1, ... nr)
Red pasiva
Y Bm h
Las tensiones Thvenin se obtienen calculando las tensiones de los nudos no lineales sin
conectar sus cargas (tensiones a circuito abierto). Ntese que las tensiones Thvenin son
igual a cero para k 1 ya que no existen generadores de armnicos en la red.
La determinacin de la matriz de impedancias Thvenin se calcula mediante el mtodo
de la inyeccin de corriente. El procedimiento consiste en inyectar una intensidad de valor
141
unidad por uno de los nudos no lineales y calcular las tensiones resultantes en cada uno de
los nudos no lineales. Las tensiones obtenidas son las impedancias Thvenin entre el nudo
en el que se inyecta la corriente y el resto. La operacin se repite para todos los nudos no
lineales hasta completar la matriz de impedancias Thvenin. Este proceso se realiza para
cada uno de los armnicos considerados, obtenindose las diferentes Zthk (k = 1, 3, 5, ...).
Una vez obtenido el equivalente Thvenin, el sistema de ecuaciones no lineales
planteado para resolver el HA est formado por:
- Las restricciones impuestas por el equivalente Thvenin,
nc
V i = E thi Z thi , j I i
k
( i = 1,...n c ;
k = 1,3,5,...)
(4.11)
j =1
142
143
Esta ltima formulacin permite, al no tener que reducirse el sistema a los nudos no
lineales, considerar la potencia debida nicamente a la frecuencia fundamental, (4.5), o
debida a las frecuencias fundamental y armnicas,
Si =
k =1, 3, ....
( )
k *
Vi Ii
( i = g + 1,...c )
(4.12)
NUDOS
TIPO
DATOS
INCGNITAS
N INCGNITAS
SLACK
V1 1
V1 h
2hm
2ag
PV
Pi, Ui
Vi k
2(hm+1)(g-1)
g+1 a c
PQ
Pi, Qi
Vi k
2( hm+1)(c-g)
c+1 a n
NLD
(2( hm+1)+r)(n-c)
144
(V ) ( i = 2,
k 2
k =1, 3, ...
... g )
(4.13)
n 1 1
S i = V Y ij V j
j =1
1
i
( i = g + 1,
... c )
(4.14)
Y i V i + Y ij V j = 0
k
( i = 1,
j =1
... c ; k = 1,3,5,..., hm )
(4.15)
donde Yik son las admitancias armnicas de los generadores y de las cargas PQ.
Las admitancias de los generadores son dato y las de las cargas se pueden calcular
a partir de (4.5) gracias a la consideracin de potencia fundamental.
- El balance de corrientes armnicas en los nudos no lineales:
n
Y
j =1
k
ij
V j = Ii
( i = c + 1,
...n ; k = 1,3,5,..., hm )
(4.16)
NUDOS
TIPO
DATOS
INCGNITAS
N INCGNITAS
SLACK
V1 1
V1 h
2hm
2ag
PV
Pi, Ui
Vi k
2(hm+1)(g-1)
g+1 a c
PQ
Pi, Qi
V i k, Y i 1
(2(hm+1)+2)(c-g)
c+1 a n
NLD
(2(hm+1)+r)(n-c)
k n k k
Pi = Re V i Y ij V j ;
k =1, 3 ...
j =1
Ui =
(V ) ( i = 2,
k 2
k =1, 3, ...
(4.17)
... g )
n k k
S i = V Y ij V j
k =1, 3 ...
j =1
k
i
( i = g + 1,
... c )
(4.18)
(V ) (Y ) ( i = g + 1,... c )
k 2
k =1, 3 ...
k *
i
(4.19)
donde Yik = fi(Yi1, k). Estas relaciones deben ser incluidas en el sistema de
ecuaciones no lineal para obtener el valor de las admitancias fundamentales.
La resolucin numrica de los sistemas de ecuaciones no lineales de las formulaciones
CHLFF y CHLFH proporciona el valor de las incgnitas del problema, siendo el
procedimiento utilizado ms comnmente en la bibliografa el mtodo de NewtonRaphson.
En un paso ms hacia delante en la mejora de esta ltima metodologa, se han propuesto
otras formulaciones completas basadas en la utilizacin de equivalentes Norton en la
modelizacin de los dispositivos no lineales y en la resolucin por separado del mtodo de
Newton-Raphson para cada armnico [62, 63].
Recientemente se ha presentado un estudio que compara y discute acerca de las
propiedades de convergencia del mtodo desacoplado del SHLF y del mtodo de Newton
unificado del CHLFF [55]. De acuerdo con este estudio, el SHLF puede presentar
problemas de convergencia debido al algoritmo de Gauss-Seidel que enlaza el CLF y el
HA, mientras que el CHLFF es ms rpido y con una convergencia ms robusta.
147
V 11= V 11 1
X1
ZT=jXT
Rectificador
Trifsico
No Controlado
R L, X L, X C , R
Carga PQ
P3, Q3
Z23=R23+jX23
U 2 , P2
X2
Figura 4.6. Sistema a estudio
i4A
RL
XL
+
i4B
i4C
RL
XL
RL
XL
uC
v4
Figura 4.7. Esquema del rectificador trifsico no controlado
Los datos para este dispositivo no lineal son las resistencias y reactancias de alterna, RL,
XL, el condensador de filtro y la resistencia de carga.
148
Los parmetros que definen el comportamiento del dispositivo son nicamente los
instantes de tiempo de inicio y final de la intensidad de alterna i, los instantes t1 y t2 en la
Figura 4.8. Ello es debido a la consideracin de simetra en el modelo utilizado para el
rectificador [66].
uC
t
i4 ( iA)
t1
t2
t1+/3
t2 +/3
(
nl (V
)
, t , t ) = i (t ) = 0
1
4
(4.20)
(
nl (V
)
, t , t ) = i (t ) = 0
1
4
,V
h
4
(4.21)
149
Asimismo, una vez caracterizados sus parmetros, la fuente de alimentacin podr ser
modelada por una fuente de corriente de valor igual a la corriente armnica que inyecta a la
red. Esta corriente puede ser calculada bajo la hiptesis de no interaccin,
( k = 1,5, 7,...)
(4.22)
( k = 1, 5, 7,... ;
h = 5, 7,...)
(4.23)
k
k
0
Y
Y
0
22
23
YBk = k
k
k
k
Y 31 Y 32 Y 33 Y 34
k
k
0
0 Y 43 Y 44
( k = 1,5, 7,)
(4.24)
siendo los elementos que forman parte de ella los expuestos a continuacin:
k
Y 11 =
k
Y 33 =
k
1
1
1
1
k
=
; Y 22 = k =
k
Z 13 R13 + jkX 13
Z 23 R23 + jkX 23
1
1
1
1
1
1
+ k + k =
+
+
k
Z 13 Z 23 Z T R13 + jkX 13 R23 + jkX 23 jkX T
k
Y 13 = Y 31 =
(4.25)
1
1
1
k
k
k
k
; Y 23 = Y 32 =
; Y 34 = Y 43 =
R13 + jkX 13
R23 + jkX 23
jkX T
150
DATOS
INCGNITAS
N INCGNITAS
V1 1
----
P2, U2
V2 1
P3, Q3
V3 1
V 4 1 , t1 , t2
X1
V1 h
2hm
X2
V2 h
2 hm
Y3 h
V3 h
2 hm
I4 h
V4 h
2 hm
VN
CLFm
NUDOS
{
Re {V (Y V + Y
Im {V (Y V + Y
1
1
3
1
31
1
1
1
32
1
3
1
31
1
1
1
32
(
nl (V
)
,t ,t ) = 0
) } P = 0
1 *
Re V 2 Y 22 V 2 + Y 23V 3
1
1
2
1
33
1
3
1
34
V +Y V +Y V
) } P = 0
1 *
V 2 + Y 33V 3 + Y 34 V 4
) } Q
1 *
4
(4.26)
=0
nl1 V 4 , t1 , t2 = 0
2
1
4
que responde a las restricciones impuestas por los nudos Slack, PV y PQ y las ecuaciones
que definen el comportamiento de la fuente de alimentacin.
A continuacin se obtienen las tensiones armnicas de los nudos del sistema aplicando
el mtodo de los nudos. La hiptesis de potencia debida nicamente a la componente
correspondiente a la frecuencia fundamental permite pasar del esquema de la Figura 4.6 al
151
1
Z13h= R13+jhX13
3
jhX1
ZTh=jhXT
I 4h
2
Z23h=R23+jhX2
Y 3h
jhX2
de la Figura 4.9, ya que las admitancias armnicas que modelan la carga PQ pueden ser
calculadas segn la siguiente expresin,
1
Y =
R3 + jhX 3
h
3
1
3
con Z = R3 + jX 3 =
(V )
1 2
3
P3 jQ3
(4.27)
h
I hB = YBm
VBh
h
h
h
0
0 V 1
Y 13m
0 Y 11m
h
0
h
h
0
Y
Y
0
V 2
22 m
23 m
= h
h
h
h
h
0 Y
Y 32 m Y 33m Y 34 m V 3
31m
h
h
h
h
I 4 0
0
Y 43m Y 44 m V 4
(4.28)
con,
h
Y 11m = Y 11 +
152
1
1
h
h
h
h
h
; Y 22 m = Y 22 +
; Y 33m = Y 33 + Y 3
jhX 1
jhX 2
(4.29)
VN
IHA
CLFm
NUDOS
DATOS
INCGNITAS
N INCGNITAS
V1 1
----
P2, U2
V2 1
P3, Q3
V3 1
V 4 1 , t1 , t2
X1
----
X2
----
Y3 h
----
V 4 h , t1 , t2
2hm+2
X1
V1 h
2 hm
X2
V2 h
2 hm
Y3 h
V3 h
2 hm
I4 h
----
Tabla 4.9. Datos e incgnitas del ejemplo para la penetracin armnica iterativa
153
Y44Rh
I 4h
h
h
h
h
h
0 Y h
0 0
0
Y 32 m Y 33m Y 34 m V 3
Y 33mR
31m
h
h
h
h
h
I 4 0
I 4 0
0
0
0
Y 43m Y 44 m V 4
h
0 V 1
h
0 V 2
h
h
Y 34 m V 3
h
h
Y 44 mR V 4
(4.30)
En este caso, el sistema slo consta de un nudo no lineal, nudo 4, as que la expresin
anterior se puede escribir como,
h
I 4 = Y 44 R V 4
(4.31)
y es la que se resuelve mediante un mtodo iterativo para cada armnico que se desea
calcular. Esto quiere decir que, asumiendo constantes las tensiones fundamentales en cada
interaccin , los ltimos valores de las tensiones armnicas (V4h) son utilizados para
actualizar el estado de la fuente de alimentacin, (t1) y (t2), y las corrientes inyectadas en
la red (I4h) = d
(V )
h +1
4
1
Y
h
44 R
I4 =
h
44 R
((V ) )
h
4
(4.32)
El tercer y ltimo paso, consiste en la aplicacin del mtodo de los nudos, tal y como se
explic en el mtodo de la penetracin armnica.
154
VN
HA
CLF
NUDOS
DATOS
INCGNITAS
N INCGNITAS
V1 1
----
P2, U2
V2 1
P3, Q3
V3 1
I4 1
V4 1
X1 , V1 1
----
X2 , V2 1
----
V3 1
----
V 1 k, t 1 , t 2
2(hm+1)+2
X1
V1 h
2hm
X2
V2 h
2hm
Y3 1
V3 h
2hm
I4 h
----
Tabla 4.10. Datos e incgnitas del ejemplo para el flujo de cargas armnico simplificado
En primer lugar se debe resolver el proceso iterativo establecido entre las formulaciones
CLF y HA para posteriormente abordar el proceso VN.
El flujo de cargas convencional, CLF, resolver el siguiente sistema de ecuaciones no
lineales:
155
V21 U 2 = 0
{
}
Re {V (Y V + Y V + Y V + Y V ) } P = 0
Im {V (Y V + Y V + Y V + Y V ) } Q = 0
Re {V (Y V + Y V ) } Re {V ( I ) } = 0
Im {V (Y V + Y V ) } Im {V ( I ) } = 0
Re V 2 (Y221 V21 + Y231 V31 ) P2 = 0
*
1
3
1 1
31 1
1 1
32 2
1 1
33 3
1 1 *
34 4
1
3
1 1
31 1
1 1
32 2
1 1
33 3
1 1 *
34 4
1
4
1 1
43 3
1 1 *
44 4
1
4
1 *
4
1
4
1 1
43 3
1 1 *
44 4
1
4
1 *
4
(4.33)
I 4 = d 1 V 4 ,V 4 , t1 , t2
(4.34)
V1
Z131
ZT1
Eth1
Zth1
Z 31
Z231
V21
Figura 4.11. Sistema a estudio para la frecuencia fundamental
1
Z13h
jhX1
Z23h
Z Th
Ethh = 0
Zthh
Z 3h
jhX2
Una vez obtenidos los equivalentes Thvenin de la red vistos desde los nudos no
lineales, el estudio del problema se reduce al planteamiento de las ecuaciones
correspondientes a los esquemas que se muestran en las Figuras 4.13 y 4.14, junto con las
ecuaciones que caracterizan el comportamiento de los dispositivos no lineales,
157
E th Z th I 4 V 4 = 0
h
Z th I 4 + V 4 = 0
(
nl (V
)
,t ,t ) = 0
(4.35)
nl1 V 4 , V4h , t1 , t2 = 0
2
1
4
, V4h
donde
h
I 4 = d h V 4 , V 4 , t1 , t2
(4.36)
4
I41
Zth1
+
Eth1
V41
4
Zthh
I4h
+
V4h
158
Por ltimo, una vez obtenidas las tensiones fundamental y armnicas de los nudos no
lineales, y los parmetros de los dispositivos no lineales, se pueden calcular las corrientes
armnicas inyectadas por stos a la red y, mediante el mtodo de los nudos, calcular las
tensiones armnicas de los nudos lineales del sistema.
NUDOS
DATOS
INCGNITAS
N INCGNITAS
V1 1
V1 h
2hm
P2, U2
V2 k
2( hm+1)
P3, Q3
V3 k
2( hm+1)
V 4 k, t 1 , t 2
2( hm+1)+2
Tabla 4.11. Datos e incgnitas del ejemplo para el flujo armnico de cargas completo
(componente fundamental)
{ (
Re {V (Y V + Y
Im {V (Y V + Y
1
) } P = 0
1 *
Re V 2 Y 22 V 2 + Y 23V 3
1
1
31
1
1
1
32
V 2 + Y 33V 3 + Y 34 V 4
1
3
1
31
1
1
1
32
V 2 + Y 33V 3 + Y 34 V 4
V1
h
h
h
h
+ Y 11V 1 + Y 13V 3 = 0
jhX 1
h
V2
h
h
h
h
+ Y 22 V 2 + Y 23V 3 = 0
jhX 2
) } P = 0
1 *
1
3
) } Q
1 *
=0
(4.37)
( h = 3,5, , hm )
( h = 3,5, , hm )
159
( h = 3,5, , hm )
Y 3 V 3 + Y 31V 1 + Y 32 V 2 + Y 33V 3 + Y 34 V 4 = 0
h
( h = 3,5, , hm )
I 4 + Y 43V 4 + Y 44V 4 = 0
(
nl (V
)
,t ,t ) = 0
(4.37 cont.)
nl1 V 4 , V4h , t1 , t2 = 0
2
1
4
, V4h
con
1
3
Z = R3 + jX 3 =
(V )
1 2
3
(4.38)
P3 jQ3
I 4 = d h V 4 , V 4 , t1 , t2
(4.39)
NUDOS
DATOS
INCGNITAS
N INCGNITAS
V1 1
V1 1
2hm
P2, U2
V2 1
2( hm+1)
P3, Q3
V3 1 , Y 3 1
2( hm+1)+2
V 1 1 , t1 , t2
2( hm+1)+2
Tabla 4.12. Datos e incgnitas del ejemplo para el flujo armnico de cargas completo
(componentes fundamental y armnicas)
160
{ (
Re {V (Y V + Y
Im {V (Y V + Y
1
) } P = 0
1 *
Re V 2 Y 22 V 2 + Y 23V 3
1
1
31
1
1
1
32
V 2 + Y 33V 3 + Y 34 V 4
1
3
1
31
1
1
1
32
V 2 + Y 33V 3 + Y 34 V 4
V1
h
h
h
h
+ Y 11V 1 + Y 13V 3 = 0
jhX 1
h
) } P = 0
3
) } Q
1 *
(
nl (V
)
,V , t , t ) = 0
(V ) (Y )
=0
( h = 3,5, , hm )
1
Y 3 V 3 + Y 31V 1 + Y 32 V 2 + Y 33V 3 + Y 34 V 4 = 0
I 4 + Y 43V 4 + Y 44V 4 = 0
( h = 3,5, , hm )
V2
h
h
h
h
+ Y 22 V 2 + Y 23V 3 = 0
jhX 2
h
1 *
1
3
( h = 3,5, , hm )
( h = 3,5, , hm )
nl1 V 4 , V4h , t1 , t2 = 0
2
S3 =
1
4
h
4
k =1,..., hm
k 2
3
k *
3
(4.40)
conduce a la obtencin de la solucin del problema armnico. Se observa que en este caso
las admitancias que modelan la carga PQ pasan a ser una incgnita ms del problema con
lo que se debe aadir una ecuacin ms correspondiente al balance de potencias en el nudo
PQ.
161
162
Convertidor
6p AC/DC
(NL1)
V11
= 10
X1 = 0.102
SN =8 MVA
UN1 = 25 KV
UN2 = 0.4 KV
cc = 5%
Wcc = 120 kW
NL1
NL2
Rectificador
Trifsico
No Controlado
(NL2)
RL = 0.01
XL = 0.1
Datos
XNL1 = 0.848
= 0
ID = 0.692
RNL2 = 0.5
XNL2 = 0.314
XC = 0.375
R = 178.75
PQ
P3 = 0.0625
Q3 = 0.0375
Parmetros
t0,
t1, t2
Los parmetros o incgnitas del convertidor de seis pulsos son los pasos por cero de la
tensin de lnea, t0, y el ngulo de retardo o solape en el establecimiento de la corriente de
lnea del lado de alterna, . stas se obtienen a partir de las condiciones impuestas, por un
lado, por los pasos por cero de las tensiones de lnea y, por otro lado, por el instante en que
la corriente del lado de alterna alcanza su valor ID [68].
Los parmetros o incgnitas del rectificador trifsico no controlado son los tiempos de
conmutacin de los diodos del puente rectificador, t1 y t2, que se obtienen a partir de las
condiciones de cambio de la conduccin del rectificador, por un lado la que impone el
inicio de la conduccin, que corresponde al instante en que la tensin de red rectificada y
la tensin del condensador se igualan, y por otro la que impone el final de la conduccin,
que corresponde al paso por cero de la onda de intensidad [34].
163
HP
TIPO
DATOS
INCGNITAS
N INCGNITAS
SLACK
V1 1
----
PQ
P2, Q2 (= 0)
V2 1
PQ
P3, Q3
V3 1
XNL1, ID,
t0 ,
ZNL2, XC, R
t1 , t2
CLFm
NUDOS
VN
NLD
SLACK
X1
V1 h
2hm
PQ
Y2h (= 0)
V2 h
2hm
PQ
Y3 h
V3 h
2hm
3
NLD
INL1 , INL2
CHLFF
NUDOS
TIPO
DATOS
INCGNITAS
N INCGNITAS
SLACK
V1 1
V1 h
2hm
PQ
P2, Q2 (= 0)
V2 k
2hm + 2
PQ
P3, Q3
3
NLD
V3 k
XNL1, ID,
t0 ,
ZNL2, XC, R
t1 , t2
2hm + 6
CHLFH
NUDOS
TIPOS
DATOS
INCGNITAS
N INCGNITAS
SLACK
V1 1
V1 h
2hm
PQ
P2, Q2 (= 0)
V2 k
2hm + 2
PQ
P3, Q3
3
NLD
V 3 k, Y 3 1
XNL1, ID,
t0 ,
ZNL2, XC, R
t1 , t2
2hm + 8
En ambos casos, las condiciones que permiten la obtencin de las incgnitas de los
dispositivos no lineales dependen, en primer lugar, de los propios parmetros de stos y, en
segundo lugar, de la tensin fundamental del nudo en el caso de la formulacin HP, y de
las tensiones fundamental y armnicas del nudo en los casos de las formulaciones CHLFF
y CHLFH.
As, los datos y las incgnitas del sistema a estudio para las tres formulaciones y
considerando hm armnicos se encuentran en la Tabla 4.13.
El nmero de ecuaciones que se debe plantear en cada una de las formulaciones, para
obtener las correspondientes incgnitas, es,
- HP: un sistema no lineal de ecuaciones con 8 incgnitas y un sistema lineal de
ecuaciones con 6 hm incgnitas.
- CHLFF: un sistema no lineal de ecuaciones con 6 hm + 8 incgnitas.
- CHLFH: un sistema no lineal de ecuaciones con 6 hm + 10 incgnitas.
Se puede apreciar la diferencia entre la formulacin HP y las otras dos formulaciones
debido a la hiptesis de no interaccin armnica. Por otro lado, la posibilidad de calcular la
admitancia
fundamental
de
las
cargas
PQ
de
acuerdo
con
la
expresin
Y31 = (P3jQ3)/(V31)2, permite una reduccin del nmero de incgnitas en el CHLFF con
respecto a CHLFH.
El problema se resuelve con el paquete informtico INTARBCN, el cual permite
analizar redes de hasta 200 nudos y 500 componentes. Dicho programa, desarrollado bajo
entorno Windows, consta de un motor de clculo programado en Fortran 4.0 y una
interfase de comunicacin con el usuario programada en Visual Basic 4.0 [52, 77, 78].
El programa incluye como elementos no lineales los desarrollados en la referencia [66]
y como elementos lineales, nudos Slack, PQ y PV, lneas (modelos longitudinal, en y
parmetros distribuidos), transformadores (de regulacin fija y variable) y equipos de
165
compensacin y filtrado [45, 51, 69, 71]. Sus prestaciones permiten resolver el problema
armnico a partir de las tres formulaciones a estudio, HP, CHLFF y CHLFH pudiendo
abordar la resolucin numrica por el mtodo de Newton o sus variantes [71].
Una vez resuelto el problema mediante las tres formulaciones anteriores, primero se
comparan los resultados obtenidos para las formulaciones HP y CHLFF en las Figuras 4.16
a 4.20.
Las Figuras 4.16 y 4.17 muestran la distorsin armnica y la distorsin armnica total
de las corrientes de los dispositivos no lineales. En estas figuras se muestran adems los
parmetros de las cargas no lineales (pasos por cero de la tensin de conmutacin, t0, y el
ngulo de retardo, , para el convertidor AC/DC, y los pasos por cero de la corriente, t1 y
t2, para el rectificador trifsico no controlado).
16
HDI (%)
12
CHLFF: t0 = -6.573 ms
= 0.713 rad
THDI = 15.38%
10
8
6
HP: t1 = 0.829 ms
t2 = 1.263 ms
THDI = 107.08%
8
7
HDI (%)
14
CHLFF: t1 = 0.932 ms
t2 = 1.820 ms
THDI = 76.67%
6
5
4
3
2
0
HP: t0 = - 6.589 ms
= 0.710 rad
THDI = 15.83%
11
13
11
13
166
3.5
3.0
4.5
4.0
3.5
THD V (%)
HDV (%)
2.5
5.0
2.0
1.5
1.0
CHLFF
3.0
2.5
2.0
1.5
1.0
0.5
0.5
0.0
11
0.0
13
NUDO N1
NUDO N2
NUDO N3
1.5
1.5
0.5
0.5
u3 (p.u.)
u3 (p.u.)
HP
0
-0.5
0
-0.5
-1
-1
-1.5
0
20
40
tiempo (ms)
60
80
-1.5
20
40
60
80
tiempo (ms)
Figura 4.20. Formas de onda de la tensin del nudo 3: HP (izquierda), CHLFF (derecha)
167
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0.0
Carga
PQ
Convertidor Rectificador
6-p AC/DC trifsico no
controlado
Total
168
Captulo 5
Conclusiones y lneas futuras de
investigacin
169
170
171
modelos
analticos
sencillos
que
permitan
determinar
el
172
Apndice
Modelizacin de los rectificadores no
controlados con filtro capacitivo
El estudio del comportamiento del rectificador monofsico y la validacin del modelo
aproximado de los rectificadores monofsico y trifsico del captulo 3 se realizan a partir
de la resolucin exacta de las ecuaciones que describen el comportamiento del dispositivo
[66].
173
i(t)
+
uC(t)
u(t)
IDC
t1 t2
- Tramo 2:
t2 t3
- Tramo 3:
t3 t4
- Tramo 4:
t4 t5
i
t1
t2
(1)
t3
(2)
t4
(3)
t5
(4)
(A.1)
174
por ejemplo, los tramos 1 y 2, y determinar los instantes de tiempo t1 y t2 (t3 = t1+T/2,
t4 = t2 + T/2, t5 = t1+T).
Durante el tramo1 el rectificador no conduce y el condensador se descarga sobre la
carga. El esquema del circuito pasa a ser el mostrado en la Figura A.3.
+
uC(1)(t)
IDC
(2)
Z i (t)
+
u(t)
+
uC(2)(t)
IDC
175
(1)
I DC
t + K1
C
(t ) = 0
(A.2)
A.1.1.1.2 Tramo 2
Para el tramo 2, las expresiones que describen el comportamiento del circuito son:
di ( )
2
u (t ) = L
uC( )
dt
2
( 2)
duC( )
= C
I DC
dt
2
(A.3)
176
Solucin homognea:
El sistema homogneo es
0=L
di (
2)
( t ) u ( 2)
C
dt
(t )
(A.4)
du ( ) ( t )
2
i ( ) ( t ) = C C
dt
2
duC( )
2)
( 2)
( t ) h = C
dt
( 2)
( 2)
= K1 C sin t + K 2
)=K
1p
( 2)
( 2)
K1 sin t + K 2
(A.5)
I(2)
UC(2)
177
( )
= I (pa) Ipa
2
j L
(A.6)
1
2
( 2)
( 2)
( 2)
= I pa( ) j
U Cpa
= U Cpa UC pa
(A.7)
( )
U Cpb
=0
2
Z
+
UC(2)
IDC
Solucin completa:
La solucin del sistema (A.3) es la suma de las soluciones homognea y particular
obtenidas:
)
K ( ) sin ( t + K ( ) ) + 2 I ( ) cos ( t + ( ) ) I
i(
2)
( t ) = K1 p
2)
)+
( )
2U Cpa
cos t + U( C)pa
2
2
pa
2
Ipa
(A.8)
DC
A.1.1.2.1 Tramo 1
Para el tramo 1, las expresiones que describen la intensidad de entrada y la tensin del
condensador son:
uC( ) ( t ) =
1
(1)
I DC
t + K1
C
(A.9)
(t ) = 0
A.1.1.2.2 Tramo 2
Para el tramo 2, las expresiones que describen el comportamiento del circuito son:
u ( t ) = Ri
( 2)
(t ) + L
di (
2)
( t ) u ( 2)
dt
(t )
(A.10)
duC( ) ( t )
i ( t ) = C
I DC
dt
2
( 2)
2)
(t ) + L
du ( ) ( t )
2
i ( ) ( t ) = C C
dt
2
di (
2)
( t ) u ( 2)
dt
(t )
(A.11)
cuya resolucin en este caso conduce a considerar dos casos, segn el signo del
discriminante del polinomio caracterstico, D. Las races del polinomio caracterstico
2 +
R
1
+
=0
L
LC
(A.12)
179
son
(A.13)
donde
=
=
R
2L
1
D
2
(A.14)
4
R
D=
L LC
Si D > 0 (caso aperidico), se tiene 1 = + , 2 = , y
uC( ) ( t ) = K1( ) e1t + K 2( ) e2t
2
duC( )
(A.15)
( 2)
( t ) h = C
dt
( 2)
( 2)
1t
= K1 C 1e K 2 C 2 e
2t
( 2 ) 1t
( 2 ) 2t
= K1 p K1 e + K 2 p K 2 e
duC( )
2)
( 2)
( t ) h = C
dt
= K1( ) e t K1 p cos t + K 2(
2
) + K ( ) C e
sin ( t + K ( ) ) )
= K1( )C e t cos t + K 2(
2)
)+ K
2)
sin t + K 2(
2)
) (A.16)
2p
I pa( ) =
( )
= I (pa) Ipa
2
R + j L
(A.17)
1
2
( 2)
( 2)
( 2)
= I pa( ) j
U Cpa
= U Cpa UC pa
(A.18)
( )
U Cpb
= I DC R
2
( )
uC( ) ( t ) = K1( ) e1t + K 2( ) e2t + 2U Cpa
cos t + U( C)pa I DC R
2
i(
2)
( )
2 I (pa) cos t + Ipa
2
(A.19)
i(
2)
2)
)+
)
sin ( t + K ( ) ) ) +
( )
2U Cpa
cos t + U( C)pa I DC R
2
(A.20)
( )
2 I (pa) cos t + Ipa
I DC
2
181
2)
( t2 ) = i (1) ( t2 ) = 0
(A.21)
(A.22)
(A.23)
(A.24)
Por otro lado, la conduccin termina cuando la onda de intensidad alcanza el valor cero:
T
2
f 2 ( t1 , t2 ) = i ( ) t3 = t1 + = 0
2
(A.25)
182
+
u(t)
i(t)
uC (t)
IDC
(1)
t1
(2)
t2
(3)
t3
(4)
t4
(5)
t5
(6)
t6
(7)
t7
(8)
t8
(9)
t9
t11
t12
t13
- Tramo 1:
t1 t2
- Tramo 2:
t2 t3
- Tramo 3:
t3 t4
- Tramo 4:
t4 t5
- Tramo 5:
t5 t6
- Tramo 6:
t6 t7
- Tramo 7:
t7 t8
- Tramo 8:
t8 t9
- Tramo 9:
t9 t10
- Tramo 10:
t10 t11
- Tramo 11:
t11 t12
- Tramo 12:
t12 t13
(A.26)
que presenta simetra 2/3 y simetra de semionda, con lo que ser suficiente el estudio del
dispositivo en dos tramos, por ejemplo, los tramos 1 y 2, y determinar los instantes de
tiempo t1 y t2 (t3 = t1 + T/6, t4 = t2 + T/6, t5 = t1 + T/3, t6 = t2 + T/3, t7 = t1 + T/2,
t8 = t2 + T/2, t9 = t1 + 2T/3, t10 = t2 + 2T/3, t11 = t1 + 5T/6, t12 = t2 + 5T/6, t13 = t1 + T).
Durante el tramo 1 el rectificador est en estado de conduccin y la tensin de
alimentacin se encuentra en su semionda positiva. En este tramo, el esquema del circuito
pasa a ser el mostrado en la Figura A.9 (se considera un rectificador ideal).
2Z
+
u(t)
i(1) (t)
+
uC(1)(t)
IDC
184
IDC
185
u ( t ) = 2 Ri ( ) ( t ) + 2 L
1
di ( ) ( t ) (1)
uC ( t )
dt
1
du ( ) ( t )
1
i ( ) ( t ) = C C
I DC
dt
1
(A.27)
(1)
duC( ) ( t )
i ( t ) = C
dt
1
(1)
(A.28)
cuya resolucin en este caso conduce a considerar dos casos, segn el signo del
discriminante del polinomio caracterstico, D. Las races del polinomio caracterstico
2 +
R
1
+
=0
L
2 LC
(A.29)
son
=
(A.30)
donde
=
=
R
2L
1
D
2
2
2
R
D=
L LC
Si D > 0 (caso aperidico), se tiene 1 = + , 2 = , y
186
(A.31)
duC( )
(A.32)
( t ) h = C
(1)
dt
(1)
(1)
1t
= K1 C 1e K 2 C 2 e
(1)
2t
(1)
1t
= K1 p K1 e + K 2 p K 2 e
2 t
duC( )
( t ) h = C
(1)
dt
))
2(R+jL)
+
UC(1)
()
U Cpa
1
()
= I (pa) Ipa
1
2 R + j 2 L
1
(1)
(1)
(1)
= I pa
j
= U Cpa UC pa
C
(A.34)
187
(A.35)
()
U Cpb
= I DC 2 R
1
2(R+jL)
+
UC(1)
IDC
Solucin completa
La solucin del sistema (A.27) es la suma de las soluciones homognea y particular
obtenidas. Para el caso aperidico resulta:
()
uC( ) ( t ) = K1( ) e1t + K 2( ) e2t + 2U Cpa
cos t + U( C) pa I DC R
1
()
()
i ( ) ( t ) = K1 p K1( ) e1t + K 2 p K 2( ) e2t + 2 I pa
cos t + Ipa
1
(A.36)
)
sin ( t + K ( ) ) ) +
()
uC( ) ( t ) = K1( ) e t cos t + K 2( ) + 2U Cpa
cos t + U( C) pa I DC R
1
i(
2)
1
2
(A.37)
()
()
I DC
2 I pa
cos t + Ipa
1
A.2.1.1.2 Tramo 2
Para el tramo 2, las expresiones que describen la intensidad de entrada y la tensin del
condensador son:
188
uC( ) ( t ) =
2
( 2)
I DC
t + K1
C
(t ) = 0
(A.38)
(A.39)
(A.40)
(A.41)
(A.42)
Por otro lado, la conduccin termina cuando la onda de intensidad alcanza el valor cero:
f 2 ( t1 , t2 ) = i ( ) ( t2 ) = 0
1
(A.43)
189
190
Bibliografa
1 Libros
[1] J. Arrillaga, D. A. Bradley, P. S. Bodger, Power System Harmonics. John Wiley &
Sons, 1983.
[2] J. Arrillaga, C. P. Arnold, Computer Analysis of Power Systems. John Wiley and Sons,
1990.
[3] A. M. Cohen, Anlisis Numrico, Ed. Revert, S.A, 1977.
[4] H. W. Dommel, Electromagnetic Transients Program Reference Manual (EMTP
Theory Book). Bonneville Power Administration, 1986.
[5] J. J. Grainger, W. D. Stevenson, Jr., Anlisis de Sistemas de Potencia. McGraw-Hill,
1996.
[6] G. T. Heydt, Electric Power Quality, Stars in a Circle Publications, 1991.
[7] C. D. McGuillem, G. R. Cooper, Continuous and Discrete Signal and System Analysis.
Third Edition. Saunders College Publications, 1991.
191
Bibliografa
[8] N. Mohan, T.M. Undeland, W.P. Calm, Power Electronics. John Wiley & Sons, 1995.
[9] A. V. Oppenheim, A. S. Willsky, Seales y Sistemas. Prentice-Hall Hispanoamericana,
1994.
Bibliografa
193
Bibliografa
Bibliografa
195
Bibliografa
196
Bibliografa
197
Bibliografa
3 Otras publicaciones
[66] O. Boix, Estudio y Modelizacin en Rgimen Permanente de Cargas No Lineales
Para el Anlisis Armnico de Redes Elctricas. Tesis Doctoral. Barcelona, 1996.
[67] J. C. Contreras, Rgimen Permanente Rpido con Especificaciones de Potencia para
el Anlisis de la Distorsin Armnica en Redes Elctricas. Tesis Doctoral. Madrid,
1999.
[68] J. Garca Mayordomo, Anlisis de los Armnicos Introducidos en las Redes
Elctricas por los Convertidores Estticos de Potencia. Tesis Doctoral. Madrid,
1986.
[69] J. Pedra, El Problema del Flujo Armnico de Cargas en Redes con Armnicos. Tesis
Doctoral. Barcelona, 1986.
[70] J. A. Rosendo Macas, Procesamiento de Formas de Onda Mediante DFT.
Aplicacin a Protecciones Digitales. Tesis Doctoral. Sevilla, 1997.
[71] L. Sainz, Estudio de la Formulacin y Resolucin del Problema del Flujo Armnico
de Cargas. Tesis Doctoral. Barcelona, 1995.
[72] J. Usaola, Rgimen Permanente de los Sistemas Elctricos de Potencia con
Elementos No Lineales mediante un Procedimiento Hbrido de Anlisis en los
Dominios del Tiempo y de la Frecuencia. Tesis Doctoral, Madrid 1990.
[73] M. Valcrcel, Anlisis del Rgimen Permanente de los Sistemas Elctricos de
198
Bibliografa
199
Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.
Alternative Proxies: