UNIVERSIDAD NACIONAL DE COLOMBIA

SEDE MEDELLÍN
FACULTAD DE MINAS
ESCUELA DE INGENIERÍA ELÉCTRICA Y
MECÁNICA

GUÍA DE PRÁCTICAS
FENÓMENOS
TRANSITORIOS
ELECTROMAGNÉTICOS

LEONARDO CARDONA C.
Profesor asistente

Medellín
Marzo de 2004
 2

1. INTRODUCCIÓN A LOS FENÓMENOS TRANSITORIOS

1.1 INTRODUCCIÓN

Un sistema eléctrico de potencia está sometido a una serie de perturbaciones

que alteran su estado normal de operación. El paso de un estado a otro
ocasionado por una perturbación se hace en forma gradual, ya que las
variables físicas como voltajes y corrientes no pueden tener cambios bruscos
debido a que las leyes de la Física no lo permiten. Todo el proceso de
cambio de un estado de operación a otro es lo que se pretende estudiar en
un curso de Fenómenos Transitorios y corresponde a lo que se llama el
régimen o estado transitorio.

El estudio del estado transitorio debe hacerse en forma cuidadosa ya que las
variables físicas pueden llegar a tomar valores extremos y exigir a los
equipos, hasta ocasionar el deterioro de los mismos o la interrupción del
suministro de energía. El conocimiento de los estados transitorios permitirá
tomar decisiones para proteger adecuadamente los equipos, lo mismo que el
aseguramiento de la calidad de la potencia eléctrica, con unos estándares de
calidad adecuados.

Los transitorios en un sistema de potencia en forma general son de tipo

eléctrico, mecánico y térmico. Los transitorios mecánicos y térmicos, a pesar
de ser muy importantes, son mucho más lentos que los de tipo eléctrico. Su
frecuencia está por debajo de la sincrónica (60 Hz). Los transitorios eléctricos
son muy rápidos y constituyen el motivo del trabajo en esta asignatura.

El sistema eléctrico, desde el punto de vista circuital, está formado por tres
clases de elementos: resistor, inductor y capacitor.

Cualquier componente de un Sistema Eléctrico de Potencia (SEP), estará

formado por la combinación de elementos RLC. En estado estacionario un
parámetro, para un determinado equipo, será dominante sobre los otros, pero

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 3

en estado transitorio hay dependencia de las características y tipo de

fenómeno transitorio. En un transformador de voltaje, por ejemplo, es claro
que la componente inductiva es dominante en un estado normal de operación
a 60 Hz y para eventos transitorios de una frecuencia relativamente baja,
pero se sabe que en eventos cuya frecuencia está por encima de unos 5
kHz, el efecto capacitivo empieza a ser importante y el transformador pierde
su precisión.

Otra característica importante es la forma como los efectos RLC se

encuentran presentes en un equipo. Se pueden presentar de manera
concentrada o distribuida. En una línea, por ejemplo, los parámetros RLC se
encuentran distribuidos a lo largo de la misma, mientras que en un reactor
para compensación de líneas, el parámetro dominante inductivo se encuentra
concentrado. El modelo para representar estos elementos dependerá de si
se puede representar de manera concentrada o distribuida. Estos
parámetros adicionalmente pueden depender de otra variable, como es el
caso de la inductancia, que representa la magnetización de un
transformador, la cual depende del voltaje aplicado.

1.2 CLASIFICACIÓN DE LOS TRANSITORIOS

Los transitorios que aparecen en un sistema de potencia se pueden clasificar

desde diferentes puntos de vista. Uno obedece a los tipos de fenómenos
físicos involucrados; para este caso, los fenómenos transitorios se clasifican
en dos categorías:

• Transitorios que resultan de la interacción entre las energías

almacenadas en capacitores y en inductores.
• Transitorios que resultan de la interacción de la energía mecánica
almacenada en las partes rotóricas de las máquinas y la energía
almacenada en los circuitos.

El interés en estas notas son los transitorios electromagnéticos

En la Figura 1.1 se ilustra esta clasificación de acuerdo con la frecuencia del

fenómeno.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 4

|_____|_____|_____|_____|_____|_____|_____|_____|_____|_____|
2 3 4 5 6 7
0.001 0.01 0.1 1 10 10 10 10 10 10 10
FRECUENCIA (Hz)

|← Fenómenos Electromecánicos |
→ ← Fenómenos Electromagnéticos → |
Control carga Estabilidad Corto Maniobra Ondas viajeras
frecuencia transitoria circuito de líneas

Estabilizadores Resonancia Voltajes TRV

Subsíncrona

