GICEP

Grupo de Investigación en Calidad de la Energía y Electrónica de potencia

150

10
0

Profesor: Armando Jaime Ustariz Farfán

0

E-Mail: ajustarizf@[Link]
-
5
0

-
10
0

-
150

0 2 4 6 8 1 12 14
0
t
[
m
s
]

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz).
CAPÍTULO 1: Conceptos Generales

1.1. Introducción

GICEP
Grupo de Investigación en Calidad de la Energía y Electrónica de potencia

150

10
0

Profesor: Armando Jaime Ustariz Farfán

0

E-Mail: ajustarizf@[Link]
-
5
0

-
10
0

-
150

0 2 4 6 8 1 12 14
0
t
[
m
s
]

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz).
 Introducción : Técnicas de Análisis

MÉTODOS DE ANÁLISIS

DOMINIO DEL TIEMPO DOMINIO DE LA FRECUENCIA

HÍBRIDOS INTELIGENTES
muy flexibles poco flexibles
flexibles flexibles
tiempos largos muy rápidos
relativamente rápidos rápidos
[1] [2]

ACERCAMIENTO INYECCIÓN INTERACCIÓN

RÁPIDO AL DE ARMÓNICOS DE ARMÓNICOS
RÉGIMEN sencillos complejos
PERMANENTE aproximados aproximados

[1] N. Watson & J. Arrillaga, “Power Systems Electromagnetic Transients Simulation”, IET Power and Energy Series, Vol. 39, 2007.
[2] E. Acha & M. Madrigal, “Power Systems Harmonics Computer Modelling and Analysis”, John Wiley & Sons, 2001.

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 3
 Introducción : Técnicas de Análisis de sistemas eléctricos

MÉTODOS DE ANÁLISIS

DOMINIO DEL TIEMPO DOMINIO DE LA FRECUENCIA

u1 (t ) R2 u2 (t ) U1 ( ) R2 U 2 ( )

i1 (t ) R1 R3 i2 (t ) I1 ( ) R1 R3 I 2 ( )

 i1 (t )   R1  R2  R12   u1 (t )   I1 ( )   R1  R2  R12  U1 ( ) 

1 1 1 1

i (t )     1 
   I ( )     1 

 
 2   R2    2   R2  
1 1 1 1
R2  R2 
R3  u 2 (t ) R3  U 2 ( )

i (t )  Y  u (t ) I ( )  Y   U ( )
Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 4
 Introducción : Técnicas de Análisis de sistemas eléctricos

MÉTODOS DE ANÁLISIS

DOMINIO DEL TIEMPO DOMINIO DE LA FRECUENCIA

u1 (t ) R2 u2 (t ) U1 ( ) R2 U 2 ( )

i1 (t ) R1 R3 i2 (t ) I1 ( ) R1 R3 I 2 ( )

u (t )  Y 1  i (t ) U ( )  Y 1  I ( )

R1  R2  R3  R1  R2  R3 

 u1 (t )   R1  R2  R3
R1R3
R1  R2  R3
  i (t )  
U1 ( )   R1  R2  R3
R1R3
R1  R2  R3
  I ( ) 
u (t )    R1R3  1  U ( )    R1R3  1
 2   R1  R2  R3
R3  R1  R2 
 i2 (t )   2   R1  R2  R3
R3  R1  R2 
  I 2 ( ) 
R1  R2  R3  R1  R2  R3 

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 5
 Introducción : Técnicas de Análisis de sistemas eléctricos

MÉTODOS DE ANÁLISIS

DOMINIO DEL TIEMPO DOMINIO DE LA FRECUENCIA

1
C u2 (t ) U1 ( ) jC U 2 ( )
u1 (t )

i1 (t ) L R i2 (t ) I1 ( ) j L R I 2 ( )

diL (t ) duC (t )
uL (t )  L iC (t )  C
 I1 ( )   R  j L   U1 ( ) 
