UNIVERSIDADE FEDERAL DO VALE DO SÃO FRANCISCO – UNIVASF

COLEGIADO DE ENGENHARIA ELÉTRICA

Disciplina: Sistemas Elétricos I
Prof.: Dr. Adeon Cecilio Pinto

Lista de Exercícios – Segunda Parte

1. Dado o sistema abaixo, assim como os dados. Pede-se:

Determinar as tensões complexas (módulo e ângulo) em todas as barras e os fluxos de
potências complexas (potências ativa e reativa) em todas as linhas.
Obs.: Escolha o método matemático.
barra 1 barra 2
0,1 + j 0,5 pu

barra 3
0,05 + j 0,2 pu
0,03 + j 0,4 pu

0,5 + j 2,0 pu

Especificação das barras (valores em pu)

BARRA P Q V  (°)
1 - - 1,1 0
2 0,8 - 1,0 -
3 -0,5 -0,2 - -

2. Dado o sistema elétrico de 4 barras abaixo, no qual são fornecidos os dados de barras e
admitâncias das linhas, em pu. Obtenha as tensões em todas as barras utilizando o método de
Newton-Raphson, considerando que o critério de parada seja P  0,001 e Q  0,001 .
 p2 = −0,5 pu
q3 = −0, 2 pu
p1 = 0, 6 pu
v1 = 1, 05 pu
Barra 1 0,5 − j 5 pu Barra 2
PV PQ

1 − j 3 pu
0,5 − j 5 pu 1 − j 3 pu

0, 2 − j 3 pu
Barra 4 Barra 3
V v4 = 1 pu PQ
 4 = 0
p3 = −0, 6 pu
q3 = −0,1 pu

Resposta: v2( ) = 0,9929 pu ; v3( ) = 0,9898 pu ; 1( ) = −0,0388 rad ; 2( ) = −0,1244 rad
2 2 2 2
e

3( 2) = −0,1138 rad .

3. A ilustração abaixo representa um sistema elétrico radial constituído por uma unidade geradora,
um sistema de transmissão, uma subestação de distribuição abaixadora (transformador) e duas
cargas, sendo uma dela alimentada em nível de transmissão (Carga 2) e outra em nível de
distribuição (Carga 3). Vale destacar que os dados do sistema são:
- Linha de transmissão: X = 0,05  km e comprimento = 100 km ;

- Transformador: Sn = 20 MVA , x = 5% e 138 kV 13,8 kV ;

- Carga 2: P2 = 10 MW e Q2 = 0 Mvar ;

- Carga 3: P3 = 5 MW e Q3 = 5 Mvar ;

Determine as tensões nas barras 2 e 3 através do fluxo de potência utilizando o Método

Desacoplado Rápido. Sabendo-se que a tensão na barra de referência (Barra 1) é mantida em
V1 = 10 pu .
 Obs.: 1) Utilizar os valores bases na barra 2 da tensão e da potência iguais a 138 kV e
10 MVA , respectivamente;
2) Utilizar quatro algarismos significativos.
Barra 1 Barra 2 Barra 3

Linha de transmissão
Gerador Subestação

Carga 3
Carga 2

4. Dado o simples de potência elementar, calcule a tensão na barra 1, aplicando o Método de

Gauss, utilizando a estimativa inicial da tensão na barra 1 igual a v1( ) = 0,90 pu . Adote o
0

nível de tolerância V =  = 0,005 pu .

V2  2 V11
Barra 2 Barra 1
E R + jX

Linha de transmissão
Gerador

Carga
P1 + jQ1

Dados do sistema:
Tensão nominal de linha: 69 kV ;

V2 = 10 pu ;

= 30 km ;

X ` = 0,555  km ;

R` = 0 ;
P1 = 10 MW ;

Q1 = 0 .

Obs.: Como sugestão adote os seguintes valores bases: Sbase = 10 MVA e Vbase = 69 kV .

Resposta: após 2 iterações obtém-se v1( ) = 0,9993− 0,035 rad pu , portanto

2

v1 = 0,0007 pu e 1 = 0,00348

5. Determine as grandezas faltosas através do estudo de fluxo de potência para o sistema
apresentado abaixo usando o Método Desacoplado, considerando um nível de tolerância
 P =  Q = 0, 001 pu .

Barra PQ Barra V  Barra PV

1 2 3
Z12 Z 23

sh sh sh sh
jB12 jB12 jB23 jB23
S1
S3
S2
jB1sh

Dados das barras

Barra
Tipo Vi esp [pu]  iesp [rad] Pi esp [pu] Qiesp [pu] Bish [pu]
i=
1 PQ - - -0,15 0,05 0,05
2 V 1,00 0,0 - - -
3 PV 1,00 - 0,20 - -

Dados das linhas

Barras
i= k= Z ik [pu] Biksh [pu]
1 2 0,03+j0,3 0,02
2 3 0,05+j0,8 0,01

Resposta: após 2 iterações obtém-se v1( ) = 1,0307− 0,0474 rad pu e 3(2) = 0,1605 rad
2

6. Repita o problema anterior empregando o Método Desacoplado Rápido, considerando um