Armónicos

Figura 1.1 Clasificación de los transitorios de acuerdo con la frecuencia

1.3 OBJETIVOS

Los objetivos básicos propuestos en el curso de fenómenos transitorios son

los siguientes:

• Comprender la naturaleza de los transitorios electromagnéticos en

sistemas de potencia.
• Conocer diversas técnicas de análisis y simulación, como también sus
aplicaciones.
• Obtener los conocimientos introductorios que resultarán indispensables
para cursos de coordinación de aislamiento.

En forma práctica, los ingenieros analistas no sólo deben conocer el

fenómeno transitorio, sino dar soluciones que apunten a mejorar la calidad de
la potencia y la confiabilidad del SEP. Para cumplir con lo anterior se deben
seguir los siguientes cuatro pasos:

• Conocer el fenómeno: Origen, formas de onda, magnitud, duración,

frecuencia de ocurrencia, evolución del fenómeno.
• Conocer la respuesta de los elementos a este tipo de excitación.
• Coordinar adecuadamente la ocurrencia del fenómeno transitorio con las
características de los elementos del sistema.
• Proteger adecuadamente, mediante el diseño, pruebas y aplicación.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 5

2. INTRODUCCIÓN AL PROGRAMA EMTP/ATP

El programa de transitorios electromagnéticos EMTP (Electro Magnetic

Transient Program), ó ATP (Alternative Transient Program), es un programa
digital utilizado para simular transitorios electromagnéticos, electromecánicos
y de sistemas de control en sistemas eléctricos polifásicos de potencia.
Inicialmente fue desarrollado junto con otros programas, como alternativa
digital del Analizador de Transitorios en Redes (TNA - Transient Network
Analyzer). Posteriormente, durante un período de quince años, se lo ha
ampliado considerablemente, resultando actualmente un programa de amplia
difusión y utilización en todo el mundo.

El programa ATP con respecto al curso de Fenómenos Transitorios

Electromagnéticos tiene una aplicación clara y es la de simular todas las
situaciones que generen esfuerzos en forma de sobrevoltajes y/o
sobrecorrientes. Estos esfuerzos exigen los aislamientos de los equipos,
generan efectos térmicos y pueden producir también esfuerzos mecánicos.

2.1 HISTORIA

El EMTP/ATP fue desarrollado en su versión inicial por el Dr. Hermann

Dommel durante la década del 60 en Alemania, quién posteriormente cedería
el programa a la Bonneville Power Administration (BPA) de los Estados
Unidos. A pesar de ser la BPA la que coordinó la distribución del programa
entre los diferentes grupos de usuarios, muchas universidades contribuyeron
al desarrollo de los diferentes modelos que constituyen lo que hoy es este
programa.

El programa ha evolucionado desde las versiones para grandes

computadoras hasta las versiones actuales para computadores personales.

Debido a que el programa fue escrito inicialmente en Fortran, la interacción

con el usuario es mediante un rígido archivo que debe cumplir ciertas normas
de dicho lenguaje. Actualmente existen programas tipo interfaz, que

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 6

permiten un trabajo más amigable para muchas de las aplicaciones. La

interfaz que se utlizará es la ATPDRAW, que es un producto de Norwegian
Electric Power Research Institute, mediante un contrato con la BPA. Esta
Interfaz ha evolucionado bastante desde una versión inicial para sistema
operativo DOS hasta las versiones actuales para ambiente Windows. Esta
interfaz cubre una buena parte de las posibilidades que tiene el programa
EMTP/ATP pero siempre se deberá tener un conocimiento básico de cómo
opera el ATP, que es realmente el programa simulador.

2.2 ESTUDIOS A REALIZAR CON EL PROGRAMA EMTP/ATP

A pesar de que el objetivo principal del programa es la obtención de la

respuesta transitoria de los sistemas eléctricos de potencia, también se
puede obtener la respuesta en estado estacionario para un sistema
alimentado con CA (respuesta fasorial).

Para la simulación del sistema de potencia el programa posee varios

modelos que deja una representación adecuada de los diferentes elementos
que lo componen.

Los modelos disponibles se pueden clasificar así:

• Modelos constituidos por elementos concentrados RLC. Estos modelos

pueden ser simples ramas RLC serie, circuitos PI polifásicos que pueden
representar líneas de transmisión o transformadores.
• Modelo de onda viajera para representar adecuadamente una línea o un
cable.
• Impedancias de tipo no lineal: Resistencias no lineales, inductancias no
lineales, resistencias variables con el tiempo.
• Suiches ideales controlados por tiempo. Suiches controlados por voltaje
para la simulación de "gaps" en flameos de aisladores. Diodos y
tiristores. Suiches de operación estadística.
• Fuentes de voltaje y corriente ideales de tipo escalón, sinusoidal, rampa,
exponencial y definidas punto a punto.
• Modelo completo para la máquina sincrónica.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 7