dt 1 1 1
dt j L
 I ( )     1  
j C  j L  U 2 ( ) 
 2   j L
1

¿…?
Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 6
CAPÍTULO 1: Conceptos Generales

1.2. Dominio del Tiempo

GICEP
Grupo de Investigación en Calidad de la Energía y Electrónica de potencia

150

10
0

Profesor: Armando Jaime Ustariz Farfán

0

E-Mail: ajustarizf@[Link]
-
5
0

-
10
0

-
150

0 2 4 6 8 1 12 14
0
t
[
m
s
]

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz).
 Dominio del tiempo: Método de Sustitución Numérica

Regla de Integración Trapezoidal:

dx  t  t
f  d  
x (t )
dx  
Sea: f t  
dt
  t t x ( t t )

f   f t 

Aplicando la regla:
f  t  t 
 f  t   f  t  t  
t     x  t   x  t  t 
 2 

 t  t   t t
x t   f t   f  t  t   x  t  t 
t

2 2
Área de un trapecio:
histórico
A  b
h  H 
2
Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 8
 Dominio del tiempo: Método de Sustitución Numérica

ikm (t ) R
1
ikm (t )  ukm (t ) Diferencia de potencial: ukm (t )  uk (t )  um (t )
uk (t ) um (t ) R

ikm (t ) L
dikm (t ) 1 1 t ikm ( t )

uk (t ) um (t ) dt
 ukm (t )
L

L t t
u km ( ) d  ikm ( t  t )
di ( )

f   d   dx  
t x (t )
t t x ( t t )

t t
Regla de Integración Trapezoidal: x  t   f t   f  t  t   x  t  t 
2 2

histórico

t t
I hL ikm (t )  ukm (t )  ukm (t  t )  ikm (t  t )
ikm (t ) 2L 2L
1 I hL
uk (t ) um (t )
RL
RL
Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 9
 Dominio del tiempo: Método de Sustitución Numérica

ikm (t ) C
dukm (t ) 1 t ukm ( t )

uk (t ) um (t )
ikm (t )  C
dt

C t t
ikm ( ) d  u km ( t t )
du ( )

f   d   dx  
t x (t )
t t x ( t t )

t t
Regla de Integración Trapezoidal: x  t   f t   f  t  t   x  t  t 
2 2

t t
ukm (t )  ikm (t )  ikm (t  t )  ukm (t  t )
2C 2C

2C 2C
I hC ikm (t )  ukm (t )  ukm (t  t )  ikm (t  t )
ikm (t ) t t
1 I hC
uk (t ) um (t )
RC
RC

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 10
 Dominio del tiempo: Método de Sustitución Numérica

Equivalente Norton
t
I hL I hL (t  t )  uL (t  t )  ikm (t  t )
ikm (t ) R 2L
uL (t  t )  ukm (t  t )  R ikm (t  t ) 
uk (t ) um (t )
RL t  t 
I hL (t  t )  ukm (t  t )  1  R  ikm (t  t )
2L  2L 
Req  R  RL

RL
iSC (t )  I hL (t  t )
R  RL

I hRL
ikm (t )
RL
I hRL (t  t )  I hL (t  t )
uk (t ) um (t ) R  RL
R  RL

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 11
 Dominio del tiempo: Ejemplo 1 – RL con fuente DC

t  0,002
u (t )
Fuente DC: Para t  0,002
iR (t ) iL (t )
iR (t )  i1 (t )
i (t ) R L
u (t )  Ri (t )  u (0, 002)  10

Para t  0,002

i1 (t )  iR (t )  iL (t )
i1 (t ) 10
R  1 1
i (t )  u (t )  iL (t )
R
L  1mH
1  diL 
iL (0, 002)  0 i (t )  L   iL (t )
R  dt 
diL R R
 iL (t )  i (t )
dt L L
Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 12
 Dominio del tiempo: Ejemplo 1 – RL con fuente DC

