UNIVERSIDAD NACIONAL DE COLOMBIA

SEDE MEDELLiN
FACIULTAID DE MIINAS  
ESCUELA DE INGENIERiA ELECTRICA Y
MECANICA

GUiA DE pRACTICAS
FEN6MENOS
TRANSITORIOS
ELECTROMAGNETICOS

ＮＭ ｾ
a
t 11>11VERS \l)Al) N
ｾｦ｡ｉＧ＠
ｾ

Sill»" ＬＮｴＧ＠
ＧｊＱｆＢ＠

Ｍ
"C\ONAL DE COLOMBIJ\
ﾷ Ｑ ｾ ﾷ ｌＮｾ＠ -
DEP'TO . DE BIBLlOTECAS
RlBUOTECA MI NAS

LEON'A RDO CARDONA C.
Profesor asistente

Medellin
Juni,Q del 2004
 GUIA DE pRACTICAS
FENOMENOS TRANSITORIOS
ELECTROMAGNETICOS

ｾ＠
Ｍ
ｾ ｬＮ＠
:'  ｾ＠
ｾ＠
LＱｾＢＧｬｭａ ｮ＠ N ACIO"iAL DE COLOMBIA 
Sl PP' ',II  ,'- ,1.1\'11

DE PTO. DE BTBLlOT£ CAS

RIBUOTECA M[NAS

LEONARDO CARDONA C.
Profesar asistente

Trabajo presentado para

Promoci6n a la Categoria
de Profesor Asociado

UNIVERSIDAD NACIONAL DE COLOMBIA

SEDE MEDELLiN
FACULTAD DE MINAS
ESCUELA DE INGENIERiA ELECTRICA Y
MECANICA
Medellin. Junia del 2004
:f 
;, 1. )/1:2 J
k
CIt--
11 ;'
J:.Iii
". 
ｾ＠
..,;i-::?",.;IJ
\.l1'l,VEltSH) . .O 1'l -"C\Ql'l .... L 010.
=­­­ _  Sl,l}\;: .... r.P . \.\.ttool
ｃ ｏ ｌ ｏ ｾ ｉ＾ ｴａＮ＠
CONTENIDO

DE\"TO . DE B1BLI01'E Ci\S

B \ SU 0 1' ECi\ M,H-lJ\S pag

INTRODUCCI6N  17
1.  INTRODUCCI6N A LOS  FEN6MENOS TRANSITORIOS  18  

1.1  INTRODUCCI6N  18  

1.2  CLASIFICACI6N  DE  LOS  TRANSITORIOS  19  

1.3  OBJETIVOS  20
2.  INTRODUCCI6N AL  PROGRAMA EMTP/ATP  21  
2.1  HISTORIA  21  
2.2  ESTUDIOS A REALIZAR  CON  EL  PROGRAMA EMTP/ATP  22
2.3  ESTRUCTURA GENERAL DE  UN ARCH/VO  PARA LA  ENTRADA DE  
DATOSALATP  24  
2.4  REGLAS GENERALES PARA LA GENERACI6N DEL ARCHIVO  DE  
ENTRADA DE  DATOS  25  
2.5  ESTRUCTURA GENERAL DEL ARCHIVO  DE  DATOS PARA UN  CASO  
SIMPLE DE  SIMULACI6N TRANSITORIA  0 DE  ESTADO  ESTACIONARIO DE  
UNA RED  ELECTRICA CON  RAMAS  SIMPLES  RLC  26  
3.  pRACTICA N01 .  CASO  SIMPLE  DE  SIMULACI6N  DE  UN  CIRCUITO RLC  
SERlE  29  
3.1  OBJETIVO  29  
3.2  DESCRIPCI6N  29  
3.3  PROCEDIMIENTO  31  
3.4  ANALISIS  DE  RESULTADOS  PARA EL CIRCUITO RLC  39  
3.5  COMPLEMENTACI6N  39  
4.  pRACTICA N02.  ONDAS VIAJERAS  43  
4.1  OBJETIVO  43  
4.2  MARCO TE6RICO  43  

