Máquina Síncrona 2018-Apuntes
Máquina Síncrona 2018-Apuntes
Máquina Síncrona 2018-Apuntes
ESTABILIDAD
Parámetros operacionales
Tratamiento de la saturación
Ecuación de "swing"
1
Introducción
2
De esta forma, la metodología adoptada en nuestro caso es ir de lo más general a lo más
particular, mediante el siguiente plan de trabajo:
3
4
Ecuaciones básicas de circuitos magnéticos
Se usará la notación:
: flujo magnético generado por una bobina
: Número de enlaces de flujo de una bobina (flujo magnético enlazado por cada espira,
multiplicado por el número de vueltas)
lk = flujo de dispersión generado por la corriente ik que enlaza sólo el devanado "k"
5
Comentarios:
A los efectos de la aplicación de estos conceptos al cálculo de los coeficientes de
inducción de la máquina síncrona se desea resaltar, en particular, que los coeficientes de
inducción son proporcionales a las permeancias del camino magnético. En particular (y
dado que las reluctancias de caminos magnéticos compuestos por materiales diversos se
agrupan como las resistencias de un circuito eléctrico): si el flujo magnético circula por
un camino ferromagnético con “gaps” de aire en serie, la permeancia del circuito
completo queda definida casi exclusivamente por la geometría de esos “gaps”
6
Esquema básico de la máquina síncrona
7
Modelo abstracto de la máquina síncrona
8
Coeficientes de inducción de la máquina síncrona
a) Estator
aa = Laa0+Laa2.cos2.
bb = Laa0+Laa2.cos2(-120)
cc = Laa0+Laa2.cos2(+120)
ab = -(Lab0+Lab2.cos2(+30))
ac = -(Lab0+Lab2.cos2(+150))
bc = -(Lab0+Lab2.cos2(-90))
Comentarios:
-Es intuitivo que estos coeficientes tengan una variación periódica con el tiempo, y de
período 180 º, dado que la geometría de entrehierros (y, por lo tanto, la permeancia de
los caminos de flujo) se repite con este período. Es razonable también imponer que los
coeficientes nunca se anulen (independientemente de la posición del rotor, siempre hay
flujos enlazados).La expresiones explícitas indicadas pueden verse como los 2 primeros
términos de los desarrollos de Fourier de estas funciones periódicas.
-Si se hace un desarrollo riguroso a partir de la geometría de la máquina, se puede ver
que algunos de los coeficientes “L” están relacionados entre sí [Alvarez].
b) Rotor
fdfd = Lfdfd kdkd = Lkdkd kqkq = Lkqkq fdkq = kdkq = 0 fdkd = Lfdkd
Comentarios: Si consideramos los flujos creados por una de las bobinas del rotor y
enlazados por otra, las permeancias involucradas no cambian con el tiempo (debido a la
geometría cilíndrica del estator),de ahí que los coeficientes de inducción resulten
constantes. Los coeficientes de inducción mutua entre devanados de eje directo y de eje
en cuadratura son obviamente nulos (por ser bobinas a 90 º)
c) Rotor-estator
afd = Lafd.cos
bfd = Lafd.cos(-120)
cfd = Lafd.cos(+120)
akd = Lakd.cos
9
bkd = Lakd.cos(-120)
ckd = Lakd.cos(+120)
akq = -Lakq.sen
bkq = -Lakq.sen(-120)
ckq = -Lakq.sen(+120)
Si la máquina tiene “p” pares de polos, las situaciones descritas se repiten toda vez que
un nuevo par de polos ocupa en el espacio la posición del par de polos anterior (por
ejemplo: si hay 2 pares de polos las geometrías de entrehierro que definen los caminos
magnéticos estator-estator se repiten cada 90 º en vez de cada 180º).Bastará, por lo
tanto, con cambiar por p en las expresiones anteriores.
10
Ecuaciones de estado instantáneas en coordenadas de fase
v = r.i+d/dt
En notación matricial:
va v fd
v
v = vb
s
v = vk d
r
(con vkd = vkq = 0) v= s
v c vk q vr
( vectores columnas de tensiones del estator, del rotor y de toda la máquina)
i a - i fd
is
is = i b ir = - i k d (observar convención de signos!!) i=
i c - i k q ir
a fd
s = b r = k d = s
c k q r
rfd
Rr = rkd
rkq
Rs
R=
Rr
11
Lss Lsr
L=
Lrs Lrr
v = d/dt - R. i
= -L. i
a) 27 de los coeficientes de la matriz L son dependientes del tiempo (no todos distintos,
debido a la simetría de los coeficientes de inducción mutua)
12
Transformada de Park
tp
T=
U3
cos -sen 1
-1
tp = cos(-120) -sen(-120) 1
cos(+120) -sen(+120) 1
Notación:
vd i d d
tp. v = v p = v q
s s
tp. i = i p = i q
s s
tp. =
s
sp = q
v 0 i 0 0
Comentario
La notación utilizada se debe a que si ma mb mc es una terna de vectores colineal con
los ejes de los devanados del estator resulta:
md = 2/3.(ma.cos + mb.cos(-120) + mc.cos(+120)) (suma de proyecciones en eje "d")
mq = -2/3.(ma.sen + mb.sen(-120) + mc.sen(+120)) (suma de proyecciones en eje
"q")
m0 = 1/3(ma + mb + mc) (componente homopolar).
v sp i sp sp
vp = ip = p =
vr ir r
13
vp = T.d/dt (T-1p) - (T.R.T-1) ip
p = -(T.L.T-1) ip
gp 0
G=
0 0
0 -1 0
gp = 1 0 0
0 0 0
p = -Lp . ip
Ld 0 0 Lafd Lakd 0
0 Lq 0 0 0 Lakq
Lp = 0 0 L0 0 0 0
3/2.Lafd 0 0 Lfdfd Lfdkd 0
3/2.Lakd 0 0 Lfdkd Lkdkd 0
0 3/2.Lakq 0 0 0 Lkqkq
Comentarios
-Más adelante se indica cómo quedan las ecuaciones de Park si uno sustituye los flujos
de la segunda ecuación en la primera (se obtienen, por lo tanto, relaciones tensión-
corriente).La formulación a través de los flujos en la forma aquí indicada es más útil
para introducir la saturación en la teoría y para establecer condiciones de continuidad
antes y después de una perturbación.
14
-Observar que si se supone w=constante, el sistema de ecuaciones diferenciales
obtenido es de coeficientes constantes. Se asume esta hipótesis en el desarrollo
posterior del modelo.
a) Ninguno de sus elementos depende del tiempo (por lo tanto: si uno asume la
aproximación de que w es constante, las ecuaciones de Park son de coeficientes
constantes)
b)Es mucho más esparsa que la matriz original
c)A diferencia de la matriz original, no es simétrica (consecuencia de que t p no se eligió
unitaria)
La transformada de Park “sustituye” los devanados del estator por 3 devanados ficticios
(d, q ,0) con las siguientes características:
a)La relación flujo- corriente para el devanado “d” está dada por: d = -Ld.id + Lafd.ifd +
Lakd.ikd ,lo cuál nos dice que el devanado “d” está “montado” sobre el rotor (dado que
los coeficientes de inducción propia y mutua no dependen del tiempo),y con su eje
sobre el eje directo (dado que interactúa con los 2 restantes devanados del eje directo, y
con ninguno del eje en cuadratura)
b)En forma análoga: el devanado “q” está “montado” sobre el eje en cuadratura del
rotor.
c)La relación flujo-corriente para el devanado “0” está dada por 0 = -L0.i0,por lo que el
devanado “0” debe verse como un devanado ficticio ortogonal al plano d-q (no
interactúa con ninguno de los otros devanados).
Si uno realiza el cálculo riguroso de los coeficientes de inducción (cosa que no hemos
hecho),se puede verificar que L0 no es más que la inductancia de dispersión “Ll” de los
devanados del estator, resultado que es intuitivamente compatible con la idea de que el
devanado “0” no interactúa con ninguno de los devanados del rotor.
La idea de “sustituir” los devanados del estator por devanados en el rotor es análoga a la
que se usa en el análisis de transformadores, cuando las reactancias de cortocircuito del
primario se refieren al lado secundario. En nuestro caso hay que tener en cuenta
adicionalmente que estamos sustituyendo devanados fijos por devanados giratorios, y el
“precio” de esta sustitución es introducir en las ecuaciones tensión-flujo el término
w.G p (la llamada “tensión de velocidad”).Este término es el único “discordante” en
las ecuaciones de Park: si no fuera por él las ecuaciones de la máquina en el dominio de
Park serían las asociadas a 6 bobinas estáticas interactuando electromagnéticamente.
15
Comentarios adicionales
Dividiendo por Lb.ib,(valores base de estator) e identificando los valores p.u con un * :
En que (por definición) Lafd*= (Lafd.ifdb)/(Lb.ib),siendo ifdb la corriente base del devanado
de campo, y en forma análoga: Lakd* = (Lakd.ikdb)/(Lb.ib).
Dividiendo por Lfdb.ifdb,(valores base de rotor) e identificando los valores p.u con un *:
Para que Lafd* = Lfda* debe cumplirse: Lfdb.ifdb2 = 3/2.Lb.ib2,y multiplicando a ambos
lados por w0: efdb.ifdb = 3/2.eb.ib .,siendo efdb y eb las correspondientes tensiones base.
Como habitualmente se eligen como corriente/tensión base en el estator los valores
nominales de pico fase-tierra, la relación anterior impone: Pfdb = Sn , siendo Pfdb la
16
potencia base en el devanado de campo, y Sn la potencia aparente nominal trifásica de la
máquina.
En forma análoga: la otras 2 condiciones de simetría (para Lakd y Lakq) imponen que las
potencias base en los 2 devanados amortiguadores sean también iguales a S n.
b) Que todas las inductancias mutuas estator-rotor en el eje directo sean iguales e
iguales a Lad (2 condiciones)
c) Que todas las inductancias mutuas estator-rotor en el eje en cuadratura sean iguales e
iguales a Laq (1 condición)
Ld 0 0 Lad Lad 0
0 Lq 0 0 0 Laq
Lp = 0 0 L0 0 0 0
Lad 0 0 Lfdfd Lad 0
Lad 0 0 Lad Lkdkd 0
0 Laq 0 0 0 Lkqkq
Comentario: La motivación para obtener una matriz de inductancias simétricas pasa por
poder obtener circuitos equivalentes (en los circuitos “físicos” las inductancias mutuas
son siempre simétricas)
17
Ecuaciones en coordenadas de Park desarrolladas
vd = d d/dt-w. q - rs.id
vq = d q /dt + w. d - rs.iq
v0 = d 0 /dt - rs.i0
vfd = d fd /dt + rfd.ifd
0 = d kd /dt + rkd.ikd
0 = d kq /dt + rkq.ikq
18
Régimen estacionario balanceado
A) Máquina en vacío
1) Condiciones de borde:
is = isp = [0], vfd = Vfd = constante
Comentario
Las relaciones tensión-corriente de los circuitos fd y kd constituyen un sistema de
ecuaciones lineal, de coeficientes constantes. Resolviéndolo y haciendo t→∞ se
obtienen las relaciones anteriores. Observar que no es riguroso imponer “de entrada”
dikd/dt = dikq/dt = 0, dado estas corrientes podrían ser, por ejemplo, senoidales en
régimen.
va = -sen.w.Lafd.Vfd/rfd
vb = -sen(-120).w.Lafd.Vfd/rfd
vc = -sen(+120).w.Lafd.Vfd/rfd
en que = w.t +(0) , y se supone (0)=0
Resulta:
va = 2.EI.cos(w.t+)
vb = 2.EI.cos(w.t+-120)
vc = 2.EI.cos(w.t++120)
con EI = 1/2.w.Lafd.ifd (la llamada “tensión proporcional a la corriente de campo” )
19
El resultado final es que la máquina genera en vacío una terna de fasores de tensión de
secuencia positiva, con módulo proporcional a la corriente de campo y fase igual a un
ángulo que también identifica la posición física del rotor.
B) Máquina en carga
1) Condiciones de borde:
ia = 2.I1.cos(wt++Φ)
ib = 2.I1.cos(wt+-120+Φ)
ic = 2.I1.cos(wt++120+Φ)
Se impone que las corrientes inyectadas al sistema formen una terna de fasores de
secuencia positiva, cuya fase es +Φ (Φ es el ángulo del fasor corriente respecto al eje
q, ver figura).
20
En notación fasorial, para la fase a: I=I1. e j (+Φ) (fasor corriente eficaz de la fase “a”)
Comentario:
Observamos que en este caso las variables transformadas son constantes (no dependen
del tiempo)..
Más en general: sea w la frecuencia asociada a la transformada de Park.
Es fácil ver que la transformada de Park transforma ternas de secuencia positiva de
cosenoides de frecuencia w’ en ternas de cosenoides de frecuencia (w-w’).
En particular: si w=w’ se obtienen constantes, y si w’=0 se obtienen cosenoides de
frecuencia w.
Vfd = rfd.ifd
vd = -rs.id+xq.iq
vq = -rs.iq-xd.id+ w.Lafd.ifd = -rs.iq-xd.id+2.EI
v0 = 0
y sustituyendo id,iq con las expresiones vistas, obtendríamos la solución del problema en
coordenadas de Park .
va = vd.cos - vq.sen
21
vb = vd.cos (-120) - vq.sen(-120)
vc = vd.cos (+120) - vq.sen(+120)
= wt+-90
lo cuál muestra (como era de esperar) que también las tensiones forman un sistema
balanceado de secuencia positiva
Métodos fasoriales
Dado que todo es balanceado, pasamos a trabajar sólo con la fase “a”:
En forma fasorial:
siendo
2.V1 = (vd2+vq2)1/2 tg = -vd/vq
2.I1 = (id2+iq2)1/2 tgΦ = -id/iq
vd = -2.V1.sen vq = 2.V1.cos
id = -2.I1.sen Φ iq =2.I1.cos Φ (relaciones ya obtenidas más arriba)
vd = -rs.id + xq.iq
vq = -rs.iq - xd.id + 2.EI
vq-j.vd = -rs.(iq-j.id) - xd.id - j.xq.iq + 2.EI
22
(vq-j.vd) e j = -rs.(iq-j.id) e j - xd.id e j - j.xq.iq e j + 2.EI e j
V = -rs.I - j.xd.Id - j.xq.Iq + EI
En condiciones de régimen ifd es constante, por lo que EI puede verse como la “fem”
interna de la máquina. Observar asimismo que es un fasor sobre el eje “q”.
Se define:
Eq = EI- j.(xd-xq).Id (Tensión atrás de la reactancia sincrónica de eje en cuadratura,
fasor sobre el eje "q")
Eq = V + rs.I + j.xq.I
23
Circuito equivalente para la máquina de polos salientes
24
Régimen transitorio balanceado aproximado
b) Asumir las derivadas de los flujos nulas no equivale a suponer los flujos constantes
en todo el estudio, sino a suponer que, una vez producida la perturbación que motiva el
estudio, las variaciones de los flujos son sumamente rápidas
Esta hipótesis es más bien discutible (la constante de tiempo de los flujos puede llegar
hasta a unos 20 ciclos [Kundur, pág. 161] ,y más bien se suele adoptar esta hipótesis
debido a que si consideráramos estas derivadas no nulas sería necesario (a efectos de
tener un modelo coherente, dado que los flujos terminales de la máquina se
“enganchan” con la red de trasmisión externa) considerar en el modelo también los
transitorios de la red de trasmisión .y no sólo las ecuaciones algebraicas de la red en
régimen.
Se ha verificado empíricamente, asimismo, que levantar esta hipótesis no altera
sustancialmente los resultados de los estudios [Kundur, pág.170].
Otro argumento habitual para mantener esta hipótesis es que, en la práctica, estos
términos son mucho más pequeños que los términos de “tensión de velocidad” (de la
forma w.) con que están sumados en las ecuaciones de Park.
d) Observar que las hipótesis 1) y 3) significan que sólo consideraremos con precisión
los transitorios del devanado de campo.
vd = -w.q - rs.id
vq = w.d - rs.iq
vfd = dfd /dt + rfd.ifd
25
d = -Ld.id + Lafd.ifd
q = -Lq.iq
fd = -3/2.Lafd.id + Lfdfd.ifd
d = -L'd id + Lafd/Lfdfd.fd
con L'd = Ld - 3/2.Lafd2/Lfdfd ( "Inductancia transitoria de eje directo" )
Se definen:
vd = -rs.id + xq.iq
vq = 2.E'q - x'd.id - rs.iq
Son ecuaciones formalmente similares a las de régimen, cambiando xd por x'd y EI por
E'q, y en forma análoga al caso de régimen, se define el fasor E’q = E’q.ejδ
dE'q/dt = (Efd.- EI)/T'd0 (forma habitual en que se encuentra en la literatura esta ecuación
diferencial).
Observar que (dado id) ésta se puede ver también como una ecuación diferencial en E’q,
dado que de las relaciones:
vq = 2.E'q - x'd.id - rs.iq
vq = 2.EI - xd.id - rs.iq
se obtiene la relación: EI = E’q + (xd-x’d).id/√2
Comentario: Se elige E’q como variable de estado para describir la dinámica del
devanado de campo porque (por ser proporcional a un flujo) no puede variar en forma
instantánea al producirse una perturbación. Esto facilita imponer condiciones iniciales
en el momento de la perturbación (del tipo: E’q (0+) = E’q (0-))
Síntesis:
vd = -rs.id + xq.id
26
vq = 2.E'q - x'd.id - rs.iq
EI = E’q + (xd-x’d).id/√2
Métodos fasoriales
Asumiendo que las corrientes son de secuencia positiva, podemos operar igual que en el
desarrollo hecho en régimen permanente:
E’q puede verse como una fem interna de la máquina en el intervalo de tiempo en que
ψfd no cambia. Si no tenemos en cuenta el regulador de tensión (AVR), este intervalo de
tiempo es habitualmente del orden de varios segundos.
Cuando sí se consideran las variaciones de ψfd (ya sea por la presencia del AVR o
porque se quiere hacer un análisis más preciso), las variaciones de E’q se calculan a
través de la ecuación diferencial vista.
Observaciones:
-Al igual que en el caso de régimen, el ángulo de fase de E'q respecto al eje sincrónico
de referencia es también el ángulo "" del eje "q" de la máquina respecto a ese eje.
-El método fasorial se aplica en este caso no en régimen (que es la aplicación habitual),
sino en el contexto de un análisis cuasiestático: los fasores van cambiando a lo largo del
tiempo, y el método fasorial se considera válido en cada intervalo de tiempo elemental
de cálculo
Una forma sencilla de hacerlo es usar el mismo E q que en los estudios de régimen,
haciendo uso de la relación fasorial E q = E’q + j(xq-x’d).Id para calcular (en cada
intervalo de cálculo) el valor de esta “fuente de tensión”
(Recordamos que Eq verifica:Eq = V + rs.I + xq.I)
Comentario
En los métodos de cálculo más antiguos se usaba habitualmente este otro modelo:
Definimos:
E' = E'q - j.(xq-x'd) Iq ("Voltaje atrás de la reactancia transitoria de eje directo")
E' = V + rs.I + j.x'd.I
E’ no es una fuente de tensión ideal (su módulo depende de la corriente),y su ángulo de
fase respecto al eje de referencia es sólo aproximadamente "" (manipulando las
ecuaciones de otra forma: sumando y restando j.xq.Id sí se obtiene, en vez de E’, un
fasor en el eje q ; pero no es habitual presentarlo de esta forma).
27
Diagrama fasorial: régimen transitorio balanceado, con E’q
28
Régimen subtransitorio balanceado aproximado
Despejando ifd , ikd de las 2 últimas y sustituyendo en la primera se obtiene una relación
de la forma:
d = -L’’d.id + K1. fd + K2. kd
Análogamente:
q = -L’’q.iq + K3. kq
2.V = (vq – j.vd ).e j = 2.(-rs .I - j.x’’d .Id – j.x’’q .Iq + E’’q + j.E’’d) , con
E’’q = ejδ f (fd, .kd) y E’’d = -j.ejδ g (kq ),siendo g y f funciones que se pueden
calcular explícitamente.
E’’2 = E’’q + j.E’’d es una “fuente de tensión” constante en el período de tiempo en que
estos flujos no varían. Este período de tiempo no puede ser mayor que las constantes de
tiempo de los devanados amortiguadores, que es de unos pocos ciclos.
Sólo tiene sentido usar este modelo para un estudio de estabilidad, por lo tanto, con un
método de resolución cuasiestático: estas fuentes de tensión van cambiando su valor a
lo largo del tiempo
Para calcular sus valores es necesario, por supuesto, tener en cuenta también las
ecuaciones diferenciales de Park restantes: dψfd/dt, dψkd/dt, dψkq/dt.
Si uno lo hiciera, aparecerían (en forma análoga a como se definió T’ d0) las constantes
de tiempo subtransitorias T’’d0 y T’’q0.
E’’= V + rs .I + j. x’’d .I
29
Parámetros operacionales de la máquina síncrona
Las ecuaciones de eje directo del rotor en el dominio de Laplace (en forma incremental
respecto a las condiciones de régimen en t=0) son:
30
T’d,T’’d : Constantes de tiempo de eje directo en cortocircuito (transitoria y
subtransitoria)
(Los nombres provienen de que son las constantes de tiempo que definen el andamiento
de la corriente de cortocircuito cuando se impone una falta trifásica en bornes de la
máquina)
Es inmediato que:
-si s0, Ld (s) Ld
-si s, Ld (s) L’’d = Ld .T’’d / T’’d0 . T’d / T’d0
-si s y no se consideran los devanados amortiguadores Ld (s) L’d = Ld .T’d / T’d0
Haciendo el desarrollo explícito se puede ver que L’d y L’’d coinciden con las ya
definidas anteriormente al estudiar el régimen transitorio y subtransitorio balanceado
[Kundur,pág. 143].
Más en general: si uno dibujara el diagrama de Bode Ld(jw) = f(w) para una máquina
típica, se encontraría que Ld(jw) Ld para frecuencias de hasta 0,02 Hz, Ld(jw) L’d
para frecuencias entre 0,2 y 2 Hz, y Ld(jw) L’’d para frecuencias >10Hz (Ver figura)
[Kundur,pág. 160].
Este comportamiento de la inductancia operacional sugiere métodos de medida de los
parámetros de la máquina basados en imponer a la máquina corrientes de estator de
frecuencia ajustable.
31
B) Parámetros de eje en cuadratura
Es inmediato que:
Observación
Si la máquina es de rotor liso, es habitual modelarla con 2 devanados amortiguadores
(reales o ficticios) en el eje en cuadratura, con constantes de tiempo bastante distintas
entre sí. En tal caso, la inductancia operacional de eje en cuadratura pasa a ser
formalmente similar a la de eje directo :
-aparece un segundo juego de constantes de tiempo T’q, T’q0
-se puede definir L’q = lim s Lq (s ) cuando no se considera el devanado
amortiguador de eje en cuadratura de constantes de tiempo más pequeñas.
32
Aplicación típica
33
Se observa en esta gráfica la variación de las reactancias de la máquina a medida que
transcurren los diversos períodos: subtransitorio (dinámica muy rápida de los
devanados amortiguadores), transitorio (dinámica del devanado de campo) y régimen
permanente.
Se ha supuesto la máquina inicialmente en vacío, de ahí que el valor inicial de la
corriente (que debe coincidir con el valor pre-perturbación, dado que los flujos en el
estator no pueden cambian instantáneamente) sea cero.
34
-la componente transitoria se amortigua con constante de tiempo 1 a 3 seg
35
Síntesis de los modelos de máquina y su interconexión con la red
Modelo más simple: fuente de tensión atrás de una reactancia, y despreciando los
efectos de saliencia. La red se representa como un conjunto de impedancias serie y
paralelo interconectadas, eventualmente con sus componentes de secuencia
(componentes simétricas).
A diferencia del anterior, ahora se consideran los efectos de saliencia, por lo que
aparece el modelado de la máquina en coordenadas de Park, y la necesidad de usar la
transformada de Park y su inversa para interactuar con la red. Las reactancias en el
modelo de máquina pueden ser sincrónicas, transitorias o subtransitorias, de acuerdo al
estudio que se esté realizando.
36
La máquina se representa ahora con las ecuaciones de Park completas, incluyendo las
que representan la dinámica del rotor .Este es el modelo más sofisticado que hemos
desarrollado, con las simplificaciones ya explicadas: w = constante y despreciando las
derivadas de los flujos del estator.
Se muestra, finalmente, este modelo aun más sofisticado; en que las simplificaciones
anteriores no son consideradas. Dado que las ecuaciones incluyen las derivadas de los
flujos del estator, es necesario también representar la dinámica de la red externa, que
pasa a ser representada por un conjunto de ecuaciones diferenciales en cada una de las
fases.
Este modelo es típico de los estudios de transitorios electromagnéticos, y no se usa
habitualmente en estudios de estabilidad.
37
Valores típicos de los parámetros de la máquina síncrona
38
TRATAMIENTO DE LA SATURACION
En estas condiciones: V√2 = vq = d = Lafd.ifd = Lad.ifd (en p.u, con w=1 en régimen, y
acudiendo a la forma de la ecuaciones p.u ya vista, que impone la simetría de la matriz
de inductancias) por lo que la curva nos permite obtener para cada punto de operación
el valor del Ksd en condiciones de vacío. Como además vd = q = 0, la curva también
puede verse como at = f(ifd) ,siendo at el flujo total a través del entrehierro (flujo total
en el estator excluyendo el flujo de dispersión).
Se asume habitualmente que (a los efectos del cálculo de Ksd) la curva también es
válida en condiciones de carga, interpretada como una curva at-i ,siendo at el flujo
magnético a través del entrehierro e i la correspondiente corriente magnetizante, que es
una combinación de las diversas corrientes id, ifd, etc.
Por lo tanto: si de la OCC deducimos una curva Ksd = f(at), asumimos que la curva
sigue siendo válida en carga.
39
el flujo leído en la OCC la suposición de que la curva OCC sigue siendo válida en
carga se escribe: Ksd = at/at0
a)El flujo de eje directo a través del entrehierro surge de las relaciones básicas flujo-
corriente:
d = -(Ll+Lad) .id + Lad.ifd + Lad.ikd, y (por definición) el flujo a través del entrehierro se
obtiene de descontar (con el signo adecuado a la orientación elegida de corrientes) el
flujo de dispersión:
ad = d + Ll.id = vq + rs.iq + Ll.id (Ecuación en p.u, con w=1 p.u.; válida en régimen
transitorio o subtransitorio bajo la hipótesis habitual d/dt d 0)
Análogamente:
aq = q + Ll.iq = -vd - rs.id + Ll.iq
Acudiendo a las definiciones fasoriales habituales:
(ad + j.aq ).ej/2 = V + (rs + j.Xl )I (siendo V,I la corriente y tensión terminales).
Si llamamos Ea = V + (rs + j.Xl )I (tensión de entrehierro) resulta Ea= at (en p.u,
valor r.m.s).
(Recordar que, por propiedades de la transformada de Park, es:
2. at=(ad2+aq2)1/2 )
Por lo tanto: para cada condición terminal V, I , calculando V + (rs + j.Xl )I se obtiene
at.
40
Estas aproximaciones analíticas permiten calcular at0 como función de at en cada
tramo, de acuerdo con lo requerido.
Comentarios
-Algunos programas asumen la curva OCC como función cuadrática del flujo de
entrehierro.
-Por lo general los programas de cálculo solicitan que se ingresen 2 puntos de la curva
OCC a efectos de poder aproximarla con la función definida.
41
ECUACION DE “SWING”
J.dwm/dt = Ta
Interesa expresar las ecuaciones en función del ángulo que, como ya vimos, tiene
también un significado en algunos modelos como ángulo de fase de la “fem” de la
máquina :
42
2.H/w0.(d2/dt2) = Pa(pu) - KD.w
Comentarios:
Sea TM el tiempo requerido para acelerar la máquina en vacío desde el reposo hasta su
velocidad nominal, aplicando el torque mecánico nominal y asumiendo KD=0 :
dw/dt = Ta.w0/2H con Ta = 1 = Sb/w0 w0 = w0/(2H) (0,TM) Ta dt = w0/(2H).TM
TM = 2H (TM se denomina “tiempo de arranque mecánico”)
Valores típicos de H:
43
LIMITES OPERATIVOS DE LOS GENERADORES SINCRONOS
Suponemos una máquina conectada a una barra de carga infinita (tensión terminal V
constante, cualquiera sea su inyección de potencia).Consideramos condiciones de
régimen, y despreciando la resistencia de estator.
Las expresiones indicadas son para las potencias por fase .Es fácil ver que también son
válidas para potencias trifásicas p.u (considerando como potencia base la nominal
trifásica, y como tensión base la nominal fase-fase).
44
Límites operativos
Se intenta definir una región en el plano P-Q en la cuál deben encontrarse los puntos de
operación de la máquina en régimen.
Se supondrá en todos los casos tensión terminal V=constante.
C)Límite de estabilidad
Se observa que P como función de es acotada, por lo que para cada EI se puede
obtener el máximo asociado a una potencia P máxima, y a su correspondiente Q. En
particular: para rotor liso máximo=90° (de lo contrario es < 90°)
Haciendo variar EI obtenemos la curva límite P máxima = f(Q),que se denomina “límite
de estabilidad teórico”.
En el caso particular de las máquinas de rotor cilíndrico este límite es la recta vertical
Q = -V2 /Xd
45
a la pequeña señal, y no debería pensarse en una curva límite individual de cada
máquina.
El límite de estabilidad aquí descrito (visto como una cota que cada máquina individual
no debería rebasar) se sigue usando habitualmente, no obstante, como complemento a
los límites que se definan en base a los estudios de sistema, (observar que es un límite
“optimista”, dado que mantener la tensión terminal constante es beneficioso para la
estabilidad de pequeña señal).
Observaciones:
-Para máquinas de rotor liso existen también limitaciones en la zona de subexcitación
vinculadas a calentamientos localizados si se trabaja con corrientes de campo muy
bajas.
-La zona de funcionamiento para corrientes de campo bajas suele limitarse también a
causa del ajuste del relé de subexcitación.
Si una máquina pierde su excitación, pasa a funcionar como un generador de inducción,
consumiendo cantidades importantes de reactiva de la red. La consecuencia práctica es
un descenso de la tensión terminal y un aumento importante de las corrientes de estator
y rotor, de ahí que suelen implementarse protecciones que evitan subexcitaciones
excesivas.
46
Límites operativos de los generadores de la Central Térmica UPM.
47
Referencias bibliográficas básicas
P.Anderson, A.Fouad : Power System Control and Stability, IEEE Press, 1994
48