nível de tolerância  P =  Q = 0, 001 pu .

Resposta: após 2 iterações obtém-se v1( ) = 1,0307− 0,0473 rad pu e 3(2) = 0,1605 rad
2

7. Determine os fluxos de potências ativa e reativa em todas as linhas do sistema abaixo, sem
utilizar nenhum método matemático, sabendo-se que v1 = v2 = v3 = 1,0 pu .
 v1 = 1, 0 pu
Barra 1 j 0,10 pu Barra 2
v2 = 1, 0 pu

j 0,10 pu j 0,10 pu

Barra 3
v3 = 1, 0 pu

Barra Demanda Demanda de Geração Geração de

i= de Ativos Reativos de Ativos Reativos
1 1,0 0,5 1,0 -
2 0,0 0,0 1,0 -
3 1,0 1,0 0,0 -
 FORMULÁRIO

Si = ( PGi − PDi ) + j ( QGi − QDi ) ; assim Pi = PGi − PDi e Qi = QGi − QDi

Pi − jQi − Yi1V1Vi* − Yi2 V2 Vi* − ... − Yin Vn Vi* ou ainda Pi − jQi = Yi1V1Vi* +Yi2 V2 Vi* + ... + Yin Vn Vi*

Pik =
1
R + X2 2 ( R V − R V V cos  + X V V sen  ) 
i
2
i k i k

Qik = 2
1
R + X2 ( X V − X V V cos  − R V V sen  )  para
i
2
i k i k
 = 1 −  2
Pki = 2
1
R + X2 ( R V − R V V cos  − X V V sen  ) 
k
2
i k i k

( )
1 2

Qki = 2 X V − X V V cos  + R V V sen 
R + X2
k i k i k

Ou da seguinte forma:

Pik = (1
R +X
2 2
R Vi − R Vi Vk cos (i −  k ) + X Vi Vk sen (i −  k )
2
)
R +X ( )
1
− X Vi Vk cos (i −  k ) − R Vi Vk sen (i −  k )
2
Qik = XV
2 2 i

Y  -Y  -Y  
11 12 1n
  
Y  -Y
=  21
 

Y22 
 - Y2n 
 

, G = Re Y   e B = ImY bus 
 
bus bus
 
    
- Yn1 - Y n2  Y nn 

1  Pi − jQi 
E(i  ( v )* −  k =1;k i Yik E k 
v +1) n
= (v)

Yii  Ei 
Pi = kni Vi Vk (Gik cos(θik ) + Bik sen(θik )) , ou
n

Pi = Pi − Vi  kni Vk ( G ik cos(θik ) + Bik sen(θik ) ) = 0 , ou ainda

n

gip = Pi − kni Vi Vk (Gik cos(θik ) + Bik sen(θik )) = 0 (estas três últimas equações são idênticas).
n

Qi = kni Vi Vk (Gik sen(θik ) − Bik cos(θik )) ou

n

Qi = Qi − Vi  k =1 Vk ( G ik sen(θik ) − Bik cos(θik ) ) = 0 , ou ainda

n

giq = Qi − kni Vi Vk (Gik sen(θik ) − Bik cos(θik )) = 0 (estas três últimas equações são idênticas).
n

 H   N´ gip P g p P g q Q g q Q
J =  J sendo: H = = ; N´= = ; J= = e L´= =
   L´  θ θ V V θ  V V

  H   N´  ΔP j 
−1
 Δθ j   ΔP j 
  −1      
− − − = −  J  * − − −  = − − − − − − −  * − − − 
 ΔV j   ΔQ j   J  L´   ΔQ j 

 θ j+1   θ j   Δθ j 
     
− − −  = − − −  + − − − 
 V j+1   V j   ΔV j 

 ΔP j    H   N´  Δθ j    H   N´.V   Δθ 
j

         −−− 
 − − −  = −  − − − − − − *
  − − −  = −  − − − − − − * 
 ΔQ 
j  J


 L´   ΔV´    J   L´.V  ( ΔV´ V ) 
j

 θ j+1   θ j   Δθ 
j

     −−− 
− − − = − − − +  
 V j+1   V j  ( ΔV j V ) 
 j
 

  H   N   ΔP j 
−1
 Δθ j 
     
 − − −  = − − − − − − −  * − − − 
 j
( ΔV V ) 
j  J  L   ΔQ j 


 ΔP j    H   N   Δθ 
j

     −−− 
 − − −  = −  − − − − − − * 
 ΔQ 
j  J
  L  ( ΔV V ) 
 j

 ΔP j    H  --   Δθ   H .Δθ j 
j
--
     −−−  
−−− 

− − − = − − − − − − − *   = − − − −
 ΔQ j   --
  L  ( ΔV V ) j   -- L. ( ΔV V ) 
j

1
Hik = Bik' =
Xik
1
H ii = Bii' = − 
Xik
Xik
Lik = Bik'' = =B
( R +Xik2 ) ik
2
ik

1
Lii = Bii'' = −  =Bii
Xik

ΔP V(v ) = −BI .Δθ(v )

ΔQ V(v ) = −BII .ΔV(v )

ΔPi  ε p e ΔQi  ε q