t 0
u (t )
Fuente DC:
i1 (t ) 10
iR (t ) iL (t )
R  1
i (t ) R L L  1mH

iL (0, 002)  0

diLn R  t L
Respuesta Natural:  iLn (t )  0  iLn (t )  Ae ; 
dt L R

iLf  t   i (t )
diL R R
 iL (t )  i (t ) Respuesta Forzada:
dt L L

 RL t
Respuesta Completa: iL (t )  iLn (t )  iLf (t )  Ae  i (t )

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 13
 Dominio del tiempo: Ejemplo 1 – RL con fuente DC

u (t ) Solución exacta:
 RL t
iL (t )  (iL (0, 002)  i(0, 002))e  i (t )
Fuente DC:
iR (t ) iL (t )

i (t )  10 R L u (t )  R  i (t )  iL (t ) 

Solución regla de integración trapezoidal:

u (t )
 
t 0
i (t )  I L  1
R  R1L u (t )
iR (t ) iL (t )
i (t )  I hL
i (t )  10 u (t ) 
R RL I hL
 1
R  R1L 
i (t )   2Lt u (t  t )  iL (t  t ) 

 R1  2Lt 
Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 14
Dominio del tiempo: Ejemplo 1 – RL con fuente DC

t 1ms

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 15
Dominio del tiempo: Ejemplo 1 – RL con fuente DC

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 16
 Dominio del tiempo: Ejemplo 2 – RLC con fuente AC

u1 (t ) C u2 (t )

iL (t ) iR (t )
Para t  0,002

u1 (t ) C u2 (t )
i1 (t ) L R i2 (t )
iL (t ) + - iR (t )
t  0,002 uc (t )
i1 (t ) L R i2 (t )

i1 (t ) 10sin t  5  R  2

i2 (t )  5sin t  45  L  1mH

f  60 Hz C  10uF iL (t )  i1 (t ) iR (t )  i2 (t )

u1 (t )  ? u2 (t )  R  i2 (t )
iL (t )  i1 (t )  iL (0, 002)  7,5 A

uC (t )  u1 (t )  u2 (t )  vC (0, 002)  7,3V 

u1 (t )   L 10sin t  95 

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 17
 Dominio del tiempo: Ejemplo 2 – RLC con fuente AC

Modelo discreto del inductor

Para t  0,002 I hC
u1 (t ) u2 (t )

t
I hL  ukm (t  t )  ikm (t  t ) RC
2L
i1 (t ) RL I hL R i2 (t )
2L
RL 
t

Modelo discreto del capacitor

I hC  
2C
ukm (t  t )  ikm (t  t ) i1 (t ) 10sin t  5  R  2
t
i2 (t )  5sin t  45  L  1mH
t
RC  f  60 Hz C  10uF vC (0, 002)  7, 4
2C
iL (0, 002)  7,5

i1 (t )  I hL  I hC   RL  RC  R1C   u1 (t ) 
u (t )  Ybus  i(t )
1 1
1
 i (t )  I   1 
 
 2   RC  
1 1
C R 
RC 
u 2 (t )

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 18
Dominio del tiempo: Ejemplo 2 – RLC con fuente AC

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 19
Dominio del tiempo: Ejemplo 2 – RLC con fuente AC

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 20
 Dominio del tiempo: Ejemplo 3 – Diferentes frecuencias

u1 (t ) C u2 (t )

iL (t ) iR (t )

i1 (t ) L R i2 (t )

t  0,002

i1 (t ) 10sin t  5  R  2

i2 (t )  5sin  5t  45  L  1mH

f  60 Hz C  10uF

iL (t )  i1 (t )  iL (0, 002)  7, 2 A

uC (t )  u1 (t )  u2 (t )  vC (0, 002)  12, 0V

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 21
Dominio del tiempo: Ejemplo 3 – Diferentes frecuencias

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 22
CAPÍTULO 1: Conceptos Generales

