PROGRAMA ESPECIAL DE TRABAJO ACADÉMICO

GUÍA PRÁCTICA DEL PROGRAMA DE SIMULACIÓN EMTP/ATP

Lina Marcela Ramírez Arbelaez

C.C. 39 192 891

Carlos Alberto Serna Machado

C.C. 15 518 343

Profesora:

Clara Rosa Rojo Cevallos

UNIVERSIDAD NACIONAL DE COLOMBIA

FACULTAD DE MINAS
SEDE MEDELLÍN

2007
PROGRAMA ESPECIAL DE TRABAJO ACADÉMICO
GUÍA PRACTICA DEL PROGRAMA DE SIMULACIÓN EMTP/ATP ii

1. INTRODUCCIÓN..........................................................................................
2. OBJETIVOS................................................................................................
2.1 OBJETIVO GENERAL
2.2 OBJETIVOS ESPECÍFICOS
3. INTRODUCCIÓN AL MANEJO DEL PROGRAMA EMTP/ATP..................................
3.1 HISTORIA
3.2 ESTUDIOS A REALIZAR CON EL PROGRAMA EMTP/ATP
3.3 Modelos
3.4 ESTRUCTURA GENERAL DE UN ARCHIVO PARA LA ENTRADA DE DATOS AL ATP
3.5 REGLAS GENERALES PARA LA GENERACIÓN DEL ARCHIVO DE ENTRADA DE
DATOS
4. ESTRUCTURA GENERAL DEL ARCHIVO DE DATOS PARA UN CASO SIMPLE DE
SIMULACIÓN TRANSITORIA O DE ESTADO ESTACIONARIO DE UNA RED ELÉCTRICA
CON RAMAS SIMPLES RLC.....................................................................................
5 MODELAMIENTO DE UN CIRCUITO RLC SERIE................................................
6. MODELAMIENTO DE ONDAS VIAJERAS..........................................................
7. ENERGIZACIÓN DE LINEA MONOFASICA SIN PÉRDIDAS CON RESISTENCIA DE
PREINSERCION.....................................................................................................
8. OBTENCIÓN DE MODELOS LEE Y CLARK CON LINE CONSTANTS.
9. Obtener el modelo de onda viajera trifásico para línea transpuesta..
10. transformación interna que utiliza el programa ATP en comparación con la
transformación de componentes simétricas.
11. Determinar la sensibilidad de los parámetros: Resistencia de pérdidas, impedancia
característica y Velocidad de Propagación, con respecto a variación de la frecuencia (10
H ≤ f ≤ 5000 Hz) y la resistividad del terreno (100 – 3000 Ωm).
12. Determinar la relación entre Zc y R/4 (con R la resistencia total de pérdidas en la
ínea).
13 MODELAMIENTO DE UN RAYO
13.1PROCEDIMIENTO
13.1.1 Condensador....................................................................................
13.1.2 Inductancia......................................................................................
13.1.3 Resistencia de prueba.......................................................................
13.1.4 Resistencia del canal.........................................................................
13.1.5 Suiche.............................................................................................
13.1.6 ATP Settings....................................................................................
13.2Simulación
Rcanal (Ω)
14. MODELAMIENTO DE LÍNEAS.........................................................................
14.1PROCEDIMIENTO
PROGRAMA ESPECIAL DE TRABAJO ACADÉMICO
GUÍA PRACTICA DEL PROGRAMA DE SIMULACIÓN EMTP/ATP iii

14.2SIMULACIÓN
15. BIBLIOGRAFIA

LISTA DE FIGURAS
Figura 1 circuito RLC ...........................................................................................
Figura 2 Parametros del elemento RLC...................................................................
Figura 3 Fuente DC..............................................................................................
Figura 4 parametros de la fuente DC.......................................................................
Figura 5 Switche..................................................................................................
Figura 6 condensador...........................................................................................
Figura 7 Parametros del condensador.....................................................................
Figura 8 Aterrizamiento del elemento.....................................................................
Figura 9 Pantalla indicadora de medicion................................................................
Figura 10 Forma de correre el programa.................................................................
Figura 11 Forma de graficar las senales de interes....................................................
Figura 12 voltaje en el capacitor............................................................................
Figura 13 Coriente en circuito serie........................................................................
Figura 14 Coriente en circuito serie con R=1...........................................................
Figura 15 Coriente en circuito serie con R=10..........................................................
Figura 16 Coriente en circuito serie con R=20.........................................................
Figura 17 Coriente en circuito serie con R=40.........................................................
Figura 18 Coriente en circuito serie con R=50.........................................................
Figura 19 Circuito RLC para la deduccion del voltaje en el capacitor...........................
Figura 20 Circuito con alimenacion senoidal

Figura 21 selección de fuente de ac........................................................................

Figura 22 Parámetros fuente senoidal.....................................................................
Figura 23 voltaje en el capacitor............................................................................
Figura 24 Voltaje en el capacitor con señal aplicada tipo escalon.................................
Figura 25 Opciones de frequency scan....................................................................
Figura 26 Frecuencia de resonancia del circuito .....................................................
Figura 27 Separación de los conductores................................................................
Figura 28 Montaje del circuito RLC.........................................................................
Figura 29 El voltaje al final del circuito...................................................................
Figura 30 Circuito PI............................................................................................
PROGRAMA ESPECIAL DE TRABAJO ACADÉMICO
GUÍA PRACTICA DEL PROGRAMA DE SIMULACIÓN EMTP/ATP iv

Figura 31 El voltaje en el capacitor 1 del circuito PI...................................................

Figura 32 El voltaje en el capacitor 2 del circuito PI..................................................
Figura 33 Circuito PI parâmetros distribuídos...........................................................
Figura 34 Onda viajera en un circuito con parámetros distribuidos.............................
Figura 35 Frecuencia de resonancia para el modelo de parámetros distribuidos............
Figura 36 Voltajes en todo el circuito......................................................................
Figura 37 voltajes en todo el circuito.....................................................................
Figura 38 falla trifásica........................................................................................
Figura 39 corriente de falla trifasica.......................................................................
Figura 40 falla monofásica.....................................................................................

Figura 41 corriente de falla monofásica

Figura 42 Montaje onda viajera.............................................................................

Figura 43 Selección de fuente tipo heidler...............................................................
Figura 44 fuente heidler.........................................................................................
Figura 45 ...........................................................................................................
Figura 46 ...........................................................................................................
Figura 47 ...........................................................................................................
Figura 48 ...........................................................................................................
Figura 49 ...........................................................................................................
Figura 50 resistencia............................................................................................
Figura 51 voltajes en el recorrido del circuito...........................................................
Figura 52 voltajes en el recorrido del circuito...........................................................
Figura 53 voltajes en el recorrido del circuito...........................................................
Figura 54 voltajes en el recorrido del circuito...........................................................
Figura 55 voltajes en el recorrido del circuito.
Figura 56 ...........................................................................................................
Figura 57 ...........................................................................................................
Figura 58 Extraída de la guía del curso por Leonardo Cardona..................................
Figura 59 ...........................................................................................................
Figura 60

Figura 61 corriente de falla monofásica

…………………………………………………………..

Figura 62 ............................................................................................................
PROGRAMA ESPECIAL DE TRABAJO ACADÉMICO
GUÍA PRACTICA DEL PROGRAMA DE SIMULACIÓN EMTP/ATP v

Figura 63 ...........................................................................................................
Figura 64.............................................................................................................
Figura 65 ...........................................................................................................
Figura 66 disposición física Torre de 230 KV para transmisión .................................
Figura 67 ...........................................................................................................
Figura 68 ...........................................................................................................
Figura 69 ...........................................................................................................
Figura 70.............................................................................................................
Figura 71.............................................................................................................
Figura 72 ...........................................................................................................
Figura 73 ...........................................................................................................
Figura 74 ...........................................................................................................
Figura 75 …………………………………………………………
Figura 76 ...........................................................................................................
Figura 77 ...........................................................................................................
Figura 78 ...........................................................................................................
Figura 79 ...........................................................................................................
Figura 80

Figura 81 Representación de la nube y el canal de descarga

Figura 82 Circuito para la simulación de la corriente de rayo......................................

Figura 83 Datos del Capacitor en el ATP...................................................................
Figura 84 Datos de la inductancia en el ATP.............................................................
Figura 85 Datos de la resistencia de prueba en el ATP...............................................
Figura 86 Datos del suiche en el ATP .....................................................................

Figura 87 Parámetro de corrida del simulador ATP

Figura 88 Forma de onda de la corriente para la resistencia del canal de 10500

Figura 89 Forma de obtener fuente tipo 13………………………………………………

Figura 90 Datos de la Fuente tipo 13 para simular una Fuente tipo rayo....................

Figura 91 Forma de onda tipo rayo simulada con la fuente tipo 13

Figura 92 Estructura típica de 230 kV.

Figura 93 Modelación de la estructura metálica

PROGRAMA ESPECIAL DE TRABAJO ACADÉMICO
GUÍA PRACTICA DEL PROGRAMA DE SIMULACIÓN EMTP/ATP vi

Figura 94 Modelo de Clarke para tramo de red monofásica

Figura 95 Forma de entrar los datos al modelo Clarke.

…………………………………………………………

Figura 96 Modelo de [Link].

Figura 97 Parametros de entrada del modelo de [Link].

Figura 98 Datos de entrada del modelo de [Link].

Figura 99 Vista de los conductores en el modelo de [Link]......................................

Figura 100 Suiche controlado por voltaje (cadena de aisladores)

Figura 101 Datos de entrada de la cadena de aisladores

Figura 102 Datos de entrada de la Fuente tipo 13, simulando descarga atmosferica de
30kA............................................................................................................
Figura 103 Tacs....................................................................................................
Figura 104 Datos de entrada de los TACS................................................................
Figura 105 Sistema completo a simular...................................................................
Figura 106 Settings de la simulación de una descarga atmosférica en un cable de
guarda..........................................................................................................

Figura 107 Voltajes en las cadenas de aisladores

Figura 108 Voltaje de soporte de la cadena de aisladores.

 1. INTRODUCCIÓN

En esta guía se mostrará la forma de realizar las simulaciones en ATP, facilitando el

ejercicio profesional y/o docente.

La guía consta de ocho capítulos enfocados a la familiarización con el programa, de

modo tal que quien la estudie este en capacidad de resolver problemas de transitorios
electromagnéticos y electromecánicos en sistemas de potencia eléctricos polifásicos.
2. OBJETIVOS

1 2.1 OBJETIVO GENERAL

Construir, desarrollar y documentar casos de simulación de sistemas de potencia

eléctrica en los distintos niveles de tensión (1,2,3 y 4) utilizando el programa ATP.

2.2 OBJETIVOS ESPECÍFICOS

 Elaborar una guía didáctica de aplicación de esta herramienta de simulación

para estudiantes y docentes.
 Presentar los diferentes sistemas con los cálculos, características y datos
necesarios, como también la forma de entrada para la simulación.
 Elaborar una guía para la construcción de los sistemas en el simulador.
 Presentar los resultados relevantes de las simulaciones para el sistema.
3. INTRODUCCIÓN AL MANEJO DEL PROGRAMA EMTP/ATP

El programa ATP (Electromagnetic Transients Program), es un programa para

computador digital, utilizado para la simulación de transitorios electromagnéticos,
electromecánicos y funciones de control en sistemas de potencia polifásicos.

3.1 HISTORIA

El ATP fue desarrollado en su versión inicial por el Dr. Hermann Dommel a finales de la
década del 60 en Alemania (inicialmente con el nombre de EMTP), quién
posteriormente vendería los derechos 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.

Del programa existen versiones para diferentes tipos de computadores. La versión

para micros es la más popular y ya poco se habla de las versiones que hubo para otro
tipo de computadores.

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 permiten un trabajo más
amigable para muchas de las aplicaciones. La interfaz que se utilizará en esta guía 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 el sistema operativo DOS hasta las versiones actuales para el
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. La interfaz ATPDRAW
normalmente se puede actualizar vía Internet; algunas direcciones de donde se puede
bajar son:

[Link]
[Link]
 3.2 ESTUDIO QUE SE PUEDEN REALIZAR CON El ATP

Entre los difeferentes estudios que se pueden realizar están los siguientes:
 Cálculo de parámetros en líneas aéreas y cables subterráneos
 Cálculo en estado estacionario sobre un sistema de potencia por maniobra
(suicheo)
 Sobrevoltajes por descargas atmosféricas
 Cálculos de cortocircuito que involucran 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

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).

3.3 MODELOS

Para la simulación del sistema de potencia, el programa posee varios modelos que dan
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 acció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.
 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).
 Máquinas eléctricas
 Simulación del control de las máquinas y del sistema de potencia

2 3.4 ESTRUCTURA GENERAL DE UN ARCHIVO PARA LA ENTRADA DE

DATOS AL ATP

Este archivo está ceñido por lo tanto a las reglas de ese lenguaje de programación
(Fortran), esto hace necesario conocer en detalle como está organizado el archivo de
entrada.

Existen dos formas generales de utilizar el programa:

 Para hacer una simulación de un sistema de potencia (o una red circuital) con el
fin de obtener la respuesta transitoria. el sistema bajo simulación puede
contener tanto elementos circuitales, como elementos de control que
interactúan con la red. Los sistemas de control Se pueden simular de manera
independiente
 Obtención de parámetros correspondientes a elementos del sistema de
potencia, como líneas aéreas, cables subterráneos, transformadores,
pararrayos, etc.
Las posibilidades de uso común del ATP como programa de simulación o de cálculo de
parámetros son las siguientes:
 Caso simple de simulación de una red eléctrica sin bloques de control.
 Caso simple de simulación de un sistema de control sin incluir la red eléctrica.
 Caso simple de simulación de red eléctrica y sistema de control en forma