• Modelo universal de la máquina eléctrica que permite representar doce

tipos de máquinas diferentes. El modelo de más interés en este grupo
de modelos es el del motor de inducción trifásico.
• Representación de los sistemas de control mediante la opción de los
TACS (Transients Analysis Control System). En las versiones más
recientes del ATP existe una opción alterna para representar la parte de
control del sistema de potencia: MODELS. En la Figura 2.1 se observan
los diferentes módulos de que dispone el programa ATP para generar
modelos. A estos módulos se les ha denominado subprogramas de
soporte y muchos de ellos se pueden manejar directamente desde la
interfaz Atpdraw.

Figura 2.1 Programas de soporte que interactuan con el EMTP/ATP

El ATP permite hacer entre otros los siguientes estudios:

• Cálculo de parámetros en líneas aéreas y cables subterráneos

• Cálculo en estado estacionario sobre un sistema de potencia
• Sobrevoltajes por maniobra (suicheo)
• Sobrevoltajes por descargas atmosféricas

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 8

• Cálculos de cortocircuito involucrando los sistemas de retorno

• Coordinación de aislamiento
• Cálculo de voltajes inducidos sobre elementos cercanos a líneas de
transmisión (otras líneas, cercas, oleoductos)
• Resonancia en líneas paralelas
• Simulación de arranque de motores
• Evaluación de armónicos
• Estudios de ferrorresonancia
• Máquinas eléctricas
• Simulación del control de las máquinas y del sistema de potencia

2.3 ESTRUCTURA GENERAL DEL ARCHIVO DE DATOS PARA UN

CASO SIMPLE DE SIMULACIÓN TRANSITORIA O DE ESTADO
ESTACIONARIO

A pesar de que el programa permite simular elementos de un sistema de

potencia con modelos muy elaborados, en este primer acercamiento al
programa ATP se utilizarán únicamente ramas simples del tipo RLC con
parámetros concentrados y fuentes simples del tipo escalón y senoidal.

Las partes del archivo de entrada de datos para este tipo de simulación en
forma general se observan en la Figura 2.2

Un archivo de entrada general al ATP se caracteriza por los siguientes

grupos de datos:

• Datos iniciales para identificar el inicio de un nuevo caso, identificación

de sí es una simulación, o es un caso de cálculo de parámetros para
algún elemento de la red (acceso a subprogramas).
• Si es un caso de simulación el siguiente grupo de datos corresponde a la
especificación de tiempo de simulación, tiempo de paso de integración,
especificación de las unidades de los parámetros para inductancias y
capacitancias, control de salida de información. También se especifica en
este bloque de datos si se requiere simulación en estado estacionario.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 9

• Datos de las ramas de la red (modelación de los componentes de la red

de potencia).
• Datos de los interruptores.
• Datos de las fuentes.
• Especificación de las condiciones iniciales si existen.
• Lista de voltajes nodales a ser tabulados y graficados.
• Bloque de datos para unas opciones gráficas y de cálculos especiales
como análisis de Fourier y tablas con resúmenes estadísticos.

1 2 3 4 5 6 7 8
1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0

B EG I N NE W D A T A CA S E
1 2 3 4 5 6 7 8
1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0

Delta T Tmax Xopt Copt

1 2 3 4 5 6 7 8

1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0

Network Steady-state Extremal Extra printout

Print Freq Plot Freq MemSave Plotted output
connectivity phasors values control

DATOS DE RAMAS
/ B RA NCH
1 2 3 4 5 6 7 8
1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0

ID BUS1 BUS2 BUS3 BUS4 R L (WL) C (WC)

DATOS DE SUICHES
/ S W I T CH
1 2 3 4 5 6 7 8
1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0

ID BUS1 BUS2 Tclose Topen Imargin

DATOS DE FUENTES
/ SOURCE
1 2 3 4 5 6 7 8

1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0

ID BUS U/I Amp f Pha A1 T1 Tsta Tstop

/ I N I T I AL
/ OU T P U T
ESPECIFICACION DE VOLTAJES NODALES A CALCULAR Y/O GRAFICAR
1 2 3 4 5 6 7 8
1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0

NODO1 NODO 2 NODO3 NODO 4 NODO 5