L-{ 3 SB ?J 0
4.3  SIMULACIONES  45  
4.3.1  Energizacion de una linea  monofasica sin  perdidas,  en  vado,  con  una  
fuente  impulso de voltaje.  45  
4.3.2  Energizacion de  linea  monofasica sin  perdidas,  en  vado  46  
4.3.3  Energizacion de  linea  monofasica sin  perdidas,  con carga  49  
4.3.4  Energizacion de  linea monofasica sin  perdidas,  con  resistencia de  
preinsercion  49  
4.3.5  Obtencion  de  modelos Lee  y Clarke con el  subprograma del ATP,  LINE  
CONSTANTS.  50  
5.  pRACTICA N°3.  FEN6MENOS TRANSITORIOS ORIGINADOS  POR  UNA  
CORRIENTE  DE  RAYO  58  
5.1  INTRODUCCI6N  58  
5.2  MODELAMIENTO DEL RAYO  59  
5.3  SISTEMA DE TRANSMISI6N A ESTUDIAR  62  
5.4  ASPECTOS  DE  MODELACI6N DE  ESTRUCTURA Y CONDUCTORES  63  
5.5  PROCEDIMIENTO Y PREGUNTAS  69  
6.  pRACTICA N°4. SOBREVOL TAJES TRANSITORIOS  EN  SISTEMAS  DE  
DISTRIBUCI6N,  DEBIDO A MANIOBRA DE  CAPACITORES  73  
6.1  ENERGlZACI6N DE  UN  BANCO  DE  CAPACITORES  73  
6.2 MAGNIFICACI6N DE VOL TAJE  POR  SUICHEO DE  BANCOS DE  
CONDENSADORES  76  
6.3  SISTEMA A ESTUDIAR  78  
6.3.1  Equivalente de cortocircuito en  la  barra de  13.2 kV  79  
6.3.2  Cable aislado de 350  KCM  para  la  salida del  circuito de  13.2 kV  79  
6.3.3  Red  aerea de  13.2 kV  79  
6.4  SIMULACIONES  80  
7.  pRACTICA N°  5.  SOBREVOL TAJES A  BAJA FRECUENCIA DEBIDOS A  
FALLAS ASIMETRICAS  83  
7.1  SISTEMA A ESTUDIAR  83  
7.2  DATOS  DE  LOS  ELEMENTOS DE  RED  84  
7.3  CALCULAR Y SIMULAR  87  
8.  pRACTICA N°6.  CARGAS  NO  LINEALES ­ PUENTE  RECTIFICADOR  
DE  SEIS  PULSOS  . 88  
8.1  INTRODUCCION  88  
8.2  DESCRIPCION  DEL  SISTEMA A ANALIZAR  89  
8.3  OBJETIVOS  90  
8.4  PUENTE  RECTIFICADOR DE  SEIS PULSOS  91  
8.5  PUENTE RECTIFICADOR CONTROLADO  DE  SEIS  PULSOS  93  
8.6  MODELO ARMONICO  98  
8.7 DISTORSION  DEBIDA  AL  PUENTE  Y  PRESENCIA  DE  
TRANSFORMADOR  102  
8.8  DISTORSION  CON  LA CONEXION  DE  BANCO CAPACITIVO  103  
8.9  NOTAS  ADICIONALES  104  
8.10  CONTENIDO  DEL INFORME  105  
BIBLIOGRAFIA  107  
ANEXOS  108  
 LlSTA DE FIGURAS

pag 

Figura  1.  Clasificaci6n de los transitorios de acuerdo con  la  frecuencia  19  

Figura  2.  Programas de soporte que interactUan con  el  EMTP/ATP  23  
Figura  3.  Formato general  para  el  programa  EMTP/ATP  28  
Figura 4.  Circuito  RLC  serie  30  
Figura  5.  Circuito RLC serie en Atpdraw  31  
Figura 6.  Caja de herramientas de la interfaz Atpdraw  32  
Figura 7.  Ventanas de datos para  rama  RLC,  suiche y fuente escal6n  32  
Figura 8.  Ventana de datos para  los "settings"  33  
Figura 9.  Caja de herramientas para  datos enteros  33  
Figura  10 .  Persiana "ATP" para  realizar el  "Make  File"  34  
Figura  11 .  Corriente en  el  circuito  RLC  serie  38  
Figura  12.  Voltajes  y corrientes en  el  circuito  RLC  serie  38  
Figura  13.  Circuito con  una  cascada de elementos RLC  40  
Figura  14.  Resultado esperado  de  la  simulaci6n con  elementos  
semidistribuidos  41  
Figura  15 .  LInea energizada con  un impulso de voltaje  45  
Figura  16.  Ventana de datos,  fuente tipo impulso de voltaje  46  
Figura  17.  Modelo de Clarke monofasico para onda viajera  46  
Figura  18.  Ventana con  la  ayuda para entrada de datos del  modelo de Clarke  47  
Figura  19.  LInea  ideal energizada con  fuente escal6n de voltaje  47  
Figura  20 .  LInea ideal energizada con  carga  49  
Figura  21 .  LInea ideal con  resistencia de  preinserci6n  49  
Figura 22.  Configuraci6n linea de 500  kV  54  
Figura  23 .  Ventana de datos para el  modelo de Clarke trifasico  55  
Figura 24.  Representaci6n de la  nube y canal de descarga  59  
Figura  25 .  Circuito para simular corriente de rayo  61  
Figura 26 .  Forma esperada de  la  corriente de  rayo  61  
Figura 27. Estructura Hpica  de 230 kV  62  
Figura  28.  Trayectoria  del  rayo  y  diagrama  para  modelaci6n  con  
ATPDRAW­ATP  63  
Figura 29 .  Modelo de  Clarke para tramo de  red  monofasica  64  
Figura 30.  Ventana de datos para  obtener modelo K.C .LEE  para  la  linea  65  
Figura 31 .  Modelo tipo "gap" para  la cadena de aisladores  65  
Figura 32 .  Dispositivo 51  de  los TACS  67  
Figura 33.  Descripci6n del  ､ｩｳｰｯｴｶｾ＠ 51  67  
Figura  34 .  Sistema completo  para  simulaci6n  68  
Figura  35. Diagrama unifilar de un  circuito de  distribuci6n con capacitor para  
suicheo  74  
Figura  36.  Circuito equivalente para energizaci6n del  banco  ｣｡ｰｩｴｶｾ＠ sin  carga  
en  red  primaria  74  
Figura  37 .  Voltaje observado en  la  barra donde esta  ubicado el  banco  75  
Figura  38 .  Voltaje observado en  un  punto del  circuito primario entre el  banco y  
la  fuente de  potencia  76  
Figura 39 .  Diagrama  unifilar circuito de distribuci6n con  capacitor para suicheo  77  
Figura 40 .  Circuito equivalente para  estudio de  magnificaci6n  de voltaje  77  
Figura  41.  Configuraci6n  de estructura de  13.2 kV,  neutro superior  79  
Figura 42 . Esquema circuital en  Atpdraw para  energizaci6n de  un  banco  
｣｡ｰｩｴｶｾ＠ 80
Figura 43 .  Resultados esperados durante la  energizaci6n  81
Figura  44.  Esquema  circuital  en  Atpdraw  para  estudio  de  magnificaci6n  de  
voltaje  82  
Figura 45 .  Diagrama unifilar de  un  sistema de distribuci6n  83  
Figura  46.  Configuraci6n de estructura de  13.2 kV, neutro inferior  85
Figura 47 .  Configuraci6n  red  de 44  kV  86  
Figura 48.  Sistema industrial a analizar  89
Figura 49.  Puente  rectificador de seis  pulsos  91
Figura 50 .  Suiche tipo  11  0  diodo  91  
Figura 51 .  Puente rectificador para simulacion  92  
Figura  52.  Puente rectificador en  Atpdraw y variables de interes  93  
Figura  53 .  Puente rectificador controlado  94  
Figura 54 .  Estrategia de control del  puente  94  
Figura 56 .  Generacion  de un  nuevo elemento con  MODELS  99  
Figura 57 .  Incorporacion del  modelo ARMO dentro de la librerfa de  MODELS  99  
Figura 58 .  Definicion de datos y nodos del  modelo ARMO  100  
Figura  59. Utilizacion de  modelo ARMO en  descomposicion de corriente del  
puente  101  
Figura 60.  Resultado de la  descomposicion armonica de  la corriente  101  
Figura 61 .  Distorsion  armonica del  voltaje antes y despues del  transformador  102  
Figura 62 .  Diagrama en  Atpdraw incluyendo el  transformador potencia  103  
Figura 63 . Formas de onda del voltaje en  baja  con  banco y sin  banco  ｣｡ｰｩｴｶｾ＠ 104  
Figura 64.  Sistema base implementado en  ATPDRAW  105  
 LISTA DE TABLAS

pag 

Tabla  1.  Archivo plano generado por la  interfaz  Atpdraw  34  

Tabla 2.  Archivo de resultados  36

Tabla 3.  Archivo plano generado por la  interfaz  Atpdraw  70

Tabla 4 .  Estrategia de control  para el  puente en  c6digo TACS  95

Tabla  5.  Modelo completo del puente controlado en c6digo TACS  96

Tabla 6.  Modelo ARMO en  lenguaje MODELS  98

 LlSTA DE ANEXOS

pag 

ANEXO A. Determinaci6n de equivalente de Thevenin trifasico a partir de los 
niveles de  cortocircuito  108 

ANEXO B.  Soluci6n de la  ecuaci6n de onda viajera  110 

 LIST A  DE T ABLAS

pag 

Tabla  1.  Archivo plano generado por la  interfaz  Atpdraw  34

Tabla  2.  Archivo de  resultados  36

Tabla 3.  Archivo plano generado por la interfaz  Atpdraw  70

Tabla 4.  Estrategia de control  para el  puente en  c6digo TACS  95

Tabla  5.  Modelo completo del  puente control ado en  c6digo TACS  96

Tabla 6.  Modelo ARMO en  lenguaje MODELS  98  

 LlSTA DE ANEXOS

pag 

ANEXO  A. Determinaci6n de equivalente de  Thevenin trifasico a partir de los 

niveles de cortocircuito  108 

ANEXO  B.  Soluci6n  de la  ecuaci6n de onda  viajera  110 

 LlSTA DE SiMBOLOS

C   Capacitancia . 

C1 Capacitancia de secuencia positiva de la  linea . 

Co Capacitancia de secuencia cera de la linea .  

De  Distancia  media  geometrica  entre  las  cargas  de  los  conductores  y sus  
respectivas imagenes.  

DELTAT  Delta de tiempo de simulaci6n  en  el  ATP.  

Dgg  Distancia media geometrica entre los cables de guarda.  

Dguardas  Distancia entre los cables de  guarda .  

DMG  Distancia  media geometrica entre fases .  

DMG fg  Distancia  media  geometrica entre  los conductores de fase  y  los cables  

de guarda .  

DMGfases­imagenes  Distancia  media  geometrica  entre  los  conductores  de fase  y  sus 

imagenes. 

DMGfases­imagenes guardas  Distancia  media geometrica entre los conductores de fase  y 

las imagenes de  los guardas. 

DMGguardas­imagenes guardas  Distancia media geometrica entre los cables de guarda  y  

sus imagenes.  

d  Distancia.  

C Permitividad electrica del  medio.  

Co   Permitividad  electrica del vacfo . 

h  

FP Factor de  potencia . 

f  Frecuencia de  la  corriente que circula  por el  conductor de fase . 

fo  Frecuencia de  la  oscilaci6n  senoidal. 

({J Diametro del conductor. 

G  Conductancia electrica . 

h  Altura  de  la  torre. 

Corriente electrica. 

Icc1 ｾ＠ Corriente de cortocircuito  monofasica. 

ｉ｣Ｓｾ＠ Corriente de cortocircuito trifasica. 

L  Inductancia. 

'A Factor de amortiguamiento de  sobretensi6n. 

).1  Permeabilidad  magnetica del  medio. 

).10  Permeabilidad  magnetica  del  vado. 

Ng  Densidad  de  rayos  en  la  zona  (Numero de descargas/km 2/ario) 

Q  Sfmbolo del  Ohmio. 

1t  Pi. 

pu  Valor en  por unidad 

R  Resistencia electrica 

Rac  Resistencia AC  del  cable.  ,t

Rc  Resistencia  de carga . 

Rdc  Resistencia  DC  del cable. ' 

RMG'  Radio  medio geometrico del  haz de conductores.  Para una  fase 
compuesta  por un  solo conductor,  equivale al  radio geometrico del conductor. 

Ro  Resistencia de secuencia cero. 

R1  Resistencia de secuencia positiva. 

rf  Radio de los conductores de fase . 

rg  Radio de los cables de guarda. 

r1  Coeficiente de reflexi6n  al  comienzo de la  linea . 

r2  Coeficiente de  reflexi6n  al final  de la  linea.  

p  Resistividad  del  suelo.  

Tmax  Tiempo de simulaci6n  en  el  ATP.  

Tf  Tiempo de frente de la  onda  

t  Tiempo  

V  Voltaje electrico  

V PK  Valor pico del  voltaje fase­tierra  de  la  fuente  

v  Voltaje electrico  0  velocidad  de propagaci6n  

Vf  Voltaje de la  fuente.  

w  Frecuencia angular de la  corriente.  

X  Reactancia  inductiva .  

Xo  Reactancia de  secuencia cero.  

X1  Reactancia de  secuencia positiva .  

(X/R) 3'1'  Relaci6n  XIR que ve  la  corriente de corto trifasica .  

(X/R) Iq> Relaci6n  XIR que ve  la corriente de corto monofasica . 

xcc  Reactancia inductiva de cortocircuito de un  transformador.  

x  Distancia .  

lAB Impedancia  de  dispersion  medida  en  el  lado  de  alta  tension  del  
transformador, con  el  lado de baja cortocircuitado  y  el  lado de  media abierto.  ­:-
I  • 

lAM Impedancia  de  dispersion  medida  en  el  lado  de  alta  tension  del 
transformador,  con  el  lado de media cortocircuitado  y el  lado de baja abierto.  ­ . 

1MB Impedancia  de  dispersion  medida  en  el  lado  de  media  tension  del 
transformador, con  el  lado de baja cortocircuitado  y el  lado de alta abierto. 

lc Impedancia caracterfstica de sobretension . 

lfg Impedancia mutua entre los cables de guarda  y  los de fase . 

19 Impedancia propia de  los cables de  guarda. 

199 Impedancia mutua entre los  cables de guarda. 

lm Impedancia  mutua entre los conductores de fase . 

ls Impedancia propia del  conductor de  fase . 

IT Impedancia de  puesta a tierra 

l1 Impedancia de secuencia  positiva de  la  linea  0 impedancia de la fuente 

de  voltaje . 

l2 Impedancia en  el  extremo de la  linea. 

 RESUMEN

EI  presente  trabajo  corresponde  a  la  documentaci6n  de  10 que  ha  sido  la 
experiencia  en  el  curso  de  Fen6menos  Transitorios  de  la  Linea  de  Profundizaci6n 
en  Coordinaci6n  de Aislamiento.  EI  trabajo es  una  gufa  detallada  de  las  practicas 
realizadas con  el  programa ATP. 

Se  hace  una  introducci6n  en  10 que  es  el  manejo  basico  del  ambiente  de 
simulaci6n  de  fen6menos  transitorios  en  sistemas  electricos  de  potencia  ATP-
ATPDRAW. En esta introducci6n se revisan las posibles aplicaciones que tiene el
programa ATP en anal isis electricos.

Se realiza una primera practica que tiene que ver con el circuito basico RLC serie
pero con una clara orientaci6n hacia los fen6menos transitorios que aparecen en
Ifneas de transmisi6n , especfficamente involucra conceptos como el de
impedancia caracterfstica .

La practica siguiente tiene que ver con el concepto de ondas viajeras en redes de
transporte de energfa electrica . Se adiciona un anexo donde esta solucionada
paso a paso la ecuaci6n de onda viajera .

Se detalla una practica sobre descargas atmosfericas en Ifneas aereas de alta

tensi6n . En esta practica se manejan desde las caracterfsticas ffsicas del rayo,
hasta el concepto de coordinaci6n de aislamiento teniendo en cuenta la longitud
de la cadena de aisladores .

Se presenta una practica sobre suicheo de condensadores en un sistema de

distribuci6n, manejando dos situaciones que pueden generar problemas de calidad
de la potencia, por los sobrevoltajes que se presentan . La primera situaci6n que
se considera, es la del suicheo de un banco situado en un circuito de distribuci6n .
La segunda situaci6n que se presenta es la de magnificaci6n de voltaje debido a la
interacci6n entre un banco situ ado en el circuito de distribuci6n y un banco
｣｡ｰｩｴｶｾ＠ situado en lado de baja tensi6n de un usuario.
Los  sobrevoltajes temporales  0  denominados de  baja  frecuencia  se  ilustran con  el 
caso  de fallas  asimetricas en  un  sistema  de  distribuci6n. Los  sobrevoltajes  que  se 
presentan se orientan hacia una adecuada especificaci6n  de  pararrayos. 

Una  practica final  se  refiere  al  comportamiento  no  lineal  de  equipos  en  un  sistema 
de  potencia .  EI  caso  especffico  es  el  debido  a  la  presencia  de  elementos  de 
electr6nica  de  potencia .  Se  considera  el  caso  de  un  rectificador  trifasico  de  gran 
potencia  en  un  sistema  industrial  donde  genera  problemas  de  distorsi6n  arm6nica 
en  el  voltaje de suministro. 
 INTRODUCCION

Un  sistema  electrico  de  potencia  esta  sometido a una  serie de  perturbaciones que 

alteran  su  estado  normal  de  operacion .  EI  paso  de  un  estado  a  otro  ocasionado 
por  una  perturbacion  se  hace  en  forma  gradual, ya  que  las  variables  ffsicas  como 
voltajes  y corrientes  no  pueden  tener cambios  bruscos  debido  a  que  las  leyes  de 
la  Ffsica  no  10 permiten.  Todo  el  proceso de  cambio de  un  estado  de  operacion  a 
otro  es  10 que  se  pretende  estudiar  en  un  curso  de  Fenomenos  Transitorios  y 
corresponde a 10 que se  llama el  regimen  transitorio. 

EI  estudio  del  estado  transitorio  debe  hacerse  en  forma  cuidadosa  ya  que  las 
variables  fisicas  pueden  lIegar  a  tomar  valores  extremos  y  exigir  a  los  equipos, 
hasta  ocasionar  el  deterioro  de  los  mismos  0 la  interrupci6n  del  suministro  de 
energfa .  EI  conocimiento  de  los  estados  transitorios  permitira  tomar  decisiones 
para  proteger  adecuadamente  los  equipos,  10 mismo  que  el  aseguramiento  de  la 
calidad  de la  potencia electrica , con  unos estandares de calidad  adecuados . 

Este  texto  ha  sido  preparado  especialmente  para  orientar el  curso  de  Transitorios 
Electromagneticos  en  Sistemas  Electricos.  Tambien  puede  ser  utilizado  en  los 
temas  de  coordinacion  de  aislamiento.  EI  trabajo  es  producto  de  una  extensa 
labor de  trabajar en  el  tema  de  coordinaci6n  de  aislamiento  y con  el  programa  de 
simulaci6n  de  transitorios  electromagneticos ATP/EMTP  desde el  ario  1988, fecha 
de  lIegada de dicho programa a la  Universidad  Nacional, Sede Medellin . 

EI  trabajo  se  dividio  en  ocho  capitulos  donde  se  trabaja  en  forma  separada  seis 
temas  con  una  propuesta  de  practica  de  simulacion  para  cada  uno  de  ellos.  En 
los  tres  primeros  capitulos  se  hace  una  introducci6n  al  tema  del  fenomeno 
transitorio  y al  programa  ATP  en  su  forma  basica  de  utilizaci6n. EI  cuarto  capitulo 
se  refiere  al  tema  de  las  ondas  viajeras.  Este  capitulo  se  complementa  con  un 
anexo  sobre  la  soluci6n  de  la  ecuaci6n  de  onda .  En  el  capitulo  cinco  se  trata  el 
tema  de  las  descargas atmosfericas  y su  impacto en  una  linea  de alta  tension .  EI 
capitulo  seis  esta  dedicado  al  tema  de  la  energizacion  de  bancos  de  capacitares 
en  un  sistema  de  distribucion  y su  impacto  en  la  calidad  de  la  potencia.  EI  tema 
de  los  sobrevoltajes  temporales  debidos  a  fallas  asimetricas  en  sistemas  de 
distribucion  es  tratado  en  el  capitulo  siete .  EI  ultimo capitulo  es  dedicado  al  tema 
de los armonicos en  un  sistema  industria\. 

17
 1. INTRODUCCION A LOS FENOMENOS TRANSITORIOS

1.1 INTRODUCCION

Los  transitorios  en  un  sistema  de  potencia  en  forma  general  son  de  tipo  eltktrico, 
mecanico  y  termico.  Los  transitorios  mecanicos  y  termicos,  a  pesar  de  ser  muy 
importantes,  son  mucho  mas  lentos  que  los  de  tipo  electrico.  Su  frecuencia  esta 
por  debajo  de  la  sincronica  (60  Hz) .  Los  transitorios  electricos  son  muy  rapidos  y 
constituyen  el  motivo de este trabajo. 

EI  sistema  electrico,  desde el  punto  de  vista  circuital,  esta  formado  por tres  clases 
de elementos:  resistor, inductor y capacitor. 

Cualquier  componente  de  un  SEp1  estara  formado  por  la  combinacion  de 
elementos  RLC.  En  estado  estacionario  un  parametro,  para  un  determinado 
equipo,  sera  dominante  sobre  los  otros,  pero  en  estado  transitorio  hay 
dependencia  de  las  caracterfsticas  y  tipo  de  fenomeno  transitorio .  En  un 
transformador  de  voltaje,  por  ejemplo,  es  claro  que  la  componente  inductiva  es 
dominante  en  un  estado  normal  de  operacion  a 60  Hz  y  para  eventos  transitorios 
de  una  frecuencia  relativamente  baja,  pero  se  sabe  que  en  eventos  cuya 
frecuencia  esta  por  encima  de  unos  5  kHz,  el  efecto  ｣｡ｰｩｴｶｾ＠ empieza  a  ser 
importante  y el  transformador  pierde su  precision. 

Otra  caracterfstica  importante  es  la  forma  como  los  efectos  RLC  se  encuentran 
presentes  en  un  equipo .  Se  pueden  encontrar  de  manera  concentrada  0 
distribuida.  En  una  linea,  por  ejemplo,  los  parametros  RLC  se  encuentran 
distribuidos a 10  largo de  la  misma,  mientras que en  un  reactor para  compensacion 
de  Ifneas,  el  parametro  dominante  inductivo se  encuentra concentrado.  EI  modelo 
para  caracterizar estos elementos depende de  si  se  puede  representar de  manera 
concentrada  0  distribuida.  Estos  parametros  adicionalmente  pueden  depender de 
otra  variable,  como  es  el  caso  de  la  inductancia,  que  representa  la  magnetizacion 
de  un  transformador,  la  cual  depende  del  voltaje  aplicado . 

1 SEP : Sistema Electrico de  Potencia 

18 
1.2 CLASIFICACION DE LOS TRANSITORIOS

Los  transitorios  que  aparecen  en  un  sistema  de  potencia  se  pueden  clasificar 
desde  diferentes  puntos  de  vista .  Uno  obedece  a  los  tipos  de  fen6menos  ffsicos 
involucrados;  para  este  caSO,  los  fen6menos  transitorios  se  clasifican  en  dos 
categorfas: 

•   Transitorios  que  resultan  de  la  interacci6n  entre  las  energfas  almacenadas  en 
capacitores  y en  inductores. 
•   Transitorios que  resultan  de  la  interacci6n  de  la energia mecimica almacenada 
en  las  partes  rot6ricas  de  las  maquinas  y  la  energfa  almacenada  en  los 
circuitos. 

EI  interes en  estas  notas son  los transitorios electromagneticos. 

En  la  figura  siguiente  se  ilustra  esta  clasificaci6n  de  acuerdo  con  la  frecuencia  del 
fen6meno. 

I  I  I  I  I  I  I 3  I 4  I 5  I 6  I7

0.001 0.01 0.1 1 10 10 2  10 10 10 10 10
FRECUENCIA (Hz) 

Iof­ Fen6menos Electromecimicos  ｾ Iof­ Fen6menos Electromagneticos  ｾ I 

Control ca rga  Estabilidad  Corto  Maniobra  Ondas viajeras  
fTecuencia  transitoria  circuito  de IIneas  

Estabilizadores   Resonancia  Voltajes TRV 

Subsincrona 

Arm6nicos 

Figura 1.  Clasificacion de los transitorios de acuerdo con  la frecuencia 

Ｑｾ＠
r--

1.3 OBJETIVOS

Los  objetivos  basicos  propuestos  en  el  curso  de  fen6menos  transitorios  son  los 
siguientes: 

•   Comprender  la  naturaleza  de  los  transitorios  electromagneticos  en  sistemas 
de  potencia . 
•   Conocer  diversas  tecnicas  de  anal isis  y simulaci6n ,  como  tambien  sus 
aplicaciones. 
•   Obtener los  conocimientos  basicos que  resultaran  indispensables  para  cursos 
de  coordinaci6n  de aislamiento. 

En  forma  practica ,  los  ingenieros  analistas  no  s610  deben  conocer  el  fen6meno 
transitorio, sino dar soluciones que apunten  a mejorar la  calidad  de  la  potencia  y la 
confiabilidad  del  SEP.  Para  cumplir con  10  anterior se  deben  seguir los  siguientes 
cuatro  pasos: 

•   Conocer el  fen6meno:  Origen,  formas  de onda , magnitud,  duraci6n,  frecuencia 

de  ocurrencia ,  evoluci6n del fen6meno. 

•   Conocer la  respuesta de  los elementos a este tipo de excitaci6n . 

•   Coordinar  adecuadamente  la  ocurrencia  del  fen6meno  transitorio  con  las 
caracterfsticas de los elementos del  sistema. 

•   Proteger adecuadamente, mediante el  diseno,  pruebas  y aplicaci6n. 

 2. INTRODUCCION AL PROGRAMA EMTP/ATP

EI  programa  ATP  con  respecto  al  curso  de  Fenomenos  Transitorios 
Electromagneticos  tiene  una  aplicacion  clara  y  es  la  de  simular  todas  las 
situaciones  que generen  esfuerzos en  forma  de  sobrevoltajes  y/o  sobrecorrientes. 
Estos esfuerzos exigen  los  aislamientos de  los equipos,  generan  efectos termicos 
y pueden  producir tambien  esfuerzos mecanicos. 

EI  programa  ATP  (Electromagnetic  Transients  Program),  es  un  programa  para 
computador  digital,  utilizado  para  la  simulacion  de  transitorios  electromagneticos, 
electromecanicos y funciones de control en  sistemas de potencia  polifasicos. 

2.1 HISTORIA

EI  ATP  fue  desarrollado en  su  version  inicial  por el  Dr.  Hermann  Dommel  a finales 
de  la  decada  del  60  en  Alemania  (inicialmente  con  el  nombre  de  EMTP),  quien 
posteriormente  venderia  los  derechos  a la  Bonneville  Power Administration  (BPA) 
de  los  Estados  Unidos.  A  pesar de  ser  la  BPA  la  que  coordino  la  distribucion  del 
programa  entre  los  diferentes  grupos  de  usuarios,  muchas  universidades 
contribuyeron  al  desarrollo  de  los  diferentes  modelos  que  constituyen  10 que  hoy 
es  este programa. 

Del  programa  existen  versiones  para  diferentes  tipos  de  computadores.  La 
version  para  micros  es  la  mas  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  interaccion  con  el 
usuario  es  mediante  un  rigido  archivo  que  debe  cumplir  ciertas  norm as  de  dicho 
lenguaje.  Actualmente  existen  program as  tipo  interfaz,  que  permiten  un  trabajo 
mas  amigable  para  muchas  de  las aplicaciones.  La  interfaz que  se  utilizara  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  version  inicial  para  el  sistema  ｯｰ･ｲ｡ｴｩｶｾ＠ DOS  hasta  las  versiones  actuales 

?1  
para  el  ambiente  Windows.  Esta  interfaz  cubre  una  buena  parte  de  las 
posibilidades  que  tiene  el  programa  EMTP/ATP  pero  siempre  se  debera  tener  un 
conocimiento  basico  de  c6mo  opera  el  ATP,  que  es  realmente  el  programa 
simulador.  La  interfaz  ATPDRAW  normalmente  se  puede  actualizar  vIa  Internet; 
algunas direcciones de donde se  puede bajar son : 

hNp:/laJpha.kisarazu. ac.jp/-aFTPlatpJatpdrawl
hNp://www.ee.mtu.edulatp

2.2 ESTUDIOS A REALIZAR CON EL PROGRAMA EMTP/ATP

A  pesar de  que  el  objetivo  principal  del  programa  es  la  obtenci6n  de  la  respuesta 
transitoria  de  los  sistemas  ehktricos  de  potencia,  tambie!n  se  puede  obtener  la 
respuesta  en  estado  estacionario  para  un  sistema  alimentado  con  CA  (respuesta 
fasorial) . 

Para  la  simulaci6n del  sistema de potencia,  el  programa  posee varios  modelos que 

dan  una  representaci6n  adecuada  de  los  diferentes  elementos  que  10  componen . 

Los  modelos disponibles se  pueden  clasificar as!: 

•   Modelos  constituidos  por  elementos  concentrados  RLC .  Estos  modelos 

pueden  ser  simples  ramas  RLC  serie,  circuitos  PI  polifasicos  que  pueden 
representar Ifneas de  transmisi6n  0 transformadores. 

•   Modelo  de  onda  viajera  para  representar  adecuadamente  una  linea  0 un 
cable . 
•   Impedancias  de  tipo  no  lineal:  Resistencias  no  lineales,  inductancias  no 
lineales, resistencias variables con  el  tiempo. 

•   Suiches  ideales  control ados  por tiempo.  Suiches control ados  por voltaje  para 

la  simulaci6n  de "gaps"  en  flameos de aisladores. Diodos  y  tiristores.  Suiches 
de  acci6n estadlstica. 
•   Fuentes  de  voltaje  y  corriente  ideales  de  tipo  escal6n,  sinusoidal ,  rampa , 
exponencial  y  definidas punto a punto . 

•   Modelo completo para la  maquina sincr6nica . 

??
 •   Modelo  universal  de  la  maquina  electrica  que  permite  representar  doce  tipos 
de  maquinas  diferentes.  EI  modelo  de  mas  interes  en  este  grupo  de  modelos 
es  el  del  motor de induccion trifasico. 

•   Representacion  de  los  sistemas  de  control  mediante  la  opcion  de  los  TACS 
(Transients Analysis  Control System).  En  las versiones mas recientes del  ATP 
existe  una  opcion  alterna  para  representar  la  parte  de  control  del  sistema  de 
potencia: MODELS.  En  la  figura  2 se  observan  los  diferentes modulos de  que 
dispone  el  programa  ATP  para  generar  modelos.  A  estos  modulos  se  les  ha 
denominado  subprogramas  de  soporte  y algunos  de  ellos  se  pueden  manejar 
directamente desde la  interfaz Atpdraw. 

Simulation part Supporting programs

ｾ ｉ＠ U NE CONSTANTS 

... ­:­ICABL.E ｃｏｈｓｔｾ ｔｓ＠

r--------------....;...;.." ... : I  S8Ift. 'fBI! SETUP
time­domain WId 
f naoquenc.y ­ domsln  ...  f  JMAR1'.  !'!ETUP 
sofubDns 
.. ICASLE ｐａｒｈｯＱｅｔｾｓ＠
repr••entatfon 0' th e :rt 
• •
NODA SEfUP 
...ctrlcal network
.. ARMN'rr p
spy ,
..  . 
:
I BC TR A" ,

ｉ＠ ｾＧｊＭ＠
ｾ i I 
I 
X  FORMER 
:  SA T U  RA 

8 H  V SO A T 

Tr.., SI8I'l1 G eneral purpose  Z.N(l"F ITTER

An atvsis of 
Qr Systems 
C t.4'llt 

TACS - Idmufation 
Ｑ Ａｴｬｾｑ･＠

MODELS
r DAT;aASiNooUl.e  ., 

Figura 2.  Programas de soporte que  interactuan con  el  EMTP/ATP 

EI  ATP permite hacer entre otros los  siguientes estudios: 

•   Calculo de  parametros en  lineas aereas y cables subterraneos 

•   Calculo en  estado estacionario sobre  un  sistema de potencia 

23

;' 
•   Sobrevoltajes por maniobra (suicheo) 

•   Sobrevoltajes por descargas atmosfericas 
•   Calculos de cortocircuito que  involucran los  sistemas de  retorno 

•   Coordinaci6n  de aislamiento 

•   Calculo  de  voltajes  inducidos  sobre  elementos  cercanos  a  Ifneas  de 
transmisi6n  (otras Ifneas, cercas , oleoductos) 

•   Resonancia en  Ifneas paralelas 

•   Simulaci6n  de  arranque de  motores 

•   Evaluaci6n  de arm6nicos 

•   Estudios de ferrorresonancia 

•   Maquinas electricas 
•   Simulaci6n  del  control de  las maquinas y del  sistema de  potencia 

2.3 ESTRUCTURA GENERAL DE UN ARCHIVO PARA LA ENTRADA DE

DATOS AL ATP

EI  programa  ATP  es  un  paquete  de  simulaci6n  digital  de  transitorios 
electromagneticos escrito en  lenguaje Fortran.  EI  programa cuando ejecuta lee los 
datos  de  un  archivo  que  debe  ser editado  previamente.  Este  archivo  esta  ceriido 
por  10 tanto  a  las  reglas  de  ese  lenguaje  de  programaci6n.  Esto  hace  necesario 
conocer en  detalle como esta organizado el  archivo de entrada . 

Existen dos formas  generales de utilizar el  programa : 

•   Para  hacer una  simulaci6n  de  un  sistema  de potencia  (0  una  red  circuital)  con 
el  fin  de  obtener  la  respuesta  transitoria .  EI  sistema  bajo  simulaci6n  puede 
contener  tanto  elementos  circuitales,  como  elementos  de  control  que 
interactuan  con  la  red . Se  pueden  simular de  manera  independiente  sistemas 
de control. 
•   Obtenci6n  de  para metros  correspondientes  a  elementos  del  sistema  de 
potencia,  como  Ifneas  aereas ,  cables  subterraneos ,  transformadores, 
pararrayos , etc. 

24
Las  posibilidades  de  uso  comun  del  ATP  como  programa  de  simulacion  0 de 
calculo  de parametros son  las siguientes: 

•   Caso simple de simulacion  de una red  electrica sin  bloques de control. 

•   Caso simple de simulacion de  un  sistema de control  sin  incluir la  red  electrica. 

•   Caso  simple  de  simulacion  de  red  electrica  y  sistema  de  control  en  forma 
conjunta . 

•   Obtencion de para metros y modelos de Ifneas de transmision . 

2.4 REGLAS GENERALES PARA LA GENERACION DEL ARCHIVO DE

ENTRADA DE DATOS

EI  archivo  de  entrada  de  datos  debe  ser  valido  para  lenguaje  Fortran .  Las  reglas 
generales son  las siguientes: 

EI  archivo debe tener un  nombre  y se  recomienda  que tenga  una  extension  *.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 
extension *.RES  0 *.LlS 

EI  archivo de entrada de datos tiene 80 columnas. 

EI  archivo de entrada  puede tener Ifneas de  comentarios,  las cuales empiezan  con 

el  caracter "C" en  la  primera  columna,  seguido  de  un  espacio en  blanco  y luego el 
comentario que se desee en  la  misma  linea. 

Los  valores  reales  con  formato  cientffico  (5.0E­2  por  ejemplo)  y  los  numeros 
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  mayorfa  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  0
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. 

25 
No  se  deben  dejar  Ifneas  en  blanco  arbitrariamente,  porque  estas  las  utiliza  el 
programa  para  identificar  cuando  ha  terminado  un  grupo  de  datos, como  se  vera 
posteriormente. 

2.5 ESTRUCTURA GENERAL DEL ARCHIVO DE DATOS PARA UN CASO

SIMPLE DE SIMULACION TRANSITORIA 0 DE ESTADO ESTACIONARIO DE
UNA RED ELECTRICA 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 
utilizaran  unicamente  ramas  simples  del  tipo  RLC  con  parametros  concentrados  y 
fuentes  simples del tipo escalon  y senoidal. 

Las  partes  del  archivo de  entrada  de  datos  para  este  tipo  de  simulacion , en  forma 
general  se  observan en  la figura  3. 

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 , identificacion de 

sf  es  una  simulacion ,  0  es  un  caso  de  calculo  de  para metros  para  algun 
elemento de  la  red  (acceso a subprogramas) . 
•   Si  es  un  caso  de  simulacion ,  el  siguiente  grupo  de  datos  corresponde  a  la 
especificacion  del  tiempo  de  simulacion ,  tiempo  de  paso  de  integracion, 
especificaci6n  de  las  unidades  de  los  para metros  LC  y  control  de  salida  de 
informaci6n.  Tambien  se  especifica  en  este  bloque  de  datos  si  es  un  caso  de 
simulacion de estado estacionario. 
•   Datos  de  las  ramas  de  la  red  (modelacion  de  los  componentes  de  la  red  de 
potencia) . 

•   Datos de  los interruptores. 

•   Datos de  las fuentes . 

•   Especificacion  de  las condiciones  iniciales,  si  existen . 
•   Lista  de voltajes nodales a ser tabulados y graficados. 
•   Bloque  de  datos  para  unas  opciones  graficas  y  de  calculos  especiales, como 
anal isis de Fourier y tablas con  resumenes estadisticos. 

26
Para  una  simulaci6n  en  particular  no  deben  aparecer  necesariamente  todos  los 
bloques  sefialados.  En  la  figura  3  se  observa  la  presentaci6n  en  detalie  de  los 
formatos  de  las  diferentes  Ifneas  de  un  archivo  para  simulaci6n  transitoria  0  de 
estado estacionario de un  caso  simple que no contiene bloques de control. 
Un  caso  simple  de  simulaci6n  se  entiende  como  aquel  donde  no  se  hacen 
estudios  de  tipo  estadfstico,  no  hay  bloques  de  control ,  no  se  utilizan 
subprogramas  para  calculos  de  parametros  de  elementos  del  sistema,  s610  se 
emplea  un  tipo  de  rama  (rama  RLC  serie  de  parametros  concentrados)  para  la 
modelaci6n  de  la  red , el  interruptor que  se  usa  es  controlado  por tiempo  en  forma 
determinfstica y se  utilizan fuentes sencilias como escal6n y cosenoidal. 

En  la  figura  3 se  observa la  forma  como estan  organizados los  diferentes datos en 

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

27
101  BUS  IU/II  Amp  Pha  Al  Tl  Tsta  Tstop 
BLANK  BRANCH 
BLANK  Svl  l TCH 
BLANK  SOURCE 
ESPECIFICACION DE  VOLTAJES  NODALES A CALCULAR  Y/O GRAFICAR 

NODOI  INO D 02  INOD03  INO D 04  INOD05 

BLANK  OUTPUT 
BLANK  PL OT 
BEGIN  NE w  DATA  CASE 
BLANK 

Figura 3. Formato general para el programa EMTP/ATP

28
 3.  PRAcTICA N01.  CASO SIMPLE DE  SIMULACION DE  UN  CIRCUITO RLC 
SERlE 

Para  realizar  una  introducci6n  en  10 que  es  la  utilizaci6n  practica  del  EMTP/ATP 
como  herramienta  de  simulaci6n ,  10 mas  apropiado  es  hacerlo  con  un  sistema  que 
sea  sencillo  (caso  del  circuito  RLC  serie),  para  el  usuario  que  se  inicia  en  su 
utilizaci6n,  y  que  a  la  vez  permita  formarse  una  idea  de  la  forma  de  utilizaci6n  del 
programa.  A partir del circuito RLC  se  pueden  dar los  primeros conceptos asociados 
con  los  fen6menos  transitorios  reales  que  se  presentan  en  un  sistema  electrico  de 
potencia . 

3.1  OBJETIVO 

EI  objetivo  de  esta  practica  es  familiarizarse  con  la  forma  del  archivo  de  entrada  de 
datos  para  un  caso  simple  de  simulaci6n  transitoria  de  un  circuito  RLC  serie. 
Adicionalmente, se  plantearan algunas situaciones en  las que se generan transitorios 
simples donde el  modelo circuital son  circuitos  RLC  simples. 

3.2  DESCRIPCION 

EI  circuito  RLC  serie  es  una  primera  aproximaci6n  a la  modelaci6n  de  componentes 
de  un  sistema  de  potencia , como una  red  de transmisi6n  de  energfa .  En  forma  real 
estos efectos,  resistivo­inductivo­capacitivo,  se  presentan  de  manera distribuida y no 
concentrada, como  es  el  caso  del  que  se  ocupa  esta  aplicaci6n.  Sin  embargo, hay 
una  serie  de  definiciones  que  se  pueden  hacer  sobre  un  circuito  RLC  serie,  que 
coinciden  cuando  una  red  de  transmisi6n  se  modela  mediante  parametros 
distribuidos.  EI  circuito a simular se observa en  la  figura 4. 

Para  observar los  estados trans itorios  de  voltaje  y de  corriente  sobre  el  capacitor se 
pueden  utilizar cualquiera de  los circuitos #1  y #2. 

EI  primer circuito incluye un  suiche conectado entre los nodos FTE y NSW,  que se 

cierra en  el  momenta que se desee.  La fuente se activa en  cualquier momenta 
mediante el  parametro TST ART. 

29
 ･
Ｌｾ＠ Ｌ ［ｾ＠
,  ; ( ｾ

ｾ ｶｾ＠
ｕＧｾｬ ｶ ｲｒ ＼［ ｉ｛ｬａ ｻＩ＠ N AC IONAl D E  CO LOMDIA 
5t.it1  "4I'}JIo.u JN 

DEPTO. DE BTBLIOTECAS
B I BLiOTECA MINAS

FTE  _i{t) 
_>  VC 
100  Ｎ＠ ｘ ｾ＠
:I:' . V NSW  T VC{t) 
Cto #1 =

i(t) 
FTE  ­­>  VC 

­:r . V T VC(t) 
Cto #2  =
Figura 4.  Circuito RLC  serie 

EI  segundo circuito  no tiene  suiche, pero se consigue el  mismo efecto al  activarse  la 

fuente en  un tiempo igual  al tiempo en  que se cerr6 el  suiche del  primer circuito. 

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

v Rt
ｉ ｬＭ ｾ＠
i(flJ = V  ｾ ａ＠ t
!  1 x e '  2L ' x  sen ＮｾＧ ｃ＠ｌ (3.1 )

Zc.p - 4,e

donde Zc es  la  impedancia caracteristica de sobretensi6n 

L
Zc = ", C (3 .2) 

Y A es el factor de amortiguamiento, 
Z
A =­=-
R
(3 .3) 

30  
r  fT· 

Cuando  el  efecto  resistivo  es  despreciable,  la  corriente  tendra  un  comportamiento 
senoidal  de  amplitud  constante  (no  hay  amortiguamiento).  La  corriente  tendra  un 
comportamiento en el tiempo de acuerdo con  la  siguiente expresi6n: 

i(t) = ­v x sen
Zc 
(1 
(LC t)
. 
(3.4) 

La  anterior  expresi6n  permite  dar una  interpretaci6n  de la  impedancia  caracterlstica 

Zc, como la  relaci6n  entre  la  amplitud  de  la  fuente  del  voltaje  de excitaci6n  y  el  valor 
pico de la  corriente senoidal del circuito serie. 

La  frecuencia de la oscilaci6n senoidal esta dada por: 

fa = 2TT I
, LC
Hz (3.5)

3.3 PROCEDIMIENTO

..l .  . ｸ
­:r  V 

Ｎ［ｷｾ＠
ｾ＠
I ｔ ＧｃｾＩ  

Figura 5. Circuito RLC  serie en Atpdraw 

Mediante  la  interfaz Atpdraw elaborar el  circuito en  forma  gratica  que se  observa  en 

la figura  5. 

Para  editar  cada  uno  de  los  componentes  se  abre  la  caja  de  herramientas  con  el 
bot6n  derecho del" mouse": 

31  
­; 

Erot>et" 3­pha... 
B.ranchu­r 
Stanch tl"""*r 
ｬＮｩｮｳ＠ ｾ＠ • 
line Qistribuled  • 
S!jiIches 
ｾＰＱＶ ｟＠

lIIiCtn. 

T[......"'fIl... 
MODELS 
Typej4 
lACS 
l.I_Specfied 
LineIC!bJe 
llverhead l ine  tFCH) 
frequency  001!!P. 

Figura 6.  Caja de herramientas de la  interfaz Atpdraw 

Editar cada  una  de  las ramas  para  formar  la  red  pas iva  RLC , los  datos  del  suiche  10 
mismo que la fuente tipo escal6n de acuerdo con  la figura 7. 

ｾ＠ ｾ Ｇ＠ ｾ＠ ｬ＠

T 

"'- f P­­­­I 
ｾ Ｇ＠ ' ... 

_._J ｾＭＢＡＮｊ＠

Figura 7.  Ventanas de datos para  rama  RLC,  suiche  y fuente escalon 

ASignar los "settings" correspondientes a la  simulaci6n de acuerdo con  la  ventana de 

datos que se observa en  la figura  8. 

32 
 )( 

S'rnul§ !on  \!ntI!go,  ｾ I I lId'«  I Becom I UM 

Jjeko T:  12E­5  Simulalforlwe ­ ｾ＠

Im"",  "'°":
1 ":02: - - - - - r.  lime dommn
r  !iaq\lenCY i Ce.n
ｾ ｰｴ＠ 10 
.(;opt ｾ＠ t  r hlllrnol\lc (MFS) 

r  _  Frequency 

ｲ Ｍ ｾｋｬ ｉ＠ t!elp  I 

Figura 8.  Ventana de datos para  los "settings" 

De  igual  manera deben asignarse los "settings" enteros: 

)( 

s..nut{IOOn  "",,9.. Ｑｾ !ｾ＠ 1110<o", IUM  1

Ｇ ｏｾ ｾ ｃｏｲｴｬＧｯ ｴ＠ Pri.·')."tII 

i>tinlll1'Q, '  P S ­:  ［Ｗ＠ Ｌｾ Ｂ｣ｮ ｲＺＡｴｊｋｩｾ＠

ｆｬｉ ｯｴｆｲ
｟ｾ＠
Ｌｾ .• r ­P' .. 810.",..,.. ",.""."
P  """",,,",_s 

- 5'Jf.!ch sQjd')r r 240. prum14 c:cntrd

I s._>t.4Y r 
Ｌ＠ ｾＢ ＢＧ ｃＢ｣ｴｹ＠ r 

ｍＮ ｾ＠ r 

ｾ＠  

Figura 9.  Caja de herramientas para datos enteros 

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

Es  importante  que  los  diferentes  parametros  que  controlan  las  simulaciones  se 
conozcan  bien.  Una  primera  fuente  de  informaci6n  es  el  HELP  que  aparece  en  la 
mayorfa de las ventanas. 

33
 _!TF
ｾ Ｎｾ｜ｐ＠
f iliUS
. ｾ＠
ｅ ｉｾ＠Ｖ
'' 'aTP  COi'Al.ij
lun PC:Ftll ｾ Ｎ Ｑ＠

nmATP IIW CI!ooIl.2
MPtaO<Y Q\-Aa . 3
i,diT.. ｾ Ｎ＠

Figura 10. Persiana "ATP" para realizar el "Make File"

Cuando se tengan completamente editados todos los componentes se realiza un

"Make File", accion que corresponde ala realizacion del archivo *.atp de simulacion .

EI archivo de simulacion *.ATP se observa en la tabla 1

Tabla 1. Archivo plano generado por la interfaz Atpdraw

BEGIN  NEW  DATA  CASE  
C  ­­ -
C  Generated  by  ATPDRAW  agost o,  jueves  17,  2000  
C  A  Bo nnevi l le  Power  Administration  program  
C  Programmed  by  H.  K.  H01dalen  at  SEfAS  ­ NORWAY  19 94 ­ 98  
C  ­ ­ ­­-
C  ｾ ｳ｣･ｬ｡ｮｯｵ＠ Data  Card  
C  dT  ><  Tmax  ><  Xopt  ><  Copt  >  
2 .E ­5  . 02 
25  1  1  1  1  o  o  1  o  
C  1  2  3  4  5 6 7   8  
C  345 67890 12 3 45678 90 123 4567 890123 45 6789 01 234 5678 9012 3 45 6789 01 2 34 567890123 4567 890  
/ BRANCH  
C  <  n  ｾ ＾ ＼＠ n  2>< ref l><re f2 ><  R  ><  L  ><  C  >  
C  <  n  : ><  n  2>< refl> <r ef2 ><  R  ><  A  ><  B  >< Leng ><><>O  
NODOI  NOD02  1 50.  1 63 .   2 
NOD02  . 7011 2  o  
/S WITC H  
C  <  n  1 ><  r.  2><  Tclose  >< Top /T de  ><  Ie  ><Vf/CLOP  ><  t yp e  >  
FTE  ;';000 1  1 .  1  
/SOURCE  
C  <  n  1 ><><  Ampl.  ><  Freq .  ><Phase /T O><  Al  ><  Tl  ><  TS TART  ><  TSTOP  >  
llFT E  0  1 00 .  ­ 1 .  1 .  
BLANK  BRANCH  
BLANK  SWI TCH  
BLANK  SOU RCS  
FTE  NOD02 

Las caracter[sticas de este archivo de simulacion son las siguientes:

• Las I[neas que empiezan con la letra "C" son simples comentarios para
documentar el archivo. EI programa en la ejecucion las ignora.

34
•   EI  tiempo  maximo  de  simulacion  es  de  20  ms.  EI  delta  de  tiempo  es  de  20  J..I.s. 
Lo  anterior da  un total de 1000 puntos de calculo. 

•   Se imprimira una tabla con  los datos de las respuestas transitorias de los voltajes 

y corrientes en  el  circuito,  con  un  intervalo  de  impresion  de  25,  para  un  total  de 
datos de  41  puntos en  el  archivo de resultados.  Este archivo de resultados tiene 
extension  *.LlS 

•   Se  solicita  una  tabla  de conexionado de  la  red,  10  mismo que los  picos  maximos 
y minimos de las variables en  el  tiempo. 

•   En  la columna 64 de la linea correspondiente a datos miscelaneos enteros (linea 

3a  del  anterior  archivo,  excluyendo  las  IIneas  de  comentarios)  se  ha  solicitado 
mediante  un  "1"  la  generacion  de  un  archivo  con  extension  *. PL4, el  cual  sirve 
para  que otras  paquetes  de graficacion  como el  PLOTXY y el  TPPLOT generen 
una grafica , con unas caracteristicas de resolucion  muy buenas. 

•   Para  la  modelacion  de  los  elementos  del  circuito  se  utilizan  dos  ramas tipo "00" 
(rama  tipo  RLC  serie).  La  primera  es  una  rama  con  valores  R,  L,  C  de  150  Q , 
163  mH  y 0.0 J..I.F.  La  segunda  rama  con  valores  de 0.0 Q,  0.0 mH  y 0.7042  J..I.F . 
Las  ramas RLC  no  puede tener los tres valores en  cero.  Cuando se asigna  cera 
a  la  capacitancia,  el  pragrama  entiende  que  es  una  rama  RL  y  simplemente 
asume  que  no  existe  el  para metra  capacitancia,  porque  de  10  contrario  un  valor 
cera de capacitancia  puede  interpretarse como  una  impedancia de  un  valor muy 
grande. 

•   Como fuente  de  excitacion  se  utiliza  un  escalon  de  100.0 voltios  de  amplitud , la 
cual  tiene un tiempo de activacion de ­1.0 segundos. 

•   No  hay condiciones iniciales (las CI  son  cera) . 

•   Como variables de salida se solicitan  los voltajes nodales correspondientes a los 

nodos  FTE  y  NOD02.  Adicionalmente,  se  habia  solicitado  el  voltaje  y  la 
corriente para la rama conectada entre nodos NOD01  y NOD02 (rama RL) 
EI  archivo  correspondiente  a  la  salida  de  resultados  de  esta  simulacion  se  observa 
en  la tabla 2. 

35  
Tabla 2.  Archivo de resultados 
ｃｯｴＡ･ＺｰＧｾｶＮ＠ Lflt.etpretat i·:.r:  ｾ Ｈ＠ wp u t  ､ｅｯＮ ｴ ｾ＠ e.rd, . ﾷ Ｇ［ｦｯｊｾ＠ J.'II  :ud  i.&!.aqe!l  o.re  sr..:ur, bel?<ol •  .1 :  £Iv  ｣､ｾＮ＠ o;ha..rA<;t4it  'Y  ＺＱｉＧｲＤ＼ｃ＠ Ｚ ｩｉｾ＠

ｑｾｚ＠ H  Ｕ＠ ｢ ＿ ｴｮ ｾ 12H ｾＶＱｒ＠ j'j 1_ 3 E.6"'.s  ＬｲＮｬ＠ ｾ H  ｾＧＶｔｦ［ＮＹＰ＠ Ｚ＠ .:'J. ｾＱＰ｡｝ｯＮＡ H!>€i a ｾＨｬＡｊ＠ Ｇ＠ ｾＶＷＱｊ｝｜ｮｈｩｴ ｍＧｾ＠
... ­ ­­­­­­­ _.  ­ ­­ _. 
,:­.(IIDItflt  ::: .. t d.  l'fI.'ta:ro  ..  ,  .  ­ ｾＧＮＢ＠ Ｚ＠ :: \ " Tf".  ,.­ Art"  , P.;.c , {{Ii 
ti­uic..:  ::.. :d  ｐｴＮｲＺｾｈ｜ＧＱ＠ r,"",  Fl1Tr· ｾｴｬｳｵ･＠ Ｎ＠ Jii;Glr,'  ｦｾ Ｇ＠ cn'i'A  -XF.
ｾ ｾｮ ﾷ＠ .ar::: .  ｴｲｾ＠ ­ J .  II: 
ｾｲｴＧ＠ ,:. ",.;j.  ｾＧｍＩｃｄﾷ＠ t .  I''::  Get1oerat.d  Coy  ｜ｮＧ ｴＮｾ＠ "g"X!lt"l) .  P  .  1,;)  'wtV_ 
":".cQlIIer....  Sf,;.! .  ft'L'MOC':)  1] .  C  A  S.,,..r .... H  .:. e  L:wec  JrdaUni.:l":.toS'. {­:­tI  ［Ｍｾ＼ｏＺＮ｜ｴｉ＠
ｾｲＮ＠ :.rd .  ｾ Ｑ Ｎｴｙｌ .......::.­ f . I:  ｲ＠ ｲ］ｮ ＡＢｉｊｃＮｏｯ ｾ＠ t ',  11  .  K.  ＢＧＬｾＺＮｲ＠ Ｌ＠ ｾＮ＠ n:: .ＮＱ  ｾ＠ ­ lii­'P.',,;r.',·  H­:" 
ｾＬ Ｇ＠ ",3 rd .  N:lt1X'O,  1 . 
ＺｊｉｲＢｾ＠
c.u'i .  NUtO:!;  E.  .  ;  ｾＬｴ .• :.U.t'\..,...s  ｯＮ ＧＺＮｾ＠ ­:.a.!":::! 
'::'=­­n'.  card.  f('J. .［Ｚｲｾ＠ ..  ,  .  j 7  x  l'IIUIX  )(  XOpt  )oo  ｾＺ＾＠

M,i!': .  ＼ｊｬ｜Ｂｾ＠ Ｎ＠ :  Ｎ＠ ｦｊｏ ｏｅＭｻｉｾ＠ c­oo ;:­r!l  . c,Qg  ｾ＠ ･ＺＮ . i._c. 

M,i:!':.  dGt •  •  1  1  ­' 
ＢＧｾ Ｎ＠ ｮＢＺＮ＠ ｾ＠ •ＮＬｾ＠ ...  ｾ ｕｍｏｃｴｉ＠ ...  11 .  l C  ｾ＠

r.­:  .• r.!.  !n..')1JCtl.  :3 .  :  ｾ＠ Ｂ Ｂｬｦ｡ｏ ::106:H,2 1 4'T'890 lZH 'lii"'1!!9Ji  ｾ＠ Ｂｾ f " ) ..,;) ｾ＠ Sc"' P;:1"",': :, ｾ＠ 1i"'l'=lJ:":' '"  5,; ­ 91) 
ｾｮＺＮ＠ ｾ＠ ｡ｴ､Ｎ＠ ｾＧｴＮｈＩＺｏ＠ :1.  ­::  <:  !)  1'><':  r.  ［ ｨ Ｇｮ Ｎｴ［＾ﾷＧＡＢｦｾ＠ ­ x.:..  y  ­
Ｚｾ "  tllcd.  Ｌ＠ Ｇ＠ ｩ＠ ｕＡＱ｛ｾＧＢ＠ L,  ｾ＠ Ｂ＠ ｾ ｬ＠ •  •  1\  ｊＮ＠ ｾ＠ l  lx'.("1·<  X  A  r< a.  Y!.ioIT,,:;'­<  Y·. ,)

:.lQ r .1V::  ｆＭｾＨＧＢＮ＠ 1..  ＮｬＧｮ ［ﾷ｣ｾ＠ ｬ＠ Ｎ＠ ＢＧｊＮｾｴ ［＠ ｾ Ｚ＠ P:/JIJ:JE'  Ot1  h<".'fr"l  Ａ＠ Ｎ＠ Ｚ＠ ｊＺｮ ｾ＠ l(. !  . 
ｾｔｩｴＡｬ＠ l:.­t­: .  . .  ｊＮＧｾｅ＠ ｾｏ＠ . . ";OO E o[J.'  Ｇ＠ Ｎ＠ ＢＺｩｾｲＭ＠ 1UlO,  . ­0 4 
8"'o.."1l  :.1Ird  " n·:hn.J  ｨｲｾｮＬＺＱｩＡ＠ Ｎ＠ I&iI,  H';'j;"  . ; :   j  IfiLAhi­'  BPJIJi ,,.. 
ｉｲＺＢｮＧｸ＠ Ｇ＠ ｾＢｩ［ｬ ｾ Ｇ＠ ＬＢ＠ Ｚ ｗｔ＠ Ｌｪ ｾ＠ 1.,.  )o<  .. ,! ,r­u' r'  X  .. )1....  :> 
(1' .  rJi'JE: ·  l'ItJ  ')  . 1)( £ ,  f  rTF  1¥)l'I­·l  1. 
ｾﾫｔ［ＱｯＨ＠ •  I  .&.l.MW ｾＺｴ ＧＺＭｈ＠
C'XGen.'__ uJ .  ｴｬｾｄﾷ＠ _2.  =
<  r.  1">­(")·  A."":j:i .  ....:  f r eq.  y:"hll,!,,,,' T):--
"  ><  Tl  x  7:)':'AP7 x  : :!"':"T ｾ＠