1.3. Dominio de la Frecuencia

GICEP
Grupo de Investigación en Calidad de la Energía y Electrónica de potencia

150

10
0

Profesor: Armando Jaime Ustariz Farfán

0

E-Mail: ajustarizf@[Link]
-
5
0

-
10
0

-
150

0 2 4 6 8 1 12 14
0
t
[
m
s
]

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz).
 Dominio de la Frecuencia: Ejemplo 4 – RLC con fuente AC

Fuente AC:
1
U1 ( ) jC U 2 ( )
 I1 ( )   j L  jC  jC  U1 ( ) 
1

 I ( )    
 
 2    j C   2
U ( ) 
1
R j  C
I1 ( ) j L R I 2 ( )


Ybus   
i1 (t )  10sin t  5  ; i2 (t )  5sin t  45  ; f  60 Hz

I1 ( ) 105; I 2 ( )  545;   2 f 
U ( )  Ybus    I ( )
R  2; L  1mH y C  10uF 1

U1 ( )  3, 7795, 2
Solución:
U 2 ( )  9,9844, 7

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 24
Dominio de la Frecuencia: Ejemplo 4 – RLC con fuente AC

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 25
 Dominio de la Frecuencia: Ejemplo 5 – Fuente AC (F1≠F2)

1
U1 ( ) jC U 2 ( )
i1 (t )  10sin t  5  ; i2 (t )  5sin  5t  45  ; f  60 Hz

I1 ( ) 105; I 2 (5 )  545;   2 f

I1 ( ) j L R I 2 ( )
R  2; L  1mH y C  10uF

 j1 L  jC  j C 
Ybus      
  j C 
1
R j  C 
 I  k    Ybus  k   U  k  
 j 51 L  j 5C  j 5C 
Ybus  5     
  j 5C 
1
R j 5 C 

1
U1 ( )  3, 7895, 0; U1 (5 )  0,37  137, 2
U  k    Ybus  k    I  k   Solución:
U 2 ( )  0, 03  175, 4; U 2 (5 )  9,9942,8

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 26
Dominio de la Frecuencia: Ejemplo 4 – RLC con fuente AC

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 27
CAPÍTULO 1: Conceptos Generales

1.4. Taller #1

GICEP
Grupo de Investigación en Calidad de la Energía y Electrónica de potencia

150

10
0

Profesor: Armando Jaime Ustariz Farfán

0

E-Mail: ajustarizf@[Link]
-
5
0

-
10
0

-
150

0 2 4 6 8 1 12 14
0
t
[
m
s
]

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz).
Taller #1: Método de Sustitución Numérica

Fuente Línea 1- 2  Capacitor Línea  2 - 3  Carga 2

RS u1 (t ) R1 L1 u2 (t ) R2 L2 u3 (t )

t  0,002

u (t )
RC

FUENTE : LINEA 1  2  CAPACITOR LINEA  2  3 CARGA

U m 13.2 kV R1  0.088 / km RC  5kW R2  0.088 / km P1  25kW
V  0 L1  0.273mH / km C  0, 47  F L2  0.273mH / km Q1  30kVAr
f  60 Hz
L(1 2)  100km L(23)  50km VN  13.2kV
RS  0.01

Sbase  100MVA

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 29
Taller #1: Método de Sustitución Numérica

Fuente Línea 1- 2  Capacitor Línea  2 - 3  Carga 2

RS u1 (t ) R1 L1 u2 (t ) R2 L2 u3 (t )

t  0,002

u (t )
RC

1. Implementar en Matlab (Tiempo y Frecuencia)

2. Implementar en Simulink (Tiempo)
3. Implementar en Matpower (Frecuencia)
4. Comparar los tres resultados

Copyright © 2015 Universidad Nacional de Colombia . All rights reserved (a.j. Ustariz). 30