B L A NK B RANCH
B L A NK SW I T CH
B L A NK SOURC E
B L A NK I N I T I AL
B L A NK OU T PU T
B L A NK P L O T
B EG I N NEW D A T A CA S E
B L A NK

Figura 2.2 Formato general para el programa EMTP/ATP

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 10

En la Figura 2.2 se observa la forma como están organizados los diferentes

datos en los campos del archivo de entrada. Estos espacios en las
diferentes columnas se deben conservar de manera rigurosa.

Para una simulación en particular no necesariamente deben aparecer todos

los bloques señalados. En la Figura se observa la presentación en detalle
los formatos de las diferentes líneas de un archivo para simulación transitoria
o de estado estacionario de un caso simple que no contiene bloques de
control.

Un caso simple de simulación se entiende como aquel donde no se hacen

estudios de tipo estadístico, no hay bloques de control, no se utilizan
subprogramas para cálculos de parámetros de elementos del sistema, solo
se utiliza un tipo de rama (rama RLC serie de parámetros concentrados) para
la modelación de la red, el interruptor que se utiliza es controlado por tiempo
en forma determinística, se utilizan fuentes sencillas como escalón y
cosenoidal.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 11

3. PRÁCTICA N°1. ANÁLISIS DE SISTEMAS SIMPLES RLC Y RL EN

ESTADO TRANSITORIO

Para realizar una introducción en lo que es la utilización práctica del EMTP/ATP

como herramienta de simulación, lo más apropiado es hacerlo con sistemas
simples RLC y RL, pero que involucren estados transitorios de interés para el
curso. A partir de circuitos RLC y RL se puede entender los primeros
conceptos asociados con los fenómenos transitorios reales que se presentan
en un sistema eléctrico de potencia.

3.1 OBJETIVO

El objetivo de esta práctica es familiarizarse con la forma de utilización del

ambiente de trabajo ATP/ATPDRAW/PLOTXY. Se pretende también conocer
el formato del archivo de entrada de datos para un caso simple de simulación
transitoria. Con los casos que se van a simular se pretende realizar una
introducción acerca de fenómeno transitorio, pero orientado hacia las
situaciones reales que se presentan en los sistemas de potencia.

3.2 DESCRIPCIÓN

El circuito RLC serie es una primera aproximación a la modelación de

componentes de un sistema de potencia, como una red de transmisión de
energía. En forma real estos efectos, resistivo-inductivo-capacitivo, se
presentan de manera distribuida y no concentrada como es el caso que nos
ocupa en esta aplicación. Sin embargo hay una serie de definiciones que se
pueden hacer sobre un circuito RLC serie que aparecen definidos cuando una
red de transmisión se modela mediante parámetros distribuidos, como es el
caso del la impedancia característica de sobretensión.

El circuito a simular se observa en la Figura 3.1. Para observar los estados

transitorios de voltaje sobre el capacitor y de corriente, se pueden utilizar
cualquiera de los dos circuitos #1 y #2.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 12

El primer circuito incluye un suiche conectado entre nodos FTE y NSW, que se
cierra en el momento que se desee. La fuente se activa en cualquier momento
mediante el parámetro Tsta, que es el tiempo de activación de la fuente.
i(t)
FTE -------> VC

V NSW vc(t)

CTO #1

i(t)
FTE -------> VC

V vc(t)
CTO #2

Figura 3.1 Circuito RLC serie

El segundo circuito no tiene suiche, pero se consigue el mismo efecto, al

activarse la fuente en un tiempo igual al que se cerró el suiche del primer
circuito.

Para el caso de una respuesta subamortiguada en el circuito, la corriente tiene

la siguiente expresión:
1
1-
V R
t × sen 4 λ2
i(t) = × e- 2L
.
.t (3.1)
1 LC
Zc 1-
4 λ2

Zc = Impedancia característica de sobretensión

L
Zc = (3.2)
C

λ = Factor de amortiguamiento;
λ = Zc/R (3.3)

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 13

Cuando el efecto resistivo es despreciable, la corriente tendrá un

comportamiento senoidal de amplitud constante (no hay amortiguamiento). La
corriente tendrá un comportamiento en el tiempo de acuerdo a la siguiente
expresión:

V 1
i(t) = × sen .t (3.4)
Zc LC

La anterior expresión permite dar una interpretación de la impedancia

característica Zc, como la relación entre la amplitud de la fuente del voltaje de
excitación y el valor pico de la corriente senoidal del circuito serie.

La frecuencia de la oscilación senoidal está dada por:

1
f0 = hz (3.5)
2π LC

3.3 PROCEDIMIENTO