ｓｾｵ｣･ Ｎ＠ :.C'")'; +02 '. CO E.. OO ｏＮＨＺｾｅＧ［＠ ｾ＠ flOSH.Q 'ilfif.  L 

Bl.r.ic C8t:l Ｎｾ､ ... ｬ｣ｏＺ＠ ＧＺｲｾＭ ｮＮｾＺＭＧ＠ "'OlLICe.x . iH...N'r. liO',;­F.c­r 

!....i.1 '. d i r.put e-i.::-I'I:;t ·, ｴｊＬＮｾ＠ .. cw Ｎ［ｾｲｪＭＺＡ､＠ ｾＡＺＮ＠ • • .::h nooe . ｃｴ｜ｾ＠ ｹ＠ ｾ Ｎ＠ ｨＮ＠ f.h:.''!'l. ｾＬｕｵＮＺ［Ｇ｟ｴｲｊ＠ ,:. f  ｾＬ ... ltl-f'I""'.1IIiI ｬ＠ ｾｲＺ .... ｡ｻｾ＠ 21"; ...' !. ｡｛ＮＴＺｬｴｾＧ Ｏｾ＠
..n1l.lrl'J:·. h"" ｾｰｬｃｩＱ＠ ace ｩ［ｲＮﾷＺｴｾ＠ f llp00!4--j .r.·t · ,!\ ;no-Hr.!Jt.o.. ｰＦＨＮＡｉｬｾｩ＠ Ｇ＠ ＺＭｍ｀［＠ ＢｩＮ Ｎｊｦｬｾ＠ Ｑ＠ ［ｷｩﾷｯＺｬｾＢ＠ "'. ir..:bOaa. a.ld.Clli1b  "'·,'.Ir::ot -
Ｑ＠ ｾｊＧＱＮＺｵ､ｩｮｴ［＠ f"'.• -:. iu'3 sa.::h.iLuy IlIa ｾｩｮＺ･ＭＱ＠ - II­:­;'.p­ UulL U. !1'. !tsaq. p::-cdu=a . ＹＮｘｾ Ｎ＠ ｲ＠ I.  ｬＧＬ［ＮｭｾｩＭ､ｦ＠ Ｚ＠ ｮＮｵ＠ r.u..i • ..I:  "L'K'O"./.X ·' .
---------,----------
FtJ::. ｾ＠ ［＠ ｾ＠ I ｾＺＬ＠ ｣ ｾ＠ .1 t !I­!; h",r.t ＺＮｕＭＧ ｾｃＺＮ＠

NOOOt  00'00.:  ｾ＠ ｾｅ＠

ｾｏＺＭＮＢｌ＠ rTS"F:AA · t/ '(.o)l

,","E  100001 '
Ｚｾ ｩｬｊＮＱ ａ＠ I  NUIX·.;. •

COrd : f  n&."'OO. (0:. f  tm.-5Cep 1'OOP outP\J!.. • r.E /""rn:

｡ｾＺＵ ｮ Ｎ Ｇ Ｎ＠ --41, 11 endin1  :-('q-... ｾＡ＠ 1 !cr  ;:lJ '.p...t. v,u·iatl "." . ｬ＠ ｾｩＮＢｈｲ＠ CUTf'UT 

Colu.c:n ｮＮｴ､ｬＹＧｾ＠ for t.h41 3 OO'f  .;A.Itput  v.tl "bl 'ltll ｦｯｬｾＮ＠ rh.. se .Ct'l di... ld"d "",on.:, ­n.  ｾ＠ po-uibh Ｎ･Ｚｾｐ＠ H · [ , ·,l"7oolS
F!r.!t. ｑ｜ｬｾｐｕ
ｉ＠ 'lOt Lot;.le:!- !!Ire ･ｬＮ｣ｴｌﾷｮｾｷＬ ＺＮ［＠ vr.>lt"lIft ｈｲｴＮＶｬ｣･ｾ＠ ( ll1=p8["  ＭＮＧＺｾｴｲ［･＠ ｉｴｬＮｲ［Ａｾ＠ ｌｾｲ＠ ·ntr. ... 1...
'i€lx t l eancn ｣ＧｊＮｲ
::<IJtp" " ·11" 1A:I\ .. '  ""C6Ｚ･ｮ ｾ＠ ｴ＠ Lｾ＠ ｬｯｐｊｾ ［ ｾ＠ Ｑ＠ ! .:'" '%l ｦｾ＠ ｴＺｾｲ＠ r.a­u  ­ .. ｾ .•  1.,:"-,,,,[ ｲｾＬ＠ ［＠
S;­ ­r  tUM  ｴｦＨＺｄＧｾＱ＠ n[  ｎ｣ＮＺｯ［ｾ＠ tTl  Ｏ ｉｃｾ＠ 1 
ｾ［｡ｯＺ＠ N.::.ool In 
S'oIltch " P"":"E ﾷ＠ ﾷ ｾＺ ｲＧﾷｴｹｬｩ＠ ｾｌ｣ＺＡＭｅ､＠ ＭｴｩｦＮｾｲ＠ t:. ｾ ［ｏＧＨ｜Ｑｬｲｊ＠ ［Ｚｏ＠ Ｂ

J.G  l  . t.  ｾＮ＠ ｏ＠ '). II  ｾ＠ Ｎ＠

· SE- J :'J ＮｾＱＹＴｈ＠ l!)O . 1" . : al,1!)t:; . • ｨ｜ＮｬＺＡＬＩｾＧｮ＠ ｴＮｪＺＭＲ｜ｊＵｩｑｾ＠ Ｑ

S'_ . 1E :  Ｍ Ｎｾ＠ Ｎ＠ ｩｪｈｴＺｦｉＢ＠ l\:I'J .  lsa . -3uE':-66 Ｎ＠ ｏ＠ ＺＧｾ＠ ｊ Ｚ＠ ｬＹ ＺＱ＠ . :JJJ80 U 99
Ｑｾ＠ .OOl!J .. ＺＵ｟ＶｾＧＸＳ＠ 1 oQ .  ｬＮｺｾ＠ Ｎ＠ b14!1a3 -. n9SB109E 0'}887098
ｉｏｾ＠ . Ｇｊｏ＠ Ｎｾ＠ Ｚ ＧＺ＠ Ｎ＠ ＵＲｴＩＺＧﾣ Ｓ＠ ｾ＠ 10-.:,. (j< . i L o·,ool Ｍ Ｎ＠ ｮＺｲｾ＠ - ."}':r:::
Ｑ＠ Ｒｾ＠ . r·..'!.  . ":P161 ＱｾＮｴ＠ Ｎ＠ 11 . .... ＧｬｾＸＩＺＢ＠ ＮＢＧｬｊｾｕ＠ ｾ＠ . ( ｾ］［＠ ｾｕ＼［ｰ＠
Hi,;  . to ll r".. ｾ＠ lOti . 11'; . ,;)1";ClH Ｎ＠ ｃｈｻ［ＢＬ ｾＰＷＡＮＱ＠ Ｎｬ ＧＮ＠ ｾＢＧＮ＠ ＢｉｊｾＡ＠
. 7S nJI!I  t. . ＸＮ［ＺｾｬＭ＿ＧＡ＠ lOO . ｾＱｅＮ＠ Ｎ＠ ･ｓｬＱＰｦｾ＠ -. 'J;£:'':'12 ".).!I314,:-r.::.
10:'1 . ｾＲＮ＠ Ｗ＠ ＰＱＧＩｮ｢＠ Ｍ Ｎ＠ ［ｾ＿ｃ ［Ｔ＠ ＭＮ＠ ｾＺ［ＮｪＶ ［ Ｂ＠ ｴ＠
ｾｏ＠
•• <  Ｎ＠ ｴＡＧＩＬ
Ｎ＠ Ｇｾｏ

ｴ＠ ｾ＠
ｾ＠ Ｔ＠ ［ＧＹＸｃｲｾＱ＠

:: .1"1(1)'016 :"0 .  ill. ｾ ｴｉｾＢ＠ . fl1l21U Ｇ ｩＬ ｾ＠ i  Ｎｻ ｜ｪＺＮｾｯｃｮ［Ｒ＠

L' ' .,)1'1" Ｎ＠ ｾｬＢ＼ｏｲＮＺＧｩ＠ . q ut,m..·

.;. 1:'- ｾ＠ ＭＬ＠ Ｎ＠ ＺＡＧＬ［ｾ＠ •. ｾＬＱｃＸｉｪｊＴ＠
.u .H.! Ｎ＠ ｾ＠ ＧＴ＠ ｕ ｾＧｊ｟＠ ,. OU!T!61
!2: . 1)01;'j
. '  · 0 ｾ＠ o;r,
jIj'.i  . 63" ｟Ｎ＠ Ｚｮ ＮＺ＠ ｾＺ ｟Ｌ＠

(Ii!,. · ,.j' S H . :1'i. Ｇ＠ Ｎ＠ ＴｌＧＺｻＬﾣ ｾ＠ ｾ＠ ＮｈｾＬＺｕＧ＠ ﾷ＠ Ｍｩ＠

,"0 .. HI  Ｍ Ｎ＠ ｊｾｌ ..  ＺｾＨＡＮ＠ lrJe._.lliC . tU Ua-c .... .:cU·-4
ｑ［ｾ＠ Ｌ＠ Ｇ［ｈ ｾ＠ . '"'H e ｾ ﾷ･ｦｩＶ＠ 100 . ｾＺＩｮ＠ ＮＨＱ＠ Ｑ＠ :HjI - .': I}I)';H ' Ｎ［＠ ｾｬｃ＠ ＧｅＭｾ＠

ＢｾｉＬＮ＠ ＮｾＱＡ＠ Ｎ＠ ＬＩｏＡＭｾＯｪＢＷ＠ r.;. ＠ｓ Ｍｬ＠ Ｎ＠ ＧＺｾ［＠ Ｎ＠ '. . Ｇｬ［ｾ .. ＮＳ ｩＢｾＺＧ＠ !·"
, I 1 , ,. '" 
ｆＱｲ Ｎ
Ｎ＠ ｃＡｾ｜ｩ［＠ . 0:: :)':583;
ＧｬＱｾ＠ 1 ..... ｾｴ Ｂ ｬＭ Ｎ＠
.. ｾ＠ Ｉ＠ Ｎ＠ H . 41" t'll U .... .I .i P!i-:. Ｎ＠ ｜ＮＧＹＡ＠
l"L':FlL ｩｵＺ［ｰｾ＠ plt>l j.)u t :. ". ;:'4"
ｾｅ＠ Ｕ［
ＮＭｉｾ＠ Ｚｴ＠
;
f il •  . 
Jct"... Ju.:t11n'·j r:,l.:.1. ｾｩｮｲ＠ ｉ［Ｚ＠ ｾ＠ 1.:-11. ,:001 ｾＡｫ＠ !i:' .. .
i:"( . l : <-14;:£ .;, 1':-(. ｾ＼Ｚｩ＠ Ｎ＠ ＼ＧｪＨＮｯ［･ｲｾ Ａ＠ .;;;: H::: , Ｎ＠ ｴＡＧＺ＠ Ｒｑ＠ ｾ＠ ｾ＠

ｓｘ Ｇ Ｌ＠ ｲｾ＠ ＧＮｉ＠ ｾ＠ ｣Ｎ｟＠ ｾＩＭｩＮｉ＠ Ｎ＠ Ｇ ｜ Ｇｬｕｩ｡ＮＺ＠ ｾ＠ .!:.:·U!o ,; . azr. a5 fClr ｾ Ｎ＠ ｨ .. ｆ Ｚ＠ ＮＺｊｩｦｬＧｾ＠ !..l!:o.-II:!.*iJ ｬｾ＠ ..U!p ... .
·J&:':"At.-ll!!' fM..­:.u..,  ｾＢ＠ Ｌ＠ ｾｬ＠ .. ': Ｎ＠ ｃ
Tum.. .;.:i ｉｍ［Ｇ Ｎ＠ ｾ＠ • :E·4 ｾ＠ ':F. ­ ' 
'JAOet·':. 1Il.1r..4.ma  ­ f,n.  ［ｾ ＡｯｪＧＺ＠

­:!,a ....  cf .u.i:r: I'  ｬｾＺＮ＠

Bl6r1k ::_d ｴＮｮｕｾｩ｝＠ ,1:1 ｾｬｲＮ Ｎ＠ cards . ｾＱａｗ＠ .:;"'1'

36
'"1 

EI  anterior archivo de  resultados esta constituido de las siguientes partes: 

•   Un  cuadra  general  que  tiene  la  forma  del  listado.  Este  listado  se  divide  en  dos 
partes; la  parte  derecha  corresponde  a una  copia  de  las  IIneas,  tal  como  fueron 
editadas  en  el  archivo  de  datos  de  simulacion;  la  parte  izquierda  corresponde  a 
la  interpretacion  que  el  programa  Ie  da  a  cada  linea.  Esta  informacion  es  muy 
valiosa  cuando  se  esta  buscando  un  error en  un  archivo  de  datos.  Esta  forma 
del  cuadro  se  repite  en  forma  intercalada  en  diferentes partes  del  ｡ｲ｣ｨｩｶｾＬ＠ hasta 
completar todas las Ifneas del caso bajo simulacion . 
•   Despues  de  la  linea "BLANK  SOURCE"  aparece  la  tabla  de  conexionado  de  la 
red .  De  esta  tabla  de  conexionado  se  excluyen  los  acoples  capacitivos  e 
inductivos y las fuentes.  Los suiches sf aparecen  en  la  tabla  de conexionado de 
la  red. 
•   Despues de la  Ifnea "BLANK OUTPUT", aparece una tabla con  los valores de las 
variables electricas en  el  tiempo.  EI  ordenamiento de  las columnas en  esta tabla 
de  datos es la siguiente: 
­ Numero de pasos de integracion . 
­ Tiempo en segundos. 
­ Diferencias de potencial  que fueron solicitados en  las diferentes ramas. 
­ Voltajes nodales especificados en  lista. 
­ Corrientes de rama . 

Valores  de  maximos  y  de  mfnimos,  10  mismo  que  los  tiempos  en  que  se  presentan 
estos valores,  para  cada  una  de  las variables  y en  el  mismo orden  de la  tabla  a que 
se  refiere el  item anterior. 

Algunos  resultados  que  se  han  obtenido  con  el  programa  graficador  "PLOTXY"  se 
observan en  las figuras  11  y 12. 

37  
---i t ·

Corriente [A] 
ＧｉＭ Ｍ ＭｾＬＮ＠

[Al 

0.12 

0.06 

0.00 
II 

­0.06 
r
...  ....  •.•.  .  •.••. I ..  •  .  .  ....... ｾ＠ .. 

­0.12­t1 ­­­y­­­­;­­­­j­­­­t­­­­;­­­­­r­­­­;­­­­,­­­­;­­­­l 
o  4  8  12  16  20 
Tiempo [ms] 

Figura 11. Corriente en  el  circuito RLC  serie 

Voltaje  y corriente [V,  mAl 
ＱＸＰｾｲＭ｟Ｎ＠

....\.  HVoitajeeh capaq:rtor . 
120 

60 _'·H. 

o 
...  j  . 

­60  ··············i·  ............. j .•. 

ＭＱＲＰＫｔｲ Ｍｲ Ｍｾ｟Ｋｲ＠
o  2  4  6  8  10 
Tiempo [ms] 

Figura 12. Voltajes  y corrientes en  el  circuito RLC  serie 

38 
3.4  ANAuSIS DE RESULTADOS PARA EL CIRCUITO RLC

•   Calcule  para  el  circuito  simulado  la  impedancia  caracterfstica,  la  frecuencia 
natural de oscilacion  y el  factor de amortiguamiento que aparece asociado con  la 
corriente transitoria. 

•   Deducir la expresion  matematica para  el  voltaje sobre el capacitor. 

•   Utilizando  el  concepto  de  impedancia  caracterfstica  explique  porque  el  voltaje 
maximo  sobre  el  condensador se  aproxima  al  doble  del  valor  de  la  fuente, para 
un  circuito  RLC  serie con  poco amortiguamiento (R«Zc) .

•   Calcule  los  parametros  RLC  totales  para  una  linea  monofasica,  con  retorno  por 
tierra,  de  una  longitud de  100  km , que tiene un  conductor de 30  mm  de diametro 
y  una  resistencia AC  de  0.06  Q/km,  situado a una  altura  sobre  el  suelo  de  20  m. 
La  resistividad  del  terreno  es  de  100  nm.  Realizar  la  simulacion  para  una 
energizacion con  una fuente  DC  de  188 kVpico.  Observar el  voltaje al  final  de  la 
linea,  considerando  un  circuito  RLC  serie  y  un  circuito  PI.  Determinar  en  que 
tiempo el  voltaje esta en ±1 0%  del valor final. 

3.5  COMPLEMENTACION

Como  complementacion  a  la  parte  basica  se  realizaran  las  simulaciones  para  los 
siguientes casos: 

•   Circuito  RLC  serie  (150  n ,  163  mH,  0.7042  ).IF),  alimentado  con  una  fuente 
senoidal  de  frecuencia  industrial  de  60  Hz,  y  una  amplitud  de  100  Vpico.  Para 
este caso se desea observar el  voltaje sobre el  condensador  y compararlo con  el 
obtenido para  la  simulacion con excitacion escalon. 

•   Para  el  circuito  anterior  y  con  una  fuente AC  de magnitud  100 Vpico, empezar a 
cambiar  la  frecuencia  desde  60  Hz  hasta  1000  Hz.  Observar  la  magnitud  en 
estado  estacionario  del  voltaje  y  compararlo  con  la  magnitud  del  voltaje  de  la 
fuente.  Obtener una  grafica que relacione : 

Magnitud vo/taje condensador estado estacionario ｦＨｾ＠

I­ .  )
­ ­ ­ ­ ­ ­ ­ ­ ­ ­ ­ ­­ ­ ­ ­ ­ ­ ­ , ==  trecuencla
Magnitud votaje fuente I  

Explicar el  resultado 

39 
•  Realizar  la  misma  simulacion  inicial  con  fuente  escalon,  pero  agregando  una 
carga  resistiva de diferentes valores alrededor de  la  impedancia caracterfstica. 
EI  valor de  la  resistencia de carga  Rc sera: 

Rc = ·Iue =  Impedancia caracterfstica de sobretension 

Rc =  j O.163/0 .7042x10 "" =481 .11D

•  Observar  el  voltaje  Vc(t) y  comparar  con  el  obtenido  inicialmente  sin  esta 
resistencia de carga. 

•  Con  el  fin  de  introducir el  concepto de  propagacion  y de parametros distribuidos, 

se  partira  un  circuito  RLC  en  10  secciones  (0 .6  Q ,  1.27  mH , 0.09  /J.F,  para cada 
tramo  que  representa  aproximadamente  10  km  de  una  linea  de  alta  tension)  de 
modo que  se  aproxime  a la  manera  ffsica  real  como  estos  efectos  aparecen  en 
las lineas.  Este tipo de modelamiento de una linea se denomina semidistribuido. 

..  •  •  •  •  •  •  ..  ..  .. 

.  .. 
!  !  !  !  !  !  !  !  !  ! 
Figura 13. Circuito con  una cascada de elementos RLC 

EI  objetivo  es  observar el  voltaje  que  se  va  propagando  por  diferentes  puntos 
de  la  linea.  La  excitacion  recomendada  para  este caso es un  escalon,  aunque 
con  una  fuente  senoidal  tambiEm  se  puede  lograr el  objetivo  que se  persigue. 
Se  recomienda  trabajar  con  un  Tmax  de  0.4  ms  y  un  Deltat  de  0.4  /J.s,  para 
analizar unicamente esta parte del transitorio. 

40  
 ｛ｖ ｾ｝＠ Ｍ ＭＬ
Voltaje 250  I  ＭＺ［ Ｍ ＺＭ ｯＭ ｾ ｲ
200 

150 

100-1 1 
· f \  \:  Ｗ＼Ｚ＾ｾ＠ ｇｾ＠ ｾ ｾ 0+  ｌ ＮＭ ｾ＠

50 

0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8 

Tiempo [ms] 

Figura 14. Resultado esperado de la simulaci6n con  elementos semidistribuidos 

Dar respuestas a las siguientes preguntas : 

(.Que  retrasos en  el  tiempo tienen  los diferentes voltajes en  cada  punto de 

la  linea con  respecto a la  fuente? 
Observar el  voltaje  a  mitad  de  la  linea  y  compararlo  con  el  voltaje  al  final 
de  la  misma . 
­ Observar  el  orden  en  que  los  voltajes  van  al  100  %  con  respecto  a  la 
fuente  y  el  orden  en  que los  voltajes van  al  200%  con  respecto a la  fuente . 
Explicar el  fen6meno . 

•  Si  se tiene  un  equ ivalente de cortocircuito en  una  barra de 220  kV 

Voltaje de  prefalla :  226 .23  kV  L O.O°

Corriente de falla  trifasica : 11872 A 
Corriente de falla  monofasica: 14771  A 

(X / ｒＩ Ｓｾ＠ =6.05 

(X / ｒＩ Ｑｾ＠ =6.44 

41 
Las  bases son  100 MVA Y 220  kV 

Donde  (X I R) 3q, ' es  la  relaci6n  XIR (reactancia  inductiva/resistencia)  que  ve  la 
corriente  de  corto  trifasica  y  (X I R) 1q, es  la  relaci6n  XIR que  ve  la  corriente  de 
corto monofasica 

Realizar la  simulaci6n  para  determinar, cuando ocurra  una  falla  trifasica  y  una falla 

monofasica  en  esta  barra ,  el  valor  de  la  corriente  dinamica.  Recordar  que  la 
corriente  dina mica  es  el  pico  maximo  de  la  corriente  que  se  puede  presentar 
durante  una  falla  y  depende  del  momento  en  que  se  presente  la  falla  y  de  la 
relaci6n  XIR. Investigar  la  relaci6n  que  propone  la  norma  IEC­909  para  calcular el 
valor de esta corriente a partir de la  relaci6n  XIR.

La  teorfa  correspondiente  al  equivalente  de  Thevenin  trifasico  se  puede  consultar 
en  el  Anexo A. 

42  
" 

4.  pRACTICA N02.  ONDAS VIAJERAS 

4.1  OBJETIVO 

EI  objetivo  en  esta  practica  sobre  onda  viajera  es  entender  la  solucion  obten ida 
analfticamente  para  la  onda  viajera  en  Ifneas  de  transmisi6n  y  darle  la 
correspondiente interpretacion fisica . 

A  partir  de  la  solucion  para  linea  ideal  de  la  ecuaci6n  del  telegrafista  se  van  a 
confirmar, mediante simulaci6n , los  resultados obtenidos de la  solucion analftica . 

4.2  MARCO TEORICO 

Los  modelos de onda viajera  para  lineas y cables se  dividen en forma  general: 

•  Modelo para  linea trifasica transpuesta  (Modelo de Clarke) . 

•  Modelo para  linea no transpuesta  (Modelo K.C . LEE) 

EI  modelo de Clarke utiliza en  los calculos  una  matriz de transformaci6n con stante, 

ya  que  para  la  linea  transpuesta  los  parametros  en  los  modos  de  propagaci6n 
siempre se  pueden evaluar de la  misma  manera. 

La  matrices  de  transformacion  10  que  hacen ,  basicamente,  es  descomponer  una 
propagaci6n  que  se  daria  por  una  linea  acoplada  (caso  trifasico  por  ejemplo)  en 
varios  mod os  de  propagaci6n .  Para  el  caso  de  la  linea  trifasica  transpuesta,  los 
parametros  para  los  tres  modos  de  propagaci6n  se  pueden  obtener a  partir de  las 
componentes  de  secuencia .  Cuando  la  linea  deja  de  ser  transpuesta  los  modos 
se  evaluan  para  cada  caso  en  forma  particular.  Para  este  caso  se  debe  calcular 
una  matriz  de  transformaci6n  particular,  que  se  hace  normal mente  con  la  opci6n 
LINE  CONSTANTS del ATP . 

Si  la  linea  es  monofasica  no  se  requiere  una  matriz  de  transformaci6n  y 
adicionalmente,  si  se  considera  el  caso  ideal  sin  perdidas ,  la  soluci6n  seria  la 
obtenida  para  la  ecuaci6n del telegrafista: 

43 
 v, (t ­ x 1v)u(t ­ x 1v) +
r Z v, (t ­ (2d ­ x)1 v)u(t ­ (2d ­ x)1 v) +
V(x,t) =  Zc  xl r rZv,(t-(2d+x)lv)u(t-(2d+x)lv)+
J
(4.1 ) 
Zl  +ZC r)rZZv,(t-(4d-x)lv)u(t-(4d - x)lv)+
r)zr} v,(t - (4d +x)lv)u(t -(4d +x)lv)+ ...

v, (t ­ x/v)u(t ­ x/v)
- r2 V, (t ­ (2d ­ x )/v )u(t ­ (2d ­ x)lv)
1 x  +r) rzv, (t ­ (2d + x )/v )u(t ­ (2d + x )/v) (4 .2) 
i(x,t) = Z  +Z 
J  c 
­ r,rzzv, (t ­ (4d ­ x)lv )u(t ­ (4d ­ x )/v)
+r,ZrZ2 v, (t ­ (4d + x )/v)u(t ­ (4d + x )/v)- .. .

