SIMULACION DE TRANSITORIOS DE LA MAQUINASINCRONICA TRIFASICA CON EL EMTP
HECTOR FABIO MARTINEZ GARAVITOJAIR AGUADO QUINTERO
i.nrvcrs¡0ü0 urr,rr0m0 de 0ccidonfo
1 5514!
CALICORPOHACION UNIVERSITARIA AUTONOMA DE OCCIDENTE
DIVISION DE INGEN IERIASPROGRAMA DE INGENIERIA ELECTRICA
l8f
SIMULACION DE TRANSITORIOS DE LA MAQUINASINCRONICA TRIFASICA CON EL EMTP
HECTOR FABIO MAHTINEZ GARAVITOJAIR AGUADO QUINTERO
Trabajo de grado presentado como requisito parcial paraoptar al título de INGENIERO ELECTRICISTA
Director:LUIS EDUARDO ARAGON
I n ge n i e ro E I éctrici sta
CALICORPORACION UN IVERSITARIA AUTONOMA DE OCCIDENTE
DIVISION DE INGENIERIASPROGRAMA DE INGENIERIA ELECTRICA
1 993
-/
/:,
Nota de AceptaciónAprobado por el comité de trabajo degrado, en cumplimiento de losrequisitos exigidos por la CorporaciónUniversitaria Autónoma de Occidente,para optar al título de lngenieroElectricista.
Presidente de jurado
)t
Cali, Mayo de 1993
AGRADECIMIENTOS
Los autores expresan sus agradecimientos:
A LUIS EDUARDO ARAGON R, 1.E., M.Sc. Profesor de la CorporaciónUniversitaria Autónoma de Occidente, Jefe de Líneas y Subestacionesde la firma consultora GERS Ltda, y Director del Trabajo por susinvaluables aportes y constante dedicacíón en todas las etapasdesarrolladas en esta investigación.
A MARTHA cEClLlA AMAYA E, LE., M.sc. Profesor de la corporaciónUniversitaria Autónoma de Occidente, Directora del programa de tngenieríaEléctrica de la Universidad del Valle y Asesor delTrabajo por losaportesy apoyo ofrecido sin los cuales no hubiese sido posible la realización deeste trabajo.
A ENRIQUE CIRO QUISPE, 1.E., Profesor de la Corporación lJniversitariaAutónoma de occidente por los aportes ofrecidos en éste proyecto.
A La Corporación Universitaria Autónoma de Occidente.
A Los profesores que en eltranscurso de nuestra carrera nos brindaronsus conocimientos.
A todas aquellas personas que en una forma u otra colaboraron en larealización del presente trabajo.
II
DEDICATORIA
Todo elesfuerzo conjugado en este trabajo se lo dedíco a mis padres, que me
brindaron todo su apoyo y que con su esfuerzo y dedicación ayudaron a
ampliar mis conocimientos:
Julio César Martínez A.
Graciela Garavito de M.
A mis hermanos fieles testigos de mis propósitos, y triunfos.
Julio César Martínez G.
María Patricia Martínez G.
Alfredo Adolfo Martínez G.
III
Hector Fabio M.
TABLA DE CONTENIDO
INTRODUCCION
,: , 1. DESCRIPICION DEL FENOMENO TRANSITOHIO
2. PARAMETROS DE LA MAQUINA SINCRONICA TRIFASICA
2.1 INDUCTANCIAS DE LA MAQUINA S/NCRONICA TRIFASICA
2.1.1 Definición de los devanados de un generador sincrónico
2.1.1.1 Discontinuidad eléctrica
2.1.1.2 Discontinuidadmagnética
2.1.2Inductancias propias del rotor
2.l.S lnductancias mutuas relativas a los devanados rotóricos
2. 1 .4 lnductancias estatóricas propias
2.1.5 lnductancias mutuas entre fases estatóricas
2.2 TRANSFORMACION DE PARK
2.2.1 Transformación de los ejes móviles a /os ejes fijos
2.2.2 Ecuaciones de enlaces de flujo
2.2.3 Ecuacíones de voltaje
Pág
1
3
12
13
13
15
15
15
16
19
23
28
29
34
38
IV
2.2.4Normalización de las ecuaciones
2.2.4.1 Cantidades base para el estator
2.2.4.2 Cantidades base para el rotor
2.2.5 Ecuaciones de par y potencia
2.2.6 Representación esquemática de los flujos
2.3 PARAMETROS TRANSITORIOS DE LA MAQUINA SINCRONICATRIFASICA
2.3.1 Reactancias de la máquina sincrónica
2.3.1.1 Reactancia longitudinal transitoria
2.3.1.2 Reactancia longitudinal subtransitoria
2.3.1.3 Reactancia transversal subtransitoria
2.3.2 Constantes de tiempo de la máquina sincrónica
2.3.2.1 Constante de tiempo subtransitoria longitudinal encortocircuito
2.3.2.2 Constante de tiempo subtransítoria transversal encortocircuito
2.3.2.3 Constante de tiempo transitoria longitudinal en coftocircuito
2.3.2.4 Constante de tiempo subtransitoria transversal encircuito abierto
2.3.2.5 Constante de tiempo transitoria longitudinal a circuitoabierto
2.3.2.6 Constante de tiempo subtransitoria longitudinal a circuitoab ie rto
V
pag
45
46
46
49
52
55
55
58
58
58
60
61
60
60
61
61
pág
3. DIAGRAMAS DE BLOQUE DE LA MAQUINA SINCRONICATRIFASICA
3.1 CORTOCIRCUITO BRUSCO DE UNA MAQUINA SINCRONICACONECTADA A UN BARRAJE INFINITO
3.1.1 Expresión normalizada para la corriente de armadura en ejesdirecto y en cuadratura
3.1.1.2 Expresión normalizada para las tensíones y flujos
3.1.1.3 Expresión normalizada para la corriente de excitación
3.1.1.4 Expresión normalizada para la ecuación de par
3.2 DIAGRAMA DE BLOQUE PARA LA MAQUINA SINCRONICAS/N SATURACION
3.2.1 Variables de entrada
3.2.1.1 Tensión de barraje infinito
3.2.1.2 Tensión de excitación
3.2.1.3 Torquemecánico
3.2.2 Variables de salida
3.2.2.1 Tensión inducida en el generador
3.2.2.2 Corriente de salida del generador
3.2.2.3 Corriente de excitación
4. APLICACION DEL EMTP PARA LA SIMULACION DETRANSITORIOS DE LA MAQUINA SINCRONICA
4.1 ESTRUCTURA DE LOS DATOS DE ENTHADA
64
65
65
69
74
76
77
79
79
79
79
79
79
80
80
81
82
VI
LISTA DE FIGURAS
pag.
FIGURA 1. Circuito que representa la característica de cortocircuito. 10
FIGURA 2. Modelo de la máquina bipolar ideal(máquina generalizada). 14
FIGURA 3. Modelización de los diferentes devanados de una máquinasincrónica. 17
FIGURA 4. Circulación del campo estatórico longitudinal,Circulación del campo estatórico transversal. 22
FIGURA 5. Descomposición de la fuerza magnetomotríz de la fase a
en sus componentes longitudinal y transversal. 22
FIGURA 6. Variación de la inductancia propia de una Íase estatóricade una máquina sincrónica. 24
FIGURA 7. Descomposición del flujo de la fase en sus componentesIongitudinal y transversal. 27
FIGURA 8. Variación de la inductancía mutua entre fases estatóricasde una máquina sincrónica. 27
FIGURA 9. Representación de la distribución espacial de la fuerzamagnetomotrí2. 31
FIGURA 10. Descomposición de una onda pulsante en fasoresespaciales 31
FIGURA 11. Proyección de la fase a,b,c en los ejes d,q,O. 33VIII
pág
FIGURA 12,
FIGUBA 13.
FIGURA 14.
FIGURA 15.
FIGURA 16.
FIGURA 17.
FIGURA 18,
FIGURA 19.
FIGURA 20.
FIGURA 21.
FIGURA 22,
FIGURA 23,
FIGURA 24.
FIGURA 25.
Diagrama esquemático de voltajes de la máquinasincrónica.
Flujos en el eje directo.
Flujos en elje transversal.
Esquema equivalente para los flujos en el eje directocaso general.
Esquema equivalente para los flujos en el eje transversal.
Esquema equivalente usual (LDd 0) para los flujosen eje directo.
Definición de Xd'.
Definición de Xd".
Definición de Xq".
Definición de Td".
Definición de Tq".
Definición de Td'.
Definición de T'd1.
Generador sincrónico conectado a un barraje infinito.
44
56
56
57
57
57
59
59
59
62
62
63
63
65
IX
RESUMEN
Este proyecto está basado en la simulación de transitorios de ta máquina
sincrónicatrifásíca, utilizando como herramienta de análisis el programa EMTP.
La idea primordiales la de mostrar el principio de funcionamiento transitorio de
la máquina sincrónica trifásica.
Este proyecto beneficia al estudiante, porque Ie permite visualizar en forma más
clara los conceptos teóricos aprendidos en sus cursos de máquinas eléctricas,
especialmente en las máquinas de corriente alterna (conversión de energía ltl).
y además, el profesor dispondrá de una valiosa herramienta para la enseñanza
del curso.
Este proyecto contará con un material bibliográfico referente a la ecuaciones
físicas, eléctricas y mecánicas que representan la máquina sincrónica trifásica.
X
INTRODUCCION
El presente proyecto contempla un análisis de los transitorios y sus efectos en
la Máquina Sincrónica Trifásica.
Para poder desarrollar la simulación de transitorios de la Máquina Sincrónica
Trifásica se tendrán en cuenta los parámetros referidos en los ejes directo,
transversal y homopolar respectivamente, en asocio con la ecuación de par y
utilizando como herramienta matemática las transformaciones de Park, para
poder obtener el modelo matemático y poder cumplir con nuestro objetivo.
Esta simulación se desarrollará utilizando como herramienta para el análisis el
programa EMTP. Se introducirán los parámetros que corresponden a uno de
los generadores trifásicos de la Central Hidroeléctrica delALTO ANCHICAYA,
diseñados y construidos por la firma CANADIAN GENERAL ELECTRIC - CGE.
2
Es importante este estudio porque elfuncionamiento de la Máquina Sincrónica
Trifásica en régimen transitorio facilita eldiseño de su sistema de protecciones,
y la selección de los diferentes controladores, reguladores, etc. Y la determinación
de su nivel de aislamiento.
Además sirve como referencia en cursos posteriores de Máquinas Eléctricas,
Sistemas Eléctricos de Potencia y Protecciones Eléctricas.
1. DESCRIPCION DEL FENOMENO TRANSITORIO
Las máquinas eléctricas en funcionamiento continuo se ven sometidas a
pertubaciones, durante pequeños intervalos de tiempo y SUS componentes
deben ser capaces de soportar los valores máximos de las tensiones y
corrientes que se originan durante estos períodos transitorios.
Las máximas tensiones y corrientes a que se ven sometidos los distintos
elementosde la máquinasincrónicase producen durante los períodostransitorios
que tienen lugar maniobras, coftocircuitos y pérdida súbita de carga. Como
consecuencia el diseño de todos los elementos de la máquina sincrónica viene
fuertemente condicionado por la consideración de los fenómenos transitorios a
que puede verse sometída.
Los fenómenos transitorios que aparecen en la máquina sincrónica son de
muchas clases y cubren un amplio margen de frecuencias: desde las oscilaciones
electromecánicas asociadas a la resonancia subsincrona, de frecuencia de
algunos Hertzios, hasta los rápidos cambios en las tensiones y corrientes,
4
produc¡dos por la aparición de una falla o simplemente por una maniobra,
capaces de crear componentes de frecuencias del orden de decenas de
Kilohertzios.
En el evento de un cortocircuito los fenómenos transitorios están influidos por
el tipo, el instante en que se produce, las fuentes de corriente implicadas, el
estado previo de carga, la ubicación y la duración delcortocircuito y la forma de
la red, así como por los parámetros característicos y el comportamiento de los
aparatos y componentes implicados.
Como la impedanciadecortocircuitode ta máquinasincrónica puedeconsiderarse
prácticamente como una reactancia inductiva, elvalor inicial de la componente
aperiódica en las condiciones de cortocircuito es casitan grande como el valor
de cresta de la corriente inicial simétrica de cortocircuito. Sin embargo la
totalidad delfenómeno de amortiguamiento no depende sólo de la componente
aperiódica, sino también, en alto grado, de la componente de la amortiguación
de la corriente simétrica de coftocircuito. La corriente simétrica de cortocircuito
se divide en las siguientes componentes:
- Componente subtransitoria
- Componente transitoria
- Componente permanente.
5
Sin embargo, para representar y calcular las relaciones de cortocircuito se
considera en la práctica una tensión constante y se supone que el fenómeno de
amortiguamiento de la corriente simétrica de cortocircuito viene provocado por
un crecimiento de las reactancias de la máquina sincrónica. De acuerdo con las
componentes citadas, las reactancias correspondientes se denominan:
- Reactancia subtransitoria X"d
- Reactancia transitoria X'd
- Reactancia síncrona Xd.
El subíndice d significa que las reactancias están referidas a una posición del
rotor para la cual coinciden los ejes de los devanados rotóricos y estatóricos,
esta posición se denomina eje longitudinal. Además de estas reactancias, la
máquina sincrónica presentaotras según elejetransversalque deben conocerse
también para poder simular con gran exactitud el comportamiento transitorio de
la máquina sincrónica.
Para determinar la variación temporal de la corriente de cortocircuito deben
conocerse también las constantes de tiempo correspondientes a las reactancias.
Estas constantes son las siguientes:
6
- Constante de tiempo subtransitoria Td"
- Constante de tiempo transitoria Td'
- Constante de tiempo aperiódica Ta.
La componente permanente (corriente estacionaria simétrica de cofto circuito
o corriente permanente de cortocircuito) tiene un valor de cresta constante.
Su valor eficaz es el cociente entre la f uerza electromotriz (nominal) fase - neutro
y la reactancia síncrona Xd (debida a la reacción de inducido) que corresponde
al régimen permanente de cortocircuito.
A la corriente estacionaria simétrica de cortocircuito se le superpone una
componente transitoria que se amortigua de forma relativamente lenta. Esta
componente es debida principalmente a las reactancias de dispersión de los
devanados estatóricos y de excitación, pues la máquina sincrónica en coftocircuito
se comporta como un transformador cortocircuitado.
La reactancia para la componente transitoria es considerablemente menor que
la correspondiente a la corriente permanente de cortocircuito. El valor eficaz de
fa corriente transitoria simétrica de cortocircuito es el cociente entre la tueza
electromotriz (nominal) fase - neutro y de la reactancia transitoria Xd'.
7
A la componente transitoria se le superpone durante algunos períodos una
componente subtransitoria. En la máquina sincrónica con devanado amortiguador
aparecen en los circuitos rotóricos unos fenómenos de compensación entre los
circuitos de amortiguación y los de excitación. Estos fenómenos provocan en los
circuitos estatóricos la componente subtransitoria de la corriente simétrica de
cortocircuito, cuyo valor inicial es igual a la amplitud máxima de la corriente
inicial simétrica de cortocircuito. Su valor eticaz es el cociente entre la fuerza
electromotriz (nominal) fase - neutro y la reactancia subtransitoria X"d.
La componente de régimen transitorio de la corriente es diferente en cada fase
porque las tensiones generadas en las fases de una máquina sincrónica
trifásica están desfasadas Fnas de otras en 120 grados eléctricos, es decir el
coftocircuito se aplica en puntos diferentes de la onda de tensión de cada fase.
En una máquina sincrónica, elflujo a través del entrehierro es mucho mayor en
el instante de ocurrir elcortocircuito. A medida que transcurre eltiempo elflujo
disminuye y la corriente del inducido decrece porque la tensión generada por el
flujo en el entrehierro determina la corriente que fluirá a través de la resistencia
y la reactancia de pérdidas del devanado del inducido.
8
En condiciones transitorias en todos los circuitos del rotor pueden existir
corrientes que varíen con eltiempo. Así mismo, las corrientes en elestator no
serán ya senoidales de amplitud constante y en consecuencia tampoco lo serán
las ondas de la componente de luerza magnetomotriz y de flujo en el espacio.
Debe tenerse en cuenta la posibilidad de que la velocidad mecánica del rotor
varíe con eltiempo, dando como consecuencia el que los flujos concatenados
con sus devanados no se mantengan constantes originándose corrientes
inducidas en todos ellos: la corriente en el devanado principal, por ejemplo, no
obedece simplemente a la ley de OHM. Todas estas corrientes, como eS
natural, repercuten sobre la excitación de la máquina sincrónica, ya que
intervienen en et flujo en el entrehierro y por consiguiente influyen sobre los
valores instantáneos de la corriente en el estator.
Los flujos concatenados con eldevanado de excitación principalse consideran
constantes. Su valor iniciales determinado por la inductancia y por la corriente
de excitación antes de producirse la falla y esta constancia deberá mantenerse
frente a la fuerza magnetomotriz desmagnetizante del estator que acompaña a
ta corriente de coñocircuito. Para contrarrestar esta fuerza magnetomotriz
desmagnetizante deberá aparecer una componente inducida de la corriente de
excitación y esta es la que determina es estado de la corriente del estator
durante el período transitorio: esto representa simplemente una excitación en
la máquina sincrónica mayor que la existente en régimen permanente con la
I
consecuencia de que la corriente en el estator durante el período transitorio es
también mayor que en régimen permanente.
La componente continua de la corriente del estator crea una componente de
campo magnético en el entrehierro fija en el espacio y por consiguiente induce
una tensión y una corriente a la frecuencia fundamental en los devanados del
rotor.
Esta componente continua no estará presente siel coftocircuito se produce en
elmomento en que elflujo concatenado con una de las fases delestatores nulo
y la onda de corriente de cortocircuito de esa fase será simétrica. Si por el
contrario, en el momento de producirse elcortocircuito los flujos concatenados
por las fases no son nulos deberán aparecer en ellas las componentes continuas
necesarias para mantenerlos constantes.
La componente continua desempeña efectivamente el mismo papel que la
componente continua transitoria en un circuito RL sencillo cuando se le aplica
súbitamente una tensión alterna. Al igual que en el circuito RL, el mayor valor
posible de la componente continua es igual al mayor valor instantáneo de la
corriente de cortocircuito durante el período subtransitorio y este máximo valor
se alcanza cuando el cortocircuito se produce en el instante en que el flujo
concatenado por una de las fases pasa por un máximo; en este caso la onda de
rrnlvcfst0J0 - ' r,0fll0 de ()ccidcnfc
Secr:ión Biblioteco
10
corr¡ente de cortocircuito queda descentrada y a un lado del eje cero. Esta
componente continua al no estar sostenida por alguna tensión en el estator va
perdiéndose a un ritmo que depende de la resistencia y de la inductancia
equivalente del circuito del estator.
La componente continua se extingue rápidamente, pero inicialmente puede
llegar a atcanzar el 50 ó 60 por ciento de la magnitud de la componente de
corriente alterna un instante después de la falla. De está forma, la corriente
inicialtotal llega a sertípicamente 1 .5 ó 1 .6 veces la magnitud de la componente
alterna sola.
La característica de cortocircuito se obtiene teniendo en cuenta el siguiente
circuito:
a
I
i iz,l
IIII
Eq
FIGURA 1. Circuito que representa la caracterÍstica de cortocircuito.
V=iR+Ldi/dt
i(t) = Vm/Z [sen(wt +cr - 0) - sen(cr - $)e-RVL]
Donde Z=1ffiml ; =Tg-l(wUR),f
11
Después del período transitorio, la corriente de falla alcanza su condición de
estado estacionario llamada comúnmente l. Su magnitud se calcula
aproximadamente por la ecuación siguiente:
| = EalX
La corriente subtransitoria está definida por la siguiente ecuación:
l" = Ea./X"
De la misma manera la corriente transitoria es:
l'= EalX'
Debe observarse que estas ecuaciones son válidas únicamente para falla
trifásica simétrica. Para fallas asimétricas el análisis se hace empleando el
método de las componentes simétricas.
2. PARAMETROS DE LA MAQUINA SINCRONICA TRIFASICA
Las máquinas rotativas de corriente alterna se diferencian entre sí por la
combinación, entre los arrollamientos del estator y el rotor.
Sin embargo el estudio de cualquier tipo de máquina se desarrolla a partir de lo
que se denomina máquina bipolar ideal, que matemáticamente será equivalente
a la real. Para desarrollar el modelo matemático de la máquina sincrónica se
parte de dos suposiciones a saber:
- Permeabilidad infinita (para despreciar los efectos de la saturación)
- Sólo se tomarán en cuenta los armónicos fundamentales de las ondas de
tensión, flujos y fuerza magnetomotriz ( no se considerarán los efectos de los
armónicos superiores).
Con base en el modelo de la máquina bipolar ideal (máquina generalizada) se
obtiene un modelo semejante al de la figura 2. Para la máquina trifásica. Para
esta transformación en la máquina sincrónica se definen como los ejes d y q,
13
el longitudinaly eltransversal respectivamente a la rueda polar, lo cual equivale
a considerar fija la rueda polar y móvil el estator. De esta manera las
inductancias mutuas (Lij) se pueden considerar constantes e independientes
del ángulo formado por los ejes elegidos y el eje definido de la rueda polar.
Una vez elegidos los ejes vemos que el arrollamiento de campo coincide con el
eje d y que hay que sustituir el arrollamiento trifásico por dos bobinas fijas a la
rueda polar.
2.1 INDUCTANCIAS DE LA MAQUINA SINCRONICA TRIFASICA
2.1.1 Definición de los Devanados de un Generador Sincrónico. El devanado
estatórico contiene tres fases idénticas desfasadas entre ellas 120 grados
eléctricos, estas tres fases son indicadas respectivamente por a, b y c. El rotor
contiene eldevanado de excitación (f) de corriente directa, distribuido siguiendo
los ejes polares. Sobre la cara de los polos se construye un devanado
amortiguador constituido por barras de cobre cortocircuitadas por dos anillos
como la jaula de un motor de inducción.
Para el análisis del compodamiento del devanado amortiguador este se
descompone en dos devanados en cortocircuito actuando respectivamente en
el eje polar d y en el eje interpolar q, magnéticamente a 90 grados el uno del otro.
14
t e¡. .nI cu¡dr¡¡ur¡
Eicd¡7á6>
FIGURA 2. Modelo de ta Máquina Bipolar ideal (Máquina Generatizada)
15
Por construcción el rotor de una máquina sincrónica presenta dos tipos de
discontinuidad:
2.1.1.1 Discontinuidad Eléctrica. En eleje de los polos o eje directo existen dos
devanados distintos, el devanado de excitación (f) y el de amortiguación
longitudinal d, mientras que en el eje interpolar o transversal sólo existe el
devanado amortiguador transversal q, cuyos parámetros característicos son
diferentes a los del devanado amortiguador longitudinal.
2.1.1.2 Discontinuidad Magnética. Principalmente en la máquina de polos
salientes debido a la permeancia diferente en el entrehierro según los ejes
directo y transversal.
2.1.2 Inductancias Propias del Rotor. Despreciando las perturbaciones locales
en ta distribución del campo magnético en el entrehierro provocadas por la
presencia de las ranuras, elcircuito magnético visto desde el rotores homogéneo
y por to tanto las inductancias propias del devanado de excitación (Lff) devanado
amortiguadoren eje directo (LDD) y devanado amortiguadoren eje en cuadratura
(Loo),son constantes.
Lff = LF (H) Ec 2.1
LDD = LD (H) Ec 2.2
LQQ = LO (H) Ec 2.3
16
2.1.9 Inductancias Mutuas Relativas a los Devanados Rotóricos. En la figura
3 se puede apreciar que la inductancia mutua entre el devanado de excitación
y eldevanado amortiguadoren eje directo (Lro= Lor), no depende de la posición
del rotor por lo tanto:
LFD = LDF = MR (H) Ec 2.4
Elcoeficiente de acoplamiento entre los ejes d y q es cero puesto que están a
90 grados magnéticos y por lo tanto las inductancias mutuas entre eldevanado
amortiguador en cuadratura y los devanados de excitación y amodiguador en
eje directo son nulas.
LFQ = LQF =0 (H)
LDQ= LQD=0(H)
Ec 2.5
Ec 2.6
Las inductancias mutuas entre los arrollamientos estatóricos y rotóricos
definen de nuevo de la figura 3.
- Para 0 = 0o el coeficiente de acoplamiento entre la fase a, y el devanado de
excitación es 1.0 y por lo tanto Laf es máximo.
- Para 0 =[lt2los ejes de la fase a y el eje directo están a 90o, el grado de
acoplamiento es nulo, por lo tanto Laf=O.
17
Axclr¡nsyers¡l ,/ Axe longiturlinrl
,^ ou rJi¡rrt
¡, b, c phrrcr dc l'enroulcmcn¡slrlo¡iouc
f en¡oulcmcnt inducleur
\o
]*
¡mortisrcur longitudinrlrmontStcu ¡ l¡¡nsvc¡s¡l
Da
FIGURA 3. Modelización de los diferentes devanados sde una Máquina Sincrónica
18
- Para 0 = n el acoplamiento vuelve a ser máximo pero en sentido contrario.
Esta variación sugiere un comportamiento cosenoidal de tal forma que:
Laf(e) = Lfa(0) = Mfcos0 (H) Ec 2.7
Ya que Ias fases b y c están 120o y 24Oo retrasadas, entonces:
Lbf(e) = Lfb(e) = Mfcos(o - 2rl3) (H) Ec 2.8
Lcf(e) = Lfc(0) = Mfcos(O+ 2d3) (H) Ec 2.9
Análogamente la inductancia mutua entre los devanados del estator y el
devanado amortiguador en eje directo puede expresarse como:
LaD(e) = LDa(0) = MDcos0 (H) Ec 2.10
LbD(e) = LDb(O) = MDcos(e -2fl3) (H) Ec 2.11
LcD(o) = LDc(o) = MDcos(e +2d3) (H) Ec 2.12
Similarmente para las inductancias mutuas entre los devanados del estatory el
devanado amortiguador en cuadratura podemos decir que:
Para 0 = 0o el acoplamiento entre la fase a y el devanado transversal es nulo.
Pára 0 = rJ2 el acoplamiento es máximo en el mismo sentido.
19
Para 0 =¡¡ el acoplamiento vuelve a ser nulo.
El comportamiento es senoidal y por lo tanto:
LaA(e) = tQa(Q) = MQseng (H) Ec 213
Lbo(e) = LQb(e) = MQsen(o -2rl3) (H) Ec 2.14
Lca(e) = LQc(O) - MQsen(e +2d3) (H) Ec 2.15
2.1.4 Inductancias Estatóricas Propias. En una máquina de rotor cilíndrico,
haciendo abstracción de la posible variación de la permeancia del circuito
rotórico según los ejes d y q debida a la ausencia de ranuras en la zona central
de los polos, el circuito magnético es prácticamente igual y entonces:
Laa = Lbb = Lcc = lh_ + los = cte (H) Ec 2.16
donde los términos:
Lh_ y Los designan respectivamente la inductancia principal (debida al flujo
total), y la inductancia de fuga.
En una máquina de polos salientes, por el contrario, la permeancia del circuito
magnético es máxima en el eje transversal q (figura 4). Sean Ad y Aq los valores
extremos de la permeancia. En la práctica las reglas adoptadas para el
dimensionamiento de las expansiones polares proporcionan una forma tal del
Universidod r,r0n0rn0 de 0ccidcnta
5er:rión 8ibli0f0(0
20
entrehierroconducente a unavariación prácticamentesinusoidalde lapermeancia
entre estos dos extremos.
Si se designa Fa como la Íuerza magnetomotriz resultante de la fase a, sus
proyecciones sobre los ejes d y q según la figura 5 serán:
Fad = Facos0 Ec 2.17
Faq - -Fasen0 Ec 2.18
Los flujos creados por dichas componentes son:
Oad = AdFad =AdFacos0 Ec 2.19
@aq - AqFaq = -AqFasen0 Ec 2.20
Los flujos totales acoplados con la fase se obtienen multiplicando por el número
de bobinas, los flujos precedentes.
Yad = FaNsKwl sAdcos0 Ec 2.21
Yaq - -FaNsKwlsAqsen0 Ec 2.22
El flujo total resultante acoplado con la fase a se obtiene proyectando ad y aq
sobre el eje de la fase a
21
Ya(0) = Yadcosg - Yaqseng
Ya(0) = FaNsKwls(Adcos20+Aqsen20) Ec 2.23
Haciendo:
Lhd_ = FaNsKwl sAd / ia
Lhq_ -FaNsKwlsAq/ia
Ec 2.24
Ec 2.25
La inductancia principal se define como la relación entre los enlazamientos de
flujo y la corriente.
Lh-(0) = Ya(0) lia= Lhd-cos20 + Lhq-sen2o
Lh_(e) = (Lhd_+Lhq_) / 2 + (Lhd_ - Lhq_)cos20/2
Lh-(e) = Lh-+Lhzcos2g Ec 2'26
donde:
Lh_ = (Lhd_ + Lhq_ )/2 Ec 2.27
Lhr=(Lhd_ - Lhq_)12
Según lo que se obtuvo, la inductancia Laa es debida al flujo que cruza el
entrehierro y no incluye el efecto delflujo de dispersión, local (esto es elflujo que
concatena la bobina a pero no cruza el entrehierro).
El flujo localtiene un circuito magnético de geometría constante , por lo tanto
Ls es independiente de , asísu contribución a la inductancia es una constante
y sólo aumenta la magnitud Lh_.
22
FIGURA 4. Circulación del campo estatóricorespectivamente.
lllongitudinal y transversal,
\A¡c rr¡nsycrs|l q
Axc pherc r
\
FIGURA 5. Descomposición de la F.M.M. de la fase a en sus componentesLongitudinal y Transversal.
\er
23
Las inductancias Lbb(e)y Lcc(e) se deducen de la ecuación anterior reemplazando
0 por (g -2112) y (0+2n12) respectivamente.
Laa(O) = Lh-+Ls + Lhrcos2g
Laa (0) = L+ Lhrcos20 = Ls+ Lmcos20 (H)
Lbb(e) = Ls + Lmcos2(0 -2?T/3) (H)
Lcc(0) = Ls + Lmcos2(0+2d3) (H)
Ec 2.28
Ec 2.29
Ec 2.30
La figura 6 muestra la variación de la inductancia propia con el ángulo 0 .
2.1.5 Inductancias Mutuas entre Fases estatóricas. En una máquina de rotor
cilíndrico las inductancias mutuas son constantes, debido a que su configuración
es igual y se tiene:
Lab = Lba = Lbc = Lcb = Lca = Lac = cte (H) Ec 2.31
El acoplamiento magnético entre las fases estatóricas es esencialmente debido
alflujo del entrehierro. Existe otro tipo de acoplamiento a nivel de cíertos flujos
de fuga (de ranura y de cabezade bobina) los cuales se tienen en cuenta en la
expresión numérica de los términos de la inductancia de fuga correspondiente.
25
En estas condiciones se obtiene fácilmente la inductancia mutua Lab entre dos
fases idénticas desfasadas 2ttl3, como
Lab = Lh_cos 2trl3 = -0,5Lh_ Ec 2.32
En la máquina de polos salientes las inductancias mutuas estatóricas dependen
de la posición del rotor. Para calcular la inductancia mutua entre las fases a y
b cuando el eje d está desfasado con respecto al eje de la fase a, se descompone
la fuerza magnetomotriz de la fase a, Fa, en sus componentes según los ejes
d y q como se hizo anteriormente. Los flujos creados por Fad y Faq están dados
por fas ecuaciones 2.21 y 2.22.
Los flujos creados por la fase a y acoplados con la fase b se obtienen
proyectando <Dad y @aq sobre el eje de la fase b ( figura 7).
Oba(O) =Oadcos(g -hú3) - <Daqsen(g -hú3)
oba(e) = FaAdcos0 x cos(O -2Tú3) + FaAqsen x sen(O -tutJg) Ec2.33
26
El flujo total mutuo acoplado con la fase b será:
Yba(0) = NsKwlsAba(g) Ec2.34
puesto que Lba(0)= Yba(0)/ia, teniendo en cuenta las inductancias Lhd-y Lhq-
definidas anteriormente, se obtiene:
Lba(e) = Lab(0) = Lhd_cos0cos(O -2rr/3) + Lhq_sen0sen(0 -21ú3)
Lba(e) = Lab(O) = -(Lhd_+Lhq_)/a + (Lhd_ - Lhq_)/2xcos(20 -2nf3)
Lba(e) = Lab(O) = - 0,5Lh- + Lhrcos(20 -27ú3)
Lba(e) = - Ms - Lmcos2(0 + n/6) (H) Ec 2.35
De la misma manera obtenemos Lbc y Lac quedando:
Lbc(O) = Lcb(O) = - 0,5Lh_+Lhrcos20
Lbc(e) = - Ms - Lmcos2(g - rJ2) (H) Ec 2.36
Lac(e) = Lca(O) = - 0,5Lh_ +Lhrcos(20 +hrl3)
Lac(O) = - Ms - Lmcos2(e +sd6) (H) Ec 2.37
La figura 8 muestra la variación de la inductancia mutua.
27
o$.o-t Phase b
t/
\ 2.,ttOrt -/
\ 0
A.rc phrse l
1\
FIGURA 7. Descomposición del flujo de la fase a en sus componentesLongitudinal y Transversal.
FIGURA 8. Variación de la Inductancia mutua entre fases estatóricas de unaMáquina Sincrónica.
rDrq
\
28
2.2 TRANSFORMACION DE PARK
Como se dijo anteriormente la máquinasincrónica que seestudiaestáconstituida
por tres devanados estatóricos, un devanado de campo y dos devanados
amortiguadores. Estos seis devanados como ya se demostró están
magnéticamente acoplados: el acoplamiento magnético es una función de la
posición del rotor; los enlazamientos de flujo asociados a cada uno de los
bobinados son también una función de la posición del rotor. La tensión
instantánea terminal de cada uno de los devanados es de la forma:
v = tlri +Y Ec 2.39
Donde representa los enlazamientos de flujo, r es la resistencia de los
devanados e ies la corriente con dirección positiva saliendo en elestator hacia
los terminales delgenerador. La notación rL es la suma de todos los términos
apropiados con sus respectivos signos.
La expresión para los voltajes de los devanados es complicada debido a la
variación de Y con la posición del rotor.
Se obtiene una gran simplificación en la descripción del modelo matemático de
la máquina sincrónica al introducir ciertas transformaciones con el objetivo de
29
eliminar la variación de las inductancias y flujos, con la posición del rotor. El
método consiste en hacer un nuevo conjunto de flujos, corrientes y voltajes en
términos de los reales mediante unatransformación que cumpla una restricción
talcomo la de que la potencia debe ser la misma en el modelo real como en el
transformado. Latransformación que cumple estos requisitos es la denominada
transformación de park.
2.2.1 Transformación de los Ejes Móviles a los ejes Fijos. Elnuevo conjunto de
variables es obtenido de la proyección de las variables reales sobre tres ejes;
uno a lo largo deleje del devanado delcampo del rotor denominado eje directo;
un segundo eje a lo largo del eje neutro del devanado de campo denominado
eje en cuadratura y un tercer eje, el estacionario, en la dirección axial de la
máquina.
La sustentación analítica de la escogencia de estos tres ejes se basa en la
distribución espacial de la circulación del campo magnético o tuerza
magnetomotriz. Las fuerzas magnetomotrices originan los flujos y los
enlazamientos de flujo asociados a los devanados. La distribución espacialde
la fuerza magnetomotriz creada por un devanado monofásico puede ser
descompuesta para sus estudio en dos fasores espaciales tales que:
S_ : corresponde a la componente debida a la distribución sinusoidalsituada en
el punto X , Y perpendicular al eje de la máquina y que corresponde
Secrión Bibliofeco
30
respectivamente al eje X al eje directo d, y el eje Y corresponde al eje en
cuadratura q.
S,: corresponde a la componente de distribución uniforme dirigida según eleje
de la máquina ó eje z.
La distribución espacial de lafuerza magnetomotriz y los vectores espaciales
asociados se muestran en las figuras 9 y 10. La componente de distribución
sinusoidal está asociada al flujo principal que atraviesa el entrehierro y la
componente de distribución uniforme se asocia a los flujos de fuga estatórico y
a un flujo axial prácticamente despreciable En máquinas de circuito magnético
simétrico.
La transformación de park se desarrolla matemáticamente de la siguiente
manera:.
Se define que el eje d del rotor está en cualquier instante de tiempo a un ángulo
radianes con respecto al eje de la fase a tal como se muestra en la figura 1 1 . Se
tomarán las corrientes de fase ia, ib, ic saliendo de los terminales del generador.
Si se proyectan estas corrientes a lo largo de los ejes d y q se tienen las
siguíentes relaciones:.
iq = 2/3f iasenO+ibsen(0 -2n /3)+icsen(0 +2n t3)
id = 2l3l ia cos0+ibcos(O -2nt3)+iccos( 0 +2n t3) Ec 2.40
31
FIGUIRA 9. Representación de la distribución espacial de la F.M.M.
FIGURA 10. Descomposición de una onda pulsante en fases espaciales.
d o¡i¡q orir
b oxi¡
ARROLLAMIENTOS:. 3 DE FASE
a,brc. CAMPO
F. 2 DE A.I'vIORTICUACION
D
a
0 (anglc électrique)
FIGURA 11. Proyección de las fases a, b, c en los ejes d, q, 0.
32
Debe recordarse que se tienen tres variables ia, ib, ic y por lo tanto se debe
obtener por la transformación otras tres. Ya se tienen id e iq y la tercera es
proporcional a la corriente de secuencia cero (i0), luego:
I*q = PI"o. Ec 2.41
Donde se definen los vectores de corriente:
Iooo = I-tabc - Ec 2.42
Y donde la transformación de park es definida como:
la
ib
ic
to
id
iq
P= Ug
11 2
cos 0
sen 0
112 1/2
cos(0 -2n 13) cos(O+2r /3)
sen(O -2rc /3) sen(0+2n /3)
Ec 2.43
Elflujo del campo principal se encuentra a lo largo del eje d del rotor e induce
una fuerza electromotriz que atrasa elflujo 90o, por lo tanto la fuerza electromotriz
de la máquina, E, está inicialmente a lo largo del eje q. Si se considera una
máquina que tiene un voltaje terminal V, para el régimen generador el fasor E
33
adelanta al fasor V. El ángulo entre E y V es el ángulo de torque de la máquina
( 6) si el fasor V está en la dirección de la referencia de la fase a.
En t=0 el fasor V está localizado en el eje de la fase a, el eje q está localizado
a un ángulo 6 y el eje d está localizado a un ángulo 0= 6 +n 12.
Para t>0 el eje de referencia está localizado a un ángulo <ont con respecto al eje
de la fase a. El eje d del rotor estará localizado a:
O=ont +6+n /2 (Rad) Ec 2.44
Donde on es la frecuencia nominal angular en Rad/Seg y 6 es el ángulo de
torque en radianes eléctricos .
Pueden escribirse expresiones similares a la ecuación 2.41 para voltajes y
flujos:
Vodq = PVabc Yodq = PYabc Ec 2.45
Si la transformación de la ecuación 2.43 es única, existe una transformación
inversa tal que:
I"o. = P '1 I*o Ec 2'46
34
La inversa de P es:
P'1 = 213
'll 2 cose sene
11 2 cos(e -2n 13) sen(O-2ru /3)
11 2 cos(O +2n 13) sen(0+2r /3)
Ec 2.47
Se nota que p'r = prlo cual signif ica que la transformación P es ortogonal. Siendo
P una transformación ortogonaltambién significa que P es una transformación
de potencia invariante lo cual se esperaba.
P = V€r ia +vb ib + vc ic = Vabct labc
P = (P '1 V0dq)t(P'1 i0dq)
P = VOdg'(P -'¡' P-1lOdq - V0dqtP P-l lOdq Ec 2.48
P = VOdqr= v0 i0 + vd id +vq iq
2.2.2 Ecuaciones de Enlaces de Flujo. Cómo ya se demostró en elcapítulo 2
la máquina tiene seis bobinas mutuamente acopladas, la expresión para los
enlaces de flujo está dada por:
Y=L¡ Ec2.49
35
Entonces se pueden escribir ecuaciones para los seis enlaces de flujo asociados
a los seis devanados, así:
Ya
Yb
Yc
YF
YD
YQ
Laa Lab Lac
Lba Lbb Lbc
Lca Lcb Lcc
LFa LFa LFb
LDa LDb LDc
LQa LQb LQc
LaF LaD LaQ
LbF LbD LbQ
LcF LcD LcQ
LFF LFD LFQ
LDF LDD LDQ
LQf LQD LQQ
Wb Ec 2.50
= inductancias estator - estator
= inductancias estator - rotor
ia
ib
ic
¡F
iD
iQ
Donde
Laa =
LaR = LRa =
Laa Lab Lac
Lba Lbb Lbc
Lca Lcb Lcc
LaF LaD LaQ
Lbf LbD LbQ
LcF LcD LcQ
LFF LFD LFQ
LDF LDD LDQ
LQF LQD LQQ
LRR = = inductancias rotor - rotor
36
Los subíndices F, D y Q hacen referencia a los enlaces, inductancias y corrientes
asociadas a los devanados de campo, amofiiguadoren eje directo y amortiguador
en eje en cuadratura.
De otro lado,
Ljk = inductancia mutua entre los devanados j y k
Ljk = Lkj en todos los casos.
Todas las inductancias varían con respecto a y de acuerdo a la ecuación 2.39
tenemos:
dY/dt=Ldi/dt+dUdti Ec 2.51
La transformación de Park nos evita la complicación que presenta la ecuación
anterior por la variación de las inductancias en el tiempo.
Haciendo una partición matricial de la ecuación 2.50 y aplicando la transformación
de Park tenemos lo siguiente:
P
0
0
u3l;:P
0
0l,.1
*"1YFDo
I
P-1
0
Laa LaR
LRa LRR
2.52
:,ll ,::Ec
37
Donde P es la matriz transformación de park y U. es la matriz unidad 3X3.
Realizando las operaciones indicadas se obtiene:.
Lo
0
0
0
Ld
0
0
0
Lq
KM, KMO
00
0
0
KMo
lo
id
iq
¡f
iD
iQ
(wb)J
0
0
0
KM, o
KMo 0
0 KMo
LF MR
MR LD
00
Ec 2.53
0
0
LQ
En esta ecuación se definen las siguientes nuevas constantes:
Ld = Ls + Ms + (3/2) Lm
L0 = Ls - 2Ms (H)
Lq = Ls + Ms -(3/2)Lm
K=312
De la ecuación 2.53 pueden extraerse importantes apreciaciones:.
- Los enlaces de flujo Yo están completamente desacoplados de los demás
circuitos puesto que la primera fila y la primera columna sólo contienen el
(H)
(H)
38
término de la diagonal. Esta observación es importante ya que las corrientes,
flujos e inductancias de secuencia cero no se presentan en condiciones de
carga balanceadas ó fallas trifásicas como en caso de que trata el presente
proyecto.
- La matriz de inductancias es una matriz de constantes, la transformación de
park remueve la variación en el tiempo como se pretendía.
- La matriz transformada es una matriz simétrica y puede ser físicamente
realizable mediante un circuito eléctrico.
2.2.3 Ecuaciones de Voltaje. Los voltajes para el generador según la ecuación
2.39 corresponde a los circuitos mostrados esquemáticamente en la figura 12.
Se omiten las inductancias mutuas para mayor claridad, pero ellas estarán
siempre presentes con los valores dados en la sección2.1 (inductancias de la
máquina síncrona). Es de notar que las corrientes del estator han sido asumidas
positivas fluyendo fuera de los terminales de la máquina cuando se comporta
como generador. Para las condiciones indicadas se puede escribir la siguiente
ecuación matricial:
V=-ri-Y+V
39
va
vb
vc
-vF
0
0
ra
0
0
0
0
0
00000
rb00000rc000
ia
¡b
ic.--IF
¡D
iQ
Ya
Yb
Yc
YF
YD
Yo
Vn
;-000000
111
111
111
rFo0h00
0
0
ro
(Voltios) Ec 2.54
La contribución del voltaje neutro a Vabc se define como:
Vn=-rn -Ln
-Rn Iabc - Ln Iabc 2.55
Si ra = rb = rc = r cofi'ro es lo usual, se puede definir:
Rabc=rUs C¿ Ec 2.56
1a
ib
ic
111
111
111
(Voltios) Ec
ia
ib
ic
Se puede entonces particionar la Ec 2.54 y reescribirla de la forma siguiente:
40
Vabc Vn
VFDQ
Donde:
I-'FDo -
\I,,- Ec 2,58
Como la ecuación2.57 es complicada por la presencia de coeficientes variantes
con el tiempo en el término Y, se aplica la transformación de park para
eliminarlos. Esto requiere que ambos miembros de la ecuación sean
premultiplicados por:
0
Y-f
YD
Yo
i-r
iD
i^
",""=lí
PO
0u3
41
Por definición para el lado izquierdo de la ecuación:
V0dq
Vroo
(Voltios) Ec 2.59
En el lado derecho, primero se calcula la caída de tensión en la resistencia:
Vabc
Vroo::l
l' ' lln"u"o ll '*| I' o lln"o.
o I le' 'll' 'l I r"u"l
l0 u3ll0 R."oll r,"o I l0 u3l l0 R,ool I 0 u3l
| 0 u3l I r,"o
I
len"o"euo ll'ooql ln"o" o ll'ooql| | | ll ll | (Voltios) Ec2.60
lo t,*l l''". llo R,ool I r,oo I
El segundo término del lado derecho es transformado así:
PO
oU3
PYabc
\y FDO
(Voltios) Ec 2.61
w
42
Para evaluar PYabc se debe recordar la definición de la Ec2.45 de la cual se
puede deducir:
Yodq=PYabc+PYabc
Entonces:
PYabc = Yodq - PYabc = yodq - P P-1yOdq (Voltios) Ec2.62
Se puede demostrar que:
P P-1YOdq = co
000
00-1010
ol
Y
Y
Finalmente eltercertérmino del lado derecho de la ecuación2.57 se transforma
de la siguiente manera:
PO Vn PVn
(Voltios) Ec 2.64
0
-(,)Yq
r¡Yd
(Voltios) Ec 2.63
lo¿q
0 0OU
43
Donde por definición q o* "r la caída de voltaje neutro a tierra en el sistema de
coordenadas 0, d, q, usando la ecuación 2.55
l'looq = PVn = - PRnP-lPlabc - PLnP-rPlabc
t'lo¿q= - PRnP'1I0dq - PLnP'llOdq
l'looq =
3rnio
0
0
3 Lnio
0
0
(Voltios) Ec 2.65
Se observa que la caída de voltaje es únicamente en la secuencia cero situación
de esperarse ya que por el neutro a tierra sólo circula dicha componente de la
corriente.
Substituyendo las ecuaciones 2.59 a2.65 dentro de la ecuación 2.57 se obtiene
lo siguiente:
V0dq Rabc 0 Iooo Y*o PP'l Yodq n0dq
Vroo o Rroo Iroo Yroo
+ +
0 0
(Voltios) Ec 2.66
44
Se nota que todos los términos en esta ecuación son conocidos dados que la
matriz de resistencias es una malriz diagonal y para condiciones balanceadas
el voltaje de secuencia cero es nulo.
to
Db
D+
v^ -0IJ
FIGURA 12. Diagrama esquemático de voltajes de la máquina sincrónica.
+
"t- Lr
Lo+
v- =0v
Lo
Io
-l¡>,F
45
2.2.4 Normalización de las Ecuaciones. Las ecuaciones de flujo, voltaje, torque
y potencia de las secciones precedentes no están expresadas en forma
conveniente para usos en ingeniería. Una dificultad estriba en la inconveniencia
numérica entre los rangos de voltaje del estator usualmente (kilovoltios), y el
rango de voltaje del campo (usualmente voltios). Este problema puede ser
resuelto normalizando a un valor base conveniente y expresando todos los
voltajes en por unidad (p.u) o porcentaje.
Hay muchas maneras de normalizar las ecuaciones y se escogerá una
conveniente para la simulación del comportamiento transitorio de la máquina;
para ello se cuenta con las siguientes guías:
- El sistema de ecuaciones de voltaje debe ser el mismo ya sea que las
ecuaciones estén en p.u o en sistema internacional (Sl) esto quiere decirque las
ecuaciones son siempre simbólicamente las mismas y no se requieren constantes
de normalización .
- El sistema de ecuaciones de potencia debe ser exactamente el mismo ya sea
que las ecuaciones estén en p.u o en Sl, esto significa que la potencia es
invariante al hacer la normalización.
- Todas las inductancias mutuas deben poder representarse por circuitos T
después de la normalización. Este requerimiento se incluye para simplificar la
simulación de las ecuaciones p.u.
46
- La mayoría de las impedancias tradicionalmente suministradas por los
fabricantes deben ser mantenidas en el sistema adoptado para conveniencia de
los usuarios. Otras impedancias p.u pueden ser derivadas fácilmente de los
datos suministrados por el fabricante.
2.2.4.1 Cantidades Base para el Estator.
Potencia base S, = S* = Voltio - amperios nominales por fase.
Voltaje base V, = V = Voltios nominales fase - neutro.
Velocidad base oe = @¡r = Radianes eléctricos/Seg nominales.
De éstos tres datos básicos se pueden derivar las demás cantidades base:.
Corriente base I, = Su / V, - SN / VN (Amperios)
Tiempo base Ts= 1 / tD, = t/co* (Segundos)
Enlaces de flujo Y, = V, T, = V¡r / ol* (Webber - Vueltas)
Inductancias base L, = V, TB / IB = V*/ I* ol* (Henrios)
Resistencia base Rr= Vr/ IB- V * / I, (Voltios)
2.2.4.2 Cantidades Base para el Rotor. Para el rotor no se pueden utilizar las
mismas cantidades bases del estator puesto que los voltios - amperios y voltios
del estator son mucho mayores que los del rotor y por lo tanto las cantidades
47
p.u resultantes serían muy reducidas. Se utiliza entonces elconcepto de iguales
enlazamientos de flujo lo cual significa que una corriente base de excitación o
corriente base de amorliguación en eje directo producirán la misma distribución
espacial de flujo en el entrehierro que la producida por la corriente base del
estator actuando en un devanado ficticio en el eje directo.
Teniendo en cuenta éste concepto y las ecuaciones de enlazamientos de flujo
presentadas en el numeral 2.2.2 se pueden deducir cantidades base del rotor.
- Devanado de campo:
Voltaje base Vr, = K, V, (voltios)
Corriente base Ir, = Ia/K, (Amperios)
Resistencia base Rr, = Kr'R, (Ohmios)
Inductancia base Lra = Kr'L, (Henrios)
- Devanado amortiguador en eje directo
Voltaje base Vo, = KoV, (Voltios)
Corriente base Iou = IB / KD (Amperios)
Resistencia base Ro, = KotR" (Ohmios)
Inductancia base Lo, = Ko'L, (Henrios)
48
- Devanado amortiguador en eje en cuadratura
Voltaje base Vou = Ko V, (Voltios)
Corriente base Ioe = IB / Ko (Amperios)
Resistencia base Ro, = Ko' R, (Ohmios)
Inductancia base Lo"= Ko'L, (Henrios)
Donde:
K, = M, /KMD Ko = M, /KMr Ko = L,o /KMo
Estos parámetros ya fueron dados en el numeral2.2.2.
Algunas de las referencias consultadas presentan indistintamente reactancia
como inductancia y esto es debido al hecho de que en condiciones cercanas a
la nominal co = co* el valor relativo en p.u sería igual a 1, por lo tanto:.
X=rol=L
En lo sucesivo y para estar de acuerdo con la bibliografía utilizada todos las
ecuaciones estarán en p.u y las inductancias tratadas como tales.
49
2.2.5 Ecuaciones de Par y Potencia. Un sistema de potencia estable es aquel
en elcual las máquinas sincrónicas cuando son perturbadas retoman su estado
originalsino hay un cambio neto de potencia, o adquieren un nuevo estado sin
salirse de sincronismo. Usualmente la perturbación causa un transitorio
oscilatorio, pero si el sistema es estable las oscilaciones serán amortiguadas.
El ángulo del rotor medido con respecto a la velocidad de sincronismo y
expresado en grados eléctricos es ta señal que permite medir la estabilidad.
La ecuación de oscilación que gobierna el movimiento del rotor de la máquina,
relacionando el torque de inercia con el resultante de los torques mecánico,
eléctrico y de amortiguación, es igual a;
J0 = Ta ( N.m) Ec 2.67
Donde J es el momento de inercia en Kg.m2 de todas las masas rotatorias
conectadas al eje, 0 es elángulo mecánico del eje en Radianes con respecto
al eje de referencia, y Ta es el torque de aceleración en (N.m) actuando sobre
el eje. Puesto que la máquina a estudiar es un generador, eltorque de entrada
Tm es mecánico y el torque de carga Te es eléctrico. Se puede escribir
entonces:
llnivsrgid,¡d !rrr¡Gmo Oo ¡U¡funlf
Sec';in 8ibri0f6f0
Ta = Tm - Te (N.m) Ec 2.68
50
Ecuación que establece una convención de signo útil, a saber: si Tm es positivo
acelera el eje, mientras que si Te es positivo lo desacelera.
Como ya se vió en este capítulo, la posición del rotor está medida a partir de la
velocidad angular de sincronismo y el eje de referencia de la fase a,
Durante el transitorio aparecen otros torques de los cuales el principal es el
torque de amortiguamiento, resultado de la interacción entre elflujo de secuencia
positiva en el entrehierro y los devanados de amortiguación. En general este
torque es beneficioso a la máquina puesto que tiende a reducir la magnitud de
la oscilación, es decir se comporta como elemento restaurador. El torque de
amortiguación es usualmente asumido proporcional a la frecuencia de
deslizamiento.
Td = Dco (N.m) EC 2.69
Donde D es una constante de amortiguación.
Se puede así completar
manera:.
la ecuación de torque acelerante de la siguiente
Ta = Tm - Te - Td (N.m) Ec 2.70
51
Como ya se demostró en la sección 2.2.1 la transformación de Park mantiene
la potencia invariante y por lo tanto la potencia de salida del generador es:
Pout = va ia + vb ib + vc ic = Vabct Iabc = vabct Iabc
Pout = vodqt lodq = vo i o + vd id + vqiq Ec 2.48
Sustituyendo elvalor de v0, vd y vq de la ecuación matricial 2.66 se obtiene:
Pout =(-r io -Yo - 3 r n io - 3 Ln io) io + (- r id - Yd - olYq) id + (- r iq - Yq+ otYd) iq
Pout = -( ioYo+idYd +iqYq)+co(iqYd -idYq) - r (io' + id2 +iq2) - 3 r n io' - 3 Ln io io
Pout = -[ r (io' + id2 + iq') + 3 r n io2l - l(ioYo + idYd +iqYq) + 3 Lnio io] +
o{iqYd - idYq) Ec 2.71
La energía eléctrica para el estator viene dada por:
dW = Pout dt = -[ r( io2 + id2 + iq' ) + 3 r n io2] dt - [ ioYo+ idYd+ iqYq + 3 Ln io ioJdt
+ <o( iqYd - idYq) dt Ec 2.72
El primer término representa la energía perdida por efecto joule en el estator,
el segundo representa la variación de la energía magnética estatórica y el
tercero representa la energía transformable en energía eléctrica, por lo tanto el
torque de la máquina es:
52
Te = dWd0 = dPout /dW = iqYd - idYq (N.m) Ec2.73
Volviendo a la ecuación de torque acelerante 2.70 y reemplazando todos los
términos por su valor se tiene:
J0 =Tm-iqYd+idYq-D0
Tm - iqYd + idYq = J0 + D0 Ec2.74
2.2.6 Representación Esquemática de los Flujos. La figura 13 esquematiza la
disposición de los arrollamientos en el eje directo y el significado de los
diferentes flujos es el siguiente:
(Dad: representa el f lujo principal acoplado con los tres arrollamientos estatóricos.
Ood, (DoD, OoF: corresponden a los flujos de fuga de los arrollamientos
estatóricos en eje directo d, devanado amortiguador en eje directo D y devanado
de excitación F.
OoDd:designa elflujo mutuo exclusivo entre elestatoryeldevanadoamofiiguador
en eje directo. En una máquina clásica estos flujos son prácticamente
despreciables.
OoDF: representa el flujo mutuo exclusivo entre el devanado amortiguador en
eje directo y la excitación.
53
OoFd:corresponde a un flujo mutuo exclusivo entre eldevanado de excitación
y los devanados estatóricos supuestos todos los demás nulos.
A estos ftujos corresponden las inductancias respectivas:
Lad, Lod, LoD, LoDd, LoDF.
Para condiciones balanceadas, las ecuaciones de flujo se obtienen de la
ecuación matricial 2.53.
Por lo tanto las inductancias propias y mutuas que figuran en estas relaciones
pueden ser explicadas a partir de las diversas inductancias pre-citadas.
Yd = Ld id + KMF iF + KMoio
Yq=Lqiq+KMoio
YD=KMoid+M"ir+LDio
YF=KMF¡d+LFir+M*io
YQ = KMoio+ LQ io
Ld=Lad+LoDd+Lod
LD=Lad+LoDd+LoDF+LoD
LF=Lad+LoDF+LoF
LDd=KMD=Lad+LoDd
Ec 2.75
Ec 2.76
Ec 2.77
Ec 2.78
Ec 2.79
Ec 2.80
Ec 2.81
Ec 2.82
Ec 2.83
54
LDF = MR = Lad + LoDF Ec 2.84
LFd=KMF=Lad Ec 2.85
En el eje transversal el esquema de flujos está ilustrado por la figura 14 en el
cual:
@aq: representa el flujo acoplado entre los arrollamientos estatórico y el
devanado amoñiguador en cuadratura.
Ooq, Qoo: designa los flujos de f uga del devanado estatórico (q) y amoñiguador
(Q), en cuadratura.
Las inductancias correspondientes son designadas por Laq, Loq y L oo, donde
de nuevo refiriéndose a las ecuaciones de flujo estarán representadas así:.
Lq = Laq + Loq Ec 2.86
LQ = Laq + LoQ Ec 2.87
LQq = KMo - Laq Ec 2.88
Introduciendo estos valores en las ecuaciones de flujo 2.75 a 2.79
Yd = Lad (id + io+ ir)+ LoDd (id + ij + Lod id Ec2.89
Yq = Laq ( iq + io) + Loq iq Ec 2.90
55
YD = Lad (id + io+ ir)+ LoDd(id+io)+LoDF(id + ir) + LoDio
YF = Lad ( id + io + ir) = LoDd ( id + ir) + LoF i,
YQ=Laq(iq+io)+LoQio
Ec 2.91
Ec 2.92
Ec 2.93
El resultado de estas relaciones es que los diferentes flujos pueden ser
calculados a partir de esquemas equivalentes obtenidos por un arreglo de las
inductancias definidas según 2.80 a 2.88 y por fuentes de flujo Yd, YD, Y'
Yq, V Yo, (figuras 15 y 16).
En el caso de una máquina clásica para la cual LoDd = 0, el esquema en eje
directo se simplifica según la figura 17.
2.3 PARAMETROS TRANSITORIOS DE LA MAQU¡NA SINCRONICA
2.3.1 Reactancias de la Máquina Sincrónica. Para elestudio de los regímenes
transitorios se definen los elementos de los esquemas equivalentes, cuyas
ecuaciones incluyen las reactancias y las constantes de tiempo. Para obtener
los esquemas equivalentes de la máquina se calculan los valores de las
reactancias a partir del esquema de la figura 16 para el eje en cuadratura y la
figura 17 para la representación en eje directo.
57
FIGURA 15. Esquemas equivalentes para los flujos en el eje directo.
v,, I
v"fro * lr.¡
b
FIGURA 16. Esquemas equivalentes para los flujos en el eje directo.
FIGURA l T.Esquema equivalente usual (LoDd=O) para los flujos en eleje directo.
t,,.l
58
2.3.1.1 Reactancia LongitudinalTransitoria (Xd'). Corresponde a la reactancia
vista en bornes de la figura 18 e incluye las reactancias de dispersión de los
devanados estatóricos y de excitación del generador.
xd'= xod + {xad(xof + xoDf)}/ xad + xof + xoDf Ec 2.93
2.3.1.2 Reactancia longitudinal Subtransitoria (Xd"). Engloba las reactancias
de dispersión de los devanados estatóricos y rotórico del generador, incluyéndose
en la dispersión rotórica la influencia del devanado o de las barras de
amortiguación y la de las piezas macizas del rotor, Se obtiene de interuenirtodos
los elementos de la figura 19.
xd"=xod+{xad [xoDf (xof +xo D) +xof xoD]]/(xad+xoDf ) (xof+xoD) +xofxoD
Ec.2.94
2.3.1.3 ReactanciaTransversalSubtransitoria(Xq").Correspondeala reactancia
vista en bornes de la figura 20 e incluye las reactancias de dispersión de los
devanados rotóricos y la reactancia de reacción de armadura en eje en
cuadratura del generador.
XQ"= Xoq + [xaqxoQ]/xaq + xoQ Ec 2.95
59
RA 18. Definición de Xd'.
trl
FlGU
l'.,xd".GURA 19. Definición de
"l
FI
x.Q
Xq".RA 20. Definición de
Ir;f
I
FIGU
60
2.3.2 Constantes de Tiempo de la Máquina Sincrónica
2.3.2.1 Constante de Tiempo Subtransitoria Longitudinal en Cofto-Circuito
(Td"). Corresponde a la componente de la corriente simétrica de coñocircuito
que se amortigua de forma muy rápida. Depende de las características de
amortiguación de los circuitos rotóricos y, especialmente, del devanado
amortiguador. Esta constante de tiempo se obtiene del circuito de la figura 21.
Td"= 1 /ro{xoD+[xadxodxof + (xad+xod) xof xoDf ]/(xad+xof +xoDf) xod+ (xof +xoDf )xadl
Td" = TdO" xd"/xd' Ec 2.96
2.3.2.2 Constante de Tiempo Subtransitoria Transversal en Corto- Circuito
(Tq"). Corresponde a la componente de la corriente simétrica de cortocircuito,
que se amortigua de forma muy rápida. Depende de las reactancias de
dispersión de los devanados rotóricos y de la reactancia de reacción de
armadura en eje en cuadratura delgenerador. Se obtiene delcircuito de la figura
22
Tq"= 1 /ro(xoQ+(xaqxoq/xaq+xoq))
Tq" = Tqo" xq"/Xq Ec 2.97
2.3.2.3 Constante de Tiempo Transitoria Longitudinal en Corto-Circuito (Td').
Corresponde a la componente de la corriente simétrica de cortocircuito que se
amortigua de forma considerablemente más lenta y se calcula desarrollando las
61
ecuaciones que gobiernan el circuito de la figura 23.
Td'= 1 /rf{xof+xoDf+(xadxod)/(xad+xod)}
Td' = TdO'xd'/xd Ec 2.98
2.3.2.4 Constante de Tiempo Subtransitoria Transversal en Circuito-Abierto
(TqO"). Incluye la reactancia de dispersión del devanado rotórico y la reactancia
de reacción de armadura en eje en cuadratura. Corresponde a la definición
siguiente:.
TqO" = 1/ro (xag + xoD) Ec 2.99
2.3.2.5 Constante de Tiempo Transitoria Longitudinala Cicuito Abierto (TdO').
Corresponde a la modificación de la fuerza electromotriz en vacío de dicho
generador al variar súbitamente la tensión de excitación, se representa de la
forma siguiente
TdO'=Tf=xflrt Ec 2.100
2.3.2.6 Constante de Tiempo Subtransitoria Longitudinal a Circuito-Abiefto
(Td0"). Corresponde a las reactancias de dispersión de los devanados rotóricos
y estatóricos del generador. Se obtiene de inte¡venir todos los elementos de la
tigura 24.
TdO"= 1/ro{xoD+[(xad+xoDf)xof/(xad+xoDf+xof)]] Ec2.101
3. DIAGRAMAS DE BLOQUE DE LA MAQUINA SINCRONICATRIFASICA
Las ecuaciones de tensión y de flujo establecidas en la sección anterior
completadas por las ecuaciones de pary de potencia definen elfuncionamiento
de la máquina en un régimen cualquiera. El sistema de ecuaciones así obtenido
presenta una no - linealidad producida por el producto rodYq y su solución debe
efectuarse por métodos numéricos o analógicos. En los regímenes de la
máquina sincrónica a velocidad constante puede presentarse:.
- Funcionamiento cuasiestacionario en régimen de corriente senoidaltalcomo
e! arranque asincrónico.
- Funcionamiento en régimen transitorio tal como el caso de corto circuito
brusco.
La linealización del sistema es posible en los casos siguientes:.
- Oscilaciones periódicas de poca amplitud de la rueda polar alrededor de su
posición de equilibrio.
65
- Variación brusca de poca amplitud de alguna cantidad eléctrica o mecánica.
3.1 CORTOCIRCUITO BRUSCO DE UNA MAQUINA SINCRONICA
CONECTADA A UN BARRAJE INFINITO.
3.1 .1 Expresión normalizada para la corriente de Armadura en Ejes Directo y en
Cuadratura. Considérese el sistema de la figura 25 donde una máquina
sincrónica está conectada a un barraje infinito a través de una línea de
transmisión la cual tiene una resistencia Re e inductancia Le. Se muestran
únicamente los voltajes y corrientes por fase, asumiendo que no existe
acoplamiento entre fases por inspección se puede escribir:.
FIGURA 25. Generador sincrónico conectado a un barraje infinito.
va
vb
vc
Vá"o
vb""
VCo
+Re Ú
ta
¡b
ic
+ LeU
ia
¡b
ic
Ec 3.1
66
En notación matricial tenemos:
Vabc = V""abc + ReUIabc + LeUIabc Ec 3.2
El cual al aplicarle la transformación de Park toma la forma:.
Vodq = PVabc = PV""abc + Relodq + Le Plabc (p.u) Ec 3.3
El primer término del lado derecho se llamará V*OdQ y su valor puede ser
determinado asumiendo que el voltaje del barraje infinito V abc es trifásico
balanceado.
V,.abc = 2Y Ec 3.4
Donde V es la magnitud r.m.s del voltaje de fase. Usando las identidades
trigonométricas y teniendo en cuenta que 0 = orut +6+n12 se puede demostrar
que:.
cos ( co*t)
cos ( rort - 120" )
cos ( ro*t + 120" )
67
V."Qdg=PV..abC=3Voo
0
- sen6
cos6
Ec 3.5
El último término del lado derecho de la ecuación 3.3 puede ser computado de
la siguiente manera, a partir de la definición de transformación de Park:.
IOdq = Plabc, derivando Odq - Plabc + Plabc entonces;.
Plabc = Iodq - Plabc = I0dq - PP'1 lOdq Ec 3.6
Donde PP-1 es una matriz conocida (ecuación 2.63), por lo tanto la ecuación 3.3
puede escribirse como:.
+RelOdq + Lel0dq -olle (p.u) Ec 3.7
Llamando 3V- = V* y rearreglando la ecuación se tiene:
Vo = Re io + Le di/dt Ec 3.8
vd = -Vxsen6 + Reid + Le did/dt + coleiq Ec 3.9
VQ = VN cos6 + Reíq + Lediq/dt - roLeid Ec 3.10
I
I'VOdq =tu-l - sen6
I cos6
0
-iq
id
68
Transformando en Laplace y teniendo en cuenta que para condiciones trifásicas
balanceados Io = 0 tenemos lo siguiente:
Vd + V*sen6 = Re Id + LeS Id + cole Iq Ec 3.11
Vq - V*cos6 = Relq + LeSIq + oleld Ec 3.12
Haciendo -V*sen6 = V*d y Vrcos8 - V*Q se puede escribir:
vd - V*d (Re + LeS) ole
Ec 3.13
Vq - V*g - cole (Re + LeS
Despejando Id e Iq y haciendo re = Le /Re se tiene lo siguiente:
Y1(S) 1 Vd . VNd
= Xe/DET Ec 3.14
- 1 Y1(S) Vq - VNq
Donde:
Y1(S) = (Re/Le) (1+S te)
Xe = ole = Ls, puesto QUO co= 1 p.u
DET(S) = Re2(1 + Ste)2 +Xe2 te = Le/Re
Id
Iq
Id
Iq
69
3.1.1.2 Expresión normalizada para las Tensiones y Flujos. De la expresión
matricial 2.66 y teniendo en cuenta condiciones trifásicas balanceadas (no
existe la secuencia cero ni la impedancia de puesta a tierra), se obtiene:.
vd=-rid+dYd/dt
Vg=-riq+dYd/dt
vF = rf if + dY.ldt
0 = Io io + dYo/dt
0 = ro io + dYo/dt
(l)Yq
olYd
Ec 3.15
Ec 3.1 6
Ec 3.17
Ec 3.18
Ec 3.19
Ec 3.20
Ec 3.21
Ec 3.22
Ec 3,23
Ec 3.24
Transformando en Laplace se tiene:.
Vd=-RId+SYd- roYq
Vq=-RIq+SYq+coYd
Vf=Rrl,+SY,
0=Rolo+SYo
0 = Rolo + SYo
De las ecuaciones de f lujo 2.89 a 2.93, con L oDd =0 y transformando en Laplace
se tiene lo siguiente:.
Univcrsidod uronomo de 0ccidente
Serc;ón Biblrotsco
70
Yd = (Lad + Lod)d + Ladlo+ LadI,
Yd = LdId + Ladlo + LadI, Ec 3.25
Yq=(Laq+Loq)lq+Laqlo
Yq + LqIq + Laqlo Ec 3.26
Yo = Ladld + (Lad+Loor+Loo)Io + (Lad+ Loor)I,
Yo = Ladld + Lolo + LorI, Ec 3'27
Yr = Ladld+(Lad + Loor)lo + (Lad+Loor+ Lor)I,
Ye = Ladld + LoFIo + Lrl, Ec 3.28
Yo=Laqlq+(Laq+Loq)lo
Yo = Laqlq + Lolo Ec 3'29
De fas ecuaciones 3.23 y 3.27 se deduce:
0 = S Ladld + (Ro + S Lj Io + SLorl, Ec 3.30
De la misma forma a partir de las ecuaciones 3.22y 3.28 se obtiene:.
Vr = SLadId + SLorlo+ (R, + SLr)I, Ec 3.31
Se obtiene así un sistema de ecuaciones con el siguiente arreglo en forma
matricial:
71
vd
0
VF
Ld Lad Lad
SLad (Ro+SLo) SLDF
SLad SLo, (R, + SLr)
Id
ID
IF
Ec 3.32
Para obtener la solución de id se resuelve la ecuación matricial 3.32 :
Yd Lad Lad
0 (Ro+ SLo) SLo,
VF SLo, (R, + SLr)
Id= Ec 3.33
Ld Lad Lad
SLad (Ro + Lj SLo,
SLad SLo, (R, + SLr)
Haciendo los cálculos y teniendo en cuenta que:.
Lo= Lad + Lod i Lor = Lad + Loo, I LD = Lad + Loor+ Loo ;
Le = Lad = Loo, + Lo, obtenemos:.
72
Yd(1 +st'dj( 1+sr"do) -(1+sroo) (Lad/Rr)V,
ld= Ec 3.34
LD(1+St'd)( 1+Sr"d)
Ld(1+St'd)(1+Std")ld (1+Stojco*T'ad
Yd=- Ec 3.35
(1+t'd)(1+Sr"do) (1+sr'do)(1+St"do)
r'do = LF / RF = Constante de tiempo transitoria longitudinal en circuito abierto.
r"do = 1 / RDI Loo + ( Lad+Loor)Lo"/ Lad+Loor+Lor] = Constante de tiempo
subtransitoria longitudinal en circuito abierto.
r'd = 1l Rr(Lo, + Loo, + (Ladlod)/(Lad+Lod) = Constante de tiempo transitoria
longitudinal en corto circuito
1
T"d = Lod +
LadLod Lor+( Lad+Lod) LorLoo,
RD (Lad+Lor+Loor)Lod+(Lor+Loor)Lad
Constante de tiempo subtransitoria longitudinal en corto circuito.
73
T'ad = Lad//R, = Constante de tiempo principal longitudinal
roD = L oo/Ro = Constante de tiempo de dispersión del devanado amortiguador
en eje directo
Yd = Ld(s)ld + G(s)V, Ec 3.36
Donde:
Ld(s) = Ld [(1+St'd)(1+Std") / (1+St'do)(1+St"d)l
G(s) = (1+ Stod)co*t'ad / (1+St'dj(1+Sr"do)
Para el eje transversal se combinan las ecuaciones 3.24 y 3.29 de la misma
forma
IQ = -lqSlaq / RQ + S(Laq + LoQ) Ec 3.37
Introduciendo este valor en la ecuación 3.26 se tiene:.
Yq=(Laq+Loq)[RQ+S(LoQ+(LaqLoq/LaqLoq))]/Ra+S(Laq+LoQ)Iq Ec 3.38
Yq = Lq [(1+ Sd'q)/(l +St"qo)]lq Ec 3.39
74
Donde:
r"go = t/Ro (Lag + LoQ): Constante de tiempo subtransitoria transversal en
circuito abierto.
r"q = 1/Ro ( LoQ + (LaqLoq)/Laq + Loq): Constante de tiempo subtransitoria
transversal en coñocircuito.
Yq = Lq(s)lq Ec 3.40
Donde:
Lq(s) =[Lq(1 +S t"q)/(1 +S r"qo]
3.1 .1 .3 Expresión normalizada para la Corriente de Excitación. Combinando las
ecuaciones 3.22 y 3.28 se obtiene:
V, = Rrl, + S/ro* [Ladld + Lorl o + Lrlrl Ec 3.41
Despejando Io de la ecuación 3.25 se obtiene
Io = 1/Lad [Yd - Ldld - Ladlrl Ec 3.42
75
Reemplazando elvalor de Yd de la ecuación 3.36 en la ecuación3.42tenemos:
ID = 1/Lad ILd(s)ld + G(s)V, - Ldld - Ladlrl Ec 3.43
Reemplazando el valor de ID dado por la ecuación 3.43 en la ecuación 3.41
teniendo en cuenta los valores de las inductancias dadas por las ecuaciones
2.80 a 2.85 se obtiene una ecuación en VF, IF e Id que expresa la dependencia
entre la corriente de campo y la tensión de excitación y de la corriente de
armadura en eje directo.
I, = [(1/Rr¡(1 +Sto)/(1+St'do)(1 +St"do)]VF+G(s)Id Ec 3.44
Donde:
ro = 1/Ro(Lad+Loor+Loo): Constante de tiempo para el devanado amortiguador
en eje directo.
I, = F(s)V, + G(s)ld Ec 3.45
Donde:
F(s) = (1 /RFX1 +Srj/(1 +St'd)(1 +St"do)
G(s): Tiene el mismo valor de la expresión de los enlaces de flujo en eje directo
de la ecuación 3.36.
76
3.1.1.4 Expresión normalizada para la Ecuación de Par. La ecuación2.74toma
después de normalizada la forma siguiente:
Tm-iqYd+idYq=M0+D0
Donde:
M = 2olrvH
H = coeficiente de inercia en p.u
D = coeficiente de amortiguamiento del generador
La ecuación 2.44 relacionaba el ángulo de giro del rotor con el ángulo de par
0=orvt + 6+n 12 Ec 2.44
Normalizando la ecuación ( o* = 1 p.u) se tiene:
0=t+6+n/2 Ec3.46
77
Derivando:
de/dt = 1+d6/dt
s0=1+s6
d2e /dt - d26/dt2
S20 =s26
Ec 3.47
Ec 3.48
Ec 3.49
Ec 3.50
De la transformación de Laplace para la ecuación 2.74 se obtiene:
Tm - IqYd+ldYq =JS20+DS0 Ec 3.51
Reemplazando S'0 y S0 por su valor en función de 6 se tiene:
Tm - IqYd+ldYq = JS26 + D(1+56)
Tm-IqYd+ IdYq- D=(JS2+DS)6 Ec3.52
3.2 DIAGRAMA DE BLOQUE PARA LA MAQUINA SIN SI\TURACION
Se hará un resumen de las ecuaciones necesarias para la elaboración del
diagrama de bloque.
78
-'l
o-l
Vd
Vq
Id
Iq
S -(1+56)
(1+SD) S
Ec 3.53
Ec 3.54
Ec 3.55
Yd
Yq
Y1(s) 1
= y2(s)
-1 Yl(s)
Donde:
Yr(s) = Xe/DET(s)
Yd Ld(s) 0
Yq 0 Lq(s)
Id
Iq
vd -V*d
Id
Iq
+
Vq -VNq
G(s)
VF
0
IF=F(s)Vr+SG(s)ld
E= [1/(JS2 + DS)][Tm - IqYd+ IdYq - D] Ec 3.56
Cómo éste es un modelo de múltiPles
una de ellas.
79
entradas y salidas estableceremos cada
3.2.1 Variables de Entrada
3.2.1.1 Tensión del Barraje Infinito V*
VNd = -VNsen6 Ec 3.57
VNq = VNcos6 Ec 3.58
3.2.1.2 Tensión de Excitación Vr. Corresponde a una señal de referencia
controlada por el respectivo regulador de tensión de la máquina.
3.2.1.3 Torque Mecánico Tm. Correspondiente a una señal de referencia
controlada por el respectivo regulador de frecuencia de la turbina.
3.2.2 Variables de Salida
3.2.2.1 Tensión lnducida en el generador Vlgv
..:LrniVe¡SIUü0 ! r1.0m0 dB
(s¡r;f¡ EibliotecO
V =WdTqT ev = tan-r VqA/d Ec 3.59
80
3.2.2.2 Corriente de Salida del Generador u
I = { (Id2 + Iq') gr = tan-t Iq / Id Ec 3.60
9.2.2.9 Corriente de Excitación Ir. Corriente que origina una fuerza magnetomotriz
de campo, la cual produce el flujo de la máquina'
4. APLICACION DEL EMTP PARA LA SIMULACION DE
TRANSITORIOS DE LA MAQUINA SINCRONICA
Para lograr el objetivo de éste proyecto, simular los fenómenos transitorios de
la máquina sincrónica, se utilizó como herramienta para elanálisis el programa
EMTP. Como aplicación específica se utilizaron los parámetros que a continuación
se listan, los cuales corresponde a uno de los generadores trifásicos de la
Central Hidroeléctrica del ALTO ANCHICAYA, diseñados y construidos por la
firma CANADIAN GENERAL ELECTRIC - CGE.
S=126MVA V=13,8kV n=450rpm f=60H2
Ra=0,2448"/" Xd=106% Xq=86% Xd',=250/o
Xq'=86% Xd"=1 8,3o/o Xq"=1 8,3o/o* Xl=7,6o/o
T'd0=6,7seg T'qo=1seg* Td0"=0,045segTqo"=0,055se9'
X0=87o Xn=o,oool% lf=1604Amp np=1opolos
* Datos asumidos para ésta máquina, tomados de la referencia Jean Chatelain.
82
Para dicha implementación se tiene en cuenta la siguiente utilización básica del
EMTP.
4.1 ESTRUCTURA DE LOS DATOS DE ENTRADA
Los datos de entrada al EMTP se suministran al programa en "tarjetas", a
continuación se sigue con la denominación de "tarjetas" que deben constituir el
fichero de entrada:.
- Tarjeta BEGIN NEW DATA CASE
- Tarjeta de datos misceláneos reales
- Tarjeta de datos misceláneos enteros
- Tarjetas de interruptores
- Tarjetas de fuentes
- Tarjeta de los parámetros de la máquina sincrónica trifásica
- Tarjetas de requerimientos de salida
- Tarjetas de requerimientos en impresora
La estructura de los datos de cada tipo de "tarjeta" se presenta en la edición que
se realiza para poderejecutarelprograma EMTP (veranexo 1)y de ésta manera
obtener la simulación que se pretende en este proyecto.
83
Los pasos que se deben de efectuar para poder utilizar el programa EMTP son
los siguientes:.
1. Encender el computador
2. Salir al directorio raíz del disco duro ( C:\ )
3. LLamar un editor de texto como el editor del Norton y crear un archivo para
la entrada de datos al EMTP. Dicha edición consta de una serie de tarjetas las
cuales se utilizan para podersimularelcaso que se requiera. Cuando setermina
la edición ésta se graba en el directorio: C:\atpsim\lec. Este directorio se crea
automáticamente cuando se instala el programa EMTP y en el aparecen los
programas tpl.exe y pcplot.exe.
4. Llamar el directorio de atpsim\lec pulsando cd\atpsim\lec
5. Pulsar enter
6. Pulsar tpl.exe
7. Pulsar enter
8. Pulsar el nombre del directorio que se creó en la edición
9. Pulsar enter
10. De ésta manera el programa EMTP corre. Cuando no aparece ningún error
el programa crea un número el cual consta de la fecha y de la hora exacta en
que comenzó a correr el programa
11. Cuando termina la ejecución del programa se pulsa la palabra STOP
12. Pulsar enter.
84
Para obtener la gráfica del caso simulado se siguen los siguientes pasos:
13. Pulsar pcplot.exe
14. Pulsar enter
15. Pulsar F1
16. Pulsar enter
17. Pulsar enter
18. Pulsar el número creado, explicado anteriormente
19. Pulsar enter
20. Pulsar enter
Aparece un cuadro donde están todas las variables que se pueden graficar el
programa pcplot.exe tiene tres formas de graficar:
1- tiempo en función
2- xy - plot
3- respuesta en frecuencia
21 . pulsar el número uno (1) que corresponde a graficar el tiempo en el eje de
las abcsisas
22. Pulsar el número de la variable que se requiera graficar. Solamente dibuja
hasta tres variables
23. Pulsar enter.
De esta manera se obtiene la gráfica que se quiere obtener
24. Terminada la gráfica se pulsa enter
25. El programa pregunta si se quiere confirmar dicha gráfica. En caso
85
afirmativo se pulsa la letra (Y)
26. Para retirarse del programa pcplot.exe se pulsa (ALT+ F10).
El computador está listo para correr otro caso.
5. CONCLUSIONES
1. Debido a la asimetría magnética de la máquina sincrónica de polos salientes,
la determinación de los circuitos equivalentes debe hacerse tanto en el eje
directo como en el eje en cuadratura.
2. Las inductancias de la máquina sincrónica varían en el tiempo y también con
la posición del rotor, esto complica el desarrollo de las ecuaciones de dicha
máquina y su solución se obtiene aplicando la transformación de Park a las
ecuaciones matriciales.
3. La inductancia transitoria del eje en cuadratura es aproximadamente igual a
la inductancia sincrona del mismo eje.
4. Para describir la parte eléctrica de la máquina sincrónica es necesario
conocer un conjunto de reactancias tanto propias como mutuas.
5. Las magnítudes que interesan en el estudio de transitorios en la máquina
sincrónica son las tensiones y corrientes del estator y la tensión de excitación,
87
pues son las variables que realmente se tienen en cuenta de bornes hacia fuera.
6. Para los casos objeto de este estudio se consideraron condiciones
balanceadas, con lo cual la parte mecánica no afecta los resultados y por
consiguiente se omiten en las formulaciones debido a que el generador gira a
velocidad sincrónica.
7. Enel modelo desarrollado para la máquina sincrónica no se consideraron los
efectos de la saturación ya que la no linealidad requiere de herramientas de
análisis como el método de los elementos finitos.
8. Antes de la falla en el generador solamente existen voltajes y corrientes de
corriente alterna, en tanto que después de la falla se encuentran corrientes
alternas y continuas.
9. La corriente subtransitoria que es el primer período de la corriente de falla es
limitada por el arrollamiento amortiguador del generador sincrónico.
10. La corriente transitoria conlleva la contribución de una componente de
corriente continua del circuito de campo en el instante del cortocircuito.
88
BIBLIOGRAFIA
ADKINS Bernard . Teoría general de las máquinas eléctricas. EdicionesUrmo, Bilbao, 1957.
ANDERSON P,M; FUAD A.A. Power system control and stability Volume l,lowa State University press, First edition , 1977.
BRANDWAJN V. EMTP Theory book, Three Phase Synchronous Machine.Agosto 1986.
CHAPMAN Stephen J. Máquinas Eléctricas. Ediciones Mc Graw Hill 1988.
CHATELAIN Jean. Machines Electriques Volume X, Editions Georgi, Premiére.1 983.
EMTP. Rule Book Julio 1987.
FITZGEHAL A.E, KINGSLEY Alexander Kusko. Teoría y análisis de las
Máquinas Eléctricas, segunda edición
GINSBERG. Desing Calculations for A-C Generators AIEE transactions Volume69, 1 950.
KIRSCHBAUM. Per -Unit Inductances of Sinchronous Machine Volume 69.AIEE transactions 1950.
KRON. Reference Frames of a Synchronous Machine AIEE transactionsVolume 69. 1950.
LIWSCHITZ Michael, WHIPPLE Clyde. Máquinas de corriente alterna,EditorialCecsa. 1970.
89
MATSCH Leander. Máquinas Electromagnéticas y electromecánicas. 1972.
Mc PHERSON. George. Introducción a Máquinas Eléctricas y
Transformadores. Editoriall Limusa. 1987.
ROEPER Richard. Conientes de cortocircuito en reded trifásicas EdicionesMarcombo - España. 1985.
STEVENSON William D. Análisis de Sistemas Eléctricos de Potencia.Ediciones Mc Graw Hill Segunda Edición 1988.
A¡IE){I] I. IDICIÍ}I{ PAftA SII'IULíiR LIHA FALLA TRIFASIC/i.
BEEIN I{EI I}ATA CASE
PFII{TED NUIIBER lllI}TH, 1I, i, ( Request nexitur precrsion for I autput colunn:
.00út00 .90t' ó0. ó0.
!1111-ll5 5 2u 'til 100 ti'|:l
c
C SECCI{}}¡ES PI [E UIIA LI}IEA
I flA[Al RJüA 9.04 41.88 ]51.28
? ilA[81 RJ0B 7.94 19.0{ -57.6 ?,2 41.5t 360.88
3 ttAccl RJüC 7.6S l{,9É -14,9á7.8{ 19.04 -$7.É 9.04 44.8S 351.28
tBLA|ll( rard endinq ¡ll branch cards
c
C IIITERRUPT|]R I}E SIITULACI{}I'| I)E FALLA DE LA FASE A A TIERRA
IIACAT .05 1ü, 1
HACBI ,ü5 10, tl'lAcc! .05 10. I
FLAlll( rard ending all sritch cerds
59 trACAr 9340.0 ú0.0 35.0
ltA[81
itAccl
TOLERAI{CES ?{'
PAfiA}iETER FITTII{T 2.0I I tót. 1. 1?ó. 13.8 1ó04,
BLAI{|( c¡rd for the 0-axrs.0ü24480 .ü7ó 1,06 0,8ó .25 0.8ó .183 .t831.7 1.0 ,ü45 .055 .08 .0t'('ü01 .00üü01
I 1.ü 38.056 .0ü01
FLAtll( card endinq rass crrdsql
ll
?l3t
FLA|I|( card endinq output specilications for this 5.1,l.
F IN ISH
BLAllf card ending all source c¿rdsñACAI llACBt }IACC1
BLAH|( card ending node voltage output requestPRINTEF PLOT
IC Titulo de la gralicer 1-7C
c
? TRAI{SITORIOS I}E LA IIABUINA SINCROIIICA TRIFASICA
'ron ,. o,o ,o nAcAr n*cBr nAC[l
BLAHI( CARI} ENDIIIE PLI}T CRRDS
BEG¡II IIE}I DATA IASE
FLAI{I(
ANEXO 2. LISTADO Y ANALISIS DE LAS GRAFICAS QUE CORRESPONDEN
A UNA FALLA TRIFASICA.
la, lb, lc: (Corrientes de las fases a,b,c en falla).
Después de que ocurre la falla, por las fases del Generador Sincrónico se
encuentran corrientes alternas y continuas, estas corrientes continuas aparecen
en el Generador Sincrónico debido a que éste es básicamente inductivo.
En 50ms tiempo en que ocurre la falla, la componente alterna de la corriente
salta a un valor muy alto pero la corriente total no cambia en ese instante.
La componente de corriente continua es suficientemente grande para que la
suma de las componentes de corriente alterna y continua justamente después.
de la falla sea igual a la corriente alterna que circulaba antes de la falla.
Se observa que la magnitud de la componente de corriente continua es
diferente en cada una de las fases, también se visualiza los tres períodos de la
componente de falla (período subtransitorio, período transitorio, período
permanente).
La corriente subtransitoria que es el primer período de la corriente de falla es
limitada por el arrollamiento amortiguador del generador sincrónico.
La corriente transitoria conlleva la contribución de una componente de corriente
continua del circuito de campo en el instante del coftocircuito.
lf (Corriente de Campo).
La componente de corriente alterna de la corriente de campo de debe a la
componente de corriente continua de la corriente de armadura y decae a la
misma velocidad, y está determinada por la resistencia de armadura.
En elestado transitorio, la componente de corriente continua de la corriente de
campo, y la componente de la corriente alterna de la corriente de armadura
decaen juntas a una determinada velocidad debido a la resistencia del devanado
de campo.
lg (Corriente de falla en el Rotor).
Señal de componente de cuadratura de falla en el rotor, señal períodica
amortiguada, no existe el fenómeno de la modulación puesto que solamente
ocurre en fallas asimétricas.
ID, lo, lo (Corrientes de Armadura en el eje d,q,O).
lo, es una corriente negativa amortiguada.
lo, es una corriente simétrica amortiguada.
lo, antes de la falla aparece una señalcon una frecuencia de 85.76H2, en 50ms
tiempo en que ocurre la falla esta decae exponencialmente.
lKD, IKQ (Corriente en el devanado amortiguador D y Q respectivamente).
La componente en cuadratura es igual a cero debido a que no hay componente
activa.
La componente en eje directo es una componente reactiva (inductiva) ésta es
la causante de la desmagnetización de la máquina.
MForce (Fuerza Magnetomotríz total en el entrehierro).
En el período subtransitorio existe una componente aperiódica que varía en
forma aperiódica, la existencia de los picos grandes siendo pequeños los
coeficientes de atenuación indica que fuerzas pedurbadoras de pequeña
magnitud provoca oscilaciones grandes en el rotor de la máquina sincrónica,
ésta f uerza magnetomotríz varía en forma exponencialy se representa en forma
de oscilaciones armónicas amortiguadas.
Alternative Transients Progran {ATP), IElt FC 11 trenslation, Copvright i?87. Use licensed only through LEC (|(.U. Leuven, Belgiue),
il¡te {dd-nth-yyi and ti¡e ol day {hh.m.ss} = 03-llar-93 t3:31:29 Naoe of disl plot file. if any, is [:33031331.p14
For infornetion, rsnsult the roprrighted ATP EIITP Rule Boot published by LEC in Julyr 1987. L¡st rajor Progrer update: llayr l9BB
Total length of "LABC0tl" tables = 48797 IIITEGER rords. '\IARDIII' List Sizes lollor : 214 211 42f 12 212{
33 594 1487 6s 135 42 50 I t6 I 34 3 5 lÉüü t&0 100 150 3599 3 400 42 42
-----------+-¡)escriptive interpretation of input deta cards. i Input d¡t¡ c¡rd images are shorn belor, all 80 colulns, chararter bv charecter
012345678012345ó7890t 23tfó7t901 23{-5ó7890123f5678t0123f5ü78901231567890123156799012345t7890
.----+-------
Co¡¡ent card, f{UltDEI} = f582, lI d¡t¡:llA0AC.l]AI
ll¿rker cerd preceding nex E|IIP data cese, |BEEIH HEt I)ATA IASE
|lidth of tire-step loop nurbers. lf=13 S=2 |PRII{TED }lUltBER HIDTHT 13, 2, { Request a¡xirul precision for I output co¡u¡ns
llisc. d¿te. 2.000E-0{ 9.0008-01 6.000E+01 i .000200 .900 60. 60.
ñisc.dat¿. t t I I 1-l 0 I 0 0 i I I I I t -t Ittl [arning. l{onzero I0PT or C0PT dilfers fror the porer frequency of 5,00000000E+01 Hz. l}id the user reellv rant this?
Printout : 5 5 20 20 100 100 i 5 5 20 20 100 100
Coment c¡rd. ilüllDEI} = 1588,
[o¡rent card. l{UllDCI) = 5589.
Co¡lent card. llUllllCD = 5596.
Coment card. llUllDCD = 5597.
IF
:C SECEIÍ]I{ES PI DE U}IA LII{EA
Ist of Pl-ckt, 9.0f0E+00 l.l90E-01 9.3tgE-07 i I ltACAl RJ0A 9.04 {4.88 J51.28
7.g1gE+00 5,051E-02-1.528E-07 9.?00E+00 r.t82E-011 2 ltACBl RJoB 7.84 19.04 -Í7,6 9.2 {f.5É 360.88
7.ó80E+00 3.9ó88-0?-3.968E-09 7,8{0E+00 5.051E-02i 3 ilACil RJoC 7,ó8 11,96 -l{.9ú7.84 19.0{ -57.6 9.04 f4,88 351.28
Co¡rent card. NUIIDID = 5593.
[oment card. llU]'ll][I} = 5594.
El¿nt card ending branches. IBR, t{T0T = 3 7 |BLAI{K IARD E]IDING ALL BRAHCI{ CARDS
lcrf
lcIC I}ITERRI]PTOR DE SIITULACIOil OE FALLA DE LA FASE A A TIEHRA
Sritch. 5,00E-02 1.¡Sl+01 0.00E+00 0.00E+00 i ttACAl .05 10,
Sritch, 5.00E-02 t.00E+01 0.00E+tr0 g,¡1rt+00 I ñAc81 .05 10.
sritch. 5,00E-02 ¡.$st+01 0,00E+00 0.t)0E+00 i lfAccl .05 10,
Blant card ending sritches. I(SITCH = 3. lBLAtlK CARI) EilDItlB ALL SflI[H CARDS
Cor¡ent c¡rd. llUltI)CD = 5É02.IF
Source, 9.34E+03 ú.(10E+01 3,50E+01 -9.99E+03 159 llACAl 9340,ú 60.0 35'0
Znd ohese ot S.ll, 9,34q6ffi+03 -85.00 I llACEl
3rd phase ol S.ll, 9.34000E+03 -205,00 I ItAECI
Epsilon. 1,00E+02 t.00E-15 1.00E-16 20 |T0LERAIiCES 20
0ptirire Park'E date. 1,000E+00 |PARAIIEIER FlfTlllB 2.0
4th S.tl. c¡rd, t I 0 ló 1.000 1,000 i t I 16l. l. 126. 13.8 ltOf'q-axis. 0,000 0,0000 0.0000 0.0000 iBLAltX CARD FOR THE 0-AIIS
5th S.ll, cerd. 0.002{ 0.07ú0 1.0600 i.002lf80 ,076 1.0É 0.8ó .25 0.8ó ,lg3 .183
tth S.ll. card, ó.7000 1.0000 0.0450 :á.7 1.0 .045 .ü55 ,08 ,000001 .000001
llass card I I t.000E+00 3.806E+01 0,000E+00 i I 1.0 38.056 .0001
Blank card teruinating S.ll. rass cards, iBLAllK CARD Ell0lll6 IíASS CARIIS
S.ñ, output request card fcr class 5.
S.l'1. output request cerd fsr class I.S,ll, output request card for class 2.
S,ll, output request card for class 3.
Blenk cerd terrinating S.ll. output requests, iBLAI{|( CARD E}lll¡llE 0UTPUT SPECIFICAIIÍII{S FtlR THIS S.ll,End of 0 S.lt.-TA[S interlece variebles, I FIlllSH
Blant card ends electric netrork sources, |BLAI{X CAR|} Et{DIl{G ALL S0URIE CAR0S
List of input elemnts that are connected to eech node. [}nly the physical connections ol rulti-phase lines ere shorn {capacitiveand inductive coupling are ignoredl. Repeated entries indic¿te parallel connections. Sritches are included, although sources
{including rotating mchinery} are oritted -- except thet U.fl, usage produces extre, internally-defined nodes 'UllIIIX".-+-----------
Fror bus nere i ll¿res ol all adjacent busses.
IIACAI |TERRA T RJOAT
RJOA i IIACAII
IIACB1 iTERRA T RJOFI
RJOB : IIACBTI
T,IAIEI :TERRA I RJOCI
5til2t3t
RJOC : ITAICII
TERRft :IERFA ITERRA ITEFRA TTERRA ITERRA ITERFA I IIAIATT I4AC[lI IIACIIT
-+-----------
Bus |t
ilACAt
Sinusoid¡l Eteady-stete phasor solutionr branch by branch,
is printed ebove the ilaginary part, the angle, or n0'.All llors ere ar¡y iror a bus,
Ihe first solutisn lrequency :Ph¿sor brench current
Rectangu I ar Pol ar
-2.3r{855508ólt 3.óf43t707f1521
2,775{435012529 130,f331220
.38191ó7203E-13 .39{771329E-13-,9992007228-r{ -rf.ü615{77
3.90373ó2541f92 3,9t86131972013
.3f113f3537ó95r r.9941982
.799360$i77E-14 .17256{0276E-13
.1529332216E-13 62.{0f5ó01
-1.8{779ó349?3? 3.6460397ó7792{
-3.r{3r28rÉt575 -12ü,{50É107
I-ragn Degrees
0pen 0pen
{}pen 0pen
0pen 0pen
InjectBd sBurce currentRectenoular Polar
3,90373ó2541f92 3.91S61319720tf
.34il343537ó95 f .99il982
-1.84ii963{9232 t.6f60397677925-3.1f3t281É1573 -120,r50úr07
end the real part, ragnitude. or oP'
6.000000008+01 Hertz.Parer llor Porer lossPend0 P¿nd0
-1612,30457t499 -1612.30{371499
-16951.79931259 -tú951.7993t2ó
.1200f783363E-9
.14128155093E-9
1.813047403f?97 1.95304740336tf
-18299.92353711 -18299.9235371
-,ó824{11úóE-t0-. t36214162E-t0
It17.35577üá372 t617.3i5776ó372-rü950.01722t08 -lÉ950,0172251
.7358281ó7{E-ll-.2381fiÉ5ó5E-9
Porer Reactive0pen 0pen
Open 0pen
0pen 0pen
Injpcted source porPr
P end 0 ll|lA and P.F.
-tat2,30f57t{99 17tr28.30073629
-1&951.79931259 -0,094ó838
1.8330t74035191 18299,92363093
-tB2?9.923537rI 0,000t0t3
1617.35577ó6372 170?7.00571559t
-16950.017?2508 0.09{9877
Bus llPh¡ssr node voltage
Rectangular Polar
7ú50.880093ó592 93f0.
5357.2039155188 35.0000000
7695.513839{5i5 9392.É152897394
5383.191613023ú 3{,983684?
814.03f637263t3 9310.
-9304.458{90177 -S5.0000000
fiJ(|B Bl6,ó99598i700{ 9387.16ó905086
-93tl.572?88ó2r -8f.0088589
-g{64,?t4i50922 93{0.39f7.2545áfó58I 155,0000000
HJ('A
ftAcct
HACFI
RJ0C -8t09.419070934 9390.59552{2733 -.2287059{38-13 .5071338ó75E-13
39i1.53?623022 154.9805057 -.4529709?4E-13 -11ó,7892997
Total netrorl loss P-loss by suming injections = 6.90f252if16{9E+00
0utput for steady-stete phasor sritch currents.ilode-X Nodp-ll
ilACAr
HACEI
ltAcct
l'lsde
ntteSource node voltage
Bectangular Pol¿r
I -rea I0pen
0pen
0pen
I -irag0pen
0pen
0pen
Solution et nsdes rith knorn voltage. ilodes that arp shortpd together by sritches are shorn es a group ol nales, rith the printed
result applying to the corposite group. The entry'llVA'is 50RT( Pllz + 0tl2 I in units of porer, rhile 'P.F.' is the
associated porer lactor.
ilACAI 7650.8800936592 9340. -2.3É4855508ótl 3,6fó3t707f1321
5357.2039155188 35.0000000 2,7754{35012527 130.{331220
llACBt 81{.r'3f637?AJl3-9304.458{80177
HACCT -8f64.914730922
3917.25{5ú{6581
9t{0.-8t.0000000
93{0,155.0000000
|IACH t Dat¡ parareters and initial conditions sl next r¡chine follar.' llACAl' unit Iflachine reartances and resistances. in IohlsJ. 0uentities l¡beled as inductances ere in f¡ct re¡ctancesl¡
É.04ú12050E+01 Lf = self-inductence of the d-¡riE field rinding;8.ó03{9127E+00 Laf = ¡utual inductance betreen field ¡nd arrature;4.97698943E+01 Lttd = field-derper rutuel inductence for the d-axis;l.É02114?lE+00 Ld = sell-inductence of the d-¡xis arlature rinding;g,ó03{9127E+00 Lald = rutu¡l inductance lror arrature to d-axis dalper¡6.38?486tBE+0t LXd = self-inductance of the d-axis darper rinding;4.59213955E+01 tg = self-inductence sf circuit I of the g-axis;ó,8i48l4tgE+00 Lag = rutu¡l inductance. q-aris circuit I to arrature;0.00000000E+00 Lgtq = rutual inductance, q-axis rircuit I to circuit 2;1.29982857E+00 Lq = self-inductance of the q-axis ¡r¡¡ture ninding;
0,000000üt)E+00 L¡kq = nutual inductance, q-axis circuit I to arrature¡3,34ó44732E+0i Ltq = self-inductence of circuit ? oi the q-axis¡
l,?0917?86E-01 Lo = zpro sequence inductance;
3.i$297714E-03 Ro = ¡ero sequence resistance;2.39370ú3EE-02 Rf = resirtance ol the d-axis field rinding;3,69997714E-03 R¡ = resistance ol the er¡aturei1.3472É082E+00 RId = resist¡ncp ol the d-¡xis darper rinding¡?.214732f5E+00 Hg = resistence of circuit I of the q-exis;0.00000000E+00 Rtq = resistance of circuit 2 of the q-axis;
llech¡nicel parareters in physiral units. as described by the colurn headings that lollor.llorent of inertie Self-darping coefficients of lass l'lutual-darping coeff. Torsionel spring constant
of rotor r¡ss Speed-deviation Absolute-speed {rith folloring rrssl trith folloring lass}
llillion fn-n/radlseclt2] ltillion [n-¡/rad/sec] ltillion ln-r/radlseclf2l ltillion [n-r/rad]1.6037169E+00 1,35ó3060E-10 O.Oggg6g6e+00 o.oqg6gggE+00 0.r1000000E+00
T¡bulation ol¡ negative sequence current. CltltlEG. CAilEE = 2.959J028314E-01 -2.0073387
zero spquencp current, CltlER, CAUER = l,033fl{B7llE-01 -3.0558f65
Tatal current injected into the netrorl at the generator bus, in phase coordinetes. For parallel rachines, this is the corbined
totel for all ol ther. The first line displays the currents as found by the phasor netrorl solution, and they rill be unbalancpd
il the netrort is, Ihe 2nd line shors only the positive-sequencp corpon€nts of these. llagnitudes ol currents are in units of
le¡psl.Phese 'a' injection Phase "b' injection Phase 'cn injection
fhqnitude Degrees |lagnitude 0egrees llagnitude 0egrees
Actual: 3.6ltit71E+00 130.133t?20 3.918ót52E+00 4.9941982 3,ó{Á0399E+00 -120.4506107
Pos. seq,. 3.72é03438+00 124,992q220 3,72ó0343E+00 4.99?!,ZZA t.t26¡343e+00 -11t.0075780
Arrature currents of generator in rotating refprence lr¡re {d-q-o coordinates}, in units ol Iarperesl,Id Iq Io
4,92038730353E+00 ú.34730901ú17E-02 -l.783525069498-01
Positive-sequence corponent ol generator ¡rr¡ture current in phase coordin¡tes, in units of Iarperes].Ar¡ature of phase 'a' Arrature of phase 'b' Arrature ol phase "c'
llagnitude Degrees llrgnitude Degrees llaqnitude Degrees
3.726034t8+00 t24.9i2422i 3.72É03{3E+00 4.9924220 3.72ó03f38+00 -115.00i3780
Field current of the generator in units of l¿rperesl. Tot¿l dc corponent
t.3287089972E+03 1.32874ú3{33E+03
Electromch¡nical torque of the generator in units of [rillion n-rl, totel dc colponentr . 5rü6992f59E-05 1 ,481f788391E-07
Critic¡l level of total air-gap llllF to begin saturation, in units of larpsl : 0.00000000000E+00 0.00000000000E+00
All rotor r¡sses heve angular velocity 0llEEA = f.7I23B89BE+01 [redians/secl. ]lechanicel angles of rotor r¡sses follor, inunits of illegreesj ¡ 15.ó2$0111
++++ Type-f9 S.ll, nurber I begins operation on segrents 0 ¡nd 0.
Card of neres for tire-step loop output. i llACAl tlACBl llACCl
Blanl card ending requests lor output variebles. llLAl{|( CARI} E}'lDIllE l{00E I/0LIA6E Í}UTPüI RE0UESI
Colu¡n headings for the 22 EIITP output variables follox, These are divided arong the 5 possible cl¡sses a¡ follors,...First 3 output variebles are electric-netrorl voltage diflerpnces (upper voltage rinus lorer voltagef;Hext 3 output veri¡bles are brench current¡ {floring fror the upper node to the lorer node};l{ext 16 output variables pertain to Type-f9 S.ll. corponents (n¡res are generated internally}¡Step Ti¡e lfACAl llACBl llACCt llACAt ltACBl ltACCl IIACH t ltA[H I
TEftRA TERRA TERRA ID IB
ilACHr |ACHI lrACHl lrACHr ltACHt frACHl ftACHl I{ACHI
IO IF IXD IÉ ¡|(O ]A IB IC
t{A[H IEFD
0.0 7650,88ü09{- . r 783i2507
-3t .80ó1423
.2E-3 7223.f$36S4
-.17ó80É09
-31,8061f23
,{E-I á759.?18ú9
-.174062153
-31 .80ó1{23
,6E-3 &25{.4?7t08-.t7012i997-31, 806l{23
.BE-r 5714.r22709-.1ó5857268
-31.8061{?3
,IE-? tt4l.ff324-. ló0295528
-st.80ótf?3.00? 1909,98082{
-. ll9f5B3ú7-31.80É1f23
.003 -1589,87263
-,0t218793-31,8061423
.004 -{8ó6.30195
,004015f939-31 .80ótf23
.008 -8261.89197
.1789356f09-3t.8061423
.otr 3828,7{9515
.0184482215
-31. B0Él{23
.016 8742.hi9724-.17tÉ922ó
-31.806r{23
,0? -7730,79412
-.0f050É839
-31 .80ú1{2J
,0{ -9358,55728
.1533{563tÉ-31 . B0ó1f23
üACH IIIFf}RCE
81{,03{ó373
1321.70899i
1329. 535947
1512. ú11679
l32B,70992r
1329.535772
2202,60b641
1328. il 176
t329,53s356
?979.998t99
1328.7t3889
t 329 . 53f9
3541. t7719{1328.717145
1329.531201
fi82,0ótt3t132S.7206f7
t329. 533{5
6962,742937
r3?8.7f29t91329.5?8üú7
87ó5.5Ár043
132S.7ó3959
1329.524il59337.316782
1328,772052
t 329 , 522{01
358. f7l0837r328. i08t9l1329.53ó054
-92e2.23984
t328.7710t71329, 522ó26
-l 525. fú285
t328.710178
t32q. 535t9t9100. ó39779
1328. ió8355
r329, 523r99
f810.47982{t328,72578{
r329.532335
ilA[H IIrAllE
-84ó4. 91471
-. (r2361f459
.1É179{8E-É
-s738.07161
-.022824232
,929t 37óE-6
-89Ét .93207
-.021 194037
,1723{238-5-9134. f3594-.0t9431095
.2399524E-5
-923Í,31304
-.01ó82É339
.3129{{iE-5-9323, 52001
-.01{053i22.3659303E-5
-8872.ó50I3
,00308518ó{
.52f0838E-5-7175.72085
.018ó829776
.3998285E-5
-f{71.05379
.02{0501i45
.57t921{E-É
7903, {20i76-. 02{243769
-.1 r29B5E-i
5{ó3,528011
.023578?706
,1s727078-5-7217,23007
-.023209767
-,23{753E-5
-É369.88051
.02183168ó
.29915348-f
4528.05772S
-.0u9222{7-.48t778E-3
0,04230,7507{9-2.10733254
0.0431.19061ú9
-,8{5872ú0f
0.0-1 522.09905
-, t065t005t0.0
-232.9t8794
.68310972ó3
0.0ró27.39{5
1.5134il311
0.0-.01 841 5031
0.0-.0i36ó{896
0,0 Ú.0
0.0 -2,5ó9{5305
0.0 0.0
0.0 -3,62738385
0.0 -3,37221101
4 .92039730{
3.9ü3736?54
4.9100353ú{
3.869297239
4.8869S4f54
3,809i6531
{.861828103
3.73102900{
f.82388r33r3.628612716
f.70323tf523.505103323
{. 527586995
2,ó1f732É82
f.2B9fi79ú5l. 35fót2009
{.201 52848
-.0950507
4,927372b24
-3.91792r9
{,2115{i{53-.398309143
4.9097307{4
3.Só83915{5
4,2{02t9837.8Bl9t 51882
4.73t7f5e3{-3.35925047
' ¡fter' efter' efter
-72473.052{
59363.3267
-26007.0799
1630?.9A059
30.3160ó952
201. r{03$1ó
-20{51 .9307
I 8003.30445
-50f80.0573
{0278.05988
.0t34730902-1,s4779ó35
.115f973129-1 . ó0ó09132
. ló8232?9r8-1. 35ó216fs
,2u5190507-1.095{05ó2
.257ó02t097-.8332ú8t li
.288ó7il81f-, 559f{5óf2
.3ót03I85{.90s5700043
,237731587f
2.05É855321
-.0ró598222
3.020if8412
-.027105317
2.22852t 579
.0756ó232{5-2.742231t4
-,1r589261
-?,57338r8
. l5s996f{552.42ü135t17
-.318779995
3.343539059
5.020000008-02
t.02000000E-02
5.02000000E-02
-2738S,3S4
-10755.0472
-32752,35{4
17638.176
2875.774ó54
-212t,7102É
27873.8007{-2783ü.545i
t262É,2r0rr-31s47. {112
IIACH I IIACH I HACH ITO EEil AI{E 1 |,EL I
rJ, u 0.0 ü.0-.ü10702312 0.0 -?.36485551
.1540ú998-{ 15,62501107 0,0
,2803{2tE-{ t5.62f01107 -,2É9031E-8
0.0 0.0 0.0-,02ó205É58 0.0 -2.75f83337
.{0833928-{ 15.6250rr07 -,É96ól3E-g
0.0 0,0 0.0-.0325¿023 0,0 -2.93099
.513{036E-{ 15.62501107 -.126952E-7
U.U
-.0393083830.0 0.00.0 -3,082Ét781
.62525f88-{ 15,É2501107 -,t97769E-7
0.0-.0f3794971
0.0 0,00.0 -3,22329768
.i00ó6{2E-4 1$,ó230t107 -,2902b2E-7
t0
l5
20
{0
60
80
100
.97ó29Í78-{ 15.ó?50110É -.796323E-7
0,0-.03ff23668
0.0 0.00.0 -3.52017998
,5770t27E-4 15,6250il06 -.126837E-t
0.0.003928ó212
0,0 0,00.0 -2,919772b7
-.f02953E-5 15.62501105 -,rtflrzE-ó0.0
.0028001891
0.0 0.00.0 1,999325939
-,65i808E-i 15.É250il03 .89105ólE-9
0,0-.009954{5f
0.0 0.00.0 3.1i2{935{l
.1856385E-4 15.62501t02 -.1t27{9E-óv.u
.01ót 1óf9r50.0 ü.00.0 -1.6010t972
-.28128fE-{ r5.ó2i0lt -,2i8llIE-80.0
-. 022f88039
0.0 0.0
,3959066E-f 15.ú2501099 -,13ó8188-6
0.0 0.0 u.u.0{70165¡2 0.0 .?813138t03
-.7737358-1 15,ú250109 -,37tbzhi-7tlsse sritch ' llACAl' to I
Close sritch " llf,CBl' to I
Close sritch ' llACCl' to '-fB5Él.02fl 5937ó.85526 -toBt7.lll4{31 .08t755 0,0 -f8608.3518
-2,479763b 15.6lfti5f9 -,036t27fó3-31092.073f 1óÍ09.23017 l7Í81.86ó774822,tb?221 0.0 -sft{1.3527-ó.719{ó98 15.58t0{661 -.031209709
1970,34?3 2t2.579{853 -2183.92983
-É35,632312 0.0 t9?5.4i8325.5425533528 15,538907?2 -.032i1fsl398bi,969727 18015,ü5093 -2788É,6029-{182 .020t2 0.0 9827.105004
5.0639f3tó7 15.f8ó93204 -.0{6S87996-s38t.93ó98 f0295.81967 -31899.7556-1708.07384 0.0 -8430.75677
200
uttilfit
300
400
500
ó00
700
sec.
sec,
sec..06 0.0
-. r5i204303-31.80ó1f23
.08 0.0-,124792ü7
-3i,90ó1{2I0,r 0.0
-.0990É2Íó
-31.80ó1f23
,t2 0,0-.078ó37936
-3t ,80út{23.14 0,0
- ,062{244f1-31.8061423
0.08t53. $55823
28t .2700397
0.0{t57.086237893,6399801
0.02549.0ó{503
1038.088729
0.04391,425f32
803.294397
0.07120.5S07S1
3i8,728óó18 3.53477i93S 15,4253083É -.0619ó804
800 ,16 ü,0 0.0 ü.0 -l?{17.7ú03 3907É,09104 -úó59.24598 -4872€.0457 -15775,{74ó
-.{){9553829 i081.2f3717 135ó.496óó9 24ú4.392ft2 0.0 -32458.68{3 39t¡65.03378 -5ó0¡.{353{
-31,80ú1423 t93.8553274 -1.{9úilllE -1,7út203{3 15,356218J1 -.0ü7ó047á1
900 .i8 0.ü 0.r) 0.0 -23{s$.32 79ú2.028A2É 15522.385t8 -225üt.87É7 -18709,795?
-.03933ó868 f809.S37078 -4?1.óúif83 2735.221ó01 0.ü -235?4.908S 7950.4f72{t 15574.39345
-31.80ó1423 628.993{704 -.óSf0f660t -t.502ú096{ lf.2S?i5&56 -,064505t99
l00a 0,2 0.0 0,0 0.0 5917.É7?531 -74{1,92338 1523.f05158 -9208.ú5207 2780.755918
-.031226{3 356{.2282f -t230.78ú4 -535.?50ó21 0,0 5S77.5ó2151 -if50.622i3 1573,00ó491
-31.80ó1423 742.99037É5 -.059295113 .J02{7821ó 11.205{8559 -.0ú{515933
il00 ,22 0,0 0,0 0.r' 15033,35376 7724,1791bt -22739.7fr07 -22732,19b2 16831.33795
-.02{788194 4ó79.f8241 -2ó1.2ó5795 -2515.97{01 0.0 t{995,13f17 7713.673123 -22708.8i02
-31.80É1423 5Br'.5795221 .5705850i24 2,802913172 15,1213718 -,071?98341
1200 ,24 0.0 0,0 0.0 -2305,0ó09 30?3r.87809 -27817,É907 -40529.ú791 7018,493877
-.0r9677388 b241,042121 904.0593636 -93f,990274 0.0 -2422,40141 30220.7718f -27798.f005
-31.90ó1423 260.9úl5lS7 1.00941829É r.8r26t208 15.0389700É -.07884{993
1300 .26 0.0 0.0 0.0 -25583.6715 30f58.03038 -f875.lBt8f -38881.0085:9754.39905
-.015ó20323 6t28,30óó81 741.478ó1f5 1535.0i95{3 0,0 -25621.3057 304f8.73533 -f827.f5672-31,80út423 191,2ó0t0üi -.6809{2383 -1.11933875 l{.9502f058 -.081t3501S
1400 .28 0.ú 0,0 0.0 -19f19.9ú17 3080,809f1 16338,30937 -231ó1.33f9 -1095É.0587
-.0123997t t737.98tÉúó -300.729209 tÉ02,780899 0,0 -t9f5É,357 3070,168843 1638É.16ó6ú
-3t,80ó1423 f91.É2{5ú7É -.{913ó7568 -t,8{79t88 r1,85957É32 -.0i8951828
1500 0.3 0.0 0.0 0,0 if80.f352ú2 -t2733,33f9 5252.258612 -t5551.0928 2022.ó436ó5
-.0099{3173 {008,525953 -7fL0i50J6 -3ó9.35i385 0.0 71t4,bil97 -127{4.3{55 5299,7197{8
-31.806t423 577.4r9t22t -,0212032i2 .168557405f lr.7t7f253i -.078697013
lú00 .3? 0.0 0.0 0,0 17886.5055É 892.0353934 -18779.3375 -23831.7173 10119.81659
-.0079t371i 4680,047832 -tf0.27ó829 -150É.3087É 0.0 17850.52232 882.3ú27293 -18732.8986
-31.80ól{23 fó2,521998 .{3tü7898ó2 l.ú1907724ó 1f.673297t9 -.082376t69
1700 .34 0.0 0.0 0,0 13ó2,15379 23702.5957 -25065,5795 -34247,97V3 4012,924ú3
-.006202692 55S3,721279 5{9.2{{883ó -527,{0153ó 0.0 1326,97581f 23É92,064 -25019,050ó
-¡1 . s061423 250, 571 537 . t571499502 .98f2671 704 14.57ó94068 -.08É121828
1800 .36 0.0 0.0 0.0 -2r38f.012 25408.{?4 -f025.31I23 -3?93t.8231 -5928.3634{
-.004923827 5{77.2É90ó5 437,t50ó5t4 935,3253171 0.0 -21{18,821S 25397.86i35 -3979,05f06
-31.80él{23 ?33.30031{8 -.3105337{8 -,óó9t021{9 tf.{7882989 -,087É52óó
1900 .38 0.0 0,0 0.0 -1753f.9522 8Ió.31459{8 r671i,82339 -23f28.i025 -6f38,835ó5
-,003908ó37 fó26.953009 -195.026{76 9{{.83i9999 0.0 -173ó9.8178 806.40?É018 ló7ó3.f3839-31.806r{23 fl4.t98{395 -,329508561 -1.03399909 r{.379ó7901 -.086104543
2000 0.{ 0,0 0,0 0.0 7f0{.33t371 -I5501,95s8 s09É.80{12 -189f7,5871 13ó4.1{73ó{
-.003102758 {188.70749ó -437.868t47 -239.549783 0.0 i369,9¡67 -15512.4189 81{2.45É97{
-31.80ó1423 {73.7383ú19 -.0062832óó .1076S20884 I+,279b6b17 -,0863il935?100 ,42 {'.0 0.0 0,0 19008.220s9 -3595.t706{ -t5ft3,ó731 -23920,9f91 6081.7S084
-.002fó3034 {583.r12316 -71,2262192 -900.t70367 0.0 l897l.0B2lS -3605,882?5 -1536S.2055
-31.80ó1f23 395.f897521 .308802759 .9432977268 lf.l7S4720t -,088f97393
2200 ,44 0.0 0.0 0.0 38{9.99t818 t9l2t.9t2l{ -2297?,7h8 -30000.561 2339.0t2608
-.ü01955208 5102,421ú7 336.4388419 -302.0tf905 0.0 3815,806079 19ilr.89899 -22927,?t)Bl
-31.80t1{?3 ?6$,6U8212,31r829{{98 .55{298895{ 1r.07591388 -.0909512ú5
23t'0 .46 0.0 0.0 0.0 -18fó6.8713 22339.108ft -3873.03ó8 -29030,82ól -35t5.370{9-,001552085 5014.198214 ?ú3.191ó{51 565.9il082É 0,0 -18500.7217 22328.695t -3827,97ú1
-31.80Á1423 2ó2.3588ó78 -.1537172t1 -.381ó7689t 13.97224138 -.ü9167512É
?t00 .{B 0.0 0.0 0,0 -túÉ92.2114 ilI.72ú68tf 16569.6ó528 -23262,1778 -378ó,13813
-.001232078 {f90.202149 -1r)9.0987 559.0335317 0.0 -16715,912 t01.285f343 1ó614.62439
-31.80á1f23 3ó7.77ú7S03 -.208303791 -,57174{50ó 13.8ó783t9{ -.09110142É
2500 0,5 0.0 0.0 0,0 ó493,045{29 -ló665.9055 t0l7?.0462i -?057ó.6313 S93,3397t22
-,?780498-3 {219.0á2É8 -2Í5.ó7914r -150.7f1991 0,0 ú{59.318325 -lÉ67ó.0És 1021ó.7f4i{-3t,80ó1423 t0É.{0713t6 .0072253852,0773t3997r t3,7b217472 -.09120ó822
2600 .52 0.r) 0.0 0.0 t9l52.3sif5 -659É.70786 -12556.{973 -2350t.5ó81 3óü1,4371ó8
-,776396E-3 {f{3,1ú5úfó -3{.22i9501 -537,S9{7S2 0.ü 19118,87522 -6607.090{8 -l25ll.7Sól-31.8061423 354.8051ó15,21S0775É26 .55ó6229527 13.65ó8618ó -,092ó21895
2700 .54 0.0 0,0 0.0 5ú26.390731 1t727.7365ó -21355.14f9 -27028,3f78 1387.3{7228
-,É163198-3 173ó.886721 207.349t952 -r74,97ó0i8 0.0 5593.t72729 15717.339{f -?1310.5132
-3t,soórf23 27b,s352299 .t877378885 .32t8477308 13.5t998f1ó -.094t8ó7082800 .5ü 0,0 0.0 0.0 -16?2t.7962 204t?.39423 -4187.41209 -2É329,8742 -2124,37586
-.{892f7E-3 {É6ó.255ó65 ló0,838270ó 3{0.92037f2 0.0 -t6219,229b 20f02.17399 -414?,94323
.58
0.É
,6?
.ú{
,bb
.ó8
0.7
7'1
,74
.76
,td
0.8
,tsI
,84
.86
.sE
0.9
3100
3700
4100
4200
43tttl
4400
-3r,80ú1423 ?7ó.36óÉ{ÍÉ0.0 0.0
-.388375E-3 4338,t359ót-31.804t423 337.28t0015
ü,0 0.0-,5083E-3 4t¿f.7lB0ól
-31.80úl{23 3ó1,02ó0194
0.0 0,0-.?4{735E-3 4286.22qi72
-31.80ó1423 327.ú304393
0.0 0.0-.19{?76E-3 f4f6.8l9SB-31.8061423 280.ó5773É
0.0 0.ü-,${228-3 {389.ó36f19
-31.80ilf23 280,2121543
0.0 0.0-.122fztE-3 4190,342931
-Jl.s06lf23 315.t502128
0.0 0,0-,971821E-4 1064,?33{99
-31.s06r423 329.04ó2ó19
0.0 0.0-.7714t2E-{ {t2{,193t66-3r.8061{23 307.5{99385
0.0 0.0-.É12395E-4 {20ó.828945-31,8061423 278.8593úil
0.0 0.0-,486132E-4 4159.4ó3729
-31.8061{23 27i.9ll?t770.0 0.0
-.3859028-{ f021.8S77{7-31.80t1{23 297.6844973
0.0 0.0-.306337E-f 3?40.37i9f1-31.80É1123 305.359t56
0.0 0.0-,2{3t77E-4 3963,9S8609-31.80ót{2t ?91,f20r933
0.0 0.0-.19303tE-4 f001,09379t-31,80ór{23 273,53{2089
0.0 0.0-.1532388-4 3960,713257
-3t,8061{2¡ 272,r84tó960.0 0,0
-.121ó44E-4 38ú6.8?0322
-31.80úl{23 283,0J9{998
0.0 0.0-.?ó56328-5 3006.48382ó
-It.8061f23 286,9013176
-.ü75924S97 -.20i6$tI250,0 -163t2.1308
-6t.ó5060i9 sil.ü863432-,124134459 -.309É39S6
0.0 5190.903108
-t{7,fl7ü7{ -93.049ó91
.015725SS2r .060115105f
0.0 18750.9t45i-l{.2521ó93 -Jzr.153311
.r{95201s56 .3343487fú8
0.0 6988.17f781
12?.lB6B72t -102,294284
.l2lé38850É .197238ó802
0.0 -1f340.0$07
r00.3059t2t 204.7380093
-,0332{1f{9 -.t034231
0.0 -16139.631-32,613t289 197,2741b24
-.0ó829179{ -,16131797{
0.0 3733,8{4213-83.t98ü69¿ -56.ó0óó078
.020955732,0492h894770.0 t80{0.5572
-3.{il27458 -t91.6É8798
.105879255 .20ó0{267É2
0.0 8090.170087
8t,8949308{ -60.29t25ú?
.08{5{31332 .126ttf33930.0 -12640.7969
É4,2$577879 122.ó523112
-.008{7911t -.0{ót907120,0 -16015,9101
-15.1{532ó5 117.5999223
-.03?t98ó99 -.07827ú68
0.0 22{5.5ü8902-45,200801 -3f.01?i9ff
.0211397987,04t89f11750.0 17150.0t?i9
?.f{820679 -lt{.3{213.077t2i8318,131528995J
0.0 9012.03977{
53,247¡7856 -35.78218t4
,0b2i37977t .085230989
0.0 -ll0f2.t74l42.ó12ó4{i9 73.33218r54
,0063268153 -,013672825
0.0 -158ü4.97il-4,82177349 70.2ó185ü05-,009216805 -.0318fre83
0,0 790,3912849
-22,809i229 -20.203876ú
.0260831224 .03É5988519
15.{42278É9 -,09f75812ú302.t¡óS7575 ló009.2f3{
0,0 -ló3t5.427513.333955ót -.09f608ó77
-tú947.9795 llú56.257850.0 -5157.óó6928
Ir.22505t5 -,09t8522-Bó53,?716ó -t0098,{873
0.0 1871i.69695
l3,tI5f598t -,09584f0ó1
t3058,007f9 -200fó.99s3
0.0 ó955.013309
13,00510ú02 -.096919252
t9t29.Á78ó2 -4790.{f394
0.0 -1f373,17ó
12.89f072t5 -,097422208
987.4f{11il 15151.3i225
0.ü -t6t72,7b?b12,78217798 -.097518Étt-ló{33.4tf9 12t98,i5501
0.0 3700.7tt838l?,6703{62t -,097839569-100?8,1707 -7943,2fr2t7
0.0 18007.5ü15
t2,537 hL22 -.09858ó743
t0834.99583 -18925.9707
0.0 8057.1099{9
L?.,44424029 -,099377711
r8192,25324 -5552.27179
ü,0 -12673.8233
12.5t02779 -,099838955
t92b,47 6b6 lt08B,ó179{0.0 -16048,9236
12,2t5i95r9 -,1000á7741
-r5É55.9961 13109.t72?8
0.0 22t2,4i2n7r2.10080813 -,10042274
-1il33,87{5 -ü016.9f3áf0.0 17117.0rú9r
lt.9852Bl36 -,10101?{9?
8S93.7382r -1790ú,5933
ü,0 8979,051285
11.969t9732 -.10tóts6ú1l7{23.5549ó -A3Bl.ó95ó7
0.0 -il07t,6653ll.il2$B58t -,10206662i
2973.4ú819ó t2890.597320.0 -15897.8{9f
ll.ó354897 -.102359783-11ó57.8118 13866,60519
0.0 757,1t7t9b21r.5t792023 -.t02722584
-22Sof.4f97 -222!,9ó235
291.70É2058 tó053,72065
-21151.2553 579.{9824&8
-tt959,3492 1t700.ó817ó
-22838.59óf 22I2.t7t8ó2-86ú3.53037 -10054,1ÉÉ9
-2{859.9383 839.0t40928
13047.ó5804 -20002.67t9
-2{3f3.54i2 -1251.9957
19il9.32591 -f7{É,1501ó
-2217J,t278 -1297.3{54
977.161041ú 15t95.ó0139
-2il26.6188 377,5215999
-16il3.7t58 127{2.9889r
-220t8.932É 1344.76074r
-10108.5129 -7898.98879
-23202.2707 519.715i647
10824.ó8735 -18881,7974
-22825.8{lú -726.63t233
tBlBl.9179 -5508.09473
-2rf57.23il -749,878723
1916.tf1032 14132,78247
-207b7.01h6 250.27360Só
-15666.30tf 13453.9U2
-212ó4.f905 825.5382ó11
-r11f4.2062 -3972.8r07
-2187{,5f0{ 531.9078072
8883.40óú32 -t7862.{58
-21584,f078 -4il.712{89t7{13.2{0Íó -6337.5752{
-20703,8058 -425.39{702
29É3.r39138 12934,71019
-20228.8óf5 171.2031f02-lf6É8.1409 13910.7?3úf
Extre¡a of output vari¡bles follor. 0rder and colurn positioning are the sare as for the preceding tire-step loop output,V¡riable raxina ¡ 9339.ó93324 9310,040987 9339.580841 19188.72739 ó4513.2{615 19814.35ó9{ 1383,4É8ó59 3Í5ó7.07196
,1790576461 8581.ó354{2 {807.21ú652 É173.{67753 0.0 19171.300{6 ó{519,0t736 t9817,62159-31.8061423 1330,010ó7 3,05278ró53 9.985r8502t 15,ó250t107 .t07gt99E-g
Tines of raxira : .0f84 .0206 ,0428 . {532 .0586 .0ó34 .067 .070ó
.0586 .0636.02f8 .058Á .0582 ,0b22 0.0 .4ó98
0,0 .0506 .0i5ó .055 0.0 .0198
Veriable rinire ; -9339.72309 -9340,06?3 -9339,52709 -52t05.76f7 -1719ó.060{ -{f03ó,9482 -78ó52.593 -{05É1,3123-.r79t3143ú 1328.i0S229 -158{.37619 -5241,5ilf9 0.0 -f2085.9538 -17398.0594 -f{003.2291
-31,80á1f23 11,78601ü31 -3.12311558 -8,31t79597 11,517920?3 -.tí?7722bTioes of niniaa : ,02t4 .045ó .0178
.01óf .0t64 .11óó
ü,0 ,1922 .t422
.üó1 .6674 .056
.r'70á 0.0 ,0ú12
.0t28 0.9 .893
,058É
,ó3{.0ó22
.05$B
Plot c¡rd ignored in search for blanl terrinator. i PRII{IER PLOT
Co¡¡ent card, flUltDCD = 5ü2{. iC[o¡¡ent cerd, lltlltl}0D = 5ó25, lC Titulo de la grafica[o¡rent card, flulll]C0 = 5É26, l[ 3-78[ompnt card. tlUltDC0 = fÉ27, itPlot card ignored in se¡rch fsr blank ter¡inator. i 2 lRAl{9lToRI0S DE LA llAoUIllA SII{CR0ilICA TRIFASICAEo¡¡ent rard. llUllDCD = 5ó29. :C
PIst c¡rd ignored in search for blanl teminator. i 144 2. 0,0 20 iACAI llAtBl iACilBlanl rard terrinating all plot cards. lFLAttt( CARII Et{DIt¡6 pL0T CARI}S
lleeory storage figures for the preceding, nor-corpleted data case,A value of '-9999' indic¿tes that no figure is available.
Size List 1. l'lurber of electric netrort nodes.Si¡e List 2. llulber of electric nettort branches.Si¿e [ist 3. llurber of d¡ta values in R. L, C tables.Size List l, Nu¡ber of electric netrork sources.Size List 5. Storage for [Y] and triangulrrized [Y]. llo, tires = 2 F¿ctors =
Size List 6. llunber sf entries in sritch t¡blp, llo. flops = 3
Size List 7. llurber of distinct ALPHAI{UllEftlC data nares plus prograr SPY v¡riebles.Size Lisi 8. History points of distributed lines,Si¿e List f. ilurber of nonlinear elerents,Eize List 10. Points of nonlinear cheractpristics,Si¿e List ll. l{urber of Type-S9 S.ll. outputs.Si¿e List 12. Total nurber ol EñIP output variables.9ize List 13. foriing space for batch/SPY plotting.Size List 14. S.ll./U.ll. connections to TACS,
Size List 15. ller TAIS char¡cter storage in brtes.Si¡e List 16. Total nu¡ber oi Type-5? S,ll, rasses,Size List 17. llurber of Type-59 Synchronous rechines,Sire List 18. Sr¿nch and sritch pomrienergy outputs.5i¡e List 19. Tot¡l floating-point TACS table space.5i¡e List 20, llon-copied recursive ronyolution d¡ta.Size List ?1, Total rodal/phese IT] ratrir storage.5i¡e List 22. Total recursive convolution history.5i¡e Lisi ?3. 6iant vectors lor renurbering, phasors.Size list 24. Peak phases of rorpensation for dat¡.Size List 25. Iotal table space for all tl.ll usage,Size List !6, Square ol lax nulber of roupled phases.
Tiling fioures characterizing central processor tCP) solution speed.Data input ti¡e {thrsugh the top of 'SUBRó,} ...,llode renuabering and phasor solutionAfter phasor solution, but belore tire-step loopIntegration of equations (tine-step looplPlotting or STATISIICS tertin¡tion overleys ,.,.
Presen tf igure
a
6
t23
17
3
I0
0
0
tó22
-9999
-9999
-9999
0
0
0
23
0
0
0
l50
-9999
-9999
CP sec
26,420
r3.290ú,530
92,220
2. f80
Prograrlirit (nare)
214 (LBUS)
T3I {LBRIICH}
424 (LDATA}
42 {LETCT)
2124 (LYüAI}
33 {LSTTCH)
59{ (LstrE7}
1487 (LPAsT)
ú3 iLilrlr{L }
I35 (LEHAR}
42 (LSfi0UTl
50 tLsrz12lr (L5rz13)
1ú {LESTAC)
r (LCTACS)
34 fLIñASS}
5 ILSYil)
5 (ilAXPE)
lú00 {LTAcsT}
aÉ0 (tFsEill
100 {rFD}I50 {LHISTI
3772 {LSr¡2313 {LCoilPI
{00 {LsPctrrlq2 (LS¡2261
I/0 sec 9ur sec
0.000 26.420
0.000 13.290
0.000 ó.530
0.000 92,?200.0ü0 2,480
Tot¡ls 140.940 0.000 lf0,9f0
Alternetive Transients Progran (ATP), IBll PC IT translation. Copyright 1987. üse licensed only through LEC i|(,U. Leuven, Belgiurl.D¿te {dd-rth-yy) end tire of day {hh.rn,ss} = 03-tlar-93 13:31:10 }lare of dist plot tile. if any, is C:330313J{,pI4For inforration, consult the copyrighted ATP EIIIP Rule Book published by LEC in July, 1987, Last ¡ejor prograr update: ltay, lgggTotrl length of 'LABC0ll' tables = 48197 Ii{TEEER rords. '|'ARD¡ll' List Siues lollor : 211 ZS4 124 42 ZlZ433 594 1487 ó3 t35 42 50 I la r 3f 3 5 t600 6ó0 100 150 3399 3 400 42 42
;;;.;;;;;;;-;;;;;;;;;;;;;;-;;-il;-;;;;-;;;;;.--l--;;;;il;;;;;;;;;-;;; ;;;;;-;;;;;;-;;; ;;;;;;;;;;;-;;;;;;;;;;0l23ffblB01234567890t231567890123{ 5É7S90123f5ú7890123{567890123f5É789012t{567890r 2345t7890
'----+-------
llarker card preceding nex E|IIP data c¿se. |BEGI}.| HEIJ DATA CA5E
[o;tent card, l'lUll[[O = 5ó]3, l[Blenl cerd ts terrinete EIIIP execution. iBLAt{l( CARD EHDI}IG PL0T CARDS
IkI
xfJ
10
U
- l(l
-tt
- jtt
- rtD
-qñ
-60{ lrlrl t$+ 3üS
1l}¡ itf,cH I - Ifi
50'i É0$
función del timpolss tenninales del
Comiente totai de la fnse a endurante una falla trifásica engenerador sincrónieo.
ú 100 .¿00 3ü0 r00
ññcH r - rE
üonniente total de la fase b endur.ante una falla tnifásica engenenadon sincnónico.
función de tiernpolos termi¡rales tlel
0Hlar93 13:31: Zs
BHlar-93 13:31:Zg
EHlan-93 13:31: Zg
l\ ililt11ililililililiill iilt f ilt iHHU+fHl']illilf il tf ttillilfiltf1ttlltiltitirl!itlillll
-.1[
-50
_ { 16;
c
ilfrcH 1 - IC t [rnsecJ
Conniente total de la fase cdnrante u¡ra falla tr.ifásieagenerador sincrónico.
en función del timrpoelr los tenritrales del
tml
50
U
-50
- 100
- 190
-2¡lCt
-
r 3]
150
100
s5*IIñLñ ¡ - ¡U
Conriente de
B}{lan9S
armaduna en eI eje 8.
BH1ar93 13:31: Zg
I
ritlt;tltlIl
iiIil[r
f..l
;'t:t.
It;
ti1i'lllir¡'l¡
li ii il ; I i\ ii i',. r I jlr¡l ¡',;í'r /:,;,i;,\ ¡,,¡,-*¡¿.i,,iri,'t¡"irTd¡'tr
r-1
-(?)nficHr-Iü
2ql'l ?00
----(?líÉcHt-Iü
Comiente de arnaduna en elConniente de arrraduna e¡r el
eje deje q
tkl
- 10
-
( lt
tc
--- ( tl nficÉr
fü
,J nfiCCI t [msec]n8c8 i
6Hlar--93 13:31:Ug
Tensionesenlasfasesa,b,cantesdepr.oducinseelfa|l@I n'votg,',t -',r r"- m0 dg OCtidcntc i
flHlar93 13:31: Zg
o
b
'l lDü 2¡10 3L1rl
10) ftftcH 1 - IF
Conriente de campo {IF}.
rlLJ
lfü0
7?rJU
100 ü
80rJ
bUU
T00
t00
U
- ?o0
- ( le)
BHlar-93 13:31:Zg
s 1ürl
llficH 1- iiFERcE
¿00
Fuerza magneturotriz total en el entrehierrode la máquina sincrónica trifásica
lI
tl1ttl
ÉHlar93 13:31: Zg
t [msecl
Falla tnifásica delconnierrte de falla
generador sincnónicoen el roton (IG).
AHEXO 3. EDIC]OIi PAfÁ SII{ULAR UI{A FALLft I'IOI{Í]FASICA
EEEI}¡ iIEH OATA IASE
PRIi'ITEI) I{UIBER HIDTH. 13, 2, { Hequest rexirur precision for I outpitt colunns
.0ü0200 ,900 ó0. ó0.
1llll-1 1
5 Í 2A 20 t('0 100
IC 5E[CIOilES PI DE UIIA LINEA
l lrA[A1 fiJt|A 9.04 44.88 351,28
? ilACBl fiJoB 7.84 19.04 -57,ú 9.2 44,5ó 3ó0.08
3 ttAcct RJ0[ 7,É8 14,9ú -14,9É7.84 19.04 -57.6 9.0{ 4{,88 351,28
;fiLAtll( cerd ending all br¡nch cerds
E
[ ¡ITTERRUPTOR DE SIIIULAEII]I{ I)E FALLA DE LA FASE A A TIERRA
ItACAt , 05 10. I
BLATX card errding all sritch cards
59 lrACAl 9t40.0 ú0.0 35.0ItACBl
ltACcl
TI]LERANCES 2i'
PARAIIETER FITTIIIE 2,0I I ló1. l. 12ü. r3.8 1ó04.
BLAI{lt card lor the 0-exis,00244S0 ,076 1.06 0.8ó ,2J 0.8É .18¡ ,lB3
6.7 l.{¡ .0{5 .055 ,0t .0('0001 .000001
t t,0 38,05ó ,0001
BLAI{K cerd ending rass cards
5ttl4ttl3t
BLAiII cerd ending output specificatians for this S.]1.
F IItI SH
BLAI{I( card ending all source cards
iA[Al ilACBI nACCt
BLAII|( card enriing node voltage output requestFfiIIITER PL[]T
II Titulo de la grafica
C
2 TRAI{SITt]RIt]S I}E LA ITAOUII{A SIilCRÍ}}IICA IRIFASICA
144 2. i'.0 20 ltACAl ltACBl llACEl
BLAII|( CARD EIII}IIIE PL[}T CARDS
BEGII¡ IIEII I)ATA cASE
BLAI{Ii
ANEXO 4. LISTADO Y ANALISIS DE LAS GRAFICAS QUE CORRESPONDEN
A UNA FALLA MONOFASICA.
la, lb, lc (Corrientes de la fase a, b, c respectivamente).
Después de que ocurre la falla, se encuentran corrientes alternas y continuas
en la fase a del Generador; en las fase b y c las corrientes son iguales a cero.
Estas corrientes continuas de la fase a aparecen en el Generador Sincrónico
debido a que éste es básicamente inductivo.
En 50ms tiempo en que ocurre la falla, la componente alterna de la corriente de
la fase a salta a un valor muy alto, pero la corriente total no cambia en ese
instante.
lf (Corriente de campo).
En 50mstiempo en que ocurre lafallaestantrabajando las barras amoftiguadoras
y porellas circula corriente continua, se inducen FEM de corriente alterna en ese
devanado. A partir de ese tiempo en adelante, trabaja el devanado de campo
por el cuál circula también corriente continua, ésta corriente de campo se
compone de una componente grande de corriente continua y de corriente
alterna.
UnivsrS¡uurJ'u'rrcfllo do
5ección Eiblioteco
lg (Coriente de falla en el rotor).
Señal de componente de cuadratura de falla en el rotor, aparece modulada en
amplitud por una señal de fracuencia menor.
lo, lo, lo (Corrientes de armadura en el eje d, q, 0).
to Corriente en eje directo, es una corriente negativa, aparece elfenómeno de
la modulación éste fenómeno desaparece en el estado transitorio.
loCorriente en eje en cuadratura, Es una corriente de doble frecuencia que la
corriente en eje homopolar (lo).
lo Corriente en eje homopolar no existe el fenómeno de la modulación.
IKD, IKQ (Corriente en el devanado amortiguador D y Q respectivamente)
La componente en eje directo es una componente reactiva (inductiva)
desmagnetiza la máquina, aparece una componente de corriente continua de
menorfrecuencia, se presenta elfenómeno de la modulación, elcuál desaparece
en el estado permanente.
MFORCE (Fuerza magnetomotríz total en el entrehierro).
En el estado subtransitorio existe una componente aperiódica que varía en
forma aperiódica, la existencia de picos grandes siendo pequeños los coefecientes
de atenuación indica que fuerzas perturbadoras de pequeña magnitud provoca
oscilaciones grandes en el rotor de la máquina sincrónica.
Esta fuerza magnetomotríz varía en forma exponencial.
Nota:
El fenómeno de la modulación se presenta por las siguientes causas:.
En fallas asimétricas, en el generador actúan sólo momentos de rotación
condicionados por las corrientes de secuencia positiva que circulan por el
estator del generador, las corrientes de secuencia cero no circulan por el
generador, en caso contrario no Se forma un campo giratorio capaz de
interacionar con el campo del rotor.
Las corrientes de secuencia negativa crean un campo giratorio que está dirigido
en sentido contrario alcampo giatorio creado porel rotor, la interacción de estos
campos provoca un momento pulsante y cambia su signo con una frecuencia
menor que la frecuencia normal; del generador.
Alternative Transienis Progran {ATP}, IBtt PC IT tr¡nslation. Capyright 1987, Use licensed only through LEÜ (K.U. Leuven, Belgiul)'
Dete {dd-nth-yyi and tine of day {hh.m.ss} = 06-Apr-93 20:37153 Hare ol disk plot file' if anl, is [:340t2037,p14
For infor¡etion, consult the copyrighted ATP EIIIP Rule Eooli published by LEü in Julyr 1987. Lest ¡ejor progrer update: tlayr 1988
Total lengih of 'LABC0ll" tables = 48797 II|IEEER rcrds. 'VAffIrlH' List Sizes lollor ¡ 214 ?34 f24 12 2124
_ __:1___:::__::1 _ i:---:::-- i1- --::-----:--.-1----1----::- --1-----:--::::---:::---y---T--::::---- :---::1----:1---::Descriptive interpretation of input data c¡rds. i Input d¡ta card ilages are shorn belor, all 80 colurns, character by character
0li345b7E0l 234567890123456i8901?3f567S90123fSói990r2345ó7890123{567890t23f5t789012345É7990
-----+-------Cor¡ent card, llUllI}CI} = 5582. !C data:liA0ACl.0AT
I'larker cerd preceding ner EIITP data case. iBEEIll t{Ef DATA CASE
Hidth ol tire-step loop nurbers. t=13 5=2 IPRIiITED HUIIBER lf¡0TH. 13, 2, { ffequest rarirur precision lor I output colurns
llisc. date. 2.000E-0{ 9.000E-úl 6.0008+0I i .000200 '900 t0. 60.
llisc,data. I I I I t-l 0 I 0 ü i I 1 I t t -l tlll [arning, llon:ero t0PT sr C0PT diflers fro¡ the porer lrequency ol 5.00000000E+01 Hz. l)id the user really sant this?
Frintout : 5 3 20 20 100 100 i 5 5 20 20 100 100
Cor¡ent card. llUfl|)CD = 5588. lC
Eoment card. tlUltDCD = 5589. iC SEC[I0]|ES Pl I}E UNA LINEA
lst ol Pl-ctt. t.0f0E+00 1.1908-01 9,318E-07 I I llA[Al RJ0A 9,01 {4.88 351.28
7.840E+00 5.031E-{)2-1.5?8E-07 9.2008+00 I.t82E-0tl 2 llACBl RJo[ 7.84 19,04 -57,6 9,2 {1.56 3t0.88
t.6gqt+00 3.968E-02-3.9ó8E-09 7.Bf0E+00 5.051E-02i 3 nAccl RJ0c 7,68 lf.96 -14.967.84 19,0{ -37.ó 9,0{ fl.88 351.28
[oment card. ]iUI'II}CD = 5593.
Cornent cerd. llUltDCD = 559{.Elank c¡rd ending branches. ISR, HTEI = 3 7 iELAI{l( CARD EtlDlll6 ALL DRAIICH CARDS
Coment card. flUllI)C0 = 559ó, iC
Coment card. l{UltDEI) = 5597, iC II{TERRUPTOR DE SlllULACIt]}¡ DE FATLA DE LA FASE A A TIERRA
$¡ritch, 5.00E-02 1.00E+01 0.00E+00 0,00E+00 i llACAt '05 10.
Conrent c¡rd. HUIIDCI) = 5599,rF
Blenk card ending sritches. (SITIH = l. iBLAil( CARI) El¡I}IllE ALL SIITCH CAR0S
Scurce. 9,34E+03 6,00E+01 3,50E+01 -9,91t+03 lf9 llA[Al 93{0.0 60,0 35.0
?nd phase of 5,11. 9.340008+03 -85.00 i llACBl
3rd phase of S.ll. 9.3{000E+01 -205,00 i ltA[Ct
[psilon. 1.00E+02 t,00E-tS 1.00E-lá 20 iT0LERAiICES 20
0ptirize Perk's dat¡. 1.000E+00 IPARAIIETER FITTIIiE 2.0
4th 5,11. card. I t 0 ló 1.000 t.000 i I t lÉ1, 1. 126. I3.8 lÉ04.q-axis. 0,000 0,0000 0,0000 0,0000 IBLAI{|( CARD Ff}R THE 0-AIIS
Sth S,ll, card. 0.0024 0.0760 1.0ú00 1.002{180 .076 1.06 0.86 .25 0,8ó .183 .183
óth 5.I1. card. 6,7000 t.0000 0.0130 ib,7 1.0 .045 ,055 .08 ,000001 .000001
ltass cerd I I 1.000E+00 3.806E+01 0,0008+00 : I t.0 38.056 .0001
Fl¡nk rard terrinating 5.11, mss cards. |BLAI{|( CARD EllDlllE IIASS CARDS
S.ll. output request card for cless 5. i 5I
S,ll. output request cerd for cless l, I llS,ll. output request card for cless 2. i 2t
S.ll. output request r¿rd for class 3. i 3lBlank c¡rd teninating 5.11, output requests, iBLAill( CAR|) EilDIllÉ 0UTPUT SPECIFICATI0I¡S FOR THIS S,ll,End ol 0 5.1'1,-TACS interlace variables. I FIt'llSH
Bl¡nl card ends electric netrork sources. lBLAilf CARD Ell0lt{6 ALL S0URCE CARI)S
List of input elelents that ¡re connected to each node. 0nly the physiral connections of rulti-phase lines erp shotn {rapacitive¿nd inductive coupling are ignoredl. Repeated entries indic¡te parallel connections. Sritches are included, although EourcPs
Y:1T:-:::::r:r-:::::r::il-:::-::::I1-:: **"ot thet u'n' usase produces extra' internarrv-derined nodes 'unlIII"
::::-:::-::::-;--11:::-:1-:11-:11:::i:-::::::'IIACA1 ITERRA I RJOAI
RJüA
fiAcBt
RJOB
flAcct
RJOC
icrr
ñACAI IRJ(¡BI
ilA[8lfRJOCI
ilA[CtlTEfiRA |TERRA TTEFRA TTENRA ITERFA ITERRA ITERRA T IIACAIT
-+-----------
Sinusoidel ste¡dy-siete phasor solution, branch by brench.
ir printed ¡bove the inaginary part, the angle, 0r u0',
Bus X Phassr node voltage
All llors are er¡y frol e bug.
The first sslution frequency =
Phasor brench current
end the real part, ragnituder Br "P'
É,00000000E+01 Hertz.Porer llor Porer lossP¡ndg P¡nd0Bus ll Rect¡nquler
llACAl 765tt.880093óf92
Rectangu I er Pol ar
-2,3t4855Íü8611 t.6ft31707f1521 -1612.30f571499 -1612,30f571499
2.7754{35012528 130.4331220 -1695L7993t259 -16951.799312ú5357.203915ÍlS8 35.0000000
RJoA 7ó95.5138394575 9392.ó152887391
5385.1916130236 34.903ó8{?
ilACBI 81{.ü3f6372ó313 9310.
-930{.4s8{80177 -85.0000000
RJoB 81ó,É99t9877004 9387,1úó90508é
-9351.572288ó21 -85.0089589
nAtcl -81ó4,914730922 93{0,
39{7.25f56{ó5Bl 155.0000000
fiJ0c -8509,11907093{ 9390,59552{2733
397r,532É23022 I5{.980Í057
Tot¿l netrorl loss P-loss by suming injections =
0utput for steady-state phasor sritch currents.
Pol ¡r
93f0.
l{ode-l( flode-lt
ilACAt
I -rea I
0pen
I -irag0pen
.3819ró7205E-t3,39f771329E-13-,999200722E-r4 -14.óó15177
3.90373É25f1f92 3,918ó131972013
,3{il3{3537ú951 1.9911982
,7993605777E-t{,172fó|0276E-13.1t2933221óE-t3 62.{045ó03
-1.8{77963f9232 3,Éf60397ó77924
-3.1f31?8tó1573 -¡20.f50ó107
-.228705943E-13 .50i{338ó75E-13
-.f5297099fE-13 -tt6.7892997
ó.90f2t?541ó{9E+00
I -ragn0pen
Degrees
0pen
.120047833É58-9
. lfl28t55093E-9
I .8530f7f03{297
-18299.92353711
-. É82{4rtúüE-t0-,43ó2t{tú28-10
1617.3S5776É372
-lú950.0172250t
,735S28167rE-u-.2381f16565E-9
Porer
0pen
I .833047403tál{-r8299 ,923537r
r6r7.35577ú6372
-tú950.0ri2251
Reartive0pen
Solution et nodes rith knotn voltage, Nodes thet ere shsrted together by sritches are shorn eE a groüp ol nales, ¡ith the printed
result applying to the corposite group, The entry'llVA'is S0RT( Pltz + 0tl2 l in units ol porerr rhile 'P.F.' is the
¡ssocieted porer fartor,l{ode Source node voltaqenele Rectanqu I ar Poler
IIACAT 7ó50.880093ú592 93{0, -2.36f855508Éil 3.óf63170741521
5357,203915518S 35.0000000 2.775{{33012527 130.1331220
ilACBI 8t4,03fó372ó313 95f0. 3,9037362541f92 3.9196lll9720lf-930{,f58f80177 -85,0000000 .r{11343537É95 {.9941982
ñACC! -Sfúf,914730922 9340, -1.8477963{9232 3,64ó0397677925
39f7.25{5Ét6581 153.0000000 -3.11Jt28161575 -l?0.f50ó107
IIACH I Data parareters end initiel conditions of next rachine lollor'' l{ACAl" unit Illachine react¿nres and resistanres, in [ohrs], 0uantities labeled es inductanres ere in fect react¡nres):
ó.04út?050E+0t Lt = self-induct¡nce of the d-axis field rinding;8.ú03491?7E+00 L¡f = rutual inductance betreen field end aruature;4.?76989{38+01 Lfkd = field-darper rutual inductance lor the d-aris¡l.ú02t1429E+00 Ld = self-indurtance ol the d-axis ar¡ature rinding;8,t03491?7E+00 Lakd = rutual inductance fror ¿rrature to d-¡xis darper¡É,382186t8E+01 L(d = sell-inductence of the d-¡xis darper rinding;4.59213955E+01 Lg = self-inductence of circuit I of the q-axis;ó.85481{lgE+00 Lag = rutuel inductance, q-axis circuit I to arnature;0.00000000E+00 Lgtq = rutuel inductance, q-rris circuit I to circuit 2;l,?9982857E+00 Lq = self-inductance of the q-axis arrature rinding¡0.000000008+00 Lakq = ¡utuel inductance, q-axis circuit 2 ta arrature;3.3464{732E+01 Lkq = self-inductance of circuit 2 ol the q-aris;1,?091728óE-01 Ls = zero spquence inductance;3,70297714E-03 Ro 3 :ero sequence resistencei
Injected sourte currentftectangular Polar
Injected sourcP porer
P end 0 llVA end P.F.
-16t2.3045i1f99 17028.300i3É29
-tó951,79931259 -0.09{ó838
1,8530t7f035I91 18299.923Á3093
-19299.92353711 0.0001013
ló17.355i766372 t7027.005715591
-tÉ950.01722508 0.09{9877
2,39I70ó38E-07 Rl = resistance of the d-exis field *indingi3.ó9997714E-03 R¿ = resiEtance of the arrature¡1,3f726082E+00 Rtd = resist¡nce of the d-¡xis d¡mper rindingi2.214732{5E+00 Rg = registance of circuit I sl the q-axis;ij.0000000üE+00 Rlq = resistance of circuit 2 of the q-arisl
llechanic¡l parareters in physical units, es described hy the colum headings that follor.llorent of inertia Sell-danping coefficients of rass llutual-darping coeff, Torsional spring constant
of rotor rass Speed-deviation Absolute-speed {rith lolloring rass) irith folloring rass}
llillion ln-r/redlsecll2] llillion ln-r/rad/secl llillion [n-r/rad/secll2J llillion [n-¡/rad]t,ó0371á9E+00 1,35630ú08-10 0,0000000E+('0 0.0000000E+00 0,0000000E+ü0
Tebulation ol¡ negative seguence current. üll{EÉ, CAIIEG = 2.959302831fE-01 -2.0073387
¡erg sequencp current, CllZER, DAZER = 1,0335148711E-0t -3.05585É5
Total current injected into the netrorl at the generator bus, in phase coordinates. For parallel rachines, this is the corbinedtotal lor ¡ll of ther, The first line displars the currents ¡s lound by the phrsor netrork solution, and they rill be unbalanced
if the netrort is. The Znd line shors only the positive-sequence corponents of these. llagnitudes of currents are in units of
Iarps].Ph¡se 'eo injection Ph¿se 'b' injection Ph¡Be 'c' injection
l'lrgnitude Degrees llagnitude llegrees llagnitude Degrees
Actual: 3.64ó3171E+00 130.1331220 3.918é132E+00 4,9?41982 t.6¡6¡399E+00 -120,450ó107
Pos. seq.: 3.7?É03{3E+00 t2{,992+220 3.72603{3E+00 1,9i24220 3.72603{3E+00 -1t5.0075780
Ar¡ature currents ol generator in rotating reference lrare {d-q-o coordinates}, in units ol l¡rperes],Id Iq Io
4.92038730355E+00 6.34730901ó17E-02 -1.783525069f9E-01
Positive-sequence corponent of generator ernature current in phase coordinates, in units of Iarperesl.Aru¡ture of phase 'a' Arrature oi phase 'b' Ar¡ature of phase 'c'
ñagnitude Degrees llagnitude llegrees ltagnitude Degrees
3.7260343E+00 12f.9924220 3.7?603138+00 f.992f220 3,72ó03ftE+00 -il5.0075780
Field current of the generator in units of larperes]. Total dc corponent
1.3287089972E+03 1.3287403433E+03
Electrorechenical torque ol the generator in units of [rillion n-rl, tot¡l dc corponent
I .5f0É992ft98-05 l. 481 f7S839lE-07
Critic¿l lpvel of total air-gap tlllF to begin saturation. in units of lerps] ¡ 0,00000000000E+00 0.0000ó000000E+00
All rotor ¡¡sses have angular velocity [}llEGA = f,71238898E+01 [redianslsecl, llechanical angles ol rotor regses lollor, inunits of Il}egrees] : 15.ú250111
++++ Type-39 S,ll, nurber I begins operation on segrcnts 0 and 0.Eard of nanes lor tire-step loop output. I llAtAl tlACBt llACCl
Elani c¡rd ending requests for output y¡ri¡bles. lBtAll|( CARD EtI}¡tlE ll0DE I/0LTAGE I}UIPUI RI0UEST
Colurn headings for the 20 EIITP output v¡riebles follor. These are divided arong the 5 possible classes as follo*s....First 3 output variables ere Electric-netrork voltage dilferences {upper voltage rinus lorer voltage};llext I output variables ere brench currents {floring frol the upper node to the lorer node};l{ext tú autput vrriebles pertain to Type-í9 S,ll. corponents {neres are generated internally};Step Tine llA[Al ltACBl ftACCt llA[At IIACH I IIACH t IIACH I tlA[H I
TERftA ¡I} IO I(' IF
ñACHI ltACHl ñACHI ilACHI lfACHl ilACHI ilACHI ñACHIII(D 16 ITO IA IB IC EFI} IIFORCF
ilACH I üACH I ilACH I IIACH IñAilE I0 6Ei¡ Ail6 ML t
0 0.0 7650.ss0094 8l{.ú3{6373 -Sfú4.91{73 0.0 {,920387304 ,0ú317i0902 -.178352507 13?8.709997-.11?3É14{39 -.010702¡12 0.0 -2,3ú481:5i 3.90373¿25{ -1,84779635 -5t.80óil23 132t.535947
.tá179{8E-ó .1540ó998-4 r5.325ut107 ü.üI .?E-3 7?25,451684 ltl:.ó11ó79 -873S.('71ói ú.0 4.91{,t)I5364 ,115f973129 -,17É80609 1328.709821
-,0228?4237 -,018415031 0.0 -2,569{5305 3.8ó9297238 -1.60608132 -il.8061423 1329.535i72
,929157áE-ó,280342tE-4 15.6250il07 -,2ó9031E-B
2 ,4E-5 6759.2t869 ?202.606ó{l -89ó1.83?07 0,0 4.88ú?84554 .1682322918 -.17{0ú2153 1328.71176
-.ü21194037 -.02ú20f658 0,0 -2,75483337 3.8095ó53t -t.35É216{3 -31.80úil23 t329.53Í356.r723f238-5 .4083392E-{ 15.ú250U07 -.ü9áúr3E-B
3 .óE-3 É2$4,427108 t879.998t99 -913f.43594 0.ü 4.8É1828103 ,2115190507 -.170329997 1328,713889
-.019f3t095 -.03?5t023 0.0 -2,9¡099 3,731029004 -1.095{0562 -31.806t{23 1329,5319
.239952fE-5 .513403ÉE-4 13.ó250t107 -,t26i12E-7{ .8E-3 571{.r22i09 35{1.17719{ -9253,3130r 0.0 4.8238S1331 .2576021097 -,1ós8572á8 1328.7171{5
-,0ró82É339 -.03930838s 0,ü -3.082Ét7Br 3,ó2861271É -.8332ó8t17 -31.80ilf23 1329.53{201
.3129f{tE-f .ó2525{88-4 15.6250tt07 -.t9t7biE-?5 .1E-? Ír41,4432f 4182,06tt31 -9323.52001 0.0 4,78323ú452.28867llBrf -,1É0295f28 1328.720ó{7
-.0tf053722 -.0{3794971 0.0 -3.223297úg 3.505103323 -.559{456t2 -31.80ú1423 1329.533f5
.3659303E-t,7006ó42E-4 t3.É2501107 -,2802t2E-7
t0 .002 1909.8808?4 É9É?.7{293i -8S72.ú50t3 0.0 4.527586995 ,3ú103385{ -.1195583ó7 132S.7{2919
.003085184{ -.0536óf896 0,0 -3.ó2738385 2.614i32ú82 .80557000t3 -3t.S061{23 132?.528óó7
.32|0B3SE-5,8762957E-4 1s,t230il0ó -.79ú3238-7
15 ,003 -1589.872ó3 87É,5.561043 -71i5.72085 0.0 f.289f779ó5 ,237731587r -,0ú219793 1328.763959
,0r8É82977t -,03{423ó68 0.0 -3,5?017998 r.355612009 2.05695532r -31,806r{23 1329.524145
.39982S5E-5,1770t27E-4 15.ú250110e -.r2ó83iE-ó20 .00{ -{8á6,30195 9337.316702 -4{n.053i9 0.0 {.20r52S{B -.01ü598222,004015{939 t3?8.772052
,02{0501t45 .0038286212 0,0 -2,918752b7 -.0950507 3.020758f12 -31.8061f23 1329.522401
.5749214E-ú -.402953E-5 15.62501t03 -.14{112E-ú
40 .ü08 -82É1.89197 358.f710837 7903.42077ó 0.i' 4,927512b21 -,0?7105317.1789336t09 13t8.70S491
-.02f2437ó9 ,0029001891 0.0 r,999325939 -3.9179219 2,22852ttr7i -3r.90ól{23 1329.336054
-.1r2985E-5 -,ó57808E-5 r3,ó2f01103 .89105É1E-9
60 .012 3828.7f9515 -9292.2398{ 54ó3,528011 0.0 {,211547f53 .07566232{5 ,0184{82215 1328.771017
.0235782i06 -,009954434 0.0 3.172491541 -,3983091f3 -2,7422311t -31.806t{23 t3?9.52262É
.t87270/E-5 .1836385E-{ t5.62t0u02 -.lf27r9E-680 ,0ró 8712,h99724 -1f25,46285 -7217.73007 0.ü f.9097307f4 -.ltÍ89261 -.176É9226 1328,710t78
-,023209767 ,01ó11ó1915 0.0 -Ló010{972 3.8ó8391545 -2.5733818 -31,80ilf23 t329.535691-.?3{i53E-3 -.281285E-4 15,6250u -.27glttE-8
100 ,02 -2730,791h2 9t00,ú39779 -t3t9,8805t ü,0 4,2{02f983i .t58996{{55 -.0f050ú839 1328,7t8355
.02185168ú -.022{88039 0,0 -r.37221t01 ,BBl9tÍ1882 2,{20135917 -3t,80ó1423 1329.523199
,2994334E-5,39590ó6E-{ 15.62301099 -.r36818E-É200 .04 -9Ir8.55728 fgt0.f79B24 4528.ü577?8 0,0 4,7347f593f -.318779?95.153345634ó 1328,725784
-.0ili?22+7 ,047016532 0,0 ,2813139403 -3,359250{7 3.3t3539059 -3t.80ót423 1329,532335-,484778E-3 -.773735E-4 t5.ó250109 -,374b26E-7
Itt Close sritch ' llACAl' to ! t ¡fter 5,02000000E-02 sec.300 ,06 0,0 -5397,57005 792b,h1617 -63220,2418 -48839.862 tá839,18792 -365t5.0t2f 57t7,09r{02
3110,78270ú -2221.93785 0.0 -63263.230f {.833877103 12.5397508 -3I.9061{23 É67.97{512r.95É2131451 6,0504572ó1 15,622835ál -.018112709
400 .08 0,0 -7129,29752 3137,7i2272 -{8{2t,7586 -23899,3679 -31537.{¿4i -27i17,0t09 3895.9rSÉ371004,091548 {823.8ff96{ 0,0 -48fú1.3705 -5,9980352 9.76f2flf8t -3l.B06tf?3 9I7.1005954-,576975?8{ -ó.60855205 15.É07621t3 -.019213618
500 0.1 0.0 ?444.133383 -ó017.2332{ 2t28,298723 -972,8t3393 1f00.083205 1211.015ú9 LiS2,b?,6992-Éf0.888?0ó -256.7938{ 0.0 2091.217303 r.3ú58ó2993 4.937538369 -31.80ó1f23 1173.ú2É835-.010025777 .2.s785241b4 I5,58850{ -.0t55t3SS2
t00 .12 0.0 75ó4.208248 -t232.269 167Jó,32ó51 -130É7.2f5ó -387r.122ó7 96f5.103159 3103.5704t981.3773{701 558,{?40382 0.0 16ó9ó,9970f f.838i3É40ó 3,972935At2 -31,8061f23 930.3080823-.095fr68{9 -.86ó5Jffál 15.5óf5fS95 -.0205218
700 ,14 0.0 2285.091812 1883.122f93 -16ü16.t552 399,0927339 13100.?Bt5f -92É5,16536 t930.5f8336-828,02f525 -1969.9É975 0,0 -1ó050,5{38 -1,047295á2 6.853935f13 -31.80ót{23 1194.798{30.1977500378 2.78{853598 15,53ü0942t -.0239918f
s00 .lt 0.0 -{958.ó73É ó815,ó+É{73 -f7088.4576 -3ór55,31¿l l3l5?,771?1 -2i205.19t7 $tóó.93á{461552,058311 -t809.19415 0.0 -47t20,3?f -ú,0{t9193 3.5829É0892 -31,S0Ál{23 597.41ü5281,úó8051709t 3.5ó895{8ót 15.50270003 -.0t583912
900 .18 0.0 -6ó94,62 310{.327023 -351É0.0595 -t77iq,902 -225t0,5154 -203t8.5885 3Éf5.522i992I9,0f02272 3427.00692 0,0 -3519f.s79 3,69{910325 -1.9r354099 -31,90ót{23 873.081902-.44ó2401{9 -4,fó85154 15.46489378 -.035749979
t00$
I lt'0
l?00
1300
t{00
$00
ló00
I 700
1800
190ü
2000
2r00
2200
2300
2r00
2500
2600
2700
2800
2900
3ú0rl
(,,1
.22
,21
.24
.28
0.¡
.32
,3f
.3ó
.38
0,4
,42
.{4
.4ó
.{8
0,5
.52
.54
.56
.58
0.ü 197i.7{5ó73 -5{20.?9548 ?3ú8,9721ú5 -4191.15323 úJt{.3872s8 539I.1795u9 2f1t.40579s-64u.i8Br'48 -974.509381 0.0 9334.755494 4.58491i309 -1,5435S3t7 -31,80Ét{23 1050.859138
,ü95{28271ó 1.29211932 15.f2500298 -,0t{13t630.0 7251.f1570ó -ó212.13357 24123,289É7 -18994.9709 -5112,3iA2ú 139t0.05f72 3707.t59249
33ó.9ó363 783.5056919 0.0 24092.2170ó -r,7582ú957 2.4É2729519 -3t,806tf23 7á4,732678i-,I04ó3477 -.98ÉEú9578 15.38r892i -,03854{t2t
0.0 2399,265008 1711.0580É9 -7ó74.8454S 3ú8.25r98 ó?90.399775 -{+{9,2{311 1982.13{489-921.f{É301 -956,346203 0.0 -7705.50272 -3,7t054349 2,97808798ó -3t,B0Alt23 I129.055793.0913ó92202 1.282275799 15,53598579 -.03884{348
0.0 -4789,{5771 6ó00.ri7f2ó5 -39298.907 -29836.4495 1187ó,39llf -227il7,918 4ó16,888332
1092.16858 -1ú60.8f152 0.0 -39330.8012 -.rg{92909{ -.?81517455 -31,80Él{23 633.751{71f.5155805399 2.897bt5t27 t3,28b8i777 -.047{3083?
0,0 -á497.71241 3135.8612 -29301.1092 -15430,7954 -18315.399{ -16935.86€8 3405.9{4r75ó5.19393831 2777,235899 0.0 -29332.353 2,3óú053873 -3,79828{ó9 -31,80¿1t23 8ó1.3388257-,rÉ7945583 -J.560273{3 15.23{41295 -.04i{0ó355
r'.0 lÉ0ó.4641t5 -5019.87797 r2751,07909 -541ó,2?013 8962.950285 7343.840876 2479,97097
-54ó.290332 -1332.17228 0,0 t2720,78266 r.159800ó63 -2.03ó9158 -3r,806t423 1009,0{162.l58fltifót 1,7b97b922J lÍ.180315Éá -.0{ó8386
0.0 7038.0f5602 -6216,22774 28286.6{098 -22{80.3037 -Í18f,21728 163t3.72748 397i.630893575,83ó70f8 S19.2ó41377 0,0 28255,7f03.A3919tó205 -.17tú323f5 -3I.806t{23 h70,t7ü67-,0915t0t77 -.98315225{ 15.123{19{3 -,051{3ü851
0,0 2623,421858 1455.478ó1f -2{08,9?187 192.4f9290ó t982.Í7ó35f -1{ü8,9557t i9ó9.4{582t-910.ú237il -315.528171 0.0 -2f39.83559 -1.7093ó571 1.1620É5531 -3t,S06t423 1092.30{785.019834ó626 .¡9f5343552 t5.06{?5ó93 -,050358919
0,0 -{558,982t9 ó{1r.205837 -3{612.6332 -25845,0292 uf96.71912 -20002,2593 {227.5464ü8
840.2150t41 -1628,33317 0.0 -34642,8504 -?.1ó5ó0882 ,08ó597f0{3 -3r.80ó1423 É65.0É70381
.{ff6t047?8 2.É0997?3?2 15.0023á3{7 -,0572ó9555
0.0 -b44b,2J7t7 330r,353147 -2ü603.2t06 -lrü87,602ó -ló035.6794 -15378,?4ú2 329r.18102841.9125828{ 2f30,3{3$94 0.0 -26ó33,7967 1,30552003ó -3.f125{92 -31.S0ó1423 839.47972ú9
-,330248551 -3.05f88499 r4,93795294 -,0t739959t
0,0 1258.É8ó{3 -46ó9,28014 t3733.19342 -5f23.813i1 9803,38744 7127,36237 2t31.6?1878-521.5ff521 -r{á6,ó701 0,0 r372r,80ó37 2,888752231 -?.7É097855'31.806tf23 989.3792177
.ts4653029ó t.920590171 1{,8i17991f -.057{ú7911
0.0 6813,358233 -ó2ó4.20344 3032S,52537 -2{329.140t -{f7f.04147 t7492.58601 {093,526733
7r9,39s373 727,2ú50363 0,0 30297,25433 ,9342i7973 -.U089952J -Jl,B06t{23 609.3775695
-.0É0592497 -.ó608655f4 14,803013?8 -,062300294
0,0 2883.3357ó5 1133,ó32357 t2{7,59298f -137.0{98É2 -99{,358252 702.1915009 l9ú2.088087-87S.3ó0ó tzB,1451681 0.0 1217,03813 -3.0985851 2.26f490905 -31.8061123 1060,5ó5ó{9
-.031t95887 -.19139252J 14,732151s9 -.060á3019{
0,0 -4274,71109 620ó,fl9702 -31419,3093 -22934.il98 1I550,15748 -181t8.5742 3930,319595
6ó5.14{3tó4 -ló51.75558 0.0 -31{50,1797 -1.91113ól.i177238302 -31.8061{23 É88.183052ó
,4107220f86 ?.48100551S lr.ú5911765 -.066478{31
0.0 -6{3ó,197i 3525,820S41 -25f75.ólt3 -148ói.2108 -t{700.3045 -I{78f.955 3259.9{4{Éó77.9306150t 2229.590S4ó 0.0 -25607.t82f 1,94929028f -3,0ó005331 -31.8061f23 806,9832121-,312702125 -2.734ó4ltl l{.58335829 -.0ó7218361
0,0 909.232729 -4330,16909 13375,97074 -f781.81902 9790.0{7311 7704,551f74 2330.98úrlB5-142.92729 -1463.t9378 0.0 13344.32893 2,3411i4309 -1.99550515 -31,80É14?3 978.526ü624
.187r2306{4 r.8773969t4 14,50ú10051 -.0ó7130il90.0 É550.5103tl -625{.0370{ 311ó9.79353 -25211.5É73 -3279.1{228 17978,2ó935 {128.6188ói
798.2052473 55ú.2511?42 0,0 31t38,30547.5192t61213,{5t?27790A -3t.80ó1f23 5ó8,7098684-.01f8{9087 -,379506399 t{.f2633529 -,07212793b
0.0 3176.0722ó{ 759.0537ó78 {104,5ó5417 -ó4f.864?Í{ -32ó?,É83ó7 2351,688012 l9ó9.215!85-831.57{401 f70,0048905 0.0 4073.25t{97 -2.59{24ó0i 2.5858á9728 -31.80ól{23 1028,89É32¿-.072S55f58 -.620078f3{ 14,34f57207 -.07025?t9?
0,0 -3938.99339 597f.ó5f443 -2895É.ll19 -2051{.028ó 11804,{19{ -1ó73á.f2t2 367ó,ó0ó889322,975232J -1i00.84597 0.0 -28987.54t2 -1,ó963f019 .9021521t83 -31.8061423 707.2ú5592i.392729f973 2.4?1988479 14.26087089 -.075300806
0,0 -ó44f.33ú87 3789,7898i8 -25f48,{ü52 -15ú04.45Ér -13758,8224 -lf7lt.5l2i 3280,4199ó8149,752É893 2090,698ó7ó 0.0 -25{80,0{79 1,2{t{3689ú -2,28093{úÉ -31,80É1423 76ó,775936$-,103É260áú -2.{9ltftfó 14.t745ó703 -.07Éfüú019
0.0 5iE.?372016 -3979.28811 l2lBl,5ó435 -382f,é99{3 915?.970t9 70I4,93{595 2205.9á632{-5S8,044599 -1370.59108 0.0 12149.É7053 2,1.46471974 -1.69387?{2 -31,80É1423 971.5tr95335
0,6
t700
.174452784{ 1.716569979 t{.ü8690007 -.07ú2üi¡5r.r¿
ú.0 62ii.óó1087 -622t,71412 3tSl5.1ó9Sl -?5490,7598 -1ió9.3637 180ú2.1ú5 41ut.950¿55837.40úB4l 3¡4,619252J 0.0 31283.2{748 .995411ólit .34{587É743 -3L9061{23 542.100914
.ü42?787il1 -.0ú9370995 13.99ó653ó9 -.081284396
0.0 3{96,8J7J0É 333.254$5i6 É583.7J8Ír3 -1t59.60922 -517{.169ó? 3783.0S?85ó 1994.7{6{53-770.137É47 . 758.03&7994 0,0 6Íi2.272t02 -2,593{871f 2,8131002ó1 -31.8061{23 995.Í288578-,r09580853 -.96178932 t3.90{5{453 -.07941928
0.0 -3569.82555 5721,097¿t7 -?6S1f.ú292 -t8273.97t5 12104.ó2288 -15i00.0078 3t39.94I77f392.438{38i -li5+.81535 ü.ü -2ó846.0483 -r.940r8s84 1.18i51886{ -3I.80út423 725.1895362
.38005É5951 2,38583Éó37 13.810ó7716 -.083777É48
0.0 -6449,318 {075,14ó4tt -25i85.792 -tóó74.0346 -1289ó.6593 -14906.3119 3332.09280ó
240.11751 l9ó5.{42óóó 0.0 -2581i,5402 1.082795{49 -2.03223ú82 -31,90É1423 721,3802873-.29581{5ó5 -2.2ú49f651 t3.il{3057 -.0853782{?
0.0 l{f.6ó85465 -3608,78ó0f 10f74.0f338 -27ó3.3?551 B0ó5.391 ó029,066955 207{,50267ó-ó42.17790s -1211.2024 0.0 l0{{2.01707 2,38859{t8É -1.77537983 -31,8061{23 9ü5.7148408
,1515812431 1,t712i0772 13,óló50315 -,084828659
0.0 5932.É9{ú23 -óró3.29098 t1020.9f903 -2530t.9097 -57.5Ét{655 17892.25989 4051.16113{Bf6.Élú8021 80,373489S5 0,0 309S9.0773 1,02ó6{5055 .1992403651 -3t.80614?3 52ú,949ó89.10Ét8198+ó,?53ó8ú5754 t3.51ó{589á -.089922213
0.0 3833.670007 -132.É01959 8897.9{2753 -2315.5{875 -6859,tÉ613 51t9,23ó$28 2042,186191-692,440{89 1013.151243 0.0 8Bt6.5ró921 -2.559óSiU 2.820525888 -31,S061{23 959.3t83{5i-,1437912f4 -1.2f710894 13.4144A847 -,088209412
0.0 -J171.98091 5{50.412182 -24770,82?L -ló0ó1.862 12333.I72É5 -1{319.9884 32ü7,2429182ó3.073É029 -1797,32121 0.0 -2480?,2t03 -1,832iúB4S l.15572ó206 -51.80ótf23 743,2508764
.3óó922990{ 2.3{f06S752 13.31087014 -.091914875
0.0 -ó431.725ó8 43ú3,9ÉBfó{ -2ó330.28{5 -l79lg.ó937 -t1924.9591 -15220.6794 3400,65899t342.7791289 182{.9il76 0,0 -2É361.9ú72 .82891717+8 -1.85189{73 -3r.80ó1f23 É72.85t7873-.28390ó833 -?.02283795 13,20485885 -,09395580t
0.0 -27ó,2907ól -3209.73095 8{19.79115t -U44.1S99 áó22.É6ó2ft 48{3.0130á4 1949,023074-ó95,14f789 -999.0últ7 0.0 8587.88848f 2,3227794,l -1,8t657467 -31.8061{23 959.4237ó18.121350?195 t.188503525 13.097{{195 -.093089047
0.0 5574,312508 -6074.41936 30f12,13707 -2f7fr.5229 t772,893303 t75f0,72083 3961,23t711830.Cs50597 -193.{20Éü4 0.ü 30380,305{ L,l7b92229t -.0ü?55023 -31.8061423 522.rór5r55.1i33653{S5 .57t3{3t074 t2.98792903 -.09gil578ú
0,0 u78,492942 -ú33.71{93 111f2,87027 -3539.8r35 -835J.2{05ú t415.395J72 211{.083878-59ó.93¿5É2 t240.751413 0.0 11il1,{44?t -2,36129526 2,707824713 -31.80ó1f23 919,Á0fÉ701-.177005873 -t,48519308 t2.8ió47i18 -.096663{li
ü,0 -27f9,3012 5155,i58917 -22701,77h? -13815,75t't 12395,39255 -13125,39ó3 2973,217294
130.5357209 -lgt4.g10ó6 0.0 -22733.r63i -1.858áló74 1.16898{3É7 -31.8061423 7&1.7t96437
,350il?799ó 2,278?343h7 12.7ó3547t -.0997205t70.0 -ó378,31992 {É41,170121 -2ó918.9168 -19210.0ó59 -10731.8655 -15560.52{t 3474,518412
4{7.9626102 1451.óó3624 0.0 -2ó9-10,5008.5019876f9t -1.ó20{1609 -31,80Ét423 623,217{577-.2ó3I9ó8f2 -t,71770123 t2.64829ó2f -,102195739
0,0 -724.873218 -27tb,8l7Lb 6113.48óil3 -880.803718 t$gÉ,720726 35It.431029 1839,20092S-739.1149ó9 -742.808741 0.0 ó091.6110{2 2.37t083t9 -2.00517t5ó -31.80át{23 951.2988118
.08547927{f .8586888301 12.53tA5{98 -.t01028ú88
Ertrena of output variables lollor, 0rder and colurn positioninq are the sere es for the preceding tire-step loop output.lleriable raxi¡a : 9339,É9332{ 9340.060987 9339.5808{1 31373.31104 3749,487732 31869.85479 18I06,2793 6308.308113
3{f2.944051 5261,71ü27 0.0 3t3¿0,07t77 17.802ó804 t8,7778705 -31.8061{23 t{14.51987{.95ú2331f51 8.31278234 15.ó2501107 .t07Bl99E-S
,É2
,É+
.bb
.ó8
0.7
11
.7{
,7ü
.78
0.8
3100
4100
.84
,8É
,88
0.94500
Ti¡es of ¡¡xira ; .0{84.041
.0ó
.020ó
.0a34
.058ó
.0428 .Í532 .0558 .0582
0.0 . fó98 ,050{ .0548
0.0 ,0{98
.0178 .0612 .0Ét .063{0,0 .0ó12 .ü50 ,0528
0.9 ,8956
,5ó98 .0ól0.0 .0562
V¿riable rinina : -9339.72309 -9340.0t93 -9339.52709 -b7427,944 -5f921.t3f -3{?S3.33?3 -38945.95ó2 l0Il.9ió0ú-98J.753132 -f593.8047ú 0.0 -67470,9498 -18,00581{7 -18,5350587 -31.S06t423 272.i2t932-.ótlós3885 -1,4749707 12.53165498 -,10418845{
Iines of rinira I .023{,t922.0b26
,045É
,0382
.0ó3É
.0É12 ,0558
0.0 ,0ól{
Plot c¡rd ignored in search for blank terrinetor, i PRINTER PL0T
IFIU
lD litulo de la qrelic¡i[ 3-78rc(U
ter¡inator. i ? TRAI{SIT0RI05 DE LA IíAüUIliA Slll[R0}llCA TfiIFASICArf
terrin¡tor. i 14{ 2. 0.ü ?0 ñACAI frACBl ilACCI
iBLAilT CARD EHDII{E PL[]T IAftDS
llerory storage figures for the precedingr nor-colpleted data case.
A v¡lue ol '-9999n indic¡tes th¡t no figure is availeble'Si¡e List 1. llurber of electric netrork nodes'
Sire List 2, l,lu¡ber of electric nettorl brenches,
Size List 3. Nulber of dat¡ values ir, R. L, I tables.Size List 4. ilurber of electric nettorl sourtes.Si¡e List 5, Storage for [Y] and triangulari¿ed [Yl. llo. tires = ? Factors = l5Size List ó, llu¡ber of entrieE in sxitch table. llo. flops = I5i¡e List 7. l{urber of distinct ALPHAllUltEftlE data n¡res plus prograr SPY veriables,Size LiEt B, History points of distributed lines,Size List 9. l{urber of nonlinear elerents.Size List 10. Points of nonlinear rh¡recteristics.Sire List 11. ilurbpr of TIpe-59 5,11. outputs.Si;e List 12. Total nu¡ber of EHTP output veriables.5i¡e List I3, llorking space for batch/SPY plotting.5i¡e List 14. S,ll./U.lf, connections to TACS.
Si¡e List lf. Hen TAES character starage in bytes.
Si¿e List 16. Iotal nu¡ber ol lype-59 S.ll, tasses.
Size List 17. ilurber of Type-59 Synchronous ¡achines.
Si¡e List lS, [rench and sritch poner/energy outputs.
Size List 19. Total llaating-point TACS teble spece.
Size List ?0. Non-rogied recursive convolution deta,
Size List ?1, Totrl rodal/phase IT] ratrix storage.
Si¡e List ??, Total recursive ronvolution histarv'Size List 23. 6iant vectors for renurbering, phasors.
Size List 2{, Peal phases of corpensation lor d¡t¡'Si¡e List 25, Tot¿l teble space for all U.ll usage,
Si¡e List ?ó. Square of lax nurber of coupled phases.
Tiring figures cherecteri:ing central protessor {CP} solution speed.
Data input tire {through the top ol 'SUBRó") ,..'l{ode renurbering and phasor solutionAlter phasor solution, but before tire-step loop
Integration ol equations (tioe-step loop)Plotting or STATISTICS tErrination overlays ,...
Co¡nent card. ttUt'lDü[ = 56??,
Csnrent rard, I'lUllDCD = f¡2.1.
Colrent card, l'lUlltlCD = 5É?4.,
Cnment card. NUI{lltD = 5615,
Plot card ignored in search lor bl¡nkCo¡¡ent card, HUIIDID = 3ó17,
Ploi c¡rd ignored in search for blaniBlank card terrinating all plot cards,
llarler card preceding ner EIITP date c¡se.Blank card to terrinate EIITP exerution,
Present Progrea
ligure lirit7 2r4ó 25ft2 424
s4257 2124
l13É 594
ü 1487
0ú30 135
tó 42
20 ¡0-9999 r-9999 tó-9999 I
0l{030l
23 lÉ00
0 óó0
0 100
0 150
15 J752
05-9999 400
-9999 42
CP sec li0 sec
20,430 0.000
il,980 0.000
6.6fü 0,000
88.820 0.000
?,470 0.ü0ü
f narel
I LBUS )
{LBR}ICH }
{ LOATA }
( LEXIT }
{LYIIAT }
{ LsHTCH }
(Ls¡ 1E7 )
{ LPASI )
tLilolrL l
{ LDHAR )
(LStt0iJT )
(LSI ¡I? )
(rsr¿r3t
{tBsTAC }
{ LCTACS }
{LrnAss}(LSYil )
(IIATPE }
{LTACSI )
(LFSEH}
{Lt0}(LHIST}
(LSI¡?3}
{Lt0lrPl(tsPcuilI
{Lsl 12ó }
Sul sec
20,430
11.980
ó, ú40
88,820
2.t20
Totals 130.290 0,000 130,290
Alternative Transients Frogral (ATPl, IBll PC II transletion, Copyright 1987, Use licensed only through LEC (f.U. Leuven, Belgiur),Date (dd-¡th-yyl and tire of day (hh,m,ss) = 06-Apr-93 20:40:27 l{are of disl plot file, if any, is [:3f062040.p14
Far infonation, consult the copyrighted ATP EllfP Rule Book published by LEE in July, 1987, Lest r¡jor prograr update; l{ay, 1988
Total length of "LAE[0ll'tables = 45197 II{TEEER rords, 'VARD¡ll' List Si¡es follor: 214 2i4 124 42 2124
33 594 1{87 63 135 42 50 I 1ó r 34 3 5 t600 úó0 t00 150 3399 3 400 42 42
-----------+-Descriptive interpretation of input data cards. i Input d¿ta card irages are shorn belor, all 80 colurns, chrracter by character
0123456180123{5679901 23156/890123{56i890123{3ó7990121{5ó7890123{5é7s9012I45ó7890123+56799['
|BEGIII I{EH OATA CASE
:FLAII(
tkI
tkI
t
I
IJfI
I ilff
t
lrll\f,ll¿q'
- l0
_ t 1l llfiCri 1
lDú 1t0 l¡tr]
..-.-.- f f¡ ¡tñccl
FaiIa Jlonofásica del ffenerndor Sincrónico (Fase a)Fen iúdo Subtransiton io.
nü Í'-i
--- ( 2l il8tEt
?SE 3r'JL1
---- ( ?¡ nficBl
a{{{thlr
l¡tltt
39C ¡l$r1
..-.--- ( 3r ilRCCIt'tficfit
Falla Hnnofásica del Generadon SincnónicnFer iúilo tr.ans itsr. is.
06-Apr'-93 28i82,.56
66-Apr93 ZBiBZJ56
rLtL ¡I¡ BÉ-fipr93 Z0:BZ:56
-'¿
-
( t)
r-_ t-! D
--- ( r! HficEl
ic0
'-,-.- ( 31 rRCCf t [nrsecJnRcR 1
Falla llonofásica del tenenador SincrónicoPeniódo permanente.
tkl
n ItJfJ
-- 11?¡ nftcH i - Ifi
tkl
0
-
¡. 5l ilficH r - ID
türl ?üü ¡tuü
--- t¡i.i ttnln ¡ - ¡E
9n0 EDo
-'-.- (l+i itftct{ i - Ic
151-r t00
.-.-.-(finREHl-I0
Falla Jlonofásica del Senerndor SincrónicoCornientes Ia, Ih, Ic.
Eñ I ñN
----t 6) nRtH t - r0
Falla Hnnofásica del generador Sincnónicoüoruientes de Armaduna e¡r los ejes d, {r 6Per iórlo Suhtrans itcr ic.
66-f,pr93 ZÜiBZi56
B6-Apr*93 Z0iEZt56
Iki
-fC'
-
I El ilnr-n ¡ _ ¡u
tttItt
. .1.
I
li
_<fu J!,u st!t)
--- .. Él ilnr-n ¡ - ¡¡
Eür] E5ú
---(ElnRcHl-IC
T'lü lzrj.-.-.- ( :1 nficH
+ $tr
I - I0
Falla llonofásica del Generador. SincrónicoConrientes de Arrnadur,a en los ejes d, {r 0Periódo Transitonio.
t.
illil:ii!:il
iliril rlli
':.l
¡itI' iilill
'litiIIIIti
-{o5LlLl
gi nstH I
55f
IO
7ü0 ?5ü
'-.-.- i ir nficH I
Faila Honofásica del ffenenndor. SincnónieeCorrientes de Armadura en lss ejes d, {r BPerióds Permanente.
B6-Apr93 Z0:fiZ:5É
BÉ-Apn93 ZBig2:56
tl E6-Apr-93 Z0:BU:56
t [msecl
Éüúú
50üD
+ 000
3000
10110
U
- 1000
- 2r-100
_ ( $l
!l ll1tl tu0 30t| +0ü
ilnLn ¡ - tt
Falla lfonofásica del ffenerado¡r SincnónicoConnie¡rte de Canps (IF).
0 lrlr.i tüü 30ú ttC'l
t -¡ ttnlfl ¡ - l¡au ¡ 1, ¡¡ñLn I - ^hu
Falla llonofásica del Generadon SincnónicoCorrientes en los devanadss ftilortiguadores.
0É-fipr-93 ZB:62:56
tkl
Ttl Stl
(itiHfi[Hr-IG
Sincrónico( IG}
ÉÉ-Apr-93 Z8:E?i56
Falla tlonofásica del Genen¡donConriente de falla etr el ¡.stsrPer ióilo Subtrans itsr io.
¿90 30rl
nficH I - IG
Falla lfonofásica del Genenadortsnriente ile falla en el rotorPeriÉ¿{s Transitnric.
Sincrónico( IG)
B6-Apr93
- (t0¡ t [msec]
tkl
7!lrl 7911
ffÉ-Apn-93 ZB:BZ.:5É
B6-Apt*93 ZBiB2:.56
ilril llil+IfillI ll ,
ftlrl
5r'ú 55t' É0'l E5ú
1¡U' ttnLn r - tu
Falla Flonofásica del SeneradonConrie¡rte de falla en el roturPeriódo Permnnente.
Sincrónico( It)
tI
l{L1Ct
1t00
100 D
¿ ü'l
lr¡JU
IJ
| ¡9-r
0 1L'10
nficH l - ñF0RcE
Falla llonofásica delFuenza llagnetmstríz
Genenailor Sincnónicotstal en el entrel¡ieruo.
t [r¡rsec]