T-cl = 0 R = 10
T-op = 1 L =150
Imar = 0 C= 0

Amp =100
Tsta =-1 C= 1
Tsto = 1

Figura 3.2 Circuito RLC serie en Atpdraw

Haciendo uso de la caja de herramientas de la interfaz Atpdraw elaborar el

circuito de la Figura 3.2 asignando los parámetros sugeridos.

Para editar cada uno de los componentes se abre la caja de herramientas con
el botón derecho del mouse:

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 14

Figura 3.3 Caja de herramientas de la interfaz Atpdraw

Editar cada una de las ramas para formar la red pasiva RLC, los datos del
suiche lo mismo que la fuente tipo escalón de acuerdo a la Figura 3.4. Se
resalta en la figura los parámetros de este par de elementos.

T-cl
T-op
Imar

Amp
Tsta
Tsto

Figura 3.4 Circuito RLC y parámetros importantes de los elementos

Asignar los "settings" correspondientes a la simulación de acuerdo a la ventana

de datos que se observa en la Figura 3.5.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 15

Figura 3.5 Ventana de datos para los "settings" de simulación

De igual manera deben asignarse los "settings" de salida de resultados

Figura 3.6 Ventana de los “settings” para control de salida de resultados

Los valores de los datos enteros que aparecen por defecto normalmente son
suficientes. Se recomienda cambiarlos cuando sea estrictamente necesario y a
medida que se adquiera destreza en el manejo del programa.

Es importante que los diferentes parámetros que controlan las simulaciones se

conozcan bien. Una primera fuente de información es el HELP que aparece en
la mayoría de las ventanas.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 16

Figura 3.7 Persiana "ATP" para realizar el "Make File"

Cuando se tenga completamente editado todos los componentes se realiza un

"Make File", acción que corresponde a la realización del archivo *.atp de
simulación

El archivo de simulación *.ATP se observa en la Tabla 1

Tabla 3.1 Archivo plano generado por la interfaz Atpdraw

BEGIN NEW DATA CASE
C --------------------------------------------------------
C Generated by ATPDRAW Marzo, Domingo 28, 2004
C A Bonneville Power Administration program
C Programmed by H. K. Høidalen at SEfAS - NORWAY 1994-2002
C --------------------------------------------------------
C dT >< Tmax >< Xopt >< Copt >
1.E-5 .05
250 1 1 1 1 0 0 1 0
C 1 2 3 4 5 6 7 8
C 345678901234567890123456789012345678901234567890123456789012345678901234567890
/BRANCH
C < n 1>< n 2><ref1><ref2>< R >< L >< C >
C < n 1>< n 2><ref1><ref2>< R >< A >< B ><Leng><><>0
R VC 10. 150. 0
VC 1. 0
/SWITCH
C < n 1>< n 2>< Tclose ><Top/Tde >< Ie ><Vf/CLOP >< type >
FTE R 1. 1
/SOURCE
C < n 1><>< Ampl. >< Freq. ><Phase/T0>< A1 >< T1 >< TSTART >< TSTOP >
11FTE 0 100. -1. 1.
/INITIAL
/OUTPUT
VC FTE
BLANK BRANCH
BLANK SWITCH
BLANK SOURCE
BLANK INITIAL
BLANK OUTPUT
BLANK PLOT
BEGIN NEW DATA CASE
BLANK

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 17

Las características de este archivo de simulación son las siguientes:

• Las líneas que empiezan con la letra "C" son simples comentarios para
documentar el archivo. El programa ATP en la ejecución las ignora.
• El tiempo máximo de simulación es de 50 mseg. El delta de tiempo es de
10 µseg. Lo anterior da un total de 5000 puntos de cálculo.
• Se va a imprimir una tabla con los datos de las respuestas transitorias de
los voltajes y corrientes en el circuito con un intervalo de impresión de 25,
para un total de datos de 21 puntos en el archivo de resultados. Este
archivo de resultados tiene extensión *.LIS
• Se solicita una tabla de conexionado de la red, lo mismo que los picos
máximos y mínimos de las variables en el tiempo.
• En la columna 64 de la línea correspondiente a datos misceláneos enteros
(línea 3a del anterior archivo, excluyendo las líneas de comentarios) se ha
solicitado mediante un "1" la generación de un archivo con extensión *.PL4,
el cual sirve para que otros paquetes de graficación como el PLOTXY
generen una gráfica, con unas características de resolución muy buenas.
• Para la modelación de los elementos del circuito se utilizan dos ramas tipo
"00" (rama tipo RLC serie). La primera es una rama con valores R, L, C de
10 Ω, 150 mH y 0.0 µF. La segunda rama con valores de 1.0 µF. Las
ramas RLC no puede tener los tres valores en cero. Cuando se asigna
cero a la capacitancia en la rama RLC el programa entiende que es una
rama RL y simplemente asume que no existe el parámetro capacitancia,
porque de lo contrario un valor cero de capacitancia puede interpretarse
como una impedancia de un valor muy grande.
• Como fuente de excitación se utiliza tipo escalón de 100.0 voltios de
amplitud, la cual tiene un tiempo de activación de –1.0 s y un tiempo de
desactivación de 1.0 s.
• No hay condiciones iniciales (las CI son cero).
• Como variables de salida se solicitan la corriente del circuito y los voltajes
de la fuente y sobre el capacitor.
El archivo correspondiente a la salida de resultados de esta simulación se
observa en la Tabla 3.1

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 18