EI  procedimiento seguido  para obtener i(x,t) y  v(x ,t) y  el  significado de cad a  uno de 

los terminos  presentes en  (4 .1)  Y (4.2),  se presentan  en  el  ANEXO  B. 

Cuando  la  linea  es  del  tipo  semi­infinita  las  ondas  reflejadas  del  lado  de  la  carga 
nunca  Ilegan  y  en  las  ecuaciones  anteriores  solamente  aparecen  los  primeros 
terminos de la serie: 

Z
V(x,t) =  C x  [v,(t-x/v)u(t-x/v)] (4.3)
Z j +Zc

i(x,t) = _  1_ _x [v, (t ­ xlv )u(t ­ xlv )] (4.4)

Zl +Zc

Lo  mismo  sucede  cuando  la  linea tiene  una  carga  resistiva  igual  a  su  impedancia 
caracterfstica . 

44  
4.3  SIMULACIONES 

Se  realizaran  una  serie  de  simulaciones,  empezando  por  IIneas  ideales  y  luego 
con  Ifneas  reales  del  sistema  de  transmisi6n  colombiano .  Para  las  Ifneas  ideales, 
los  parametros  se  pueden  estimar  mediante  calculos  manuales,  pero  para  lineas 
reales  se  usa  el  subprograma  LINE  CONSTANTS  para  que  determine  los 
para metros  correspondientes  a  los  modelos  de  onda  viajera.  Este  subprograma 
hace  parte del  programa ATP. 

4.3.1  Energizacion  de  una  linea  monofasica  sin  perdidas,  en  vacio,  con  una 
fuente impulso de voltaje. 

Esta  primera simulaci6n se realizara considerando una  linea ideal  sin  perdidas.  La 

linea se  divide en  dos tramos de 25  km  para efectos,  unicamente,  de  determinar 10 
que Ie  sucederfa a un  observador que este situado en  la  mitad  de  la  linea.  EI 
dielectrico es el  vado. 

Longitud  =50  km 

Impedancia caracterfstica  =500  n 
Velocidad de  propagaci6n  = 3*10 5  km/s

Tmax  =2.5 ms 
ｾｴ＠ =0.15  ).ls 

vr 
Figu ra  15.  Linea energizada con un impulso de vottaje 

La  fuente  de  voltaje  a considerar es  una  de  tipo  impulso,  de  8 ).ls  de  frente  y  una 
constante de tiempo de cola  de  20  ).ls  (Ver figura  16) 

45 
 ­ 
2robO $; & 3-p he.ue

a,..,ce, Linger
Bre nch .tionline Gl"
Un"  Lump e d 
U n"  QlsttlbObOd 