conjunta.
 Obtención de parámetros y modelos de líneas de transmisión.

3.5 Reglas generales para la generación del archivo de entrada de datos

El archivo de entrada de datos debe ser válido para lenguaje Fortran. Las reglas
generales son las siguientes:

 El archivo debe tener un nombre y se recomienda que tenga una extensión *.ATP
que es la que por defecto asigna la interfaz ATPDRAW, y para darle nombre al
archivo de resultados se recomienda igualmente el mismo nombre, pero con
extensión *.RES o *.LIS
 El archivo de entrada de datos tiene 80 columnas.
 El archivo de entrada puede tener líneas de comentarios, las cuales empiezan con
el carácter “C” en la primera columna, seguido de un espacio en blanco y luego el
comentario que se desee en la misma línea.
 Los valores reales con formato científico (5.0E-2 por ejemplo) y los números
enteros, deben estar ajustados a la derecha en el campo. Los valores reales con
formato tipo “F” pueden ir en cualquier lugar dentro del campo, siempre y cuando
se les coloque el punto. La mayoría de los datos en el ATP de tipo real se piden con
formato tipo “E”. Para este tipo de formato el dato se puede colocar como tal o
como si fuera en formato tipo “F”, pero colocando el punto. Ejemplo, si un valor de
entrada es 10E-02, este dato se puede colocar como 0.1 con la ventaja de que no
tiene que estar ajustado a la derecha.
 No se deben dejar líneas en blanco arbitrariamente, porque éstas las utiliza el
programa para identificar cuando ha terminado un grupo de datos.
4 Estructura general del archivo de datos para un caso simple de simulación
transitoria o de estado estacionario de una red eléctrica con ramas simples
RLC

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.

Un archivo de entrada general al ATP se caracteriza por los siguientes grupos de datos:
 Datos iniciales para identificar el comienzo 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 del tiempo de simulación, tiempo de paso de integración,
especificación de las unidades de los parámetros LC y control de salida de
información. También se especifica en este bloque de datos si es un caso de
simulación de estado estacionario.
 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.
Para una simulación en particular no deben aparecer necesariamente todos los bloques
señalados.

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, sólo se emplea un tipo de rama (rama RLC
serie de parámetros concentrados) para la modelación de la red, el interruptor que se
usa es controlado por tiempo en forma determinística y se utilizan fuentes sencillas
como escalón y cosenoidal.

APLICACIÓNES

5 Modelamiento de un circuito RLC serie

Objetivo de la práctica

Familiarizar el programa EMTP/ATP mediante el análisis de circuitos simples RLC, RL

en estado transitorio.

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

transmisión con parámetros distribuidos.

En la figura 1 que se muestra a continuación se observa un circuito RLC serie con los
estados transitorios de voltaje y corriente en el capacitor.

VC
FTE
U

Figura 1 Circuito RLC

Los parámetros del circuito son los siguientes

R = 50 Ω
L = 150mH
C = 1uF
 Figura 1 Parámetros del el elemento RLC

Luego para sacar la fuente DC propuesta en el montaje de la figura No 1 haremos lo

siguiente, hacemos clic derecho y aparcera el pantallazo mostrado en la figura 3.

Figura 2. Fuente DC
Luego de obtener en nuestras pantallas la fuente DC hacemos doble click y
parametrizamos nuestra fuente:

Figura 4. Parámetros de la fuente DC

Ahora miremos el switche, este elemento nos define los tiempos de cierre y de
apertura. En la figura 5 se mostrara la forma de obtener este elemento

Figura 5. Switche
La figura 5 nos muestra las opciones de medida sobre este elemento, como es notorio
se observa medida de corriente, voltaje, corriente y voltaje y energía.
Por ultimo tenemos el condensador, de la misma manera como los anteriores casos,
click derecho y desplegamos la ventana siguiente

Figura 6 Condensador

Figura 7. Parámetros del Condensador

Luego aterrizamos el condensador dando doble click izquierdo sobre el y chuleamos
donde dice tierra como nos muestra la figura siguiente

Figura.8 Aterrizamiento del elemento

Lo más común es llamar por un nombre (tiene que ser en mayúscula) los nodos que
necesitemos medir, en este caso es el nodo del capacitor, lo hacemos dando doble
click derecho en el nodo y aparecerá la siguiente ventana (Figura 9).

Figura [Link] indicadora de medición.

Para la simulación es necesario cerciorarse de que el montaje esté en las condiciones

requeridas, que no presente ningún tipo de error, como sobre poner líneas, tierras etc.

Después del montaje del circuito es necesario correr el programa mediante la tecla F2
ó como se muestra la siguiente figura (10)
 Figura [Link] de correr el programa.

Luego de correr el programa para obtener las señales se hace lo siguiente:

Figura 11. Forma de graficar las señales de interés.

Las señales voltaje y corriente se pueden observar en las figuras 12 y 13

respectivamente
200
[V] 0.3
[A]
160
0.2

120
0.1

80
0.0

40
-0.1

0
0.00 0.02 0.04 0.06 0.08 [s] 0.10 -0.2
0.00 0.02 0.04 0.06 0.08 [s] 0.10
(file trabajo.pl4; x-var t) v:VC
(file trabajo.pl4; x-var t) c:FTE -XX0013

Figura 12 Voltaje en el capacitor

Figura 13 Corriente en circuito serie

5.1.1 Impedancia característica, frecuencia de oscilación y factor de

amortiguamiento asociado con la corriente transitoria.

 Calcular para el circuito simulado la impedancia característica, frecuencia natural

de oscilación y el factor de amortiguamiento que aparece asociado con la corriente
transitoria.

Impedancia característica

L
Zc 
C ( 1)

Donde L = 150mH y C = 1uF

Reemplazando en 1 los valores se obtiene:

150 E  3
Zc  Por lo tanto Z c  387.29833
1E  6
Frecuencia natural de oscilación
1
f0  Hz ( 2)
2 LC
Al remplazar valores se obtiene:
1
f0  Hz = 410.93Hz
2 150 E  9

Factor de amortiguamiento
Zc
 ( 3)
R
Donde:
387.29833
 , por lo tanto   7.7459
50

5.1.2 Relación entre la amplitud del voltaje del a fuente y el primer pico de
corriente.

 Calcular la relación entre la amplitud del voltaje de la fuente y le primer pico de

corriente. Comprar esta relación con la impedancia característica del circuito.

Tabla 1. Valores específicos

En la anterior tabla es notorio que el valor del primer pico de corriente es
aproximadamente de 0.231 A y la amplitud de la fuente de voltaje es de 100 la
relación entre estos dos valores es de V/I = 432.9
Los valores de la impedancia característica calculada anterior mente son aproximados
a los calculados por medio de la relación V/I siendo esta última mayor, teniendo en
cuenta:
 La impedancia característica es la impedancia vista por la fuente sin tener
en cuenta el efecto resistivo.
 En la relación V/I se tiene en cuenta el efecto resistivo.
Si embargo los valores obtenidos son aproximados y se pueden considerar como
iguales.

5.1.3 Cuando el circuito RLC presenta poco amortiguamiento (R<<<Zc), el

valor de voltaje máximo en el condensador, se aproxima al doble del valor de
la fuente.

 Utilizar el concepto de impedancia característica explicar porqué el voltaje

máximo sobre el condensador se aproxima al doble del valor de la fuente para
un circuito RLC con poco amortiguamiento (R << ZC). Asumir R de 1 Ω.

El voltaje a lo largo de una línea está formado por dos componentes, incidente y
reflejado significa esto que en cualquier punto de la línea el voltaje es la suma de dos
componentes cuando en una línea no hay perdidas (R = 0), para esto tenemos la
expresión de voltaje en el condensador, suponiendo R = 0 :

VLC  1 1 
VC   cos t   V (4)
ZCC  LC LC 


y reemplazando Zc obtenemos

1
VC  V cos t  V ( 5)
LC
Se deriva V con respecto al tiempo y se iguala a cero para hallar el máximo
 dVc 1 1
V sen t  0 ( 6)
dt LC LC
1
sen t 0 (7)
LC
1
t  2 n n  0,1,2,3,...
LC
t  LC 2 n 8)

Reemplazando para estos tiempos en (5) obtenemos

Vc  V cos 2 n  V
cos 2 n  1
Vc  2Vc

La impedancia característica aparece cuando en una red hay fenómenos transitorios y

esta permite despreciar los efectos resistivos de la red simulada por un circuito RLC, y
como la resistencia se supone cero por lo dicho anteriormente lo que pasa que el flujo
de la corriente no encuentra oposición en entonces el voltaje en el condensador
alcaza a ser el “máximo” y se considera ideal mente el doble de la fuente esto ocurre
por el fenómeno de reflexión.

5.1.4 Variación del primer pico de corriente con respecto a la variación del
valor de resistencia en el circuito (1Ω<R<50Ω).

Determinar la variación del primer pico de corriente con respecto a la resistencia

del circuito para 1Ω < R < 50Ω. ¿Cuál es la razón de esta variación?

 Para R = 1 Ω. , se obtiene la figura 14

 0.3
[A]
0.2

0.1

0.0

-0.1

-0.2

-0.3
0.00 0.02 0.04 0.06 0.08 [s] 0.10
(file trabajo.pl4; x-var t) c:FTE -XX0013

Ilustración 14 Corriente en circuito serie para R= 1 Ω.

Tabla 2. valores obtenidos

Como se observa el valor pico de la corriente es de 0.25A aproximadamente

 Para R = 10 Ω. Se realiza el mismo procedimiento

0.3
[A]

0.2

0.1

0.0

-0.1

-0.2

-0.3
0.00 0.02 0.04 0.06 0.08 [s] 0.10
(file trabajo.pl4; x-var t) c:FTE -XX0013

Figura 15. Corriente en circuito serie para R= 10 Ω

 Tabla 3. valores obtenidos

Para este caso se observa que el pico máximo de la corriente es de 0.24949 A

 Para R = 20 Ω. Se realiza el mismo procedimiento

0.3
[A]

0.2

0.1

0.0

-0.1

-0.2

-0.3
0.00 0.02 0.04 0.06 0.08 [s] 0.10
(file trabajo.pl4; x-var t) c:FTE -XX0013

Figura 16 Corriente en circuito serie para R= 20 Ω.

Tabla 4. valores obtenidos

En este caso el valor pico es de 0.244 aprx

 Para R = 40 Ω.
 0.3
[A]
0.2

0.1

0.0

-0.1

-0.2

-0.3
0.00 0.02 0.04 0.06 0.08 [s ] 0.10
(file trabajo.pl4; x-var t) c:FTE -XX0013

Figura 17. Corriente en circuito serie para R= 40 Ω.

Tabla 5. Valores obtenidos

En este caso se observa que el valor pico es de 0.23548

 Para R =50 Ω.

0.3
[A]

0.2

0.1

0.0

-0.1

-0.2
0.00 0.02 0.04 0.06 0.08 [s] 0.10
(file trabajo.pl4; x-var t) c:FTE -XX0013

Figura 18. Corriente en circuito serie para R= 50Ω.

 Tabla 6 Valores obtenidos

En este caso se observa que el valor pico es de 0.2314

Es notorio observar que la diferencia del pico máximo entre el primer caso y el último
caso es muy pequeña no es muy significante, lo que quiere decir esto que la línea no
depende de la resistencia si no de la resistencia característica, esto ocurre en los
primeros ciclos. Esto permite verificar que para una línea infinita no se tiene en cuenta
los efectos resistivos.

5.1.5. Expresión matemática del voltaje del condensador

Figura 19. Circuito RLC para la deducción del voltaje en el condensador

di 1
V (t )  Ri (t )  L 
dt C  i(t )dt 9)

Con: V (t )  U  Cte

Derivando con respecto al tiempo.

dv di d 2i 1
 R  L 2  i (t ) 10)
dt dt dt C
 d 2 i R di 1
2
  i (t )  0 11)
dt L dt LC

Ecuación característica.

R 1
2    0 12)
L LC

R 1 R2
1, 2     2 13)
2L LC 4 L

R 1  R 2C 
1, 2    1   14)
2L LC  4 L2 

De acuerdo a la naturaleza de estas raíces obtenemos la forma de respuesta del

circuito, para el circuito de segundo orden, como el planteado estas raíces son
complejas, por lo tanto la respuesta tiene una atenuación exponencial.

La solución de la ecuación es de la forma:

     
 c cos 1 1  R C  t  c sen 1 1  R C  t 
R 2 2
 t
i (t )  e 2L
15)
 1  LC  4 L  
2
 LC  4 L2   
    

Aplicamos las condiciones iniciales para la solución anterior obtenemos las constantes
C1 y C2:

Con i (0)  0 , C1=0 entonces:

   R 2C  
   
R  1  4 L2  
 t  
i (t )  e 2L
c 2 sen t  16)
 LC 
 
 

di (t ) V
Con t = 0 
dt L

Reemplazando C2 obtenemos:

 R    
 R  2 L t   1  R C     2 L t 1  R C    1  R 2 C    
2 R 2
di(t )
 c2   e sen 1 t e 1 cos  1  t 17)
dt  2L   LC  4 L2    LC  4 L2    LC  4 L2    
        

V  1  R 2C  
 c2  1   18)
L  LC  4 L2  
 
V
c2 
 1  R 2C   19)
L 1  
 LC  4 L2  
 
V
c2 
 L R 2C  
 1  2  
 C 4L  
 