Tabla 3.2 Archivo de resultados

--------------------------------------------------+--------------------------------------------------------------------------------
Descriptive interpretation of input data cards. | Input data card images are shown below, all 80 columns, character by character
0 1 2 3 4 5 6 7 8
012345678901234567890123456789012345678901234567890123456789012345678901234567890
--------------------------------------------------+--------------------------------------------------------------------------------
Comment card. KOMPAR > 0. |C data:C:\ATP\DRAW37\ATP\RLC.ATP
Marker card preceding new EMTP data case. |BEGIN NEW DATA CASE
Comment card. KOMPAR > 0. |C --------------------------------------------------------
Comment card. KOMPAR > 0. |C Generated by ATPDRAW Marzo, Domingo 28, 2004
Comment card. KOMPAR > 0. |C A Bonneville Power Administration program
Comment card. KOMPAR > 0. |C Programmed by H. K. Høidalen at SEfAS - NORWAY 1994-2002
Comment card. KOMPAR > 0. |C --------------------------------------------------------
Comment card. KOMPAR > 0. |C dT >< Tmax >< Xopt >< Copt >
Misc. data. 1.000E-05 5.000E-02 0.000E+00 | 1.E-5 .05
Misc. data. 250 1 1 1 1 0 0 1 0 0 | 250 1 1 1 1 0 0 1 0
Comment card. KOMPAR > 0. |C 1 2 3 4 5 6 7 8
Comment card. KOMPAR > 0. |C 345678901234567890123456789012345678901234567890123456789012345678901234567890
Comment card. KOMPAR > 0. |C < n 1>< n 2><ref1><ref2>< R >< L >< C >
Comment card. KOMPAR > 0. |C < n 1>< n 2><ref1><ref2>< R >< A >< B ><Leng><><>0
Series R-L-C. 1.000E+01 1.500E-01 0.000E+00 | R VC 10. 150. 0
Series R-L-C. 0.000E+00 0.000E+00 1.000E-06 | VC 1. 0
Blank card ending branches. IBR, NTOT = 2 3 |BLANK BRANCH
Comment card. KOMPAR > 0. |C < n 1>< n 2>< Tclose ><Top/Tde >< Ie ><Vf/CLOP >< type >
Switch. 0.00E+00 1.00E+00 0.00E+00 0.00E+00 | FTE R 1. 1
Blank card ending switches. KSWTCH = 1. |BLANK SWITCH
Comment card. KOMPAR > 0. |C < n 1><>< Ampl. >< Freq. ><Phase/T0>< A1 >< T1 >< TSTART >< TSTOP >
Source. 1.00E+02 0.00E+00 0.00E+00 -1.00E+00 |11FTE 0 100. -1. 1.
Blank card ends electric sources. KCONST = 1 |BLANK SOURCE

List of input elements that are connected to each node. Only the physical connections of multi-phase lines are shown (capacitive
and inductive coupling are ignored). Repeated entries indicate parallel connections. Switches are included, although sources
(including rotating machinery) are omitted -- except that U.M. usage produces extra, internally-defined nodes "UMXXXX".
--------------+------------------------------
From bus name | Names of all adjacent busses.
--------------+------------------------------
R |VC *FTE *
VC |TERRA *R *
FTE |R *
TERRA |VC *
--------------+------------------------------
Card of names for time-step loop output. | VC FTE
Blank card ending requests for output variables. |BLANK OUTPUT