ｓｬＡ ･ｨ ･ｾ＠
­. 
Q C lVpa  11 
MAchine. Bftm p type  12  
1llop R­Rmnp type 1 J  
T (e n5fcrme,..
AC typ.. 14  
MODELS 
Type l!.  
IACS  
SUIU'>'P e  l S  

:IACS . ouree 'If-@-

ｾＢＧ Ｑ＠
r_
-,
ｾｲｓ ｐ Ｂ｣ｩｨ･ ､＠ AC;Ph .type Ｑｾ r
Llne/Ceble
ａ ｾ ｕｮ ｧ ｲｯ ｵｮ､Ｘ
Overh....d Una (PCH)
ｾ＠
D C !.Lngro un d gd QI< IjoIo
Ereq.uancy comp o 

Figura 16. Ventana de datos, fuente tipo impulso de voltaje

ｾ＠
Lo que se pretende con esta simulaci6n es que se entienda la soluci6n dada a la
propagaci6n de una onda viajera planteada en las ecuaciones 4.1 y 4.2. La
energizaci6n se realiza con una funci6n impulso, para que los terminos de la
soluci6n a la ecuaci6n de onda aparezcan graficamente en forma separada .

Se consideraran diferentes situaciones de impedancia de la fuente y de

impedancia de la carga.

4.3.2 Energizacion de linea monofilsica sin perdidas, en vacio

Para este caso se utiliza el modelo monofasico de Clarke. Se selecciona de

acuerdo con la figura 17.

Erob.. &ｾｨ･ ｳ ･＠

Jl.rcndlUn. ", 
eotn bute. ,

61Mdl l:ionlln"o,
1.Q£b. 
I.ine Lumped
RJI 

ｾ ･Ｎ＠
Untmn. p. Tone. {KC1.nl  '   Z ph . ...  c::=:::::::>  IA  500

!!.  JE5 
J ph"".  I 100
Sources
ｾｰｨｍｴｬ＠
Mp,ne..  IUNE
6 ph mutu. 1 
T (lInaftumo .. apha,. 
hjOOELS
Type 94 
lACS 
\Jle, Spoc:msd 
QIIerto.ed u n. (PCH) 
Freqo8ncy co mp o 

Figura 17 . Modelo de Clarke monofasico para onda viajera

46 
--

Para entrar los datos de la linea hay varias opciones que dependen de la opcion
ILlNE, esta puede tomar los valores 0, 1 0 2, Para comprender estas opciones se
puede consultar la ventana del HELP (Ver figura 18),

Itlame  lINEZT_l  ­ Distributed  parameters ,  sigle  phase 

Card  BRANCH 
Data  R/  l=  Reslstance  pr o  length  in  (Ohm/  leng th] 
A=  I1INE=0 '  A=inductance  I'  in  (mH/  lengthJ  if  Xopt  =0 
A=I'  in  (ohm /  length]  if  Xopt . =power  frequency . 
I1INE=1  A=Modal  s urge  impedance  in  (ohm]  Z =sqrt ( L '/C') 
IIINE=2  AcModal  surge  impedance  in  [ohm]  Z=sqrt(L'/C') 
B=  I1INE=O  B=capacitance  C'  in  Ｈｾｆ Ｏ ｬ･ｮｧｴｨ｝＠ if  Copt . =O 
B=C'  in  ｛ｾｍｨｯ Ｏ ｬ･ｮｧｴｨ｝＠ if  Copt . =power  frequency . 
IIINE=l  B=propagation  velocity  in  [length/sec . )  v=l/sqrt(C'*L') 
I1INE=2  B=modal  travel  time  in  (sec]  T=sqrt ( L'*C')*length 
1=  length  of  line  (>0  for  transposed  lines) 
IIINE=  Takes  values  from  0  to  2 
Xopt  and  Cop t  .  are  set  in  menu .  ATPISettings/s imulation . 
Node  From=  Start  node  of  line 
To=  End  node  of  line 
On ly  branch  voltage  output  is  reliable . 
RuleBook :  IV . D . l 

Ｉ＠ ｾ . .i
.!l

Figura 18, Ventana con la ayuda para entrada de datos del modelo de Clarke

Los para metros que pide el modele son los siguientes:

• Resistencia de perdidas (Q/unidad de longitud)

• Impedancia caracterfstica (Q para ILlNE=1)

• Velocidad de propagacion (unidad de longitud/segundos para ILlNE=1)

• Distancia del tramo de red (unidad de longitud)

• Variable de control de unidades ILiNE. Para el caso de las definiciones

anteriores esta variable debe tener valor 1

Nodo.-,A Nodo_B

ｾ＠

T
Figura 19, Linea ideal energizada con fuente escalon de voltaje

47
.  " 

Para  la  simulaci6n ,  implementar  el  circuito  de  la  figura  19.  La  energizaci6n  se 
hara , inicialmente , con  una fuente constante y para  linea ideal  sin  perdidas. 

Una  precauci6n  importante en  las  simulaciones  que se  hagan  con  los  modelos de 
onda viajera es la de que la  variable  DELTAT (delta de tiempo de simulaci6n)  debe 
ser menor que el  tiempo de viaje de  la  linea  mas corta .  Para  el  caso de la  anterior 
linea , que tiene  una velocidad  de  propagaci6n  de 300 .000  km/s y una  distancia  de 
100 km: 

,.  distancia 100  333  33

DELTAT maximo ｾ＠ =  = .   )J.s
velocidad 300000

Se podrfa tomar un  DEL TAT de 10 )J.s. 

Con  el  fin  de determinar el  tiempo  de  simulaci6n  (TMAX)  se  debe tener en  cuenta 
la  duraci6n  del  periodo de la  oscilaci6n  de  voltaje  al  final  de  la  linea , que esta  del 
orden  de  cuatro  veces  el  tiempo  de  viaje,  y considerando  que se  quiere  realizar la 
simulaci6n  durante  10  periodos: 

TMAX = 10x 4 x 333 .33 x l O-6 =1 3.3 ms 

Es  bueno chequear el  numero de pasos de calculo que tendria  la  simulaci6n : 

,  3 
' d d 'l 1 TMAX
=  13.3 x l O­6  = 1 3 30  pasos
Numero e pasos e ca  cu 0  = ­ ­
DELTAT 10 x l O­

Este  chequeo  se  realiza  porque  el  numero  de  pas os  pod ria  dar  lugar  a  una 
simulaci6n  muy  pesada  y  habria  que  afectar,  adecuadamente ,  los  "settings" 
enteros.  Recordar que la  variable "Plot Freq"  de los "settings" enteros debe ser un 
numero entero  impar. 

48  
 4.3.3 Energizaci6n de linea monofasica sin perdidas, con carga

1" 
.. 
Nodo_A 
ｾ Ｎ＠
Nodo_B 
,f 
Figura 20 .  Linea ideal energizada con carga 

Energizar nuevamente la  linea  pero considerando una carga  resistiva  con  los 

siguientes valores: 
1  . 

Rcarga de  1000,  100 Y 500  Q 

Para  cada  caso,  calcular  el  coeficiente  de  reflexi6n  y  compararlo  con  el  tipo  de 
oscilaci6n  que se  presenta . 

4.3.4 Energizaci6n de linea monofasica sin perdidas, con resistencia de

preinserci6n

T 
Figura 21. Linea ideal con  resistencia de  preinserci6n 

Una  manera  de  controlar la  magnitud  del  sobrevoltaje  que  se  presenta  al  extremo 
de  una  linea durante su  energizaci6n  es  mediante  una  resistencia  de  preinserci6n. 
EI  objetivo  de  esta  resistencia  es  hacer  que  la  onda  viajera  que  regresa  hacia  la 
fuente yea  un  coeficiente de  reflexi6n  muy  pequeno .  Considerar el  sistema que se 
ilustra  en  la  figura  21 .  A  la  resistencia  de  preinserci6n  para  efectos  de  la 
simulaci6n  se  Ie deben controlar dos  parametros basicos : 

•  Valor de la  resistencia 

49 
:1  

•  Tiempo que debe permanecer conectada en  serie con  la  linea. 

Este par de parametros se  deben escoger y justificar adecuadamente de acuerdo 
con  el  sistema  bajo simulaci6n . 

4.3.5 Obtenci6n de modelos Lee y Clarke con el subprograma del ATP, LINE
CONSTANTS.

Los  modelos  de  onda  viajera  para  linea  transpuesta  y  no  transpuesta  son  los 
modelos  mas  utilizados  para  representar  adecuadamente  una  red  de  transmisi6n 
de energia electrica  en  diferentes estudios de  un  sistema de potencia . 

Para  estudio  de  transitorios  y  considerando  una  linea  trifasica  se  hace  una 
descomposici6n  en  modos  de  propagaci6n .  Lo  anterior  es  semejante  a  la 
descomposici6n  que  se  hace  en  redes  de  secuencia  para  estudios  de  estado 
estacionario.  Cuando  la  linea  es  transpuesta  y  los  transitorios  a  estudiar  no  son 
de  alta  frecuencia ,  los  modos  de  propagaci6n  se  pueden  obtener  a  partir  de  la 
teorfa  basica de impedancias de secuencia . 

Para  la  practica ,  se  van  a  calcular  los  modos  de  propagaci6n  a  partir  de  los 
parametros de secuencia. 

Rmodo =  Rsecuencia (Q/km) (4 .5) 

fTr7ciUCiBnCi a de secuencia (H ) (0)

(4.6)
Zcaracteristica  del modo =\1  Capacitancia de secuencia (F)

1
Vel. de prop. = ­ (km/s) (4 .7) 
,jlnduc de sec. (H / km) x Capac de sec.(Flkm )

De  manera  analoga  a  como  se  obtienen  los  parametros  de  secuencia,  para  los 
modos de  propagaci6n  se  consiguen tres tipos, donde el  primero es  el  modo tierra
(ground)  y  tienen  influencia  las  caracteristicas  del  suelo,  los  dos  siguientes  son 
iguales entre sf  y  no tienen  influencia  del  suelo. 

La  formulaci6n  basica  para  el  calculo  de  los  para metros  de  secuencia ,  de  una 
linea con  dos cables de guarda [8] , es: 

50 
'  :1 

2 
6  Z'g
Z  =  (Zs  +  2  ZM)­ Z  +  Zgg (4 .8)
o  9 

=
Zs  + 2Zm  Rae + (3 x  0 .0592) + j3 x  wp In[  Oe  ]  [.Q/km] (4 .9) 
2 TT ｾ RMG' . DMG;,

r­­:  -
De = 658,86 .. P /  f  [m] (4 . 10) 

Z,g = 0 .0592 + j wJ.1

2 TT
ｾｊ＠
In(DMG,g [Q/km] (4 .11 )

Z99 =0.0592 + j wJ.1 In[De J [Dlkm] (4.12) 

2 TT Ogg

Zg = Rae + 0.0592  + j wJ.1

2TT
In(. RMG'
De )  [Q/km] (4 .13)

­ In ( DMG
Z]  =Rae + j WJ.1 J Q/km ff
(4 .14)
2TT RMG'

Donde, 

Zs: es la  impedancia  propia del conductor de  fase . 

ZM es  la  impedancia mutua entre los conductores de fase. 

Zfg: es la  impedancia mutua entre los cables de guarda  y los  de fase . 

Zg: es  la  impedancia propia de los cables de guarda . 

Zgg: es  la  impedancia mutua entre los  cables de guarda . 

51 
,  Ii ___ n'  T 

Rae: Resistencia AC del conductor (Qlkm) 

w : es  la  frecuencia angular de la  corriente ( rad/s) 

J.1: Es  la  permeabilidad  magnetica del  medio. 

p: es la  resistividad  del  suelo  (Om)  

f  es  la  frecuencia  de la  corriente  (Hz)  

De: es  la  distancia  media  geometrica  entre  las  corrientes  de  los  conductores  de  
fase  y  sus respectivas  imagenes. 

oMG: es  la  distancia  media geometrica entre fases 

RMG ': es  el  radio  medio  geometrico  del  conductor.  Dato  que  normal mente  se 
obtiene de  las tablas de fabricantes de conductores.  

oMG fg: es  la  distancia  media  geometrica  entre  los  conductores  de  fase  y  los  
cables de guarda .  

o gg : es  la  distancia  media  geometrica entre los cables de guarda. 

Z1: es  la  impedancia de secuencia  positiva de la  If nea . 

Matematicamente  una  distancia  media  geometrica  (DMG)  entre  un  grupo  de 
elementos de un  conjunto  con  otro grupo de elementos de otro conjunto, se define 
como  la  rafz  n­esima  de  todas  las  distancias  posibles,  entre  cada  uno  de  los 
elementos del  primer conjunto con  los elementos del segundo conjunto . 

Para  el  calculo  de  las  capacitancias  se  pueden  utilizar  las  siguientes  expresiones 
[8] , 
2m: 55.55 nF
(4.15)
C[  =  oMG oMG km
In  In ­
r, r,

52 
: I  

1  2TT£
C  = ­x- (4 . 16) 
o  3 

ｾ＠
Ｒ

In  OMG ,aseS-imilgeneS guardas

In  ｏｍｇＧ ｾ ｡ｳＭｩｭｴｊｧｮ･ｓ＠ _ 
[ 
OMG,aseS- gUardaS

J J{,. OMG;' In  OMGguardas - imagenesguardas

7J { g . 0 guardas

Donde ,  

C1  :  es  la  capacitancia de secuencia positiva de la  linea.  

Co  : es  la  capacitancia de secuencia cero de la  linea.  

€  :  es  la  permitividad  electrica del  medio.  

{f :  es el  radio  de los conductores de fase . 

{g  :  es el  radio de los cables de guarda.  

Oguardas:  es la  distancia entre los cables de guarda.  

OMG fases-imagenes : es  la  distancia  media  geometrica entre  los  conductores  de fase  y  
sus imagenes. 

es  la  distancia  media  geometrica  entre  los  conductores  de 

OMGfases-imagenes gua rdas:
fase  y  las  imagenes de los guardas. 

OMG guardas-imagenes guardas: es  la  distancia  media  geometrica  entre  los  cables  de 
guarda  y sus imagenes.  

Realizar  los  calculos  de  los  modos  de  propagacion  para  una  linea  de  500  kV,  la  
cual  tiene la  siguiente configuracion .  

53  
:/ 

­ 25m 
­ i 

ｉｴＭＧｾ＠
\/'*'rl'" 
ｾＭＢ＠ ! 
• 
46cm 
­ 43m 

33  m 

, 
Figura 22 .  Configuraci6n  linea de 500  kV 

En  el  preinforme  deben  aparecer  los  calculos  para  estos  modos  de  propagaci6n , 
de  manera  que se  confronten  con  los  que se  obtengan  mediante el  modelo que  se 
genere con  el  programa  ATP.  Hacer el  calculo  para  frecuencias  de  60,  100,  500  Y 
1000  Hz,  y  para  una  resistividad  del  terreno de  300  Om. 

Los  datos ffsicos  del  conductor y  del  cable de guarda  son : 

Conductor  Cable de guarda 

Diametro  25.17  mm  8.71  mm 

Rdc  0.089120/km 1.8440/km

Considerar para el  cable de guarda  la  conexi6n  tradicional ; es  decir,  continuo  y 

aterrado en  cada  estructura . 

54 
Los  resultados que se obtengan sirven  para alimentar la  ventana de datos del 
modelo, que en  el  programa se  llama  modelo de  Clarke (3  fases) . 

frobes 8.  11lha.e 

!l.rencll Li'iG'"  
firencll /:iOnii"ear  

­
ｬｰ ｾ ｡ ｳｑ＠

apha. 8 
ｾＡ ｉ ｜ｕ ｴｦ ｾ＠

lOourcs.  
MJ!cI1  netl 
T[allsiorm9f3 
i 

ｾ
ｾ phas. 
6 ph mulu81  
ph"••  
D
MODELS 
ｔ ｹｰ ･ｾ＠
IACS 
,  nml  ll ""I'I'  I Jr. ,  \  ,. 
).i.e. Sped1led  ｾ ｢ ｵ ｴ･ ｳ＠ I 
L,no/Qabl • 
.Q)ieme" d Line  (?CH)  r.:I"
DA .. ­ •. ­. :­｜ｬ ＺＭ ｾ
T:­:A­­
.  ' ＺＭｊ ［Ｚｅ ］ Ｚ ＭＧ ｩｬ＠ ｉｎｏｄ ｾ＠ JpHASE  1-
f requency comp o  ｾ _ _  0  ­.J  ｾ＠ J  I 
.Rl19_  _   0tJT1_  .3  1 
A·  51lO 
, 500 
300000 000 
j ＱＰ Ｐ Ｐ Ｐ Ｐ ｾｏ＠

; 5r11 
IlJNE  ­.:J 
G.IOUp  No: 10   LoII.1 

Cort,ment:, f
Output   .­ r  H,d_ 
(""  CUrrlln!  r  ','ott"'l"  (""  C\(rr&Voli  r  eow•• &E",,'IlY  , r 
ｲ＠ ｾ •. l 
V

ｾ＠   Qt(  ｾｲＮ｣ｯ ｬ＠ ｴＺＡ Ｎ ｾ＠ II 

Figura 23.  Ventana de datos para el  modelo de Clarke trifasico 

Consultar en  el  Help la  definici6n de cada  uno de estos parametros: 

Name  LIN EZT  3  ­ Dls trl bu t ed  p a r ameters  (Clark e )  3  pha s e .  Trans p osed 
Card  BRANC H 
Data   R/ l +=  Res i sta n ce  pr o  l e n g th  in  [Ohm / l e n gth ) .  PO S .  seq ue n c e. 
R/IO­ Res is t ance  p ro  l e n gth  in  [Ohm/le ng t hl  .  Zero  s eq uenc e . 
(A:  ILI NE=O:  A=induct an c e  L'  I n  [rnH/l eng t h]  if  Xo p t.= O 
(  A=L '  i  n  [ ohm/ l e n gthj  i  f  ;.:o pt .= p O\.Je r  f requ en c y. 
I LI NE=l  :  A=Modal  s ur g e  i mp edanc e  i n  ( o hm]  Z= sqrt (L ' / C'  ) 
IL IN E=2:  A=Modal  su rge  lmped a nc e  l  D  [ o hm ]  Z=s q r t  ( L '/C ' ) 
A+  Po s  .  s equence.  
AO  Ze ro  seq u ence.  
( B :  I LI NE =O:  B=capaclt a n c e  C '  i n  ｛ｾ ｆＯ ｬ ･ｮｧｴｨｊ＠ I  f  Cop t  . =O 
(  B=C'  In  ｛ ｾ ｯＯｬ･ ｮ ｧ ｴｨｊ＠ i  f  Co p t. =p owe r  frequ enc y. 
{  I LI NE=l :  B=p rop a g a tlo n  v e lo c l t y  i n  (len g th / sec . )  v=l/ sq r t {C ' *L') 
f  I LI NE=2:  B=moda l  tra v el  t i me  in  [ sec )  T= s q r t  ( L "C') * le n gth 
x o p t  .  a nd  Co p t  .  is  s e t  i n  me n u:  ATPIS e t t i n g s/ S imu la t io n . 

55  
 B+  Pos.  sequence. 
8 0  Zero  sequence. 
1  length  of  line  (>0  for  transposed  lines) 
I LINE%  Takes  values  f rom  0  to  2 
IPUNCH:  O:  lumped ­ resistive  mode ll ing.  G ' :Q 
%1:  distor t ion l ess­mode  mod elli ng.  G'%R ' *C'/ L ' 
Node  I N1=  Sta r t  node .  Letter  'A',  ' 8 '  or  ' c ' 
is  added  to  nodename  for  node  1,  2  and  3 . 
OUT1%  End  node.  Let ter  'A ',  ' B'  or  'C' 
lS  added  to  n ode name  for  node  1 ,  2  and  3. 

Consultar  la  opcion  LINE  CONSTANTS  del  ATP  con  el  fin  de  entender  las 
diferentes  opciones  que  tiene  este  subprograma .  A  partir  de  la  version  2.0,  del 
programa  ATPDRAW,  los  datos  de  la  configuracion  ffsica  de  la  linea  se  pueden 
entrar  directamente  mediante  la  ventana  de  datos.  Esta  opcion  corresponde  a 
LINE/CABLE  de  la  caja de  herramientas. 

Lo  que se  puede obtener a partir de esta opcion , que ofrece la  interfaz ATPDRAW, 

es  una  corrida  del  subprograma  LINE  CONSTANTS  para  generar  modelos  de 
Ifneas.  Si  se  quieren  resultados , como  matrices  de  impedancias  y capacitancias, 
parametros  de  secuencia ,  etc. ,  se  debe  utilizar  la  interfaz  LCC ,  que  es  un 
programa independiente del ATPDRAW. 

Para esta linea se  requiere obtener la  siguiente informacion : 

•   EI  modelo  de  onda  viajera  trifasico  para  linea  transpuesta.  De  acuerdo con  el 
formato  para  el  modelo  de  CLARKE ,  leer  los  valores  de  resistencia  de 
perdidas ,  impedancia  caracteristica ,  velocidad  de  propagaci6n  para  los 
diferentes  modos  de  propagacion.  Para  el  caso  de  linea  trifasica  son  tres  (el 
concepto  es  equivalente  a  una  descomposicion  en  redes  de  secuencia) .  La 
longitud  de  linea  es  de  209  km .  En  la  obtencion  de  estos  modelos  es  muy 
importante  la  frecuencia ,  que  debe  ser  muy  aproximada  a  la  real ,  cuando 
ocurra  el  transitorio.  Suponer  una  resistividad  del  terreno  de  100  Qm,  y  una 
velocidad  de  propagacion  igual  a la  de  la  luz  para  el  calculo de  los  parametros 
(esto  ultimo equivale a un  calculo que se  debe hacer de la frecuencia) . 

•   En  el  modelo  de  Clarke ,  el  programa  ATP  utiliza  una  transformacion  interna  _ 
que  tiene .  Investigar  cual  es  dicha  transformacion  y  porque  no  es  valido 
utilizar la  transformaci6n de componentes simetricas . 

•   Determinar  la  sensibilidad  de  los  parametros:  Resistencia  de  perdidas,  
impedancia  caracterfstica  y  velocidad  de  propagacion ,  con  respecto  a  la  

56 
 variacion  de  la  frecuencia  (10Hz  :$  f  :$  5000  Hz)  y  la  resistividad  del  terreno 
(100­ 3000  Qm) .  Estas  sensibilidades  se  deben  presentar  en  forma  tabular  y 
en  forma  grafica. 

•   Determinar  la  relacion  entre  Zc  Y RI4 (R es  la  resistencia  de  perdidas  total  de 
la  linea)  con  el  fin  de  verificar  la  necesidad  de  considerar  tramos  mas  cortos 
de  linea  en  la  modelacion.  Se  debe  cumplir  que  Zc»RI4 para  garantizar  un 
funcionamiento  adecuado  de  estos  model os  de  onda  viajera.  Justificar  la 
relaciones  entre  las  resistencias  de  perdidas ,  impedancias  caracteristicas  y 
velocidades de  propagacion  de los  dos modos  de  propagacion  para  una  linea 
transpuesta . 

•   La  misma  informacion  del  primer  item,  pero  considerando  linea  no 
transpuesta.  Para  este caso,  aparece una  matriz de transformacion  diferente a 
la  de Clarke,  que es calculada para  cada  configuracion de linea . 

•   Realizar  una  simulacion  energizando  una  sola  fase  con  una  fuente  tipo 
escalon  (identificada  como  fuente  DC  tipo  11  en  el  ATP) .  Determinar  la 
frecuencia  real  de  la  oscilacion  del  voltaje  al  final  de  la  linea  en  la  fase 
energizada.  Una  vez  determinada  la  frecuencia  real  de  oscilacion  del  voltaje, 
cambiar  la  fuente  ｰｾｲ＠ una  del  tipo  sinusoidal  (identificada  como  fuente AC  tipo 
14  en  el  ATP),  y  observar  nuevamente  el  voltaje  al  final  de  la  linea, pero  esta 
vez  interesa  observar el  resultado  cuando  el  voltaje  se  estabilice .  Justificar el 
resultado.  GCual  serfa  la  longitud de linea  para que en  la  anterior fuente no  se 
presente  un  sobrevoltaje en  estado estacionario  por encima  del  10%?  A  partir 
de  los  resultados  obtenidos  dar  una  posible  justificacion  de  la  frecuencia 
comercial  de 60  Hz. 

•   Energizar  la  linea  transpuesta  con  una  fuente  trifasica  sinusoidal, 
considerando  diferentes  tiempos  de  cierre  de  los  interruptores.  Determinar  la 
condicion  mas  crftica  en  los  anteriores  tiempos ,  tomando  como  serial  de 
criticidad  el  maximo sobrevoltaje al  final  de  la  linea en  cualquiera de las fases. 

57  
 5.  PRAcTICA N°3.  FENOMENOS TRANSITORIOS ORIGINADOS POR UNA
CORRIENTE DE RA YO

5.1  INTRODUCCION

En  esta  practica,  que  hace  enfasis  sobre  el  efecto  de  las  descargas  atmosfericas 
en  lineas  de  transmision ,  se  estudiara  el  efecto  del  rayo  sobre  el  aislamiento  de 
una  linea.  Se  pretende,  ademas  de  estudiar  la  parte  fisica  involucrada  en  la 
interaccion  del  rayo  con  la  linea ,  obtener  algunas  cifras  que  midan  el 
comportamiento  del  aislamiento  de  la  linea  frente  a  un  rayo  de  determinadas 
caracteristicas.  EI  parametro  que  clasicamente  ha  medido  dicho  comportamiento 
es  el  numero  de  salidas  de  la  linea  par  cada  100  km  de  longitud  y  por  cada  ario, 
parametro  que  constituye  un  criterio  de  diserio  del  aislamiento,  0  una  medida  del 
desemperio para  una  linea que ya  este construida . 

EI  rayo  es  un  fenomeno  que  origina  una  sobretension  de  origen  externo  0 
simplemente de origen  atmosferico. 

Las  sobretensiones  de  origen  atmosferico  son  independientes  de  los  voltajes  de 
operacion  de  los  sistemas  de  potencia , a  diferencia  de  las  que  se  presentan  por 
maniobra ,  que  si  dependen  del  nivel  de  tension  nominal  de  los  sistemas.  La 
importancia  relativa  de  las  sobretensiones  atmosfericas  crece ,  a  medida  que 
disminuye el  nivel de tension  de  los sistemas. 

Los  sistemas  de  transmision  se  diserian  normalmente  con  un  apantallamiento 
perfecto  0  con  una  probabilidad  muy  pequeria  de  falla  del  mismo.  Las  descargas 
atmosfericas  que  caerfan  en  forma  directa  serian  de  muy  poca  probabilidad .  EI 
flameo  inverso ,  ocasionado  par  una  descarga  que  cae  sobre  la  torre  0  sobre  el 
cable  de  guarda ,  es  la  causa  mas  comun  de  salida  de  la  linea.  Se  denomina 
flameo  inverso ,  porque  este  se  presenta  del  neutro  hacia  la  fase ,  debido  al 
aumento del  voltaje en  el  mismo. 

58  
5.2 MODELAMIENTO DEL RA YO

Una  nube  cargada  se  puede  representar  circuitalmente  como  un  condensador. 
Para  efectos  de  calculos  aproximados  una  nube  se  puede  aproximar 
geometricamente a  una  placa  circular, paralela  a la  superficie  del  suelo . Para  esta 
aproximaci6n se  utiliza la  ecuaci6n: 

c =  Eo x  A _  Eo x rrr 2 
d - - d
F
donde  Co  =  8.854 x  10­ 12 ­ (5.1 ) 
m 

d 
c 

Figura 24.  Representacion de la  nube  y canal de descarga 

EI  canal  de  la  descarga  se  comporta  como  un  hila  delgado , de  aproximadamente 
dos  mm  de  diametro,  pero  debido  al  intenso  campo  electrico  que  se  origina ,  el 
efecto  corona  hace  que  la  onda  de  corriente  de  descarga  adopte  un  diametro 
mayor,  que  para  efectos  de  calculo  se  puede  suponer  en  unos  20  cm .  La 
inductancia del canal de descarga se puede calcular con  la  siguiente expresi6n : 

59 
 L = 110
2rr 
[In  ­1] 
2d
r
x  d +  110 x  d
8rr 
(5 .2)  
_ H
donde  fl o = 47r x l0  7  _ 
m
Para  la  practica, 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  n 
Tmax  = 100  ｾｳ＠
Oeltat  = 0.1  ｾｳ＠

Calcular el  valor de  la  resistencia  del  canal , de  tal  manera  que  se  obtenga  10  mas 
proximo  a  una  forma  de  impulso  de  corriente  de  1.2/50  ｾｳ･ｧ Ｎ＠ Esta  resistencia , 
fisicamente,  es  de  un  valor  alto  y  es  la  que  asegura  que,  final mente,  un  rave  se 
pueda  representar circuitalmente como una fuente de corriente. 

Para  realizar  la  simulacion  de  la  corriente  de  rayo  se  usa  cualquiera  de  los  dos 
circuitos de la figura  25 . 

En  el  primer  circuito,  el  suiche  S1  cierra  en  tiempo  negativo  y  se  abre  en  t=0,
mientras  que  el  suiche  S2  se  debe  cerrar a  partir  de  t=0 y permanecer cerrado  el 
resto de tiempo de la  simulacion . 

En  el  segundo  circuito  se  utiliza  la  opcion  de  condensador  con  condicion  inicial. 
Este  metodo  tiene  restricciones  en  algunas  maquinas  cuando  hay  varios 
condensadores con  condicion  inicial. 

60  
 S 1  R canal  L canal  S2 

50MV 
"1"  T  ­I  Rprueba 

R canal  L canal  S2 

ｾ＠
ｾ＠ R prueba 

£ 
Figura 25.  Circuito para simular corriente de rayo 

Como  resultado  se  espera obtener una forma  de corriente como la  de la  figura  26. 

Corriente  [A] 
ＵＰ Ｇ ＢＭｯｾＬ＠

4000 

3000  .... ­ . 

.. - . 
Vt ' 
2000 

1000  •••••• ••••• ••• oj.

o I 
o  10 20 30 40 50
Tiempo [us] 

Figura 26. Forma esperada de  la  corriente de rayo 

Una  vez  que  se  obtenga  la  forma  esperada  de  la  corriente  de  rayo,  verificar  la 
sensibilidad  con  respecto  a  la  resistencia  de  prueba,  para  determinar  el  range  de 
valores de la  carga  en  el  que  el  rayo  realmente  se  pueda  simular como  una  fuente 
de corriente. 

61 
5.3  SISTEMA DE TRANSMISION A ESTUDIAR ! I "t­ ......,  1'1  11 1".

DEPTO .  D E  B I RLTOT E CJ\S 

131  BLT OTE CA  M I NAS 
Para  este  trabajo  se  considerara  una  linea  de  230  kV,  con  dos  cables  de  guarda 
como  se observa en  la figura  27 . 

8.0m 

6.0m 

6.0m 

20 .0m 

Figura 27.  Estructura tipica  de  230  kV 

DATOS DEL CONDUCTOR DE  FASE 
Rdc =  0.05086  Q/km
0  =  30.98  mm 

DATOS DEL CABLE DE  GUARDA 
Rdc =  1.9014  Q/km
0  =9 .78mm 

Donde 0  es  el  diametro del  conductor. 

62 
5.4 ASPECTOS DE MODELACION DE ESTRUCTURA Y CONDUCTORES

La  modelacion  de  la  estructura  metalica  para  un  estudio  de  descargas 
atmosfericas se  puede hacer como aparece en  la  figura  28. 

C  t e  R  It  Y  0 

RS T 
ABC 

Figura 28.  Trayectoria del  rayo y diagrama  para modelaci6n con  ATPDRAW­ATP 

En  la  modelacion  de  la  estructura  impactada  por el  rayo  se  consideran  los  tramos 
mas importantes.  Si  se desea,  se  puede detallar aun  mas , incluyendo  por ejemplo 
las  crucetas  en  las  fases,  pero  la  precision  que  se  logra  al  detallar  no  es 
importante. 

Cada  uno  de  los  elementos  metalicos  se  puede  modelar como  un  tramo  de  linea 
ideal  sin  perdidas ,  con  una  impedancia  caracteristica  dada  por  (estructuras  de 
forma  conica) : 

ZT =301n[ ＲＨｨｾ ｾ＠ l] 

r'  Q  (5 .3) 

Donde,  h es  la  altura  de  la  torre  y  r el  radio  en  la  base,  suponiendo  un  cono 
eq uivalente. 

63 
 Como valor tfpico para  Zr se  puede asumir un  valor de  200  n. 

La  velocidad de propagacion  es  aproximadamente la de  la  luz. 

La  torre  se  divide  de  acuerdo  con  los  tramos  que  se  necesiten ,  con  la  misma 
impedancia  caracterlstica  de  200  n  y  de  acuerdo  con  la  distancia  que  se 
considere.  EI  modelo  a  considerar  es  el  de  CLARKE  para  el  caso  monofasico  y 
sin  perdidas, como se observa en  la figura  29. 

­1 
001*­_ 
ir.ft:hLN.- 'I  
,  •
VI0J£  NOEIE 
FreM
PH..s E 

s_ 
l llC  To 
ｂ Ｌ｡ｲ､ｴｬｦ ｾ＠ • 
.iOC1QfJOC3
boo u...­o  '  I  Ｍ ｾ

11<*.... _ \!noel  '   l_ Iu.€

 
i
J _ 
s ......  j  ...... 
6.011 __ 
"""""' 
ｔ ｬ Ｂｾ＠ j - Ｂ Ｂ＾ｉｎｯ Ｚ ｾ＠ I.eIlOt  r
!!!OOELS 
ｔ ＬＮｾ＠

1'£5 
11_ 5" -
 
CO£rir.-.,-r I 

ｏ ｾ ｵｴ  

r CluTer'  r ｾ･ ｱ＠ (" C¥'.....,.oa r ｾ Ｈ ＡｗｉｴＹｙ

r.
r 
U _u.1I'Oi1
f_"""" 
ｾ＠ ［ ｾ＠ ee:­ ｾ＠  

Figura 29 .  Modelo de Clarke para tramo de red  monofasica 

Los  conductores  de  fase ,  conjuntamente  con  los  cables  de  guarda ,  se  modelan 
como  8  fases  acopladas  entre  sf,  mediante  el  modelo  de  K.C.LEE ,  a  una 
frecuencia  que  puede ser de  100  kHz  (/0  que  se  denom ina el  modelo  Bergeron  en 
el  ATPDRAW).  Suponer un  valor del  vane de 0.5 km 

Para  obtener  el  correspondiente  modele  se  puede  utilizar  la  opcion  LINE/CABLE 
de  la  caja  de  herramientas  del  ATPDRAW.  EI  modelo  a  utilizar  es  el  de  onda 
viajera  (Bergeron)  para  linea no transpuesta. 

64  
 " 
Ｇ ､ ｯｬ＠
ｾＢＬ
ｬｴ ｾ＠

M
..  · 1ftN!  t  I""  ', Ｇ｟ＮｪｾＬ＠
ｓｯ｡ｮｾｬｌ､ ｟＠

l'u 'll  !  \11'[1\ .,1,.01. ,) ... , '  .  I  d .  ... ",,·.H nI  •• ,. 

r  T_1IOUd 
r  ｾ ｩＮＢ ＢＬﾷ Ｎｧ＠  
Ml>tIeI  12"'8  I
>1 Sg  _  
Ｑ ｾＧ ［ＹＧ Ｐ＼ｊｮ［＠
,  I' ;.
ｲＮ ｾ＠
ICI" 
>1  SO.. Mml ".. ....  ("  Ａｮ ｾ I  ｾ､ｴ＠ ｾ

1549  0050086  H  32  ｾ ＮＳｊ ｾ＠ ＢＭ ｉｾ＠ ｾ Ｎ＠ ｾ＠

I<IodoI  1 549  005 0086  ­4 1  26
T)'p' 
ＱＮｳｾＹ＠ 0050086  <1  26 
"  IiIon;oftn  
r  Pi   1549  0050086  ­41  20 

ｾ＠
r.- 1.549  0 050086  41  20 
r  t!lod.  0 489  1 90 H  ­4  36 

I" ...  0489  ;  90 14  36 

ｾ Ｇ ｉ＠ X  ­:;u; Y  lUI 
ｾ ｾ｣｡ｬ＠ import 

eddraw  Qel.le 1801 row 

W<  !;encal  I  ｾ ｾ ａｳ＠ I  R"n eW  Iｾ ｾ t!e lp 

Figura 30.  Ventana de datos para obtener modelo  K.C .LEE  para la  linea 

La  cadena de aisladores  se  modela  como  un  suiche tipo "gap" que se  controla  por 

voltaje (figura 31) . 

Figura 31.  Modelo tipo "gap" para la  cadena de aisladores 

65 
 Para que se  pueda observar perfectamente el  voltaje que aparece sobre la  cadena 
de  aisladores,  se  coloca  un  voltaje  de  fiameo  bien  alto  con  respecto  al  que 
normal mente aparecerfa sobre la  misma  (5  MV,  por ejemplo) 

La  resistencia  de puesta  a tierra  se  modela como  una  del  tipo  concentrado,  con  el 
valor correspondiente  (20  n , en  forma  tipica) 

La  corriente  de  rayo  se  puede  simular con  una  fuente  rampa  de  doble  pendiente 
(identificada  como  fuente  tipo  13  en  el  ATP) ,  0  con  una  fuente  tipo  impulso,  bien 
sea  doble exponencial  0  de  un  exponencial  (identificada como fuente tipo  15  en  el 
ｾ＠
ATP) .  Se  recomienda  para  esta  practica  la  fuente  tipo  13  por  la  facilidad  de 
manejo. 

De  manera  adicional  se  puede  simular  el  comportamiento  del  aislamiento  de 
acuerdo  con  la  formulaci6n  que  propone  la  IEEE,  donde  la  resistencia  que  ofrece 
el  aislamiento depende del tiempo: 

Vjlameo = (0.4 + ｾ［Ｉ＠ x  d [MV] (5.4) 

Donde ,  t es el  tiempo en  J.1S 

d distancia longitudinal de la  cadena en  m 

Esta  respuesta  del  aislamiento se  ha simulado mediante la  opci6n de los TACS del 

ATP que permite implementar f6rmulas matematicas y bloques de control: 

Vjlameo =  (0.4 + ( t  ＰＧｾＩ＠10  ＰＷＵ＠ .  Jx # azsladores xlong por aislador(m) x 10

x 
6
[V] ｈ ｌｯＮ ｾＱｉ＠ >-
Para  controlar  el  hecho de  que para  valores  iniciales  muy  pequerios  de tiempo  el 
voltaje  de  flameo  tiende  a  infinito, se  utiliza  el  dispositivo  tipo  51  de  los  TACS ,  el 
cual  hace  que  la  serial  que  se  coloque  a  la  entrada ,  solamente  pasa  a  la  salida 
cuando se  de una  serial  de habilitaci6n . 

66 
 GAP15P 

Senal de  habil itaci6n 

Figura 32 .  Dispositivo 51  de  los TACS 

;: 

EI  funcionamiento del dispositivo 51  de los TACS  se  describe en  la  figura  33 

IN ·  •  ;:;'f  . . .  OUT 

IN •  .  )6.1  OUT  =:> 
D_sign  •  " 
/ ｾ＠
ｾ
T_oold  D_sign 
1tC  0.5  t  (,L S) 
T_oold  +  Lrold 

Figura 33 .  Descripci6n  del  dispositivo 51 

EI  suiche esta normal m ente abierto  y  cierra  si  se cumple la  siguiente condici6n : 

!D_sign!?: T _hold + I_h old (5 .5) 

D_sign es  una  serial  de entrada  que en  este  caso  es  una  del  tipo  escal6n  unitario 
retrasada  un  tiempo de 0.5 )..ls. 

T_hold es  una  serial  que  se  suma  con  el  valor  fijo ,  t_hold, para  determinar  la 
condici6n  de cierre del  suiche.  Para  el  caso  particular la  serial  T_hold no existe  y 
el  valor fijo  Chold es de 0.5.  La  condici6n  de cierre sera: 

!D sign! ?: 0.5 (5 .6) 

EI  sistema , finalmente , puede tener la  presentaci6n  de la  figura 34 

67 
 4.5 m  4.5 m 

QeRayo 

ｾ＠
ｾ • •  ｾ RST  ｶ＠
AOC 
ｾ ｉＢ＠ ­ ｾ＠ Q\P15P 

oc 

ｾ ｟Ｌｾ＠
Figura 34.  Sistema completo para  simulaci6n 

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  mas  influyen  son  las  que 
quedan  mas  cerca  de  la  torre  impactada  por  el  rayo .  EI  hecho  de  que  en  la 
realidad  sean  muchos  vanos  se  simula  con  una  linea  de  longitud  infinita  (para  que 
no  regresen  las  reflexiones) ,  0  mediante  elementos  resistivos  acoplados  con 
valores  de  resistencia  igual  a  las  impedancias  caracteristicas  de  la  linea .  Para  el 
caso  de  este  trabajo  se  pueden  considerar  tramos  de  linea  de  una  distancia 
suficientemente  grande,  para  que  las  reflexiones  regresen  fuera  del  tiempo  de 
anal isis,  10  cual  equivale  en  la  practica  a  considerar  infinita  la  linea .  Para  el 
sistema  se  han  considerado  tramos  de  linea  de  2.5  km  cuyas  reflexiones 
regresarian en  un tiempo aprox. de  16.7 !­Is. 

Hay que tener especial  cuidado con  el  delta  de tiempo  de  simulaci6n , ya  que debe 

ser mas  pequeno que el  tiempo de viaje  mas corto  (llt < 4. 5m/3x1 0
8 
m/s  = 15  ns) . 
EI  tiempo  de  simulaci6n,  Tmax,  debe  estar entre  10  Y 20  Jl.s  para  que  no  sea  muy 
pesada  la  simulaci6n  y,  adicionalmente,  en  este  tiempo  los  peores  sobrevoltajes 
sobre la  cadena  ya  se  han  presentado. 

68 
5.5  PROCEDIMIENTO Y PREGUNTAS

•   Interesa en  primera  instancia establecer la  diferencia entre  los  efectos sobre el 

aislamiento  (cadena  de  aisladores)  de  una  descarga  directa  en  la  fase  y  una 
descarga en  los cables de  guarda  0  en  la  torre. Para  este caso se  impactara la 
linea  primero  en  la  punta  de  la  torre  y  luego  en  la  fase .  Para  ambos  casos, 
determinar el  voltaje  que  aparece  en  la  cadena  de  aisladores  y  determinar  la 
influencia  del  valor  de  puesta  a  tierra  y  de  la  variaci6n  de  la  corriente .  GSe 
puede  controlar  el  sobrevoltaje  en  las  cadenas  de  aisladores  cuando  el  rayo 
impacta la  fase?  GPor que? 

•   Considerar la  repuesta  del  aislamiento en  el  tiempo y determinar que magnitud 

de  corriente  de  rayo  puede  romper  el  aislamiento  de  una  cadena  de  15 
aisladores.  En  el  ensayo  se  considera  la  simulaci6n  de  la  respuesta  del 
aislamiento en  el  tiempo  mediante los  TACS . 

•   Calcular  el  sobrevoltaje  que  aparece  sobre  cada  una  de  las  cadenas  de 
aisladores  para  una  corriente  de  rayo  de  30  kA,  2/50 )..!S  Y resistencia  de 
puesta  a  tierra  de  20  0  .  GEn  cual  cadena  aparece  mayor  sobrevoltaje  y  por 
que?  GQue  signo tienen esos voltajes,  cual  es  la  raz6n  de ese signo? 

•   EI  concepto de acople se  entiende como  el  voltaje que aparece inducido sobre 

las fases, debido  al  voltaje  generado en  los  cables  de  guarda  por  la  descarga 
atmosferica.  Normalmente,  se  define  en  forma  porcentual  con  respecto  al 
voltaje  que  aparecen  en  los  cables  de  guarda.  GPara  la  misma  situaci6n  del 
ftem  anterior, cual  es el  valor de estos acoples en  %? 

•   Identificada  la  cadena donde aparece el  mayor sobrevoltaje y partiendo de una 

condici6n  base  de  Tf de  1.2 )..!seg,  Rtierra  de  20  0  Y una  corriente  de  rayo  de 
30  kA, hacer sensibilidad  con  respecto  a  la  corriente  de  rayo  (1 kA a  200  kA) , 
resistencia  de  puesta  a  tierra  de  las  torres  (50 a  10000), frente  de  onda  (0.5 
)..!S  a 4 )..!s) .  Realizar un  ordenamiento de  estos  parametros  de  acuerdo con  su 
influencia  en  la  magnitud  del  sobrevoltaje  y  dar  las  razones  de  dicho 
ordenamiento. 

•   Para  una  corriente de  rayo fija  de  1 kA,  lIenar la  tabla  3 

69  
Tabla 3.  Archivo  plano generado por la  interfaz  Atpdraw 

Icritica (kA)  para 
Rtierra  (0)  Tfrente  Ｈｾｳ･ｧＩ＠ Vcadena  (kV/kA)  diferentes cadenas 
12  13  14  15  16  17 
20   1.2 
5   2.0 
10   2.0 
20   2.0 
50   2.0 
100   2.0 
200   2.0 
1000   2.0 
20   4.0 

La  corriente  critica  se  obtiene dividendo  la  tensi6n  nominal  soportada a/  impulso 
atmosferico  (BIL) ,  total  de  la  cadena ,  para  impulsos  ｮ･ｧ｡ｴｩｶｾｳＬ＠ por  el 
sobrevoltaje  que  aparece  sobre  la  misma  con  un  kA  de  corriente  de  rayo . 
Considere  que  el  BI L  de  cada  aislador  es  de  120  kV  en  seco  y  80  kV  en 
humedo. 

Por ejemplo, si  para  un  Tf de  1.2  ｾｳ＠ y una  Rtierra de  20  0  el  sobrevoltaje sobre 
la  cadena  fue  de  20.3  kV/kA  (corriente  de  rayo  unitaria  de  1kA),  la  corriente 
critica  de  flameD  para  una  cadena  de  15  aisladores seria  de  15 x  120/20.3 kA, es 
decir  88.7  kA. 

•   Para  cada  una  de  las  corrientes  crfticas  de  flameo , calcular la  probabilidad  de 
que  en  realidad  estas corrientes sean  superadas .  Existen diferentes funciones 
de  probabilidad  para  las  corrientes  de  rayo.  Una  de  uso  frecuente  es  la 
propuesta por Anderson­Eriksson: 

70 
 P(J"210)=
1+   ­
1J2.  (5 .7) 
( 31 

La  funci6n  de probabilidad  anterior dice que posibilidad tiene una corriente crftica 

de  s er  superada.  Para  el  caso  del  ejemplo,  la  probabilidad  de  que  la  corriente 
sea superada serra de 6.1 % 

•   Tal  como  se  mencion6,  una  manera  de evaluar el  desempeno  del  aislamiento 
frente  a  un  esfuerzo  originado  por  descargas  atmosfericas  es  mediante  el 
fndice :  Numero de salidas/100 km/ano. Para  calcular  el  numero  de  salidas  se 
debe  conocer  el  numero  de  descargas  que  caen  por cada  km  de  longitud  de 
linea .  Una  vez  conocido  este  numero  de  descargas  se  multiplica  por  la 
probabilidad  de  que  superen  la  corriente  critica.  Adicionalmente,  hay  que 
considerar  que  no  todas  las  descargas  caen  sobre  la  torre,  sino  que  un 
porcentaje  cae  a  mitad  del  vano,  donde  no  se  produciria  f1ameo  porque  la 
distancia en  aire entre  conductores es  mayor que la  longitud  en  aislamiento de 
cadena  y  por  el  fen6meno  de  predescarga.  Se  considera  que  de  todas  las 
descargas  que  caen,  el  60%  10  hace  sobre  la  estructura  0  dentro  del  area  de 
influencia  de  la  misma  y  que  el  40%  cae  a  mitad  del  vano .  Para  calcular  el 
numero de descargas  que caen  por km  de linea  y  por ano se  puede  utilizar la 
siguiente relaci6n  [10] : 

Numera de descargas / km / ana = (4h + b )Ng (5.8) 

Donde,  h =  altura del cable de guarda (km) 

b =separaci6n  horizontal entre cables de guarda (km) 
Ng = Densidad de rayos en  la  zona (No descargas/km 2/ano) 
Se puede suponer un  valor tlpico para  Ng  de  10  para nuestra regi6n. 

Para  las  dimensiones  de  la  torre,  el  numero  de  descargas  serra  de  1.52 
descargas/km/ano.  Para  100  km  serfa  de  152  descargas/100  km/ano.  Y  el 
numero  de  salidas,  si  se  considera  una  cadena  de  15  aisladores,  serfa  de 
152xO.061 xO.6, es decir, 5.56  salidas/100 km/ano. 

71 
•  Si  se  quiere  manejar el  criterio  de que  el  numero  de  salidas  por  cada  100  km 
por  ario  este  en  el  rango  de  1  a  3,  cual  serra  el  valor  de  la  resistencia  de 
puesta a tierra  y el  numero de aisladores a garantizar? 
•  A  partir  de  una  condici6n  base  como  punto  de  partida  (Rt=200) ,  deducir  un 
criterio  basico  de  diserio que diga:  cuantos  ohmios  hay  que  rebajar a partir de 
un  valor de  referencia  de  20  n  de  puesta  a tierra , para  que sea  equivalente a 
colocar un aislador mas a partir de  15 aisladores. 

72  
 6.  PRAcTICA N°4.  SOBREVOL TAJES TRANSITORIOS EN SISTEMAS DE
DISTRIBUCION, DEBIDO A MANIOBRA DE CAPACITORES

En  los  sistemas  de  potencia,  una  alternativa  para  realizar  control  del  voltaje, 
compensar  reactivos  y  dar  soporte  de  voltaje ,  es  la  instalaci6n  de  bancos  de 
condensadores . Esta alternativa es  de  bajo costo , comparada con  otras  soluciones 
como  generadores  sincr6nicos .  La  conexi6n  se  realiza  en  forma  de  derivaci6n 
(shunt)  y se  debe conmutar de acuerdo con  la  potencia reactiva requerida. 

La  conexi6n  y  desconexi6n  de  los  bancos  capacitivos  generan  transitorios,  que 
afectan  el  comportamiento de  la  red  al  originar variaciones  del  voltaje  que  pueden 
ocasionar la desconexi6n de cargas electr6nicas sensitivas  0  afectar el  aislamiento 
de  equipos .  Los  equipos  que  mas  sufren  por estas  variaciones  de  voltaje  son  los 
de tipo electr6nico,  como variadores de velocidad  (ASD) . 

6.1  ENERGIZACION DE UN BANCO DE CAPACITORES

Cuando  se  energiza  un  banco  de  capacitores  se  presenta  una  oscilaci6n  del 
voltaje  en  mayor  0  menor grado, dependiendo  del  valor instantaneo  que tenga  la 
misma  fuente al  momenta de  la  energizaci6n . 

La  oscilaci6n  se  presenta  debido  a  que  un  capacitor  no  puede  cambiar 
instantaneamente  el  valor  de  su  voltaje .  En  la  barra  donde  esta  instalado  el 
banco,  el  voltaje  instantaneamente  forma  un  valor  de  0  voltios ,  sigue  un  voltaje 
n3pido  de  recuperaci6n  con  un  pico  de  voltaje  que  puede  alcanzar,  en  forma 
te6rica ,  un  valor  de  2.0  por  unidad .  Para  sistemas  reales ,  debido  al 
amortiguamiento que se  presenta , los  sobrevoltajes  pueden  alcanzar valores entre 
1.1y1 .6pu . 

73  
 EI capacitor es suicheado 
aqui 

ｾＨ＠
ｾ Ｍｪ
)>­­­­­­­t­I ｾ Ｍｴ Ｍｴｲｩ＠

T 
Transitono abservado en algun punta de 
carga 

Figura 35. Diagrama unifilar de un  circuito de distribuci6n con  capacitor para suicheo 

En  la  figura  35  se  observa  el  diagrama  unifilar de  un  circuito de  distribuci6n , donde 
esta  instalado  un  banco  de  condensadores  y  se  realiza  el  correspondiente 
suicheo.  Se  desea  observar  los  voltajes  durante  el  suicheo  en  el  banco  y  en  un 
punto de  la  red  primaria, entre el  banco  y la fuente de  potencia . 

Para  la  energizaci6n  del  banco  ｣｡ｰｩｴｶｾ Ｌ＠ el  circuito  equivalente  es  el  que 
aparece  en  la  figura  36 .  La  inductancia  L1 representa  la  impedancia  de  la  fuente 
mas  la  del  transformador  de  potencia  de  la  subestaci6n .  La  resistencia  ｒｾ＠ y  la 
inductancia  L2 representa  la  impedancia  de  la  red  primaria .  C  corresponde  a  la 
capacitancia del  banco. 

Ll R L 

'U/ 
, 

I e 
­

Figura 36 .  Circuito equivalente para energizaci6n del  banco capacitivo sin  carga en  red 

primaria 

74 
La  frecuencia  de esta oscilaci6n  esta  dada  por: 

1 
Is =  ＲＱ｛ ｾ ｌＬ＠ x  C 
(6 .1 )

EI  valor pico de  corriente esta dado por la expresi6n : 

I  ­ VpK  (6.2) 
PK­ ｾ Ｎ  ｳ＠ ＧＮ
＠ｾ ｾ ｃ＠  
Donde:  Ls es  la  inductancia  equivalente  (L 1+L 2 ) desde  el  banco  hasta  la  fuente 
(H)  
C es la  capacitancia del  banco (F)  
VPK es el  valor pico del  voltaje fase­tierra  de  la fuente(V)  

Si  la  energizaci6n  se  realiza  cuando  el  voltaje  esta  en  el  valor maximo,  las formas 
de los voltajes en  el  banco  y en  un  punto de la  red  primaria , situado entre  el  banco 
y la  fuente  de potencia , se observan en  las figuras  37  y 38  respectivamente . 
Voltaje  [pu] 
2.0  I 

1.5 

1.0 

0.5 

0.0 

­0 .5 

­ 1.0  I 
o 
"­'  "4 '  . '­' 
I
10  20  30  40  50 
Tiempo [ms] 

Figura 37.  Voltaje observado en  la  barra donde esta  ubicado el  banco 

75 
 Vollaje [pu] 
1.5  I 

1.0 

0.5-

0.0 

-0.5 .; 

­1.0  I  ><  ¥ '-£

o  10  20 30 40 50
Tiempo  [ms] 

Figura 38. Voltaje observado en  un  punto del  circuito  primario entre el  banco  y la fuente 

de potencia 

La  frecuencia  de  oscilaci6n  es  del  orden  de  (300­1000)  Hz.  La  magnitud  de  la 
oscilaci6n  depende  de  varios  factores  relacionados  con  las  caracteristicas  de  la 
fuente  y  de  la  red  de distribuci6n, 10 mismo que con  la  capacidad  del  banco. 

6.2 MAGNIFICACION DE VOL TAJE POR SUICHEO DE BANCOS DE

CONDENSADORES

EI  fen6meno  de  magnificaci6n  de  voltaje  se  presenta  cuando  se  energiza  un 
banco  ｣｡ｰｩｴｶｾ＠ en  un  circuito de distribuci6n  y la  frecuencia  natural  de  oscilaci6n 
excita  un  banco situado  en  el  lade de  baja  de  un  usuario  (figura  39).  La  magnitud 
de  esta  oscilaci6n  es  mucho  mas  severa  que  una  energizacion  normal  de  un 
banco  ｣｡ｰｩｴｶｾＮ＠ EI  pico  de  la  oscilaci6n  puede IIegar a valores  de 4.0 en  pu.  Para 
que ocurra  la  magnificaci6n  del  voltaje se  debe cumplir: 

1;  ': :', /2

1 1
con  .f..  = 2 
7r ｾ＠
IL­'I  xC '  I  '  Y 12 = 27r "'Vｾ  D2  X C 
2

76
 Ll L2
TrBnsformador 
de  usuario 

ｾ＠
Allm,ofaOO, d,  Ｎｾ ｾ＠ • C".,. 

Subestaci6n 
Banco  de  Banco de 

T
condensadores 

e 
, 
c' T 
condensadores del 
usuario en  el  lado de 
b.ja "  ..11)0

Figura 39 .  Diagrama unifilar circuito de distribuci6n  con  capacitor para suicheo 

Ll  L2

"v 
Tel  T C2 

Figura 40.  Circuito equivalente para estudio de  magnificaci6n de voltaje 

La  frecuencia  natural de oscilaci6n  del  banco en  alta se  puede calcular can : 

MVA corto circuito lado de alta x;;

II =
MVAR banco capacitivo de alIa a 
(6 .3)

y la  del  banco de baja ten si6n  se  calcula como: 

77 
 f ==  K VA tramformador x  j,
2  ｾ X CC( pU) transformador x KV AR  banco capacitivo de baJa 0 
(6.4) 

En  la  anterior  aproximaci6n  se  esta  suponiendo  que  el  transformador  de 
distribuci6n  se encuentra muy cerca del banco situado en  la red  de distribuci6n .  Si 
la  longitud  de  la  red  primaria  es  apreciable,  se  debe tener en  cuenta  en  el  calculo 
de la  inductancia  L2 . Para este  caso , la frecuencia  natural  de oscilaci6n  del  banco 
de baja  sera  menor. 

ｾ］ ｘｨ＠ (6.5) 
f2 =  l(lOOO XXred (O ) + Xcc(p u) transfo rmador Ix 
KVAR banco
Vlinea\ k V ) K VA transformador )

Donde  fo  es  la  frecuencia  fundamental  y  Xcc(pu)  es  la  reactancia  de cortocircuito 
del  transformador en  por unidad . 

Una  expresi6n  aproximada,  que  permite  verificar  si  existe  posibilidad  de 
magnificaci6n  de  voltaje  en  el  lado  de  baja  por  suicheo  del  banco  de  alta ,  es  la 
siguiente: 

KVAR banco de baja MVAR banco capacitivo de alta x  (6 .6) 

KVA transformador MVA corto circuito /ado de alta Xcc(pu ) transformador

La  expresi6n anterior es  muy practica  porque  la  relaci6n  entre los  KVAR del  banco 

y  los  KVA  del  transformador  debe  ser  menor  que  la  unidad  y  rapidamente  se 
verificarfa,  para  todo  un  circuito  de  distribuci6n,  si  existe  posibi lidad  real  de 
magnificaci6n  para  un  determinado  banco  que  se  haya  instalado  en  la  red  de 
distribuci6n. 

6.3  SISTEMA A ESTUDIAR

EI  sistema  a  estudiar es  un  circuito  tfpico  de  distribuci6n , de  13.2  kV,  compuesto 
de  cable  aislado  y  red  aerea.  Un  banco  ｣｡ｰｩｴｶｾ＠ se  energiza  en  el  lado de alta, 
cerca  de la  subestaci6n  de potencia  de distribuci6n .  Tambien  se  coloca  un  banco 
de  correcci6n  de  factor de potencia en  el  lado  de baja  de  un  usuario y  se  analizan 
las  condiciones  bajo  las  cuales se  puede presentar magnificaci6n de voltaje. 

78 
6.3.1 Equivalente de cortocircuito en la barra de 13.2 kV.

Como fuente  se  tomara el  equivalente de  cortocircuito en  la  barra de  13.2 kV 

Voltaje de  prefalla  = 13.2 kV 

ｉ｣Ｓ ｾ＠ = 13774 A ,  XlR  = 15 

Icc1 ｾ＠ = 17907 A ,  ｘｬｒｾ ＼ｸＺ＠

Donde  ｉ｣Ｓ ｾ＠ es  la  corriente de cortocircuito trifasica e  Icc1 ｾ＠ es  la  corriente de 

cortocircuito  monofasica. 

6.3.2 Cable aislado de 350 KCM para la salida del circuito de 13.2 kV

Z1  =0.1455 + jO .1413 nlkm 

Zo  =  0.4795 + jO .1295 nlkm 
C1  = Co  = 170 nF/km 
En  forma tipica se  considerara  una longitud de  1 km  de cable aislado. 

6.3.3 Red aerea de 13.2 kV

14 m
­. 114 m  -

065m

11  n1 

', '

Ｇ ＱＭ ｾ［Ｌ ｴ ＧＱ ＭＡｪ Ｉ＠ ｾＧＺ Ｍ［ ＭＢ ［Ｇｾ＠ ｲＭＺＮｔｾＧ ﾷ＠ ＬｆＮｖｾ ＢｖＯ ｟ ｾｪ

Figura 41. Configuraci6n de estructura de 13.2 kV , neutro superior

79
Conductor de fase

Tipo =  Penguin  (4/0 ACSR) 

Rdc (50  DC)  = 0.2979  Q/km
Diametro =  14.31  mm 

Conductor neutro

Tipo  = Raven  (1/0 ACSR)  superior 

Rdc  (50  DC)  =0.5091  nlkm 
Diametro =  10.11  mm 

Resistividad  del terreno = 250 Q*m 

Se  considerara  una  longitud  de  red  aerea  de  5  km , con  cargas  concentradas cada 
km  de 500  KVA, fp  de 0.9 

6.4 SIMULACIONES

•   Energizacion  de  un  banco  capacitivo  de  500  KVAR  instalado  al  final  de  la  red 
aerea . 

RAEREA  RAEREA  RAEREA  RAEREA  RAEREA 

SOO  tWV  FP=0.9  ­=

Figura 42. Esquema circuital en  Atpdraw para  energizacion de un  banco capacitivo 

80  
De  la  energizacion  con  carga  conectada  se  esperan  resultados  como  el  de  la 
figura 43. 

······· 
VoltaJes  [kV] 

20 

E 
­- - -- ..  .-- - ｾ＠ . - ..  ...
15  •• - 10'r - '" .\ ,. •• ;,. , -•
. . . __ ｾＬ＠ rOO ., Ｎ＠ Ｌ＠ ｾ＠ Ｎ｟＠ ｾ＠ • • _ \ _

­- . ｾ＠ - .:. - - : - . -' - - . - _. r
, . ,. .
10  ::....'"!-.:._..._ - - . . ..ｾ＠

-5

-10

· 15
6 10 20 30 40 50
Tiempo [ms]

Figura 43. Resultados esperados durante la energizaci6n

• Determinar el sobrevoltaje maximo en pu que puede aparecer. i,Como varia el

sobrevoltaje con respecto a la carga que se encuentre conectada a 10 largo del
circuito de 13.2 kV? (carga de 0-5 MW).

• i,De que depende la frecuencia de oscilacion y como varia con respecto a la

carga conectada?

• i,Como varia el sobrevoltaje a 10 largo del circuito?

• i,Si se ubica el banco a la salida del cable aislado , como varian los
sobrevoltajes en magnitud y frecuencia de oscilacion?

81
•   Con  el  fin  de  observar el  fenomeno  de  magnificacion  de  voltaje  se  simulara  el 
sistema de la figura  44 . 

RAEREA  RAE REA  RAE REA  RAEREA  RAEREA 

873  KVAR  $  500 KW 

FP=0.9 
500  KW 
FP=0.9 
500 KW 
FP=0.9 
500  KW 
FP=O.9 
60  !WAR

+ 

Figura 44.  Esquema circuital  en  Atpdraw para estudio de magnificacion de voftaje 

•   AI  final  de  la  red  primaria  se  conecta  un  banco  ｣｡ｰｩｴｶｾ＠ de  50  KVAR  por 
medio  de  un  transformador  trifasico  de  500  KVA,  X=5% ,  13200/207 V, 
conexion  L\- Y  aterrada .  AI  final  del  primer  km  de  red  aerea  hay  un  banco 
｣｡ｰｩｴｶｾ＠ de  873  KVAR ,  el  cual  se  energizara  cuando  el  banco  de  baja  se 
encuentra en  operacion normal. 

•   Determinar  los  sobrevoltajes  que  aparecen  en  los  bancos  capacitivos  para 
diferentes  condiciones  de  carga  de  la  red  primaria.  GCual  banco  presenta 
mayor  factor  de  sobrevoltaje?  GCual  es  la  razon  de  la  diferencia?  GComo 
varian  los  sobrevoltajes  con  respecto  al  aumento  de  longitud  de cable  en  baja 
tension ,  tomando  como  referencia  la  condicion  resonante?  GCual  sera  el 
efecto de aumentar la  longitud de red  de  13.2 KV?  Explicar las  diferencias. 

•   Para  la  red  primaria  sin  carga  y un  banco fijo  de  50  KVAR  en  el  lade  de baja, 
determinar los  KVAR  capacitivos en  el  lado  de  13.2 kV que crean  condiciones 
de  magnificacion  de  voltaje  para  diferentes  puntos  de  ubicacion  a  10  largo  del 
circuito. 

82  
 7.  pRACTICA N°  5.  SOBREVOLTAJES A BAJA FRECUENCIA DEBIDOS A
FALLAS ASIMETRICAS

Cuando  hay  una  falla  en  un  sistema  electrico  de  potencia  se  presenta  una 
condici6n  transitoria  que  se  amortigua  rapidamente ,  quedando  una  condici6n  de 
estado  estacionaria,  pero  donde  las  variables  electricas  no  presentan  el  mismo 
valor  que  antes  de  la  falla.  Esta  condici6n  se  mantiene  hasta  que  la  falla  sea 
despejada. 

Para  fallas  de tipo asimetrico  10  que ocurre  es  que las  fases  no  falladas  presentan 
un  sobrevoltaje , dependiendo de  las relaciones de  XoIX1 Y RoIX1.

La  magnitud  de  estas  sobretensiones  temporales  depende  basicamente  del  grado 
de  aterramiento del sistema en  el  punto de  la falla. 

7.1  SISTEMA A ESTUDIAR

Para  esta  practica  se  analizara  un  sistema  tfpico  de  distribuci6n  en  una 
subestaci6n , con  niveles  de  voltaje  de  110,  44  Y 13.2  kV,  como  se  observa  en  la 
figura  45. 

Este  sistema  esta  constituido  por  una  red  de  44  kV  Y una  red  de  13.2  kV, como 
ilustra el  siguiente diagrama unifilar 

A 
ｾ＠

ｾ＠ 10 KV 

>  I 13.2 KV 

>.. 
iJ.  "n­

44  KV 

Figura 45.  Oiagrama unifilar de  un  sistema de  distribuci6n 

83 
7.2 DATOS DE LOS ELEMENTOS DE RED

Como fuente se tomara la  barra de  110 kV 

Equivalente de cortocircuito en la barra de 110 kV

Voltaje de prefalla =  112.08  kV 

Voltaje base =  110 kV 
Potencia  base  = 100 MVA 
Z1  =0 .0143 + jO.0583  pu 
Zo  =0.0163 + jO .0856 pu 

Transformador de potencia

60120/60 MVA  
110/46.9/15  kV  
ZAB = 17.85% (Base 60 MVA)  
ZAM =9.87%  (Base 20  MVA)  
2MB = 2.86%  (Base 20  MVA)  

Cable aislado de 350 KCM para salida de circuito de 13.2 kV

Z1  = 0.1455 + jO .1413 

Zo  =0.4795 + jO .1295 

Cable aislado de 500 KCM para salida de circuito de 44 kV (pantalla

electrostatica)

Z1  =  0.104 + jO .248 

Zo  =  1.074 + jO .845 

Red aerea de 13.2 kV

Conductor de fase 

Tipo""  Penguin (4/0 ACS R) 

Rdc (50  °C)  =  0.2979 Q/km 
Rac (50  °C)  =  0.3616 

84 
Diametro;;;;;  14.31  mm 

Conductor neutro

Tipo;;;;;  Raven  (I/O  ACSR) superior 

ｾ＠
Rde (50  °C)  =  0.5091  O/km 
ｾ ｾ ［ｾＭｪｊ
ｾ .­¥­.. ｾ＠

ｕ Ｂ ｉ ｖｆ ｒ ｾｉＩ ａｏ＠ ｎ ａｃ ｊＨ Ｉｎ ｾ＠ ｛ ｟＠ DE'  C n l OMllIA 

Rae (50  °C)  = 0.6031  ntkm  "" I", '0' .  I "'" 

Diametro ;;;;;  10.11  mm  DEPTO .  DE  RTBLTOT E CAS 

RIRLIOT ECA  MI NAS 

Resistividad  del terreno = 250 Om 

0.8m 

Ｍ ﾷｾ＠ ［ Ｑ ＱＮＴ＠ m 

0.65  m 

:  #  . 

11 m

"  ,',' 

Figura 46. Configuraci6n de estructura de  13.2 kV , neutro inferior 

85 
Red aerea de 44 kV

Conductor de fase 

Tipo =  Waxwing  (266.8 ACSR) 

Rdc  (50  °C)  =  0.2383  n/km 
Rac (50  °C)  =  0.2432  nlkm 
Diametro  = 15.47  mm 

Conductor neutro 

Tipo =  Raven  (1/0 ACSR) superior 

Rdc (50  °C)  =0.5091  n/km 
Rac (50  °C)  =0.6031  nlkm 
D1ametro  = 10.11  mm 
Resistividad  del  terreno  =  250  n*m 

1.5 m 

0.85  m 

.:1" 
ｾＮ ＧＭｾ＠

r •  • \"

"  m  ｾ＠ ｾ｜＠
.•. .. ｾ ｾ

Ｌｾ .. 
'­'­r.x;:""0V  Ｏ＠ Ｂ＠ ｖ ｾＯＮＧｖ ﾷ ﾷｖｻＬＭＺ［ ｾ ｖ＠ · Z ' ,\.Iy.v ''(  ｊ ｾ ﾷ ｖ＠

Figura 47. Configuraci6n  red  de 44  kV 

86 
Las  saltdas en  cable son de  1000 m, tanto para  13.2 como para  44 kV. 

EI  sistema  de  44  es  en  delta  aislado;  eI  de  13.2  kV,  y  eI  de  110  kV  estan  en 
estrella s61idamente aterrados en  la  subestaci6n  yen cada  estructura . 

Los circuitos tienen cada  uno una distancia de  10  km . 

7.3  CALCULAR Y SIMULAR

•   Si  se  presenta  una  falla  monofasica  a  tierra  al  final  de  la  red  de  13.2  kV, 
calcular  los  sobrevoltajes  que  se  presentarian  en  las  fases  sanas  en  el  punta 
de falla.  Calcular los voltajes que aparecerian en  la  red  de 44  kV. 

•   Si  se  presenta  una  falla  monofasica  en  la  red  de  44  kV,  calcular  los 
sobrevoltajes que aparecerian  en  la  red  de 44 kV Y en  la  de  13.2  kV. 

•   Cual  serra  el  efecto en  los dos anteriores casos si  la  red  de 44  kV se aterra por 
medio  de  un  transformador  zigzag  de  864  kVA,  2.3%  de  impedancia  de 
secuencia cero y 57  n  de resistencia  de puesta a tierra? 

•   De  acuerdo  con  el  comportamiento  observado  del  sistema  de  44  kV  y de 
13.2  kV, dimensionar adecuadamente  los  pararrayos  desde  el  punto  de  vista 
de  operaci6n  en  estado  estacionario.  Determinar  el  voltaje  nominal  del 
pararrayos y el  MCOV en  ambos niveles de tensi6n . 

87  
 8.  pRACTICA N°S.  CARGAS NO LINEALES - PUENTE RECTIFICADOR DE
SEIS PULSOS

8.1  INTRODUCCION

EI  servicio  de  energia  electrica  se  debe  prestar  en  unas  condiciones  de  calid ad 
para  la  potencia  adecuada.  La  manera  basica  de  medir esta  calidad  es  mediante 
la  forma  de  onda  del  voltaje,  que  debe  estar  exenta  de  distorsiones  en  forma 
permanente  y  en  estado  transitorio .  En  la  practica ,  la  completa  forma  sinusoidal 
del  voltaje  no  se  puede  conseguir,  pero  la  distorsion  que  tenga  la  forma  de  onda 
del  voltaje  debe  estar  ajustada  a  unos  topes  que  se  definen  en  cada  pais  de 
acuerdo con  una  normativa. 

La  distorsion  del  voltaje  ocurre  debido  a  la  presencia  de  cargas  de  tipo  no  lineal, 
que a pesar de  suministrarles  un  voltaje sinusoidal,  la  corriente que  requieren  para 
su  funcionamiento  es  no  sinusoidal.  Esta  corriente,  con  componentes  armonicas, 
al  circular  por  elementos  del  sistema  de  potencia  que  ofrecen  oposicion  a  la 
corriente  (impedancias),  generan caidas de tension  igualmente distorsionadas que 
afectan  la  forma  de  onda  del  voltaje, en  mayor grado  mientras  mas  cerca  se  este 
de  1a  carga inyectora de armonicos de corriente. 

Generalmente,  los  armonicos  de  tipo  impar  son  los  de  interes  en  un  sistema 
electrico de  potencia .  Esto es  propio de  equipos que desde el  punto  de  vista  de la 
corriente  alterna  son  bilaterales,  es  decir,  no  importa  la  polaridad.  Este  tipo  de 
equipos solamente producen  armonicos impares. 

Existen  algunos equipos que si producen armonicos pares : 

•  Rectificador de media onda 

•  EI  transformador, al  producir la  corriente inrush  durante  la  energizacion 

•  Hornos de arco, durante periodos de arco inestable 

88 
La  presencia  de  un  arm6nico  par  generalmente  indica  la  falla  de  un  equipo  0  la 
presencia  de  una  componente  DC  que  puede  contribuir  nocivamente  a  la 
saturaci6n de los transformadores. 

Las fuentes de armonicos, en  forma general , son tres: 

•   Dispositivos ferromagneticos.  Transformadores  y  motof6s 

•   Dispositivos  de  arco  electrico .  Lamparas  de  arco:  fluorescente,  mercurio, 
sodio­vapor, etc.).  Hornos de arco. 

•   Convertidores  de  electronica  de  potencia.  Drivers  de  velocidad  variable  para 
motores  (Adjustable  Speed  Drivers  0  ASD)  Fuentes  ininterrumpibles  de 
potencia  (UPS),  rectificadores,  fuentes  de  potencia de computadores. 

Estos  elementos  se  encuentran  conectados  generalmente  en  forma  shunt  en  el 
sistema  de  potencia  y,  mas  exactamente,  en  el  sistema  de  consumo,  es  decir, 
directamente  en  la  carga,  en  los  puntos  de  baja  tension .  Los  elementos serie  del 
sistema  de  distribuci6n  son  en  su  mayorfa  !ineales.  Los  transformadores  son 
dispositivos  que  tienen  dos  componentes:  una  parte  lineal  en  serie,  que  es  la 
impedancia de dispersi6n  (leakage),  y  una  parte shunt,  que es  la  magnetizaci6n, la 
cuat es del tipo no linea! per sus caracterfsticas de saturaci6n .

8.2  DESCRiPCiON DEL SiSTEMA A ANALiZAR

Ei  sistema a considerar en  esie trabajo es  iiustrado en  la  figura  48 

Rectificador de 6
Fuente 

ｓ･｣＠
=Vn  44 kV 
ﾷ ｾ＠ 266.7 MVA 
(Icc  30 =3.5 kA) 
6  '{ 
)( 
[  let1 pulsos
3.5MW
FP  =.  0.88 

ｾ Ｍ ｈｾ ｾ＠
) 
Punta de  comun  Carga lineal 
acop!e  (peC) 
44000/480 V 
5MVA 
X  =7 .6  ｾ Ｏｏ＠
L  2.0MW
2 .23 MVAR

Figura 48.  Sistema industriai a anaiizar 

89 
Ei  sistema es el  de una planta industrial que se al imenta de un  circuito de 44  kVen 
delta  aislada.  EI  punto  de  comun  acople  (PCC),  en  este  caso,  serfa  la  barra  de 
44  kV en  la  subestaci6n  de  la  planta .  La  planta  basica  la  forman  un  transformador 
de  potencia  de  5  MVA 44000/480  V,  una  carga  e!ectr6nica  tipo  puente  rectificador 
de  6  pulsos  de  3.5  MW,  FP  de  0.88  y  el  resto  de  la  planta  que  puede 
representarse  por una carga  lineal  de 2.0 MW y 2.23  MVAR. 

Durante  condiciones  de  demanda  maxima,  el  iacror  de  potencia  de  la  pianta 
alcanza  el  valor  de  0.8  y  se  desea  compensar  este  valor  a  0.93 ,  mediante  la 
instalaci6n de  un  banco de condensadores al  100 % de operaci6n. 

8.3  OBJETIVOS

Famiiiarizarse  can  ios  modeios  impiementados  y  verificar  ios  caicuios  reaiizados 

para el  sistema base (equivalente del  sistema,  transformador y carga  lineal) . 

Con  el  sistema  sin  compensacion  reactiva  (sistema  base)  observar las  formas  de 
onda de la  tensi6n  y la  corriente en  44  kV y 480V. 

Ca!cular  los  reactivos  necesarios  para  obtener  un  FP  =  0.93 .  Modelar  el  sistema 
compensado y observar las formas de onda de la  tensi6n  y la  corriente . 

Calcuiar  ei  contenido  arm6nico  y  THD  de  ias  ondas  de  tensi6n  y  corriente 
mediante  la  subrutina  Fourier  0  mediante  el  uso  de  los  MODELS  (para  el  caso 
actual  se  suministran  dos  modelos para  caleulo el  de arm6nieos). 

Verifiear  si  existen  proolemas  de  resonancia  arm6niea  y  determinaria  can  ei 
caleulo  de  la  impedaneia  arm6niea  en  barras  de  480  V,  mediante  la  subrutina 
Frequency Scan. 

Disenar  el  banco  de  eondansadores  como  fiitro  de  absorci6n ,  para  disminuir  ei 
impaeto de  la  carga  no lineal  dentro de la  planta  y haeia el  sistema . 

Con  e!  sistema  compensado  y  sintonizado,  observar  las  formas  de  onda  de  la 
tensi6n  y la  corriente y caleular su  contenido arm6nico y distorsi6n total  THO. 

90 
En  cada  uno  de  ios  items  se  debe  verificar  en  ei  punto  de  comun  acople  con  ia 
planta  (nivel  de 44  kV) . Los  indices de distorsi6n,  con  los limites establecidos en  la 
norma IEEE  519  de  1992. 

Debe  verificarse,  mediante  caicuios  manuaies  (similar  a  io  hecho  en  ciase) ,  el 
efecto  de  la  resonancia  generada  por  los  condensadores  y  el  impacto  del  filtro 
disenado a la  frecuencia arm6nica del  problema . 

8.4  PUENTE  RECTIFiCADOR DE  SEiS PULSOS 

La  carga  no  lineai  eiectr6nica Hpica  es  ei  puente rectificador de  seis pulsos , como 
se  observa en  la  figura  49. 

SJs!emaAC 
ｮ ｾ＠

Figura 49.  Puente rectificador de seis puisos 

Ei  eiemento basico del  puente es ei  diodo,  que dentro de ia  librer[a del  programa 

ATP es el  suiche tipo  11. 

Erobe s  &  3­phase 

Branch Lin.. ",  
Branch tio n lin e er  
L ine L umpe d  
tJ n @Oistrib uted   ｉ ｾＮ＠ ｾ ­­­=­r­ rN Oce  JJ!ltiASE  ,..u ·-
I 

,Sourc e !;  
S'Ni'tch ti me ｾｯ ｮ ｴｲ Ｎ＠
S witch lime J­ph  ==­__  ­"i (J 
o  ""T 

M A.chin Qs  ｾ ｣ｨ＠ voltage c ontY.  

T{MSformers 
ｾ｡ｬｶ･＠ (lype 11)
MODELS 
T rfac (typ e  1 2) 
｟ＢＧ ｾ＠ '­­1>.. 1 
ｔｹｰ ＮＡｉ ｾ＠
ｃｯｾ ｲ Ｌ ｬ［Ｇｲ＠
T AGS  switch (typ e  1 3) 
r­ H it.
I ACS 
M e asurin g 
""'""" 
r ｃＧｊｴｲｉ ｾ＠ ("" ＧｌｯＡｴ ｾ＠ ,..  ｣ＮｷＧ ｉｬｾＧ＠ r E:c­"'f.En. rw  r 
Use, S pecified 
ｌ ｭ･Ｏｾｬ｡＠ SlAlistic swi1ch 
OVerhe ad Lin .. (PCH)  Systlimalic ..... i tch  ｾ＠ tlIolp 
E requ en cy comp o 

Figura 50.  Suiche tipo  11  0  diodo 

91 
Para  efecios  de  simuiacion ,  al  diodo  se  Ie  coiocan  elementos  que  normaimente 
estan  presentes  en  un  sistema  real  (figura  51).  En  paralelo  can  cad a  diodo  se 
coloca  una  rama  tipo  "snubber"  (amortiguador) ,  formada  par  la  serie  de  una 
resistencia  y  un  capacitor,  tfpicamente  de  vaiores  oe1200  n  y  0. 1  iJ.F ,  Y  una 
resistencia  serie, que  debe  ser  pequena  can  respecto  al  valor  de  resistencia,  can 
la  cual  se  simula  la  carga  que  se  coloca  al  rectificador.  En  este  caso,  se  coloco 
una  resistencia  serie de  0.001  Q  (compa rar can  respecto  a !a resistencia  de  carga 
de 0.093  Q) . 

Sistema AC 
GaIga 

Figura 51.  Puente rectificador para  simuiaci6n 

EI  puente en  ia  herramienta Atpdraw puede quedar como se observa  en  ia  figura 

52 .  Los  datos  mas  relevantes  de  simulacion  son  los  siguientes: 

Tmax = 50 ms
.1t =  1  )..1s 
Plot Freq  =  3 

92  
 i!!l 
,/I 
/.\ 

ｾ＠

, , 

Ｗ3.5 Ｇ＠ ＰｾＱ ｖ Ｕ＠ ｾ ｦ
00 
­35 
­7 .0  .......,,) ' 10 ­3 
o  10  20  30  40  50 
­3.5 
­7 .0 +1ＭＮＬｲＢＧｾ＠
o  10  20  30  40  50 

Figura 52.  Puente rectificador en  Atpdraw  y variables de interes 

Para  este puente se  debe real izar minimo io  siguiente: 

•  Identificar  las  formas  de  onda  de  voltaje  DC  para  carga  resistiva .  Deducir  la 
expresi6n para  calcular el  valor DC  de  este voltaje . 

•  Identificar la  forma de  la  corriente  AC , para carga  resistiva. 

•  Deduci r los arm6nicos para  la  corriente AC , supon iendo carga  RL. 

8.5  PUENTE RECTiFiCADOR  CONTROLADO DE  SEIS  PULSOS 

Cuando  ei  puente  5e  encuentra  formado  por  diodos,  su  conmutaci6n  se  hace  en 
forma  natural,  pero  si  se  quiere  un  control  del  voltaje  DC  hay  que  controlar  el 
angulo  de  disparo  de  los  tiristores,  mediante  una  estrategia  de  control  de  los 
mismos. 

93 
 $,,,,\8"'01  AC  Ca rg;. 

Control de IIngulo de disparo 

Figura 53.  Puente redificador controiado 

Ei  puente  controlado  no  existe  directamente  como  modeio  en  ei  ATP.  Lo  que  se 
debe  es  desarrollar  un  modelo  mediante  alguna  de  las  dos  herramientas 
disponibles para  realizar funciones de control , que son  los TACS  0  los MODELS. 

Gille 1 

G,,1112 

Gala3 

ｾ｡ｴＴ＠

G"til5 

ＢＧＺＱｉｾＶ＠
L-.----

\ t1 ­=  ·1oe •
I 
\ t1 -= 60

Figura 54. Estrategia de control de! puente

94 
La estrategia de control se puede implementar can la opci6n TACS. EI c6digo
correspondiente se observa en el cuadro siguiente.

Tabla 4. Estiategia de contiOl para el puente en c6digo TAC S

ll DLYbUD  . UU2 7 777 7 8 

90 REFPOS 
9 0 REFNEG 
98VAC  =REFPOS­RE: FNEG 
l 
98RAMPl  58+VN ITY 
(')O r """,,,,·nl
­ I.,. .... ....... . r ｾ＠
­   =( r.A!'-1P l Ｎ ｉｬｊｾ ｾｅ＠ ｬ＠ ＲＮ･＠ ｃｾ＠ Ｎｐｾｎ｛［Ｌ＠ ｬｾｔＧｴ
120.00  0.0  1 . OVAC 
I
9 8DCMP1  5 4 +COMP l  5. 0£­ 3 
. NOT .  DCMP l  . AND .  COMPI 
98 PU LS l 
ｾ｢ｆ Ｎｵ ｾ Ｇ＠ ｾＴ＠ Ｚ＠ ｆ ｕｌｾ ｬ＠ DLYu0C  I
98PULS3  54+ PU LS2  DLY6 0D 
98PULS4  5 4+PULS 3 
­ DL Y6 0D 
Ｇ Ｚｈｬｐｕｌｓｾ＠ ｾ ＴＫ ｐｕ ｌ ｓＴ＠ D[,YoUD 
98PULS6  54 + PULS 5  DLY60D 
98GATE: l  PULS l  . OR.  Put,S;':> 
98GATE2  PUL S2  . OR.  pULS3 
98GATE3  PULS3 
­ .OR.  PULS4 
r\? IT !"" C
ｧ ･ ｇｲ ｾ ｔｃ Ｑ＠ PU LS1 .OR . r '"'..I..•...J.J

98GATE5  PUL S5  . OR.  PULS6 

98GATE6  PULS6 .OR.  PULS1 
­

Tado el puente formado par los tiristores , sus ramas snubber, resistencias serie y
la estrategia de control implementada en TACS , se puede agrupar en un solo
modelo al cual se Ie asigna un icono, como aparece en el cuadro siguiente.

95
Tabia 5.  fvlodeio compieto del  puente controiado en  c6digo TACS 

T ACS 
l  lDL Y60D  . 002777 778 
Q O RF.F'POS 
9 0RE FNE G 
98VAC  ­REF POS­ REFNEG 
ＹＸ ｈｊ ｾｉＮ ...""!P l  5e ..L.   T.JN !'!'Y  l.2C  . 'JO  0. 0  1  .  CI .....A ': 
98COMP l  =(RAMPl  ­ AN GLE_/ lSO)  . AND .  UN IT Y 
98DCMP l  S 4+ COHPl  S.OE ­3 
.92 P I.)"L3 :  - Ｎ＠ ＡＺｾ ｉ＠ ｔＮ＠ DCI.. .. ::? l  Ｎｾｬ ﾷ ｩｄ＠ Ｎ＠ ｃ ｏｾ ..l ==: 
98PULS2  ｓＴｾ ｐｕｌｓｬ＠ DLY6 00 
ＹＸ＠ ｾ ｕ ｌ ｓＳ＠ 5 4+ P ULS2  DLY6 0D 
9B PULS 4  S 4 + P UL S3  i)LY b \.iD 
9 8 PULS 5 ­ S4 +PUL S 4  DI.Y 6 0D 
98PULS 6  5 4+P UL S 5  DL Y60 D 
':J 8 GATE l  PUl,S  .OR .  PULs 2 
9 8GATE 2  PUL S2  . OR .  PU LS3 
98GATE3  PULS 3  .OR .  PULS 4 
98GATE 4  PUL S4  . OR.  PULS 5 
9 8 GATE S  PUL S 5  R .  PULS6 
9 8 GATE6  PUL S6  .OR .  PUL s  l 
I BRANC H
$ V I NTAGE,O 
POS  MIDl  Rsnub  Cs nub 
pas  MID3  POS  MIo l 
pa s  MID S  POS  MI Dl 
u  A..t·l1 i..i4  r 'U ::;  £-1'1 1.01
U  BM I D6  POS  MIDl  
U _ _CMID2  P OS  MI Dl  
HIDl  U  A  kserie  
MI  D3­U­­  BMI Dl  ­ U ­ ­ ­A  
MIDS  U  CMI DI ­ U ­ ­­ A  
MID4  NEG  ID l  ­ U ­ ­­A  
MI D6  NEG  MIDI ­ U ­ ­ A  
MJ!'2  N Er:  MI!'l  Tj  .".  
ｉｓ ｗｉ ｔｃ ｾ＠ ­ ­ ­ ­ -
l  lMIDl  POS  GATE l 
l11­..1 1 r ::  i: "(_'c.  Sh.7E:3 
l  l MIDS  pas  GATES 
IlMID4  U  A  GAT E4 
1.J.. t.fID6  U  D  ｇａｔｾｴＭ
ll MID2 ­ U­­C  GA TE2 
$EOF  User­supp l ied  h e ader  c ards  fo llow.  09­ Feb ­ 02  03 . 59.12 
ARG , U  A , U  ｾ Ｌ ｵ＠ C ,  POS  , NEG  ,REFP OS ,  REFNEG 
ARG , ANGLE  , Rsnub  , Cs nu b _  , Rseri e  ­-
Nm1,ANGLE  , Rsnub  , C snub  ,  Rserie 
DUM,  PUL S 1  ,P ULS 2  ,  PUL S 3  ,  PU L S 4  ,PUL S 5  ,PULS 6  ,HI Dl  ,MI 02  ,M ID 3 
UM,  GAT El  ,  GAT E2  ,  GAT E3  ,  GATE4  ,GATES  ,  GAT E 6­, VAC _ _ ,RAMP  ,  COM P l_ 
DUM,  DCMPl  ,DLY60D , MID'!  , MID':>  ,MI D 6 

La  ligura  55  ind ica  la  utiiizaci6n  del  modelo.  Los  parametros  de  este  modelo  son 
los siguientes : Angulo de  disparo , resistencia  de  la  rama  Snubber,  capacitancia  de 
la  rama  Snubber,  resistencia  en  serie con  los tiristores. 

96  
 Mas 

Menos 
It 

O:Il.  ｉ ｒｩｆｐｏｾ＠
.R.EFNf):;: . , 

ｾ ｯｵｰｎ＠ 10  Loll"·  r 

ｾ ｲｬＭ Ｍｾ Ｍ Ｍ

­ U ....Ｂ､ｉｯ ｾ＠ ""'7  ­ ­ r  H!iIJ 

'_do   ICI. "  \",,,d,"'"'USPI  ­ ­I ｐr  ｾ OId J­phG... 

Ｍ ＢＧ Ｇ＠
r 

Qonat.  tio", 

Figura 55. Utiiizaci6n  dei  modeio  de  puente  reGtificador  controiado  y  su  ventana  de 
dialogo 

Los  datos  mas  reievantes de la simuiaci6n  son  ios siguienies : 

Tmax  = 50  ms 

M  = 5  X  10­7 s 
Print Freq  =  25000 
Piot Freq :::  15 

Para  el  puente controiado se  debe reaiizar como  minima io  siguiente: 

•   identificar  ias  formas  de  onda  de  voltaje  DC ,  para  carga  resistiva  y  para 
diferentes angulos de disparo. 

•   Identificar la  forma de la  corriente AC , para  carga  resistiva , 

•   identificar  ei  desrase  existente  entre  ei  voltaje  de  rase  AC  de  ia  fuente  y  ia 
fundamental  de corriente.  Interesa en  este caso  identificar en  forma  cualitativa 
y cuantitativa  ei concepto  de  factor  de  poten cia  de  desplazamiento  y  ei  factor 
de  potencia verdadero , introducido por el  puente controlado . 

97 
B.6 MODELO ARMONICO 

Este  modelo  es  simpiemente  de  caiculo  y  io  que  hace  es  evaiuar  cualquier 
arm6nico de voltaje  0  de corriente  y el  correspondiente valor RMS . 

La  farrna  como  se  in:p!ementa  65te  iTlodelo  con  !a  herramienta  rv'ODELS  59  ind ica 
en  la  tabla  siguiente:  

Tabla 6.  Modelo ARMO en  !enguaje MODELS  

I  N  PUT  1  n 
D.A"rA 
ｾｔＧ Ｇ＠［ Ｃ］ＭＡＢＧＺｊ＠
'· .,.2'.... P
ｾｯｲＨｬ＠ 1  f'
ｲ ｏＺ＾ａｉｾｴ］＠

i  n  t:.  ｾｴ＠
.;  . . . .,  :"""  
t.Se  
1C03  
pro39Cl 

n  
l::. n  
C  r t    
.l.r"lt:;.l  
r.­n"" 
ＮＺｴｬｮｾ＠

O  LTT PfJT 
H I S T Or,y 
rrn.e> ,
.i.nr l 
;Lnt:.::d 
V"  ｣ｭＭｾ
ＡｆＢｊ｟ｾ
,  U  1: 'L,.'1 
ﾷｴ＠
［＠  

  0, 
U, 
IN IT  "'­1. 1,  0  .  0  
1,­ n  0  . 0  
1  n t .  ::l  0  _  0  
i  n r. .2  :  ­ 0  _ U  
..3 o g   0  . 0  
E  N  DINIT  
E  X EC  
Ｍ ｾＭ］＠ ｃ＠ ｬＢＺｩ＠ ｌ｣Ｚｬ＠ ｌＮＺｊ＠ ｣ｊ＠ ｦｾ＠ ｾ ｊ __  ｲ＠ ＮＩＰＧＺＢＱｾ［＠ Ｚｬ＠ u ­.:  ｌＮｩ＠ Ｇｾ Ｇ＠ ｬＮＭ＾＠ ＼ＭＮＩ＠ ｾＨ＠ ｌ｟Ｎｲ＠ Ｌ］＠ G  ­ T  
Lb  ｾ Ｍ ｾｲ ｵ ｮ｣Ｈ Ｍ ｦｲ｡ｱ ＩＯ＠ ｃ｣ｾ ｱ＠

Ｍ Ｍ Ｑ［ＮＸｮｾｲＳ｣ＺｩＶ ｪ＠ ､ｾ＠ ＡＮＺｊｕＢＦｾ＠ cl\w.!c  r: rt=cJ  V  __  I:...  ｲｵ､ｾ

: : : :  r  Ｈ ｾ ＾＠ 1  /  :I:L6·:­d ­1:..L"tl'''''::.<:jt..  t::. Ai'''::D  ＭＺｬ ＺＮ＼ｾ＠ Ｏ＠ 1::.1:.04 )  "rllDt.:t  I  ....  ｾ ......  .e  :.  ­1. _  u 
ｅｌｓｾ＠ ｄ｡ｾｰ＠ O . 
"ND I F  

｣ ｾ ｮＧ ｾ Ｍ Ｓ＠ ｬ ｾ ｌ ｾ Ｍ ｾ｣ｱ Ｍ ｮ Ｍ ｾＩ Ｍ ｾ Ｍ ］ｲｾ＠  
.L ｣ＮＺ＼ＧＢ ｾ＠ Ｚ｣ｊ＠ CCJ.3  (  Z  ....  3  Ｑ＠ ｾ＠Ｔ Ｖ Ｍ ｔ ｊＮ ･ ｱＭｲＮＱ＠ ....  L:)  ­ :2.  ­­ J:  .r.:oq  
ｰｲｾａ･ｮ＠ ' n  ­ fs=.  
ｲＧ ｴ ｾｏｑｕ＠ ｾｮ Ｍ ｌ ｣］ ｮ＠  

­ ­ I n t : e Q . t ­ ; ; \ c i  <;"n  
T';'P.LACF  (  tnt""­ :t  /  FJr ,_­.  s­­=­n)  l  /  (  ­L  "'  )  
T..... A.PLACF!  ( j  n t Z  /  p.L"JC0C)  ｾ＠ / ( l  s )  

­ ­ M  LJ. C  ;_:  t:.  :t  ＮＧＢ ｾ＠ Y  l:..­c..;rr .,J. "­.J  c)  d(.· ｩｮ ｴ＠ Ｇ ｾｃＧＢｊｲＮ＠ Ｎ ］ｴ＠ 1  
:r  F  Ｚｴ ｾＮ ］ｵｬ｜ｰ＠ ﾷ＠ ｾ＠ >  0  .  0  THP.N  ' 0   i  n t  J.  
on  int:.2  

ir"lt2  0  . 0  
E  N  DIF  
-Yl ＮｴＺＧＢｌＨ｡ｲ＠ ｾ＠ Ｂ＠ ...  r .... rl ­ ...  ?)  
r:::­"  ＭＺＧｾ＠ "':"::=: :  Ｌ ＧＢｴＮＭｾ＠ ＭＺｾ＠ ........  -
ｅｬＮｾ ｾ＠ ｉ ＺＧ＠ s n ...;;  :  u'Lb.n.2  (  bn, .,  "::an)  
RNDIF  
on  rnlr:, :  ­ e  ll ­ :CO:  .i.."'l  ( 2  ­ 3  _  J...  /l  1.  r.., - ｾ＠ rec <.:::r  ­ n  ­ t.­ .  r.1{ I'J)  
Ｍ ｾ Ｚ＠ ｾＺＮ［＠ ］＠ ＺＧ Ｍ Ｏ＠ ］ＺＭＢＮ ＺＨ＠ ｾ＠ ｾ＠

FND8 X  8C 
P­ N  I)MO Dc.;L 

Para generar un  modeio nuevo se  empieza como aparece en  la  figura  :iU . 

98 
 Qblects  Tools  Window  tjelp 

Edi! ｑｯｭｰ ｾＺ ｾ ［ Ｇ＠ lJ  [1'_  ｾ Ｍ［ ｰ

fji3;S EC  'iPMiLT­ ­ ­ ­ ­ ­ ­ ­ ­­ ­ ｾ＠

J..[s er Specified  ｾ］Ｍｆ Ｍ ﾭ T!.r'\.tT

OCJ'T'FU":' 
­.  t­ld• •   o  f  Inpyt:  v15rtaol .  . 

­ ｎ ｾＮ＠ o f   ＴＩ＠ ｾｉＬ＠ ｵＧｴＮ ＧｴＮ＠ V'4Irl .. to  t • • 

V «r­:..bte  1'1 4_ '  ｟ｾＬＭＮｴ＠ tfy  .... ,  Cp

New sup­fil"  DATA  ｉｾＢＬＮ ｶ＠ Ot  ､ＭｬｾＮ＠ v .• ,.­,.bl ... . 

EditIACS 
Edit s up­file ｉｖａ Ｎ ｾ＠ ­­ Hoe\•  •  ..,(  ,oc; Il1 ·  c u<:. pv ,,:  'N'l l obi 

HC'!'O R,{  ­­ j . Co. ,, :<: v ol" • •  of  v o. r1 41 b: • •  OD d.  .  x pr• •• ,:>r..  ( DFL.T : '1 }  

1)lIT  ­­ hill , .. I  ...... .  I  n , ,  

Editmod­1ile 
It;tJDUf  I  

!!:X!!:c  ­ ­ ＡＺｚＮ Ｍ］＠ Ｇ ｾ＠ ｴＮ＠ ｬＡＩ ｊｬ＠  

",:nfXlY1  
u  ;:lHOD!l 

l.:J 
­­' IiiJ""" 
ｾ＠

Figura 56.  Generaci6n de un  nuevo elemento con  MODELS 

En  el  editor  nativo  del  Atpdraw  aparecen  ias  partes  con  las  que  se  estructurarfa 
en  MODELS  un  modelo  general.  En  esta  misma  ventana  se  puede  editar  0  en 
otro  editor  de  archivos  pianos  y  se  guarda  con  la  correspondiente  extensi6n 
*.tv10D  en  el  directorio  ｃ Ｚ ｜Ｌｾｴｰ､｡Ｇ［ｊｲｶＱｯ Ｎ＠ La  fOima  con­tO  56  incorpora  e!  mode!o 
ARMO  dentro  del  ATPDRAW con  la  herramienta  MODELS  se  ilustra  en  la  figuras 
57  y  58 . 

2robes &  3­phase 

§ranch linear  fiuscaren: 
Branch ｾｯｮｬｩ･｡ｲ＠
ｾ ａｲ｣Ｓｰ＠ ｾａｓｄＹＴ＠ ｾｬｭＳｰ＠ ｾｒ･ｧｵｬ ｡＠
]..ine Lumped  GCopiade TAPTRAFO  2Jlm3p  2lREGULA 
dlArc3p 
line .o.is\ributed 
ｾ＠ GOetault  Glnd94  BRegulv

ｓｾｴ｣ｨ･ｳ＠
G  armonlco  I2l Default  ｾ ｉｎＰＹ Ｔ＠ .] Regulv 
0  armon ico  GOist1  GMedPO  Grele sup 
ｾｯｵｲ｣･ｳ＠
eJAsd6p  eJ OISTl  2] MedPO  2lrelesup
2lAsd6p  G  Flash_l  8Modos  8Rms 
MJ:!chines  .JJ Flosh_'  2J Modos  2JRms 
­­­­­ GAsd94 
Transformers  <4 
.!.1 
t::!ombre de archive:  IArrno  6brir 
ｔｹｰ･ｾ Ｔ＠ ｾ

IACS  Iipo de archives :  I M­o­d­e­If­il­er­­.mod­­

r ­ ­ ­­­­­­­­­ij­
­ ;"'.s­up) ....­. Cancelar 
!.tser Specified 
ｕｮ･Ｏ ｾ｡ ｢ｬ･＠
Qverhead Line (PCH) 
frequency compo 

Figura 57.  Incorporacion del modeio ARiviO  dentro de ia  iibreria de MODELS 

99  
 ­i ')   Model successfully ｩ ､･ｮ ｾ ｦｩ ･､＠ I ­ SUlflde.'d de.t!l­ ­-
"V  Input­I. Ou\put=2. De.te.· 2 
Edilfile?  .r  High pre cisi an  j 
r  Output enable 
ｉｮ ｾ Ｎ ＺＱ ｟ＭＬ＠ t!lo  Cencel I 
IMo del  ｾ＠

12±l l' 
Ｍｾ＠
f3:B 
ｾ ｾ＠
.­Standard data 
r  High  ｰ ｲ ･ ｣ｩ ｓ ｾ＠  
erma 
vrms 
'q"" '
32
• 
"  " , 
I r  Output ene.ble 
TVOE   IModel  iJ 
S l!V9  I Save ｾ＠
｣ＺＮｾ Ｌ ｾ＠
I  ｾ Ｚ＠ ｆ ｾ＠ I i t］ｾ Ｎｬ＠ ｟ ｉ＠ ｌ _ｾ｟ ｄ ｟ . . _' _t!=el=p ｾｉ

ｾ ｾ＠
ｾ ｓ ｡ ｶ･ｾ＠ I
•  ｉ ＨＭ .  . ｾ｝＠ .__..,  t!el p

Figura 58. Definicion de datos  y nodos dei  rnode io ARivl0 

Utiiizando el  modelo ARiviO  determinar: 

•   La  descomposicion arm6nica de ia corriente del  puente 

09  Calcuiar ei  factor de distorsi6n  para  !a  corriente. 

100  
 11  15  17  111 

t. 

"1
Mas 

Ｌｾ＠ 1:\·­:..- L-( ..

Fuente  ｾＭｴＩ Ｑ＠ Puent, "

. _ 
+ Menos 

. D

t.. 

Figura 59. Utilizaci6n de modelo ARMO en descomposici6n de corriente del puente 

Corrriente (kA) 
7 ill l .... ｟ ｾ＠ ...&&iC 7"1

525  _ .. __ _ . 

3BJ 

1.75 

am
­1.75 

­3 BJ 

­:525 

-7 m I  '>.E - , , ., '=-"" i  "10-J I 

15.00 2222 V .77 33.33 33.00 44.44 tllDO

tiempo ems) 

Figuia 60. Resuitado de ia  descomposicion armonica de ia  corriente 

·1 rH
I V  I 
8.7  DiSTORSION DEBIDA AL PUENTE Y PRESENCIA DE  TRANSFORMADOR 

El  verdadero  probiema  por  ia  presencia  de  corrientes  armonicas,  es  que  estas 
circulan  hacia  la  fuente  por  elementos  de  la  red  que  tienen  impedancia 
(transformadores  y  lineas) ,  generando  distorsiones  que  afectan  la  forma  de  onda 
del  'voltaje  en  if,ayor grado  mientras  mas alejado  se  este de  18  fuente .  En  la  figure 
61  el  transformador de  potencia  se  ha  referido  al  lade  de  baja , de  manera  que  se 
pueda observar la  distorsi6n del  voltaje antes  y  despues del  transformador . 
..., 
300  .,  I I
1 11  I  1\ 
,.. , lL 
200 30 

If 
Ｌ＠ ｾ＠ \ 
·'00 
"  \  \. 
."" 
. 
ｾ＠

·300  /  7 
­400  ' i  5ll  '"  . 
rV 
23 81  ｾ ＧＮ＠ Ｘｊ＠
J ,. ,
.6'00 
... ｐｔﾣ ｐｾｽＢ ＭＧ ｗＢＧＮ Ｇ ｾＮ＠ ＢＧ＠ Ij  ｾ ｐ ｔｅ＠ Mas

. - ﾷ Ｈｴｲ［ＮＭﾱｬ＼ ］Ｍｬｩｾ＠ Menos

'DO 

lOO 

"'"
"'"
.,,,, 
· 200 

· 300 

....,
l i50  13'&7  Ｚｗ ｾＸＳ＠ .., 00 
ｾ PTEJTftAFO_REF _WA.plo&: X'I'f( I) .." ... _a..­

Figura 61. Distorsi6n arm6nica del voltaje antes y despues del transformador

El transformador se puede representar en forma real conectado a sus respectivos

niveles de tensi6n nominales y teniendo en cuenta la forma de conexi6n de sus
devanados primario y secundario. (Ver figura 62)

102 
 Mas 
u 
ｾ＠
ＫＮＭ
'I  ,.  LIN. 
Itt... ( A')
_ )L  ＬＭＮＧ ｾ ［Ｎ Ｚ Ｍ Ｍｾ Ｈ Ｍ Ｎ ＭＫ ＭＫ＠
Menos 
Fuente 

Figura 62.  Diagrama en  Atpdraw inciuyendo 6i  transformador potencia 

8.8  DiSTORSION CON LA CONEXION DE BANCO CAPACITiVO

Para  ilevar  ei  factor  de  potencia  a  un  vaior  de  0,93,  ei  banco  capacitivo  puede 
generar  problemas  de  resonancia  a  frecuencias  cercanas  a  los  arm6nicos 
presentes  inyectados  por  el  puente,  Se  trata  de  identificar  el  arm6nico  que  esta 
presentando  amplificaci6n,  para  impiementai  medidas  iemediaies  como  son  los 
filtros tipo shunt 

103 
Vo ltaJes  en  el  banco  [V] 

voltaje  con  banco 

,,,,
10  15  20 25  .;)u  35  40  45  50 
Tiempo  [ms] 

Figura 53. Formas de onda del voitaje en  baja con  banco  y sin  banco capacitivo 

8.9 NOTAS ADiCiONALES

Adjunto  a  este  documento  Se  envlan  ios  archivos  base  denominados 

ARMBASE.ADP y ZWBASE.ADP, que  corresponden  al  sistema  industrial 
modelado,  y su  correspondiente archivo  para  utilizar la  opci6n  Frecuency  Scan  en 
480 V. 

Ei  archivo corresponde ai  sistema  industrial en  condiciones de demanda  maxima  y 

sin  la  instalaci6n  del  banco  de  condensadores  que  se  utilizara  para  mejorar  el 
factor de  potencia  de la  planta. 

Ei  sistema  fUe  implementado  con  ia  interiaz  A TPDRAVv.  Para  efectos  de 
simulaci6n ,  se  incluyen  los  respectivos  modelos  del  equivalente  del  sistema , 
transformador de  5 MVA,  carga  lineal y puente  rectificador de 6 pulsos,  de acuerdo 
con e! disgrama unifi!ai dado en la figure 64.

104 
 MIs 
!' 
,..,.. 
M!nt:e 

L 
2 WMI h  '1  1­.'1  2. ZJM.IAA

Figura ＶＴ ｾ＠ Sistema base imp!ementado en ATPDRA'/J

B.l0 CONTENiDO DEL iNFORME

Ei  informe finai  debera ser estructurado de  ia  siguiente forma : 

'1.  introduccion 
2.  Descripcion  del  sistema 
3.  Formulas empleadas 
4. Ca!culos reallzados
5.  Comportamiento del  sistema actual 
5.1  Grafico 
5.1 Formas de onda y espectro arm6nico en barras de 4e,o \/
5.2  Niveles de distorsion armonica 
5.3  impedancia armonica en  barras de 480\:' 
6.0  Compensacion  reactiva 
6.1  Reactivos a compensar 
6.2  Sistema compensade (grance) 
6.3  Formas de onda  y contenido espectral  en  barras de 480 V 
6.4  Niveles de distorsion 
6.5  !mpedancia armonica 

-105
7.  Diserio dei fi itro 
7.1  Metodologfa 
7.2  Calculo 
7.3  Sisterna  sintonizado 
7.4  Impedancia arm6nica 
7.5  Formas de onda  y contenido espectral en  barras de  480 V 
7.6 ｾＴ＠ Ａ ｶ･Ａ ･ ｳ＠ de distorsi6n arm6n ica
8.  Conclusiones 
9.  Bibiiograffa 
10.  Anexos  (adjuntar solo archivo *.atp caso  base) 

106 
 BIBUOGRAFIA

1.   CARDONA  CORREA,  LEONARDO.  Teorfa  y  Prcktica  con  el  ATP . 

Medellin , Colombia .  Universidad  Nacional  de  Colombia ,  Faculta  de Minas, 
-1995. 270 p.

2.   BEATY  H.  VvAYNE,  DUGAN  ROGER  C.  Y  McGRANAGHAN  ivlARK  F, 
Electric Power System  Quality.  New York ,  1996.  264  p, 

3.   CRESPO ,  FRANCISCO.  Sobretensiones  en  las  Redes  de  Alta  Tensi6n . 
Madrid ,  Asociaci6n  de  Investigaci6n  Industrial  Electrica  (Asinel),  1975. 
222  p. 

4.   F1LHO ,  JORGE  AMON .  PEREIRA,  MARCO  POLO ,  Atp  Alternative 

Transients  Program,  Curso  Basico  sobre  a  UtilizaC;ao  do  Atp.  Sao  Paulo, 
BrasiL  1994. 

5.   FURNAS ,  Curso  sobre  Transltonos  Liectromagneticos,  Coordinaci6n  de 

Aislamiento y Programa  EMTP.  Medellin,  Colombia.  1985. 

6.   FURNAS.  Transitorios  eieciricos  e  coordinac;ao  de  isolamienio,  apiicac;ao 

em  sistemas  de  potencia  de  alta­tensao,  Rio  de  Janeiro.  Universidad 
Federai  Fiuminense.  1987. 

7.   GREEN\NOOD,  ALLAN.  Electrical  Transients  in  Power  Systems ,  New 

York ,  Jhon Wiley  & Sons, Inc,  1991 , 750  p. 

8.   HERMANN  VV.  DOMMEL ,  Electromagnetic  Transients  Program  Reference 

Manual , EMTP Theory Book . Boneville Power Administration , BPA.  1986, 

9.   VvESTiNHOUSE  ELECTRiC  CORPORATION .  Electromagnetic  Transients 

Program  (EMTP) , Primer.  Electric Power Research  Institute,  EPRI.  1985, 

10.   \NEST1NHOUSE  ELECTRiC  CORPORATiON.  Electromagnetic  Transients 

Program  (EMTP) ,  Aplication  Guide,  Electric  Power  Research  Institute, 
EPRi.  1986, 

ＮＨＧ｜ｾ

I \J {
ANEXO A. DETERMINACION DEL EQUiVALENTE DE THEVENIN TRiFAslCO
A PARTIR DE LOS NIVELES DE CORTOCIRCUITO

Si  se  quiere  determinar  ei  equivalente  de  Thevenin  trifasico  a  partir  de  ia 
informacion  de  los  datos  de  cortocircuito,  hay  que  tener  en  cuenta  que  las 
expresiones aproximadas  para  sistemas dominantemente reactivos  pueden  no  ser 
aplicabtes , especial mente  cuando e! equiva!ente corresponde a nive!es de vo!taJe
de  distribucion,  donde  la  componente  resistiva  empieza  a  tener  importancia 
ralativa  con  respecto a la  componente reactiva. 

La  informacion  necesaria  para  determinar  ei  equivaiente  de  Thevenin  trifasico  en 
algun  punto de la  red  es  la siguiente: 

hX.:3rP =Corriente de cortocircuito trlfasica 
Iccup =  Corriente de cortocircuito monofasica 
(X I R)3¢ =  Relacion  XIR que ve  la  corriente de corto trifasica 
(X / R)1¢ :::  Reiacion  XiR que ve  ia  corriente de corto  monofasica 

A  partir  de  los  datos  y  considerando  ei  voitaje  nominal,  ias  expresiones  que 
permiten calcular los parametros de secuencia son  las siguientes: 

F 
R ­ r====
' L (1 ) 
I  ­ ｾ＠ Ｏ Ｓ＠ x  Icc3r/J x  ­,./1+ (X / R)23¢

(X / R)3¢ XVL
X= =
I  ­/"3x  Iccvp x  .i1+ (X / R)23¢ (2) 

ｾ = _  ­.,'3  xV,. LXV

ｾ＠ ,­ 2 
!cC1¢ x  ,11+ (X I R)\£ 3 x  !cC3¢ Y  •  .  }  + (X / R) 3¢-
(3)

_ (X i R)1¢-'3xVL 2x(X I R/3¢ XVI.

Xo ­ = Ｎ＠ Ｍ Ｍ］ＬＮｾ (4 )
f ccup x . 1+ eX / R)2!"./3 x lccJ¢ x \ 1+ (X / R/3¢

108
Las  expresiones  para  X 1 y  X o,  cuando  las  reiaciones  XiF? son  altas, tienden  a  ias 
formas  clasicas : 

T,'
YL
(5) 
Xl =  3X /CC3¢ 

ｪＯｾ 2 x  ｪＯｾ＠

, 3 x ­­ ­ (6) 
Xo =  f CC1¢ 3 X  Jcc.1rjJ

Para  las  aproximaciones,  se  considera  que  reiaciones  altas  son  las  que  estan  por 
encima de 3.16,  10  que da  una  precision del  95% . 

E!  mode!o  mas  adecuado  y  practico  pam  representar  un  equiva!ente  de  Thevenin 
trifasico  en  el  ATP  es  el  de  tipo  de  elementos  concentrados  RL,  en  la  opcion  de 
componentes simetricas (Ver figura  1) 

.Erobes  & 3­phase 

!2ranch Linear 

RLC Pi­equiv. 1  .1 
Lins Qistributed  •  RL Qoupled 51  • 
r  ｾ＠ 'L.  (  I,  :­.  • 

Sn:itches 
.6.  ph  ｾ＠
Sources 
ｍｾ｣ｨｩｮ･ｳ＠
)It 

Alribute;. 1
Transformers 
ｄ ｾｬＮａ＠ VIltiJE  NODE ' PHASE  IllAME 
MODELS  B2..­..  a  ｾ ｟ Ｑ＠ _­.­1 3. .  .• 
La  : H2  Ot,!T1
_  J_ __ 
ｔｹｰ･ｾＴ＠
R·  ｩ＠ ｾ＠
pes  L·  11  0  X Ｚ［＠ Ｇ＠ Ｏ Ｓ＠ ｸｬ ｾ＠ ｟＠ ｾ＠ Ｒ＠ ｸ＠ ｖ ｾ＠
­ Q JcClIP ,.r3 X  JcC 3¢
!.!ser Specified 
_ V_L__
Une/Qable 
Ｓ＠ ｸ＠ ｬｾ＠
Qv erhead Line (PCH)  llroup No  [0  Lalla\:  I  .
Erequency comp o  ｃ ｯ ｭＮｾｲｬＭ Ｍ Ｍ

r  Hill_ 
r  · 

ｾ＠ IlK  !;en",,1  tjolp 

Figuia 1.  Ventana de  datos para ei  equivaiente de Thevenin trifasico 

109 
 ﾷｾ＠
ｾ＠ ­/ ..... ｾＮ＠

­ .. ·Yn"" , \/' ｾＮＺ＠ .. .  (­_ ·hll\.\f. 

U·" . 
j)f 
.  IN 
C UU)M RlA

DEPTO.  DE  RIRLTOT ECAS 

BIBLlOTFCA  MINAS 

ANEXO B. SOLUCION DE LA ECUACiON DE ONDA VIAJERA

Una  linea  de  transmisi6n  eS  fundamentaimente,  una  red  RLCG  de  elementos 
distribuidos,  donde  la  energfa  se  va  propagando  de  elemento  en  elemento, 
generandose  10  que  se  llama  una  onda  viajera.  Una  primera  aproximaci6n  al 
fenorrleno de la propagaci6n es representar !a Hnea por una serle de elernentos
RLCG  en  cascada. 

Para  determinar  en  forma  exacta  ei  comportamiento  de  ia  linea  en  regimen 
transitorio,  se  plantean  las  ecuaciones  de  la  telegraffa ;  ecuaciones  que  se 
deducen a partir de  un  elemento de linea como  el  de la  figura  1. 

z, 
L
ｾ＠
r.(I  •, 

+ +
R .m L m ­
i ＨｾＮｴＩ＠ :;(x+6x.:)
ｾ＠

+ ｾ＠ ｾ ＭＮＺｮ＠
-
!:J(t}

Vf(l) vdt) ,,(x, I)  GLlx,,"  >1< Cilx vlx+iltl) ｾ ｾ Z, 

< m ｾ＠
I
x =O x =d
Figura 1. Circuito equivalente de linea con pan3rnetros distribuidos

ｒ ｾ＠ Resistencia por unidad de iongitud

L -7 Inductancia per unidad de !ongitud

G-7 Conductancia por unidad de longitud

ｃｾ＠ Capacitancia por unidad de longitud

En el eiemento de ifnea se piantean ecuaciones de voltaje y de cornente de

acuerdo con las leyes de Kirchoff.

i 10
De ia  ley de voltajes  : 

v(x, t) ­ v(x Ｋｾ Ｌ＠ t) = R.llx.i( x, t) + L.llx. ai(x,t) (1 )

at

Dividiendo  por  t1X  en ambos miembros de  ia  ecuaci6n  Cl )

v (x , t ) ­ v (x + ill<, t) = R.i(x t) + L _a i (x, t ) (2) 

ｾ＠ ' at

Tomando i[mites cuando  ｾ ｾ＠ 0 

av(x,t) = -R.i ( x ,l ) ­ L---at

ai(x, t )
(3)
ax

De ia  ley de corrientes : 

i ( x ,t) ­ i (x + Llx,t ) =  G. Llx. v (x + Llx,t ) + C. Llx. ov( x + LU, t) (4)
ot

Dividiendo por  ｾ en ambos  miembros de la  ecuaci6n (4) 

i(x+Lll,)
i(x,t) - fu  ｾ＠ ,
[\  = G .v( x + fu,t ) + C ov\x + Lix,t) (5)
ot

Tomando limites cuando ｾ ｾ＠ 0 

a
oi (xx, t ) = -G.v (x , t) ­ C 3v( x, t ) (6)
at

Las  ecuaciones  (3)  y  (6)  se  deben  resolver  simuitaneamente .  Existen  diferentes 
metodos para resolver las ecuaciones de la  Ifnea : 

•  Ei  metodo  de  D'Al embei l 

: Et metoda clasico

11 1
•  EI  metodo de la  transformada  doble de Laplace 
Los  dos  primeros  metodos  son  los  mas  utilizados  porque dan  una  solucion  simple 
y parten  de un  supuesto de  solucion  que cumpla  con  las condiciones  iniciales y de 
borde del  problema . 

EI  metodo  de  la  transformada  de  Laplace  es  mas  elaborado  en  el  proceso  de 
encontrar  la  solucion , pero  tiene  la  ventaja  de  que  la  interpretacion  de  la  solucion 
final  es  mas  sencilia .  Otra  ventaja  es  que  la  suposicion  de  casos  particulares  de 
solucion ,  como  la  de  una  linea  sin  perdidas,  solo  hay  que  hacerla  al  final , 
perm itiendo un  paso facil  hacia  los  modelos de IIneas que consideran  perdidas . 

En  este documento se  utiliza el  metodo de la  doble transformada  de  Laplace. 

SOLUCION MEDIANTE LA TRANSFORMADA DE LAPLACE

Aplicando  la  transformada  de  Laplace  a  las  ecuaciones  (3)  y  (6) ,  con  respecto  al 
tiempo : 

:£  [c{f(t)] = sF(s) ­ f( t = 0)  (7)

dt

dV(x,s) = -RJ(x,s) ­ sLJ(x,s) + Li(x,t = 0)  (8)

dx

En  la  ecuacion  (8)  i(x, t=0) = O.  Esta ecuacion  queda: 

dV(x,s) = -(R + sL)J(x,s) (9)

dx

dV(x,s) = -GV(x,s) ­ sCV(x,s) + Cv(x,t = 0)  (10) 

dx

En  la  ecuacion  (10)  el  voltaje en  t=O  es  cero,  ya  que la  onda  no se  ha  propagado , y 
la  ecuacion  se  reduce a: 

112 
 ､ｬＨｸＬ ｾ = -(G + sC)V(x,s) (11 )
dx:

Oerivando las ecuaciones (9)  y  (11)  con  respecto a x: 

2
d V(x,s) =  (sL + R)(sC + G)V(x,s)
dx 2 (12) 

2
d lex ,s) = (sL + R)(sC + G)I(x,s) ( 13) 
dx 2

Aplicando la  transformada  de  Laplace a la  ecuacion  (12)  con  respecto  a  la  variable 
de posicion x.  ;£ [f(x)] =  F(p)

:f [d 2 f 2ex)]  = p 2 F(P) _ pj(x =  0) _  dj(x = 0)  (14 ) 

dx dx

p 2 V(p,s) ­ p vex = O,s)  + (sL + R)I(x = O,s) = (sL + R)(sC + G) V(P,s) (15)

Oefiniendo  : 

y2 =  (sL + R)(sC + G)  (16) 

Oespejando el  voltaje  V(p, s)  : 

p .  sL +R
(p,s) = ­ 4  ­ i­ ­ VI(S) ­ ­ 2  2- J I(S) (17) 
p ­ y (s) p ­ y (s)

Aplicando  transformada  inversa  de  Laplace  a  la  ecuacion  (17)  con  respecto  a  la 
variable  x. c'l - 1[(F(p, s)]  ­) F(x, s) 

:L ＭＱｬ ｾ
p  ­a  
2]  =  cosh(ax);  :L ­1 [2  1­2]  = ｾa Ｎ＠ senh(ax) 
p  ­ a 
(18)

Y definiendo, 

113 
 Zc(s) =  , I sL + Rｾ＠ =  Impedaneia Caraeterlstiea
' (19)
sC+G

Vex,s) = V/s) cosh (y(s)x )- Ze(s)l /s)senh(y(s)x)  (20)

Con el mismo procedi miento se lIega a la ecuaci6n para la corri ente :

l(x,s) = l/s)cosh(y(s)x)­ VI(S) senh (y(s)x) (21 )

Ze(s)

Las ecuaciones (20) y (2 1) al utilizar las expresiones del cos h y senh , toman la
siguiente forma:

VI (s) + Z  (.\')11 (s) -"(s) ,  VI(S) ­ Ze(s)l l(S) '(s)x

Vex s) =  e e ' .  +  e' (22)
' 2   2 

1(x s) = 1I ( s) + ( VI  (s ) / Z c (s )) e -)'(s)x + 1I  (s) ­ (V; (s) I Z e ( 5 )) e Y(S)X (23)
' 2   2 

Definiendo :

V1 (s)+ Zc(s)11 (s)

A =  ｾＭＮＺＧ   (24)
2 

B =  VI  (s) ­ ZcC.5)11 (s) (25)