V
c2 
  R 2C  
Zc  1  
  4 L  


La corriente que circula en el circuito es de la forma

  2 
 1 R 
V R  4 L 2

i (t )  e 2 L sen t  20)
 R 2C   LC 
Zc 1  
2   
 4L   

Zc
Ahora, si  , entonces:
R

 
 1 1 
V 4
sen t  21)
R
i (t )  e 2L

 1  LC
Zc 1    
 4   

Ya obtenido la ecuación para la corriente, procedemos ha hallar la ecuación para el

voltaje en el condensador de la siguiente forma:

1
VC 
C  i (t ) dt

E integrando la corriente a través de ayuda de tablas de integrales tenemos la

siguiente ecuación solución:

 1 1 1 
R  t 1 1 1 
V e 2L
 R   
VC   sen 4 t 4 cos 4 t
1  1   4L LC LC LC 
Z C C 1 1   22)
4   R   
2
4   
 4 L2   LC 
   
 

5.1.6 Cambiar la fuente del circuito por una senoidal con frecuencia de 60
Hz.

 Cambiar la fuente de DC por una de tipo sinusoidal de frecuencia industrial 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.

El montaje utilizado para esta simulación es el siguiente:

AMP = 100Vpico, Tsta = -1, Tsto = 1 , f=60 Hz

FUENTE VC
U

Figura 20 Circuito con alimentación senosoidal

Para obtener la fuente sinusoidal haremos lo siguiente

Figura 21 .selección de fuente de ac.

 Figura [Link]ámetros fuente senoidal.

200
[V]
150

100

50

-50

-100

-150

-200
0.00 0.02 0.04 0.06 0.08 [s] 0.10
(f ile linea.pl4; x-var t) v:VC

Figura [Link] en el capacitor

 900
[V]

560

220

-120

-460

-800
0.00 0.02 0.04 0.06 0.08 [s] 0.10
(file linea.pl4; x-var t) v:VC

Figura [Link] en el capacitor con señal aplicada tipo escalon

En los dos casos anteriores se nota que existe un fenómeno transitorio en los primeros
ciclos a demás esta se va atenuando

En comportamiento transitorio, la onda tipo escalón producida por la fuente DC se

amortigua según su componente DC hasta alcanzar el valor constante de 100 V,
siguiendo a la fuente; mientras que en la respuesta a la onda senoidal el capacitor
oscile de forma estable ante el sistema en 100 V.

5.2 Frecuencia de resonancia utilizando Frequency scan del ATP.

 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?.

Mediante el uso de la opción de frequency scan en el ATP, la cual se encuentra en los

settings de la simulación, obtenemos la siguiente gráfica:
 Figura 25. Opciones de frequency scan

Figura [Link] de resonancia del circuito

Lo anterior muestra que la frecuencia de resonancia del circuito es 410. Si la fuente

tuviera la frecuencia de resonancia del circuito este fenómeno causaría una oscilación
crítica y nuestro sistema, la línea, y todos los elementos de esta, comenzarían a
dañarse, cuando los elementos de la línea entran en resonancia (L=C) en un sistema,
habrá un intercambio de energía, causando esto perjuicios para el sistema de potencia
por las fallas ocasionadas en la construcción y en el aislamiento en los mismos
elementos que conforman el circuito.

5.2. Cálculo de los parámetros de una línea

Calcular los parámetros RLC 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 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. 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 de la figura 3.10 ¿Qué resulta ser un modelo adecuado para modelar la línea
en un estudio de energización?
Para la simulación se recomiendan los siguientes parámetros:

Figura [Link]ón de los conductores

Para encontrar los parámetros RLC de la línea se utilizan las siguientes fórmulas:
 DMG 
Z1  RAC  j Ln 23)
2 RMG km
 55.55 nF
C1 
DGM km 24)
Ln
RMG
Donde

DMG  3 Dab  Dac  Dbc = 7.5 m 25)

1

RMG  r  e 4 =0.015 m 26)

De lo anterior tenemos que:

  
Z1   0.06  j 0.48   200 km  12  j 97.486 
 km 

nF
C  8.593  200km  1.719 F
km
R  12 

5.2 Circuito serie RLC

Figura 28. Montaje del circuito RLC.

El voltaje al final del circuito se muestra en la siguiente grafica.

 400
[kV]
350

300

250

200

150

100

50

0
0.0 0.4 0.8 1.2 1.6 [ms] 2.0
(file linea.pl4; x-var t) v:VC

Figura 29 El voltaje al final del circuito

5.2.1 Circuito PI

Figura 30. Circuito PI

 400
[kV]
300

200

100

-100

-200

-300

-400
0.0 0.4 0.8 1.2 1.6 [ms] 2.0
El (file linea.pl4; x-var t) v:VC1

Figura 31 El voltaje en el capacitor 1 del circuito PI

En las dos gráficas anteriores se nota que las dos señales son muy similares solo
varían en magnitudes y a demás el tiempo para alcanzar el valor máximo de la
segunda señal es un poco mas rápido.

200
[kV]
150

100

50

-50

-100

-150

-200
0.0 0.4 0.8 1.2 1.6 [ms] 2.0
(file linea.pl4; x-var t) v:VC2

Figura 32. El voltaje en el capacitor 2 del circuito PI

5.2.2 Modelo PI con parámetros semidistribuidos

Para el modelo con parámetros semidistribuidos se debe tener en cuenta las partes en
que es dividida la línea de transmisión, para este caso los 200 km de línea los se
dividen en 20 partes y a cada elemento determinado (R, L, C) deben ser
multiplicado por 10 Km. :

R  0.06  10   0.6 

XL
X  0.487  10   4.8   L   12.732  mH 

C  8.593  10 nF  85.593 nF  0.0856 nF

Al final se obtiene que:

Los parámetros Obtenidos son :
R = 0.6 Ω
Xl = 4.87 Ω entonces L = 12.732 mH
C = 0.08595 Uf

Figura 33. Circuito PI parâmetros distribuídos

 500
[kV]

400

300

200

100

0
0.0 0.4 0.8 1.2 1.6 [ms] 2.0
(file Practica1.pl4; x-var t) v:VC

Figura 34 .Onda viajera en un circuito con parámetros distribuidos

En las líneas de transmisión el voltaje es una onda que viaja a través de esta a la
velocidad de la luz (300000 m/s), la mejor configuración para representar la línea
será aquella en la que el voltaje alcance el valor máximo de voltaje en la energización
en el menor tiempo, y de información de los transitorios que ocurren en el cierre de un
alinea de transmisión. Y se llego también a concluir que el modelo más adecuado para
la energización de la línea es el modelo con parámetros semidistribuidos.

5.2.3 Frecuencia de resonancia para el modelo de parámetros distribuidos.

Utilizando el Frequency Scan en el ATP:

 12
[MV]

10

0
0 200 400 600 800 [s] 1000
(file Practica1.pl4; x-var t) v:VC

Figura 35. Frecuencia de resonancia para el modelo de parámetros

distribuidos.

Tabla [Link] de la Frecuencia de resonancia

La frecuencia de resonancia es de 370 Hz. Si esta frecuencia se presenta en la red se
podrían ocasionar daños en los elementos y en su aislamiento debido a la elevación
de corriente y de voltaje

5.2.4 Para el circuito de parámetros distribuidos determinar voltaje en la

línea cada 50km.

 Para el circuito anterior observar los voltajes cada 50 Km. para un tiempo de
simulación de 0.8 ms En este ítem se quiere introducir el concepto de propagación..

Los voltajes obtenidos a 50 km, 100 km, 150 km y 200 km se presentan en la

siguiente figura 36:
 500
[kV]

400

300

200

100

0
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 [ms] 0.8
(file Practica1PICDA.pl4; x-var t) v:KMIV v:KMI v:KMII v:KMIII

Figura [Link] en todo el circuito

 Dar respuestas a las siguientes preguntas de acuerdo al resultado esperado en 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?.

Los tiempos tomados para cada punto de la línea se muestran en la siguiente tabla:

Longitud [km] Tiempo [ms]

50 0.1
100 0.235
150 0.41
200 0.54
Tabla 23. Tiempos en los voltajes de la línea
Observe que la onda de voltaje se va propagando a lo largo de la línea, y la magnitud
se va incrementando con la longitud, esto se debe al efecto Ferranti o efecto capacitivo
de la línea. El cual consiste cuando la línea esta en vacío esta ve impedancia infinita y
la onda se refleja y hasta se pude duplicar esta magnitud de la onda

 Observar el voltaje a mitad de la línea y compararlo con el voltaje al final de la

misma.

El voltaje a mitad de la línea (km = 100) y al final de la línea (km = 200) lo podemos
observar en la figura del numeral 9, como se mencionaba debido a la propagacion de
la onda de voltaje, se demora mas tiempo en incrementarse este al final de la línea,
además la magnitud del voltaje al final de la línea es mucho mayor, cuando la onda
llega al final de la línea sufre un fenómeno de reflexión, esta reflexión interfiere
constructivamente con la onda que llega y por ende aumenta su amplitud.

 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 Tmáx de 1.5 ms. Explicar el fenómeno.

El resultado de la simulación es el siguiente:

 250

[V]

200

150

100

50

0
0.0 0.3 0.6 0.9 1.2 [ms] 1.5
(file ModeloPropagación.pl4; x-var t) v:FTE v:50 v:100 v:150 v:200

Figura 37 voltajes en todo el circuito

Cuando la onda de voltaje llega al final de la línea y se refleja hacia la fuente, esta en
esta refracción ayuda constructivamente con las ondas que van de la fuente al final de
la línea a duplicar su amplitud.

5.3.5 Carga en el extremo de l a línea y determinación de los voltajes

cada 50km.

Para el mismo circuito anterior colocar una resistencia de carga al final de la línea, de
manera que se obtengan los resultados 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?

En una línea infinita no existe la reflexión ya que esta tiene un lado Terminal y evita la
propagación de la onda y a demás cuando se conecta una carga a la línea, esta se
comporta como un elemento infinito donde no hay opción para la reflexión de la onda
a y cuando se conecta una Zc = Zcaracterística el voltaje presenta la misma magnitud
durante toda la línea por lo tanto se comporta como una línea infinita.
 En conclusión el valor de la impedancia de carga Terminal ZT que se debe conectar
para evitar ondas reflejadas debe ser igual a la ZC de la línea de transmisión.

En la figura anterior observamos que el voltaje al final de la línea cae con respecto a
los demás esto se debe a que la impedancia de carga es menor que la impedancia
característica.

5.3 Calculo de corto circuito

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

Voltaje de prefalla 226.23 0.0º kV .
Corriente de falla trifásica 11872 A.
Corriente de falla monofásica 14771 A.a
La siguiente información es obtenida por el sistema de distribución de EEPPM, datos
dados por la Empresa.

X
   6.05
 R 3

X
   6.44
 R 1

Las bases son 100 MVA y 220 kV.

X
   Relación X/R que ve la corriente de corto trifásica. 27)
 R 3

X
   Relación X/R que ve la corriente de corto monofásica. 28)
 R 1

Las bases son 100 MVA y 220 kV.

X
   Relación X/R que ve la corriente de corto trifásica.
 R 3
X
   Relación X/R que ve la corriente de corto monofásica.
 R 1

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 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.

Calculamos las impedancias de secuencia.

VL
X1 
3  I CC 3
29)

3  VL 2  VL
X0  
I CC 1 3  I CC 3
30)
Donde
X 1  11 .00 
X 0  4.52 

5.3.1 Cortocircuito trifásico

 Ilustración 38 falla trifásica

La corriente de falla es la siguiente

20

[kA]

15

10

-5

-10

-15

-20
0 4 8 12 16 [ms] 20
(file ultimo.pl4; x-var t) c:X0007A-

Ilustración 39 corriente de falla trifasica

El valor de la corriente de corto tiene un valor aproximado de 17000 A. valor que

corresponde al corriente dinámica para la falla trifásica en la barra de 200kV.

5.3.2 Cortocircuito monofásico

 Probes &

Ilustración 40 falla monofásica

30
[kA]
20

10

-10

-20

-30
0 4 8 12 16 [ms] 20
(file ULTIMO.pl4; x-var t) c:X0001A-

Ilustración 41 corriente de falla monofásica

La anterior gráfica nos muestra la corriente durante la falla, cuyo valor pico es
aproximadamente 20900 A, que corresponden al valor de corriente dinámica para la
falla monofásica en la barra de 200 kV.
5.3.3 Aplicación de la norma IEC 600909

Esta norma presenta diferentes recomendaciones para el análisis de corto circuito, se

sugiere el cálculo de la corriente inicial simétrica trifásica

I 
´''
K y de la corriente pico

I 
´
P de la siguiente manera :

cU m
I K'' 
3Z 31) I P  x  2  I K'' 32)

Suponiendo que la falla fue balanceada (Falla Trifásica), el corto circuito fue cercano a
la fuente y que la red es radial (Líneas de Transmisión)

Donde:

c: es una constante que para el nivel de voltaje con el cual se simuló corresponde a
valores entre 1.1 y 1.00. Para el cálculo se tomó el valor de c=1.

Um: Valor fase tierra del voltaje equivalente.

Z: Impedancia de secuencia positiva.

 X
3 R
x: factor multiplicativo que depende la relación (R/X) x  1.02  0.98  e 33)

Reemplazamos en las ecuaciones anteriores con los datos de nuestro problema y

obtenemos que:

1  220 kV
cU m 3
I 
''
K   6.6667  kA
3Z 3  11