Column headings for the 3 EMTP output variables follow. These are divided among the 5 possible classes as follows ....
First 2 output variables are electric-network voltage differences (upper voltage minus lower voltage);
Next 1 output variables are branch currents (flowing from the upper node to the lower node);
Step Time VC FTE FTE
R
*** Switch "FTE " to "R " closed after 0.00000000E+00 sec.
0 0.0 0.0 0.0 0.0
250 .0025 8.94555411 100. .037379781
500 .005 19.5264118 100. .070603312
750 .0075 31.1544467 100. .09638356
1000 .01 43.2829353 100. .115007267
1250 .0125 55.4324609 100. .126954719
1500 .015 67.1964968 100. .132854097
1750 .0175 78.2437117 100. .133437921
2000 .02 88.3174269 100. .129502721
2250 .0225 97.232677 100. .121872771
2500 .025 104.871343 100. .111368468
2750 .0275 111.175812 100. .098779696
3000 .03 116.141601 100. .084844289
3250 .0325 119.809355 100. .070231527
3500 .035 122.256566 100. .055530444
3750 .0375 123.589346 100. .041242593
4000 .04 123.934514 100. .027778824
4250 .0425 123.432206 100. .015459568
4500 .045 122.229177 100. .00451808
4750 .0475 120.472921 100. -.00489392
% % % % % % Final time step, PLTFIL dumps plot data to ".PL4" disk file.
Done dumping plot points to C-like disk file.
5000 .05 118.306651 100. -.01269878

Extrema of output variables follow. Order and column positioning are the same as for the preceding time-step loop output.
Variable maxima : 196.016314 100. .253034843
Times of maxima : .00122 .1E-4 .61E-3
Variable minima : 0.0 0.0 -.24296238
Times of minima : 0.0 0.0 .00183

Blank card terminating all plot cards. |BLANK PLOT

--------------------------------------------------+--------------------------------------------------------------------------------

El anterior archivo de resultados está constituido de las siguientes partes:

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 19

• Un cuadro general que tiene la forma del Listado. En este listado se divide
en dos partes; la parte derecha corresponde a una copia de las líneas, tal
como fueron editadas en el archivo de datos de simulación; la parte
izquierda corresponde a la interpretación que el programa le da a cada
línea. Esta información es muy valiosa cuando se está buscando un error
en un archivo de datos. Esta forma del cuadro se repite en forma
intercalada en diferentes partes del archivo, hasta completar todas las
líneas del caso bajo simulación.
• Después de la línea "BLANK SOURCE" aparece la tabla de conexionado
de la red. De esta tabla de conexionado se excluye los acoples capacitivos
e inductivos y las fuentes. Los suiches si aparecen en la tabla de
conexionado de la red.
• Después de la línea "BLANK OUTPUT", aparece una tabla con los valores
de las variables eléctricas en el tiempo. El ordenamiento de las columnas
en esta tabla de datos es la siguiente:
- Número de paso de integración.
- Tiempo en segundos.
- Diferencias de potencial que fueron solicitados en las diferentes ramas.
- Voltajes nodales especificados en lista.
- Corrientes de rama.

• Valores de máximos y de mínimos, lo mismo que los tiempos en que se

presentan estos valores, para cada una de las variables y en el mismo
orden de la tabla a que se refiere el ítem anterior.

Algunos resultados que se han obtenido con el programa graficador “PLOTXY"

se observan en la Figuras 3.8 y 3.9.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 20

0.3
[A]

0.2 Corriente del circuito RLC serie

0.1

0.0

-0.1

-0.2

-0.3
0 10 20 30 40 [ms] 50
(file RLC.pl4; x-var t) c:FTE -R

Figura 3.8 Corriente en el circuito RLC serie

200

Voltaje en capacitor
150 Voltaje en la fuente

100

50
Corriente del circuito

-50
0 10 20 30 40 [ms] 50
(file RLC.pl4; x-var t) v:VC v:FTE c:FTE -R
factors: 1 1 1 200 200
offsets: 0 0 0 0 0

Figura 3.9 Voltajes y corriente en el circuito RLC serie

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 21

3.4 ANÁLISIS DE RESULTADOS PARA EL CIRCUITO RLC

• Calcule para el circuito simulado impedancia característica, frecuencia