2 

Las ecuaciones (22) y (23) quedan de la siguiente forma:

Vex,s) =  A e -;(SlX + Be )('<)x (26 )

l(x,s) = _ A _ e -iiS)x _  ｾ ･ ｽＨＧ ｬｸ＠ (27)

Z e(S) Zc (s)

114
Para  resolver,  finalmente  las  ecuaciones  para  el voltaje  y  la  corriente  en  cualquier 
punta  de  la  linea,  se deben,  calcu lar los coeficientes  A y B  de  las ecuaciones  (26) 
y  (27).  Para  evaluar  esos  coeficientes  se  plantean  las  siguientes  ecuacio nes de 
borde : 

Vj(s) =  Z JI (s) + Vl (s) (28) 

V2 (s)= 2 2 12 (s)  (29) 

Vls')= A+B (30) 

1 (s) =A-B (31 )

1
2 C 

V (s) =  Ae -}{.<)ti + Be )\s)tl (32) 

2

12(S )= Ae - )\s)tl _Be )'S )tl

(33) 
?c(s )

De (28) , (30)  Y (31) : 

VI(s ) =  21  + Zc A + Zc ­ Ｇ ｾ＠ ｬ＠ B (34) 

2 c  Zc

De  (29),  (32)  Y (33): 

B _  Z2 ­ Z  2i
(35) 
­ 2  + Z C x  A x  e- ( 5)t1