3 1 
x  1.02  0.98  e 6.05
 1.61686
 I P  x  2  I K''  1.61686  2  6.6667kA  15.243  kA

Así mismo la norma IEC 60909 para estudios de corto donde la falla fue desbalanceada

(Falla Monofásica) sugiere el cálculo de la corriente inicial simétrica trifásica I 

´''
K y de

la corriente pico
I 
´
P como sigue:

3  c U m
I K'' 
2  Z (1)  Z ( 0) I P  x  2  I K''
34)

Usando los datos del problema obtenemos que:

3  1  220
kV
3  c U m 3
I  ''
K   8.2956  kA
2  Z (1)  Z ( 0 ) 2  11  4.52


3 1 
x  1.02  0.98  e 6.44
 1.6351

I P  x  2  I K''  1.6351  2  8.2956kA  19.182  kA

ahora comparamos los resultados obtenidos mediante la simulación en ATP con los
calculados según la norma IEC 60909:

Simulación ATP Norma IEC 60909

Falla
Corriente pico [kA] Corriente pico [kA]
Monofásica 20.895 19.182
Trifásica 16.792 15.243
Tabla 23
De acuerdo con la siguiente tabla, podemos apreciar que la variación entre los
resultados obtenidos para la corriente de falla a través norma IEC 60909 y de la
simulación en ATP no muestran mucha diferencia entre ellas, de hecho, el error entre
ambos datos para cada falla es del 9%, lo cual nos da a pensar que la simulación de
este evento puede ser llevada acabo con alto grado de precisión mediante el uso de el
ATP, así como mediante los cálculos simplificaos propuestos por la norma IEC 60909.

[Link] DE ONDAS VIAJERAS

Como sabemos la onda de voltaje en un Sistema de Potencia se propaga a la velocidad

de la luz (300000Km/s) por lo tanto dicha onda tarda algún tiempo en propagarse por
la línea de transmisión, a la vez como el voltaje se propaga sobre una onda senoidal,
sufre todos los fenómenos que se presentan en las ondas que se propagan en un
medio entre estos fenómenos esta el de reflexión, este fenómeno lo podemos estudiar
con la ecuación del telegrafista. Una línea de transmisión fundamentalmente es una
red RLC de elementos distribuidos, donde la energía se va propagando de elemento en
elemento generándose lo que se llama una Onda Viajera. Una primera aproximación
al fenómeno de la propagación es representar la línea por una serie de elementos RLC
en cascada. Para determinar en forma exacta el comportamiento de la línea en
régimen transitorio, se plantea las ecuaciones del telegrafista.
Con el siguiente montaje se realizaran simulaciones que permitirá comprender el
comportamiento de las ondas viajeras en líneas de transmisión y la incidencia del
coeficiente de reflexión para su análisis en cualquier punto de la línea.

Figura 42. Montaje onda viajera

Longitud = 50 km
Impedancia característica = 500 Ω
Velocidad de propagación = 300000 km/seg.

En nuestro caso se analizará el comportamiento de la onda en el principio, mitad y

final de la línea de transmisión.
La solución de la ecuación de onda para nuestro caso es la siguiente:

  x  x 
V f  t   u  t    
  v  v 
   2d  x     2d  x   
 r2V f  t   u t     
   v    v  
 
zc   2d  x     2d  x  
V  x, t   *  r1r2V f  t     u t       34)
z1  zc    v    v  
 
 r r 2V  t   4d  x   u  t   4 d  x    
 1 2 f       
  v    v 
 
 2 2   4d  x     4d  x   
 r1 r2 V f  t   v
u t  
v
   ...
       

Simulaciones

1. La primera simulación corresponde a una resistencia de fuente Rf = 0.1  y una

resistencia de carga Rc = 500  , V f  100

Z1  Z c
r1  35), coeficiente de reflexión de una onda viajera de la línea hacia la
Z1  Z c
fuente.

Z 2  Zc
r2  36), coeficiente de reflexión de una onda viajera de la línea hacia la
Z2  Zc
carga.

0.1  500
r1   0.9996
0.1  500
 500  500
r2  0.
500  500

Si llamamos  el tiempo que tarda la señal a la mitad de la línea se tiene que:

25 Km
   83.3333 s
300000 Km s

Según la solución de la ecuación de onda se tiene que:

- Para X=0

500  x x 
V  x, t   V f (t  )u (t  ) 
0.1  500  v v 
V  0, t   99.98  tu (t ) 

d
- Para X= con d= 50 Km
2
500  25 25 
V  25, t   V f (t  )u (t  )
0.1  500  3E 5 3E 5 

V  25, t   99.98  (t  8.33E  5)u (t  8.33E  5) 

- Para X=d con d= 50 Km

500  50 50 
V  50, t   V f (t  )u (t  )
0.1  500  3E 5 3E 5 

V  50, t   99.98  (t  1.666E  4)u (t  1.666 E  4) 

Ahora para la simulación obtendremos cada uno de los componentes requeridos

 Figura 43 Selección de fuente tipo heidler

De la anterior obtenemos la fuente tipo (figura 3)

Figura 44 fuente heidler

Luego obtendremos la línea utilizada en la simulación.

 Figura 45

Figura 46
De lo anterior se obtiene el elemento a continuación.

Figura 47
 Figura 48

Luego para obtener la resistencia

Figura 49
 Figura 50 resistencia

100
[V]

80

60

40

20

0
0.0 0.5 1.0 1.5 2.0 [m s] 2.5
(f ile juan.pl4; x-var t) v:FTE v:LNS v:COLA

Figura 51 voltajes en el recorrido del circuito

FTE= fuente.
LNS= mitad de la línea.
COLA= final de la línea.

- Después de haberse aplicado el impulso de voltaje en t = 0, a los 25 Km de la

fuente aparece una señal “idéntica” (igual amplitud y forma) a la de la
alimentación después de transcurridos 83.55 s.
- Después de haberse aplicado el impulso de voltaje en t = 0, a los 50 Km de la
fuente aparece una señal “idéntica” (igual amplitud y forma) a la de la
alimentación después de transcurridos 166.5 s.
- Después de presentarse este impulso de voltaje en los dos puntos de observación
no se observan más perturbaciones (el voltaje permanece en cero).
Cuando en un extremo de una línea de transmisión se produce una excitación de
voltaje la señal se propaga a través de todos los elementos que forman a la línea a
una velocidad finita de tal manera que la señal tarda un tiempo en ir de la fuente a
otro punto relativamente alejado de ella. Cuando esta señal llega al extremo final de
la línea encuentra un obstáculo físico y sufre lo que ha llamado una “reflexión” cuya
magnitud depende del coeficiente de reflexión que se tenga en la carga y varía del

100% al 100%. En caso de que el coeficiente de reflexión sea diferente de cero se

crea una “onda reflejada” en dirección hacia la fuente que dependiendo de lo que se
tenga como fuente se reflejará nuevamente para crear una oscilación en la red. La
oscilación se atenuará o se mantendrá dependiendo de las características de la línea.

Como el coeficiente de reflexión r2=0, al final de la red no se refleja nada. Esto quiere
decir que la onda al encontrar una impedancia de carga igual al a impedancia propia
de la línea se comporta como si la carga de 500  hubiera absorbido todo el “impacto”
de la onda propagada.
Esta última explicación aclara porque no se vuelve a presentar ninguna otra
perturbación a los 25 y 50 km.
La verificación de los tiempos de propagación se presenta en la siguiente tabla:

Tiempo de propagación Tiempo de propagación

teórico (s) simulación (s)
X = 25 83.33 83.55
Km
X = 50 166.66 166.5
Km
Tabla 24

2. Rf = 0.1  , y Rc = 250  , V f  100

 0.1  500
r1   0.9996
0.1  500

250  500
r2   0.3333
250  500

d
Para X= con d= 50 Km
2
V  25, t   99.98  (t  8.33E  5)u (t  8.33E  5)   33.32  (t  2.5 E  4)u (t  2.5 E  4) 
33.3  (t  4.1E  4)u (t  4.1E  4) 

Para X=d con d= 50 Km

V  50, t   66.6  (t  1.66E  4)u (t  1.66 E  4)  22.21 (t  5E  4)u (t  5E  4) 

100

[V]

75

50

25

-25

-50
0.0 0.3 0.6 0.9 1.2 [ms] 1.5
(file juan.pl4; x-var t) v:FTE v:COLA v:LNS

Figura 52 voltajes en el recorrido del circuito

Se analizarán las primeras manifestaciones de la onda a lo largo de la línea.

- x = 25 km

En este caso y de manera similar a como se explicó antes, la señal de voltaje llega sin
ninguna deformación a la mitad de línea en un tiempo dado por , cuando la señal llega
a la carga a los 2 encuentra un coeficiente de reflexión de –0.333333. Esto indica que
de la señal incidente se refleja –33.33%, o sea del impulso de 100 v pico se refleja un
“impulsito” invertido de 33.33 voltios. Así, esta onda reflejada llega a la mitad de la
línea a los 3 para mostrar en este punto un impulso de -33.33 voltios pico. Este
impulso reflejado llega a la fuente a los 4 encontrando un coeficiente de reflexión de
–1. Esto indica que este impulso de refleja totalmente en la carga invirtiéndose y
llegando nuevamente a la mitad de la línea los 5 pero con un pico positivo. El proceso
se produce sucesivamente hasta que el impulso se atenúa completamente.
Comparemos estos resultados en la siguiente tabla:

Resultados Teóricos Resultados Simulación

Tiempo de llegada Tiempo de llegada
Amplitud de la Amplitud de la
de la señal de de la señal de
señal (voltios) señal (Voltios)
voltaje (s) voltaje (s)
83.33 100 83.40 100
250.00 -33.32 250.35 -33.292
416.66 33.33 416.25 33.296
Tabla 25.

x = 50 km

En t = 2 llega la señal al final de la línea en donde encuentra un coeficiente de

reflexión de –0.333, de tal manera que se refleja un impulso de valor pico –33.33
voltios. este impulso reflejado se resta con la señal incidente para dar un impulso de
voltaje en la carga de amplitud 66.66 voltios pico.
En t = 4 llega la señal a la fuente en donde encuentra un coeficiente de reflexión de –
1 de tal manera que este impulso reflejado de –33.33 voltios se refleja totalmente y se
invierte para dirigirse nuevamente hacía la carga.
En t = 6 llega la señal de 33.33 voltios pico al final de la línea y de la cual se refleja el
-33.33%. Al igual que antes la reflejada se resta a la incidente para dar origen a un
impulso en la carga de valor pico 22.2222 voltios.
Los resultados se resumen en la siguiente tabla:

Tabla 26
Resultados Teóricos Resultados Simulación
Tiempo de llegada Tiempo de llegada
Amplitud de la Amplitud de la
de la señal de de la señal de
señal (voltios) señal (Voltios)
voltaje (s) voltaje (s)
166.66 66.666 166.50 66.645
500.00 22.222 499.65 22.203

Rf = 0.1  , y Rc = 750  , V f  100

0.1  500
r1   0.9996
0.1  500

750  500
r2   0.2
750  500

d
Para X= con d= 50 Km
2
99.98  (t  8.33E  5)u (t  8.33E  5)   19.99  (t  2.5 E  4)u (t  2.5 E  4)   
V  25, t    
19.9  (t  4.1E  4)u (t  4.1E  4)  

Para X=d con d= 50 Km

V  50, t   119.97  (t  1.66 E  4)u (t  1.66 E  4)   23.98  (t  5 E  4)u (t  5 E  4) 

 comportamiento con carga de 750 ohms
120
[V]
100

80

60

40

20

-20

-40
0.0 0.3 0.6 0.9 1.2 [ms] 1.5
(f ile juan.pl4; x-var t) v:FTE v:COLA v:LNS

Figura 53 voltajes en el recorrido del circuito

- x = 25 Km
En este caso tenemos que:
En t =  llega la señal a la mitad de la línea sin sufrir ninguna deformación y
manteniendo la amplitud constante.
En t = 2 llega la señal al final de la línea en donde encuentra un coeficiente de
reflexión de 0.2, o sea que al llegar al final de la línea se refleja un impulso de voltaje
cuya amplitud corresponde al 20 % de la onda incidente para originar un impulso de
voltaje de 20 voltios pico.
En t = 3 llega la señal reflejada. Un impulso de voltaje de 20 voltios pico.
En t = 4 llega la señal reflejada a la carga en donde encuentra un coeficiente de
reflexión de –1, o sea que al llegar se refleja totalmente invertida.
En t = 5 llega la señal reflejada de la fuente a través de un impulso de voltaje de –20
voltios pico de amplitud.

Los resultados se expresan en la siguiente tabla:

Resultados Teóricos Resultados Simulación

Tiempo de llegada Tiempo de llegada
Amplitud de la Amplitud de la
de la señal de de la señal de
señal (voltios) señal (Voltios)
voltaje (s) voltaje (s)
83.33 99.98 83.25 100.00
250.00 19.99 249.90 20.023
416.66 -19.95 416.40 -19.997
Tabla 27

- x = 50 Km
En t = 2 llega la señal al final de la línea en donde encuentra un coeficiente de
reflexión de 0.2. esto implica, tal como se ha explicado antes que la señal reflejada se
sume con la incidente para dar origen a un impulso de voltaje parecido al de la fuente
pero con amplitud de 120 voltios.
En t = 4 llega la señal a la fuente en donde encuentra un coeficiente de reflexión de –
1. El impulso de voltaje de amplitud igual a 20 voltios se refleja invertido en la fuente.
En t = 6 llega la señal reflejada de la fuente a la carga en donde nuevamente se
refleja solo el 20 % de la señal. Esto da origen a un impulso de voltaje de amplitud –
24 voltios.
Los resultados se expresan en la siguiente tabla:

Resultados Teóricos Resultados Simulación

Tiempo de llegada Tiempo de llegada
Amplitud de la Amplitud de la
de la señal de de la señal de
señal (voltios) señal (Voltios)
voltaje (s) voltaje (s)
166.66 119.97 166.50 119.96
500.00 -23.98 499.95 -23.98
Tabla 28

Rf = 0.1  , y Rc = 1 M  , V f  100

0.1  500
r1   0.9996
0.1  500

1E 6  500
r2   0.9990
1E 6  500
 d
Para X= con d= 50 Km
2
99.98  (t  8.33E  5)u (t  8.33E  5)   99.88  (t  2.5 E  4)u (t  2.5 E  4)   
V  25, t    
99.84  (t  4.1E  4)u (t  4.1E  4)  

Para X=d con d= 50 Km

V  50, t   199.86  (t  1.66E  4)u (t  1.66 E  4)   199.58  (t  5E  4)u(t  5E  4) 

comportamiento con carga de 1M ohm

200
[V]
150

100

50

-50

-100

-150

-200
0.0 0.3 0.6 0.9 1.2 [ms] 1.5
(f ile juan.pl4; x-var t) v:FTE v:COLA v:LNS

Figura 54 voltajes en el recorrido del circuito

En este caso se tiene que:

- x = 25 Km

En t =  llega la señal a la mitad de la línea sin sufrir ninguna deformación y

manteniendo la amplitud constante.
En t = 2 llega la señal al final de la línea en donde encuentra un coeficiente de
reflexión de 1, o sea que al llegar al final de la línea se refleja un impulso de voltaje de
amplitud igual al incidente, o sea de 100 voltios.
En t = 3 llega la señal reflejada. Un impulso de voltaje de 100 voltios pico.
En t = 4 llega la señal reflejada a la carga en donde encuentra un coeficiente de
reflexión de –1, o sea que al llegar se refleja totalmente invertida.
En t = 5 llega la señal reflejada de la fuente a través de un impulso de voltaje de –
100 voltios pico de amplitud.
En t = 6 llega la señal a larga nuevamente en donde se refleja toda la señal incidente.
En t = 7 llega la señal reflejada de la carga a través de un impulso de voltaje de -100
voltios pico.

Los resultados se expresan en la siguiente tabla:

Resultados Teóricos Resultados Simulación

Tiempo de llegada Tiempo de llegada
Amplitud de la Amplitud de la
de la señal de de la señal de
señal (voltios) señal (Voltios)
voltaje (s) voltaje (s)
83.33 99.98 83.25 99.975
250.00 99.88 249.75 99.897
416.66 -99.84 416.40 -99.787
Tabla 29

- x = 50 Km

En t = 2 llega la señal al final de la línea en donde encuentra un coeficiente de

reflexión de 1, esto implica, tal como se ha explicado antes que la señal reflejada se
sume con la incidente para dar origen a un impulso de voltaje de dos veces el de la
fuente con amplitud de 200 voltios.
En t = 4 llega la señal a la fuente en donde encuentra un coeficiente de reflexión de –
1. El impulso de voltaje de amplitud igual a 100 voltios se refleja invertido en la
fuente.
En t = 6 llega la señal reflejada de la fuente a la carga en donde nuevamente se
refleja totalmente. Esto da origen a un impulso de voltaje de amplitud –200 voltios.
O sea que teóricamente en este caso lo que tendremos es un impulso de voltaje de
amplitud  200 Voltios que oscila indefinidamente en la línea de transmisión de un
extremo a otro.

Los resultados se expresan en la siguiente tabla:

Resultados Teóricos Resultados Simulación

Tiempo de llegada Tiempo de llegada
Amplitud de la Amplitud de la
de la señal de de la señal de
señal (voltios) señal (Voltios)
voltaje (s) voltaje (s)
166.66 199.86 166.50 199.83
500.00 -199.58 499.80 -199.52
Tabla 30

En este punto vale la pena resaltar que el anterior análisis es solo un resultado teórico
pues es imposible que en una línea real después de haber realizado una pequeña
excitación en un extremo de ella se tenga una señal oscilante indefinidamente a menos
que exista una fuente que permanentemente se encuentra excitando a la red.

3. Rf = 500  , y Rc = 1 M  , V f  100

500  500
r1  0
500  500

1E 6  500
r2   0.9990
1E 6  500

500
Amplitud  100  50
500  500
d
Para X= con d= 50 Km
2
V  25, t   50  (t  8.33E  5)u (t  8.33E  5)   49.95  (t  2.5 E  4)u (t  2.5 E  4) 
Para X=d con d= 50 Km

V  50, t   99.95  (t  1.66 E  4)u (t  1.66 E  4) 

comportamiento con carga de 1M ohm y R de fuente de 500 ohms

100
[V]

80

60

40

20

0
0.0 0.3 0.6 0.9 1.2 [ms] 1.5
(file juan.pl4; x-var t) v:FTE v:COLA v:LNS

Figura 55 voltajes en el recorrido del circuito

- x = 25 Km

En t =  y tal como lo hemos planteado anteriormente, llega la señal a la mitad de la

línea sin sufrir ninguna deformación y manteniendo la amplitud constante. (50 voltios)
En t = 2 llega la señal al final de la línea en donde encuentra un coeficiente de
reflexión de 1, por tanto al final de la línea se refleja un impulso de voltaje de amplitud
igual al incidente, o sea de 50 voltios. En este caso se dice que la señal incidente se
refleja totalmente.
En t = 3 llega la señal reflejada a la mitad de la línea. Un impulso de voltaje de
50voltios pico.
En t = 4 llega la señal reflejada a la carga en donde encuentra un coeficiente de
reflexión de 0, lo cual indica que no hay reflexión en la carga. O sea, es como si en la
fuente se absorbiera todo el “impacto” de la señal incidente.
A partir de este momento desaparece de la línea cualquier señal.

Los resultados se pueden resumir en la siguiente tabla:

Resultados Teóricos Resultados Simulación

Tiempo de llegada Tiempo de llegada

Amplitud de la Amplitud de la
de la señal de de la señal de
señal (voltios) señal (Voltios)
voltaje (s) voltaje (s)
83.33 50.00 83.25 49.998
250.00 49.95 249.75 49.958
Tabla 31.

x = 50 Km

En t = 2 llega la señal al final de la línea en donde encuentra un coeficiente de

reflexión de 1, esto implica, tal como se ha explicado antes que la señal reflejada se
sume con la incidente para dar origen a un impulso de voltaje de dos veces el de la
fuente con amplitud de 100 voltios.
En t = 4 llega la señal a la fuente en donde encuentra un coeficiente de reflexión de
0, por tanto pasa lo mencionado anteriormente que la señal incidente se “disipa”
completamente y desaparece de la línea.

Los resultados se expresan en la siguiente tabla:

Resultados Teóricos Resultados Simulación

Tiempo de llegada Tiempo de llegada
Amplitud de la Amplitud de la
de la señal de de la señal de
señal (voltios) señal (Voltios)
voltaje (s) voltaje (s)
166.66 99.95 166.50 99.93
Tabla 32

6 ENERGIZACION DE LINEA MONOFASICA SIN PÉRDIDAS EN VACIO.

 U

Figura 56

Se energiza la línea con una fuente DC que con la ayuda de un swiche (Cerrando en t
= 0) simula una función escalón. El modelo con que representaremos a la línea
monofásica sin pérdidas es el mismo con el que se viene trabajando: Modelo de Clarke
monofásico.
Para asignar los setings en el ATP se debe tener en cuenta la duración de viaje de la
señal de voltaje en la línea de 100 Km.

Figura 58 (Extraída de la guía del curso por Leonardo Cardona)

Lo que se desea es observar como se comporta el voltaje al final de la línea.

Mediante simulación en ATP se obtiene el siguiente resultado:
 200
[V]

160

120

80

40

0
0 3 6 9 12 [ms] 15
(f ile EnergiVacio.pl4; x-v ar t) v :INI v :FIN

Figura 59

La anterior gráfica se muestra las señales de voltaje en la fuente de alimentación y en

el final de la línea. Podemos observar que la señal de voltaje no aparece
instantáneamente en el extremo final de la línea después de cerrar el circuito, sino que
existe un tiempo de propagación durante el cual la señal “va” de la fuente al extremo
de la red.
Este tiempo estará dado por:
100 Km
t propagación     333.33us
300000 Km
s
Efectivamente se puede ver que en la simulación la señal aparece en el extremo de la
línea en t=333.33 us.

Al analizar la respuesta oscilatoria de acuerdo la teoría de Onda Viajera, es necesario

definir los coeficientes de reflexión.

Z 2  Z c   500
r2   1
Z 2  Z c   500

Z 1  Z c 0  500
r1    1
Z 1  Z c 0  500

El tiempo que tarda la onda de voltaje en propagarse desde la fuente hasta el final de
la línea es  = 333.33 us.
Cuando t =  la señal llega al final de la línea en donde encuentra un coeficiente de
reflexión de 1. Esto indica que en este punto la onda de voltaje incidente sufre una
reflexión total, de tal manera que esta última se suma con la señal reflejada para
formar una señal de doble amplitud con respecto a la de la fuente. Luego en el
extremo final de la línea se tiene un sobrevoltaje por lo antes mencionado.
Sucesivamente la onda reflejada se dirige en dirección de la fuente. Finalmente en
t = 2 llega a la fuente en donde encuentra un coeficiente de reflexión de -1. Esto
indica que la señal se refleja nuevamente pero restándose con la amplitud de la fuente
Esta onda de amplitud 100 voltios al llegar a la carga en t = 3 se refleja totalmente
para originar una señal de amplitud –200 que se suma con los 200 voltios ya presentes
en ese punto. Luego en la carga en el intervalo 3 < t < 5 se tienen 0 voltios.
Los resultados se resumen en la siguiente tabla:

Resultados Teóricos Resultados Simulación

Tiempo de Amplitud de la Tiempo de Amplitud de la
llegada de la señal (voltios) llegada de la señal (Voltios)
señal de voltaje señal de voltaje
(s) (s)
333.33 200.00 333.3 200.00
999.99 0.00 999.9 0.00
Tabla 33.

7. ENERGIZACION DE LINEA MONOFASICA SIN PÉRDIDAS CON CARGA

Figura 60
CON R=1000

Calculemos los coeficientes de reflexion:

Z 2  Z c 1000  500 1
r2     0.33333
Z 2  Z c 1000  500 3

Z 1  Z c 0  500
r1    1
Z 1  Z c 0  500
Simulando el circuito de la figura 10 en ATP tenemos.

160
[V]
140

120

100

80

60

40

20

0
0 3 6 9 12 [ms] 15
(f ile EnergiCarga.pl4; x-v ar t) v :INI v :FIN

Figura 61

De acuerdo con la teoría de Onda Viajera.

La onda de voltaje demora un tiempo en propagarse hasta la carga, cuando esta llega
a la carga encuentra un coeficiente de reflexión de 0.3333, quiere decir que la onda
incidente se refleja en la carga el 33.33%. Esta onda reflexión se superpone con la
incidente para y forma una onda de 133.33 voltios de amplitud. Sucesivamente la
señal reflejada de 33.33 voltios se dirige hasta la fuente allí el coeficiente de reflexión
es –1, la onda se refleja de forma invertida. Estos –33.33 voltios llegan a la carga en t
= 3 punto en donde se refleja el 33.33% de la señal que tiene un valor de 133.33V,
esto da origen a una señal de voltaje de amplitud igual a –43.99 voltios. Estos se
suman con 133.33 ya existentes para dar origen en la carga a un voltaje de 88.88
voltios. Así sucesivamente las reflexiones se van haciendo mas pequeñas hasta
estabilizarse en el valor de la fuente.
 Resultados Simulación
Resultados Teóricos

Tiempo de Amplitud de la Tiempo de Amplitud de la

llegada de la señal (voltios) llegada de la señal (Voltios)
señal de voltaje señal de voltaje
(s) (s)
333.33 133.33 330 133.33
999.99 88.88 1030 88.88
1666.65 103.703 1660 103.70
Tabla 34.

CON R=100

Tenemos:

Z 2  Z c 100  500 2 Z 1  Z c 0  500

r2      0.66666 r1    1
Z 2  Z c 100  500 3 Z 1  Z c 0  500
 100
[V]

80

60

40

20

0
0 3 6 9 12 [ms] 15
(f ile EnergiCarga.pl4; x-v ar t) v :INI v :FIN

Figura 62
La onda de voltaje llega al final de la línea después de transcurrir un tiempo , allí se
refleja el -66.66% del valor de el voltaje de la fuente, luego la señal incidente se
superpone con la reflejada para formar una señal en la carga de 33.333 voltios. Estos
-66.66 voltios reflejados llegan a la fuente en t = 2 en donde encuentra un
coeficiente de reflexión de –1. Esto permite que en la carga se refleje una señal de
66.66 voltios. La onda relejada de 66.66 voltios de amplitud llega a la carga en t = 3
donde se produce una reflexión de –44.44 voltios; esta onda reflejada se suma con la incidente y
con el voltaje ya presente en el punto para dar origen a una señal de voltaje de 55.55 voltios de
amplitud.. El proceso se repite sucesivamente hasta que la onda de voltaje se llega al valor de la
fuente.
Resultados Simulación
Resultados Teóricos

Tiempo de Amplitud de la Tiempo de Amplitud de la