natural de oscilación y factor de amortiguamiento que aparece asociado
con la corriente transitoria.
• Calcule la relación entre la amplitud del voltaje de la fuente y el primer pico
de corriente. Compare esta relación con la impedancia característica del
circuito.
• Utilizando el concepto de impedancia característica explique porqué el
voltaje máximo sobre el condensador se aproxima al doble del valor de la
fuente para circuito RLC serie con poco amortiguamiento (R<<Zc). Asumir
R de 1 Ω.
• Determine la variación del primer pico de corriente con respecto a la
resistencia del circuito, para 1≤ Ω R ≤ 50 Ω. ¿Cuál será la razón de esta
variación?.
• Deducir la expresión matemática para el voltaje sobre el capacitor.
• Cambiar la fuente DC por una de tipo sinusoidal de frecuencia industrial de
60 Hz, y una amplitud de 100 Vpico. Para este caso se desea observar el
voltaje sobre el condensador y compararlo con el obtenido para la
simulación con excitación escalón.
• Determinar la frecuencia de resonancia serie del circuito utilizando la
opción Frequency Scan. ¿Qué sucede si la fuente sinusoidal tiene esta
frecuencia de resonancia?.
Calcule los parámetros RLC de secuencia positiva para una línea de 200
km de longitud, con un conductor de fase de 30 mm de diámetro, una
resistencia AC de 0.06 Ω/km, una separación entre fases es de 7.5 m,
una disposición de conductores de fase en delta. Realizar la simulación
para una energización con una fuente AC de 230 kV
( Amp = 230000 × 2
3 ). Observar el voltaje al final de la línea considerando
un circuito RLC serie, un circuito PI y una cascada de 20 elementos PI
(Modelo de elementos semidistribuidos). De acuerdo al resultado
observado en la Figura 3.10. ¿cuál resulta ser un modelo adecuado para
modelar la línea en un estudio de energización?. Para esta simulación se
recomienda los siguientes parámetros:

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 22

Delta T = 1e-6 s, Tmax = 0.002 s.

500
[kV]

400

300
Fuente

200

100

0
0.0 0.4 0.8 1.2 1.6 [ms] 2.0
(file Noname.pl4; x-var t) v:FTE v:VCRLC v:VCPI v:VC20 v:VC20

Figura 3.10 Resultado esperado de la simulación con elementos concentrados

y semidistribuidos

• Para el circuito anterior y utilizando modelo de parámetros semidistribuidos

obtener la frecuencia de resonancia de la línea por cualquier método. ¿Qué
significado físico tiene esta frecuencia?.

• En este ítem se quiere introducir el concepto de propagación. Para la

mismo circuito anterior observar los voltajes cada 50 km para un tiempo de
simulación de 0.8 ms.

250

[V]

200

150

100

50

0
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 [ms] 0.8
(file Noname.pl4; x-var t) v:FTE v:V50 v:V100 v:V150 v:V200 v:V200

Figura 3.11 Resultado esperado para el efecto de propagación

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 23

Dar respuestas a las siguientes preguntas de acuerdo al resultado

esperado de la Figura 3.11:

- ¿Qué retrasos en el tiempo tienen los diferentes voltajes en cada

punto de la línea con respecto a la fuente?
- ¿Observar el voltaje a mitad de la línea y compararlo con el voltaje al
final de la misma?
- Observar el orden en que los voltajes van al 100 % con respecto a la
fuente y el orden en que los voltajes van al 200% con respecto a la
fuente. Para este caso tomar un Tmax de 1.5 ms. Explicar el
fenómeno.

• Para la mismo circuito anterior colocar una resistencia de carga al final de

la línea, de manera que se obtengan los resultados observados en la
Figura 3.12. Identificar la relación de esta resistencia con las características
eléctricas del circuito. ¿Cuál será la razón física de que no se presente
sobrevoltaje al final de la línea?

150
[V]

120

90

60

30

0
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 [ms] 1.6
(file Noname.pl4; x-var t) v:FTE v:V50 v:V100 v:V150 v:V200 v:V200

Figura 3.12 Resultado esperado para el efecto de propagación

• Se tiene un equivalente de cortocircuito en una barra de 220 kV

Voltaje de prefalla 226.23 kV ∠0.0°

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA

 24

Corriente de falla trifásica: 11872 A

Corriente de falla monofásica 14771 A
( X / R ) 3φ = 6.05
( X / R )1φ = 6.44
Las bases son 100 MVA y 220 kV

( X / R ) 3φ = Relación X/R que ve la corriente de corto trifásica

( X / R )1φ = Relación X/R que ve la corriente de corto monofásica

Realizar la simulación para determinar cuando ocurra una falla trifásica y

una falla monofásica en esta barra el valor de la corriente dinámica.
Recordar que la corriente dinámica es el pico máximo de la corriente que
se puede presentar durante una falla y depende del momento en que se
presente la falla y de la relación X/R. Investigar la relación que propone
la norma IEC-909 para calcular el valor de esta corriente a partir de la
relación X/R.

La teoría correspondiente al equivalente de Thevenin trifásico se puede

consultar en el Anexo 1.

LEONARDO CARDONA C. UNIVERSIDAD NACIONAL DE COLOMBIA