2 ｾ｣＠

Resolviendo  las  ecuaciones  (34)  y (35) se  lIega  a  las expresiones fi nales para  A y 

B. 

A =  V[(s) x  ­ c  (36) 

115 
 -Z){s)d
Zc r2 e
B = VI(S) X  x  -1----=-- ---=-
2r{-:-:s)d--;- (37)
J ZI  + Zc -rl r2 e

Donde, 

rl = 
z ­z 
Ie,
. 
Coejiciente de reflex ion af comienzo de fa linea (38) 
Zl  + Zc

Z  ­ Z 
r 2 =  2  c , Coejiciente de r eflex io n af jinal de la l inea (39) 
Z 2 + Zc

Reemplazando A y Ben (26) : 

Vex,s) = VI(s) Zc
Zl  +zc
x [e -i'(S)X + r 2e- 1.. )( 2d-x)
l - r l r z e-1.. S )2d
S
1  (40)

Con  el  uso  del  teorema  del  binomio  se  puede  descomponer  una  parte  de  la 
ecuaci6n  (40)  en  una serie infinita, de la  siguiente forma ; ya  que  r l r 2e- y(s)Zd < 1 

1  2 3 n
l- r r e- ){S)Zd = I- a =1+a + a + a +  ..  a  +  (41 )
l 2