llegada de la señal (voltios) llegada de la señal (Voltios)
señal de voltaje señal de voltaje
(s) (s)
333.3 33.3 330 33.35
999.3 55.5 990 55.55
1667.0 70.4 1670 70.35
2333.3 80.3 2310 80.28
2999.9 86.8 2980 86.85
Tabla 35.
CON R=500

Tenemos

Z 2  Z c 500  500
r2   0
Z 2  Z c 500  500

Z 1  Z c 0  500
Y : r1    1
Z 1  Z c 0  500

100
[V]

80

60

40

20

0
0 3 6 9 12 [ms] 15
(f ile EnergiCarga.pl4; x-v ar t) v :INI v :FIN

Figura 63

Cuando la onda se propaga hasta la carga, allí se encuentra un coeficiente de reflexión

cero, quiere decir que no se presenta el fenómeno de reflexión, esto ocurre cuando la
Zcarga = Zcaracteristica

7. ENERGIZACIÓN DE LINEA MONOFASICA SIN PÉRDIDAS CON RESISTENCIA

DE PREINSERCION
Una manera de controlar la magnitud del sobrevoltaje que se presenta al final de una
línea durante su energización es mediante una resistencia de preinserción. El objetivo
de esta resistencia es hacer que la onda viajera que regresa de la línea hacia la fuente
vea un coeficiente de reflexión muy pequeño.

Figura 64

A la resistencia de preinserción para efectos de la simulación se le debe controlar dos

parámetros básicos:
 Valor de la resistencia.
 Tiempo en que debe permanecer conectada en serie con la línea.

Para la selección de esta resistencia se debe conocer primero que todo tres elementos
importantes de la línea:
- Velocidad de propagación.
- Longitud de la línea.
- Impedancia característica.
Para este caso :
Datos

v  300000 Km L  100 Km Z C  500 

s

Conociendo esto, definimos:

Con la resistencia de preinserción buscamos evitar sobrevoltaje al final de la línea
durante el proceso de energización. Una línea abierta se tiene un coeficiente de
reflexión en la carga de 1, por lo tanto en el extremo será dos veces el voltaje
incidente. Luego para que el voltaje en el extremo final de la línea sea el nominal debe
asegurarse que el voltaje que vea la línea en el mismo momento de energizarse sea la
mitad del de la fuente. Esto se consigue con el “divisor de onda” que se establece con
la resistencia de inserción, tal que:
Z1
Voltaje incidente  V fuente 37)
Z1  Z C

Luego se debe garantizar que: Z 1  Z C

Además de esto también se debe de garantizar que cuando la onda relejada al final de
la línea llegue a la fuente nuevamente, encuentre un coeficiente de reflexión cero (0).

Z1  Z C
r1   0 , si Z1  Z C .
Z1  Z C
Por lo tanto el valor de la resistencia de preincersion debe ser igual a la impedancia
característica. Ahora evaluemos el tiempo que debe permanecer esta resistencia en la
línea.

La resistencia de preincersion se debe desconectar después de 2. Como  =

333.33μs, 2 =666.66 μs. Abramos el interruptor en t = 670μs

La simulación en ATP arroja los siguientes resultados.

 120
[V]
100

80

60

40

20

0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 [ms] 3.5
(file pablo1.pl4; x-var t) v:TRES v:DOS

Figura 65

8. OBTENCIÓN DE MODELOS LEE Y CLARK CON LINE CONSTANTS.

Los modelos de onda viajera para línea transpuesta y no transpuesta son los modelos
más utilizados para representar adecuadamente una red de transmisión de energía
eléctrica para diferentes estudios de un sistema de potencia.
Para la práctica se van a calcular los modos de propagación a partir de los parámetros
de secuencia.
Para ello utilizaremos las siguientes ecuaciones:


Rmod o  RSecuencia 
Km
 (38)

LSecuencia  H 
Z C Modo    (39)
C Secuencia (F )

1  Km 
Vel . Pr opa. 
Ind . de sec uencia H  Km
  Capac. de sec uencia F Km  Seg.  (40

)
Lo primero que se debe hacer es calcular las impedancias de secuencia para la línea de
transmisión real que se va a estudiar.
Sea una línea de transmisión de 500 Kv (del Sistema Interconectado Nacional) que
tiene la siguiente configuración:

Figura 66 disposición física Torre de 230 KV para transmisión

Se busca calcular los elementos de secuencia para la red.

Los datos físicos para el conductor de fase y del cable de guarda son:

Conductor del cable de

Conductor de guarda
fase

25.17 8.71
Diámetro (mm)

0.08912 1.844
Rdc (/Km)

Tabla 36.

Donde se usan las siguientes ecuaciones:

Z0 = Impedancia de secuencia cero.
Z1 = Impedancia de secuencia positiva (que es igual a la de secuencia negativa)
C0 = Capacitancia de secuencia cero.
C1 = Capacitancia de secuencia positiva.

Las siguientes ecuaciones resumen el cálculo de parámetros en líneas de

transmisión.

w DMG
Z 1 = Zs  Zm  Rac  j ln  / km 41)
2 RMG 

3 Z 2fg
Z0 = ( ZS +2 ZM )- 42)
Zg

w De
Zs = Rac  0.0592  j ln  / km 43)
2 RMG fase
 w De
Zm = 0.0592  j ln  / km 44)
2 DMG fases

w De
Zs  2 Zm = Rac  3 * 0.0592  j 3  ln  / km 45)
2 3
RMG  . DMG 2
De = 658.86  / f metros 46)

w De
Zfg = 0.0592  j ln  / km 47)
2 DMGfg

w De
Zg = Rac  0.0592  j ln  / km 48)
2 
RMGcg

2  55.55 nF
C1 = 
DMG DMG km 49)
ln ln
rf rf
1 2.
C0 =  2
3  DMG fases imagen guarda 
 ln 
DMG fases imágenes  DMG  50)
 
fases  guarda
ln
3 r f .DMG 2 DMG guarda imagen guarda
ln
rg

Se presentaran cálculos para frecuencias de 60, 100, 500 y 1000 Hz y con una
resistividad del terreno  = 300 *m.

Se debe recordar para estos cálculos que   1,25664E - 06 H  m

Radio del conductor equivalente al haz de conductores de las fases del sistema de
500kV está dado por:

RMG '  n
nrA n 1 51)
 Frecuencia Zfg (Ω/Km) Zgg (Ω/Km) Zg (Ω/Km)
(Hz)

60 0.06+0,34i 0.0592+0,3074i 1,9032+0,979i

100 0.06+0,53i 0.0592+0,4802i 1,9032+1,599i
500 0.06+2,165i 0.0592+1,895i 1,9032+7,49i
1000 0.06+3,89i 0.0592+3,355i 1,9032+14,55i
Tabla 37.

Impedancias y capacitancias de secuencia:

Frecuencia Z0 (Ω/Km) Z1 (Ω/Km) C1 C0

(Hz) (nF/Km) (nF/Km)

60 0,4497+1,1161i 0.08912+0,332i 7,7889 5,453

100 0,5760+1,6476i 0.08912+0,554i 7,7889 5,453
500 0,7095+6,42536i 0.08912+2,7702i 7,7889 5,453
1000 0,664+12,308i 0.08912+5,54032 7,7889 5,453
i
Tabla 38.
Se obtiene que:

R1 R0 L1(H) C1(F) L0(H) C0(F) Zc1 Zc0 Vel1 Vel0

Frec. (Ohm/Km) (Ohm/Km) (Ohm) (Ohm) (Km/sg) (Km/sg)
(Hz)

60 0,089 0,449676 0,18429 1,628E- 0,6187 1,14E- 336,46 736,82 381581,3 248878,47
06 06 5
100 0,089 0,576031 0,18429 1,628E- 0,548 1,14E- 336,46 693,44 381581,3 264446,57
06 06 5
500 0,089 0,70947 0,18429 1,628E- 0,42746 1,14E- 336,46 612,42 381581,3 299433,45
06 06 5
1000 0,089 0,664530 0,18429 1,628E- 0,4094 1,14E- 336,46 599,35 381581,3 305960,29
06 06 5
Tabla 39.
La tabla 16 muestra en resumen los modos de propagación para el sistema de
transmisión en estudio.
Ahora se hará el cálculo de estos parámetros usando la aplicacion Line Constant que
calcula estos valores para la línea en estudio (totalmente transpuesta). Se hará la
obtención para una resistividad del terreno de 300 *m y para el barrido de
frecuencias anterior.
El modelo usado en el ATP es el dado por la siguiente gráfica:

Figura 67

Se obtienen los siguientes resultados:

Frecuencia R1 Zc1 Vel1 R0 Zc0 Vel0

(Hz) (Ω/Km) (Ω) (Km/s) (Ω/Km) (Ω) (Km/s)
60 0.02396 258.406 293548.0 0.382217 628.731 210925.0
100 0.02628 258.135 293856.0 0.561509 585.392 226541.0
500 0.04621 256.640 295568.0 1.01812 519.244 255401.0
1000 0.06036 256.137 296148.0 1.39648 510.809 259618.0
Tabla 40.

Los cálculos que arrojó el ATP se diferencian un poco de los obtenidos analíticamente
debe ser por causa de la aproximación Lewis que se hace en las ecuaciones analíticas
que se presentaron (que es válida para baja frecuencia). La impedancia característica
del modo positivo no varia con la frecuencia para los valores calculados analíticamente
lo que es lógico si se tiene en cuenta que es el cociente de dos “parámetros” de la
línea que dependen solo de la configuración. Para el caso de los valores arrojados por
el ATP tenemos que estos varían con la frecuencia.

10 Obtener el modelo de onda viajera trifásico para línea transpuesta.

El estudio de ondas viajeras arrojan resultados diferentes si la línea es transpuesta o
no, se usan modelos diferentes para cada caso (Clark y Lee). Para ello se usaran el
aplicativo Line Constants que incluye los dos modelos dependiendo si se hace una
selección transpuesta o no transpuesta.
Este aplicativo es muy importante conocer la frecuencia del fenómeno, o sea “la
frecuencia natural” de la línea. Para este cálculo debemos recordar que:

1 1 1
f     358.8517 Hz  360 Hz
4 L  Km 209 Km
4
v 
Km
s
4
300000 Km
s

Así teniendo en cuenta esta sugerencia procedemos a hallar lo valores de impedancia

característica, velocidad de propagación y resistencia de pérdidas para los diferentes
modos de propagación.

R1 Zc1 Vel1 R0 Zc0 Vel0

(Ω/Km) (Ω) (Km/s) (Ω/Km) (Ω) (Km/s)

0.0406595 256.913 295253.0 0.846 513.251 258383.0

Tabla 41

Los modelos de Onda Viajera de Clarke y de Lee lo que hacen es descomponer las
respuesta transitoria del fenómeno en Modos de Propagación (modos aire y modos
tierra) de tal manera que de un sistema donde existen acoples se obtengan tres
sistemas (para este caso) desacoplados entre sí que simplifican enormemente las
ecuaciones de solución del problema

Respuesta de voltaje al final de la línea de 500kV en el momento de la energizacion.

 Figura 68

De esta gráfica se tiene que al igual como se ha observado antes tenemos que la señal
de voltaje no aparece al final de la línea inmediatamente se energiza la red sino
después de un tiempo de propagación de 707.8 μs.
Obsérvese que al inicio de la energizacion se produce un sobrevoltaje que se va
atenuando con el tiempo, esto se debe a la propagación y reflexión de la onda a lo
largo de la línea

10 transformación interna que utiliza el programa ATP en comparación con

la transformación de componentes simétricas.

El programa ATP tiene como herramienta de análisis de algunos problemas transitorios

los modelos de Clarke y de Lee. El primero usado para sistemas totalmente
transpuestos y el segundo para sistemas sin transponer.
El modelo de Clarke usa para el análisis la descomposición del modelo real en dos
modos de propagación: modos tierra y modos aire que resultan ser análogos a las
descomposición en componentes simétricas que se hace de un sistema trifásico
desbalanceado.
El modelo CLARKE considera

V   T  V 
fase 0 y  I   T   I 
fase 0 51)
Tenemos

V   T  V 
0
1
fase y  I   T   I  52)
0
1
fase

V0 

Con: V   
0  V
V 


52)
 

La matriz T  se conoce como matriz de transformación. Esta matriz en su forma

general y para un sistema de M fases se puede escribir como:

53)

Pero en nuestro caso, como es un sistema trifásico se puede escribir como:

54)
Se sabe de la teoría de componentes simétricas que la matriz de transformación tiene
entradas imaginarias (fasores) que solo es aplicable para el estudio en estado estable
(no transitorio) de sistemas n-fásicos. Por estas dos razones es que no se puede
aplicar la matriz de transformación de componentes simétricas para el análisis de
modos de propagación en comportamientos transitorios.

11 Determinar la sensibilidad de los parámetros: Resistencia de pérdidas,

impedancia característica y Velocidad de Propagación, con respecto a
variación de la frecuencia (10 H ≤ f ≤ 5000 Hz) y la resistividad del terreno
(100 – 3000 Ωm).

Teniendo en cuenta el aplicativo LCC con el que se viene trabajando lo que se pretende
es mirar como varían los parámetros: Resistencia de pérdidas, impedancia
característica y Velocidad de Propagación para variaciones de frecuencia (manteniendo
ρ = cte. = 100 Ω/m) y como varían con la resistividad del terreno manteniendo la
frecuencia f = cte. = 60 Hz.

Tenemos
 ρ = 100 Ωm

Frecuencia R1 Zc1 Vel1 R0 Zc0 Vel0

