La Maquina Sincronica Muy Buen Apunte Teorico
La Maquina Sincronica Muy Buen Apunte Teorico
La Maquina Sincronica Muy Buen Apunte Teorico
La Máquina Sincrónica
1. Introducción
1
(a) Guri en Venezuela (b) Diablo Cayon en California
2
(a) Estator de la máquina sincrónica (b) Rotor de polos salientes
tensión variable. El principal inconveniente para su uso como motor es que no desarrolla par
de arranque, pero si se incluye en el rotor de la máquina un devanado auxiliar de jaula de
ardilla, es posible obtener par de aceleración como motor de inducción hasta una velocidad
cercana a la de sincronismo, y excitar en el momento apropiado la bobina del campo, con la
nalidad de sincronizar la máquina a la red mediante los pares transitorios adicionales que
se obtienen durante este proceso. Si la fuente de alimentación puede reducir la frecuencia
angular de las tensiones o corrientes de armadura a valores muy bajos, la máquina es capaz
de sincronizarse a esa red y posteriormente ser acelerada al mismo tiempo que se incrementa
paulatinamente la frecuencia de la fuente. Como la construcción de fuentes de gran potencia
controladas en frecuencia es hoy día factible mediante puentes inversores con interruptores
estáticos, es posible que en el futuro esta máquina incremente notablemente su importancia
como accionamiento industrial, e incluso desplace a las máquinas de corriente continua.
3
Durante la operación de la máquina sincrónica en régimen permanente, la velocidad mecá-
nica del rotor es igual a la velocidad angular del campo magnético rotatorio producido por
el estator. En estas condiciones, sobre los conductores o bobinas del campo no se induce
fuerza electromotriz. Para producir fuerza magnetomotriz en el rotor es necesario inyectar
corriente en esta bobina mediante una fuente externa. De esta forma se obtienen dos campo
magnéticos rotatorios que giran a la misma velocidad, uno producido por el estator y otro
por el rotor. Estos campos interactúan produciendo par eléctrico medio y se realiza el pro-
ceso de conversión electromecánica de energía. De acuerdo con la expresión ??, la condición
necesaria, pero no suciente, para que el par medio de la máquina sea diferente de cero es:
ωe = p · ω m (1)
donde:
Te = k · Fr Fe sin δ (2)
donde:
Las fuerzas magnetomotrices del estator Fe , y del rotor Fr tienen una amplitud constante
y para que en la expresión 2 el par medio resulte constante, es necesario que el ángulo δ
entre las dos fuerzas magnetomotrices sea constante en el tiempo durante la operación en
4
(a) Modelo elemental demostrativo
5
régimen permanente. Para lograr esto, las dos fuerzas magnetomotrices deben girar a la
misma velocidad angular.
Cuando la máquina sincrónica se encuentra desequilibrada, el campo magnético rotatorio
producido por las bobinas del estator es elíptico. Este campo se puede descomponer en dos
campos magnéticos rotatorios circulares de sentidos contrarrotativos. Para que sea posible
la producción de par eléctrico medio en estas condiciones, es necesario que la velocidad del
rotor esté sincronizada con uno de los dos campos magnéticos contrarrotativos. El campo
que está fuera de sincronismo y gira en el sentido contrario del rotor, produce par eléctrico
transitorio, pero su valor medio es cero.
Si se cortocircuita la bobina de campo en el rotor de la máquina sincrónica, es posible en
ciertos casos acelerar el rotor como si fuera un motor de inducción con rotor devanado. En el
campo se inducen fuerzas electromotrices con la frecuencia del deslizamiento cuando el campo
magnético rotatorio del estator corta a los conductores del campo. La fuerza electromotriz
inducida en el rotor fuerza la circulación de corrientes por este devanado. Aun cuando el par
eléctrico puede ser muy reducido, en algunas ocasiones este método puede ser utilizado para
arrancar en la máquina sincrónica sin cargas mecánicas acopladas.
d dθ d
[vabc,f ] = [Rabc,f ] [iabc,f ] + [Labc,f ] [iabc,f ] + [Labc,f ] [iabc,f ] =
dt dt dt
= [Rabc,f ] [iabc,f ] + [Labc,f ] p [iabc,f ] + θ̇ · [τabc,f ] [iabc,f ] (5)
El sistema de ecuaciones diferenciales 5 representa el comportamiento dinámico de las bobi-
nas de la máquina sincrónica en coordenadas primitivas. Este sistema se expresa en forma
canónica como:
n h i o
p [iabc,f ] = [Labc,f ]−1 [vabc,f ] − [Rabc,f ] + θ̇ · [τabc,f ] [iabc,f ] (6)
La matriz de inductancia [Labc,f ] depende de la posición relativa θ del rotor con respecto
al estator, por esta razón la matriz de transición de estado también depende de la posición
6
angular del rotor. Si la velocidad de la máquina es constante, la posición angular del rotor
es:
θ = θ 0 + ωm t (7)
La solución del sistema 6 puede obtenerse mediante métodos numéricos de integración, utili-
zando algoritmos tales como Euler, Runge-Kutta o Adams entre muchos otros. El principal
inconveniente que se presenta es la necesidad de evaluar e invertir la matriz de inductancias
de la máquina en cada paso de integración, debido a la dependencia de esta matriz con la
posición angular del rotor. Los computadores personales actuales son capaces de resolver
este problema, aun cuando en el pasado estos cálculos representaba grandes dicultades.
Por este motivo durante varias décadas se desarrollaron transformaciones de coordenadas
que simplican el problema, aceleran notablemente los cálculos y permiten interpretar más
fácilmente el comportamiento dinámico y estático de la máquina sincrónica.
Durante los períodos transitorios, la velocidad angular de la máquina cambia y la posición
angular del rotor es una nueva variable de estado que debe ser evaluada para determinar su
dependencia temporal. En este caso es necesario incorporar una ecuación adicional al sistema
6 para determinar el comportamiento dinámico del eje mecánico de la máquina:
1
[iabc,f ]t [τabc,f ] [iabc,f ] − Tm = J θ̈ + αθ̇ (8)
2
Esta expresión representa el balance del par eléctrico y mecánico en el eje del rotor. El par
acelerante es igual al par eléctrico del convertidor, menos el par resistente opuesto por la
carga y por las pérdidas mecánicas. La ecuación diferencial 8 puede ser expresada mediante
dos ecuaciones diferenciales de primer orden:
( ³ ´
ω̇m = 1 1
[iabc,f ]t [τabc,f ] [iabc,f ] − Tm − αθ̇
J 2 (9)
θ̇ = ωm
donde:
El sistema de seis ecuaciones diferenciales formado por las cuatro ecuaciones del sistema 6,
y las dos ecuaciones mecánicas representadas por la expresión 9, denen el comportamiento
dinámico y transitorio completo de la máquina sincrónica de la gura 3. Este sistema de
ecuaciones diferenciales es no lineal y los coecientes son variables en el tiempo, por este
motivo es necesario recurrir a técnicas numéricas para evaluar el comportamiento de la
máquina o simplicar el problema mediante la técnica de transformación de coordenadas.
En la matriz de inductancia de la máquina sincrónica, se encuentra toda la información
necesaria para determinar su comportamiento. En la matriz de inductancia se resume la
información sobre la disposición geométrica de las bobinas, sus acoplamientos, números de
7
vueltas y reluctancias de los diferentes caminos magnéticos. Una vez conocida la matriz
de inductancias se puede evaluar la matriz de par calculando la derivada parcial de esta
matriz con respecto a la posición angular del rotor. La matriz de inductancias de la máquina
sincrónica esquematizada en la gura 3 posee la siguiente estructura:
· ¸
[Lee (θ)] [Ler (θ)]
[Labc,f (θ)] = (10)
[Lre (θ)] Lf
Laa (θ) Mab (θ) Mac (θ) Maf (θ)
[Lee (θ)] = Mba (θ) Lbb (θ) Mbc (θ) ; [Lef (θ)] = [Lf e (θ)]t = Mbf (θ)
Mca (θ) Mcb (θ) Mcc (θ) Mcf (θ)
donde:
Cada una de las inductancias de la máquina sincrónica se puede representar como una
función del ángulo θ. Esta función es periódica porque se repite nuevamente cada vez que el
rotor realiza un giro completo. Esta propiedad permite representar estas funciones mediante
expansiones en series de Fourier, con el ángulo θ como variable. Si la pieza polar se diseña
convenientemente7 , es posible representar las inductancias de la máquina con un número
reducido de los términos de la serie. La expresión de la matriz de inductancias más simple
consiste en considerar términos dependientes hasta en 2θ, para las inductancias estator-
estator y términos en θ para las inductancias estator-rotor.
La inductancia del rotor Lf , es independiente de la posición θ del rotor debido a que el
estator de la máquina es aproximadamente liso8 . El resto de las inductancias propias y
mutuas depende de la posición angular θ, si el rotor de la máquina es de polos salientes.
Las permeanzas de los caminos magnéticos de las bobinas del estator y de los acoplamientos
estator-rotor son dependientes de la posición angular θ. Cuando la pieza polar del rotor se
encuentra alineada con una de las bobinas del estator, el camino magnético posee la máxima
permeanza. Si la pieza polar se encuentran en cuadratura con la bobina, el entrehierro es
muy grande y disminuye la permeanza. La variación de la permeanza depende del ángulo
2θ porque una bobina alineada con el polo norte del rotor tiene el mismo camino magnético
cuando el alineamiento ocurre con el polo sur. Estas inductancias se pueden representar
aproximadamente mediante las siguientes funciones:
8
4π
Lcc (θ) = L1e + M2e cos 2(θ − ) + ··· (13)
3
π
Mab (θ) = Mba (θ) = −M1e − M2e cos 2(θ + ) + · · · (14)
6
π
Mac (θ) = Mca (θ) = −M1e − M2e cos 2(θ − ) + · · · (15)
6
π
Mbc (θ) = Mcb (θ) = −M1e − M2e cos 2(θ − ) + · · · (16)
2
donde9 :
r
3 3 3
Ld ≡ (L1e + M2e ) ; Lq ≡ (L1e − M2e ) ; Ldf ≡ Mef (17)
2 2 2
Ld + Lq L d − Lq
L1e = ; M2e = (18)
3 3
L1e
M1e ≃ (19)
2
En lo que se reere a los acoplamientos mutuos estator-rotor la funcionalidad de las induc-
tancias es diferente porque al girar el rotor 180o , la bobina del campo invierte su polaridad.
Las inductancias del estator varían entre un valor máximo y un mínimo, siempre positivo
respecto a la posición angular del rotor. Sin embargo, los acoplamientos mutuos estator-rotor
varían entre un valor máximo positivo hasta un valor máximo negativo, que en valor abso-
luto son idénticos, cuando el rotor de la máquina gira 180o . Las inductancias mutuas entre
el estator y el rotor pueden ser aproximadas mediante las siguientes funciones:
Maf (θ) = Mf a (θ) = Mef cos θ + · · · (20)
2π
Mbf (θ) = Mf b (θ) = Mef cos(θ − ) + ··· (21)
3
4π
Mcf (θ) = Mf c (θ) = Mef cos(θ − ) + ··· (22)
3
Si el rotor de la máquina sincrónica es liso, todas las inductancias del estator son indepen-
dientes de la posición del rotor. En esta situación la matriz de inductancias [Labc,f (θ)], se
expresa de la siguiente forma:
L1e M1e M1e Mef cos θ
M1e L1e M1e Mef cos(θ − 2π )
(23)
[Labc,f (θ)] = 3
4π
M1e M1e L1e Mef cos(θ − 3 )
2π 4π
Mef cos θ Mef cos(θ − 3
) M ef cos(θ − 3
) Lf
Aun para el caso de una máquina sincrónica de rotor liso, la solución del sistema de ecuaciones
diferenciales que determina el comportamiento de la máquina sincrónica requiere el uso de
9 En este caso la aproximación se debe a que la dispersión de la bobina no está siendo considerada. La
dispersión puede ser considerada en el modelo como si estuviese conectada externamente a los bornes de la
máquina.
9
métodos numéricos, debido a la dependencia de las inductancias mutuas entre el estator y
el campo, con la posición θ del rotor. El modelo de la máquina sincrónica de rotor liso o de
polos salientes se puede obtener mediante transformaciones del sistema de coordenadas. La
transformación a vectores espaciales y la transformación a coordenadas dq0, introducidas en
el tema II.
donde:
r
2¡
va + αvb + α2 vc
¢
ve = (27)
3
r
2£
1 α α2 [Re ] [ie ] = Re ie
¤
(28)
3
r
2£
1 α α2 [Lee ] [ie ] =
¤
3
r
2£
1 α α2 · · · ·
¤
=
3
− cos 2(θ + π6 ) − cos 2(θ − π6 )
L1e −M1e −M1e cos 2θ
··· · −M1e L1e −M1e + M2e − cos 2(θ + π ) cos 2(θ − 2π ) − cos 2(θ − π ) [ie ] =
6 3 2
−M1e −M1e L1e − cos 2(θ − π6 ) − cos 2(θ − π2 ) cos 2(θ − 4π )
3
3
(L1e + M1e ) ie + M2e ej2θ i∗e =
2
1 1
= (Ld + Lq ) ie + (Ld − Lq )ej2θ i∗e (29)
2 2
r r
2£ 3
1 α α2 [Lef ] if = Mef ejθ if = Ldf ejθ if
¤
(30)
3 2
10
Reemplazando las deniciones de los vectores espaciales en la ecuación 25, se obtiene:
ejθ i∗e + e−jθ ie
½ · ¸ ¾
vf = Rf if + p Ldf + Lf if (31)
2
1 1
Te =[iabc,f ]t [τabc,f ] [iabc,f ] = [ie ]t [τee ] [ie ] + [ie ]t [τef ] if =
2 2
−j π3 j π3 π π
ej 3 e−j 3
−1 −1
e e
M2e
π 4π π 4π
=j [ie ]t ej2θ e−j 3 −e−j 3 e−jπ − e−j2θ ej 3 −ej 3 ejπ [ie ] + · · ·
2 π 8π π 8π
ej 3 e−jπ −e−j 3 e−j 3 ejπ −ej 3
1 1
Mef t
jθ −j 2π −jθ j 2π
··· + j [ie ] e e 3 −e e 3 i =
f
2 −j 4π j 4π
e e
3 3
r
3 Mef 3 © jθ ∗
= M2e (e−jθ ie )2 − (ejθ i∗e )2 + j e ie − e−jθ ie if =
© ª ª
4j 2 2
1
= (Ld − Lq )ℑm (e−jθ ie )2 + Ldf ℑm e−jθ ie if (32)
© ª © ª
2
Las expresiones 26, 31 y 32 modelan la máquina sincrónica utilizando vectores espaciales.
La principal ventaja de esta transformación consiste en la reducción de las tres ecuaciones
del estator a una sola en variable compleja. Por otra parte, aun cuando la dependencia
angular en θ se mantiene en este sistema de coordenadas, las correspondientes expresiones
han sido simplicadas convenientemente al utilizar los términos e±jθ . En la expresión 32
correspondiente al par eléctrico pueden observarse dos componentes: el par de reluctancia y
el par producido entre las fuerzas magnetomotrices del estator y del campo.
idq
e ≡ id + jiq = ie e
−jθ
(34)
11
Derivando la expresión 34 se obtiene la relación siguiente:
1 2
Te = (Ld − Lq )ℑm (idq + Ldf ℑm e−jθ ie if =
© ª © ª
e )
2
= (Ld − Lq ) id iq + Ldf iq if = λd iq − λq id = λdq dq
e × ie (40)
donde:
λd = Ld id + Ldf if ,
λq = Lq iq ,
λf = Lf if + Ldf id ,
λdq
e = λd + jλq .
12
6. Transformación de Park
En la máquina sincrónica, el campo magnético rotatorio producido por las fuerzas mag-
netomotrices de los devanados estatóricos, gira a la velocidad sincrónica ωe . El rotor de la
máquina también gira a la velocidad sincrónica ωr = ωe . Por esta razón es conveniente referir
las ecuaciones diferenciales que denen el comportamiento de la máquina a un sistema de
coordenadas solidario con el rotor. De acuerdo con estos lineamientos se denen los siguientes
ejes magnéticos:
Eje d : Gira con respecto al estator a la velocidad del rotor, y en todo momento se
encuentra colineal con el eje magnético del campo.
Eje q : Rota con respecto al estator a la velocidad del rotor, y en todo momento se
encuentra en cuadratura con el eje magnético del campo.
Eje 0 : Fijo en el estator y se encuentra desacoplado magnéticamente del resto de los
ejes de la máquina.
Eje f : Solidario con el sistema rotórico y colineal con el eje magnético de la bobina de
campo.
Aun cuando los ejes d y q giran a igual velocidad que el rotor, estos ejes representan variables
del estator. El eje 0 es necesario para permitir que la transformación de coordenadas sea bi-
direccional, es decir, se pueda transformar de variables primitivas a variables dq0 y viceversa.
El eje 0 tiene una estrecha relación con las variables de secuencia cero de la transformación
de componentes simétricas. En la práctica el eje 0 permite representar ujos de dispersión
que no están acoplados con otras bobina de la máquina. En la gura 3(b) se ha representado
el sistema de coordenadas dq0 − f .
La matriz de transformación de coordenadas dq0 − f a coordenadas primitivas se dene
mediante la relación:
[iabc,f ] = [A] [idq0,f ] (42)
Si la transformación anterior se escoge de tal forma que la matriz [A] sea hermitiana10 ,
la transformación es conservativa en potencia. Cuando la matriz es hermitiana y real, se
obtiene:
[idq0,f ] = [A]−1 [iabc,f ] = [A]t [iabc,f ] (43)
La matriz de transformación [A] se puede obtener multiplicando la transformación de coor-
denadas primitivas a coordenadas ortogonales αβ011 , por la transformación de coordenadas
αβ0 a coordenadas dq012 :
√1
ia
r 1 0 2 iα
√
2 1
ib = −2
3 √1 iβ (44)
2√ 2
3
ic − 21 − 23 √1 i0
2
10 Inversa de la matriz de transformación [A] igual a su traspuesta conjugada.
11 Transformación de Clark.
12 Rotación en θ introducida en el tema II.
13
iα cos θ − sin θ 0 id
iβ = sin θ cos θ 0 iq (45)
i0 0 0 1 i0
r √1
cos θ − sin θ
ia 2 id
ib = 2 2π 2π
¡ ¢ ¡ ¢ 1
cos ¡θ − 3 ¢ − sin ¡θ − 3 ¢ √2 iq (46)
3
ic cos θ − 4π 3
− sin θ − 4π 3
√1
2
i0
cos ¡θ − 2π cos ¡θ − 4π
¡ ¢ ¡ ¢
cos θ ¢ 0
id 3
2π
¢ 3
4π
ia
− sin θ − sin θ − 3 − sin θ − 3 0
r
= 2
iq ib
√1 √1 1
√ 0 (47)
i0 3
2 2 2 q ic
if 3 if
0 0 0 2
donde:
[Rdq0,f ] = [A]t [Rabc,f ] [A] (50)
d
[Gdq0,f ] = [τdq0,f ] + [Hdq0,f ] = [A]t [τabc,f ] [A] + [A]t [Rabc,f ] [A] (52)
dθ
Por otra parte, la ecuación dinámica del movimiento se puede expresar de la siguiente forma:
1
J θ̈ + ρθ̇ = [idq0,f ]t [τdq0,f ] [idq0,f ] − Tm (53)
2
14
Figura 4: Modelo en coordenadas dq0 − f de la máquina sincrónica
En un sistema trifásico sin neutro no circula corriente de secuencia cero, pero cuando las tres
corrientes de fase encuentran un camino de retorno, es necesario considerar esta componente.
La componente de secuencia cero representa la circulación de corrientes iguales y en fase por
las bobinas de la máquina. Estas corrientes no producen magnetización debido a que la
suma de las fuerzas magnetomotrices de las tres bobinas es cero. Sin embargo, los ujos de
dispersión si poseen componente de secuencia cero. En el modelo de la máquina no existe
acoplamiento magnético de esta secuencia con el resto de las bobinas. Esta componente no
puede producir par eléctrico, pero inuye en las pérdidas de la máquina y en las fuerzas
electromotrices sobre las bobinas. En la expresión 54 no aparecen fuerzas electromotrices de
generación sobre la bobina de campo. Esto se debe a que el sistema de coordenadas dq0 es
solidario al eje f del campo. Los ujos de las bobinas d y q no cruzan tangencialmente a los
conductores del campo. Sin embargo, en este eje pueden aparecer fuerzas electromotrices por
transformación, debido a que el ujo de la bobina del eje directo atraviesa el devanado de
campo. Por el contrario, el eje cuadratura no puede producir ningún efecto sobre el campo
debido a que se encuentra permanentemente en una posición ortogonal.
En la gura se presenta el modelo en coordenadas dq0 − f que satisface las ecuaciones 54. En
la máquina real, las corrientes id e iq no circulan por ningún devanado físico, para determinar
las corrientes reales es necesario aplicar la transformación inversa de coordenadas dq0 − f a
coordenadas primitivas.
15
Cada pareja de escobillas separa las capas de corriente de las bobinas equivalentes. La fuer-
za electromotriz de todos los conductores que forman cada una de las bobina se obtiene
en bornes de las escobillas. Cuando por un par de escobillas se inyecta una corriente, esta
circula entrando a los conductores a la derecha del eje que dene la posición de estas esco-
billas, y saliendo en los conductores a la izquierda. Esta conguración produce una fuerza
magnetomotriz orientada en el eje de las escobillas tal como se muestra en la gura 4.
Las fuerzas electromotrices de generación que aparecen sobre los conductores se recolectan
en los circuitos que se encuentra en cuadratura con el ujo que las producen. El campo y
la bobina del eje d producen generación sobre la bobina del eje q , y la bobina del eje q
produce generación sobre el eje d, pero sobre la bobina de campo no se produce generación
por que este devanado no es cortado por el ujo de los demás ejes. En el sistema de referencia
utilizado, las fuerzas electromotrices de generación aparecen adelantadas 90o con respecto a
los ujos que las producen. Si en las bobinas primitivas se inyecta un sistema balanceado de
corrientes trifásicas, se obtienen las siguientes corrientes en el sistema de coordenadas dq0:
cos ¡θ − 2π 4π
r ¡ ¢ ¡ ¢
id cos θ 3 ¢
cos θ − 3 ¢ √ cos(ωt + α)
iq = 2 − sin θ − sin θ − 2π − sin θ − 4π
¡
3 3
2Ie cos(ωt + α − 2π ) =
3
3 √1 √1 √1 4π
i0 2 2 2
cos(ωt + α − 3
)
id √ cos(θ − ωt − α)
iq = 3Ie − sin(θ − ωt − α) (55)
i0 0
Si la posición angular θ del rotor se sincroniza13 con la variación angular de las corrientes en
la expresión 55, las corrientes en las coordenadas dq0 son independientes del tiempo. En esta
condición, los términos que dependen de las derivadas de las corrientes se anulan. Corrientes
constantes en el tiempo en este sistema de coordenadas, producen fuerzas magnetomotrices
constantes en las bobinas dq0 transformadas. Como la transformación está sincronizada con
la velocidad angular de las corrientes durante el régimen permanente, el campo magnético
producido por las bobinas d y q , gira con la misma velocidad y como resultado se obtiene
el mismo campo magnético rotatorio de la máquina sincrónica en coordenadas primitivas,
excitada mediante un sistema trifásico balanceado de corrientes.
El par electromagnético de la máquina está determinado por la interacción entre fuerzas
magnetomotrices no alineadas. Por una parte la fuerza magnetomotriz del campo produce
par al interactuar con el ujo de la bobina que representa al eje q . La fuerza magnetomotriz
del eje d produce par en su interacción con la fuerza magnetomotriz del enrollado cuadratura.
Exactamente igual pero con sentido contrario, la fuerza magnetomotriz del eje q produce par
con la fuerza magnetomotriz del eje d. Si la reluctancia de los caminos magnéticos d y q son
iguales, estos dos pares se neutralizan. Cuando la reluctancia del eje d es menor que la del
eje q , el par que produce la fuerza magnetomotriz del eje d sobre el eje q es mayor que en la
dirección contraria y se produce un par neto resultante debido a la variación de reluctancia
entre los dos ejes. Desde otro punto de vista se puede interpretar que la pieza polar intenta
13 θ(t) = ωt + θ0 .
16
alinearse con la fuerza electromotriz resultante en la máquina. Si la máquina posee un rotor
cilíndrico, este par es nulo. En la ecuación 40, el par eléctrico se divide en dos componentes,
la primera es proporcional al producto de la corriente de campo if por la corriente de la
bobina cuadratura iq y la segunda componente depende del producto de las corrientes id e
iq . Esta última componente se anula si la inductancia Ld es igual a la inductancia Lq . La
inductancia Ld está denida por la permeanza del eje directo, mientras que la inductancia
Lq está denida por la permeanza del eje cuadratura.
En la gura 4 se han representado las fuerzas magnetomotrices en coordenadas dq0. Se
observa que sobre la pieza polar aparecerá un par eléctrico que intentará alinear el rotor con
la fuerza magnetomotriz total. Cuando se analizan las fuerzas electromotrices de generación
en el sistema de ecuaciones 54 se observan dos términos similares, el primero depende de la
inductancia Lq , que es proporcional a la permeanza del camino cuadratura y determina la
generación sobre el eje directo, el segundo término depende de Ld , y es proporcional a la
permeanza del camino directo y determina parte de la generación sobre el eje cuadratura.
7. Régimen permanente
17
8. Diagrama fasorial
De acuerdo con esta expresión, el fasor que representa el valor efectivo de la tensión en la
fase a del estator de la máquina sincrónica, en régimen permanente es:
vd vq
Ve = Vd + Vq = Vd + jVq = √ + j √ (64)
3 3
Con un razonamiento similar se obtiene la siguiente expresión para las corrientes en régimen
permanente:
id iq
Ie = Id + Iq = Id + jIq = √ + j √ (65)
3 3
Reemplazando las expresiones 64 y 65 en las ecuaciones 58 y 59, se obtienen las siguientes
relaciones fasorales:
Vd = Re Id + jXq Iq (66)
1
Vq = Re Iq + jXd Id + j √ ef = Re Iq + jXd Id + Ef (67)
3
Ve = Vd + Vq = Re (Id + Iq ) + jXd Id + jXq Iq + Ef ⇒
Ve = Re Ie + jXd Id + jXq Iq + Ef (68)
En las expresiones 66 a 68, los fasores con subíndice d están orientados según la dirección
del eje directo, y los fasores con subíndice q , apuntan en la dirección del eje cuadratura. El
fasor Ef se orienta en la dirección del eje q debido a que representa la fuerza electromotriz
producida por la corriente del campo if sobre el eje q . En la ecuación 67 se observa que
el fasor Ef se obtiene multiplicando por j 15 la fuerza electromotriz ef producida por el
campo y dividiendo este resultado por el factor √13 . Todas las magnitudes de los fasores de
las expresiones anteriores se han denido en términos de valores efectivos,√ por esta razón
no aparece en la denición de cada uno de los términos el coeciente 2. En la ecuación
14 En el sistema trifásico balanceado se tiene: v0 = va + vb + vc = 0.
15 Dirección del eje cuadratura.
18
Figura 5: Diagrama fasorial de la máquina sincrónica de polos salientes en convención motor
fasorial 68 aparecen los términos jXd Id y jXq Iq , aun cuando aparentan ser caídas de tensión
reactivas, en realidad representan fuerzas electromotrices de generación. Es necesario recordar
que el operador imaginario j , produce una rotación de 90o . Como el fasor Xd Id está dirigido
según el eje directo, el fasor jXd Id se orienta según la dirección del eje cuadratura. En
otras palabras, el ujo producido por la bobina del eje directo de la máquina, corta a los
conductores jos del estator e induce fuerza electromotriz de generación en el eje cuadratura.
De forma semejante el término Xq Iq representa un fasor con dirección cuadratura, jXq Iq rota
90o y el fasor resultante apunta en la dirección negativa del eje directo. En la gura 5 se
representa el diagrama fasorial de la máquina sincrónica en régimen permanente.
Si el rotor de la máquina sincrónica es liso, las reactancias directa y cuadratura son iguales,
en este caso se dene una sola reactancia denominada reactancia sincrónica Xs . Para la
máquina sincrónica de rotor liso la ecuación fasorial 68 se simplica cuando se agrupan los
términos de generación:
Ve = (Re + jXs ) Ie + Ef (69)
Las relaciones anteriores están escritas en la convención motor. En otras palabras, las corrien-
tes que circulan por las bobinas de la máquina entran por su punto de polaridad relativa. En
la convención motor una potencia positiva indica que la máquina consume potencia eléctrica.
Si la potencia es negativa, la máquina genera potencia eléctrica. Las máquinas sincrónicas
son utilizadas con mucha frecuencia como generadores y es ventajoso en estos casos utilizar
la convención generador en lugar de la convención motor para describir su comportamiento.
En la convención generador las corrientes de armadura salen por el punto de polaridad de
cada bobina. En ambas convenciones, la dirección de referencia de la corriente de campo se
dene entrando por el punto de polaridad relativa por que este eje eléctrico es pasivo y en
general consume potencia eléctrica. El cambio de convención se realiza invirtiendo el sentido
de circulación de las corrientes de los ejes directo y cuadratura, para este n se cambia el
signo de las corrientes Ie , Id e Iq , en las ecuaciones 66, 67 y 68. La fuerza electromotriz que
produce el campo no cambia de signo en la nueva convención, debido a que la corriente de
campo if mantiene la misma referencia en las dos convenciones. De esta forma, la ecuación de
la máquina sincrónica de polos salientes en régimen permanente y en convención generador
19
Figura 6: Diagrama fasorial de la máquina sincrónica de polos salientes en la convención
generador
20
9. Potencia y par eléctrico
Para calcular del par eléctrico se puede utilizar las expresiones 40 o 61. Sin embargo, las
variables independientes de esta ecuación son cticias, por esta razón es conveniente expresar
el par y la potencia eléctrica mediante variables asociadas con el diagrama fasorial. Las
máquinas sincrónicas tienen rendimientos muy altos, particularmente cuando son de gran
potencia. En una máquina sincrónica típica, la potencia mecánica en el eje es prácticamente
igual a la potencia eléctrica en bornes de la máquina. Empleando esta aproximación es
posible desarrollar expresiones del par y de la potencia eléctrica dependientes de variables
mesurables en la práctica. Con estas condiciones se tiene:
Pm = Tm · ωm ≈ Pe = Te · ωe (78)
La potencia eléctrica se determina de la siguiente forma:
Pe (t) = va ia + vb ib + vc ic = vd id + vq iq + v0 i0 (79)
En régimen permanente equilibrado, las corrientes y las tensiones en coordenadas transfor-
madas son independientes del tiempo. La corriente y la tensión de secuencia cero son nulas.
La potencia eléctrica se calcula como:
√ √ √ √
Pe (t) = vd id + vq iq = 3Vd 3Id + 3Vq 3Iq = 3(Vd Id + Vq Iq ) (80)
Despreciando la caída de tensión en la resistencia Re en el diagrama fasorial representado en
la gura 6, se deducen las siguientes relaciones:
Ef − Ve cos δ
Ve cos δ + Xd Id = Ef ⇒ Id = (81)
Xd
Ve sin δ
Ve sin δ = Xq Iq ⇒ Iq = (82)
Xq
Vd = Ve sin δ (83)
Vq = Ve cos δ (84)
Reemplazando las ecuaciones 81 a 84 en la expresión 80 se obtiene el siguiente resultado:
Ef Ve Xd − Xq 2
Pe = 3 sin δ + 3 V sin 2δ (85)
Xd 2Xd Xq e
21
(a) Potencia activa (b) Potencia reactiva
ujo del campo y la tensión de armadura. El ángulo de carga dene el estado o punto de
operación de la máquina, es análogo a la variable deslizamiento en el caso de la máquina de
inducción. En la gura 7(a) se presenta el gráco potencia eléctrica con respecto al ángulo
de carga para una máquina sincrónica típica, indicando las dos componentes de la potencia
eléctrica y la potencia eléctrica total.
La potencia aparente en el estator de la máquina sincrónica se calcula de la siguiente forma:
Se = 3Ve · I∗e = 3(Vd + jVq )(Id − jIq ) =
22
Figura 8: Variación de la potencia eléctrica con el ángulo de carga y punto de máxima
potencia
∂Pe Ef Ve Xd − Xq 2
(δmax ) = cos δmax + V cos 2δmax = 0 (89)
∂δ Xd Xd Xq e
Para las máquinas sincrónicas de rotor liso, las reactancias directa y cuadratura son iguales,
y en este caso se obtiene a partir de la expresión 89:
π Ef Ve
δmax = arc cos(0) = ⇒ Pe max = (93)
2 Xs
23
10. Valores nominales de la máquina sincrónica
Tanto los transformadores, las máquinas de corriente continua, las máquinas de inducción
como las máquinas sincrónicas poseen un conjunto de valores nominales o datos de placa
que determinan sus límites de operación en régimen permanente. Los valores nominales
están denidos por los límites térmicos del convertidor. Si estos límites de temperatura son
excedidos durante períodos prolongados de tiempo, se produce un rápido envejecimiento del
aislamiento de la máquina y un deterioro prematuro de la misma. Los datos de placa de las
máquinas sincrónicas son los siguientes:
Está determinada por las pérdidas Joule de las bobinas del estator, y por las caracte-
rísticas térmicas de la máquina. Las características físicas y químicas de los materiales
utilizados en el aislamiento de los devanados de armadura, denen la máxima tempera-
tura de operación que permiten estos materiales sin que se produzca una degradación
prematura de los mismos. El sistema de intercambio de calor juega un papel impor-
tante porque a medida que aumenta la capacidad de extracción de calor, es posible
incrementar las pérdidas sin aumentar la temperatura máxima de las bobinas. El fa-
bricante dene durante la etapa de diseño, y posteriormente en el banco de pruebas, el
valor máximo de la corriente de armadura que no excede la temperatura máxima del
aislamiento con el sistema de refrigeración utilizado por la máquina. Durante la opera-
ción de la máquina este valor puede ser sobrepasado por un cierto tiempo, aun cuando
exceder la corriente nominal de armadura durante periodos prolongados de tiempo o
en régimen permanente, reduce la vida útil de la máquina.
24
(a) Factor de potencia unitario (b) Condensador sincrónico
Figura 9: Operación de la máquina con potencia aparente nominal con cos φ unitario y
puramente inductivo
encuentran entre 1, 0 y 1, 4 W
m2
b
para los materiales ferromagnéticos de grano no orien-
tado utilizados habitualmente en la construcción de las máquinas eléctricas rotativas.
25
Figura 10: Máquina sincrónica de polos salientes en el punto nominal de operación
26
6. La velocidad nominal:
2πf
ωn = ωsin = (96)
p
A partir del modelo 54 que dene el comportamiento dinámico de las corrientes de la máquina
sincrónica en convención motor, se puede modelar la máquina mediante cuatro circuitos
eléctricos acoplados por términos de generación y transformación, mediante transformadores
y fuentes de tensión dependientes de corriente. En la gura 11(a) se presenta el circuito
equivalente de la máquina sincrónica de polos salientes. Durante la operación equilibrada
en régimen permanente, las corrientes id , iq e if son constantes en el tiempo, y la corriente
i0 es nula. De esta forma, el circuito desacoplado correspondiente a la secuencia cero no
tiene inuencia, y las inductancias del resto de los circuitos no producen caída de tensión.
La corriente del campo se puede calcular evaluando el cociente entre la tensión aplicada al
campo y la resistencia de esta bobina. Los dos circuitos restantes, correspondientes al eje
directo y cuadratura, están congurados tan sólo mediante resistencias y fuentes de tensión
dependientes de corrientes que circulan por otros circuitos. Asociando con el eje real la
polaridad positiva de la corriente y tensión del circuito correspondiente al eje directo, y con el
eje imaginario la del circuito cuadratura, se representa en la gura 11(b) el circuito equivalente
fasorial de la máquina sincrónica de polos salientes en régimen permanente equilibrado.
Si la máquina sincrónica es de rotor liso, las reactancias del eje directo y del eje cuadratura
son iguales y se denomina entonces reactancia sincrónica Xs . Para máquinas sincrónicas de
rotor liso se deduce de la gura 11(b), y de las expresiones 33 y 34, la siguiente relación:
Ve = Vd + jVq = (Re + jXs )(Id + jXq ) + jEf = (Re + jXs )Ie + Ef (97)
27
(a) Modelo transitorio de polos salientes
28
12. Análisis de la máquina sincrónica considerando la sa-
turación
29
Figura 12: Característica de magnetización en vacío de la máquina sincrónica
30
Figura 13: Diagrama fasorial de fuerzas magnetomotrices y enlaces de ujo de la máquina
sincrónica
Las permeanzas de los caminos magnéticos del eje directo y del campo son iguales, por lo
tanto es posible utilizar la característica de magnetización para corregir la reactancia del
eje directo de la máquina sincrónica en un punto de operación determinado. La reactancia
del eje cuadratura no se satura por que su permeanza es muy reducida. Este hecho permite
identicar la posición del eje cuadratura y determinar de esta forma la magnitud del fasor
Eeq . Si esta componente de la fuerza electromotriz fuese generada mediante el ujo producido
por la corriente de campo, se necesitaría un determinado valor de corriente de campo en la
máquina no saturada y una magnitud mayor en la máquina saturada. La proporción entre
estas dos corrientes es dependiente directamente de la variación de permeanza existente entre
el caso lineal y el saturado. La reactancia de magnetización del eje directo se puede expresar
como:
Xmd = ωLd = ωNe2 ℘md = ωNe2 ℘df (100)
En las dos expresiones anteriores se observa que las permeanzas de los caminos magnéticos
son iguales, por esta razón se utiliza la característica de magnetización o curva de vacío para
evaluar la permeanza del camino magnético. En la gura 13 se observa que para producir
una determina fuerza electromotriz Eq es necesaria la corriente if 0 para una máquina lineal
y if sat para la máquina saturada. Mediante la expresión 101 se deduce la siguiente relación
31
entre la permeanza saturada y lineal de la máquina para el camino magnético que une al
campo con la bobina del eje directo:
1 1
Eeq = √ ωNe Nf ℘df lin if 0 = √ ωNe Nf ℘df sat if sat (102)
3 3
De la expresión 102 se puede deducir el siguiente resultado:
℘df lin if
= sat = s (103)
℘df sat if 0
donde:
La expresión anterior indica que la reactancia de magnetización saturada Xmd sat es menor
que la reactancia lineal Xmd lin . debido a que el grado de saturación s es siempre mayor o
igual a la unidad según la denición 103. La reactancia saturada total del eje directo se
calcula añadiendo a la reactancia saturada de magnetización del eje directo, la reactancia de
fuga de la máquina. Esta última reactancia es independiente del grado de saturación de la
máquina:
Xmd lin Xd lin − Xσ 1 s−1
Xd sat = Xmdsat + Xσ = + Xσ = + Xσ = Xd lin + Xσ (105)
s s s s
por lo tanto:
1 if 0 Eeq
Ef sat = √ ωLdf lin if = if (107)
3 if sat if sat
La ecuación 107 indica que para el punto de operación dado, la característica que determina
la fuerza electromotriz producida por el campo en la condición de carga es una recta que pasa
por el origen y por el punto de intersección entre la fuerza electromotriz Eeq y la corriente
if sat , tal como se ilustra en la gura 12.
32
Si se conocen las magnitud de la tensión y corriente de armadura, el ángulo del factor de
potencia, las reactancias lineales de la máquina y la característica de magnetización en vacío,
se puede obtener el diagrama fasorial completo y la corriente de magnetización necesaria
para mantener el punto de operación. En primer lugar se procede a ubicar la posición del eje
cuadratura, multiplicando la corriente de armadura por la reactancia cuadratura y sumando
este término a la tensión de armadura.
El segundo paso consiste en calcular el fasor Ee sumando a la tensión de armadura la caída
en la reactancia de fuga. Posteriormente se proyecta ortogonalmente este fasor con respecto
al eje cuadratura para obtener la magnitud de Eeq . Con esta fuerza electromotriz se entra en
la curva de vacío y se calcula el grado de saturación correspondiente al punto de operación
de la máquina. El grado de saturación s, determinado mediante la característica de vacío
se utiliza para corregir la reactancia del eje directo mediante la expresión 105, y con esta
reactancia se calcula la fuerza electromotriz producida por el campo.
Para determinar la corriente de campo necesaria para mantener el punto de operación es
necesario utilizar la ecuación 107. Si la máquina cambia su punto de operación, es necesario
recalcular el nuevo grado de saturación y evaluar la reactancia del eje directo y la corriente
del campo. Las guras 12 y 13 muestra el procedimiento descrito anteriormente.
El análisis de la saturación de la máquina sincrónica de rotor liso presenta un inconveniente
práctico. En la máquina de polos salientes, la reactancia del eje cuadratura no se satura
debido a que en este eje la reluctancia es muy grande. En los rotores lisos, este hecho es dife-
rente y la reactancia de eje cuadratura también se satura. Si esta reactancia está saturada y
el grado de saturación depende del punto de operación de la máquina, no es posible evaluar
directamente el valor de esta reactancia sin determinar previamente la posición del eje cua-
dratura. Como esta posición está indeterminada, se utiliza un proceso iterativo para localizar
el eje. En principio se puede considerar que la máquina no está saturada en el eje cuadratura
y calcular con esta aproximación la proyección de la fuerza electromotriz Ee sobre el eje di-
recto. Con esta proyección se determina el grado de saturación del eje cuadratura mediante
el procedimiento descrito anteriormente para el eje directo. Esta primera aproximación al
grado de saturación se utiliza para corregir la reactancia del eje cuadratura y recalcular la
posición de este eje. Con la nueva posición se repiten todos los pasos anteriores hasta que
el grado de saturación en la iteracción anterior y la actual converjan a un valor del error
inferior al que se ha especicado previamente. A partir de este punto, el diagrama fasorial se
determina de igual forma que en el caso de la máquina sincrónica de polos salientes, debido
a que la posición de los ejes de la máquina queda completamente denida. En presencia de
saturación las reactancias del eje directo y cuadratura de la máquina sincrónica de rotor liso
son diferentes.
En el pasado se utilizaba con frecuencia la curva de saturación en carga para analizar los
fenómenos de saturación de la máquina sincrónica. Esta característica representa la tensión en
bornes de la máquina en función de la corriente de campo. Se determina para las condiciones
de mayor caída de tensión en las reactancias, es decir con la corriente de armadura nominal
y carga puramente inductiva. En este caso las caídas en las reactancias están en fase con
la tensión de armadura. La característica de saturación en carga se obtiene a partir de
la característica de excitación en vacío restando la caída en la reactancia de dispersión,
33
Figura 14: Característica de saturación en carga - triángulo de Potier
Al igual que en las otras máquinas eléctricas y en el sistema eléctrico de potencia es con-
veniente utilizar el sistema adimensional de unidades, también conocido como sistema en
por unidad, para cuanticar las diferentes variables que denen el comportamiento de la
máquina sincrónica. Como sucede en los transformadores, la máquina sincrónica posee un
rendimiento elevado en el rango de operación industrial y es conveniente denir la potencia
eléctrica aparente de la máquina como potencia base del sistema adimensional de unidades.
En muchas ocasiones es aconsejable utilizar como base de potencia la potencia aparente
monofásica de la armadura de la máquina, especialmente cuando se denen las bases del
rotor, con la nalidad de centrar las magnitudes en por unidad de las variables y parámetros
asociados con el rotor. Estas variables pueden obtener valores extremos, porque el circuito de
campo consume una potencia que es varios órdenes de magnitud inferior a la de armadura.
Además de jar las bases de potencia, tensión y corriente en los convertidores electromecá-
nicos es necesario denir las bases de par, velocidad, tiempo, impedancia e inductancia. El
par base se dene como el cociente entre la potencia base y la velocidad base de la máquina:
SB
TB = (108)
ωB
34
La velocidad base de la máquina en general se escoge igual a la velocidad sincrónica:
ωe f
ωB = ωs = = 2π (109)
p p
donde p es el número de pares de polos del convertidor y f es la frecuencia de las corrientes
de armadura. Al denir la velocidad base de la máquina, el tiempo base y el ángulo base se
encuentran relacionados; al escoger uno, el otro queda automáticamente denido:
αB αB 1 rad
ωB = ⇒ tB = = (110)
tB ωB ωs
Para convertir las inductancias al sistema adimensional de unidades es necesario denir la
impedancia e inductancia base del circuito, de acuerdo con las expresiones 109 y 110 se tiene:
V
ZB B
VB2
LB =
I
= B = (111)
ωB ωB ωB SB
Con las bases denidas anteriormente, la inductancia e impedancia adimensional es:
L ωB L Z
Lpu = = = = Zpu (112)
LB ωB LB ZB
Las inductancias propias de las bobinas de la máquina sincrónica están constituidas por una
componente de magnetización y otra de dispersión:
Ld = Lmd + Lσd = Lmd (1 + σd ) (113)
Lq = Lmq + Lσq = Lmq (1 + σq ) (114)
Lf = Lmf + Lσf = Lmf (1 + σf ) (115)
donde:
Lσd Lσq Lσf
σd = ; σq = ; σf = (116)
Lmd Lmd Lmf
Las inductancias de magnetización del eje directo y del campo, están relacionadas con la
inductancia mutua entre el eje directo y el campo porque poseen el mismo camino magnético.
Los valores físicos de estas inductancias se diferencian debido a que cada una de ellas está
denida por diferentes números de vueltas:
Lmd = Nd2 ℘md ; Lmf = Nf2 ℘mf ; Ldf = Nd Nf ℘df (117)
donde:
℘md = ℘mf = ℘df (118)
Si los números de vueltas de las bobinas del campo y del eje directo son iguales, las tres
inductancias de la expresión 117 tienen el mismo valor. Deniendo las bases del sistema
adimensional de unidades, es posible obtener esta simplicación. Cuando se dene la potencia
base de la máquina para el circuito de armadura y la tensión base de cada bobina como su
tensión nominal, en el estator la corriente bases es igual a la corriente nominal. En el rotor
35
esta situación es diferente, la tensión nominal de la bobina de campo y la potencia aparente
de armadura no producen una corriente base del mismo orden de magnitud que la corriente
nominal del campo. Los valores en por unidad son extremadamente grandes o pequeños, y
se pierden las ventajas del sistema adimensional de unidades. Una solución posible para este
problema consiste en denir los enlaces de ujo base de cada bobina:
λBmd = Lmd IB = Ldf IBF (119)
λBF m = Ldf IB = Lmf IBF (120)
Multiplicando las expresiones anteriores se obtiene la siguiente relación entre las corrientes
base en la armadura y el campo:
s
IBF Lmd Nd
Lmd IB2 = Lmf IBF
2
⇒ = = (121)
IB Lmf Nf
36
Cuadro 1: Rango típico de los valores de las inductancias de la máquina sincrónica de polos
salientes
Inductancia Rango en pu
Ldf = Lmd = Lmf 0, 7 ∼ 1, 1
Lmq 0, 5 ∼ 0, 7
Lσd ≈ Lσq = σd Ldf (0, 1 ∼ 0, 2) Ldf
Lσf = σf Ldf (0, 2 ∼ 0, 3) Ldf
Ld = (1 + σd )Ldf (1, 1 ∼ 1, 2) Ldf
Lf = (1 + σf )Ldf (1, 2 ∼ 1, 3) Ldf
Lq = (1 + σq )Lmq (1, 1 ∼ 1, 2) Lmq
′ L2df
Ld = Ld − Lf
(0, 27 ∼ 0, 43)Ldf
′ L2df
Lf = L f − Ld
(0, 29 ∼ 0, 47)Ldf
En el sistema adimensional que utiliza las bases de enlace de ujo y potencia aparente
monofásica, las tres inductancias calculadas en las expresiones 126, 127 y 128 son iguales. El
valor en por unidad de estas inductancias depende de la permeanza del camino magnético
mutuo estator-rotor ℘df . Valores típicos adimensionales de las inductancias propias, mutuas
y de dispersión en las máquinas sincrónicas convencionales se presentan en la tabla 1.
37
(a) Disposición física de las bobinas (b) Modelo en coordenadas ortogonales
38
Figura 16: Circuito equivalente de la máquina sincrónica con devanado amortiguador
39
Algorithm 1 Análisis transitorio de la máquina sincrónica con devanados amortiguadores
// ANALISIS TRANSITORIO DE MAQUINAS SINCRONICAS CON DEVANADOS AMORTIGUADORES
// p[i]= -[L]^-1*([R]+w*[G]) + [L]^-1*[v]
// [A] = -[L]^-1*([R]+w*[G])
w = 1.0 // velocidad sincrónica w = 1
rd=0.01;rq=0.01;rad=0.03;raq=0.03;rf=0.01;
R = diag([rd rq rad raq rf],0); // matriz de resistencias
Ld =1.0;Lq =0.6;Lad=0.90;Laq=0.50;Lf=1.0;
Ldad=0.85;Lqaq=0.40;Ldaq=0.40;Lqad=0.80;Lfad=0.70;Ldf= 0.8;
// matriz de inductancias
L = [Ld,0,Ldad,0,Ldf;0,Lq,0,Lqaq,0;Ldad,0,Lad,0,Lfad;0,Lqaq,0,Laq,0;Ldf,0,Lfad,0,Lf];
// matriz de generación
G = w*[0,-Lq,0,-Ldaq,0;Ld,0,Lqad,0,Ldf;0,0,0,0,0;0,0,0,0,0;0,0,0,0,0];
Lin = inv(L); A = -Lin*(R+G); Ain = inv(A); // matriz de transición de estado
[V,g] = spec(A); // autovalores y autovectores
io = [0;0;0;0;1]; // condiciones iniciales
t = 0:.25:160; // tiempo (tbase 1/377 s)
e1=exp(g(1,1)*t);e2=exp(g(2,2)*t);e3=exp(g(3,3)*t);e4=exp(g(4,4)*t);e5=exp(g(5,5)*t);
ip = -Ain*Lin*[0;0;0;0;rf*1]; // solución particular p[i]=0
ih = io-ip; // condiciones homogéneas
k= inv(V)*ih; // coeficientes homogéneos
ke = diag(k,0)*[e1;e2;e3;e4;e5]; i = V*ke; // solución homogénea
for m=1:length(t); // solución completa
i(1,m)=i(1,m)+ip(1); // id(t)
i(2,m)=i(2,m)+ip(2); // iq(t)
i(3,m)=i(3,m)+ip(3); // ida(t)
i(4,m)=i(4,m)+ip(4); // iqa(t)
i(5,m)=i(5,m)+ip(5); // if(t)
end;
scf(0)
plot(t,i(1,:),t,i(2,:),t,i(5,:)) // Figura 17(a)
xgrid
scf(1)
plot(t,i(3,:),t,i(4,:)) // Figura 17(b)
xgrid
scf(2)
ia = (i(1,:).*cos(w*t) - i(2,:).*sin(w*t)); // Transformación de Park
plot(t,ia,t,i(5,:)) // Figura 17(c)
xgrid
Referencias
[1] B. Adkins & R. G. Harley, "The General Theory of Alternating Current Machines,"
Chapman and Hall, London 1975.
[2] P. Anderson & A. Fouad, "Porwer System Control and Stability," The Iowa State Uni-
versity Press, US 1977.
40
(a) Corrientes directa, cuadratura y de campo
42