INDICE DE TRANSITORIOS ELECTROMAGNÉTICOS

1. NOCIONES FUNDAMENTALES.

1.1.¿PORQUÉ ESTUDIAR LOS TRANSITÓRIOS?

1.2. AMBITO DE LOS FENÓMENOS TRANSITÓRIOS

2. DEFINICIONES

2.1. PARAMETROS BÁSICOS

2.2. COMPORTAMIENTO FÍSICO

2.2.1. SOLUCIÓN MATEMATICA

2.2.1.1.RESPUESTAS EXPONENCIALES
2.2.1.2.RESPUESTAS OSCILATORIAS

2.2.2. INTERPRETACIÓN DE RESPUESTAS TRANSITÓRIAS

3. MÉTODOS DE TRANSFORMADA DE LAPLACE PARA SOLUCIÓN DE

ECUACIONES DIFERENCIALES

3.1. MÉTODO OPERACIONAL

3.2. SOLUCIÓN DE CIRCUITOS “R-L”
3.3. SOLUCIÓN DE CIRCUITOS “R-C”
3.4. SOLUCIÓN DE CIRCUITOS “L-C”
3.5.ESTUDIO DE LOS CIRCUITOS “R – L – C”

4. SIMPLES TRANSITÓRIOS DE CONEXIÓN.

4.1. INTERPRETACION FISICA

4.2. SOBRETENSIONES EN LA CARGA
4.3. RESISTENCIA DE ATENUACIÓN
4.4. SOBREVOLTAJE TRANSITÓRIO

4.4.1. INTERRUPCIÓN DE CORRIENTE ASIMÉTRICA

4.4.2.RECUPERACIÓN DE LA TENSIÓN DE ARCO EN EL
INTERRUPTOR LUEGO DE APERTURA
4.4.3. SWITCHING DE CONDENSADORES
4.4.4. TRANSITÓRIO DE DOBLE FRECUENCIA

5. ONDAS NAVEGANTES Y OTROS TRANSITÓRIOS EN LÍNEAS DE

TRANSMISIÓN

5.1. ASPECTOS GENERALES

5.2. CONCEPTO GENERAL
5.3. PROPAGACIÓN DE LAS ONDAS
5.4. ECUACIONES DE LINEA DE TRANSMISIÓN
5.5. SOLUCIÓN DE LAS ECUACIONES DE LÍNEA.
 5.6. RELACIÓN DE ONDAS TENSIÓN Y CORRIENTE
5.7. ONDAS DIRECTAS Y REVERSAS
5.8. INTERACCIÓN DE ONDAS VIAJERAS
5.9. REFLEXIÓN Y TRANSMISIÓN DE ONDAS
5.10. REFLEXIÓN Y REFRACCIÓN DE ONDAS

5.10.1. EJEMPLO 1
5.10.2. EJEMPLO 2

5.11. EFECTOS DE LAS TERMINACIONES EN LÍNEAS

5.11.1. TERMINAL ABIERTO

5.11.2. TERMINAL EN CORTO CIRCUITO
5.11.3. TERMINAL CON RESISTENCIA
5.11.4. TERMINAL CON CAPACITANCIA
5.11.5. TERMINAL CON REACTOR

5.12. RESUMEN EN DOMINIO DEL ESPACIO Y TIEMPO

5.13. ATENUACIÓN Y DISTORSIÓN DE ONDAS
5.14. DIAGRAMAS DE LATTICE

6. TEORÍA ELECTROMAGNÉTICA DE LOS EFECTOS TRASITÓRIOS

7. PRINCIPIOS DE MODELACIÓN

7.1. INTRODUCCIÓN
7.2. ELEMENTOS LINEARES CONCENTRADOS Y DESACOPLADOS

7.2.1. Resistance R

7.2.1.1. Análisis del Error

7.2.1.2. Ejemplo de Resistencia en Redes

7.2.2. Inductancia Propia “L”

7.2.2.1. Error Análisis

7.2.2.2. Amortiguamiento de la "Oscilación Numérica" con
Resistencia en Paralelo
7.2.2.3. Razones Físicas para Resistencia en Paralelo

7.2.3. Capacitance C

7.2.3.1. Análisis del Error

7.2.3.2. Amortiguamiento de "Oscilaciones Numéricas" con
Resistencia Serie
7.2.3.3. Razones Físicas de la Resistencia Serie
7.2.3.4. Example for Network with Capacitances
 7.2.4. Series Connection of R, L, C

8. INTERRUPTORES

8.1. Conceptos Básicos

8.2. PARTES
8.3. ACCESORIOS
8.4. PARÁMETROS DE LOS INTERRUPTORES
8.5. INTERRUPCIÓN DEL ARCO
8.6. PRUEBAS DE INTERRUPTORES
8.7. RESISTENCIA AMORTIGUAMIENTO DE TVR
8.7.1. RESISTENCIA AMORTIGUAMIENTO DE TVR

9. MANIOBRAS DE BANCO DE CONDENSADORES

9.1. ENERGIZACIÓN DE BANCO DE CONDENSADORES

9.2. LIMITAR INRUSH DE BCO.CONDENSADORES
9.3. ENERGIZACIÓN DE VARIOS BCOS. DE CONDENSADORES
9.4. DESCONEXIÓN DE BANCOS DE CONDENSADORES
9.5. EJEMPLOS DE APLICACIÓN

10. EFECTOS ARMÓNICOS

 TRANSITORIOS
ELECTROMAGNÉTICOS
1. NOCIONES FUNDAMENTALES.

1.1. ¿PORQUÉ ESTUDIAR LOS TRANSITÓRIOS?

¾ El crecimiento de la demanda eléctrica exige mayor confiabilidad y

continuidad de los suministros;

¾ Las implicancias económicas y ambientales que restringen las nuevas

construcciones del Sector Eléctrico obligan a un máximo aprovechamiento de
las instalaciones existentes, por lo que pueden dar origen a problemas
transitorios.

¾ El crecimiento poco planificado del sistema Eléctrico Peruano esta

originando la aparición de fenómenos dinámicos y transitorios que afectan la
operación. Siendo, el condicionante para la aparición de Transitorios de alta
frecuencia los niveles elevados de tensión de operación;

¾ Por ello, el previo conocimiento y dominio de los fenómenos transitorios,

nos ayudarán a prevenir y dar solución a los mismos.

1.2. AMBITO DE LOS FENÓMENOS TRANSITÓRIOS

1h 10 1min 10 1s 100 10 1ms 100 10 1us 100 10 1

0.1 1 10 100Hz 1 10 100kHz 1 10 100MHz

Acción de los Controles de las Transitórios Electromagnéticos

Turbinas

Transitórios
Electromecánicos
Transitórios
Maniobras

Aplicación de Fallas
Descargas
Oscilación Torsional
Atmosféricas

Control Carga- Transitórios

Frecuencia Rápidos

Resonancia Subsíncrona RSS

Efecto Corona
Estabilidad
Transitória
 TRANSITÓRIOS
SOBRECORRIENTE

SOBRETENSIONES

Transitórias Sostenidas
A frecuencia fundamental ó
Altas frecuencias y
múltiplos y no son amortiguadas
rápidamente amortiguadas

Temporárias
A frecuencia fundamental ó múltiplos y un poco
Amortiguadas.

2. DEFINICIONES

• SOBRETENSIONES
Súbitos incrementos de tensión a altas frecuencias ó frecuencia industrial
originados por impactos en el sistema.
• SOBRECORRIENTES
Gran crecimiento de corriente que resulta de fallas en el sistema y que ayuda a la
determinación de los esfuerzos mecánicos y térmicos dentro de un equipo.
• ONDAS ANORMALES
Distorsiones que aparecen en el sistema durante ciertas condiciones operativas.

2.1. PARAMETROS BÁSICOS

diL
vL = L
vR = RiR dt

R iR =
vR
R L
iL =
1 t
L∫0
vL dt + iL (0)

1 2
EnergíaL = LiL
2
 1 t
C ∫0
vC = iCdt + vC (0)

C iC = C
dvC
dt

1
EnergíaC = CvC2
2

2.2. COMPORTAMIENTO FÍSICO

Elementos almacenadores de Energía

∆ iL
vL = L
∆t
L vL finito→∆iL = 0 en ∆t = 0
iL(0−) = iL(0) = iL(0+)
iL vL

El estado de Energía no se puede alterar instantáneamente

∆ vC
iC = C
∆t

C iL finito → ∆vC = 0 en ∆t = 0
vC (0 − ) = vC (0) = vC (0 + )

iC vC
2.2.1. SOLUCIÓN MATEMATICA
R
dV1
V = R.C * + V1
dt

V C V1
−t
V1 = V − A.e R .C

Resolviendo: Tensión
V = R * I + V1
V
1
C∫
V1 = I .dt

dV1 V1(0)
I =C*
dt

t Tiempo

2.2.1.1. RESPUESTAS EXPONENCIALES

Combinaciones “R-L”

t
eτ, τ = RL
−

ó
τ = cons tan te det iempo

Combinaciones “R-C”
t
eτ, τ = RC
−

τ = cons tan te det iempo

ó
 2.2.1.2. RESPUESTAS OSCILATORIAS

Combinaciones “L-C”

senω0 .t ó cos ω0 .t

1
donde : ω0 =
ó LC

ω0 = frecuencia natural de oscilación

t t
e τ .senω .t e τ . cos ω .t
− −

0 ó 0

ó 1
donde : ω0 =
LC

2.2.2. INTERPRETACIÓN DE RESPUESTAS TRANSITÓRIAS

ω = Frecuenciafundamental

ω0 = Frecuencia no − amortiguada
3. MÉTODOS DE TRANSFORMADA DE LAPLACE PARA SOLUCIÓN DE
ECUACIONES DIFERENCIALES

3.1. MÉTODO OPERACIONAL

Modelación del
problema en el domínio

Obtención de
Condiciones la solución
Iniciales operacional

Respuesta (tensión ó
corriente) en el domínio del
tiempo

TRANSFORMADAS DE LAPLACE

d2y
k2
dt 2
+ k1
dy
dt
+ k0 = f (t ) [ ]
k 2 S 2Y ( S ) − Sy (0) − y ' (0) + k1 [SY ( S ) − y (0)] +
k0
S
= F (S )
 Aplicación
k súbita de
y(t ) = k Y (S ) = corriente
S
continua

k Rampa o
Y (S ) =
y(t ) = kt S2 estiramiento
s pequeños
y(t ) = k1t + k2 k k de un
Y ( S ) = 12 + 2 senoide
S S

−t
k Señales de
Y (S ) =
y(t ) = ke τ S +1 τ ocasion

ω
Y (S ) = k
S2 +ω2 Fuente de
y(t ) = k senωt alimentación

S
Y (S ) = k
S2 +ω2 Fuente de
y(t ) = k cosωt alimentación

ω
Y ( S ) = k1
S2 +ω2
Fuente de
y(t ) = k sen(ωt + θ )
S
+ k2 2 2 alimentación
S +ω
 3.2. SOLUCIÓN DE CIRCUITOS “R-L”

Ecuaciones el el Tiempo: V = Ri + L
di
dt
Transformada de Laplace:
t=0

R V
= RI ( S ) + L[ SI ( S ) − i (0)]
S
V i (t ) = ?

v L (t ) = ?
Condiciones Iniciales:
L i(0) = iL (0) = 0

V
Solución Operacional: I (S ) =
S ( SL + R )

Conceptos:
Transformada Inversa de Laplace:

L → em t = 0 →
L ⎡ t ⎤ t
V ⎢ −
τ ⎥; −
i (t ) = 1− e di V τ
L R⎢ ⎥ v (t ) = L = e
L → em t = ∞ → ⎣⎢ ⎦⎥ L dt R

L
donde : τ = ;
R

SOLUCIÓN DE CIRCUITOS “R-L”

V ⎡ −
τ
t ⎤
i (t ) = ⎢1 − e ⎥;
R ⎢⎣ ⎥⎦
 3.3. SOLUCIÓN DE CIRCUITOS “R-C”
Ecuaciones el el Tiempo:

1
∫
t
V = Ri + 0
idt + vC ( 0 )
C
t =0
Derivando: di 1
R + i = 0
dt RC
V i (t ) = ? Transformada de Laplace:
vC (t ) = ? 1
C [ SI ( S ) − i (0)] + I (S ) = 0
RC

Condiciones Iniciais: V − V0
i ( 0) =
R
Solución Operacional:

Conceptos: ⎛ V − V0 ⎞ ⎛ 1 ⎞
I (S ) = ⎜ ⎟ ⎜S + ⎟
⎝ R ⎠ ⎝ RC ⎠
C → em t = 0 → C
Transformada Inversa de Laplace:
C → em t = ∞ → C V − V0 −τt V −τt
i (t ) = e ; vC (t ) = (V0 − V ) e +V
R R

donde : τ = RC;

SOLUCIÓN DE CIRCUITOS “R-C”

V − V0 −τt
i (t ) = e ;
R
 3.4. SOLUCIÓN DE CIRCUITOS “L-C”

Ecuaciones el el Tiempo:
t=0 vL (t ) = ? di 1 t
dt C ∫0
L V =L + idt + vC (0)

V i(t ) = ? d 2i 1
vC (t ) = ?
Derivando: +
dt 2 LC
i=0
C
Transformada de Laplace:
1
[S2I (S) − Si(0) − i' (0)]+ I (S) = 0
LC
V − V0
Condiciones Iniciales: i (0) = 0 e i ' (0) =
L
V − V0 ω0 1
Solución Operacional: I (S ) = siendo ω0 =
L S 2 + ω0
2
LC

Transformada Inversa de Laplace:

V − V0 1 1 1
i (t ) = sen t vL (t ) = (V − V0 ) cos t vC (t ) = V − (V − V0 ) cos t
LC LC LC LC

SOLUCIÓN DE CIRCUITOS “L-C”

1
vC (t ) = V − (V − V0 ) cos t
LC
 3.5. ESTUDIO DE LOS CIRCUITOS “R – L – C”

CIRCUITO SERIE: CIRCUITO PARALELO:

L
R L C
fuentes de alimentación:
C
- tensión contínua
- tensión rampa
- tensión alterna

Estos circuitos se consideran importantes, porque las redes eléctricas reales en su

configuración contienen éstos componentes que envuelven los problemas
transitórios.

SERIE
PARALELO

d 2 Ψ R dΨ Ψ d 2Φ 1 dΦ Φ
+( ) + = F(t) 2
+( ) + = F(t)
dt 2 L dt LC dt RC dt LC

Ψ = Voltaje de un componente . Φ = Corriente de un ramal

donde : Tp = R.C donde : TS =

L
R
Si multiplicaramos las dos constantes de tiempo tenemos:

TP .TS = L.C = T 2
Sea : Z 0 = ( L / C )1/ 2 impedancia de onda

1/2
R ⎛C ⎞
Definiendo " η " como : η = = R *⎜ ⎟
Z ⎝ L ⎠
o

TP R 2 .C
=η
2
Luego : =
TS L
De la solución común establecida con la relación anterior, calcularemos el
comportamiento de la corriente por el ramal del inductor para el caso particular del
circuito paralelo.

Y del cual obtendremos las siguientes gráficas.

IL(t)

0 t

η = 1 / 2 : Criticamente amortiguamiento

IL(t)

0 t

η 〈 1 / 2 : Sobre amortiguamiento
IL(t)

SOLUCIÓN GENERALIZADA

Del ejemplo anterior se puede plantear una solución generalizada que tendrá la
siguiente solución:

de esta ecuación separamos la parte variable en una función de “η” y la

graficaremos: V (0) 2ηε −t '/ 2η t'
I L (t ) = *T * * sen(4η 2 −1)1/ 2
L (4η −1)
2 1/ 2
2η
2η t'
f1 (η ) = * ε −t '/ 2η * sen(4η 2 − 1)1/ 2
(4η − 1)
2 1/ 2
2η
 FAMILIA DE SOLUCIONES

4. SIMPLES TRANSITÓRIOS DE CONEXIÓN.

4.1. INTERPRETACION FISICA

Interpretación Física:
Red Eléctrica
P+jQ

Circuito Equivalente:
PROPOSITO:
Investigar los efectos de la tensión en los terminales de la carga y entre los contactos del
interruptor.
Por ejemplo, en los rechazos de carga
 4.2. SOBRETENSIONES EN LA CARGA

Para el cálculo de la tensión transitória en

v( t ) terminales de la carga “R-L-C”, debe tomarse
en cuenta que, al interrumpirse la corriente
i( t ) i(0)= 0, la tensión estará passando por un valor
v(0) = Vo.
En estas condiciones, utilizando la
metodología anterior se obtiene para la tensión
entre los terminales de la carga la expresión η
V0
>½.
Luego de esta maniobra, se presentan
sobretensiones transitórias (primeros ciclos) y
sostenidas (luego de estabilizarse).

_
v(t)
v(t) =
V0
⎧
−t '

[
v(t ) = e ⎨cos (4η − 1) .(t 2η ) −
_
2η 2 12 '
] [
sen (4η 2 − 1) .(t ' 2η ) ⎫
12

⎬
]
⎩ (4η 2 − 1)1 2 ⎭

4.3. RESISTENCIA DE ATENUACIÓN

L Para un sistema de
345kV,25000MVA.
Capacidad del INT. 40Ka

R 345 2
VT = 2 * = 562kV
3
C La magnitud de corriente
inyectada que interesa se da
en el primer medio ciclo de
inicio de la falla esto es porque la frecuencia natural del circuito es usualmente mucho
mayorque el de la red. Por eso la corriente que interesa se considera una rampa:

V
I= *t
L
Cuando el sist. Tiene neutro a tierra el sobrevoltaje se reduce al 70% osea 393.4kV que
según las graficas anteriores corresponde a η=1,8.

X L = 345 / 40000 3 = 5Ω
L = 13.2mH

CÁLCULO RESISTENCIA DE ATENUACIÓN

Para un nivel de tensión de 345kV, tiene una capacitancia aproximada de 25000 pF.
Luego:

1/ 2 1/ 2
⎛L⎞ ⎛ 1.32 *10−2 ⎞
Z0 = ⎜ ⎟ =⎜ −8
⎟ = 725 Ω
⎝C ⎠ ⎝ 2.5 *10 ⎠

Luego, para :η = 1,8

R = η * Z 0 = 1,8 * 725 = 1300 Ω

4.4. SOBREVOLTAJE TRANSITÓRIO

4.4.1. INTERRUPCIÓN DE CORRIENTE ASIMÉTRICA

4.4.2.RECUPERACIÓN DE LA TENSIÓN DE ARCO EN EL INTERRUPTOR
LUEGO DE APERTURA

4.4.3. SWITCHING DE CONDENSADORES

 4.4.4. TRANSITÓRIO DE DOBLE FRECUENCIA

5. ONDAS NAVEGANTES Y OTROS TRANSITÓRIOS EN LÍNEAS DE

TRANSMISIÓN

A B

A B

…..
 5.1. ASPECTOS GENERALES

Una onda eletromagnética se propaga i v

en un medio (línea, cable, etc…)
através de una forma tal que:
Z0
A B

Siendo: v
= Z0
i

Donde: L Impedancia característica de la línea

Z0 =
C
Para líneas: µ0 ⎛ D ⎞ πξ0
L= ln⎜ ⎟ [ H / m] C= [ F / m]
π ⎝r⎠ ; ⎛D⎞
ln⎜ ⎟
⎝r⎠

kg .m Hr Coulomb2 .seg2 F
µ O = 4π . 10 − 7 2
= 4π . 10 − 7 ε O = 8.854x10−12 = 8.854x10−12
3
Coulumb m kg.m m

250 hasta 400 Ω - lineas aéreas

Z0
Por tanto:
40 hasta 80 Ω - Cables Aislados
50 hasta 100 Ω - Máquinas Rotativas

5.2. CONCEPTO GENERAL

Matemáticamente, una onda viajera puede representarse simbólicamente por: V(x,t) =

F(t-x/vt)
Donde, V(x,t) es la tensión en función del tiempo y su posición en el conductor, definida
por una distancia “x” a partir de un punto de referencia.
Esta onda viajera se desplaza en el conductor con una velocidad de propagación igual a:
Generalizando, para cualquier posición a lo largo del conductor X=Xn, una onda definida
como V(Xn,t) = F(t-Xn/Xt) no aparece en ese punto hasta después de transcurrido un
tiempo t=Xn/vt después de una inicio de la onda en x=0.
 5.3. PROPAGACIÓN DE LAS ONDAS

1° Para líneas aéreas conduce a:

• Velocidad de propagación = 300.000 km/s

2° Para cables aislados µ cable ≈ µ0 e ξcable ≈ 3 a 5 veces ξ0
Así adoptando un valor medio para : ξcable ≈ 4ξ0

• Velocidad de propagación = 150.000 km/s

3° Para Máquinas Rotativas, se desplaza por las bobinas según :

• Velocidad de propagación = 30,000 km/s

5.4. ECUACIONES DE LINEA DE TRANSMISIÓN

Circuito básico de una Línea:

Expandiendo con serie de Taylor y reacomodando las ecuaciones se llega a la siguiente

relación simplificada:
Luego, desarrollando en ecuaciones de segundo orden, se expresiones que describen el
comportamiento de la Tensión y corriente a lo largo de la línea.

Considerando una línea ideal sin pérdidas y haciendo que R = G = 0, se tienen la

ecuación simplificada siguiente

5.5. SOLUCIÓN DE LAS ECUACIONES DE LÍNEA.

Resolviendo con la transformada de Laplace las ecuaciones diferenciales anteriores se

llega a la siguiente solución:

∂ 2V ( x, s ) donde : V ( x , s ) es la transforma da
= s 2 LC .V ( x, s ) de Laplace de V ( x , t )
∂x 2

Donde, la parte exponencial es función de la unidad de longitud y las constantes A(s) y

B(s) son función del tiempo. Por propiedad de Laplace se tiene que:

V ( x, s) = A( s ).e − s LC . x
+ B( s ).e s LC . x
− − − − − −(5)

I ( x, s) = D( s).e − s LC . x
+ E ( s).e s LC . x
− − − − − −(6)

Aplicando esta propiedad a las ecuaciones (5) y (6), se tiene:

Estos resultados pueden ser interpretados como dos ondas viajeras de tensión ó corriente
que viajan a lo largo del conductor, una en dirección positiva “+x” y la otra en la dirección
reversa “-x” desplazándose con la velocidad :

V ( x, t ) = A(t − LC .x).u (t − LC .x) + B (t + LC .x).u (t + LC .x)

I ( x, t ) = D (t − LC .x).u (t − LC .x) + E (t + LC .x).u (t + LC .x)
Estos resultados pueden ser interpretados como dos ondas viajeras de tensión ó corriente
que viajan a lo largo del conductor, una en dirección positiva “+x” y la otra en la dirección
reversa “-x” desplazándose con la velocidad:

vt = 1 / LC

5.6. RELACIÓN DE ONDAS TENSIÓN Y CORRIENTE

En las ecuaciones (1) y (2) desarrolladas, aplicaremos las transformadas de Laplace y

obtendremos:
∂V ( x, s ) ∂I ( x, s )
− = sL. I ( x, s ) − = sC. V ( x, s )
∂x ∂x
Luego, aplicando la solución de las ecuaciones de línea se llegan a las siguientes
relaciones:

V ( x, s) = s LC A( s).e − s LC . x
− s LC B( s).e s LC . x

I ( x, s) = sL D( s ).e − s LC . x
+ sL E ( s ).e s LC . x

Igualando estos coeficientes, a los de las ecuaciones (5) y (6) llegamos a lo siguiente:

V ( x, s ) = A( s ).e − s LC . x
− B( s ).e s LC . x
V ( x, t ) = F1 ( x − vt .t ) + F2 ( x + vt .t )

A( s ) − s B(s ) s
I ( x, s ) = .e LC . x
+ .e LC . x
1
Zo Zo I ( x, t ) = [F1 ( x − vt .t ) + F2 ( x + vt .t )]
Zo
1 1 1
donde : D(s) = A( s ) ; E (s) = B( s) ; Z O =
LC LC LC

5.7. ONDAS DIRECTAS Y REVERSAS

• Ondas Directas
(polaridades de las tensiones y corrientes coincidentes):
• Ondas Reversas
(polaridades de las tensiones y corrientes contrárias):

5.8. INTERACCIÓN DE ONDAS VIAJERAS

• Ondas de Tensión positivas y opuestas

• Ondas de Tensión superpuestas

• Ondas de Tensión continuando su propagación

5.9. REFLEXIÓN Y TRANSMISIÓN DE ONDAS

Sean dos conductores empalmados, de diferentes características por los que viajan las
ondas de tensión y corriente.
La relación de la onda incidente::
V1
I1 =
ZA
Y las relaciones de las ondas reflejadas y transmitidas son:
V2 V3
I2 = − I3 =
ZA ZB

Además:
V1 + V2 = V3 I1 + I 2 = I 3

Se puede establecer:
V1 V2 V3
− =
Z A Z A ZB
Con estas ecuaciones se puede calcular las siguientes expresiones:

⎛ Z + ZA ⎞ ⎛ Z − ZA ⎞
V2 ⎜⎜ B ⎟⎟ = V1 ⎜⎜ B ⎟⎟
⎝ Z B .Z A ⎠ ⎝ Z A .Z B ⎠

⎛ Z − ZA ⎞
V2 = ⎜⎜ B ⎟⎟.V1 − −(7)
⎝ B
Z + Z A ⎠

⎛ Z + Z A ⎞ 2V1
V3 ⎜⎜ B ⎟⎟ =
⎝ Z B .Z A ⎠ Z A
⎛ 2Z B ⎞
V3 = ⎜⎜ ⎟⎟.V1 − −(8)
⎝ ZB + Z A ⎠
La potencia transferida por la línea:
V12
I1V1 = watts − −(9)
ZA

5.10. REFLEXIÓN Y REFRACCIÓN DE ONDAS

De las ecuaciones (7) y (8) se obtienen:

 ⎛ Z − ZA ⎞ ⎛ 2Z B ⎞
Coef . Re flexión = ⎜⎜ B ⎟⎟ Coef .Transmisión = ⎜⎜ ⎟⎟
⎝ ZB + Z A ⎠ ⎝ ZB + Z A ⎠

Efectos de los puntos de unión de líneas de Transmisión en Tensión y corriente

5.10.1. EJEMPLO 1

• Z2 > Z1 - onda propagándose de A hacia B:

kT
Z1 = 60Ω kR Z 2 = 400Ω
B C
A
2 × 400
kT = = 1,74
400 + 60

400 − 60
kR = = 0,74
400 + 60

• La onda transmitida es de 174% de la

onda incidente;

• La onda reflejadada es de 74% de la

onda incidente y se suma a esta;
 5.10.2. EJEMPLO 2

• Z2 < Z1 - onda propagándose de A hacia B:

kT
Z1 = 300Ω kR Z 2 = 50Ω
B C
A
2 × 50
kT = = 0,286
300 + 50

50 − 300
kR = = −0,714
300 + 50

• La onda transmitida es de 28,6%

de la onda incidente;

• La onda reflejada es negativa e

igual a 71,4% de la onda incidente.
Esto resulta en una onda total de
28,6%;

5.11. EFECTOS DE LAS TERMINACIONES EN LÍNEAS

Una situación típica frecuentemente encontrada es:

v1
v1 v3
Donde:
v1 - onda incidente
v2 - onda reflejada Z1
v2 B
v3 - onda transmitida A Z2
Siendo: La Máxima transferencia se produce
cuando: V3 = V1
•v2 = kR.v1
Coef .Transmisió n : k T = 1
• v3 = kT.v1 ⎛ 2Z 2 ⎞
1 = ⎜⎜ ⎟⎟.
•k T = 2Z2 /(Z1+Z2) ⎝ Z 2 + Z1 ⎠

•k R = (Z2 - Z1) /(Z1+Z2) Z 2 + Z1 = 2 Z 2

Z1 = Z 2
Osea cuando las impedancias de la línea y la
carga son iguales

5.11.1. TERMINAL ABIERTO

ANTES
Z2 → ∝ (circuito abierto)-onda Circuito
V abierto
propagándose de A hacia B.

A B X
Z1 I
V Z2 → ∞

DESPUÉS
2V
kT = 2 kR = 1
V

V
V
De los índices se deduce que la onda
reflejada es el 100% de la onda incidente y X
se suma a ella duplicando su tensión final, I
mientras que la corriente tendrá valor I
negativo que se resta a la incidente
resultando finalmente una corriente cero.
-I

5.11.2. TERMINAL EN CORTO CIRCUITO

Z2 = 0 (corto circuito)
onda propagándose de A hacia B.
 A B

Z1
Z2 = 0
V

kT = 0 k R = −1

De los índices se deduce que la onda reflejada es el 100% negativa de la onda

incidente y se resta a ella dando una tensión final de cero, mientras que la corriente
tendrá valor positivo y se suma a la incidente resultando finalmente una corriente
doble.

ANTES Corto
V Circuito

2V DESPUÉS

V
V

-V
I

5.11.3. TERMINAL CON RESISTENCIA

Caso Particular, cuando: R = Z1.

A B

Z1
V R

kT = 1 kR = 0

Para la condición R=Z1, toda la tensión y corriente se transmiten a la carga. En

forma general, para otros valores de “R” se cumplirá que:

⎛ R − Z1 ⎞ ⎛ 2R ⎞
Coef . Re flexión = ⎜⎜ ⎟⎟ Coef .Transmisión = ⎜⎜ ⎟⎟
⎝ R + Z 1 ⎠ ⎝ R + Z1 ⎠

2V B
ANTES

V
V

I
I

DESPUÉS
V

V V

I1 I3

5.11.4. TERMINAL CON CAPACITANCIA

 Zc → -j 1/ωC (circuito capacitivo)

A B
Z1
V C1

V1 ⎛ (1 / C1s ) − Z1 ⎞ V2 ( t ) = V1 (1 − 2.e −α .t )
V2 ( s ) = ⎜ ⎟→
s ⎜⎝ (1 / C1s ) + Z1 ⎟⎠

V1 ⎛ 2 / C1s ⎞
V3 ( s ) = ⎜ ⎟→ V3( t ) = V1.(2 − 2.e −α .t )
s ⎜⎝ (1 / C1s) + Z1 ⎟⎠

1
Siendo : α=
C.Z1
La onda de tensión se incrementará, mientras que la onda de corriente tenserá a cero,
después que se cargue el condensador.

5.11.5. TERMINAL CON REACTOR

Z2 = jωL → circuito inductivo.

Z1
L
V

V1 ⎛ ( L .s ) − Z 1 ⎞
⎟⎟ → V2 ( t ) = −V1 (1 − 2.e )
− β .t
V2 (s) = ⎜⎜
s ⎝ ( Ls ) + Z 1 ⎠
V ⎛ 2 L .s ⎞
V (s) = 1 ⎜ ⎟ → V3( t ) = 2.V1.e − β .t
3 s ⎜ ( Ls ) + Z ⎟
⎝ 1 ⎠
Z1
siendo : β=
L

Inicialmente, la tensión se incrementa transitoriamente tendiendo a estabilizarse luego,

mientras que la corriente sufre una disminución instantánea y luego se incrementa.
 5.12. RESUMEN EN DOMINIO DEL ESPACIO Y TIEMPO

5.13. ATENUACIÓN Y DISTORSIÓN DE ONDAS

En una línea real, existen pérdidas longitudinales y transversales, que originan que las
ondas viajeras sufran pérdidas de amplitud y distorsiones en frente de onda.
Estos procesos se pueden expresar por las siguientes relaciones exponenciales:

⎛ G x ⎞ ⎛ G ⎞ ⎛ R x ⎞ ⎛ R ⎞
⎜− ( )⎟ ⎜ − ( ).t ⎟ ⎜− ( )⎟ ⎜ − ( ).t ⎟
V = VO e ⎝ C v ⎠
= VO e ⎝ C ⎠ I = I O e⎝ L v ⎠
= I O e⎝ L ⎠

1/ 2
V ⎛L⎞
La relación tensión corriente, se expresa por: = ZO = ⎜ ⎟
I ⎝C ⎠
⎡ G R⎤
VO ⎢⎣ −( C ) + L ⎥⎦ .t
.e = ZO
IO
⎡
Para un caso particular, en que: G R⎤
⎢
⎢ − ( ) + ⎥⎥ = 0
⎢⎣ C L ⎥⎦

R G
Se obtiene que: =
L C
ATENUACIÓN Y DISTORSIÓN DE ONDAS

La reducción de la amplitud
máxima de la onda puede ser 1 ⎛ R G ⎞
G β= ⎜
2 ⎝ L
−
C⎠
obtenida de la ecuación:
V = VO .e −α .t
β = Coeficiente
donde: R
Distorsión
1⎛R G⎞
α= ⎜ + ⎟
2⎝ L C ⎠
1 2 3 4

Onda incidente de Distorsión por

Frente escarpado Efecto Corona

Atenuación del Atenuación de la

Frente por Pérdidas Amplitud por Pérdidas
Longitudinales Transversales
α = Coeficiente de Ped L = I 2 .R PedT = V 2 .G

atenuación

5.14. DIAGRAMAS DE LATTICE

Antes de aparecer los sistemas computarizados en el mercado, Bewley propuso un método

de representación Gráfica para realizar los estudios de sobretensiones, ésta era una
adecuación para el estudio de ondas viajeras através de diagramas espacio-tiempo.
Para apreciar la aplicación de esta metodología, mostraremos a continuación un sistema
simple en el cual la propagación de ondas viajeras se inicia con el cierre del interruptor
“S”. Éstas se propagan hasta el terminal resistivo con una impedancia distinta,
prodiciéndose una transmisón y una reflexión la que retorna a la fuente, repitiéndose el
proceso.
 Z
S T=
R
V

0
V1
T V’1
2T
V2
3T
V’2
t
4T
V3
5T V’3
6T

Sean las siguientes ecuaciones formuladas a partir del diagrama anterior:

V1 = Z .I 1 (onda incidente)

En el punto de conexión de la línea con la carga resistiva se tendrá:

VT = R.I T

VR = α1.V1 (tensión reflejada )

VT = V1' = β1.V1 (tensión transmitida )

R−Z
Coef . de Re flexión = α1 =
R+Z
2R
Coef . de Transmisión = β1 =
R+Z
Luego, al retornar la onda al punto de la fuente, se tendrá:

Z FUENTE − Z
Coef . de Re flexión = α 2 = = −1
Z FUENTE + Z
2 Z FUENTE
Coef . de Transmisión = β 2 = =0
Z FUENTE + Z

Los coeficientes de reflexión y transmisión se relacionan así:

β1 = 1 + α1
β2 = 1+ α2
Con éstas relaciones establecidas y aplicadas al circuito simple se podrá trazar el diagrama
de Lattice, teniendo en cuenta que el tiempo de tránsito de las ondas en la línea esta dado
por la relación:

T =(L.C)1/2

Donde: L y C son la inductacia y capacitancia de toda la línea.

DIAGRAMAS DE LATTICE

Sean las siguientes relaciones:

V1 = V = Onda incidente en el terminal remoto.

V’1 = α1.V1 = α1.V = Reflexión de V1 en el terminal de la carga.
α1 = Coeficiente de reflexión en el terminal de la carga.
V2 = α2.V’1 = α1. α2.V1 = Reflexión de V’1 en el terminal de la fuente.
α2 = Coeficiente de reflexión en el terminal de la fuente.
V’2 = α1.V2 = α21. α2.V1 = Reflexión de V2 en el terminal de la carga.
En cualquier punto la tensión de la onda se obtiene de la suma de las ondas incidentes y
reflejadas que pasan por ese punto. En el terminal de la carga, con relación a la tensión se
tiene:

t<T Tensión = 0
T < t < 3T Tensión = V+ α1 .V
3T < t < 5T Tensión = β1.(1+ α1. α2 ).V
5T < t < 7T Tensión = β1.(1+ α1. α2+ α21. α22 ).V
Corto Circuito

Z
S T = (L.C)1/2
R=0
V

En el terminal de la carga R=0, y se

tiene que:
α1 = - 1
V1 = - V1 = - V
Esta onda de tensión estará asociada
a una onda de corriente de sentido
opuesto, propagándose en dirección
reversa (negativo).
Así se tiene, que la onda de corriente
reflejada la primera vez, duplica el
valor de la onda incidente

6. TEORÍA ELECTROMAGNÉTICA DE LOS EFECTOS TRASITÓRIOS

6.1.

Los Fenómenos de transitorios electromagnéticos son un conjunto de interacciones de

efectos electrostáticos y electromagnéticos, de cambios de estados e intercambios de
energía.

Es por ello que resulta importante tener una sólida base con los conceptos fundamentales
de la teoría Electrostática, así como la formación y efectos de los campos
electromagnéticos.

Por ejemplo, una línea muy larga en 345 kV al encontrase operando, almacena una energía
electrostática transversal, mediante el condensador que se forma por conductor tensionado
y la tierra.
Un efecto estático aparece en una línea desenergizada debido a la fricción del aire con la
superficie del conductor que le da una carga respecto del suelo.

El efecto electrostático se pone de manifiesto cuando se induce tensión en una línea

telefónica que pasa cercana a las líneas de alta tensión.

Los efectos electromagnéticos que aparecen en los elementos conversores de energía,

como son los generadores y motores eléctricos.

El Campo Eléctrico aparece entre dos cuerpos cargados con diferentes signos. La
diferencia de potencial que aparece es proporcional a la carga.

Q = C.V ó V= Q/C

El Campo Eléctrico aparece entre dos cuerpos cargados con diferentes signos. La
diferencia de potencial que aparece es proporcional a la carga.

Q = C.V ó V= Q/C

ATENUACIÓN Y DISTORSIÓN DE ONDAS

En una línea real, existen pérdidas longitudinales y transversales, que originan que las
ondas viajeras sufran pérdidas de amplitud y distorsiones en frente de onda.

7. PRINCIPIOS DE MODELACIÓN

7.1. INTRODUCCIÓN

En este resumen se discute los grandes métodos matemáticos de solución que son
empleados en los Programas de Transitorios Electromagnéticos para la simulación digital de
los fenómenos de transitorios. En el desarrollo de los métodos, se adecua la solución que es
de propósito general del programa, y éstos son discutidos a continuación.
 Para analizar problemas específicos, se podrán comparar otros métodos que son
competitivos; por ejemplo, el método de la transformada de Fourier puede ser preferible
para estudios de atenuación y distorsión de ondas a lo largo de una línea en casos donde el
tiempo de estudio es corto, pero no es lo mismo cuando la onda viene de un lugar remoto
hasta el final de la línea.
El programa mediante cálculos computacionales, no pueden ser simulados fenómenos
continuos, solamente discretos en intervalos de tiempo de paso de iteración “∆t”. Esto
genera pequeños errores de truncación que se van acumulando paso a paso y pueden causar
divergencias con la solución verdadera. Los métodos usados son numéricamente estables y
aborda éstos tipos de error de diseño de los algoritmos.

Los programas diseñados para Transitorios Electromagnéticos pueden resolver algunas

redes que consisten de interconexiones de resistencias, inductancias, capacitancias, de
simple y multifases de circuitos, líneas con parámetros distribuidos y otros elementos. En
la Figura 1.1 se muestran detalles de una gran red en una región alrededor del nodo 1.
Suponiendo que las tensiones y corrientes tendrán que ser calculadas instantáneamente en
los instantes 0, ∆t, 2∆t, etc., siendo el paso de iteración “t-∆t” y la solución más cercana en
“t”. En algún instante de tiempo, la suma de las corrientes que fluirá desde el nodo 1 hasta
cada uno de los ramales, será igual a la suma de las corrientes inyectadas, como:

i 12 (t) + i 13 (t) + i 14 (t) + i 15 (t) = i 1 (t) (1.1)

Fig. 1.1 – Detalles de una gran red alrededor del nodo N°1.

Los nodos de voltaje son usados como variables de estado en el Programa, éstas son
necesariamente expresados en ramales de corrientes, i12, etc., como funciones de los nodos
de voltaje.

Para el caso de una resistencia tenemos la siguiente relación:

1 (1.2)
i12 (t) = {v 1 (t) - v 2 (t)}
R
Para una inductancia, una simple relación es obtenida para reemplazar la ecuación
diferencial:
di
v=L
dt
con la siguiente ecuación:
v(t) + v(t - ∆t) i(t) - i(t - ∆t)
=L
2 ∆t
Esta ecuación puede ser escrita, para el caso de la Fig. 1.1, como:

∆t (1.3a)
i 13 (t) = { v 1 (t) - v 3 (t)} + hist 13 (t - ∆t)
2L

con el término “hist13“ conocido desde el valor previo al paso de tiempo.

∆t
hist13 (t - ∆t) = i13 (t - ∆t) + { v1 (t - ∆t) - v 3 (t - ∆t)} (1.3b)
2L

La derivada para la ecuación del ramal capacitivo es análogo, e igual a:

2C
i14 (t) = { v1 (t) - v4 (t)} + hist14 (t - ∆t) (1.4a)
∆t
con el “hist14” siendo conocido desde el valor previo al paso de tiempo
2C
hist14 (t - ∆t) = - i14 (t - ∆t) - { v1 (t - ∆t) - v4 (t - ∆t)} (1.4b)
∆t

Para la línea de transmisión de la Fig. 1.1, entre los nodos 1 y 5, las pérdidas serán ignoradas
en esta solución y las ecuaciones de onda son:
∂v ∂i
- = L′
∂x ∂t

∂i ∂v
- = C′
∂x ∂t

donde:
L', C' = inductancia y capacitancia por unidad de longitud (1),
x = distancia desde el inicio al final,

Tenemos ya conocida la solución debido a d'Alembert:

i = F(x - ct) - f(x + ct)

v = Z.F(x - ct) + Z.f(x + ct) (1.5a)

con:

1
El prima usado en L', C' es para distinguir estos parámetros distribuidos de la ecuación (1.5a); en que la corriente es
multiplicada por Z y sumada al voltaje, de los parámetros concentrados L, C.
F(x - ct) y f(x + ct) = funciones de la expresión compuesta (x – ct) y (x + ct),
Z = impedancia de onda Z = √L'/C' (constante),
c = velocidad de propagación de la onda (constante).

v + Zi = 2ZF(x - ct) (1.5b)

En esta expresión compuesta (v + Z.i) no cambia si (x – ct) tampoco cambia. Por ejemplo,
imaginemos un observador ficticio viajando en línea con la velocidad de onda “c”. La
distancia recorrida por este observador será x = x0 + ct (x0 = localización en punto de
inicio), con x - ct constante.
Sí, x – c.t es constante, el valor de v + Z.i vista por el observador será necesariamente
constante, con un tiempo de viaje: τ = longitud de línea / c,

Por ejemplo para la Fig. 1.1, sea un observador que sale del nodo 5 en un tiempo t - τ verá
el valor de v5(t - τ) + Zi51(t - τ), y al llegar al nodo 1 (luego del tiempo de viaje τ) verá el
valor v1(t) - Zi15(t) (el signo negativo es porque i15 tiene la dirección opuesta de i51); pero
después este valor visto por el observador será constante,
1
i15 (t) = v1 (t) + hist15 (t - τ) (1.6a)
Z

Donde los términos hist15 son conocidos de los valores previos calculados:

1
hist15 (t - τ) = - v5 (t - τ) - i 51 (t - τ) (1.6b)
Z

Otro ejemplo: sea el paso ∆t = 100 µs y τ = 1 ms (tiempo total); en las ecuaciones (1.6)
pueden ser calculados los datos históricos de la línea y almacenados en espacios de tiempo
iguales a “τ”, donde los valores de la ecuación (1.6b) son primeramente calculados en 10
pasos de time. Las ecuaciones (1.6) son soluciones exactas de las pérdidas en la línea si
“∆t” es un integrador de los múltiples “τ”; si no lo es usará la interpolación linear
incurriendo en errores. Las pérdidas pueden ser representadas con suficiente precisión
insertando una resistencia de parámetro concentrado en algunos tramos a lo largo de la línea,
así como también presentar sofisticados métodos con parámetros dependientes de la
frecuencia los que serán detallados más adelante en el Tema de modelación de líneas.

Si las ecuaciones (1.2), (1.3a), (1.4a) and (1.6a) son reemplazadas en la ecuación (1.1), para
el nodo 1 se tendrá:
i1 (t) - hist13 (t - ∆t) - hist14 (t - ∆t) - hist15 (t - τ) (1.7)
Para algún tipo de red con ”n” nodos, un sistema de “n” ecuaciones puede ser formado, con
la siguiente relación:
[G] [v(t)] = [i(t)] - [hist] (2) (1.8a)

con: [G] = matriz simétrica de conductancia nodal de “n x n”,

[v(t)] = vector voltaje de “n” nodos,
[i(t)] = vector corriente de “n” fuentes, y
[hist] = vector de “n” términos históricos conocidos.

En una red, algunos nodos tienen voltajes conocidos debido a que se encuentran conectados
a fuentes de voltaje y otros nodos conectados a tierra; por eso la ecuación anterior puede ser
particionada en nodos “A” con voltajes desconocidos y nodos “B” con voltajes conocidos,
de la cual podemos determinar los voltajes desconocidos [vA(t)] resolviendo la ecuación:

[G AA ] [ v A (t)] = [i A (t)] - [hist A] - [G AB] [ vB (t)] (1.8b)

Haciendo una reseña comparativa de los métodos de solución mencionamos lo siguiente:

ANTES:
En los inicios de simulación de transitorios, los programas fueros diseñados para arrancar el
caso desde la condición inicial de cero. En algunos casos, los términos históricos hist13,
hist14 y hist15 en la ecuación (1.7) son simplemente cero. Pero había en algunos casos, en
que la simulación transitoria se iniciaba desde la condición inicial AC de estado estable a la
frecuencia de 50 ó 60 Hz. Originalmente, las condiciones iniciales de estado estable eran
leídas de las salidas de programas auxiliares que daban esta solución (subrutinas) los cuales
tenían modelada la red, y a partir de estas salidas recién se realizaba la primera iteración
transitoria y esta información se iba almacenando como información histórica, proceso que
era tedioso y que tomaba mas tiempo de simulación.

ACTUALIDAD:
Los métodos computaciones actuales empleados por los Programas de Transitorios proceden
de la siguiente forma:
Se construyen las matrices [GAA] y [GAB], luego [GAA] es triangularizada con eliminación
por ordenación y aprovechando la esparcidad. En cada paso de tiempo de iteración, el
vector de la derecha de la ecuación (1.8b) es "ensamblado" desde los términos históricos
conocidos y de las fuentes de tensión y corriente conocidos. Luego, el sistema de ecuaciones
lineares es resuelto para [vA(t)], usando la información contenida en la matriz de

(2)Los corchetes son usados para indicar que es una matriz y cantidad vectorial.
conductancia triangularizada. En este proceso de solución repetitiva se aprovecha al
máximo la simetría de la matriz. Luego, para el proceso en el próximo paso de iteración
toma los datos históricos “hist” de las ecuaciones (1.3b), (1.4b) y (1.6b) y los actualiza para
usarlos en los siguientes pasos de iteración.

Luego, continuando con la formulación planteada tenemos:

La solución AC de estado estable desarrollada para el caso de la Fig. 1.1, usando las
ecuaciones de nodo de la ecuación (1.1), tendremos lo siguiente:

I12 + I13 + I14 + I15 = I1 (1.9)

Donde la corriente “I” es un fasor de cantidades complejas ahora ⏐I⏐⋅ejα; y para elementos
de parámetros concentrados tenemos las siguientes ecuaciones:

1
I12 = ( V1 - V 2 ) (1.10)
R
para la inductancia,

1
I13 = ( V1 - V 3 ) (1.11)
jωL

y para la capacitancia ,
I14 = jωC ( V1 - V 4 ) (1.12)

Para una línea con parámetros distribuidos R', L', G', C', la solución de estado estable es:
⎡ 1 - Yseries ⎤
⎢ Yseries + 2 Yshunt ⎥
⎡ I15 ⎤ ⎢ ⎥ ⎡ V1 ⎤ (1.13)
⎢ ⎥=⎢ ⎥ ⎢⎢ ⎥
⎢⎣ I 51 ⎥⎦ ⎢ - Yseries 1 ⎥ ⎣ V 5 ⎥
⎦
⎢⎣ Yseries + Yshunt ⎥
2 ⎦

Si la representación del circuito “π” equivalente de la Fig.1.2 es usada con:

Fig. 1.2 Circuito Equivalente “π” para solución

AC de estado estable de líneas de
Transmisión.
 1 s e nh( γ_ )
Yseries = , con Zseries = (R ′ + jωL ′)
Zseries γ_

⎛ γ_ ⎞
tanh⎜ ⎟
1 1 ⎝ 2 ⎠
Yshunt = (G ′ + jωC′) (1.14)
2 2 γ _
2

Agrupando tenemos:

1
Yseries + Yshunt = cosh( γ_ ) • Yseries
2

1
Yseries + Yshunt = cos(ω_ L ′C′ ) • Yseries
2

donde “γ” es la constante de propagación,

γ = (R ′ + jωL ′) (G ′ + jωC′) (1.15)

Para el caso de pequeñas pérdidas se considera: R' = 0 y G' = 0, y en la ecuación (1.14) se

tiene:

s e n (ω_ L ′C′ )
Zseries = _ • jωL ′ •
ω_ L ′C′

⎛ ω_ ⎞
tan ⎜ L ′C′ ⎟
1 _ ⎝ 2 ⎠
Yshunt = • jωC′ • (1.16)
2 2 ω_
L ′C′
2

Donde: si el valor de “ω_” es pequeño, típicamente ≤ 100 km a 60 Hz para líneas cortas, la

relación senh(x) / x y tanh(x/2) / x/2 en la ecuación (1.14) y el sen(x) / x
, tan(x/2) / x/2 en la ecuación (1.16) todos serán considerados igual a 1.0.
De esta simplificación del circuito “π”, se tiene:

Zseries = (R ′ + jωL ′)

1 1
Yshunt = (G ′ + jωC′) (1.17)
2 2

Con el circuito equivalente “ π” de la Fig. 1.2, las ecuaciones de los ramales para pequeñas
pérdidas en la línea finalmente será:
 1
I15 = ( Yseries + Yshunt ) V1 - Yseries V 5 (1.18)
2

Ahora, escribiendo las ecuaciones para el nodo 1, desde las ecuaciones (1.10), (1.11), (1.12)
y (1.18) dentro de la ecuación (1.9), tendremos:

⎛1 1 1 ⎞ 1 1
⎜⎜ + + jωC + Y series + Y shunt ⎟⎟ V 1 - V 2 - V 3 - jωC V 4 - Y series V 5 = I 1 (1.19)
⎝ R jωL 2 ⎠ R jωL

Para algún tipo de red con “n” nodos, puede ser formado un sistema de “n” ecuaciones,
siendo su forma general:
[Y] [V] = [I] (1.20)

con:
[Y] = Matriz de admitancia simétrica nodal, con elementos complejos,
[V] = vector voltaje de “n” nodos (valor del fasor complejo),
[I] = vector de fuente de corriente (valor del fasor complejo).

Donde los términos de la izquierda de la ecuación (1.20) es particionada dentro en una red
con “A” nodos con voltajes desconocidos y en otra de “B” nodos con voltajes conocidos.
Los voltajes desconocidos son hallados al resolver el sistema linear, de la siguiente ecuación
algebraica:

[Y AA ] [V A ] = [I A] - [Y AB] [V B ] (1.21))
 7.2. ELEMENTOS LINEARES CONCENTRADOS Y DESACOPLADOS

Entre los elementos lineares concentrados desacoplados tenemos la resistencia “R”,

inductancia propia “L” y Capacitancia “C”, son usualmente parte de circuitos
equivalentes, como pueden representar generadores, transformadores, cortas secciones de
líneas de transmisión, u otros componentes de un sistema eléctrico de potencia, ó ellos
pueden representar un componente por si mismo.

7.2.1. Resistance R

Los elementos resistivos son usados para representar algunos efectos entre los cuales
tenemos:
(a) Resistencias de cierre y apertura de interruptores,
(b) Resistencias de pie de torres (como una aproximación gruesa de un cálculo
complicado, de la impedancia de tierra en función de la frecuencia),
(c) Resistencia de tierra de transformadores y generadores con neutro a tierra,
(d) "Mediciones" – ramales en los cuales se necesita medir las corrientes ó voltajes que
no pueden ser calculadas por el programa de la topología directa,
(e) Como partes equivalentes de redes, por ejemplo en paralelo con inductancias lo que
produce un amortiguamiento para su respuesta en frecuencia.
(f) Para la representación de líneas largas en estudios de descargas atmosféricas si es
estudiado la reflexión desde el punto terminal de la línea durante el tiempo tmax
estudiado.

Ejemplo:

Fig. 2.1 – Circuito equivalente de una línea larga sin

reflexiones desde el terminal

La ecuación de una Resistencia Concentrada entre los nodos “k” y “m”,

i km (t) =
1
(vk (t) - vm (t)) (2.1)
R
 Tiene una solución exacta de un programa si tiene un valor razonable, pero no será lo
mismo si tiene un valor muy pequeño.

7.2.1.1. Análisis del Error

Valor de “R” muy Grande, son valores aceptados y no degradan la solución de la red
completa. En el caso limite, R = ∞, es reciproco al valor 1/R y simplemente dará la
pérdidas en [G] de la ecuación (1.8).
Este valor no tiene ninguna influencia en la solución general de la red, siendo un
límite práctico de este valor, R < 1038
Resistencias muy grandes pueden ser usadas para reemplazar circuitos serie con
elementos R-L simétricos en circuitos multifásicos cuando nos interesa observar
únicamente el efecto de acoplamiento capacitivo de una línea.

Fig. 2.2 - Conversión de un circuito “π” de una línea corta con sólo el
acoplamiento capacitivo.

Valores muy pequeños “R”, crea problemas de precisión, valores muy pequeños de “R” crea
valores de conductancias muy grandes 1/R en la matriz [Y] de estado estable y en la matriz
[G] para las soluciones transitorias.
Valores pequeños de resistencias pueden usarse para unir interruptores ó pequeños tramos
de líneas de longitud despreciable, también para simular corto circuitos con característica
resistiva. Actualmente, los softwares existentes ya han superado este problema.

7.2.1.2. Ejemplo de Resistencia en Redes

El propósito de resistencia en las redes es limitador; un simple caso se muestra en la
Fig.2.3. Se asume un voltaje DC desde la fuentes con impedancia de fuente despreciable
es conectado a través de un interruptor con una resistencia de inserción de cierre de valor
“Rclose”

Fig. 2.3 - Energización de una línea larga (Z = impedancia de onda, τ = propagación)

• Contactos cerrados (típicamente de 8 a 10 ms),

• Sistema trifásico,
• La línea con bajas pérdidas,
• Multiples reflexiones.
• Fuentes de impedancia cero.
• Para determinar factor de polos en cierre de interruptores.
En un sistema típico, la máxima sobretensión de maniobra está entre 2.4 a 2.8 p.u. sin
resistencias de inserción de cierre (versus 2.0 p.u. in Fig. 2.3), los cuales se pueden
reducir a 1.5 ó 2.2 p.u. con resistencia de inserción de cierre.

7.2.2. Inductancia Propia “L”

Circuitos magnéticamente acoplados prevalecen en los sistemas de potencia, arranque de
generadores, arranque de transformadores, el acoplamiento magnético de conductores
trifásicos de las líneas. También, hay casos como inductancias propias acopladas, como
las presentadas en la Fig.2.4
(a) Reactor shunt monofásico y reactor de neutro de compensación.
(b) Circuitos de descarga de capacitores serie.
(c) En equipos de HVDC de una estación conversora, para flitros en AC,
(d) Inductancia de circuitos Thevenin equivalente,
 (a) Conexión trifásica (b) Esquema de compensación
normal con cuatro inductores.

Fig. 2.4 – Conexión Shunt de un Reactor.

(e) Inductancia monofásica de circuitos trifásicos conectados en “π” que
representan circuitos balanceados de impedancia de secuencia positiva ó
secuencia cero,
(f) Parte de circuitos de carga (Fig. 2.7), y su modelación en las altas
frecuencias.

Fig. 2.5 – Circuito Thevenin equivalente

con Zpos = Zneg = Zzero

Fig. 2.6 – Representación típica de redes de secuencia positiva

Fig. 2.7 – Modelo de carga para

estudios de armónicos.

(g) Modelos de pararrayos para simulación de la característica dinámica,

 (h) Partes de circuitos electrónicos.

La ecuación de la inductancia propia “L” entre los nodos “k” y “m” es resuelta con
precisión en las condiciones de estado estable con la ecuación (1.11). Sólo debe tenerse
cuidado de que no sean extremadamente pequeña por las mismas razones del caso de la
resistencia. Para simulaciones transitórios la ecuación diferencial exacta será.

di km
vk - vm = L (2.2)
dt

esta es reemplazada por la ecuación diferencial aproximada:

v k (t) - v m (t) + v k (t - ∆t) - vm (t - ∆t) i (t) - i km (t - ∆t)

= L km (2.3)
2 ∆t

La misma ecuación diferencial es obtenida si la regla trapezoidal de integración es aplicada

a la integral.

i km (t) = i km (t - ∆t) + ∫t -∆t [v k (u) - v m (u)] du

1 t
(2.4)
L

dando:
∆t
i km (t) = i km (t - ∆t) + { v k (t) - v m (t) + v k (t - ∆t) - vm (t - ∆t)} (2.5)
2L

Las ecuaciones (2.3) y (2.5) pueden ser reescritas dentro de la ecuación del ramal.
∆t
i km (t) = {vk (t) - vm (t)}+ hist km (t - ∆t) (2.6)
2L

con el término histórico “histkm(t - ∆t)” conocido desde la solución precedida al paso de
tiempo de iteración,
∆t
hist km (t - ∆t) = i km (t - ∆t) + {vk (t - ∆t) - vm (t - ∆t)} (2.7)
2L

Esta ecuación (2.6) puede ser convenientemente representada por una resistencia
equivalente Requiv = 2L/∆t, en paralelo con una fuente de corriente conocida histkm(t - ∆t),
como se muestra en la Fig. 2.8. Primeramente, todos los nodos de voltajes tendrán en
particular un mismo paso de tiempo en el instante “t”, el término histórico de la ecuación
(2.7) será dado por:
 Fig. 2.8 – Circuito resistivo Equivalente
para la solución transitória de
impedancia concentrada.

Cada ramal de inductancia usará en el siguiente tiempo el paso “t + ∆t”. Las corrientes en
los ramales primeramente serán calculadas desde la ecuación (2.6), ó alternativamente si las
ecuaciones son combinadas, usará la formula:

∆t
hist km (t) = {vkm (t) - vm (t)}+ hist km (t - ∆t) (2.8)
L

7.2.2.1. Error Análisis

La ecuación diferencial (2.2) es resuelta en forma aproximada, por lo que es importante
tener el conocimiento acerca del error causado por la aplicación de la regla trapezoidal
de integración.

La regla trapezoidal es numéricamente estable, y afortunadamente hay siempre una

explicación física al error. Para llegar a esta explicación se considerará la
representación de una línea corta que tiene una inductancia de parámetro distribuido L'
por una inductancia de parámetro concentrado por la longitud de la línea.
L′ _ = L (2.9)

Fig. 2.9 – Línea cortocircuitada en los extremos reemplazada por una

 Inductancia concentrada con Z = 2L/∆t y τ = ∆t/2

Siendo, L' un valor conocido, el próximo parámetro a ser determinado será el tiempo
de propagación “τ” por unidad de longitud. Donde:
τ = (L ′ ) (C′ )
(2.10)
El corto tiempo de propagación, será por el menor valor de la capacitancia parásita C’;
donde el tiempo para la simulación transitoria con paso de iteración ∆t es:
∆t
τ= (2.11)
2
Con este valor, la condición en el terminal 1 en “t – ∆t” llegará al final en “t - ∆t/2” y
dará una reflexión al terminal 1 en un tiempo “t”; con lo cual su representación será:
2L ∆t
Z= , τ= (2.12)
∆t 2

El cambio del tamaño del paso de iteración de uno pequeño a otro más grande crea
imprecisiones en los cálculos con más capacitancias parásitas que estarán dadas por la
relación:
(∆t )2
C′ = (2.14)
4L
Pero cuando las ecuaciones de onda sean resueltas con mayor precisión la solución
contendrá menos ruido.

La solución obtenida con la regla trapezoidal contra la solución exacta dado por los
circuito π en cascada es fácilmente susceptible de error como una función de la
frecuencia.

Asumiendo que una inductancia “L” es conectada a una fuente de voltaje de

frecuencia angular “ω”, con una resistencia “R” para propósitos de amortiguamiento;
la simulación transitoria normalmente dará la correcta solución de estado estable. Esta
solución a una frecuencia “ω”, también es obtenida de un circuito π equivalente
exacto. Cortocircuitando el terminal 5; la impedancia de entrada será:
1
Z input = (2.15)
1
Y series + Y shunt
2
ó con la ecuación (1.16),
⎛ ∆t ⎞
tan⎜ ω ⎟
⎝ 2⎠
Z input = jωL •
∆t
ω
2 (2.16)
El radio equivalente de la inductancia aparente resultante de un circuito π en cascada
ó desde la regla trapezoidal y la inductancia exacta tenemos:
 ⎛ ∆t ⎞
tan ⎜ ω ⎟
L trapezoida l ⎝ 2 ⎠
= (2.17)
L ∆t
ω
2
El error de fase es cero sobre el rango completo de frecuencias. Entonces el sistema de
potencia es básicamente operado como voltaje constantes de la red, esto marca la
sensibilidad en algunos casos al asumir este voltaje VL(jω) a través de la inductancia
es mas ó menos fijo, y la corriente IL(jω) le sigue. Si nosotros comparamos esta
corriente del circuito π equivalente en cascada ó la solución de la regla trapezoidal,
con la corriente de la solución exacta de inductancia de parámetro concentrado,
entonces obtendremos la respuesta en frecuencia de la Fig. 2.10 (el recíproco de la
ecuación 2.17), donde la relación es mostrada como una función de la frecuencia de
Nyquist.
1
f Nyquist = (2.18)
2 ∆t

Fig. 2.10 – Relación de Amplitud Itrapezoidal/Iexact a través de una

inductancia como una función de la frecuencia.

Porque del error en la ecuación (2.17), hay una pequeña discrepancia entre la
condición inicial con la ecuación (1.11), y la respuesta de la frecuencia/potencia en
paso de tiempo. Para 60 Hz, este error debe ser 0.012% con ∆t = 100 µs, ó 1.2% con
∆t = 1 ms. Si una perfecta igualdad es deseada, puede tener dos soluciones de estado
estable; uno idéntica por inicialización, y la otra respuesta de estado estable en una ó
más frecuencias.
 Valores Grandes de “L” son aceptados como (ωL)2 ó 2L/∆t ; para obtener el flujo:
flux = ∫ v.dt a través de un ramal una gran inductancia puede adicionarse en paralelo
y la corriente de salida puede obtenerse. La necesidad para este requerimiento puede
presentarse si el ploteo de un flujo de corriente es requerido para una inductancia no
linear; con L = 1010 H.

Valores Pequeños de “L” de “ωL” ó “2L/∆t” origina imprecisiones al igual que las
pequeñas resistencias.

7.2.2.2. Amortiguamiento de la "Oscilación Numérica" con Resistencia

en Paralelo
Mientras la regla trapezoidal filtra las corrientes de alta frecuencia en inductancias
conectadas a fuentes de voltaje, desafortunadamente siempre se amplifican los
voltajes de alta frecuencia a través de inductancias en situaciones donde corrientes
estan saliendo dentro de ellos. En el primer caso, la regla trapezoidal trabaja como
integrador, para el cual es eficiente, en el segundo caso trabaja como diferenciador
para el cual es deficiente. El problema mostrado origina las “Oscilaciones Numéricas”
en casos donde la derivada de la corriente cambia abruptamente.
Por ejemplo, cuando la corriente es interrumpida en un interruptor (Fig. 2.11), la
solución exacta para “VL” es mostrada como una línea contínua con un repentino salto
a cero en el instante de interrupción, donde la solución numérica de los programas se
muestra con una línea punteada.

Fig. 2.11 - Voltage after current interruption

Donde:
v L (t) =
2L
{i(t) - i(t - ∆t)} - v L (t - ∆t) (2.19)
∆t
y asumiendo esta la solución de voltaje será correcta seguido a la interrupción de la
corriente, seguirá esta vL(t) = -vL(t - ∆t) en los puntos 2, 3, 4,... como son las corriente en t -
∆t and t para cero; entonces, la solución para VL será oscilante alrededor de cero con la
amplitud del valor de pre-interrupción.
Hay casos donde el repentino salto deberá ser una respuesta inaceptable, y debería
indicarse un inapropiado modelaje del sistema real. Un ejemplo sería el cálculo de un
transitorio del voltaje de restablecimiento, donde algún interruptor se reinicializará si el
voltaje será máximo con una razón infinita inmediatamente después de la interrupción de
una corriente. Para el cálculo de un transitório de un voltaje de de restablecimiento, la
solución debería incluir la capacitancia propia de dispersión desde el nodo 1 a tierra y
posiblemente desde 1 a 2 y de 2 a tierra.

Hay casos donde al usuario no le interesa los detalles de la rapidez de cambio del
voltaje, y podría estar contento de aceptar resultados con un salto súbito. Un típico ejemplo
podría ser el repentino cambio de voltaje causado por la saturación de un transformador con
el modelo de inductancia de dos pendientes para la no linearidad, como observamos en la
Fig. 2.12.

Fig. 2.12 – Salto de Voltaje causado por la saturación de un transformador.

El "Manual de referancias de respuestas” recomienda ser reinicializado las variables

en el instante del salto de valor. Esto debería ser fácilmente representado si las ecuaciones
son escritas en variables de estado desde [dx/dt] = [A][x]. Con las ecuaciones nodales
usadas en los programas de Transitórios, la reinicialización será muy eficaz.
Kulicke muestra su método de amortiguamiento que soluciona este proble y que
además tiene una base física que hace sólido su planteamiento y que será analizado en el
punto 2.2.3.
 V. Brandwajn [16] y F. Alvarado [17] describieron un método para amortiguar las
“Oscilaciones Numéricas”, que es mediante el uso de una resistencia amortiguadora en
paralelo con la inductancia. Entonces, para una determinada corriente inyectada en el
circuito de la Fig. 2.13, la solución de la regla trapezoidal del circuito paralelo será:

2L
Rp -
v(t) =
1
{i(t) - i(t - ∆t)} - ∆t v(t - ∆t) (2.20)
∆t 1 2L
+ Rp+
2L R p ∆t

Fig. 2.13 – Amortiguamiento Paralelo

En el primer término de la ecuación (2.20) desaparecerá la oscilación y en el Segundo

término (derecha) tendremos el causante de la oscilación numérica:

v(t) = - α • v(t - ∆t)

con:
2L
Rp -
α= ∆t (2.21)
2L
Rp+
∆t

Donde el término “α “ es el llamado factor de oscilación y será amortiguado si: α < 1; como
se muestra en la Fig. 2.14 para Rp = 10⋅(2L/∆t) ó α = 9/11, y para Rp =
2.(2L/∆t) ó α = 1/3. La oscilación deberá disiparse en un paso de tiempo de iteración de
Rp = 2L/∆t ó “α = 0” (caso críticamente amortiguado)3.

3El caso de críticamente amortiguado con la regla trapezoidal con R = 2L/∆t es identico con el métodod de
p
backward Euler, como se muestra en el Anexo 1.9.
 (a) (b)

Fig. 2.14 – (a)Término muy oscilante. (b)Término optimamente amortiguado

Sí “Rp“ es muy grande, el efecto de amortiguamiento es pequeño; pero sí tiene un

valor pequeño de “2L/∆t” (valor ideal para amortiguamiento), introduciría un gran error
dentro de la representación de la inductancia. En la Fig. 2.15 mostramos el error de
magnitud y fase de la impedancia para: Rp = 4.(2L/∆t) y Rp = 8.(2L/∆t), como un error de
magnitud desde la ecuación (2.17) el cual existe alrededor de la inductacia con la regla
trapezoidal. Es interesante observar también, que el error de magnitud con una resistencia
en paralelo es mucho mejor.
 Fig. 2.15 – Error de Fase y Magnitud con resistencia en paralelo[16]

Además, la resistencia en paralelo no tiene efecto determinante en la magnitud de la

respuesta en frecuencia.

Desde un punto de vista puramente numérico, un buen punto de equilibrio entre un

amortiguamiento razonable (Rp lo más bajo posible) y un aceptable error de fase (Rp lo más
alto posible) los valores intermedios estan dentro de los siguientes rangos:

2L 2L
5.4 ≤ R p ≤ 9.4 de acuerdo con Brandwagn [16] (2.22)
∆t ∆t
ó:
 20 2L
Rp = • de acuerdo con Alvarado [17] (2.23)
3 ∆t

El error introducido dentro de la conexión en paralelo en la Fig. 2.13 introducido por

la regla trapezoidal con se ve en la Fig. 2.16, en el cual

Rseries +jXseries = Rp(jX)/(Rp + jX)

Con:
X = ωL,
y con solución para la regla trapezoidal ωLtrapezoidal desde la ecuación (2.17).

Dentro de los programas de transitórios electromagnéticos actuales, la resistencia de

amortiguamiento en paralelo con una inductacia es automáticamente insertada internamente.
Por ejemplo podemos mencionar que el programa SYSCAP para análisis de circuitos
electrónicos de la Rockwell International Corp., emplea técnicas similares a la mostrada,
donde tiene que Rp y Rs de la Fig. 2.18 se contruye dentro del modelo del inductor, con
valores de Rs = 0.1 Ω y Rp = 1012 Ω [22, p. 715].

Esta empresa menciona que, la posibilidad de oscilación numérica se presentaba en

los casos en donde la constante de tiempo del modelo del inductor Fig. 2.18, es más pequeña
comparados con “∆t” [22, p. 773]. Por lo tanto, es mejor escoger un paso de iteración
“∆t” menor que la constante de tiempo de una inductancia para disminuir la probabilidad de
oscilación numérica.
Fig. 2.16 – Resistencias serie aparente e inductancias serie para la conexión en paralelo
de la Fig. 2.13, con Rp = (20/3) • (2L/∆t). La region de la derecha de
f/fNyquist = 0.5 es de poco interes práctico porque el muestreo debería
también bajar con frecuencias adecuadas.
Fig. 2.17 Relación “L/R” de la impedancia de corto circuito de transformadores
típicos [18]. Gráfico tomado de la CIGRE.

Fig. 2.18 – Circuito equivalente para la impedancia

de corto circuito de un transformador.
 7.2.2.3. Razones Físicas para Resistencia en Paralelo

Hay algunos casos en las cuales la inductancia debe tener una impedancia en
paralelo por razones físicas. A continuación damos algunos ejemplos de aplicación en la
cual el amortiguamiento de la oscilación ó ruido numérico se da por el modelo físico de
equipo a representar.
Desde luego, que éstos ejemplos no cubre todas las aplicaciones pero serán las mas
representativas que se mostrarán a continuación:

(a) Impedancia de Corto Circuito de transformadores

La impedancia de corto circuito de un transformador no tiene la constante de
relación L/R; esta relación L/R decrece con un incremento en frecuencia, como se
muestra en la Fig. 2.17. Si nosotros tomaramos la curva de un transformador de 100
MVA y asumimos L = 1H (ó mH, ó p.u.) con un ∆t = 100µs, donde el valor
de Rp = 163,000 Ω (ó mΩ, ó p.u.) será producto de la propia relación L/R a 1 kHz,
estando este valor entre el rango de los límites de 108,000 Ω y 188,000 Ω
recomendados en la ecuación (2.22) mencionada anteriormente. La resistencia serie
de Rs = 9.4 Ω (ó mΩ, ó p.u.) puede ser adicionando obteniendo la corrección de la
relación L/R a 50 Hz, el cual empleará el circuito equivalente de la Fig. 2.18 para la
impedancia de corto circuito del transformador. Con Zinput desde la ecuación (2.16),
la relación Ltrapezoidal/R de este circuito equivalente es mostrada con línea punteada en
la Fig. 2.19, el cual es una buena aproximación a la curva experimental (con línea
solida), y mucho mejor a la constante de relación L/R sin Rp4. Esto es interesante
porque también el grupo de trabajo del CIGRE en problemas de interface
recomienda el mismo circuito equivalente de la figura Fig. 2.18 para el análisis de
armónicos [9], con:
SN Rp
13 < 2
< 30 (2.24)
VN
y
2
90 < V N < 110 (2.25)
S N Rs

Donde:
SN= Potencia y VN= voltaje del transformador.

4In the EMTP, the L/R-ratio without R would actually increase with frequency, since L
p trapezoidal of Eq. (2.17)
increases with frequency.
 Fig. 2.19 – Relación L/R de la impedancia de corto circuito de un
transformador de 100 MVA (línea punteada desde el
circuito equivalente de la Fig. 2.18 con “L” resuelto por la
regla trapezoidal, línea solida desde [18]). Representación
recomendada por la CIGRE.

Si asumiéramos Xshort-circuit = 0.05 ó 0.10 p.u. a 50 Hz, en la ecuación (2.24) se

convierte con ωL = (0.05 to 0.10) ⋅ VN2/SN,
(40,841 a 81,681) ⋅ L < Rp < (94,248 a 188,496) ⋅ L (2.26)
ó con un paso típico de: ∆t = 100 µs,
(2.04 to 4.08) 2L/∆t < Rp < (4.07 to 9.4) 2L/∆t,
con el gran número para bajas reactancias de corto circuito de 0.05 p.u. El valor de
Rp esta ene el mismo rango de las ecuaciones (2.22) y (2.24)
La ecuación (2.25) implica ua relación L/R a 50 Hz de 0.014 a 0.018 para una
reactancia de corto circuito de 0.05 p.u., ó de 0.028 a 0.035 para una reactancia de
corto circuit de 0.10 p.u., el cual es menor para valores a 50 Hz en la Fig. 2.17.

(b) Impedancia de Magnetización de Transformadores y Núcleo de Fierro de los

Reactores.
Resistencias en paralelo son adicionadas en la inductacia de magnetización de
transformadores para representar una gruesa aproximación de la Histéresis y las
pérdidas por corrientes de Eddy . En forma similar, el circuito equivalente de la Fig.
2.18 es recomendado para reactores con núcleo de fierro, con Rs representando las
 pérdidas en el bobinado “I2R” y “Rp “ representando las pérdidas en el núcleo de
fierro.

(c) Generador Síncrono

El Grupo de Trabajo del CIGRE sobre problemas de Interfase recomienda una
resistencia en paralelo con la inductancia de secuencia negativa (el cual es
prácticamente idéntica a (L"d + L"q). Una curva típica para la relación L"d/R de
generadores, similar a la Fig. 2.19 fue publicada, y deberían ser usados con valores
razonables de Rp.

(d) Circuito π nominal

Conexiones en cascada de un circuito π nominal son usadas para representar una
línea de transmisión en un analizador de transitorio de redes. La eliminación de las
oscilaciones numéricas que son causadas por la aproximación de las inductancias
concentradas a los parámetros distribuidos, estan superadas al adicionar una
resistencia en paralelo (ver Fig. 2.20). Valores típicos son Rp = 5 Zsurge, el cual
deberá ser calculado de Rp = 5. (2L/∆t) en el tramo de línea representada, con la
inductancia en la Fig. 2.9.

.
Fig. 2.20 –Resistencia de amortiguamiento en un circuito π nominal.

(e) Impedancia de Fuentes

El circuito equivalente Thevenin de la Fig. 2.5 es obviamente una primera
aproximación del resto del sistema a frecuencia diferente desde la relación
potencia/frecuencia.
Para hacer la respuesta en frecuencia de un circuito más realista, resistencias de
amortiguamiento son conectados en paralelo con los ramales inductivos “L”.
 7.2.3. Capacitance C
Los elementos capacitivos son usados para representar, algunos components del sistema
como:
(a) Condensadores serie y shunt,
(b) Condensadores shunt en representación de circuitos π nominal de líneas de
transmisión,
(c) Equipos de estaciones conversoras en HVDC, el cual como parte de circuitos
snubber y filtros,
(d) Capacitancias de dispersion de transformadores, generadores, etc., especialmente en
transitórios de voltaje de restablecimiento y estudios de descargas atmosféricas,
donde la impedancia “ωL” se vuelve tan alta a las altas frecuencias y las
capacitancias de dispersión de valor “1/ωC” son dominates al encontrase en paralelo
con la impedancia de transformador,
(e) Capacitancia de los transformadores de tensión y divisores de voltaje,
(f) Partes de generadores de onda.

La ecuación de la capacitancia concentrada entre los nodos “k” y “m” es resuelta

exactamente en la solución de estado estable AC con la ecuación (1.12).. Se debe tener
cuidado que “ωC” no debe ser extremadamente grande, por las mismas rezones que como se
explico para resistencia pequeña en la sección 2.1.1.
Para la simulación transitória, la ecuación diferencial exacta es:
d( v k - v m )
i km = C (2.27)
dt
luego, es reemplazada por aproximación a la ecuación diferencial:

i km (t) + i km (t - ∆t) { v (t) - vm (t)} - { v k (t - ∆t) - vm (t - ∆t)}

=C k (2.28)
2 ∆t
donde la ecuación para un nodo con varios ramales, tendremos:
2C
i km (t) = { v k (t) - v m (t)} + hist km (t - ∆t) (2.29)
∆t
con el "término histórico" histkm(t - ∆t) conocido desde la solución previa al time step,
2C
hist km (t - ∆t) = - i km (t - ∆t) - { v k (t - ∆t) - v m (t - ∆t)} (2.30)
∆t
Luego, analogo a la inductancia, un resultado identico se obtendrá desde un integrador de la
ecuación (2.27) con la regla trapezoidal. La ecuación (2.29) puede ser representada como
una resistencia equivalente Requiv = ∆t/2C, en parallelo con una fuente de corriente
conocida histkm(t - ∆t), como se muestra en la Fig. 2.22.
 Fig. 2.22 – Circuito equivalente resistivo para una
solución transitória de capacitancia concentrada.

En un instante “t”, el término histórico de la ecuación (2.30) será para cada ramal capacitivo
por uso del próximo paso “t + ∆t”, primero hallará la corriente desde la ecuación (2.29). Asi
se tiene:
C
hist km (t) = - { vk (t) - v m (t)} - hist km (t - ∆t) (2.31)
∆t
El cual, es el mismo de la ecuación (2.8) para la inductancia si seguimos por una señal
reversa.

Fig. 2.23 – La capacitancia de parámetro concentrado

reemplaza a toda la porción de línea de parámetros
distibuídos con Z = ∆t/2C y τ = ∆t/2.

7.2.3.1. Análisis del Error

Se efectuara un análisis de error similar al de la inductancia; para una representación
física del error, se presenta una porción de línea en la Fig. 2.23, en el cual la capacitancia
concentrada es reemplazada por una línea abierta sin pérdidas.. El parámetro obtenido es
razonablemente equivalente al total de la capacitancia distribuída en la línea, igual a la
capacitancia concentrada,
C′ = C (2.32)
Con C' conocido, el próximo parámetro a ser determinado es el tiempo de viaje “τ”.
La ecuación (2.10) mostrada en el corto tiempo de viaje, el menor será el valor de la
"parasita". Para pasos de iteración “∆t”, el corto tirmpo de viaje de la onda será:
 ∆t
τ= (2.33)
2
Con las ecuaciones (2.32) y (2.33) la impedancia de onda será: Z = ∆t/2C.

Igualmente al caso de la inductancia, existe un error de cálculo que en el dominio de

la frecuencia se incrementará. Asumiendo que la capacitancia C esta conectada a una fuente
con frecuencia angular “ω”, a traves de alguna red con amortiguamiento; la solución de
estado estable es conocida desde el circuito π equivalente exacto de la Fig. 1.2, con terminal
5 abierto,
2
1 Y series
Y input = ( Y series + Y shunt ) - (2.34)
2 1
Y series + Y shunt
2
después de algunos arreglos la ecuación (1.16) será,
⎛ ∆t ⎞
tan ⎜ ω ⎟
⎝ 2⎠
Y input = jωC • (2.35)
∆t
ω
2
Luego el error análogo a la ecuación (2.16) aplicado para la capacitancia, será:
⎛ ∆t ⎞
tan ⎜ ω ⎟
C trapezoidal
= ⎝ 2⎠ (2.36)
C ∆t
ω
2
El error es cero sobre el rango de frecuencia entero. Si nosotros forzamos una
corriente Ic(jω) dentro de la capacitancia, el voltaje a través de la porción de línea,
comparado con la solución exacta tendrá la respuesta en frecuencia de la Fig. 2.24, el cual es
idenco con Fig. 2.10 si la relación de corriente es reemplazado por la relación de voltajes. La
regla trapezoidal filtra los voltajes de alta frecuencia.
 Fig. 2.24 – Relación de Amplitud Vtrapezoidal/Vexact de la capacitancia
como una función de la frecuencia.

Hay pequeñas discrepancias entre la condición inicial con la ecuación (1.12), y la

respuesta a potencia/frecuencia desde la ecuación (2.35) en el paso de iteración del lazo (a
60 Hz, 0.012% error con ∆t = 100 µs ó 1.2% con ∆t = 1 ms)

Valores muy pequeños de “C” son aceptables como grandes impedancias (1/ωC)2 ó
∆t/2C que pueden ser calculados sin mayores problemas. Valores muy grandes de “C”
crea problemas de precisión iguales al de pequeñas resistencia, pero en la práctica tienen
poca probabilidad de ocurrir.

7.2.3.2. Amortiguamiento de "Oscilaciones Numéricas" con Resistencia

Serie
La regla trapezoidal filtra los voltajes de alta frecuenciaa través de la capacitancia
para una determinada corriente de inyección, sin embargo amplifica las corrientes las
corrientes de alta frecuencia para un determinado voltaje a través de “C”.
La oscilación numérica discutada en para la inductancia en la sección 2.2.2 se aplica
para las corrientes capacitivas si hay un cambio brusco en dVc/dt; por lo que los problemas
de oscilación numérica también estan presentes en las capacitancias.

En forma anláloga a la inductancia, estas oscilaciones numéricas deben ser

controladas también con resistencias de amortiguamiento, las que serán conectadas en serie
de valor “RS” (ver Fig. 2.25). Usando la teoría de Alvarado, la soución de la regla
trapezoidal para un impulso de voltaje aplicado al circuito de la Fig. 2.25 será:
 ∆t
- Rs
1
i(t) = {v(t) - v(t - ∆t)} - 2C i(t - ∆t)
∆t ∆t
+ Rs + Rs
2C 2C

Fig. 2.25 – Resistencia Serie Amortiguamiento.

Después que el impulso de voltaje “V” cae a cero, el segundo término de la derecha es el
que causa la oscilación numérica,

i(t) = - α • i(t - ∆t)

∆t
- RS
con : α= 2C (2.38)
∆t
+ RS
2C

en analogía a la ecuación (2.23), un razonable valor de resistencia de amortiguamiento será

el siguiente:
∆t
R s = 0.15 (2.39)
2C

7.2.3.3. Razones Físicas de la Resistencia Serie

G.W.A. Dummer, sugiere el circuito equivalente de la Fig. 2.26, y con un “Rs”
dominante para las muy altas frecuencias, donde “Rp” es dominante en muy bajas
frecuencias,

Fig. 2.26 – Circuito equivalente para

un capacitor con pérdidas.
 Fig. 2.27 – Factor de pérdidas. Curva obtenida por Springer-Verlag y A.W. Roth
Un condensador puede ser sujeto a corto circuitos often have series resistors built in.
Similarmenete, la protección de sobrevoltaje de condensadores serie con espinterómetro
(Fig. 2.28) incluye limitadores de corriente con elementos “R-L” en el circuito de descarga,
con un típico "zumbido" durante la descarga a una frecuencia de 400 Hz.

Fig. 2.28 – Espinterómetro como protección de capacitors serie.

7.2.3.4. Example for Network with Capacitances

Modifiquemos el estudio de la corriente de avería de la sección 2,2,4 para un caso en el cual

la línea de la transmisión serie-se compense con los condensadores (Fig. 2.29). Más aun
asumamos que el Ls en la sección 2,2,4 representó la reactancia neta Xnet = ωL - 1/ωC at 60
Hz, para hacer ambos resultados directamente comparables. Con Rs = 0.18 p.u., Xs =
1.0833 p.u., ωC = 2.695 p.u. y ∆t = 100 µs, la corriente de falla de fig. 2,30 se obtiene (los
datos tomados de [ 74 ], con la conexión de la localización de falla al bus infinito a la
izquierda apagado). Para los propósitos de la comparación, la corriente de falla con la
reactancia neta representada por Ls, según lo hecho en la sección 2,2,4, se demuestra
también; diferencia apreciable de la solución más exacta con el modelo del circuito de fig.
2,29. Esta diferencia tiene consecuencias para la exactitud de las simulaciones de la
estabilidad, puesto que los reactancias netos se utilizan prácticamente siempre en estudios de
la estabilidad. La fig. 2,31 compara las curvas del oscilación obtenidas con una reactancia
neta y una representación del L-c para un caso similar al modelo de la prueba patrón de
IEEE para la resonancia subsincrona estudia [ 21 ].

Fig. 2.29 – Sistema monofásico con falla a tierra en un

sistema con capacitor en serie

Fig. 2.30 – Corriente de falla en un sistema de trabajo de compensador serie de la

Fig. 2.29 (linea sin simbolos). Por comparacion resulta de la Fig. 2.21
con la reactancia del mostrado de la siguiente manera (linea sin
simbolo)
 Fig. 2.31 – Curvas de Oscilación con representación “R-L” y “R-L-C”.

7.2.4. Series Connection of R, L, C

Si son combinados los elementos R, L, C con frecuencia en pares como conexiones de

serie de R-l, de R-c, o del L-c, o como conexión de serie de los tres elementos R-l-c,
entonces llega a ser más eficiente tratar la conexión de serie como solo rama, de tal modo
reduciendo el número de nodos y de ecuaciones nodales. Esto se ha puesto en ejecucio'n
en el EMTP para la conexión de serie de R-l-c (fig. 2,32). Para la solución de estado
estacionario, la ecuación del rama está simplemente

Fig. 2.32 – Conexión Serie de R, L, C

 1
I km = (V k - V m )
R + j( ωL - 1/ωC)

La derivada de la ecuación del ramal de la simulación transitória, sumando la caída de

tensión a través de R, L, y C:
v k - v m = v R + v L + vc

con la caída de tension expresado como una función de la corriente con la ecuación (2.1),

(2.6) y (2.29):

⎛ 2L ∆t ⎞ 2L ∆t
vk (t) - v m (t) = ⎜ R + + ⎟ i km (t) - hist L (t - ∆t) - hist C (t - ∆t) (2.40)
⎝ ∆t 2C ⎠ ∆t 2C

Después reemplazando el término histórico histL y histC con las expresiones de las
ecuaciones (2.7) y (2.30), la ecuación de este ramal será:

i km (t) = G series { v k (t) - v m (t)} + hist series (t - ∆t) (2.41a)

con:
1
G series = (2.41b)
2L ∆t
R+ +
∆t 2C

y la combinación de términos históricos:

⎧⎛ 2L ∆t ⎞ ⎫
hist series (t - ∆t) = G series ⎨⎜ -R- ⎟i(t - ∆t) + v k (t - ∆t) - v m (t - ∆t) - 2 v C (t - ∆t)⎬. (2.42)
⎩⎝ ∆t 2C ⎠ ⎭

For updating this history term, la nueva corriente es primeramente calculado desde la
ecuación (2.41a), y el nuevo voltaje en el capacitor “vC“ será:

∆t
vC (t) = vC (t - ∆t) + {i(t) + i(t - ∆t)}
2C
La ecuación (2.42) no esta solamente dada por la expression histórica de la combinación de
términos, pero será usado por los programas de transitórios electromagnéticos.

8. INTERRUPTORES

8.1. Conceptos Básicos

El interruptor es un dispositivo destinado al cierre y apertura de la

continuidad de un circuito eléctrico bajo carga, en estado normal y de corto
circuito.

Su operación determina el estado de confiabilidad que se puede tener en un

sistema de potencia.

El interruptor debe ser capaz de interrumpir corrientes de diferentes

intensidades y factores de potencia, pasando de corrientes capacitivas hasta
las corrientes inductivas.

8.2. PARTES

Parte Activa
conformada por las cámaras de extinción que soportan los contactos fijos y
el mecanismo de operación que soportan los contactos móviles.
 Parte Pasiva

formada por una estructura que soporta uno ó tres depósitos de aceite, en
casos que el interruptor sea de aceite, en el que se aloja la parte activa.

- Protege eléctrica y mecánicamente el interruptor,

- Soporte para el gabinete de control,
- Ofrece puntos para levantamiento y transporte.

8.3. ACCESORIOS

1.- Conectores
2.- Placa de datos
 3.-Gabinete con dispositivos de control, protección y medición, como:
compreseora, resorte, bobinas de cierre ó de disparo, resistencia de
Calefacción, y tros.

EL ACCIONAMIENTO DE LOS DISPOSITIVOS DE CONTROL

PUEDEN SER DE TIPO NEUMÁTICO, ELECTROHIDRAÚLICO Y
DE RESORTE, DEPENDIENDO DEL NIVEL DE TENSIÓN AL
CUAL VAYAN A OPERAR.

8.4. PARÁMETROS DE LOS INTERRUPTORES

1.-TENSIÓN NOMINAL
Es el valor eficaz de la tensión entre fases del sistema en que se instala el
interruptor.
2.-TENSIÓN MÁXIMA
Es el valor máximo de la tensión para el cual está diseñado y el límite
máximo al que puede operar.
3.-CORRIENTE NOMINAL
Es el valor eficaz de la corriente máxima que puede circular
contínuamente a través de los contactos sin exceder los límites de
elevación de tenperatura.
4.-CORRIENTE DE CORTO CIRCUITO INICIAL
Es el valor pico de la primera semionda de corriente, comprendida en
ella la componente transitória.

5.-CORRIENTE DE CORTO CIRCUITO

 Es el valor eficaz de la corriente máxima de cortocircuito que pueden
interrumpir las cámaras de extinción del arco. Las unidades son
Kiloamperios, pero comunmente se dan en Potencia de cortocircuito
(MVA).
6.-TENSIÓN DE RESTABLECIMIENTO
Es el valor eficaz de la tensión máxima de la primera semionda de la
componente alterna, que aparece entre los contactos del interruptor
después de la extinción de la corriente. Tiene una influencia muy
importante en la capacidad de apertura del interruptor y presenta una
frecuencia que es del orden de miles de Hertz, de acuerdo con los
parámetros eléctricos del sistema en la zona de operación.

Esta tensión tiene dos componentes: una a la frecuencia nominal del

sistema y la otra superpuesta que oscila a la frecuencia natural del
sistema.
7.-RESISTENCIA DE CONTACTO
Es el contacto metálico de área muy pequeña por el que pasa la corriente
nominal ó de falla y que origina un calentamiento por esta superficie.
La Resistencia de Contacto varía de acuerdo con la fórmula:
R = K.ρ.(D/F)

R = Reistencia de contacto en Ω.
K = Constante dada por el fabricante del interruptor
ρ = Resistividad del metal del contacto.
D = Dureza del metal del contacto.
F = Fuerza que mantiene cerrado el contacto.
 8.-CAMARA DE EXTINCIÓN DEL ARCO
Es la parte principal de cualquier interruptor, en donde al separarse los
contactos se transforma en calor la energía que circula por el cicuito del
que se trate.
Dichas cámaras deben soportar los esfuerzos electrodinámicos de las
corrientes de cortocircuito, así como los esfuerzos dieléctricos.
El fenómeno de interrupción se iniciarse con la separación de los
contactos, apareciendo un arco a través de un fluído, que lo transforma
en PLASMA.
La corriente se interrumpe durante el paso de la onda por cero y
posteriormente aparece entre los contactos una tensión transitória de
restablecimiento.

8.5. INTERRUPCIÓN DEL ARCO

Aparecen fenómenos:
1.-Altas Temperaturas debido al plasma creado por el arco;
2.-Altas presiones debido al punto anterior,
3.-Flujos turbulentos del gas que adquieren velocidades variables entre 100
y 100 m/s y que producen el soplado del arco, su alargamiento y extinción.
4.-Masa metálicas en movimiento(contacto móvil) que se aceleran en
pocos milésimos de segundo hasta llegar a velocidades del orden de 10 m/s.
5.-Esfuerzos mecánicos debidos a la corriente de cortocircuito.
6.-Esfuerzos dieléctricos debidos a la tensión de restablecimiento.
8.6. PRUEBAS DE INTERRUPTORES
1.-Falla en los Terminales; En este caso la oscilación de la tensión se
amotigua por la resistencia propia del circuito de potencia y “fo”,
dependiendo de los valores de Inductancia y Capacitancia del lado de la
fuente:
fo = 1/2√πLC
2.-Falla kilométrica; efectuada en una línea Corta, considerada como crítica
para los interruptores, cuando sucede a los 3 a 5 km de distancia. La tensión
de restablecimiento esta dada por la diferencia de tensión entre el lado de la
fuente y el lado de la línea, con una frecuencia de oscilación del doble de
fundamental.
3.-Apertura en Oposición de Fases; se produce en el caso de tener una fase
mal conectada ó en pérdida de estabilidad de sistemas.El corto circuito
producido llega ocasiona una sobretensión hasta tres veces la tensión fase-
neutro con una frecuencia el doble de la fundamental.
4.-Apertura de Pequeñas corrientes Inductivas; La apertura de estas
corrientes puede provocar la llamada “Falla Evolutiva” que en un
interruptor puede ser severa.
5.-Falla Evolutiva; ocurre cuando al abrir un interruptor aparece una
sobretensión que ocasiona el contorneo de aisladores externos.
8.7. RESISTENCIA AMORTIGUAMIENTO DE TVR

o En interruptores con varias cámaras de extinción, los resistores ayudan

a distribuir la tensión de restablecimento transitória de forma uniforme
entre todas las unidades de interrupción;
o Los resistores reducen la severidad de la tensión de restablecimento
transitória, en los instantes de la interrupción;

8.7.1. RESISTENCIA AMORTIGUAMIENTO DE TVR

o Utilizando princípios conocidos es Vm

posible obtener la siguiente expresión LC
operacional para el TRV: V (S ) =
⎛ 2 S 1 ⎞
S ⎜⎜ S + + 2 ⎟⎟
⎝ TP T ⎠
o Donde:
Vm - Tensión de pico fase-neutro;
TP = RC;
R - Resistencia de Inserción;
C - Capacitância de Fuga;
T2 = LC = 1/(ω02) = T0/(2π);
`L - Indutancia del alimentador;
ω0 - frecuencia angular del circuito oscilador;
T0 - período del oscilador;
Introduciéndo una nueva variable η (adimensional) tal que:

R R TP
η= = =
En la ecuación L Z 0 T operacional para o TRV habrá las
tres hipóteses C para a variável η :
9. MANIOBRAS DE BANCO DE CONDENSADORES

9.1. ENERGIZACIÓN DE BANCO DE CONDENSADORES

9.2. LIMITAR INRUSH DE BCO.CONDENSADORES
9.3. ENERGIZACIÓN DE VARIOS BCOS. DE CONDENSADORES
9.4. DESCONEXIÓN DE BANCOS DE CONDENSADORES
9.5. EJEMPLOS DE APLICACIÓN
10. EFECTOS ARMÓNICOS