(Hz) (Ohm/Km) (Ohm) (Km/sg (Ohm/Km) (Ohm) (Km/sg)
10 0,022335 258,613 293313 0,0621507 705,006 188105
60 0,023978 258,403 293551 0,339583 608,665 217879
100 0,02628 258,135 293856,0 0,561509 585,392 226541,0
300 0,0378392 257,090 295051 0,796516 517,327 256347
500 0,04621 256,640 295568,0 1,01812 519,244 255401,0
1000 0,06036 256,137 296148,0 1,39648 510,809 259618,0
1500 0,0733177 255,902 296420 1,72254 494221 268332
2000 0,0839495 255,776 296566 2,10795 491,013 270085
2500 0,0939642 255,688 296668 2,48870 488,515 271466
3000 0,103586 255,622 296745 2,86357 486,467 272609
3500 0,112934 255,568 296807 3,23213 484,733 273584
4000 0,122078 255,524 296859 3,59429 483,231 274434
4500 0,131062 255,486 296903 3,95011 481,909 275187
5000 0,13991 255,453 296942 4,29972 480,730 275862
Tabla 42
 R1 en función de la frecuencia

0,16

0,14

Resistencia (Ohm/Km)
0,12

0,1

0,08

0,06

0,04

0,02

0
0 1000 2000 3000 4000 5000 6000

Frecuencia (Hz)

Figura 68

R0 en función de la frecuencia

4,5
Resistencia (Ohm/Km)

3,5

2,5

1,5

0,5

0
0 1000 2000 3000 4000 5000 6000

Frecuencia (Hz)

Figura 69
 Zc1 en función de la frecuencia

259

258,5

258

257,5

257

256,5

256

255,5

255
0 1000 2000 3000 4000 5000 6000

F re c ue nc ia ( Hz )

Figura 70

Zc0 en función de la frecuencia

75 0

70 0

65 0

60 0

55 0

50 0

45 0
0 5 00 10 00 15 00 20 00 25 00 30 00 35 00 4 000 45 00 5 000

F r e c u e n c i a ( Hz )

Figura 71

V1 en función de la frecuencia

30 000 0

29 900 0
Velocidad (Km/s)

29 800 0
29 700 0

29 600 0

29 500 0

29 400 0

29 300 0
29 200 0

29 100 0

29 000 0
0 5 00 1 000 1 500 2 000 25 00 30 00 3 500 4 000 4 500 50 00

F r e c u e n c i a ( Hz )

Figura 72
 Vo en función de la frecuencia

280000
270000

Velocidad (Km/s)
260000

250000
240000

230000

220000
210000

200000

190000
180000
0 500 1000 1500 2000 2500 3000 3500 4000 4500 5000

F re c ue nc ia ( Hz )

Figura73

 f = 60 Hz
R1 Zc1 Vel1 R0 Zc0 Vel0
ρ (Ωm)
(Ohm/Km) (Ohm) (Km/sg (Ohm/Km) (Ohm) (Km/sg)
100 0,0239780 258,403 293551 0,339583 608,665 217879
300 0,0239634 258,406 293548 0,382217 628,731 210925
500 0,0239598 258,407 293547 0,402696 637,584 207996
800 0,0239575 258,407 293547 0,421890 645,443 205464
1200 0,0239561 258,407 293546 0,438704 651,996 203399
1500 0,0239556 258,407 293546 0,448052 655,513 202308
1800 0,0239551 258,408 293546 0,455739 658,339 201439
2100 0,0239548 258,408 293546 0,462270 660,695 200721
2400 0,0239546 258,408 293546 0,467950 662,711 200110
2700 0,0239544 258,408 293546 0,472978 664,470 199581
3000 0,0239543 258,408 293546 0,477490 666,029 199113
Tabla 43
 R1 en función de la resistividad del
terreno
0, 02398

0, 0 239 75

0, 02397

0, 0 239 65

0, 02396

0, 0 239 55

0, 02395
0 50 0 10 00 1500 2 000 2 500 30 00

F r e c u e n c i a ( Hz )

Figura 74

R0 en función de la resistividad del terreno

0, 6

0, 5

0, 4

0, 3

0, 2

0, 1

0 50 0 1000 1 500 2 000 2 500 30 00

F r e c u e n c i a ( Hz )

Figura 74

Zc1 en función de la resistividad del

terreno

2 60
25 9, 5

2 59
25 8, 5
2 58
25 7, 5
2 57

25 6, 5
2 56

25 5, 5
2 55

0 5 00 10 00 15 00 2 000 2 500 3000

F r e c u e n c i a ( Hz )

Figura 75
 C o m p o r t a m i e n t o d e Zc 0 e n f u n c i ó n d e l a
re s ist iv ida d de l t e rre n o

670

660

650

640

630

620

6 10

600

0 500 10 0 0 15 0 0 2000 2500 3000

Fr e c u e n c i a ( Hz )

Figura 76

V1 en función de la resistividad del terreno

295000
294500
294000
293500
293000
292500
292000
291500
291000
290500
290000
0 500 1000 1500 2000 2500 3000

Fr e c u e n c i a ( Hz )

Figura 77

Com portam iento de V0 en función de la

resistividad del terreno

220000

215000

210000

205000

200000

195000

190000
0 500 1000 1500 2000 2500 3000

F re c ue nc ia ( Hz )

Figura 78

En este caso cuando la frecuencia de operación se mantiene constante y solo se varía

la resistividad del terreno se tienen los siguientes resultados:
 Los parámetros en el modo de propagación aire se mantienen más o menos
constantes.
  La resistencia de pérdidas y la impedancia característica en el modo de
propagación tierra aumentan con el aumento de la resistividad del terreno.
 La velocidad de propagación disminuye con el aumento de la resistividad del
terreno.

12 Determinar la relación entre Zc y R/4 (con R la resistencia total de

pérdidas en la línea).

Frecuencia Zc1 R1/4 10R1/4 Zc0 R0/4 10R0/4

(Hz) (Ohm) (Ohm) (Ohm) (Ohm) (Ohm) (Ohm)
10 258,613 1,16700375 11,6700375 705,006 3,24737408 32,4737408
60 258,403 1,2528505 12,528505 608,665 17,7432118 177,432118
100 258,135 1,37313 13,7313 585,392 29,3388453 293,388453
300 257,090 1,9770982 19,770982 517,327 41,617961 416,17961
500 256,640 2,4144725 24,144725 519,244 53,19677 531,9677
1000 256,137 3,15381 31,5381 510,809 72,96608 729,6608
1500 255,902 3,83084983 38,3084983 494221 90,002715 900,02715
2000 255,776 4,38636138 43,8636138 491,013 110,140388 1101,40388
2500 255,688 4,90962945 49,0962945 488,515 130,034575 1300,34575
3000 255,622 5,4123685 54,123685 486,467 149,621533 1496,21533
3500 255,568 5,9008015 59,008015 484,733 168,878793 1688,78793
4000 255,524 6,3785755 63,785755 483,231 187,801653 1878,01653
4500 255,486 6,8479895 68,479895 481,909 206,393248 2063,93248
5000 255,453 7,3102975 73,102975 480,730 224,66037 2246,6037
Tabla 44
Zc1 R1/4 10 R1/4 Zc0 R0/4 10R0/4
ρ (Ω*m)
(Ohm) (Ohm) (Ohm) (Ohm) (Ohm) (Ohm)
100 258,403 1,2528505 12,528505 608,665 17,7432118 177,432118
300 258,406 1,25208765 12,5208765 628,731 19,9708383 199,708383
500 258,407 1,25189955 12,5189955 637,584 21,040866 210,40866
800 258,407 1,25177938 12,5177938 645,443 22,0437525 220,437525
1200 258,407 1,25170623 12,5170623 651,996 22,922284 229,22284
1500 258,407 1,2516801 12,516801 655,513 23,410717 234,10717
1800 258,408 1,25165398 12,5165398 658,339 23,8123628 238,123628
2100 258,408 1,2516383 12,516383 660,695 24,1536075 241,536075
2400 258,408 1,25162785 12,5162785 662,711 24,4503875 244,503875
2700 258,408 1,2516174 12,516174 664,470 24,7131005 247,131005
3000 258,408 1,25161218 12,5161218 666,029 24,9488525 249,488525
Tabla 45
 En el estudio de líneas de transmisión se tiene que una línea que en principio está
representada por parámetros distribuidos, se puede concentrar en ciertos tramos
de red.
En estas condiciones se tendría que en un tramo:
R
Z C equivalente  Z C  55)
4
Donde Zc es la impedancia característica de toda la línea de transmisión. Para que
esta impedancia característica sea la misma en todos los tramos, se debe
garantizar que el efecto de R/4 sea mínimo, esto es que: Zc >> R/4 (por lo menos
Zc > 10(R/4)). Si esto no se cumple se deben tomar tramos de red más pequeños
para el estudio de la línea.
De los anteriores resultados se observa que para el modo aire se cumple esta
relación para cualquier valor de frecuencia y cualquier valores de resistividad del
terreno dentro del rango que se está estudiando.
Para el modo cero esta relación deja de ser válida para frecuencias por encima de
los 300 Hz. Esto sugiere entonces que se tomen tramos concentrados más
pequeños para el estudio de la red.

Obtener el modelo de Onda Viajera triásico para una línea no transpuesta.

Figura 79
Ojo valores como se ingresan

Los resultados arrojados por el ATP son

The transformation matrix was calculated at 3.60000000E+02 Hz.

$VINTAGE, 1
-1IN___AOUT__A 8.41461E-01 5.09889E+02 2.58703E+05- 2.09000E+02 1
3
-2IN___BOUT__B 5.05505E-02 2.81925E+02 2.95737E+05- 2.09000E+02 1
3
-3IN___COUT__C 3.07708E-02 2.34289E+02 2.97800E+05- 2.09000E+02 1
3
$VINTAGE, 0
0.60551764 -0.70710678 -0.40941809
-0.00188920 0.00000000 -0.00147788
0.51640631 0.00000000 0.81532061
0.00376219 0.00000000 -0.00126038
0.60551764 0.70710678 -0.40941809
-0.00188920 0.00000000 -0.00147788

Obsérvese los datos mas importantes resumidos en la siguiente tabla

R0 0.841461
(Ohm/Km)
Zc0 (Ohm) 509.889
V0 (Km/s) 258703.0
R1 0.0505505
(Ohm/Km)
Zc1 (Ohm) 281.925
V1 (Km/s) 295737.0
R2 0.0307708
(Ohm/Km)
Zc2 (Ohm) 234.289
V2 (Km/s) 297800.0
Tabla 46

Y se tiene una matriz de transformación dada por:

 0.60 5517 64 - 0
 - 0 .001 8892 0 0.

 0.51 6406 31 0.
Tv  
 0.00 3762 19 0.
 0.60 5517 64 0.


 - 0 .001 8892 0 0.

55)

Tenemos

    Tv  1  Z fase Y fase  Tv  con Y   jwC 

fase fase 56)
  Realizar una simulación energizando una sola fase con una fuente tipo escalón
(11). Determinar la frecuencia real de la oscilación del voltaje al final de la línea
en la fase energizada. Una vez determinada la frecuencia real de oscilación del
voltaje, cambiar la fuente por una del tipo sinusoidal (tipo 14) y observar
nuevamente el voltaje al final de la línea, pero esta vez interesa observar el
resultado cuando el voltaje se estabilice. Justificar el resultado. ¿Cuál sería la
longitud de línea que para la anterior fuente no se presente un sobrevoltaje en
estado estacionario por encima del 10%? A partir de los resultados obtenidos
dar una posible justificación a la frecuencia comercial de 60 Hz.

Figura 80

De esta gráfica se tiene que la frecuencia de oscilación del voltaje al final de la línea
está dado por:

f  330 Hz (56)
Cambiemos ahora la fuente por una senoidal triásica de 330 Hz

Figura 81

De donde se tiene en estado estable un voltaje pico de:

V pico  6815.3 Kv  9.638 p.u (57)

Luego se tiene un sobrevoltaje de aproximadamente 10 veces el voltaje de

alimentación.
Una manera de impedir este sobrevoltaje es cambiar la frecuencia del transitorio de
voltaje al final de la línea.

[Link] DEL RAYO

3 PROCEDIMIENTO

Una nube cargada se puede representar circuitalmente como un capacitor, con una
determinada carga dependiendo de su geometría. Para efectos de cálculos una nube se
puede aproximar geométricamente a una placa circular paralela a la superficie del
suelo.
Para esta aproximación se utiliza la Expresión :

 0  A  0  R 2
C  
d d

Expresión 58 Representación circuital de una nube cargada

Donde,
C = Capacitancía de la nube (F)
 0 = 8.854x10-12 (F/m)
R = Radio de la nube (m)
D = Altura promedio de la nube al suelo (m)

En la Figura se puede ver la representación de la nube y del canal de descarga.

Figura 81 Representación de la nube y el canal de descarga

El canal de la descarga se comporta como un hilo delgado de aproximadamente 2 mm

de diámetro, pero debido al gran campo eléctrico que se origina, el efecto corona hace
que la onda de corriente de descarga vea un diámetro mayor, que para efectos de
cálculo se puede asumir de unos 20 cm. La inductancia del canal de descarga se puede
calcular con la Expresión .
 0  2d  0
L ln r  1  d  8  d
2  

Expresión 59 Inductancia del canal de descarga.

Donde,
L = Inductancia del canal de descarga (H)
 0 = Permeabilidad magnética del aire = 4  107  (H/m)
r = radio del canal de descarga (m)
d = altura promedio de la nube (m)

Para efectos prácticos se considerará:

Radio de la nube = 1 km