1 - }(s)Zd =1+ r r e-1.. S)2d +r Zr 2e- ,(S)4d +r 3 3 -i(,)6d

-
l z
r e (42)
r r
1 l 2
e I  2  I  2  +

Reemplazando (42)  en  (40) , 

Vex,s) = VI(s) ｾ＠ z  c  x 
[e-)(*Z _+r 2e- / <)( Zd-xJ +r r Ze- )(s)( 2dH) +] l
(43) 
J  + Zc r r e  )«)(4d-x) + r Z Z -)(s)( 4d+x)
I  Z I  rz e  + 

La  ecuaci6n  para corriente (27) , despues de reemplazar las expresiones obtenidas 

para A y B queda: 

l ex,s) = )/f(s) ｾ Ｎ＠
I+ Z c
x ｛･ Ｍ ｩＨ ｓｾ＠ ­ r2e-)'(,f 2d- x)
l_rr e- i {.,) Zd
X
1  (44)
I  2 

116 
Reemplazando (42)  en  (44)  se obtiene la  ecuaci6n final  para la  corriente  I(x,s). 

I(x,s) ==  Vf(s) - _ 1_  x  [ e-V(S)X ­ r2 e - V(S)(2d-X) + r,r2 e- V(S)(2d+x) ­] 

(45)
Z,  + Zc  r,r22 e - V(S)(4d-X) + r, 2 r2 2 e - V(S}(4 d +x) ­ ...

SOLUCION DE LA ECUACION DE ONDA EN EL DOMINIO DEL TIEMPO PARA

UNA LiNEA IDEAL SIN PERDIDAS

La  transformada  inversa para  las ecuaciones  (43)  y (45)  solo se puede calcular en 

casos  particulares.  Interesa en  esta  oportunidad  resolver el  caso  de  una  linea  sin 
perdidas  (R=O,  G=O). 

Para  el  caso  de  la  linea  ideal,  la  impedancia  caracterfstica  Zc no  depende  de  s  y 
se convierte en  una constante. 

It 
Zc(S)== oJ SL+R == 
sC+G  tc  (46)

La  expresi6n  e-r(s)x se convierte en  un  termino de la forma  e - r.r ,

e -r (S)x ==  e - ,,I( SL+R)(SC +G ) x :::;  e-s I "Lcx

(47) 

La  transformada  inversa para el  termino generico que aparece es : 

­ l
:f f(s)e- ST ]= f(t ­ r)u(t ­ r)  (48) 

La  soluci6n  anterior  se  interpreta  como  una  onda  desplazada  en  el  tiempo  que 
conserva la  misma forma a  10  largo de  la  linea . 

EI  termino  . lLex corresponde  al  tiempo  de  propagaci6n  de  la  onda  hasta  la 
distancia  x.  Es  usual  expresarlo  en  funci6n  de  la  velocidad  de  propagaci6n ,  en 
lugar del tiempo de propagaci6n . 

117 
 v  = ｾ Le
1 
= Velocidad de propagaCion[Unidades de distanCia]
segundos (49)

La  soluci6n  para  el  voltaje,  despues  de  aplicar  la  transformada  inversa  a  la 
ecuaci6n  (43),  con  las consideraciones de  linea ideal , es: 

(1  + ( 2  )v, (t­ r )u(t ­ r)+

Z  (1(2(1 +(2)v ,(t-3r)u(t-3r) +
V 2 (x,t) =  _ c _._ X  2  2  (50) 
Z1  +Zc (1  ( 2  (1+(2)v,(t-5r)u(t .,- 5r)+

3  3 
(1  ( 2  (1  + (2)V, (t ­ 7r )u(t ­ 7r )+ ... 

v, (t ­ xlv )u(t ­ x/v)

­ ( 2 v, (t ­ (2d ­ x)/ v)u(t ­ (2d ­ x)/ v)

1 x  +(1(2V,(t - (2d +x)/v)u(t -(2d + x)/v) (51 ) 

i(x,t) =  Z  + Z 
1  c 
- (1(22V,(t ­ (4d ­ x )/v)u(t ­ (4d ­ x)/v)
+ (1 2 (22  v,(t -(4d +x)/v)u(t -(4d +x)/v)- ...

Las  anteriores  ecuaciones  dan  la  soluci6n  para  el  voltaje  y  la  corriente  en 
cualquier  punto  de  la  linea  e  instante.  Una  soluci6n  de  interes  practico  es  el 
voltaje  al  fina l de  la  linea ; para este caso  la  ecuaci6n (50)  se  convierte en: 

(1  + r2 )v f (t ­ r )u(t ­ r )+ 

r1r2 (l + r2 )vf (t ­ 3r )u(t ­ 3r) +
Zc
v 2 (x,t) =  2 1 +Zc
x  I  2  2 
r1 r2 (l+rJv f ( ) 
t-5r ut( 
- 5r ) +
(52)  

r1\ 3
2 (1  + r2 )v f (t ­7r )u(t ­7r ) + ... 

Algunas  situaciones  particulares  ayudan  a entender la  soluci6n  de  la  ecuaci6n  de 
onda .  Una  de  elias  es  considerar  una  linea  de  longitud  muy  grande  y  tratar  de 
darle  una  interpretaci6n  para  este caso a las ecuaciones (50)  y (51). 

Con  las  ecuaciones  (50)  y  (51)  se  quiere  evaluar  el  voltaje  y  la  corriente  en  un 
punto  de  la  linea  a una  distancia  x, donde  se  cumple que  la  distancia  d,  es  mucho 

118 
mayor  que  x  (d»  x).  Para  este  caso,  las  ecuaciones  para  voltaje  y  corriente 
quedarfan  reducidas a: 

vex,!) =  Zc  x[v/ (I ­ x' v} u(t ­ xlv)] (53) 

ZI  +Zc

i( x,t) =  1  X[V/ (/ -X 'V)U(t - x'v)] (54) 

ZI  + Zc

De  acuerdo  con  las  ecuaciones  (53)  y  (54),  el  voltaje  aparece  como  una  copia  a 
escala  con  respecto  al voltaje de la fuen te, un  tiempo despues de  la  energizacion, 
dado  por  xlv;  es  decir,  el  tiempo  que  demora  la  onda  viajera  en  recorrer  una 
distancia x a una velocidad  v, 

119