Altura de la nube sobre el suelo = 5 km
Radio del canal de descarga = 10 cm
Voltaje de carga de la nube = 50 MV
Resistencia de prueba = 50 Ω
Tiempo máximo de simulación = 100 mseg
Delta de tempo = 0.1 mseg

El primer paso es calcular el valor de la resistencia del canal de tal manera que se
obtenga lo más próximo a una forma de impulso de corriente 1.2/50 mseg. Esta
resistencia físicamente es de un valor alto y es la que asegura que finalmente un rayo
se pueda representar circuitalmente como una fuente de corriente.
Se debe obtener un valor inicial para R de tal manera que la respuesta de corriente del
circuito RLC serie, sea críticamente amortiguado. El valor que se obtenga para R es la
suma de la resistencia del canal de la descarga y la resistencia de prueba.

En el circuito se utiliza la opción de condensador con condición inicial y es la manera

más sencilla de implementar esta simulación.
 Figura 82 Circuito para la simulación de la corriente de rayo

Para encontrar esta resistencia se deben hacer varias simulaciones hasta obtener la
que aproximadamente tuviera las condiciones requeridas (forma de impulso de
corriente 1.2/50 mseg).

Los parámetros que se ingresaron al circuito de la Figura , son los siguientes:

13.1.1Condensador
El condensador de la nube fue calculado con la Expresión y la forma como se ingresan
los datos al ATP se pueden ver en la Figura 81.

Figura 81 Datos del Capacitor en el ATP

13.1.2 Inductancia
La inductancia del canal fue calculada con la Expresión y la forma como se ingresan
los datos al ATP se pueden ver en la Figura .
 Figura 84 Datos de la inductancia en el ATP

13.1.3 Resistencia de prueba

La resistencia de prueba, para nuestro caso se considera de 50 Ω y la forma como se
ingresan los datos al ATP se pueden ver en la Figura 82

Figura 82 Datos de la resistencia de prueba en el ATP

13.1.4Resistencia del canal

La resistencia del canal se debe variar hasta obtener la forma de onda tipo rayo, para
nuestro caso se variará de 9500 Ω a 11500  Ω y la forma como se ingresan los datos al
ATP se pueden ver en la Figura 82

13.1.5Suiche
El suiche se cerrará en t = 0 s y se abrirá en t = 1 s; en este elemento se medirá
corriente y la forma como se ingresan los datos al ATP se pueden ver en Figura 83.
 Figura 83 Datos del suiche en el ATP

13.1.6 ATP Settings

Finalmente se ingresan los parámetros de corrida del simulador, los cuales se pueden
ver en la Figura 84

Figura 84 Parámetro de corrida del simulador ATP

4 Simulación

Después de tener el circuito con los parámetros adecuados procedemos a simular,

variando la resistencia del canal cada 100 Ω.
Finalmente los datos obtenido en las simulaciones se pueden apreciar en la Tabla :

Tabla 47 Resultados de la simulación variando la resistencia del canal

Rcanal (Ω) Imaximo (A) I 50% (A)

9500 4908.9 2092.3
9600 4865.7 2090.3
9700 4817.3 2088.1
9800 4772.7 2085.7
9900 4728.8 2083.2
1000 4685.7 2080.5
10500 4294.8 2046.3
10750 4294.8 2046.3
11000 4294.8 2046.3
11250 4294.8 2046.3
11500 3961.7 2002.9

Con los datos de la Tabla , se puede decir que la resistencia del canal mas apropiada
es de 10500 Ω, ya que la corriente a los 5 mseg tiende a ser mas cercana a la
corriente pico.

Figura 85 Forma de onda de la corriente para la resistencia del canal de

10500 Ω.

De la Figura 85 se observa un frente de onda con un tiempo de 4.2 µS, y una corriente

pico de 4294.8 A, además un tiempo de cola de 50 µS donde se observa
aproximadamente el 50 % del valor de la corriente pico (2046.3 A), razón por la cual
se considera que el valor de resistencia de canal que se utilizo para llegar a este
resultado 10500 Ω fue adecuado. Luego el valor de RTOTAL es de 10550 Ω.
El programa EMTP/ATP tiene una fuente tipo rayo, ya establecida que se comporta de
igual manera que el circuito de la Figura , para constatar esto se mostrará la fuente y
la forma de onda de corriente que esta produce.

En la Figura 86, se mostrará como obtener la fuente 13 del ATP, que simula una fuente
tipo rayo.

Figura 86 Forma de obtener fuente tipo 13

La forma de ingresar los datos a la fuente tipo 13, para que simula una fuente tipo
rayo, se muestra en la Figura 7.
 Figura 7 Datos de la Fuente tipo 13 para simular una Fuente tipo rayo

La fuente tiene una corriente pico de 4300 A con un tiempo de frente de onda de
1.2 µS, tiempo de cola de 50 µS donde se da el 50 % del valor de la corriente pico
(2150 A), lo anterior se puede observar gráficamente en la Figura 8.

5000
[A]
4000

3000

2000

1000

0
0 10 20 30 40 [us] 50
(f ile Ray o230.pl4; x-v ar t) c:FTE -X

Figura 8 Forma de onda tipo rayo simulada con la fuente tipo 13

 14 MODELAMIENTO DE LÍNEAS

5 PROCEDIMIENTO

Para este modelamiento se va a considerar una línea de 230 kV con dos cables de
guarda tal como se observa en la Figura .

Figura 92 Estructura típica de 230 kV.

 Datos del conductor de fase

Rdc = 0.05086 Ω/km
Diámetro = 30.98 mm

 Datos del cable de guarda

Rdc = 1.9014 Ω/km
Diámetro = 9.78 mm

La modelación de la estructura metálica se puede hacer tal como aparece en la Figura

9.
 Figura 9 Modelación de la estructura metálica

En la modelación de la estructura se consideran los tramos más importantes. Si se

desea se puede detallar aún más incluyendo por ejemplo las crucetas en las fases,
pero la precisión que se logra al detallar no es importante.
Cada uno de los elementos metálicos se puede modelar como un tramo de línea ideal
sin pérdidas con una impedancia característica (estructuras de forma cónica), como se
muestra en Expresión :

 2( h 2  r 2 ) 
Z T  30 ln  
 r2 

Expresión 60 Impedancia característica de la torre

Donde,

h es la altura de la torre (m)

r el radio en la base, suponiendo un cono equivalente (m)
Como valor típico para ZT se puede asumir un valor de 200 Ω.
La velocidad de propagación es aproximadamente la de la luz.
La torre se divide de acuerdo a los tramos que se necesiten con la misma impedancia
característica de 200 Ω y de acuerdo a la distancia que se considere. El modelo a
considerar es el de CLARKE para el caso monofásico y sin pérdidas, tal como se
observa en la Figura .

Figura 94 Modelo de Clarke para tramo de red monofásica

La forma de entrar los datos al modelo Clarke para tramo de red monofásico se puede
apreciar en la Figura 10.
 Figura 10 Forma de entrar los datos al modelo Clarke.

Los conductores de fase conjuntamente con los cables de guarda se modelan como 8
fases acopladas entre si, mediante el modelo de [Link], ver Figura 11, a una
frecuencia que puede ser de 100 kHz (lo que se denomina el modelo Bergerón en el
ATPDRAW).
Asumir un valor del vano de 0.5 km.
Para obtener el correspondiente modelo se puede utilizar la opción LINE/CABLE de la
caja de herramientas del ATPDRAW. El modelo a utilizar es el de onda viajera
(Bergerón) para línea no transpuesta, para ingresar estos parámetros ver la Figura 12.
 Figura 11 Modelo de [Link].

Figura 12 Parametros de entrada del modelo de [Link].

Para configurar la geometría de los conductores, se ingresan el radio del conductor, la

resistncia DC y las distancias entre conductores y la altura del conductor al suelo en la
ventana de datos del modelo de [Link] como se puede ver en la Figura 13.
 Figura 13 Datos de entrada del modelo de [Link].

Para verificar si la disposición de los conductores es correcta abre la ventana view

model, como se puede ver en la Figura 14.
 Figura 14 Vista de los conductores en el modelo de [Link].

La cadena de aisladores se modela como un suiche tipo “gap” que es controlado por
voltaje (Figura 15).

Figura 15 Suiche controlado por voltaje (cadena de aisladores)

La forma de ingresar los datos al suiche controlado por voltaje se puede ver en la
Figura 16, donde se puede ver que la cadena de aisladores tiene un voltaje de flameo
de 5000000 V y el tiempo de respuesta es de 1E-8 s.

Figura 16 Datos de entrada de la cadena de aisladores

La resistencia de puesta a tierra se modela como una del tipo concentrado, con el valor
correspondiente (20 ohm en forma típica)
Para retomar la forma de onda tipo rayo, simularemos una descarga atmosférica en el
cable de guarda de la línea, esta corriente de rayo se simulará con una fuente rampa
de doble pendiente tipo 13 (ver Figura 17)

Figura 17 Datos de entrada de la Fuente tipo 13, simulando descarga

atmosferica de 30kA.

De manera adicional se puede simular el comportamiento del aislamiento de acuerdo a

la formulación que propone la IEEE, donde la resistencia que ofrece el aislamiento
depende del tiempo (ver Expresión 1).

 0. 71 
Vflameo   0. 4  0. 75   d MV
 t 

Expresión 1 Comportamiento del aislamiento de acuerdo con la IEEE.

Donde,

t es el tiempo (useg)
d es la distancia longitudinal de la cadena (m)

Esta respuesta del aislamiento se ha simulado mediante la opción de los TACS del ATP
(Ver Figura 19 y Figura 19) que permite implementar fórmulas matemáticas y bloques
de control, la formula que se ingresa al TAC se puede ver en la Expresión 2

 0. 71 
Vflameo   0. 4  # aisladores  long por aislador ( m )  106 V
 ( t * 106 ) 0. 75 

Expresión 2 Formula ingresada a los TACS

Figura 18 TACS
 Figura 19 Datos de entrada de los TACS

Finalmente el sistema a simular es (ver Figura 20):

Figura 20 Sistema completo a simular

Se puede observar que hay dos vanos de 500 m a lado y lado de la torre fallada. Es
suficiente esto, porque las puestas a tierra que más influyen son las que quedan más
cerca de la torre impactada por el rayo.

El hecho que en forma real sean muchos vanos se simula con una línea de longitud
infinita (para que no regresen las reflexiones), o mediante elementos resistivos
acoplados con valores de resistencia igual a las impedancias características de la línea.
Para nuestro caso se puede considerar tramos de línea de una distancia
suficientemente grande para que las reflexiones regresen fuera del tiempo de análisis,
lo cual equivale en la práctica a considerar línea infinita.

Para el sistema se ha considerado tramos de línea de 2.5 km cuyas reflexiones

regresarían en un tiempo aprox. de 16.7 µseg.

Hay que tener especial cuidado con el delta de tiempo de simulación, ya que debe ser
más pequeño que el tiempo de viaje más corto (_t < 4.5m/3x108 m/s = 15 nseg).
El tiempo de simulación, Tmáx, debe estar entre entre 10 y 20 mseg para que no sea
muy pesada la simulación y adicionalmente en este tiempo los peores sobrevoltajes
sobre la cadena ya se han presentado (ver Figura 21).

Figura 21 Settings de la simulación de una descarga atmosférica en un cable

de guarda.

6 SIMULACIÓN

Los voltajes en las cadenas de aisladores se muestran en la Figura 22

 100
[kV]

-100

-300

-500

-700
0 3 6 9 12 [us] 15
(f ile Ray o230.pl4; x-v ar t) v :CTO1_B-MEDIO v :CTO1_C-INF
v :CTO2_A-SUP v :CTO2_B-MEDIO v :CTO2_C-INF
v :CTO1_A-SUP

Figura 22 Voltajes en las cadenas de aisladores

En la Figura 22 se aprecia cómo, ante una descarga de 3kA, con forma de onda
1.2/50s, los voltajes sobre las cadenas de aisladores se elevan, en el punto más
crítico, hasta valores alrededor de 683.42kV, para una resistencia de puesta a tierra de
25.
Los voltajes que aparecen sobre las fases son inducidos por la circulación de corriente
a través del cable de guarda y son negativos porque el voltaje de referencia (el neutro)
se subió. El nivel del sobrevoltaje depende de la distancia (para el caso de descarga
sobre el neutro) entre el guarda y los conductores de fase, de modo que mientras
mayor sea la diferencia de voltaje entre ellos, mayor es el sobrevoltaje que se
presenta. Es por lo anterior que, a pesar de que la fase C es la que más retirada se
encuentra del guarda, es en ella donde se aprecia el mayor sobrevoltaje.
Además se puede apreciar que la cadena de aisladores protegen bien la torre ya que,
los sobrevoltajes en la cadena son mucho menores a los voltajes que puede soportar
(ver Figura 23).
 4
*10 6

-2

-4
0 3 6 9 12 [us] 15
(f ile Ray o230.pl4; x-v ar t) v :CTO1_A v :CTO1_B v :CTO1_C
v :CTO2_A v :CTO2_B v :CTO2_C t: GAP15N t: GAP15P

Figura 23 Voltaje de soporte de la cadena de aisladores.

BIBLIOGRAFIA

 CARDONA, Leonardo de Jesús. Guía de prácticas Fenómenos

Transitorios Electromagnéticos. Universidad Nacional de Colombia.
2004. Pag: 25-39.
 GAVIRIA, Alvaro. Notas de clase de Teoría Electromagnética.2001.
Universidad Nacional de Colombia.

 THEORY BOOK. Cap. 4. Pag: 4_1 – 4_86.