Apostila Curto Circuito - Arthur Damaso PDF

SISTEMAS ELTRICOS
DE POTNCIA
Notas de Aula
Prof. Alvaro Augusto W. de Almeida

Universidade Tecnolgica Federal do Paran UTFPR
Departamento Acadmico de Eletrotcnica DAELT
[email protected]
Verso 13/08/2014
Sistemas Eltricos de Potncia Notas de Aula Verso 13/08/2014 12:05 2
SUMRIO
1. GLOSSRIOS ............................................................................................................................................... 6
2. INTRODUO ............................................................................................................................................. 8
4. O SISTEMA POR UNIDADE .......................................................................................................................11
4.1. INTRODUO .............................................................................................................................................. 11
4.2. DEFINIO DE PU ....................................................................................................................................... 11
4.3. MUDANA DE BASE ..................................................................................................................................... 15
4.4. TRANSFORMADOR DE DOIS ENROLAMENTOS .................................................................................................. 17
4.5. TRANSFORMADOR DE TRS ENROLAMENTOS .................................................................................................. 21
4.6. TRANSFORMADOR COM TAP FORA DO VALOR NOMINAL ................................................................................... 23
4.7. MODELOS DE GERADORES SNCRONOS ........................................................................................................... 32
4.8. MODELOS DE LINHAS DE TRANSMISSO ......................................................................................................... 36
4.8.1. LINHA CURTA .............................................................................................................................................. 36
4.8.2. LINHA MDIA .............................................................................................................................................. 36
4.8.3. LINHA LONGA .............................................................................................................................................. 37
4.8 MODELOS DE CARGAS .................................................................................................................................. 38
4.9 INTRODUO AOS ESTUDOS DE CURTO-CIRCUITO. ........................................................................................... 42
4.10 EXERCCIOS ................................................................................................................................................ 49
5. COMPONENTES SIMTRICAS ..................................................................................................................54
5.1. INTRODUO .............................................................................................................................................. 54
5.2. O TEOREMA DE FORTESCUE ......................................................................................................................... 54
5.3. POTNCIA COMPLEXA .................................................................................................................................. 61
5.4. IMPEDNCIAS DE SEQUNCIA........................................................................................................................ 62
5.5. IMPEDNCIAS DE SEQUNCIA DOS COMPONENTES DE UM SEP.......................................................................... 64
5.4.1 LINHAS DE TRANSMISSO............................................................................................................................. 65
5.4.2 GERADORES SNCRONOS .............................................................................................................................. 66
5.4.3 TRANSFORMADORES DE DOIS ENROLAMENTOS ............................................................................................... 67
5.4.4 TRANSFORMADORES DE TRS ENROLAMENTOS .............................................................................................. 72
5.5 EXERCCIOS ................................................................................................................................................ 85
6. CLCULO DE CURTO-CIRCUITO ..............................................................................................................89
6.1. INTRODUO .............................................................................................................................................. 89
6.2. CURTO-CIRCUITO FASE-TERRA ..................................................................................................................... 90
6.3. CURTO-CIRCUITO FASE-FASE ........................................................................................................................ 92
6.4. CURTO-CIRCUITO FASE-FASE-TERRA ............................................................................................................. 94
6.5. MTODO DA MATRIZ IMPEDNCIA DE BARRA ............................................................................................... 102
6.7. EXERCCIOS .............................................................................................................................................. 109
7. FLUXO DE POTNCIA .............................................................................................................................112
7.1. INTRODUO ............................................................................................................................................ 112
7.2. MTODO DE GAUSS ................................................................................................................................... 115
7.3. MTODO DE GAUSS-SEIDEL ....................................................................................................................... 121
7.4. MTODO DE NEWTON-RAPHSON ................................................................................................................ 122
7.5. MTODO DESACOPLADO RPIDO ................................................................................................................ 137
7.6. EXERCCIOS .............................................................................................................................................. 141
8. ESTABILIDADE ESTTICA E TRANSITRIA .........................................................................................142
8.1 ESTABILIDADE EM REGIME PERMANENTE .................................................................................................... 142
8.2 MQUINA DE POLOS LISOS EM REGIME PERMANENTE .................................................................................... 143
8.3 CURVA DE CAPABILIDADE E CURVAS V ...................................................................................................... 149
Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

8.4 MQUINA DE POLOS SALIENTES EM REGIME PERMANENTE ............................................................................ 151

8.5 DINMICA DA MQUINA SNCRONA LIGADA AO BARRAMENTO INFINITO .......................................................... 159
8.6 EQUAO GERAL DA OSCILAO.................................................................................................................. 159
8.7 ANLISE LINEARIZADA MQUINA DE POLOS LISOS ..................................................................................... 162
8.8 MTODO DAS REAS .................................................................................................................................. 170
8.8.1 TOMADA SBITA DE CARGA.......................................................................................................................... 172
8.8.2 NGULO CRTICO PARA UM GERADOR E UMA LINHA ......................................................................................... 175
8.8.3 NGULO CRTICO PARA UM GERADOR E DUAS LINHAS....................................................................................... 179
9. ESTABILIDADE EM SISTEMAS MULTIMQUINAS ...............................................................................189
9.1 INTRODUO ............................................................................................................................................ 189
10. ESTABILIDADE TRANSITRIA NO DOMNIO DO TEMPO ..............................................................190
10.3 INTRODUO ............................................................................................................................................ 190
11. ESTABILIDADE DE TENSO ..............................................................................................................191
11.1 INTRODUO ............................................................................................................................................ 191
12. OPERAO ECONMICA DE SISTEMAS ...........................................................................................191
12.1 INTRODUO ............................................................................................................................................ 191
13. REFERNCIAS .....................................................................................................................................192

UM POUCO DE HISTRIA
O curso de Engenharia Eltrica da UTFPR, anteriormente denominado Engenharia Eltrica,
nfase Eletrotcnica, tem no momento mais de trinta anos de existncia e uma posio consoli-
dada junto comunidade acadmica paranaense. Nossa Universidade Tecnolgica Federal do
Paran completou seu primeiro centenrio em 2009 e, embora tenha ainda pouca experincia
como universidade, posio que assumiu apenas em 2005, tambm uma instituio consolidada
no Paran e no Brasil.
Durante toda a existncia do nosso curso de Engenharia, a disciplina de Sistemas Eltri-

cos de Potncia tem sido obrigatria. Anteriormente, nas grades antigas, este curso era oferecido
em apenas um semestre, com 60 horas semanais, no qual se abordavam basicamente curto-
circuito e fluxo de potncia, assim como os conceitos tericos necessrios ao desenvolvimento
de tais contedos. Assuntos como estabilidade de sistemas eram abordados de maneira introdu-
tria na disciplina de Gerao de Energia e havia pouco espao para estudos sobre transitrios,
fluxo de potncia timo e outros.
A disciplina Sistemas de Potncia 1 da grade atual tem basicamente a mesma ementa da

antiga Sistemas Eltricos de Potncia: modelagem de sistemas de potncia, sistema por unida-
de, componentes simtricas, curto-circuito e fluxo de potncia. J a disciplina de Sistemas de
Potncia 2 trata de controle de potncia ativa, reativa, tenso e frequncia, estabilidade esttica e
transitria e mtodos de anlise do problema da estabilidade. Adicionalmente, as disciplinas de
Planejamento de Sistemas, Proteo de Sistemas, Linhas de Transmisso, Subestaes, Comerci-
alizao de Energia Eltrica, dentre outras, possibilitam que o estudante possa concentrar seus
estudos, caso este seja seu objetivo, na rea de Sistemas de Potncia.
Longe de pretenderem substituir a literatura existente na rea, as presentes Notas de Aula

tm o objetivo de facilitar a introduo do aluno nessa rea fascinante dos Sistemas Eltricos de
Potncia, to essencial a um pas dotado de um imenso sistema eltrico interligado, como o
caso do Brasil.
Prof. Alvaro Augusto W. de Almeida

UTFPR, Curitiba, 2012

Observaes:
1) Todas as figuras deste trabalho, exceto a Figura 3.1, foram elaboradas pelo autor usando o
GNU Image Manipulation Program GIMP 2.6, disponvel em www.gimp.org.
2) Todas as fotografias so de domnio comum.
3) No momento (13/08/14) este um trabalho em progresso. Em caso de constatao de erros, o
autor agradece notificaes enviadas pelo e-mail [email protected].

1. GLOSSRIOS
Glossrio de smbolos usados como subndices ou sobre-ndices
Smbolo Indicao dada pelo ndice
0 Componente de sequncia zero
1 Componente de sequncia positiva
2 Componente de sequncia negativa
012 Sistema de sequncia (equilibrado)
a Fase a
b Fase b, valor base
c Fase c, perdas no ncleo (core)
abc Sistema original (desequilibrado)
ca Circuito aberto
cc Curto-circuito
d Componente de eixo direto
ef Valor eficaz
elt Grandeza eltrica
g Entreferro, componente de entreferro
h Ordem de um harmnico
i Entrada (input)
q Componente de eixo em quadratura
Tenso ou corrente de linha
max Valor mximo
mec Grandeza mecnica
mit Mquina de induo trifsica
mst Mquina sncrona trifsica
mim Mquina de induo monofsica
min Valor mnimo
msm Mquina sncrona monofsica
mdc Mquina de corrente contnua
m Grandeza magntica, magnetizao
n Valor nominal
n Componente normal
o Sada (output)
pu Por unidade (valor por unidade)
q Componente de eixo em quadratura
r Componente radial, rotor
rb Rotor bloqueado
s Saturado, sncrono, sncrona
T Total
Componente tangencial
Perdas hmicas, componente de perdas hmicas

Glossrio de smbolos gerais

Smbolo Unidade Descrio
2
a m/s Acelerao
A m2 rea da seo reta
B T Induo magntica
C F Capacitncia
e V Fora eletromotriz instantnea
E V Fora eletromotriz eficaz
f Hz Frequncia
fp Fator de potncia
F A-e/m Fora magnetomotriz
H A-e Intensidade magntica
I A Corrente eltrica
m Comprimento
L H Indutncia
N rpm Rotao, velocidade
Ns rpm Velocidade sncrona
p Nmero de polos
P W Potncia ativa
q Nmero de fases
Q var Potncia reativa
r Resistncia eltrica
r m Raio
s Escorregamento
S VA Potncia aparente
t S Tempo, intervalo de tempo
T Nm Torque, conjugado ou binrio
V V Tenso nos terminais
x Reatncia
xL Reatncia indutiva
xC Reatncia capacitiva
Z Impedncia
Graus, rad ngulo de carga
Graus, rad ngulo do fator de potncia
Wb Fluxo magntico por polo
Rendimento, eficincia
H/m Permeabilidade magntica
m Resistividade eltrica
Rad/s Velocidade angular ou frequncia angular
s Rad/s Velocidade angular sncrona


2. INTRODUO
Em agosto de 2010 o Sistema Interligado Brasileiro (SIN) era composto por 2.271 empreendi-
mentos de gerao em operao, totalizando uma potncia instalada de 110.224 MW. Desta po-
tncia, 69,24% correspondem a Usinas Hidreltricas (UHEs), 25,15% correspondem a Usinas
Termeltricas Convencionais (UTEs), 1,82% a Usinas Termeltricas Nucleares (UTNs) e o res-
tante a Pequenas Centrais Hidreltricas (PCHs), e Centrais Elicas (EOL).
A necessidade de alimentar os grandes centros consumidores do Sudeste a partir da ener-

gia produzida em regies remotas do pas tornou necessria a construo de uma extensa rede de
transmisso, conforme ilustrada na Figura 3.1.
Figura 3.1
Integrao eletroenergtica do Sistema Interligado Nacional (SIN)
Fonte: ONS, http://www.ons.com.br/conheca_sistema/mapas_sin.aspx
A interligao do SIN feita por meio da Rede Bsica, redefinida em 1998 por meio da
Resoluo ANEEL 245/1998. A Rede Bsica dos sistemas eltricos interligados constituda por
todas as linhas de transmisso em tenses iguais ou superiores a 230 kV e subestaes que con-

tenham equipamentos em tenso igual ou superior a 230 kV, integrantes de concesses de servi-
os pblicos de energia eltrica. Excepcionalmente, linhas e subestaes em tenses inferiores a
230 kV podem fazer parte da Rede Bsica, desde que autorizadas pelo Operador Nacional do
Sistema (ONS).
O planejamento e operao de sistemas eltricos, do porte do sistema brasileiro ou no,

demandam estudos bastante complexos, tais como os estudos de previso de carga, de fluxo de
potncia, de curto-circuito e de estabilidade. Modernamente todos esses estudos so feitos por
meio de computadores, alguns deles em tempo real. A finalidade das disciplinas de Sistemas de
Potncia 1 e 2 fornecer uma introduo a tais assuntos.
Fluxo de potncia, tambm conhecido como fluxo de carga, um problema matemti-

co, cujo objetivo determinar as tenses e potncias em todos os barramentos de um sistema
eltrico. Dessa forma podemos dimensionar linhas de transmisso, transformadores e demais
equipamentos que faro parte do sistema, bem como oper-los corretamente, de modo a manter
os padres adequados de tenso e frequncia. Um mtodo elementar para soluo de fluxo de
potncia o Gauss-Seidel, que, embora didtico, apresenta as desvantagens de no convergir
sempre e de no se prestar a sistemas com mais do que algumas dezenas de barras. Nesses casos,
outros mtodos, como o Newton-Raphson, devem ser utilizados. O mtodo desacoplado rpi-
do, que uma simplificao do Newton-Raphson, tambm pode ser utilizado em alguns casos.
A operao correta dos sistemas tambm depende do conhecimento dos nveis de curto-
circuito em cada barramento, de modo que sistemas adequados de proteo possam ser dimensi-
onados. Em linhas gerais, o problema de curto-circuito nada mais do que um problema de fluxo
de potncia no qual uma das barras submetida a condies de curto, ou seja, forada a manter
tenso nula ou quase nula. O curto, mais apropriadamente denominado falta, pode ser simtrico,
como nos casos dos curtos trifsico e trifsico-terra, ou assimtrico, como nos casos dos curtos
fase-terra, fase-fase ou fase-fase-terra. Sendo um problema de fluxo de potncia em condies
excepcionais, poderamos em princpio usar os mtodos de fluxo de potncia para resolver pro-
blemas de curto-circuito. Conduto, no caso dos curtos assimtricos, o problema se torna mais
complexo, pois as correntes em cada uma das fases sero diferentes. Felizmente, em situaes de
curto podemos fazer algumas simplificaes no sistema e podemos tambm usar o mtodo das
componentes simtricas, o que nos permitir conhecer correntes e potncias de curto em cada
uma das barras do sistema.

Fluxo de potncia e curto-circuito formam o contedo bsico de Sistemas de Potncia 1,

juntamente com os contedos auxiliares (modelos de equipamentos, sistema por unidade e com-
ponentes simtricas).
Finalmente, estudos de estabilidade tm o objetivo de determinar os limites operacionais

de geradores sncronos operando em sistemas multimquinas. Tais estudos podem ser divididos
em estabilidade angular, estabilidade em frequncia e estabilidade de tenso e, na atual gra-
de do nosso curso, so abordados na disciplina de Sistemas de Potncia 2.

4. O SISTEMA POR UNIDADE

4.1. Introduo
Sistemas Eltricos de Potncia (SEPs) so geralmente formados por vrios transformadores ele-
vadores e abaixadores. Em decorrncia disso, haver vrias tenses e correntes nominais em ca-
da lado de cada um dos transformadores, tornando os clculos bastante trabalhosos e complexos.
Assim, em vez de usarmos as unidades convencionais, como volts, amperes e ohms, conveni-
ente introduzirmos um sistema de unidades, denominado sistema pu (por unidade), no qual,
como veremos, todas as relaes de transformao de todos os transformadores se tornam unit-
rias,facilitando enormemente os clculos.
4.2. Definio de PU
Um valor em pu nada mais do que o valor original de uma grandeza qualquer, tal como tenso,
corrente, impedncia, etc., escrito em relao a um valor base da mesma grandeza. Sendo Vreal o
valor da grandeza original e Vbase o valor base, o valor expresso em pu ser
Vreal Definio de um valor

V pu . (4.1) em pu.
Vbase
Um valor expresso em pu igual a um centsimo do mesmo valor, quando expresso de

forma percentual. Da mesma forma que percentuais, valores em pu so adimensionais. Todavia,
costumamos anexar a partcula pu ao final dos valores, de modo a evitar confuso.
Quando expressamos valores finais, tanto faz usar pu ou %. Nos clculos, contudo, o sis-
tema pu mais adequado. A razo que dois valores percentuais, quando multiplicados, devem
ser divididos por 100 para resultar em um novo valor percentual. Por outro lado, a multiplicao
de dois valores em pu j fornece o novo valor tambm em pu.
____

Exemplo 4.1. Um transformador de tenso nominal primria igual a 13,8 kV opera momentane-
amente em 14 kV. Expresse a tenso de operao em pu e em percentual, na base do equipamen-
to.
Soluo. O valor em pu
Vreal 14 kV
V pu 1,0145 pu ,
Vbase 13,8 kV
enquanto o valor percentual
V % 100 V pu 100 1,0145 101,45% .
As principais vantagens do sistema por unidade so:
1) Os fabricantes de equipamentos tais como geradores, motores e transformadores costu-

mam fornecer reatncias e impedncias j em pu ou em %, expressas nas bases nominais
dos equipamentos.
2) Equipamentos semelhantes (mesma tenso, mesma potncia, etc.) tm impedncias seme-

lhantes quando expressas em pu. Isso facilita os clculos para substituio de equipamen-
tos e para expanso e reformulao de redes.
3) O uso do fator 3 minimizado nos clculos trifsicos em pu.
4) Como veremos, a impedncia de transformadores, quando expressa em pu, independen-

te do lado (alta, mdia, baixa tenso) que tomamos como referncia. Alm disso, a impe-
dncia de transformadores torna-se independente do tipo de ligao (delta-estrela, delta-
delta, estrela-estrela, etc.).
5) Em pu mais fcil identificar quando os valores de grandezas como tenses e potncias

se afastam dos valores nominais. Por exemplo, as tenses em qualquer barramento podem
variar em 5% em relao tenso nominal. Logo, as tenses mnima e mxima permiti-
das sero respectivamente iguais a 0,95 pu e 1,05 pu em relao tenso nominal, seja
qual for esta.
6) Caso a tenso seja 1 pu, a potncia aparente e a corrente em pu sero numericamente
iguais, por causa do cancelamento do fator 3 , como segue
3VI
S pu V pu I pu .
3Vb I b

Em princpio, h um grande grau de arbitrariedade na escolha do valor base para determi-

nada grandeza. Em sistemas de potncia, entretanto, estamos geralmente mais interessados em
quatro grandezas inter-relacionadas, o que far com que as respectivas bases sejam tambm in-
ter-relacionadas. So elas:
1) Tenso eltrica V (em kV).

2) Potncia aparente S (em MVA).
3) Corrente eltrica I (em A ou kA).
4) Impedncia Z (em ).
Escolhendo-se as bases para duas das grandezas acima, as bases para as outras duas se-
guem diretamente.
Geralmente iremos escolher as bases para tenso (Vb) e para potncia (Sb), calculando as
bases para impedncia (Zb) e corrente (Ib). Em circuitos trifsicos, que o caso usual, teremos
Zb
Vb 2 . (4.2)
Impedncia-base em
funo de Vb e Sb.
Sb
Sb Corrente-base em fun-
Ib . (4.3) o de Vb e Sb.
3Vb
Observaes:
1) A potncia-base nica e uma s para todos os barramentos do sistema em anlise.
2) As bases de tenso, corrente e impedncia transformam-se de acordo com as relaes de

transformao usuais dos transformadores.
3) Linhas de transmisso e impedncias em srie e em paralelo no afetam as bases de ten-

so, corrente e impedncia. Apenas transformadores afetam tais bases.
O exemplo a seguir esclarece essas caractersticas das bases das diversas grandezas.

Exemplo 4.2. Converta para pu as impedncias do sistema abaixo e determine as bases de tenso
e de impedncia em cada barramento. Considere que a potncia-base 20 MVA e que a tenso-
base no primeiro barramento 13,8 kV.
Figura 4.1
Sistema para o Exemplo 4.2
Soluo. A tenso-base na barra 1 Vb1 13,8 kV . A tenso-base na barra 2 pode ser obtida
considerando-se a relao de transformao do transformador, ou seja
138 kV
Vb2 kT12 Vb1 13,8 Vb2 138 kV
13,8 kV
A tenso-base na barra 3 igual tenso-base na barra 2, pois linhas de transmisso no
afetam as bases de tenso:
Vb3 138 kV
As impedncias-base podem ser obtidas a partir da potncia-base e das tenses-base
Z b1
V 13,8 kV
b1
2 2
Zb1 9,522
Sb 20 MVA
Z b2
V 138 kV
b2
2 2
Zb2 952,2
Sb 20 MVA
Z b3 Zb2 Zb3 952,2
As reatncias do gerador e do transformador podem ser facilmente convertidas para pu
j10%
xG1 xG1 j0,10 pu
100
j12%
xT12 xT12 j 0,12 pu
100

A reatncia em pu da linha de transmisso pode ser obtida dividindo-se a reatncia em

ohms pela impedncia-base nas barras 2 e 3
j80
xLT12 xLT12 j 0,084 pu
952,2
4.3. Mudana de base

As impedncias de equipamentos tais como geradores, motores e transformadores so geralmen-
te expressas pelo fabricante nas respectivas bases nominais. Contudo, as bases do sistema em
anlise geralmente so diferentes das bases dos equipamentos, sendo necessrio transformar de
uma para outra e vice-versa. Sejam inicialmente as variveis abaixo:
Sbv = potncia-base velha (equipamento).
Sbn = potncia-base nova (sistema).
Vbv = tenso-base velha (equipamento).
Vbn = tenso-base nova (sistema).
Z = impedncia original do equipamento, em ohms.
Z vpu = impedncia em pu na base velha.
Z npu = impedncia em pu na base nova.
Retomando a definio de pu, podemos agora escrever
Z Z
Z vpu e Zn
pu
.
Z bv Z bn
Igualando Z nas expresses acima, vm
Z Zvpu Zbv Z npu Zbn .

Queremos obter a impedncia em pu na base nova em funo da impedncia em pu na
base antiga. Logo, devemos escrever
Z bv
Z npu Z vpu .
Z bn
Substituindo Z bv Vbv / Sbv e Z bn Vbn / Sbn , teremos

2 2

2 Mudana de bases de
S V
Z n
pu
Z bn bv .
v
pu
(4.4) uma impedncia em pu.
Sbv Vbn
____
Exemplo 4.3. Considerando, no sistema abaixo, que a potncia-base 50 MVA e que a tenso-
base na barra 1 15 kV, converta todas as impedncias para pu, nas bases do sistema.
Figura 4.2
Soluo. A tenso-base na barra 1, Vb1 , foi arbitrada em 15 kV. A tenso-base na barra 2 pode
ser obtida a partir de Vb1 , ou seja
125 kV
Vb2 15 Vb2 135,87 kV
13,8 kV
A tenso-base na barra 3 igual tenso-base na barra 2
Vb3 135,87 kV
A tenso-base na barra 4 pode ser calculada da mesma maneira
6,6 kV
Vb4 135,87 Vb4 6,50 kV
138 kV
A nica impedncia-base que interessa a das barras 2 e 3, pois somente nesse trecho
temos impedncias em ohms que devem ser convertidas para pu
Z b2 Z b3
V 135,87 kV
b2
2 2
Zb2 369,21
Sb 50 MVA

As reatncias de G1, T12 e T34 podem agora ser expressas em pu e transformadas para as
bases novas (do sistema)
2
50 15
xG1 j 0,08 xG1 j 0,1333 pu
30 15
2
50 13,8
xT12 j 0,10 xT12 j 0,0846 pu
50 15
2
50 138
xT34 j 0,12 xT34 j 0,1547 pu
40 135,87
A reatncia da linha de transmisso pode finalmente ser calculada como
j100 j100 j100

xLT23 x 0,2709 pu
Zb2 Vb 2 / Sb 135,87 2 / 50 LT23
2
A carga na barra 4 tambm pode ser escrita em pu
25 MVA
S4 S4 0,50 pu
50 MVA
Sabendo agora que todos os elementos do sistema podem ser representados por meio de
suas impedncias, podemos desenhar o diagrama da Figura 4.3.
Figura 4.3
Diagrama de reatncias para o Exemplo 4.3
4.4. Transformador de dois enrolamentos

Podemos agora mostrar que a impedncia em pu de um transformador de dois enrolamentos a
mesma, independentemente do lado que se tome como, referncia. Considere inicialmente o mo-
delo de circuito equivalente de um transformador genrico de dois enrolamentos, como mostrado
na Figura 4.4, no qual os parmetros do secundrio foram referidos ao primrio por meio da re-
lao de espiras k=N1/N2.

Figura 4.4
Circuito equivalente por fase de um
transformador de dois enrolamentos
O circuito ilustrado para uma fase apenas, pois os circuitos para as demais fases so
idnticos, a menos das defasagens adequadas de tenses e correntes. Os parmetros do circuito
equivalente, em /fase, so:
r1 = resistncia eltrica do primrio.
k 2 r2 = resistncia eltrica do secundrio referida ao primrio.

x1 = reatncia de disperso do primrio.
k 2 x2 = reatncia de disperso do secundrio referida ao primrio.

rc = resistncia eltrica correspondente s perdas no ncleo (histerese e Foucault).
xm = reatncia de magnetizao.
O procedimento matemtico de se referir as impedncias do secundrio ao primrio per-

mite substituir o acoplamento magntico do transformador por um acoplamento eltrico, mais
fcil de ser tratado.
Em transformadores de potncia, que sempre o nosso caso, a corrente de excitao I
desprezvel frente corrente do primrio I1 . Sendo assim, e desde que o transformador esteja
prximo condio nominal, o ramo de excitao pode ser removido. O circuito equivalente
simplificado resultante mostrado na Figura 4.5.
Uma segunda simplificao possvel, pois transformadores de potncia so construdos

com condutores de seo reta elevada e, logo, de baixa resistncia eltrica. Assim, as resistncias
r1 e r2 podem ser desprezadas frente s reatncias x1 e x2 , resultando no circuito equivalente
mostrado no Figura 4.6.

O transformador de potncia pode assim ser representado por uma nica reatncia referi-
da ao primrio. Contudo, essa mesma reatncia pode tambm ser referida ao secundrio, resul-
tando em xT x1 / k 2 x2 .
As reatncias referidas ao primrio e ao secundrio sero tanto mais diferentes entre si

quando maior for o valor da relao de transformao k.
Figura 4.5
Circuito equivalente simplificado de
um transformador de dois enrolamentos
Figura 4.6
Circuito equivalente simplificado final de
um transformador de dois enrolamentos
A reatncia xT=x 1 + k2x2 pode ser obtida por meio do ensaio de curto-circuito, tambm
conhecido como ensaio de corrente nominal.
Podemos agora mostrar que, quando expressa em pu, xT independe de que lado tomamos
como referncia. Sejam inicialmente:
x A = reatncia prpria do lado de alta tenso.

xB = reatncia prpria do lado de baixa tenso.
xTA = reatncia total, referida ao lado de alta tenso.
xTB = reatncia total, referida ao lado de baixa tenso.
A reatncia total referida ao lado de alta ser

xTA x A k 2 xB ,
a qual, convertida para pu, poder ser escrita como

x A k 2 xB
x pu
TA .
Z bA

Sabendo ainda que Z bA VbA / Sb e k VbA / VbB , vem
2
xTpuA x Apu
V / V x
bA bB
2
B
V / S 2 ,
bA b
ou,
xB
xTpuA x Apu x Apu xBpu .
V / S 2
(4.5)
bB b
De forma semelhante, podemos escrever a reatncia referida ao lado de baixa como
xA
xTB xB ,
k2
a qual, convertida para pu, poder ser escrita como
xA
2
xB
xTpuB k .
Z bB
ou
pu

V / V x
bB bA
2
A
xBpu ,
V / S
xTB 2
bB b
ou, ainda
xA
xTpuB xBpu x Apu xBpu .
V / S 2
(4.6)
bA b
Comparando (4.5) e (4.6), vem que
Reatncia total, em pu,

xTpuA xTpuB x Apu xBpu . (4.7) de um transformador.
Sabendo que as reatncias de um transformador, em pu, so iguais, independente do lado

ao qual forem referidas, segue tambm que a relao de transformao k, em pu, unitria

Relao de tenses de
xTpuA
k pu
1. (4.8) um transformador, em
xTpuB pu.
Essa provavelmente a maior vantagem do uso do sistema pu, pois podemos tratar trans-
formadores como meras impedncias, sem nos preocuparmos com referncias a enrolamentos e
fatores de transformao.
4.5. Transformador de trs enrolamentos

Transformadores de trs enrolamentos so bastante comuns em sistemas de potncia e podem ser
representados em diagramas unifilares por meio do smbolo unifilar da Figura 4.7(a). Para fins
de clculos, contudo, deveremos adotar a representao da Figura 4.7(b), onde:
xam = reatncia de disperso entre os terminais de alta e de mdia tenso, com o ter-
minal de baixa tenso aberto.
xab
= reatncia de disperso entre os terminais de alta e de baixa tenso, com o ter-
minal de mdia tenso aberto.
xm b
= reatncia de disperso entre os terminais de mdia e de baixa tenso, com o
terminal de alta tenso aberto.
O modelo resultante uma espcie de delta, mas devemos salientar que h pouco em co-
mum entre este delta e as ligaes homnimas comuns em circuitos trifsicos. Assim, no po-
demos usar as transformaes Y estudadas em circuitos eltricos.

Figura 4.7
(a) Smbolo unifilar de um transformador
de trs enrolamentos; (b) modelo em delta
de um transformador de trs enrolamentos
Para facilitar os clculos e evitar a circulao de correntes fictcias, podemos converter o

modelo delta para um modelo estrela, conforme a Figura 4.8.
Figura 4.8
Modelo em estrela de um
transformador de trs enrolamentos
Tomando os enrolamentos aos pares, sempre com o terceiro a vazio, podemos escrever
xam xa xm . (4.9)
xab xa xb (4.10)
xmb xm xb (4.11)

Resolvendo o sistema acima, teremos
xa 1
2 xam xab xmb . (4.12) Reatncias de um mode-
lo Y para um transfor-
xb 1
2 xab xmb xam . (4.13)
mador de trs enrola-
mentos.
xm 1
2 xam xmb xab . (4.14)
4.6. Transformador com tap fora do valor nominal

Muitas vezes os transformadores operam fora da tenso nominal, por meio de taps (derivaes),
e, assim, precisamos desenvolver um modelo para esses casos. Iniciamos definindo uma varivel
auxiliar
A
a , (4.15)
B
onde
Tenso do lado de AT
A , (4.16)
Tenso nominal do lado de AT
Tenso do lado de BT
B , (4.17)
Tenso nominal do lado de BT
O transformador fora do tap nominal pode agora ser modelado como na Figura 4.9, ou se-
ja, um transformador ideal de relao a : 1 em srie como uma admitncia yT , que representa o
transformador quando operando no tap nominal.
Entre as barras a e r, que correspondem ao transformador ideal, podemos escrever
Sa Sr ,
ou,
Va Ia* Vr Ir* Vr Iab
*
,
ou, ainda,
V *
Va Ia* a Iab .
a
Finalmente,

I
Ia ab* . (4.18)
a
Figura 4.9
Modelo inicial para o transformador
com tap fora do valor nominal
Podemos tambm escrever a corrente Iab em funo das tenses nas barras r e b, ou seja
V
Iab a * Ia Vr Vb yT a Vb yT .
a
ou,
y y
Ia T2 Va T* Vb . (4.19)
a a
Da mesma forma, podemos escrever a seguinte relao para a corrente no lado de baixa
V
Ib Iab Vb Vr yT Vb a yT ,
a
ou,
y
Ib T Va yT Vb . (4.20)
a
Escrevendo (4.19) e (4.20) sob forma matricial, teremos
Ia yT / a 2 yT / a * Va
(4.21)
I b yT / a yT Vb

A equao (4.15) acima formalmente idntica equao matricial de um circuito equi-

valente , conforme mostrado na Figura 4.10, e que pode ser escrita como
Ia Yaa Yab Va
(4.22)
I b Yba Ybb Vb
Igualando (4.21) e (4.22), teremos
Yaa yT / a 2
Ybb yT
(4.23)
Y y / a *
ab T
Yba yT / a
Lembrando que uma das propriedades dos elementos da matriz admitncia nodal que
Yab Yba , segue-se que devemos ter a a , ou seja, a deve ser um numero real, o que significa
*
que, como sabemos, os taps do transformador apenas alteram o mdulo da tenso, mas no o
ngulo de fase. Note tambm que, na nossa notao, Y representa um elemento da matriz admi-
tncia nodal [Y ] , enquanto y representa uma admitncia fsica do circuito.
Figura 4.10
Modelo para o transformador
com tap fora do valor nominal
Escrevendo as equaes nodais para o sistema da Figura 4.10, teremos
Ia Va y a (Va Vb ) y ab .
Ib Vb yb (Vb Va ) y ab
ou,
Ia Va y a y ab Vb y ab .

Ib Va y ab Vb yb y ab
Escrevendo as equaes acima sob forma matricial, teremos
Ia y a y ab y ab Va Yaa Yab Va
(4.24)
I b yab y b y ab Vb Yba Ybb Vb
As regras de formao da matriz admitncia nodal podem ser escritas como:
Yaa y a y ab
Ybb y b y ab (4.25)
Y Y y
ab ba ab
Em resumo, os elementos Yii so iguais soma de todas as admitncias que se ligam ao
n i, enquanto a admitncia Yij Yji igual ao recproco da admitncia fsica que liga os ns i e
j.
Comparando as equaes (4.23) e (4.25), e considerando tambm que a a * a , tere-

mos
yT / a 2 y a y ab
yT y b y ab
(4.26)
yT / a y ab
yT / a y ab
Das equaes (4.26), segue-se que
yT 1 a Admitncias de um
y a . (4.27) transformador com tap
a2
fora do valor nominal.
yT a 1
yb (4.28)
a
yT
y ab (4.29)
a

Note que, se tivermos a=1, ou seja, se ambos os taps do transformador estiverem na ten-
so nominal, teremos y a yb 0 e y ab yT 1 / xT , e voltaremos ao modelo original de um
transformador de potncia de dois enrolamentos.
____
Exemplo 4.4. Para o sistema da Figura 4.11, pede-se: (a) considerando que a potncia-base 100
MVA e que a tenso-base 15 kV no barramento 1, converta os parmetros do sistema abaixo
para pu; (b) apresente os resultados em diagrama unifilar, na forma retangular.
Figura 4.11
Soluo. Primeiramente devemos calcular as tenses-base em cada um dos barramentos.

Vb1 15 kV . A tenso-base na barra 2 pode ser obtida a partir da relao de transformao do
transformador 1-2, que um elevador de tenso
138 kV
Vb2 kT12 Vb1 15 Vb2 138 kV
15 kV
Sabendo que no h queda de tenso-base em uma linha de transmisso, as tenses-base
nas barras 2 e 3 sero iguais
Vb3 138 kV
As tenses-base nas barras 4 e 5 so calculadas a partir das relaes de transformao do

transformador de trs enrolamentos 3-4-5, que um abaixador de tenso
69 kV
Vb4 kT34 Vb3 138 Vb2 41,4 kV
230 kV

13,8 kV
Vb5 kT35 Vb3 138 Vb2 8,28 kV
230 kV
Vb6 Vb4 Vb6 41,4 kV
Finalmente, a tenso-base na barra 7 decorre da relao de transformador do transforma-

dor abaixador 6-7
11 kV
Vb7 kT67 Vb6 41,4 Vb7 10,517 kV
3 25 kV
O clculo da reatncia do gerador 1 um caso de mudana de base. Aplicando a relao

(4.4), teremos
2
S V
Z n
pu
Z bn bv .
v
pu
Sbv Vbn
2
100 13,8
xG1 j 0,1 xG1 j 0,1058 pu
80 15
Da mesma forma, teremos a seguinte relao para o transformador 1-2
2
100 15
xT12 j 0,11 xT12 j 0,1222 pu
90 15
A reatncia da linha de transmisso 2-3 j est em pu, mas est expressa nas bases 230
kV e 50 MVA. Logo, devemos fazer uma mudana de bases
2
100 230
xLT23 j 0,03 xLT23 j 0,1667 pu
50 138
As reatncias do transformador 3-4-5 j esto nas tenses-base corretas, bastando mudar

as bases de potncia
2
100 230
xam j 0,13 xam j 0,4012 pu
90 138
2
100 230
xab j 0,15 xab j 0,8333 pu
50 138
2
100 69
xm b j 0,11 xmb j 0,3395 pu
90 41,4

As reatncias acima correspondem ao modelo delta da Figura 4.7(b). Devemos ento

convert-las para o modelo estrela da Figura 4.8
xa 1
2 xam xab xmb 1 2 j 0,4012 j 0,8333 j 0,3395 xa j 0,8950 pu
xb 1
2 xab xmb xam 1 2 j 0,8333 j 0,3395 j 0,4012 xb j 0,7716 pu
xm 1
2 xam xmb xab 1 2 j 0,4012 j 0,3395 j 0,8333 xm j 0,0926 pu
A reatncia da linha 4-6 est expressa em ohms. Para convert-la para pu devemos dividi-
la pela impedncia-base do trecho 4-6, ou seja
j 20 j 20
xLT46 xLT46 j1,1669 pu
(Vb 4 ) / Sb (41,4) 2 / 100
2
O clculo da reatncia do transformador 6-7 exige algum cuidado, pois se trata de um

banco trifsico com trs unidades monofsicas. Assim, as tenses dadas so de fase e a potncia
monofsica. Logo, teremos
2
100 25 3
xT67 j 0,08 xT67 j 0,2917 pu
3 10 41,4
Finalmente, a reatncia-base do gerador 7 ser

2
100 15
xG7 j 0,12 xG7 j1,2205 pu
20 10,517
Figura 4.12
Diagrama de reatncias para o Exemplo 4.4.
Todas as reatncias esto em pu.

Exemplo 4.5. Para o sistema da Figura 4.13, sabendo que a tenso na barra 5 1,0 pu e conside-
rando Sb=50 MVA e Vb1=13,8 kV, pede-se: (a) a corrente na barra 5, em pu e em amperes; (b) a
tenso na barra 1, em pu e em volts.
Soluo. Fazendo Vb1=13,8 kV, todas as tenses-base j so iguais s respectivas tenses nomi-
nais. Alm disso, as reatncias do gerador e dos transformadores j esto nas tenses-base corre-
tas. Basta reescrev-las para a nova potncia-base. Logo
2
50 13,8
xG1 j 0,1 xG1 j 0,0667 pu
75 13,8
2
50 13,8
xT12 j 0,08 xT12 j 0,0444 pu
90 13,8
2
50 138
xT34 j 0,12 xT34 j 0,1000 pu
60 138
Figura 4.13
As reatncias das linhas podem ser convertidas para pu dividindo-as pelas respectivas
impedncias-base
j 50 j 50
xLT23 xLT23 j 0,1313 pu
(Vb 2 ) / Sb (138) 2 / 50
2
j 20 j 20
xLT45 xLT45 j 0,2100 pu
(Vb 4 ) / Sb (69) 2 / 50
2
Devemos converter para pu tambm as potncias nas barras 4 e 5, dividindo-as pela po-
tncia-base
20
S4pu S4pu 0,4 pu
50
30
P5pu P5pu 0,6 pu
50

O diagrama unifilar simplificado resultante mostrado na Figura 4.14.
Figura 4.14
Diagrama de reatncias para o
Exemplo 4.5. Todos os valores em pu
A corrente na barra 5 pode ser obtida a partir da tenso e da potncia nessa barra, ou seja
P5pu V5pu I 5pu cos45
ou
1,0 I 5pu 0,9 0,6 I 5pu 0,667 pu
A corrente-base na barra 5
Sb 50 10 6
I b5 418,37 A
3Vb5 3 69 10 3
Logo, a corrente em amperes na barra 5 ser
I 5 I 5pu I b5 0,667 418,37 I 5 278,93 A
Para calcular a tenso na barra 1 devemos antes calcular a tenso na barra 4
V4pu V5pu jx45 I5pu 1,0 j0,21 0,667 25,84 V4pu 1,06856,7756 pu
A corrente na barra 4 pode ser obtida a partir da tenso e da potncia nessa barra, ou seja
S 4pu V4pu I 4pu
ou
1,0685 I 4pu 0,4 I 4pu 0,3744 pu
A corrente entre as barras 1 e 4 ser a soma de I4 e I5 , ou seja
I45pu I4pu I5pu 0,3744 18,195 0,667 25,84 I45pu 1,0393 23,093 pu

Finalmente, a tenso na barra 1 ser
V1 pu V4pu jx14 I45pu 1,06856,7756 j0,2757 1,0393 23,093
ou
V1 pu 1,06856,7756 0,286566,907 V1 pu 1,236118,33 pu
Sabendo que a tenso-base na barra 1 13,8 kV, a tenso em volts na barra 1 ser
V1 13,8 1,236118,33 V1 17,05818,33 kV
O Exemplo 4.5 ilustra um clculo elementar de fluxo de potncia, no qual desejamos

calcular a tenso e a potncia em cada um dos barramentos. A situao seria muito mais compli-
cada se, em vez de termos a tenso e a potncia na barra 5, desejando a tenso na barra 1, o in-
verso acontecesse, ou seja, se tivssemos a tenso na barra 1 e potncia na barra 5, desejando a
tenso na barra 5. Ao escrevermos as equaes do circuito, perceberamos que o sistema de
equaes resultantes seria no linear. Com o aumento do nmero de barras, a soluo analtica
do sistema seria muito difcil ou mesmo impossvel. Nesse caso, mtodos mais genricos e pode-
rosos devem ser desenvolvidos, como veremos no captulo 7.
4.7. Modelos de geradores sncronos

Um gerador sncrono composto por dois circuitos acoplados magneticamente. O primeiro a
armadura trifsica, localizada no estator e responsvel pela transferncia de potncia eltrica AC
entre a mquina e o sistema de potncia ao qual ela se conecta. O segundo circuito o campo,
localizado no rotor e alimentado com corrente contnua, de modo a produzir um fluxo magntico
constante. Sendo N f o nmero de espiras por fase da armadura, f1 a frequncia das correntes da
armadura, 2 o fluxo magntico por polo produzido pelo rotor, a fora eletromotriz E f induzi-
da em cada fase da armadura a vazio ser
Fora eletromotriz indu-

E f 2f1 N f 2 k w , (4.30) zida em cada fase de
uma armadura a vazio.
onde k1w , ainda, o fator de enrolamento da armadura, tipicamente maior do que 0,85 e menor
ou igual a 1,0.

Quando alimenta uma carga qualquer, de maneira isolada ou conectado ao sistema, a ten-
so nos terminais do gerador ser V1 E f , indicando a presena de uma impedncia interna,
usualmente representada em srie. Contudo, por causa do desacoplamento eltrico entre campo e
armadura, o gerador sncrono uma fonte de corrente quase ideal, podendo ser representado ini-
cialmente como na Figura 4.15, onde xm a reatncia de magnetizao, x1 a reatncia de dis-
perso da armadura, r1 a resistncia hmica da armadura e rc a resistncia de perdas no n-
cleo (histerese e Foucault). Todos os parmetros so expressos em ohms por fase.
Figura 4.15
Modelo inicial de um gerador
sncrono trifsico
possvel fazer algumas simplificaes no circuito da Figura 4.15. Nos geradores co-
muns em sistemas de potncia, sempre da classe MVA, os condutores da armadura tm bitola
larga a ponto da resistncia r1 ser desprezvel. As perdas no ncleo tambm so desprezveis, o
que significa que a resistncia rc muito grande em comparao com xm, e podemos fazer
rc // xm xm . O resultado o circuito da Figura 4.16, que consiste de um equivalente Norton em
srie com uma reatncia de disperso jx1.
Figura 4.16
Modelo intermedirio de
um gerador sncrono trifsico

Finalmente, o equivalente Norton pode ser convertido em um equivalente Thvenin, no

qual E f jxm I f e xd xm x1 denominada reatncia sncrona de eixo direto. O circuito
equivalente final, mostrado na Figura 4.17, adequado a geradores sncronos de polos lisos, que
geralmente o caso de turbogeradores. Para geradores de polos salientes, que geralmente o
caso de hidrogeradores, algumas modificaes devem ser introduzidas, as quais sero objeto do
captulo 9.
Figura 4.17
Modelo de circuito equivalente de
um gerador sncrono de polos lisos
Considerando que, em um gerador, o sentido da corrente de armadura I1 da mquina

para a carga, a equao fasorial correspondente pode ser escrita como
Equao fasorial de um
E f V1 jxd I1 . (4.31) gerador de polos lisos
em regime permanente.
A nica modificao necessria para transformar o gerador descrito pela equao (4.31)
em um motor sncrono a mudana do sentido da corrente, resultando na seguinte equao
Equao fasorial de um
E f V1 jxd I1 . (4.32) motor de polos lisos em
regime permanente.
As equaes (4.31) e (4.32) descrevem bastante bem o comportamento da mquina sn-

crona de polos lisos funcionando em regime permanente. No caso de geradores funcionando em
regime transitrio deveremos introduzir correes nas reatncias sncronas.

Vamos supor que um gerador sncrono esteja funcionando a vazio quando um curto-
circuito trifsico ocorre. Vamos supor tambm, por simplicidade, que o curto ocorre exatamente
quando a tenso alternada do gerador instantaneamente nula. Por causa do carter indutivo do
gerador, a corrente no atingir imediatamente um valor de regime constante, mas se comportar
como mostrado na Figura 4.18. A envoltria da senoide uma exponencial mais complexa do
que o usual, pois sua taxa de decaimento no constante. Para evitar a dificuldade de se traba-
lhar com uma quantidade muito grande de constantes de tempo, costumamos definir trs pero-
dos de tempo, cada um deles caracterizado por uma reatncia sncrona:
1) Perodo subtransitrio: corresponde aos primeiros ciclos aps o curto, durante os quais
a corrente decai muito rapidamente; caracterizado pela reatncia subtransitria de eixo
direto, xd ' ' .
2) Perodo transitrio: corresponde ao perodo aps o perodo subtransitrio e antes da cor-
rente ter se estabilizado, durante o qual a corrente decai mais lentamente; caracterizado
pela reatncia transitria de eixo direto, xd ' .
3) Perodo de regime permanente: corresponde ao perodo aps a corrente ter se estabili-
zado; caracterizado pela reatncia sncrona de eixo direto usual, xd .
Figura 4.18
Corrente de armadura de um gerador sncrono
em curto-circuito trifsico simtrico
A Tabela 4.1 mostra os valores tpicos das reatncias de algumas mquinas sncronas.
Note que a relao entre as reatncias sncrona xd e subtransitria xd ' ' pode chegar a 11 vezes
no caso do gerador de polos salientes. Como veremos no captulo 5, essa diferena torna bastante
crtica a escolha do perodo no qual devemos calcular as correntes de curto-circuito.

A corrente de curto da Figura 4.18, denominada corrente de curto simtrica, um caso

particular de um caso mais geral, o das correntes de curto assimtricas, as quais tm uma com-
ponente contnua que as desloca para cima ou para baixo. Uma corrente assimtrica corresponde
a uma corrente simtrica mais uma componente contnua que decai exponencialmente.
Tabela 4.1 Reatncias tpicas de mquinas sncronas
Reatncia Gerador de Gerador de Motor de

polos lisos polos salientes polos salientes
Sncrona, xd (pu) 1,10 1,10 1,10
Transitria, xd (pu) 0,20 0,35 0,50
Subtransitria, xd (pu) 0,10 0,23 0,35
4.8. Modelos de linhas de transmisso

Ao contrrio do que acontece com as redes de distribuio, as linhas de transmisso trifsicas,
quando em regime, operam geralmente de maneira equilibrada, o que permite a classificao de
tais equipamentos em trs tipos bsicos: linhas curtas, linhas mdias e linhas longas.
4.8.1. Linha curta

Linhas de transmisso curtas so aquelas de comprimento inferior a 80 km. Nesse caso, adota-
do um modelo simplificado que nada mais do que uma impedncia Z LT rLT jxLT por fase,
representado de maneira unifilar como na Figura 4.19. Neste modelo, rLT a resistncia hmica,
responsvel pelas perdas por efeito Joule, e xLT a reatncia indutiva da linha. Ambos os par-
metros so especificados em ohms por fase.
Figura 4.19
Modelo de uma linha de transmisso curta
4.8.2. Linha mdia

Linhas cujo comprimento superior a 80 km, mas inferior a 240 km so denominadas linhas
mdias. Nesse caso as capacitncias entre a linha e o terra no podem ser desprezadas e devere-
mos usar o modelo T, conforme representado na Figura 4.20, ou o modelo , conforme repre-

sentado na Figura 4.21. Em ambos o termo jBc ( jxc ) 1 representa a susceptncia total da
linha.
Figura 4.20
Modelo T de uma linha
de transmisso mdia
Note que a nica diferena entre os modelos e T uma distribuio diferente da impe-
dncia srie e da susceptncia paralela ao longo do trecho em questo. Quando a capacitncia em
paralelo for desprezvel, o que significa Bc , ambos os modelos se reduzem ao modelo de
linha curta. Daremos sempre preferncia ao modelo e, quando nada for mencionado, este o
modelo que deve ser usado.
Figura 4.21
Modelo de uma linha
de transmisso mdia
4.8.3. Linha longa

Linhas de comprimento superior a 240 km so consideradas longas, caso no qual o modelo com-
pleto da linha de transmisso deve ser usado. Neste modelo as impedncias srie e susceptncias
paralelas so consideradas uniformemente distribudas ao longo da linha. Considerando que z e
b so, respectivamente, a impedncia e a susceptncia por unidade de comprimento, e que l o

comprimento total da linha, podemos escrever equaes diferenciais parciais para a linha, as
quais, uma vez resolvidas, resultam em
z senh(l ) Parmetros de uma

Z eq , (4.33) linha de transmisso
longa.
2b tanh( l / 2)
B eq , (4.34)

onde o parmetro zb denominado constante de propagao. Depois de calculados, os
parmetros Z eq e B eq podem ser inseridos em um circuito equivalente T, como na Figura 4.20,
ou , como na Figura 4.21.
Em nossas simulaes as impedncias e susceptncias da linha de transmisso sempre se-

ro parmetros conhecidos. Assim, no faz muita diferena se o modelo a ser utilizado para
linha mdia ou linha longa. Caso a linha seja longa, simplesmente consideraremos que algum j
calculou Z eq e B eq para ns.
4.8 Modelos de cargas

Dentre os vrios parmetros de um SEP a carga dos consumidores a de determinao mais dif-
cil. Considerando que o valor da carga varia de segundo a segundo e que existem milhes de
consumidores, cada um absorvendo energia de acordo com sua exigncia individual, a determi-
nao das exigncias futuras um problema estatstico. A curva de carga de um dado barra-
mento de distribuio, ilustrada de forma genrica na Figura 4.22, decorre de hbitos de consu-
mo, temperatura, nvel de renda, forma de tarifao, etc.

Figura 4.22
Curva de carga tpica de um
barramento de distribuio
A carga total do sistema pode, grosso modo, ser repartida entre usurios industriais e re-
sidenciais. A potncia consumida pelos consumidores industriais varia de um tero nas horas de
pico at metade nas horas de carga mnima. Uma diferena muito importante entre os dois tipos
de consumidores que nos industriais existe uma porcentagem elevada de motores de induo
(cerca de 60 por cento), enquanto nos consumidores residenciais predominam as cargas de aque-
cimento e iluminao.
No Brasil a tarifa dos consumidores residenciais monmia, ou seja, existe apenas uma
tarifa, especificada em R$/kWh, que incluiu simultaneamente demanda e energia. J consumido-
res industriais so geralmente tarifados por meio de uma tarifa binmia, do tipo horo-sazonal.
Nesse tipo de tarifa a demanda cobrada em R$/kW, com valores diferentes para perodos de
ponta e fora de ponta. A energia cobrada em R$/MWh, com valores tambm diferentes para
perodos mido (dezembro a abril) e seco (maio a novembro). O horrio de ponta, no Brasil,
definido como o perodo de trs horas consecutivas, de escolha da distribuidora, compreendido
entres as 17h e as 22h.
Em pases mais desenvolvidos, nos quais existe algum tipo de Gerenciamento pelo Lado
da Demanda (GLD), existem tambm tarifas binmias para consumidores residenciais. Nesse
caso o consumidor paga mais caro, em R$/kWh no horrio de ponta, e mais barato, tambm em
R$/kWh, no horrio fora de ponta. A finalidade incentivar a migrao do consumo residencial
do horrio de ponta para o horrio fora de ponta, reduzindo a necessidade de investimentos em
distribuio para atendimento ao horrio de ponta. Uma maneira relativamente fcil de implantar
a GLD em um pas como o Brasil seria, por exemplo, pr-aquecer a gua durante o perodo fora
de ponta, armazenando-a em reservatrios trmicos especiais, para utilizao no horrio de pon-

ta, seja para o banho, seja para outro tipo de uso. Contudo, enquanto a energia tiver o mesmo
preo dentro e fora da ponta, esse tipo de GLD no teria sentido econmico para consumidores
residenciais. A ANEEL pretendia implantar em 2014 uma tarifa residencial denominada tarifa
branca, formada por trs componentes, todas em R$/kWh: uma componente reduzida, no hor-
rio fora de ponta, uma componente elevada, no horrio de ponta, e uma componente intermedia-
ria, uma hora antes do horrio de ponta e uma hora depois. Conduto, por causa de dificuldades
de implantao, a tarifa branca foi deixada para 2015.
Alguns conceitos importantes no estudo das cargas de SEPs so os seguintes:
1) Demanda mxima: valor mdio da carga durante o intervalo de tempo de meia hora
em que a demanda mxima.
2) Fator de carga: relao entre a demanda mdia e a mxima em um determinado in-
tervalo de tempo. O fator de carga ideal deve ser elevado. Caso seja unitrio, significa
que todas as unidades geradoras esto sendo utilizadas a plena carga durante o pero-
do considerado. Seu valor varia com a natureza da carga; sendo baixo para cargas de
iluminao (cerca de 12 %) e elevado para cargas industriais.
3) Fator de diversidade: relao entre a soma das demandas mximas individuais dos
consumidores e a demanda mxima do sistema. Este fator mede a diversificao da
carga e diz respeito capacidade de gerao e transmisso instalada. No caso da de-
manda mxima de todos os consumidores ocorrer simultaneamente, isto , fator de
diversidade unitrio, dever-se-o instalar muitos outros geradores. Felizmente, este
fator muito maior que a unidade, especialmente para consumidores residenciais. Em
um sistema de quatro consumidores o fator de diversidade poderia ser elevado, com
os consumidores absorvendo energia como na Figura 4.23.

Figura 4.23
Representao dos extremos do fator de diversidade
de uma instalao com dois consumidores
Em estudos de fluxo de potncia o ideal seria realizar um estudo para cada hora da curva
de carga da Figura 4.22. Isso, contudo, exigiria um esforo computacional muito grande, alm de
exigir uma previso de cargas muito complexa. Por outro lado, o modelo de dois patamares (pon-
ta e fora de ponta) adotado no nvel de distribuio (tenses inferiores a 230 kV), pouco descri-
tivo para estudos de sistemas de transmisso (tenses iguais ou superiores a 230 kV). Assim, em
estudos de transmisso geralmente adotamos o modelo de trs patamares (cargas mdia, leve e
pesada) da Rede Bsica brasileira, conforme mostrado na Tabela 4.2.
Tabela 4.2 Definio dos patamares de carga da Rede Bsica brasileira
Patamar Sem Horrio de Vero Com Horrio de Vero

de carga
2 feira Domingos 2 feira Domingos
sbado e feriados sbado e feriados
Leve 0h s 6h59 0h s 16h59 0h s 6h59 0h s 17h59
22h s 23h59 23h s 23h59
Mdia 7h s 17h59 17h s 21h59 7h s 18h59 18h s 22h59
21h s 23h59 22h s 23h59
Pesada 18h s 20h59 19h s 21h59
Para nossos fins, as cargas sero usualmente representadas em MVA ou MW, juntamente
com o fator de potncia, em um dos patamares da Tabela 4.2. No diagrama unifilar as cargas
sero representadas por meio de setas, como na Figura 4.29, indicando potncia absorvida, ou
por meio de impedncias, como na Figura 4.31.

4.9 Introduo aos estudos de curto-circuito.

Estudos de curto-circuito so necessrios no s em sistemas de potncia, mas tambm em sis-
temas industriais, e tm os seguintes objetivos gerais:
1) Ajustar rels de proteo e selecionar fusveis.

2) Selecionar os disjuntores que iro interromper as correntes de curto.
3) Estimar as consequncias das correntes de curto sobre cabos, transformadores, secciona-
doras, cabos para-raios, barramentos e outros equipamentos eltricos.
4) Determinar sobretenses em vrios pontos do sistema.
5) Permitir o dimensionamento de malhas de terra e de cabos para-raios.
6) Determinar as impedncias corretas dos transformadores de fora.
Os tipos de curto-circuito em um sistema trifsico so listados na Tabela 4.3 abaixo, jun-

tamente com as frequncias tpicas de ocorrncia.
Tabela 4.3 Tipos de faltas e estatsticas
Tipo de falta 69 kV 138 kV 230 kV

38,6% 36,7% 47,0%
Fase-terra
11,8% 10,0% 8,0%
Bifsico
(fase-fase)
25,5% 12,7% 5,0%
Bifsico-terra
(Fase-fase-terra)
6,3% 2,0% 0,6%
Trifsico
1,1% 0,7% 1,4%
Trifsico-terra
16,7% 37,9% 38,0%
Causa desconhecida
As causas dos curto-circuitos so diversas. Em linhas de transmisso as causas mais co-

muns so quedas de rvores, vendavais, descargas atmosfricas e vandalismo. No perodo seco,
quando as queimadas se tornam comuns, o ar pode se ionizar, provocando uma falta fase-fase
resultando em desligamento de sistemas. Em transformadores e geradores as faltas so menos
comuns e se devem a erros de operao e manuteno inadequada.
Em sistemas de potncia, compostos por geradores, transformadores, linhas e demais

equipamentos sempre equilibrados, os curtos trifsico e trifsico-terra resultam em corrente de
neutro nulo, sendo denominados faltas simtricas, por as correntes de curto so iguais em todas
as fases. O mesmo no acontece com os curtos fase-terra, fase-fase e fase-fase-terra, que produ-

zem correntes de curto diferentes em cada uma das fases, sendo denominados faltas assimtri-
cas.
A rigor, tanto faltas simtricas quanto assimtricas deveriam ser calculadas a partir das
tcnicas de fluxo de potncia, que sero vistas a partir do captulo 7, fazendo-se a impedncia de
curto igual a zero. Contudo, em sistemas de pequeno porte e em casos nos quais no se exige
muita preciso, podemos desenvolver uma metodologia simplificada, partindo das seguintes con-
sideraes:
1) A tenso pr-falta de todos os geradores igual a 1,0 pu. Sabendo que a tenso dos
geradores de um sistema de potencio pode variar entre 0,95 pu e 1,05 pu, a tenso
mais provvel de operao dos geradores 1,0 pu, onde a tenso-base a tenso no-
minal do gerador.
2) As cargas so desprezveis durante o curto, pois, sabendo que o sistema de pequeno
porte (poucas barras), a ocorrncia de um curto-circuito desvia das cargas toda a po-
tncia produzida pelos geradores.
3) As capacitncias em paralelo de linhas de transmisso tambm so desprezveis, pelo
mesmo motivo anterior.
A corrente trifsica (ou trifsica-terra) de curto-circuito franco, ou seja, sem impedncia

de curto, em uma determinada barra do sistema, pode agora ser determinada reduzindo-se o sis-
tema a um equivalente Thvenin cujas respectivas tenso e impedncia so Vth 1,00 e Zth .
A corrente de curto ser, portanto
Corrente trifsica de
V 1,00 curto-circuito franco em
I pu
cc 3 th , (4.34)

Z th Z th um sistema de potncia
de pequeno porte.
onde Zth a impedncia de Thvenin vista da barra onde ocorre o curto-circuito. Caso o curto se
d atravs de uma impedncia de falta Z f , basta adicion-la a Z , ou seja
th
Corrente trifsica de
Vth 1,00 curto-circuito atravs de
I ccpu3 , (4.35) uma impedncia em um
Zth Z f Zth Z f sistema de potncia de
pequeno porte.

O estudo das faltas assimtricas um pouco mais complexo, exigindo tcnicas especiais
que sero descritas no captulo 5, juntamente com vrios outros conceitos de curto-circuito.
Exemplo 4.6. Para o sistema da Figura 4.24, calcule a corrente trifsica de curto-circuito na bar-
ra 3, em pu e em amperes. Considere que a potncia-base 50 MVA e que a tenso-base na barra
3 69 kV.
Figura 4.24
Soluo. Inicialmente, substitumos os geradores por suas respectivas impedncias internas, des-
prezamos as cargas e isolamos a barra na qual desejamos calcular a falta. O resultado o dia-
grama de reatncias da Figura 4.25.
Figura 4.25
A impedncia equivalente de Thvenin, vista da barra 3, pode agora ser calculada
Zth j 0,10 j 0,20 j 0,15 //( j 0,42 // j 0,42) j 0,10
ou,
Zth j 0,45 // j 0,21 j 0,10 Zth j 0,1836 pu
Considerando as simplificaes feitas anteriormente, a corrente trifsica de curto-circuito

na barra 3 ser

1,0 pu 1,0 pu
Iccpu3 Iccpu3 j 5,448 pu

Z th j 0,1836
Para converter a corrente de curto para amperes, precisamos antes calcular a corrente-
base, que ser
Sb 50 10 6
I b3 I b3 418,37 A
3Vb3 3 69 103
Assim,
Icc3 Iccpu3 I b3 j 5,448 418,37 Icc3 j 2.279,29 A
____
Exemplo 4.7. Para o sistema da Figura 4.26, calcule a corrente trifsica de curto-circuito nas
barras 1 e 7. Utilize as bases de 60 MVA e 69 kV na barra 2.
Figura 4.26
Soluo. A Figura 4.27 ilustra o diagrama de reatncias resultante aps a converso para pu nas
bases indicadas, j com as cargas desprezadas e os geradores substitudos por suas respectivas
reatncias internas.
Quando o curto ocorre na barra 1, as barras 3 e 7 so flutuantes, pois as cargas nelas so

desprezadas. A impedncia equivalente de Thvenin ser ento a reatncia de j0,4 pu do gerador
1 em paralelo com a reatncia equivalente direita da barra 1, com as barras 3 e 7 abertas, ou
seja
Zth j 0,4 // j 0,16 j 0,45 j 0,105 j 0,1008 j 0,063 // j 0,0504 0,0756 ,

Zth j 0,4 // j 0,715 j 0,1638 // j 0,126 Zth j 0,2651 pu
A corrente de curto na barra 1, em pu, ser
1,0 pu 1,0 pu
Iccpu3 Iccpu3 j3,7719 pu
Z th j 0,2651
A respectiva corrente de curto em amperes ser
S 60 106
Icc3 Iccpu3 b j3,7719 j3,7719 502,04 Icc3 j1.893,66,3 A
3Vb3 69 10 3
3
Figura 4.27
O curto na barra 7 um pouco mais complicado, pois apenas a barra 3 ser flutuante e o
diagrama de reatncias resultante formar um delta entre as barras 2, 4 e 6, como mostrado na
Figura 4.28.
Figura 4.28
Sistema do Exempo 4.7 com curto na barra 7

A maneira mais fcil de calcular a reatncia equivalente na barra 7 transformar o delta

entre as barras 2, 4 e 6 para um estrela. Note que essa transformao nada tem a ver com a
transformao delta-estrela do transformador de trs enrolamentos. Usando as frmulas
tradicionais de transformao delta-estrela, teremos
x24 x26 j 0,1638 j 0,0504

x2 x2 j0,02849 pu
x24 x26 x46 j 0,1638 j 0,0504 j 0,0756
x24 x46 j 0,1638 j 0,0756

x4 x4 j0,04273 pu
x24 x26 x46 j 0,1638 j 0,0504 j 0,0756
x26 x46 j 0,0504 j 0,0756

x6 x6 j 0,01315 pu
x24 x26 x46 j 0,1638 j 0,0504 j 0,0756
A Figura 4.29 ilustra o circuito resultante com uma estrela entre as barras 2, 4 e 6.
Figura 4.29
Diagrama resultante para curto na barra 7
Agora fcil calcular a reatncia equivalente na barra 7
Zth j 0,195 j 0,01315 j 0,02849 j 0,56 // j 0,04273 0,555 ,
Zth j 0,20815 j 0,58849 // j 0,59773 Zth j 0,50479 pu
A corrente de curto na barra 7, em pu e em amperes, sero respectivamente
1,0 pu 1,0 pu
Iccpu3 Iccpu3 j1,981 pu

Z th j 0,50479
A respectiva corrente de curto em amperes ser

S 60 106
Icc3 Iccpu3 b j1,981 j1,981 2.309,4 Icc3 j 4.574,92 A
3Vb 7 15 10 3
3

4.10 Exerccios
4.10.1 Descreva algumas vantagens de se usar o sistema por unidade (pu) em vez das unidades
convencionais (volts, amperes, etc.).
4.10.2. Mostre como escrever em pu a impedncia de um equipamento, expressa originalmente

nas bases Vb1 e Sb1, quando integrado a um sistema cujas bases so Vb2 e Sb2.
4.10.3. Dois transformadores esto conectados em srie. Um deles especificado para 15 MVA,
69 kV/125 kV, X=10%. O outro, para 10 MVA, 13,8 kV/69 kV, X=8%. Determine a rea-
tncia de cada transformador e a reatncia total, em pu, nas bases de 30 MVA e 138 kV.
4.10.4 Trs transformadores monofsicos, 5 MVA, 8/2,2 kV, tm reatncia de disperso de 6% e

podem ser conectados de vrias formas de modo a suprir trs cargas resistivas idnticas
de 5 . Vrias conexes dos transformadores e cargas so ilustradas na Tabela 4.4 abai-
xo. Complete a tabela, usando potncia-base trifsica de 15 MVA (No se esquea de
mostrar os clculos!).
Tabela 4.4
Caso Conexo dos Carga Tenso-base Carga Impedncia

Transformadores conectada ao (kV, linha) R total vista do
secundrio lado de alta
Primrio Secundrio Alta Baixa (pu) (pu) ()
Y Y Y ? ? ? ? ?
1
2
Y Y ? ? ? ? ?
3
Y Y ? ? ? ? ?
4
Y ? ? ? ? ?
5 Y Y ? ? ? ? ?
6 Y ? ? ? ? ?
7 Y ? ? ? ? ?
8 ? ? ? ? ?
4.10.5 Trs geradores, cujos parmetros so listados na Tabela 4.5, esto conectados a um bar-
ramento comum de 13,8 kV. Determine a reatncia equivalente, resultante da ligao em
paralelo dos trs geradores, nas bases de 100 MVA e 15 kV.

Tabela 4.5
Gerador Potncia Tenso Reatncia

(MVA) (kV)
1 20 13,2 24%
2 50 13,8 20%
3 80 13,5 12%
4.10.6 Os geradores do exerccio 4.5 so conectados a um transformador de 160 MVA, 13,8/225

kV, 60 Hz, X=10%, o qual, por sua vez, conectado a uma linha de transmisso de 50
km de comprimento, cuja resistncia por fase 0,12 /km e cuja indutncia por fase
1,25 mH/km. Considerando que a linha de transmisso esteja a vazio, calcule a corrente
trifsica de curto-circuito em ambas as extremidades dela, em pu e em amperes.
4.10.7 Para o sistema da Figura 4.30, utilizando potncia-base de 50 MVA e tenso-base igual a
13,8 kV no barramento 1, pede-se: (a) converta para pu os valores de todos os parme-
tros; (b) calcule a tenso no barramento 1 de modo que a tenso no barramento 5 seja
0,95 pu.
Figura 4.30
Sistema para o Exerccio 4.11.7
4.10.8 Um gerador sncrono trifsico, 60 Hz, 50 MVA, 30 kV, tem reatncia sncrona igual a 9
por fase. A resistncia de armadura desprezvel. O gerador est entregando potncia
nominal com fator de potncia de 0,8 em atraso, sob tenso nominal, a um barramento in-
finito. Pede-se: (a) Determine a tenso interna do gerador e o ngulo de carga ; (b) com
a tenso interna mantida constante no valor do item anterior, a potncia de entrada do ge-
rador reduzida a 25 MW; determine a corrente e o fator de potncia.
4.10.9 Um transformador trifsico de dois enrolamento especificado para 60 kVA, 240/1200

V, 60 Hz. O rendimento deste transformador 96% e ocorre quando opera sob carga no-
minal e fator de potncia 0,8 em atraso. Determine as perdas no ferro e as perdas no co-
bre do transformador para rendimento mximo.

4.10.10 Para o sistema da Figura 4.31, pede-se: (a) converta todos os parmetros para pu,
usando potncia-base de 100 MVA e tenso-base de 15 kV na barra 9; (b) obtenha a
matriz admitncia nodal do sistema, em pu (Adaptado de ARLEI, 1998).
Figura 4.31
4.10.11 Considere o circuito equivalente de um transformador de dois enrolamentos, com

reatncias e resistncias dos lados primrio e secundrio, resistncia de perdas no fer-
ro e reatncia de magnetizao. Todas as grandezas esto referidas ao primrio. Pe-
de-se: (a) mostre que, para transformadores de potncia, o circuito equivalente pode
ser reduzido a uma nica reatncia; (b) mostre que, quando expressa em pu, a reatn-
cia do transformador de potncia tem o mesmo valor, independente de estar referida
ao primrio ou ao secundrio.
4.10.12 Uma linha de transmisso trifsica, 225 kV, tem comprimento de 40 km. A resistn-
cia por fase 0,15 /km e a indutncia por fase 1,326 mH/km. As capacitncias
em paralelo so desprezveis. Um transformador trifsico conectado a um dos lados
da linha, e uma carga de 380 MVA, com fator de potncia de 0,9 em atraso, sob 225
kV, conectada ao outro lado. Usando o modelo de linha curta, determine a tenso e
a potncia do lado do transformador.
4.10.13 Para o sistema da Figura 4.32, considere que a potncia-base 100 MVA e que a
tenso-base 15 kV na barra 4. Pede-se: (a) desenhe o diagrama unifilar para o sis-
tema, representando o transformador (que est fora do tap nominal) como um mode-
lo e incluindo as susceptncias das linhas; (b) converta todas as impedncias e ad-
mitncias para ohms e siemens, respectivamente; (c) desprezando cargas e suscep-
tncias em paralelo, calcule as impedncias Thvenin equivalentes nas barras um,

trs e quatro; (d) calcule as correntes trifsicas de curto-circuito nas barras um, trs e
quatro (Adaptado de ARLEI, 1998).
Figura 4.32
4.10.14 Repita o exerccio 4.13, considerando que a capacitncia em paralelo 0,0112

F/km, que o comprimento da linha 100 km e usando: (a) o modelo da linha m-
dia; (b) o modelo T da linha mdia.
4.10.15 Uma linha de transmisso trifsica, 345 kV, 130 km, tem impedncia srie por fase
igual a z=0,036+j0,3/km. A admitncia em paralelo por fase y=j4.2210-6 S/km.
Um dos lados da linha ligado a uma subestao e absorve 400 A, sob fator de po-
tncia 0,95 atrasado e 345 kV. Usando o modelo da linha mdia, determine a tenso,
corrente, potncia e fator de potncia do lado da carga.
4.10.16 Uma linha de transmisso de 500 kV tem comprimento de 250 km. A impedncia
srie por fase z0.045j0.4 /km e a admitncia em paralelo y=j410-6 S/km.
Determine os parmetros do modelo de linha longa para esta linha.
4.10.17 Uma linha de transmisso de 200 km conecta uma usina geradora a um sistema de
distribuio. Os parmetros da linha, por quilmetro e por fase, so: R=0,1; L=1,25
mH; C=0,01 F. Considerando que a tenso e a corrente do lado do sistema de dis-

tribuio so, respectivamente, 1320 kV e 164 37 amperes, calcule a tenso e

a corrente do lado da usina.
4.10.18 Um transformador trifsico de trs enrolamentos com tenses 132/33/6,6 kV tem as

seguintes reatncias em pu, medidas entre os enrolamentos de Alta, Mdia e Baixa
tenses e referidas a 30 MVA, 132 kV: xam = 0,15 , xab = 0,09 , xmb = 0,08 . O enro-
lamento de 6,6 kV alimenta uma carga equilibrada com corrente de 2.000 A e fator
de potncia 0,8 em atraso. O enrolamento de 33 kV alimenta um reator de j50,0
/fase conectado em estrela. Calcule a tenso no enrolamento de alta tenso para que
a tenso no enrolamento de baixa tenso seja de 6,6 kV.
4.10.19 Dois geradores com potncias individuais de 80 MVA esto conectados como ilus-
trado na Figura 4.33. A tenso nominal de cada gerador 11 kV e a reatncia sncro-
na de cada um deles 12%. Os transformadores elevadores so especificados para
11/66 kV, 90 MVA, e tm reatncia de disperso igual a 10%. Um curto-circuito tri-
fsico ocorre na barra 6 em um momento no qual nenhuma corrente circula atravs
da linha 2-5. Determine a corrente trifsica de curto-circuito na barra 6, em amperes.
Figura 4.33
4.10.20 Um gerador ligado por meio de um transformador a um motor sncrono. Quando

expressas na mesma base, as reatncias do gerador e do motor so 0,15 pu e 0,35 pu,
respectivamente, e a reatncia de disperso do transformador 0,10 pu. Uma falta
trifsica ocorre nos terminais do motor quando a tenso nos terminais do gerador
0,9 pu e a corrente de sada do gerador 1,0 pu com fator de potncia 0,8 adiantado.
Determine a corrente em pu da falta, no gerador e no motor. Use a tenso terminal do
gerador como fasor de referncia.

5. COMPONENTES SIMTRICAS
5.1. Introduo
O clculo de curto-circuitos assimtricos (fase-terra, fase-fase, fase-fase-terra) poderia, em prin-
cpio, ser realizado por meio das ferramentas convencionais de anlise de circuitos polifsicos
(malhas, ns, equivalentes, etc.). Contudo, o esforo computacional envolvido aumentaria com a
terceira potncia do nmero de barras do sistema, tornando a tarefa impossvel a partir de algu-
mas poucas barras. Felizmente, um teorema enunciado por Charles L. Fortescue em 1918 possi-
bilita a simplificao da anlise de faltas assimtricas, como veremos a seguir.
Nota biogrfica: Charles LeGeyt Fortescue (1876 1936) foi um engenheiro eltrico nascido
em York Factory, um entreposto comercial que funcionou at 1957 no noroeste da provncia ca-
nadense de Manitoba. Em 1898, Fortescue tornou-se um dos primeiros engenheiros graduados
pela Queens University, localizada em Ontario, Canad. Aps sua formatura, Fortescue ingres-
sou na Westinghouse Corporation, nos Estados Unidos, onde permaneceu durante toda sua vida
profissional, vindo a trabalhar com transformadores de alta tenso e problemas a eles relaciona-
dos. Em 1918, Fortescue publicou o artigo Method of symmetrical co-ordinates applied to the
solution of polyphase networks (AIEE Transactions, vol. 37, p. 1027-1140), dando origem ao
estudo das componentes simtricas e reduzindo enormemente o esforo computacional envolvi-
do nos clculos de curto-circuitos assimtricos. Em 1939, o IEEE (Institute of Electrical and
Electronics Engineers) criou uma bolsa de mestrado em homenagem a Fortescue, concedida
anualmente .
5.2. O teorema de Fortescue

O teorema de Fortescue aplica-se a sistemas polifsicos e foi originalmente enunciado da se-
guinte forma: qualquer sistema de N fasores desequilibrados, sendo N um nmero primo, pode
ser escrito como a soma de N conjuntos de fasores equilibrados.
O teorema de Fortescue pode ser escrito para tenses ou correntes. Para N=3 (sistema tri-
fsico), os trs conjuntos de fasores equilibrados so conhecidos como sequncias e definidos
da maneira a seguir.

1) Sequncia positiva
A sequncia positiva definida como a sequncia de fases do sistema em anlise, ou seja, a

sequncia de fases dos geradores conectados ao sistema. Usaremos o sobre-ndice 1 para re-
present-la, mas outros ndices usuais so + e abc.
A defasagem entre duas fases quaisquer da sequncia positiva sempre 120 e os mdu-
los das correntes (ou tenses) so iguais entre si. Denotando as trs fases por a, b e c, teremos
Ia1 Ib1 Ic1 . (5.1)

2) Sequncia negativa
A sequncia negativa gira no sentido inverso ao da sequncia positiva, tambm com ngulos de
120 entre duas fases quaisquer. Usaremos o sobre-ndice 2 para represent-la. Outros ndices
usuais so e cba. Da mesma forma que na sequncia positiva, teremos
Ia2 Ib2 Ic2 . (5.2)
3) Sequncia zero
Uma terceira sequncia, ou sistema de fasores, necessria para satisfazer o teorema de Fortes-
cue. Nesta sequncia, denominada zero e usualmente representada pelo sobre-ndice 0, os
fasores no giram, permanecendo paralelos entre si. Da mesma forma que nas sequncias anteri-
ores, teremos
Ia0 Ib0 Ic0 . (5.3)
A Figura 5.1 ilustra as trs sequncias em termos de seus fasores.
Figura 5.1
Sequncias de fase: (a) positiva;
(b) negativa; (c) zero

Intuitivamente, podemos imaginar as trs componentes de Fortescue como sendo produ-

zidas por geradores comuns. A sequncia positiva seria produzida por um campo girante trifsi-
co, girando no sentido da sequncia de fases do sistema. A sequncia negativa seria produzida
por um campo girante trifsico, mas girando no sentido oposto sequncia de fases do sistema.
J a sequncia zero seria produzida por um campo pulsante, no girante.
Dado um sistema de correntes desequilibradas Ia , Ib e Ic , o teorema de Fortescue pode

ser agora escrito em funo das componentes de sequncia positiva, negativa e zero
Ia Ia0 Ia1 Ia2 (5.4)

0 1 2
I b I b I b I b
0 1 2
I c I c I c I c
O sistema de equaes (5.4) pode ser simplificado introduzindo-se o operador unitrio a ,

definido como
a 1120 . (5.5) O operador unitrio a .
fcil verificar que o operador a tem as seguintes propriedades:
a 1120 1 240 (5.6)

a 1240 1 120
2
(5.7)
a 2 a * (5.8)
(a 2 )* a (5.9)
a 3 1360 a (5.10)
a 4 a 2 a 2 1 240 a (5.11)
1 a a 2 0 (5.12)
a a 2 a 3 0 (5.13)
Tomando a fase a como referncia, podemos agora escrever as correntes de sequncia

positiva da seguinte forma

Ia1 10 Ia1 (5.14)

1
I b 1240 Ia a Ia
1 2 1
1 1 1
I c 1120 I a aI a
De maneira semelhante, as correntes de sequncia negativa podem ser escritas como
Ia2 10 Ia2 (5.15)

2
I b 1120 Ia aIa
2 2
2 2 2 2
I c 1240 I a a I a
As correntes de sequncia zero podem ser escritas de maneira ainda mais simples
Ia0 Ib0 Ic0 . (5.16)
Substituindo as relaes (5.14), (5.15) e (5.16) na relao (5.4), teremos que
Ia I a0 Ia1 Ia2 (5.16a)

Ib Ia0 a 2 Ia1 aIa2 (5.16b)
I I0 aI1 a 2 I2
c a a a
(5.16c)
Escrevendo a relao acima em forma matricial, teremos
Ia 1 1 1 Ia0 (5.17)

I b 1 a
2
a Ia1
I 1 a a 2 I 2
c a
ou, em notao mais compacta

I A I ,
abc 012 (5.18)
onde
1 1 1 (5.19)

A 1 a 2 a
1 a a 2
O sobre-ndice abc denota o sistema desequilibrado original e o sobre-ndice 012 denota
o sistema de sequncia.

A matriz de transformao A tem algumas propriedades interessantes. Primeiro, pode-

mos notar que ela simtrica, ou seja,
A AT , (5.20)
onde o sobre-ndice T denota a matriz transposta. Alm disso, podemos verificar que
AT A* 3I , (5.21)
onde I a matriz-identidade. Este resultado ser til mais tarde. Finalmente, a matriz A
invertvel, com inversa dada por
1 1 1 (5.22)
A 1 a a 2
1 1
3
1 a 2 a
Pr-multiplicando a relao (5.18) por A , podemos agora obter as componentes de se-

1
quncia em funo das componentes do sistema abc original
I A I ,
012 1 abc (5.23)
ou
Ia0 1 1 1 Ia
Sistema de sequncia
1 1 012 escrito em termos do
I a 1 a a 2 Ib (5.24) sistema abc original.
I 2 3 1 a 2 a Ic
a
Note, da relao acima, que
I
Ia0 Ia Ib Ic n ,
1
3 3
onde In a corrente de neutro. Assim, s haver corrente de sequncia zero em circuitos nos
quais houver caminho para a corrente de neutro. Quando tal caminho no existir, como o
caso de conexes delta, a corrente de sequncia zero ser nula.

Exemplo 5.1. Usando a relao (5.24), calcule as correntes de sequncia para um sistema abc
equilibrado.
Soluo. Em um sistema totalmente equilibrado, teremos, por exemplo, Ia I a 0 ,
Ib I b 120 , Ic I c120 , onde I a I b I c . Aplicando (5.24), vem
Ia0 1 1 1 I a 0
1 1
I a 1 a a I a 120
2
I 2 3 1 a 2 a I a 120
a
ou,
0 I a
I a 3 1 1 120 1120 0

1 I a
I a 1 1120 1 120 1240 1120 I a
3
2 I a
I a 3 1 1240 1 120 1120 1120 0

Os resultados acima indicam que um sistema equilibrado com os ngulos 0, 120,
+120 (sistema de sequncia positiva) tem apenas componente de sequncia positiva. Se os n-
gulos fossem 0, +120, 120, caracterizando um sistema de sequncia negativa, apenas a com-
ponente de sequncia negativa existiria. Em ambos os casos a componente de sequncia zero
seria nula.
____
Exemplo 5.2. Em um sistema desequilibrado circulam as correntes Ia 80 , Ib 6 90 e
Ib 16143,1 . Calcule as correntes de sequncia e desenhe os diagramas fasoriais para cada

uma delas.
Soluo. De acordo com a relao (5.24), teremos
Ia0 1 1 1 80
1 1
I a 1 a a 6 90
2
I 2 3 1 a 2 a 16143,1
a
ou,

Ia0 1 1 1 80
1 1
I a 1 1120 1 120 6 90
I 2 3 1 1 120 1120 16143,1
a
ou, ainda,
Ia0 80 6 90 16143,1 2143,05

(5.25)
1 1 9,8118,38
I a 80 630 1623,1
I 2 3

80 6 210 16 263,1
4,3 86,08
a
De acordo com 5.14, teremos
Ib1 Ia1240 9,8118,38 240 9,81258,38 (5.26)

1 1
I c I a 120 9,8118,38 120 9,81138,38
E, da relao (5.15), teremos
Ib2 4,3 86,08 120 4,333,92 (5.27)

2
I c 4,3 86,08 240 4,3153,92
A Figura 5.2 ilustra o diagrama fasorial completo, mostrando a composio das correntes
de sequncia a partir das correntes do sistema abc original.
Figura 5.2
Diagrama fasorial mostrando a composio de um
sistema desequilibrado a partir de trs sistemas equilibrados

5.3. Potncia complexa

A potncia complexa trifsica do sistema original abc pode ser escrita, em pu, como

S3 V abc Iabc .
* (5.25)
Esse resultado pode ser escrito tambm em forma matricial

S3 V abc I abc ,
T * (5.26)
bastando que se defina os seguintes vetores-coluna
V a I a (5.27)

V abc V b , e I abc Ib
V c Ic

De (5.18) , sabemos que
I A I .
abc 012
Esse resultado vale tambm para tenses:
V A V .
abc 012
Assim, a relao (5.26) pode ser escrita como

S3 A V 012 A I ,
T 012 * (5.28)
ou,

S3 V 012 A
T
A I ,
T * 012 * (5.27)
Usando a relao (5.21), podemos finalmente escrever
Potncia complexa trif-

S3 3 V I ,
012 T 012 *
(5.29)
sica escrita em funo
dos componentes de
sequncia.
onde

V 0 I0 (5.30)

V 012 V 1 , e I012 I1 .
V 2 I 2

A relao (5.28) pode ser escrita em forma explcita
S3 3V0 I0* 3V1I1* 3V2 I2* . (5.31)
Decorre que a potncia total a soma das potncias de cada sequncia. Assim, cada um
dos trs circuitos de sequncia absorve uma parte da potncia total absorvida pelo circuito abc
original.
5.4. Impedncias de sequncia

Considere agora o circuito da Figura 5.3, que ilustra uma carga trifsica equilibrada, com impe-
dncia srie Z s por fase, ligada em estrela aterrada por uma impedncia de neutro Z n e alimenta-
da por uma fonte trifsica cujas tenses de fase so Va , Vb e Vc . As fases esto acopladas entre si
por meio de impedncias mtuas Z m , as quais podem ser resultantes de capacitncias ou indu-
tncias entre os condutores das linhas.
Figura 5.3
Carga trifsica equilibrada
com impedncias mtuas
As tenses de fase Va , Vb e Vc podem ser escritas como

Va Z s Ia Z m Ib Z m Ic Z n In

Vb Z s Ib Z m Ia Z m Ic Z n In

Vc Z s I c Z m I a Z m I b Z n I n
Lembrando que In Ia Ib Ic e reordenando os termos das equaes acima, vem
Va Z s Z n Ia Z m Z n Ib Z m Z n Ic

Vb Z m Z n Ia Z s Z n Ib Z m Z n Ic
Vc Z m Z n I a Z m Z n I b Z s Z n I c

ou, em forma matricial,
Va Z s Z n Z m Z n Z m Z n Ia (5.32)

Vb Z m Z n Z s Z n Z m Z n Ib
V Z Z Z m Z n Z s Z n Ic
c m n

A equao (5.32) pode ser escrita em forma mais compacta utilizando-se a notao matri-
cial
V Z I ,
abc abc abc (5.33)
onde
Z s Z n Z m Z n Z m Z n (5.34)
Z
abc
Z m Z n Z s Z n Z m Z n

Z m Z n Z Z Z s Z n
m n
a matriz-impedncia do sistema abc original.
Nossa inteno obter a equao de sequncia correspondente (5.33). Substituindo a re-

lao (5.18) na (5.33), teremos
A V Z A I .
012 abc 012 (5.35)
Pr-mutiplicando ambos os lados de (5.35) por A 1

, vem
V A Z A I .
012 1 abc 012 (5.36)


Sabendo que V 012 tem dimenso de volts e que I012 tem dimenso de amperes, ento,
por fora da lei de Ohm, o termo A Z A dever ter dimenso de ohms, sendo deno-
1 abc
minado matriz-impedncia de sequncia
Z A Z A .
012 1 abc (5.37)
Substituindo as relaes (5.19), (5.22) e (5.34) em (5.37), teremos, aps um calculo direto
Z s 3Z n 2Z m 0 0 Matriz-impedncia de
Z
012
0 Z s Z m 0 .

(5.38)
sequncia.
0 0 Z s Z m

A relao (5.37) deixa claro que as componentes simtricas funcionam como um mtodo
de diagonalizao da matriz-impedncia. A consequncia eltrica desse fato ainda mais inte-
ressante. Por exemplo, substituindo (5.38) em (5.36), podemos escrever
V 0 Z s 3Z n 2Z m 0 0 I
0 (5.39)
1
V 0 Z s Z m
0 I1
V 2 0 0 Z s Z m I 2

ou, de forma mais explcita
V 0 Z s 3Z n 2Z m I0 (5.40)
1
V Z s Z m I
1
V Z s Z m I
2 2
Assim, tenses de uma sequncia produziro correntes desta sequncia apenas. Em outras
palavras, os circuitos de sequncia so eletricamente desacoplados entre si.
5.5. Impedncias de sequncia dos componentes de um SEP

Um sistema equilibrado que opera alimentando cargas tambm equilibradas s contm compo-
nentes de sequncia positiva. Logo, as impedncias de sequncia positiva dos diversos compo-
nentes do circuito, tais como geradores, linhas de transmisso e transformadores, so as respecti-

vas impedncias j conhecidas. Contudo, precisamos analisar ainda a representao das impe-
dncias de sequncias negativa e zero de tais equipamentos.
5.4.1 Linhas de transmisso

As impedncias de sequncias positiva e negativa de uma linha de transmisso dependem apenas
da geometria desta e, logo, so idnticas, valendo a relao (5.41) abaixo.
J a reatncia de sequncia zero de uma linha de transmisso muito maior, por causa da
diferente distribuio de fluxos magnticos produzida pelas trs correntes em fase. Sendo Dn a
distncia entre a linha e o neutro, D a distncia entre as trs linhas e o comprimento da linha,
conforme mostrado na Figura 5.4, a reatncia de sequncia zero dada pela relao (5.42) a se-
guir.
Z LT
2
Z LT
1
, (), (5.41) Impedncias de sequn-
cias de uma linha de
transmisso. Note que a
D
0
xLT x1LT
1,2f ln n , (m/km). (5.42) reatncia em (5.42)
D dada em m/km.
3 3
10 10
Figura 5.4
Corte de uma linha de transmisso para
clculo da reatncia dada por (5.42)

5.4.2 Geradores sncronos

Como vimos na seo 4.7, o gerador sncrono caracterizado por trs diferentes reatncias: a
reatncia sncrona de eixo direto xd (correspondente ao funcionamento em regime), a reatncia
subtransitria de eixo direto xd ' ' (correspondente ao funcionamento no perodo subtransitrio) e
a reatncia transitria de eixo direto xd ' (correspondente ao funcionamento no perodo transit-

rio). Para uma reviso da definio de cada um desses perodos, consulte a Figura 4.18. A rea-
tncia de sequncia positiva do gerador ser ento igual a xd , xd ' ou xd ' ' , dependendo do
perodo no qual desejarmos calcular o curto-circuito.
Sabendo que a sequncia negativa gira no sentido contrrio da positiva, a reatncia de se-
quncia negativa do gerador dever ser calculada com o dobro da frequncia de operao. Uma
frmula prtica considerar que tal reatncia aproximadamente igual reatncia subtransitria
de eixo direto.
No caso da sequncia zero, as correntes giram junto com o campo girante. Logo, haver
apenas fluxo disperso, no fluxo magnetizante. A reatncia de sequncia zero ser portanto apro-
ximadamente igual reatncia de disperso da armadura.
x1g xd ou x1g xd ' ou x1g xd ' ' (), Impedncias de sequn-

(5.43) cias de um gerador sn-
dependendo do perodo desejado. crono.
xg2 xd ' ' , (), (5.44)
xg0 x , (). (5.45)
Alm de observarmos os valores de (5.43), (5.44) 4 (5.45), devemos tambm observar

que o tipo de conexo do gerador determinar o circuito a ser utilizado para geradores nos casos
das sequncias negativa e zero. Por exemplo, apenas a sequncia positiva gera tenso a vazio,
pois corresponde sequncia de fases do sistema. A sequncia negativa no gera tenso, de mo-
do que o respectivo circuito equivalente deve ter a f.e.m. Ef substituda por um curto-circuito.
Devemos nos lembrar tambm de que, no caso da sequncia zero, haver circulao de
corrente somente quando houver conexo ao terra. Assim, nos casos de conexo delta e estrela
aberta o circuito equivalente ser tambm aberto para sequncia zero. Finalmente, quando a co-

nexo for estrela, mas aterrada atravs de uma impedncia Zn, o circuito ser fechado para se-
quncia zero, mas a impedncia aparecer multiplicada por trs, conforme a equao (5.38).
Todos os circuitos equivalentes para geradores sncronos esto resumidos na Tabela 5.1.
5.4.3 Transformadores de dois enrolamentos

Como vimos na seo 4.4, transformadores de dois enrolamentos so representados de maneira
simplificada, considerando-se novamente apenas a impedncia de disperso, a qual ser igual s
impedncias para as sequncias positiva, negativa e zero.
Impedncias de sequn-
ZT1 ZT2 ZT0 Z , (). (5.46) cia para transformado-
res.
Tabela 5.1 Circuitos equivalentes para geradores sncronos
Conexo Sequncia Positiva Sequncia Negativa Sequncia Zero

Da mesma forma que no caso dos geradores, o tipo de conexo dos transformadores in-
fluenciar os circuitos para as sequncias negativa e zero, conforme mostrado na Tabela 5.2.
Tabela 5.2 Circuitos equivalentes para transformadores de dois enrolamentos

Tabela 5.2 Circuitos equivalentes para transformadores de dois enrolamentos (cont.)

No que diz respeito sequncia zero, as regras gerais, tanto para geradores quanto para
transformadores, so as seguintes:
1) A conexo estrela aterrada deixa passar corrente de sequncia zero, sem restries.
2) A conexo estrela aterrada por impedncia Zn deixa passar corrente de sequncia zero,
mas devemos adicionar a parcela 3Zn impedncia de sequncia zero do transforma-
dor.
3) A conexo estrela sem aterramento bloqueia completamente a passagem da corrente
de sequncia zero.
4) A conexo delta bloqueia a passagem da corrente de sequncia zero que sairia do
transformador; se o outro lado estiver ligado em estrela aterrada ou estrela aterrada
por impedncia, a corrente de sequncia zero ser desviada para o terra.
Os alunos que se deparam pela primeira vez com componentes simtricas geralmente en-
tendem como bastante naturais os circuitos para sequncia zero de transformadores conectados
em estrela, estrela aterrada e estrela aterrada por impedncia, mas veem como reservas a conexo
delta. Para melhorar a compreenso, devemos nos lembrar de que o transformador funciona por
compensao de fora magnetomotriz. Por exemplo, sendo IA a corrente no lado de alta e IB a
corrente no lado de baixa, devemos ter
Impedncias de sequn-
N A IA N B IB , (5.47) cia para transformado-
res.
onde N A e N B so os nmeros de espiras dos lados de alta e de baixa tenso, respectivamente.

Para que haja corrente de um lado, deve haver corrente do outro lado tambm. E, embora a co-
nexo delta no deixe passar corrente de sequncia zero, esta corrente circula dentro do delta.
Assim, em um transformador cujo lado de baixa(ou de alta) est ligado em delta, haver corrente
de sequncia zero no lado de alta (ou de baixa) se este estiver ligado em estrela aterrada solida-

mente aterrada ou aterrada por impedncia. As figuras 5.5 e 5.6 ilustram essa situao para a
conexo estrela aterrada-delta. Uma fonte de tenso monofsica foi ligada ao lado de alta, de
maneira a se simular a sequncia zero, o mesmo acontecendo com a impedncia de carga do lado
de baixa. A figura 5.6 deixa claro que h circulao de corrente no lado de baixa, por dentro do
delta. Logo, haver corrente tambm no lado de alta. Contudo, essa corrente no circula pela
carga e, assim, o circuito equivalente do lado de baixa aberto. No lado de alta o nico caminho
para a corrente para o terra, conforme ilustrado.
Figura 5.5
Transformador trifsico abaixador, conectado em estrela aterrada-
delta, ligado de maneira a simular a sequncia zero
A situao se torna um pouco mais complicada se tivermos uma fonte no lado em delta e
uma carga no lado em estrela aterrada. Contudo, podemos invocar a simetria implcita na relao
(5.47) e argumentar que o circuito equivalente da Figura 5.6 vale tambm para esse caso. Situa-
o semelhante ocorre no caso da conexo delta-delta (ltima linha da Tabela 5.2). Do ponto de
vista eltrico seria indiferente representarmos um circuito aberto de ambos os lados, com uma
impedncia duplamente aterrada no meio, ou representarmos apenas um circuito aberto, remo-
vendo a impedncia. Entretanto, novamente por razes de simetria, preferimos a primeira repre-
sentao.

Figura 5.6
Simplificao do circuito da Figura 5.5
5.4.4 Transformadores de trs enrolamentos

As regras para impedncias de sequncia de transformadores de trs enrolamentos decorrem das
regras j vistas para transformadores de dois enrolamentos. A Tabela 5.3 ilustra as impedncias
para algumas das ligaes mais comuns, na qual o modelo estrela utilizado.
Tabela 5.3 Circuitos equivalentes para transformadores de trs enrolamentos

Tabela 5.3 Circuitos equivalentes para transformadores de trs enrolamentos (cont.)

Exemplo 5.3. Um gerador sncrono trifsico, 25 MVA, 11 kV, tem reatncia subtransitria de
20%, reatncia de disperso de 1% e alimenta dois motores por meio de uma linha de transmis-
so e transformadores, conforme a Figura 5.7. Os motores so especificados para 15 MVA e 7,5
MVA, respectivamente, e ambos tm reatncia subtransitria de 25%, reatncia de disperso de
2% e tenso nominal de 10 kV. Os transformadores so ambos especificados para 30 MVA,
10,8/121 kV, com reatncia de disperso de 10% cada. A reatncia srie da linha 100 . Dese-
nhe os diagramas de sequncia positiva, negativa e zero para o perodo subtransitrio. Considere
que as bases do sistema so iguais aos dados nominais do gerador e que as impedncias de neu-
tro do gerador e do motor 2 so ambas iguais a 0,1 pu, j nas bases do gerador.
Figura 5.7
Soluo. Devemos antes escrever todas as reatncias nas bases do gerador: Sb=25 MVA, Vb1 =11
kV . As reatncias do gerador j esto na base correta, logo
xGd1 ' ' j 0,20 pu
xG1 j 0,01 pu
A reatncia do transformador 1-2, de acordo com a relao 4.4, ser

2
25 10,8
xT12 j 0,1 xT12 j 0,0803 pu
30 11,0
Para convertermos as demais reatncias, precisamos das tenses-base nas barra 3 e 4
121
Vb3 11 Vb3 123,24 kV
10,8
10,8
Vb3 123,24 Vb3 11 kV
121
As reatncias do transformador 3-4 e dos motores sero, respectivamente
2
25 121
xT34 j 0,1 xT34 j 0,0803 pu
30 123,24
2
25 10
xMd1 ' ' j 0,25 xMd1 ' ' j 0,3444 pu
15 11
2
25 10
xM1 j 0,02 xM 1 j 0,0275 pu
15 11
2
25 10
xMd2 ' ' j 0,25 xMd 2 ' ' j 0,689 pu
7,5 11
2
25 10
xM 2 j 0,02 xM 2 j 0,0551 pu
7,5 11
Finalmente, a reatncia da linha de transmisso 2-3 ser
j100
xLT23 xLT23 j 0,1646 pu
(123,24) 2 / 25
Da relao (5.43), sabemos que as reatncias de sequncia positiva dos geradores e moto-
res sero iguais s respectivas reatncias subtransitrias. Alm disso, as reatncias de sequncia
positiva dos transformadores so iguais s respectivas reatncias de transmisso e a reatncias de
sequncia positiva da linha de transmisso ser igual reatncia prpria da linha. Assim, o dia-
grama de sequncia positiva pode ser desenhado conforme a Figura 5.8 abaixo.

Figura 5.8
Circuito de sequncia positiva para o Exemplo 5.3
Sabendo que devemos desenhar o circuito de sequncia negativa para o perodo subtran-
sitrio, os valores das reatncias de sequncia negativa do gerador e dos motores so iguais s
respectivas reatncias de sequncia positiva. As reatncias dos transformadores e dos geradores
tambm permanecem as mesmas. Assim, o diagrama de sequncia negativa pode ser desenhado
conforme a Figura 5.9 abaixo.
Figura 5.9
Circuito de sequncia negativa para o Exemplo 5.3
O diagrama para sequncia zero mostrado na Figura 5.10 abaixo. Note a interrupo do
circuito nas barras 2 e 3, por causa das ligaes delta. Note tambm que as reatncias do gerador
e dos motores foram substitudas pelas respectivas reatncias de disperso, conforme recomen-
dado pela relao 5.45. Alm disso, as reatncias de neutro do gerador e do motor 2 aparecem
multiplicadas por trs, conforme a relao (5.39).

Figura 5.10
Circuito de sequncia zero para o Exemplo 5.3
____
Exemplo 5.4. Para o sistema da Figura 5.11, com os dados da Tabela 5.4, pede-se: (a) desenhe
os diagramas para as sequncias positiva, negativa e zero; (b) calcule as impedncias equivalen-
tes de Thvenin na barra 5 para as sequncias positiva, negativa e zero. As bases so Sb=120
MVA e Vb1=13,8 kV. Considere que a reatncia de neutro do gerador da barra 9 j0,4 pu, j
convertida para a base nova (Adaptado de ARLEI, 1998).
Figura 5.11
Sistema para o exemplo 5.4
Tabela 5.4 Impedncias originais do Exemplo 5.4
Equipamento x1 x2 x0
Gerador 1 11% 9% 1,6%
Transformador 1-2 8% 8% 8%
Linha 2-3 5+j11 5+j11 8+j20

Transformador 3-4-5, xab 10% 10% 10%

Transformador 3-4-5, xam 8% 8% 8%
Transformador 3-4-5, xmb 14% 14% 14%
Gerador 4 10% 8% 1,5%
Linha 5-6 3+j10 3+j10 7+j28

Linha 7-8 16+j40 16+j40 30+j90

Gerador 9 15% 12% 2%
O primeiro passo converter as reatncias e impedncias para a nova base. A Tabela 5.5
ilustra os resultados para as trs sequncias e as Figuras 5.12, 5.13 e 5.14 ilustram, respectiva-
mente, os circuitos para as sequncias positiva, negativa e zero. Note que as reatncias do trans-
formador 3-4-5 j foram convertidas para o modelo estrela, conforme as relaes (4.12), (4.13) e
(4.14).
Tabela 5.5 Impedncias convertidas para as bases novas
Equipamento x1 x2 x0
Gerador 1 11% 9% 1,6%
Transformador 1-2 7,38% 7,38% 7,38%
Linha 2-3 3,15% + j6,93% 3,15% + j6,93% 5,04% + j12,6%
Transformador 3-4-5, xa (5) 2,67% 2,67% 2,67%
Transformador 3-4-5, xb (4) 10,67% 10,67% 10,67%
Transformador 3-4-5, xm (3) 8% 8% 8%
Gerador 4 60% 48% 9%
Linha 5-6 0,68% + j2,27% 0,68% + j2,27% 1,59% + j6,35%
Transformador 6-7 28,8% 28,8% 28,8%
Linha 7-8 40,33% + j100,82% 40,33% + j100,82% 75,615 + j226,84%
Gerador 9 60% 48% 8%
Os circuitos para as sequncias positiva, negativa e zero so mostrados nas Figuras 5.12,
5.13 e 5.14, respectivamente.

Figura 5.12
Circuito de sequncia positiva para o exemplo 5.4
Figura 5.13
Circuito de sequncia negativa para o exemplo 5.4

Figura 5.14
Circuito de sequncia zero para o exemplo 5.4
A impedncia equivalente de Thvenin, vista da barra 5, para sequncia positiva, ser
Z th1 0,0068 j 0,0277 j 0,288 0,4033 j1,008 j 0,33 j 0,6 //

j 0,0267 j 0,08 0,0315 j 0,0693 j 0,0738 j 0,11 // j 0,1067 j 0,6
ou,
Z th1 0,4101 j 2,2537 // j 0,0267 0,0315 j 0,3331 // j 0,7067
Zth 0,01589 j0,22818 pu

1
Para sequncia negativa, teremos
Z th2 0,0068 j 0,0277 j 0,288 0,4033 j1,008 j 0,33 j 0,48 //

j 0,0267 j 0,08 0,0315 j 0,0693 j 0,0738 j 0,09 // j 0,1067 j 0,48
ou,
Z th2 0,4101 j 2,1337 // j 0,0267 0,0315 j 0,3131 // j 0,5867
Zth 0,01475 j 0,20902 pu

2
Finalmente, a impedncia de sequncia zero ser
Z th0 0,0156 j 0,0635 j 0,288 // j 0,0267
Zth 0,00008 j 0,02482 pu

0

Observe que, no caso da sequncia zero, toda a impedncia esquerda da reatncia de

j0,0267 do circuito original desconsiderada, por estar aterrada em ambas as extremidades.
O clculo de impedncias de sequncia ser importante no prximo captulo, quando

formos calcular correntes de curto-circuito assimtricos.
____
Exemplo 5.5. Desenhe os circuitos de sequncia negativa e zero para o sistema de potncia da
Figura 5.15 abaixo. Os valores-base so 50 MVA e 138 kV na barra 2.
Figura 5.15
Considere que as reatncias de sequncia negativa das mquinas sncronas so iguais s

respectivas reatncias subtransitrias e que as reatncias de sequncia zero das linhas de trans-
misso so iguais a 300% das respectivas reatncias de sequncia positiva.
Soluo. As tenses-base podem ser determinadas rapidamente por inspeo:
Vb2 Vb3 138 kV ,
Vb1 Vb4 13,2 kV ,
Vb5 Vb8 138 kV ,
Vb7 6,9 kV ,

Vb6 Vb9 138 kV ,
Vb4 13,2 kV .
As reatncias de sequncia negativa das linhas de transmisso podem ser escritas como
j 40
2
xLT xLT
2
j 0,105 pu
23
138 / 50
2 23
j 20
2
xLT xLT
2
xLT
2
xLT
2
j 0,0525 pu
58 69
138 / 50
2 58 69
Os transformadores esto todos na nova tenso-base, restando ajustar para a nova potn-
cia-base
50
xT212 xT215 xT234 xT246 j 0,1 xT2 xT2 xT2 xT2 j 0,2 pu
25 12 15 34 46
50
xT278 xT279 j 0,1 xT2 xT2 j 0,3333 pu
15 78 79
Finalmente, as reatncias de sequncia negativa dos geradores so

2
50 13,8
x x j 0,15 x g1 x g4 j 0,3279 pu
2 2 2 2
g1 g4
25 13,2
50
xg27 j 0,2 xg2 j 0,3333 pu
30 7
O diagrama de reatncias para sequncia negativa mostrado na Figura 5.16.

Figura 5.16
Diagrama de sequncia negativa para o Exemplo 5.5
As reatncias de sequncia zero das linhas e geradores so
0
xLT23
4 xLT
2
23
4 j 0,105 xLT
2
23
j 0,42 pu
0
xLT58
xLT
0
69
4 xLT
0
69
4 j 0,0525 xLT
0
58
xLT
0
69
j 0,21 pu
2
50 13,8
x g01 x g04 j 0,08 xg1 x g4 j 0,1749 pu
0 0
25 13,2
50
xg07 j 0,08 xg0 j 0,1333 pu
30 7
As reatncias de neutro de sequncia zero dos geradores so

2
50 13,8
x j 0,05 xng1 xng4 j 0,1093 pu
0 0 0 0
x ng1 ng4
25 13,2
50
0
xng j 0,05 xng
0
j 0,0833 pu
7
30 7
Finalmente, as reatncias de sequncia zero dos transformadores so iguais s respectivas

reatncias de sequncia negativa, calculadas anteriormente.

No caso do diagrama de sequncia zero devemos tomar cuidado de com as ligaes dos
transformadores, bem como adicionar as reatncias de neutro dos geradores, multiplicadas por
trs, em srie com as respectivas reatncias de sequncia zero. O resultado mostrado na Figura
5.17.
Figura 5.17
Diagrama de sequncia zero para o Exemplo 5.5

5.5 Exerccios
5.6.1. Enuncie o teorema de Fortescue e descreva suas vantagens no clculo de faltas assi-
mtricas (fase-terra, fase-fase, etc.). Esse teorema traria alguma vantagem no clculo
de faltas simtricas, tais como a trifsica e a trifsica-terra?
5.6.2. Seja um sistema eltrico cujas tenses em determinada barra so, em kV:
Va 13,87 , Vb 10,2100 , Vc 4,5 90 . Pede-se: (a) determine as tenses
de sequncia para as fases a, b e c; (b) desenhe o diagrama fasorial completo, ilus-
trando como as tenses do sistema original desequilibrado so formadas a partir das
somas adequadas das tenses de sequencia.
5.6.3. Seja um sistema eltrico cujas tenses de sequencias, para a fase a, so:
Va1 13,87 , Va2 10,2100 , Va0 4,5 90 . Pede-se: (a) determine as tenses
das fases a, b e c do sistema original desequilibrado; (b) desenhe o diagrama fasorial
completo, ilustrando como as tenses do sistema original desequilibrado so forma-
das a partir das somas adequadas das tenses de sequencia.
5.6.4. Esboce os diagramas de sequncia para transformadores de dois enrolamentos liga-

dos em: (a) delta-delta; (b) estrela aterrada-delta; (c) estrela aterrada estrela aterra-
da; (d) estrela-estrela aterrada.
5.6.5. Um transformador de dois enrolamentos est ligado em delta-estrela aterrada. Expli-

que porque, no caso da sequncia zero, haver corrente em ambos os enrolamentos,
mas no haver corrente para fora do lado em delta.
5.6.6. Considere um sistema composto por: (a) gerador trifsico, 13,8 kV, 50 MVA,
x=10%; (b) transformador de dois enrolamentos, 15 kV/69kV, 70 MVA, ligado em
estrela-estrela, com os dois neutros aterrados, x=8%; (c) linha de transmisso com
x=20 ohms; (d) carga de 20 MVA, com fator de potncia unitrio. Todos os elemen-
tos esto ligados em srie, na sequncia gerador, transformador, linha, carga. A po-
tncia base 100 MVA e a tenso base 15 kV na barra do gerador. Pede-se: (a)
converta os valores para pu; (b) calcule todas as impedncias equivalentes de Thve-
nin (sequncias positiva, negativa e zero) na barra de carga.

5.6.7. As tenses entre as fases de um sistema trifsico so Vab=218 V, Vbc=154,1 V e

Vca=154,1 V. A sequncia de fases positiva. Um conjunto de impedncias
Z an 270 , Z bn 3545 , Z cn 270 ligado em estrela s trs fases a,
b e c, na ordem indicada pelos subndices. Pede-se: (a) desenhe um esquema do cir-
cuito trifsico resultante; (b) determine as correntes de linha Ian , Ibn e Icn pelo m-
todo das componentes simtricas.
5.6.8. Considere o sistema da Figura 5.18, cujas impedncias esto representadas na Tabela
5.6. (a) Desenhe os diagramas de reatncias para as sequncias positiva, negativa e
zero, com todos os parmetros representados; (b) calcule as impedncias para as trs
sequncias nas trs barras (uma de cada vez).
Figura 5.18
Tabela 5.6 Dados do Exerccio 5.6.8
Equipamento x1 (pu) x2 (pu) x0 (pu)

G1 0,15 0,15 0,05
G2 0,15 0,15 0,05
T1 0,12 0,12 0,12
T2 0,12 0,12 0,12
L12 0,25 0,25 0,73
L13 0,15 0,15 0,4
L23 0,13 0,13 0,3

5.6.9. Um gerador trifsico de 30 MVA, 13,8 kV, possui uma reatncia subtransitria de
15%. Ele alimenta dois motores atravs de uma LT com dois trafos nas extremida-
des, conforme diagrama unifilar. Os valores nominais dos motores so 20 e 10 MVA,
ambos com 20% de reatncia subtransitria. Os trafos trifsicos so ambos de 35
MVA 13,2 - 115Y (kV), com reatncia de disperso de 10%. A reatncia em srie
da LT 80 . Faa o diagrama de reatncias com todos os valores em pu. Escolha os
valores nominais do gerador como base do circuito do prprio gerador.
Figura 5.19
5.6.10. Desenhe o diagrama de sequncia zero para o sistema da Figura 5.20 abaixo.
Figura 5.20
5.6.11. (a) Desenhe os circuitos de sequncia negativa e de sequncia zero para o sistema de
potncia da Figura 5.21. Expresse os valores de todas as reatncias em pu nas bases
30 MVA e 6.9 kV na barra 1. Os neutros dos geradores das barras 1 e 5 esto ligados
terra por meio de reatores limitadores de corrente com reatncia de 5%, cada qual

tendo como bases os valores dos respectivos geradores; (b) calcule as impedncias
equivalentes de Thevnin, para as sequencias negativa e zero, na barra 3.
Figura 5.21
5.6.12. Para o sistema do Exemplo 5.5, calcule as impedncias equivalentes de Thevnin,

para as sequencias negativa e zero, na barra 5.

6. CLCULO DE CURTO-CIRCUITO
6.1. Introduo
Neste captulo abordaremos o clculo de faltas assimtricas, quais sejam: fase-terra, fase-fase e
fase-fase-terra. Os conceitos introdutrios vistos na seo 4.10 continuam vlidos, mas agora
devemos aplicar o mtodo das componentes simtricas aos problemas em questo.
Inicialmente precisamos aprender a escrever as equaes para as trs sequncias de um

gerador a vazio. Sendo E a a tenso interna de fase, Z S a impedncia sncrona por fase e Va a
tenso de fase nos terminais da fase a de um gerador trifsico, os diagramas para as sequncias
positiva, negativa e zero podem ser representados como na Figura 6.1.
Figura 6.1
Circuitos de sequncia de um gerador a vazio:
(a) sequncia positiva; (b) sequncia negativa; (c) sequncia zero
As seguintes equaes podem ser abstradas dos circuitos da Figura 6.1 acima, os quais
so semelhantes aos circuitos da segunda linha da Tabela 5.1:
Va1 E a Z S1 Ia1 , (6.1)
Va2 Z S2 Ia2 , (6.2)
Va0 Z S0 Ia0 . (6.3)
As equaes (6.2) e (6.2) podem parecer um pouco estranhas, pois descrevem correntes
circulando e tenses terminais sem que haja fems internas. Contudo, devemos nos lembrar de que
tais equaes decorrem do teorema de Fortescue e, assim, apenas a soma das trs equaes acima
tem significado fsico. Note tambm que, como j vimos na Tabela 5.1, apenas o circuito de se-
quncia positiva apresenta fem interna ( E a ) no nula.

6.2. Curto-circuito fase-terra

Considere o circuito da Figura 6.2, que consiste de um gerador trifsico cuja fase a foi conectada
ao terra por meio de uma impedncia de falta Z f . O termo Z S representa a impedaimpedncia
sncrona do gerador em srie com qualquer impedncia a ele conectada, como impedncias de
linhas de transmisso, transformadores, etc. Por simplicidade, as fases b e c, que estariam ligadas
a cargas, so consideradas abertas, pois, conforme vimos no item 4.10, as cargas so considera-
das desprezveis durante um curto. Podemos ento escrever as seguintes condies de contorno
para o curto-circuito fase-terra
Va Z f Ia , (6.4)
Ib 0 , (6.5)
Ic 0 . (6.6)
Figura 6.2
Gerador com a fase a em curto com o terra
Escrevendo a transformao de Fortescue (5.24), com as correntes dadas por (6.5) e (6.6),
vem
Ia0 1 1 1 Ia (6.7)
1 1
I a 1 a a 2 0 ,
I 2 3 1 a 2 a 0
a
ou,

I (6.8)
Ia0 Ia1 Ia2 a .
3
Lembremos ainda que o teorema de Fortescue permite escrevermos a tenso na fase a

como
Va Va1 Va2 Va0 . (6.9)
Substituindo as relaes (6.1), (6.2) e (6.3) em (6.9), teremos
Va E a Z S1 Ia1 Z S2 Ia2 Z S0 Ia0 . (6.10)
A relao (6.8) nos garante que, na presente situao, as correntes de sequncia so

iguais. Logo, podemos escrever (6.10) como

Va E a Ia0 Z S1 Z S2 Z S0 . (6.11)
Substituindo (6.4) e (6.8) na relao acima, vem
I

(6.12)
Z f Ia E a a Z S1 Z S2 Z S0 .
3
Rearranjando os termos de (6.12) teremos
E a (6.13)
Ia0 ,
Z S1 Z S2 Z S0 3Z f
ou, finalmente, considerando que Ia 3Ia0
3E a Corrente de curto-
Iccft Ia (6.14) circuito fase-terra.
Z S1 Z S2 Z S0 3Z f
Conforme vimos na seo 4.10, E a a tenso pr-falta, considerada igual a 1,0 pu na

nossa formulao simplificada.
A corrente Ia0 , dada por (6.13), resulta de uma tenso E a aplicada a uma impedncia total
Z S1 Z S2 Z S0 3Z f . A Figura 6.3 mostra um circuito mnemnico1, que, por sua vez, ilustra a
1
Mnemnico vem de Mnemosine, a deusa grega da memria, e um elemento grfico ou verbal cuja finalidade
auxiliar a memorizao frmulas, listas ou outras informaes. Em outras palavras, trata-se de um macete.

relao (6.13) de maneira grfica. Note que este circuito nada mais do que a ligao em srie
dos circuitos da Figura 6.1, ligados ainda em srie a uma impedncia 3Z f .
Figura 6.3
Circuito mnemnico para o curto-circuito fase-terra
6.3. Curto-circuito fase-fase

Considere agora o circuito da Figura 6.4, no qual as fases b e c foram curto-circuitadas por meio
de uma impedncia de falta Z f e a fase a foi deixada em aberto. As condies de contorno so
agora
Vb Vc Z f Ib , (6.15)
Ib Ic Iccff , (6.16)
Ia 0 . (6.17)
Aplicando novamente a transformao de Fortescue (5.24) a essas correntes, teremos
Ia0 1 1 1 0 (6.18)
1 1
I a 1 a a 2 Ib ,
I 2 3 1 a 2 a Ib
a
ou,
Ia0 0 , (6.19)
1

Ia1 a a 2 Ib ,
3
(6.20)
1

Ia2 a 2 a Ib .
3
(6.21)

Figura 6.4
Gerador com as fase b e c em curto por meio de impedncia
De acordo com o teorema de Fortescue expresso por (5.16a), (5.16b) e (5.16c), as tenses
nas fases a, b e c podem ser escritas como
Va Va0 Va1 Va2 (6.22a)

Vb Va0 a 2Va1 aVa2 (6.22b)
V V 0 aV 1 a 2V 2
c a a a
(6.22c)
A tenso sobre a impedncia Z f ser ento

Vb Vc Z f Ib a 2 a Va1 Va2 , (6.23)
Substituindo (6.1) e (6.2) na relao acima, teremos

Z f Ib a 2 a E a Z S1 Ia1 Z S1 Ia2 . (6.24)
Das relaes (6.20) e (6.21) vem que Ia2 Ia1 . Logo, a relao (6.24) pode ser escrita
como

Z f Ib a 2 a E a Ia1 Z S1 Z S2 . (6.25)

De (6.2), temos ainda que Ib 3Ia1 / a a 2 , o que nos permite escrever (6.25) como

3Ia1 Z f

(6.26)
a a 2

a 2 a E a Ia1 Z S1 Z S2 .

Isolando Ia1 na relao acima, e levando em considerao que a a 2 a 2 a 3 , te-
remos
E a (6.27)
Ia1 .
Z Z S2 Z f
1
S
Substituindo (6.20) em (6.27) teremos finalmente a corrente de curto-circuito fase-fase
j 3E a Corrente de curto-
Iccff Ib (6.28) circuito fase-fase.
Z S1 Z S2 Z f
O circuito mnemnico para o curto-circuito fase-fase mostrado na Figura 6.5.
Figura 6.5
Circuito mnemnico para o curto fase-fase
6.4. Curto-circuito fase-fase-terra

Considere finalmente o circuito da Figura 6.6, no qual as fases b e c foram curto-circuitadas dire-
tamente e conectadas ao terra por meio de uma impedncia de falta Z f . A fase a foi deixada em
aberto. As condies de contorno so
Vb Vc Z f Ib Ic , (6.29)
Iccfft Ib Ic , (6.30)
Ia 0 . (6.31)

Figura 6.6
Gerador com as fase b e c em curto para o terra
Substituindo Vb Vc em (6.22b) e (6.22c), teremos
Va1 Va2 . (6.32)
Substituindo (5.16b) e (5.16c) em (6.29), teremos

Vb Z f Ib Ic Z f Ia0 a 2 Ia1 aIa2 Ia0 aIa1 a 2 Ia2 , (6.33)
ou,

Vb Z f 2 Ia0 a 2 a Ia1 a a 2 Ia2 , (6.34)
ou, ainda, considerando que a 2 a 1 ,

Vb Z f 2 Ia0 Ia1 Ia2 . (6.35)
Sabendo que Ia 0 , e considerando (5.16a), teremos
Ia0 Ia1 Ia2 0 . (6.36)
Substituindo (6.36) em (6.35), vem
Vb 3Z f Ia0 . (6.37)
Substituindo Va Va em (6.22b), teremos

1 2


Vb Va0 a 2 a Va1 , (6.38)
ou,
Vb Va0 Va1 . (6.39)
Igualando (6.39) e (6.37) e substituindo as relaes (6.1) e (6.3), vem que
3Z f Ia0 Z S0 Ia0 E a Z S1 Ia1 . (6.40)
Resolvendo para Ia0 , teremos finalmente
Z 1 I1 E a Corrente de sequncia
Ia0 S a (6.41) zero para o curto fase-
3Z f Z S0 fase-terra.
Comparando (6.29), (6.30) e (6.37), vem que
Iccfft 3Ia0 (6.42)

Corrente de curto fase-
fase-terra.
Agora falta apenas determinarmos Ia1 , valor que dever ser usado para a determinao de
Ia0 em (6.41). Considerando novamente que Va1 Va2 e igualando as relaes (6.1) e (6.2), tere-
mos inicialmente
1 1
I 2 Z S I a Ea , (6.43)
Z S2
a
Substituindo (6.41) e (6.43) em (6.36), vem
Z S1 Ia1 E a 1 Z S1 Ia1 E a
Ia 0. (6.44)
3Z f Z S0 Z S2
Agora basta isolar Ia1 , conforme o processo de clculo a seguir:
1 1 1 1
I 1 Ea Z S I a Ea Z S I a , (6.45)
3Z f Z S0 Z S2
a

2 1 1 0 1 1
I 1 Z S Ea Z S I a 3Z f Z S Ea Z S I a ,
a

Z S2 3Z f Z S0 (6.46)
2 1 2 1 1 1 0 0 1 1
I 1 Z S Ea Z S Z S I a 3Z f Ea 3Z f Z S I a Z S Ea Z S Z S I a ,
a
Z 2 3Z Z 0S f S (6.47)
E a
Ia1 1 2 0

Z S Z S Z S 3Z f 3Z f Z S2 Z S0 Z S2
,
(6.48)
Z 2 Z 0 3Z
S S f
E a
Ia1
Z
Z S 2S S 0
1
2 Z 0 3Z ,
f (6.49)

Z Z 3Z
S S f
Finalmente, notando que o segundo termo do denominador representa Z S em paralelo

2
com Z S 3Z f , podemos escrever:

0
E a Corrente de sequncia
Ia1

Z S1 Z S2 // 3Z f Z S0 (6.50) positiva para o curto
fase-fase-terra.
Para calcular a corrente de curto fase-fase-terra, devemos calcular inicialmente Ia1 , depois
Ia0 e depois Iccfft , de acordo com as relaes (6.50), (6.41) e (6.42), respectivamente.
O processo de clculo das correntes de curto pode parecer um pouco tedioso, mas tudo se
resume ao clculo das impedncias Z S1 , Z S2 , e Z S0 , identificadas com as impedncias equivalentes
de Thvenin das sequncias positiva, negativa e zero, respectivamente, e da aplicao de algu-
mas frmulas prontas. A Figura 6.7 ilustra o circuito mnemnico para o caso em questo.

Figura 6.7
Circuito mnemnico para o curto fase-fase-terra
____
Exemplo 6.1. Calcule as correntes de curto fase-terra, fase-fase e fase-fase-terra na barra 5 para
o sistema da Figura 6.7. Considere que as impedncias de neutro dos geradores so iguais a j0,1
pu e que a impedncia de falta igual a j0,2 pu. As demais impedncias so dadas na Tabela 6.1.
Figura 6.8
Tabela 6.1 Dados do Exemplo 6.1

G1 0,12 0,10 0,015
G3 0,12 0,10 0,015
T12 0,10 0,10 0,10
T34 0,10 0,10 0,10
L25 0,25 0,25 0,50
L24 0,15 0,15 0,30
L45 0,13 0,13 0,20

Soluo. Inicialmente devemos desenhar o circuito de sequncia positiva, ilustrado na Figura

6.9.
Figura 6.9
A impedncia de sequncia positiva nada mais do que a impedncia de Thvenin vista

da barra 5. Para determin-la, devemos substituir os geradores por suas respectivas impedncias
internas e transformar o delta entre as barras 2, 4 e 5 para uma estrela, conforme mostrado na
Figura 6.10.
Figura 6.10
Diagrama de reatncias de sequncia positiva para o Exemplo 6.1
As reatncias entre as barras 2, 4 e 5, convertidas para estrela, so:
x24 x25 j 0,15 j 0,25

x2 x2 j0,0708 pu
x24 x25 x45 j 0,15 j 0,25 j 0,13
x24 x45 j 0,15 j 0,13

x4 x4 j0,0368 pu
x24 x25 x45 j 0,15 j 0,25 j 0,13

x25 x45 j 0,25 j 0,13

x5 x5 j 0,0613 pu
x24 x25 x45 j 0,15 j 0,25 j 0,13
A impedncia de Thvenin de sequncia positiva na barra 5 ser
Zth1 5 j 0,0613 j 0,0708 j 0,10 j 0,12 // j 0,0368 j 0,10 j 0,12
Zth5 j 0,1977 pu
1
O circuito de sequncia negativa semelhante, exceto pelos valores das reatncias dos
geradores, conforme mostrado na Figura 6.11.
Figura 6.11
Diagrama de reatncias de sequncia negativa para o Exemplo 6.1
A impedncia de Thvenin de sequncia negativa na barra 5 ser
Zth2 5 j 0,0613 j 0,0708 j 0,10 j 0,10 // j 0,0368 j 0,10 j 0,10
Zth5 j 0,1876 pu
2
O diagrama de sequncia zero exige ateno triplicada. Em primeiro lugar, por causa dos
valores diferentes das reatncias dos geradores e das linhas de transmisso. Em segundo, por
causa das ligaes dos transformadores, especialmente o lado ligado em delta. Em terceiro, por
causa das reatncias de neutro dos geradores, que agora devem ser multiplicadas por trs e inclu-
das em srie com as respectivas reatncias de sequncia zero. O diagrama resultante mostrado
na Figura 6.12.
No que se refere ao clculo da impedncia de Thvenin de sequncia zero, agora no

necessrio converter o delta para estrela, pois a barra 2 apenas uma conexo entre as reatncias

de j0,50 e j0,30, que agora esto em srie. Logo, a impedncia de Thvenin de sequncia zero
ser
Zth0 5 j 0,30 j 0,015 j 0,10 j 0,20 // j 0,50 j 0,30
Zth5 j 0,575 pu
0
Figura 6.12
Diagrama de reatncias de sequncia zero para o Exemplo 6.1
As correntes de curto podem ser agora facilmente calculadas. De acordo com a relao
(4.35), a corrente de curto trifsico ser
1,00 1,00
I cc3

Icc3 j 2,5145 pu
Z th5 Z f
1
j 0,1977 j 0,2
A corrente de curto fase-terra ser, de acordo com (6.14)
3,00 3,00
Iccft
Z 1
th 5 Z th5 Z th5 3Z f
2 0
j 0,1977 j 0,1876 j 0,575 3 j 0,2
Iccft j1,9227 pu
De acordo com (6.28), a corrente de curto fase-fase ser
j 3 j 3
Iccff Iccff 2,9593 pu
Z 1
th 5
2
Z th 5 Z f j 0,1977 j 0,1876 j 0,2
De acordo com (6.50), a corrente de sequncia positiva para o curto fase-fase-terra ser
1,00 1,00
Ia1
Z 1
th 5
2
th 5
Z // 3Z f Z th5
0

j 0,1977 j 0,1876 // 3 j 0,2 j 0,575

Ia j 2,7819 pu
1
A corrente de sequncia zero para o curto fase-fase-terra ser, de acordo com (6.41)
I 0 Z th5 I a 1,0 j 0,1977 j 2,7819 1,0 Ia0 j 0,383 pu

1 1
3Z f Z th0 5 3 j 0,2 j 0,575
a
Finalmente, de acordo com (6.42), a corrente de curto fase-fase-terra ser
Iccfft 3 j 0,383 Iccfft j1,149 pu
6.5. Mtodo da matriz impedncia de barra

O mtodo das impedncias de Thvenin pode ser um pouco demorado caso desejemos calcular
as correntes de curto em mais de uma barra e, alm disso, de difcil generalizao e implemen-
tao computacional. Felizmente, o mtodo da matriz impedncia de barra possibilita tal genera-
lizao, bem como maior rapidez.
Seja um sistema de potncia de n barras, cuja matriz admitncia nodal pode ser escrita
como
Y11 Y12 Y1n

Y22 Y2 n

Y
Y 21

. (6.51)

Y1n Y2 n Ynn
Como sabemos das aulas de circuitos, a matriz admitncia nodal deve satisfazer a seguin-
te equao
I Y V ,

onde I e V so os vetores corrente e tenso respectivamente.
Relembrando das relaes (4.25), os elementos Yii da matriz admitncia nodal so iguais
soma de todas as admitncias que se ligam ao n i, enquanto as admitncias Yij Yji so iguais
ao recproco das respectivas admitncias fsicas que ligam os ns i e j. Logo, a matriz admitncia

nodal Y pode ser rapidamente construda a partir de regras simples. A matriz impedncia de
barra Z correspondente ser

Z11 Z12 Z1n

Z 22 Z 2 n

Z Y 21
1 Z

. (6.52)

Z1n Z 2 n Z nn
Aps obtermos a matriz impedncia de barra, a impedncia equivalente de Thve-

nin para curto-circuito na barra k igual ao elemento Z kk da matriz Z .
Em vrias circunstncias mais fcil inveter a matriz admitncia nodal, processo para o
qual existem algoritmos prontos e bastante conhecidos, do que calcular diretamente a impedncia
de Thvenin em cada barra do sistema.
____
Exemplo 6.2. Repita o Exemplo 6.1 pelo mtodo da matriz impedncia de barra.
Soluo. Inicialmente devemos redesenhar o circuito da figura 6.9, porm invertendo as impe-
dncias de modo a obtermos admitncias. O resultado mostrado na Figura 6.13.
Figura 6.13

De acordo com as regras (4.25), os elementos de Y 1 da diagonal principal sero
Y111 j8,333 j10,0 j18,333 ,
Y221 j10,0 j 4,0 j 6,667 j 20,667 ,
Y331 j8,333 j10,0 j18,333 ,
Y441 j10,0 j 7,692 j 6,667 j 24,359 ,

Y551 j 4,0 j 7,692 j11,692 .

De acordo com as mesmas regras (4.25), os elementos de Y fora da diagonal principal
sero
Y121 j10,0 j10,0 Y211 ,
Y131 0 Y311 (as barras 1 e 3 no se ligam diretamente) ,
Y141 0 Y411 ,
Y151 0 Y511 ,
Y231 0 Y321 ,
Y241 j 6,667 Y421 ,
Y251 j 4,0 Y521 ,
Y341 j10,0 Y431 ,
Y351 0 Y531 ,
Y451 j 7,692 Y541 .

A matriz admitncia nodal Y 1 ser
j18,333 j10,0 0 0 0
j10,0 j 20,667 0 j 6,667 j 4,0

Y 1 0 0 j18,333 j10,0 0 .

0 j 6,667 j10,0 j 24,359 j 7,692
0 j 4,0 0 j 7,692 j11,692

A inverso de Y 1 pode ser realizada por meio de um software numrico, como o Ma-
tLab ou o SciLab2, ou ainda por meio de uma calculadora cientfica, resultando em
2
O SciLab uma das alternativas freeware ao MatLab. Veja www.scilab.org .

j0,0937 j0,0718 j0,0263 j0,0482 j0,0563

j0,0718 j0,1316 j0,0482 j0,0884 j0,1032

1 1 1
Z Y j0,0263 j0,0482 j0,0937 j0,0718 j0,0637 .

j0,0482 j0,0884 j0,0718 j0,1316 j0,1168
j0,0563 j0,1032 j0,0637 j0,1168 j0,1977
Como pode ser observado, a matriz impedncia de barra simtrica, da mesma forma que
a matriz admitncia nodal. Contudo, enquanto a matriz admitncia nodal esparsa (tem muitos
elementos nulos), a matriz impedncia de barra cheia (tem poucos elementos nulos). Da matriz
Z tambm fica evidente que impedncia de curto na barra 5, para sequncia positiva, ser
1
Z55
1
j 0,1977 pu ,
que igual impedncia Z th1 5 obtida anteriormente no Exemplo 6.1.
O diagrama de admitncias para a sequncia negativa mostrado na Figura 6.14 abaixo.
Figura 6.14
Circuito de sequncia negativa para o Exemplo 6.2
Considerando que apenas as admitncias ligadas s barras 1 e 3 mudam, apenas os se-

guintes elementos da matriz admitncia nodal mudaro
Y112 j10 j10,0 j 20,0 ,
Y332 j10,0 j10,0 j 20,0 .

A matriz admitncia nodal Y 2 ser

j 20,0 j10,0 0 0 0
j10,0 j 20,667 0 j 6,667 j 4,0

Y 2 0 0 j 20,0 j10,0 0 .

0 j 6,667 j10,0 j 24,359 j 7,692
0 j 4,0 0 j 7,692 j11,692
A matriz impedncia de barra correspondente ser
j0,0803 j0,0606 j0,0197 j0,0394 j0,0467

j0,0606 j0,1212 j0,0394 j0,0788 j0,0933

1
Z 2 Y 2 j0,0197 j0,0394 j0,0803 j0,0606 j0,0533 .

j0,0394 j0,0788 j0,0606 j0,1212 j0,1067
j0,0467 j0,0933 j0,0533 j0,1067 j0,1876
A impedncia de curto na barra 5, para sequncia negativa, ser
Z552 j 0,1876 pu ,
que mais uma vez igual impedncia Z th2 5 obtida anteriormente no Exemplo 6.1.
O diagrama de admitncias para a sequncia zero mostrado na Figura 6.15 abaixo.
Figura 6.15
Circuito de sequncia zero para o Exemplo 6.2
Note, em relao ao circuito da Figura 6.15, que devemos primeiro adicionar a impedn-
cia de cada gerador e sua respectiva impedncia de neutro, e depois invert-las para obter a ad-
mitncia. Esse procedimento deve ser realizado sempre que tivermos algum equipamento, gera-
dor, motor ou transformador, aterrado por meio de impedncia. Feita tal considerao, os ele-
mentos da matriz admitncia nodal so

Y110 j3,1746 j10,0 j13,1746 ,
Y220 j 2,0 j3,333 j5,333 ,
Y330 j3,1746 j10,0 j13,1746 ,
Y440 j10,0 j3,333 j5,0 j18,333 ,
Y550 j 2,0 j5,0 j 7,0 .
Y120 0 Y210 (as barras 1 e 2 esto agora desconectadas),
Y130 0 Y310 ,
Y140 0 Y410 ,
Y150 0 Y510 ,
Y230 0 Y320 ,
Y240 j3,333 Y420 ,
Y250 j 2,0 Y520 ,
Y340 j10,0 Y430 ,
Y350 0 Y530 ,
Y450 j5,0 Y540 .

A matriz admitncia nodal Y 0 agora ser
j13,1746 0 0 0 0
0 j5,333 0 j3,333 j 2,0

Y 0 0 0 j13,1746 j10,0 0 .

0 j3,333 j10,0 j18,333 j5,0
0 j 2,0 0 j5,0 j 7,0
A matriz impedncia de barra correspondente ser

j0,0759 0 0 0 0
0 j0,6250 j0,3150 j0,4150 j0,4750

Z Y
0 0 1
0 j0,3150 j0,3150 j0,3150 j0,3150 .

0 j0,4150 j0,3150 j0,4150 j0,4150
0 j0,4750 j0,3150 j0,4150 j0,5750
Novamente, como esperado, teremos
Z550 j 0,5750 pu ,
que mais uma vez coincide com a impedncia Z th0 5 obtida anteriormente no Exemplo 6.1. As
correntes de curto-circuito na barra 5 podem ser agora facilmente calculadas, conforme vimos no
Exemplo 6.1. Alm disso, agora podemos calcular as correntes de curto em qualquer uma das
outras barras.
Devemos ter um cuidado especial em problemas que envolvam transformadores de trs

enrolamentos, como mostra o Exemplo 6.3 a seguir.
____
Exemplo 6.3. Calcule a correntes de curto-circuito trifsica e fase-terra na barra 3 do sistema da

Figura 6.16, pelos mtodos de Thvenin e da matriz impedncia. As reatncias dos geradores so
ambas iguais a j0,1 pu e as reatncias do transformador so: x12=j0,12 pu, x13=j0,13 pu e
x23=j0,14 pu. A impedncia de falta nula e, para os fins do presente problema, as reatncias de
sequncia positiva, negativa e zero de cada equipamento podem ser consideradas iguais.
Figura 6.16
Soluo. O primeiro passo desenhar o diagrama de reatncias de sequncia positiva, conforme

mostrado na Figura 6.17.

Figura 6.17
Diagrama de reatncias de sequncia positiva para o Exemplo 6.3
xxx
6.7. Exerccios
Observaes: (1) quando as impedncias de falta e de neutro no forem mencionadas, conside-
re-as nulas; (2) quando no mencionado, as impedncias esto especificadas nas bases dos res-
pectivos equipamentos; (3) quando no mencionado, os equipamentos esto ligados em estrela
solidamente aterrada.
6.7.1. Explique porque, em um curto-circuito fase-fase, no necessrio conhecer a impedncia

de sequncia zero.
6.7.2. (a) Quais as vantagens e desvantagens do clculo de curto-circuito por meio das impe-
dncias de Thvenin?; (b) quais as vantagens e desvantagens do clculo de curto-circuito
por meio da matriz impedncia de barra?
6.7.3. (a) Descreva as regras para a construo da matriz admitncia nodal de um sistema de
potncia; (b) como possvel obter a matriz impedncia de barra a partir da matriz admi-
tncia nodal?
6.7.4. Desenhe os circuitos mnemnicos para os curto-circuitos trifsico, fase-terra, fase-fase e

fase-fase-terra.
6.7.5. Um gerador encontra-se ligado a um transformador atravs de um disjuntor e tem os se-

guintes valores nominais: 7,5 MVA; 6,9 kV, com reatncias xd ' ' = 9%, xd ' = 15%, xd =
100%. O gerador opera em vazio e com tenso nominal quando ocorre um curto-circuito
trifsico entre o disjuntor e o transformador. Determine: (a) a corrente permanente de
curto-circuito no disjuntor; (b) a corrente inicial eficaz simtrica no disjuntor.

6.7.6. O transformador trifsico ligado ao gerador descrito no exerccio anterior tem os seguin-
tes valores nominais: 7,5 MVA; 6,9/15 kV, xd = 10%. Se ocorrer um curto-circuito trif-
sico no lado da alta tenso do transformador com tenso nominal e em vazio, determine:
(a) a corrente inicial eficaz simtrica no lado de alta tenso do transformador; (b) a cor-
rente inicial eficaz simtrica no lado da baixa tenso do transformador.
6.7.7. Um gerador de 60Hz, 650 kVA, 480 V, xd ' ' = 0,08 pu, alimenta uma carga puramente
resistiva de 500 kW, sob 480 V. A carga ligada diretamente aos terminais do gerador.
Se todas as trs fases da carga forem simultaneamente curto-circuitadas, determine a cor-
rente inicial eficaz simtrica no gerador, em pu, nas bases do gerador.
6.7.8. Um turbogerador, especificado para 60 Hz, 10 MVA e 13,8 kV, ligado em estrela soli-
damente aterrada e funciona sob tenso nominal, em vazio. As reatncias so xd ' ' = x2 =
0,15 pu e x0 = 0,05 pu. Determine a relao entre a corrente subtransitria para uma falta
fase-terra e a corrente subtransitria para uma falta trifsica simtrica.
6.7.9. No exerccio anterior, determine a relao entre a corrente subtransitria para uma falta
fase-fase e a corrente subtransitria para uma falta trifsica simtrica.
6.7.10. No Exerccio 5.6.11, considere que as impedncias de sequncia negativa e positiva so

iguais. Calcule as correntes de curto trifsico e fase-terra, na barra 3, pelo mtodo das
impedncias de Thvenin. Considere que a reatncia de falta nula.
6.7.11. Calcule todas as correntes de curto-circuito para o sistema do Exerccio 5.6.6, na barra
de carga, pelo mtodo das impedncias de Thvenin. Considere que a reatncia de falta
j0,2 pu.
6.7.12. Calcule todas as correntes de curto-circuito para o sistema do Exerccio 5.6.8, na barra 3,
pelos mtodos das impedncias de Thvenin e da matriz impedncia de barra. Considere
que a reatncia de falta nula.
6.7.13. Para o sistema da Figura 6.18, com os dados da Tabela 6.2, determine as correntes de
curto-circuito trifsica e fase-terra na barra 5.

Figura 6.18
Tabela 6.2 Dados do Exerccio 6.7.13

G1 0,14 0,12 0,05
G5 0,14 0,12 0,05
G8 0,14 0,12 0,05
T12 0,12 0,12 0,12
T78 0,12 0,12 0,12
T34 0,11 0,11 0,11
T46 0,11 0,11 0,11
T45 0,11 0,11 0,11
L23 0,12 0,12 0,3
L27 0,15 0,15 0,3
L45 0,12 0,12 0,25
L67 0,12 0,12 0,25
6.7.14. Calcule as correntes de curto-circuito fase-fase e fase-fase-terra para o sistema do Exerc-

cio 5.6.9, na barra 4, pelo mtodo das impedncias de Thvenin. Considere que a reatn-
cia de falta nula.

7. FLUXO DE POTNCIA
7.1. Introduo
Fluxo de potncia, tambm conhecido como fluxo de carga, um problema matemtico cujo
objetivo determinar as tenses e potncias em todos os barramentos de um sistema eltrico.
Como vimos no Exemplo 4.5, nos casos de real interesse o problema do fluxo de potncia no
linear e deve ser resolvido por meio de mtodos numricos, como Gauss-Seidel, Newton-
Raphson e outros. Iniciamos classificando os barramentos de um SEP, que podem ser divididos
em trs tipos bsicos:
1) Barramento de Carga (PQ): As potncias ativa e reativa (P e Q) so conhecidas,

mas o mdulo (V) e o ngulo da tenso () so desconhecidos.
2) Barramento de Gerao (PV): A potncia ativa (P) e o mdulo da tenso (V) so
conhecidos, mas a potncia reativa (Q) e o ngulo da tenso () so desconhecidos.
3) Barramento de referncia (V): O mdulo (V) e o ngulo da tenso () so conhe-
cidos, mas as potncias ativa e reativa (P e Q) so desconhecidas. Tambm denomi-
nado barramento flutuante, oscilante ou swing.
H apenas um barramento de referncia por sistema, mas pode haver vrios barramentos
de carga e de gerao. O fato de um barramento ter geradores no significa que seja um barra-
mento de gerao, ou seja, se um barramento tiver potncia ativa e mdulo da tenso conhecidos,
ele ser considerado PV, ou seja, de gerao. Ademais, barramentos que sirvam apenas de cone-
xo entre outros barramentos, sem geradores ou cargas conectados, cujas tenses e ngulos se-
jam desconhecidos, sero considerados barramentos PQ, com P=0 e S=0.
Em um problema de Fluxo de Potncia a configurao do sistema suposta inaltervel. O

objetivo determinar mdulos e ngulos de todas as tenses em todos os barramentos. Tendo as
tenses e as impedncias do sistema, podemos a seguir determinar o fluxo de potncia em cada
linha ou transformador, bem como as perdas no sistema.
Os problemas de fluxo de potncia no so lineares e geralmente apresentam uma com-

plexidade que s permite a soluo numrica. Em alguns casos didticos mais simples, contudo,
a soluo analtica possvel, como ilustrado no exemplo a seguir.

Exemplo 7.1. Dado o sistema da Figura 7.1, determine V3 , sendo V1 1,00 pu e
Sb 10 MVA .
Figura 7.1
Soluo. Inicialmente, a queda de tenso na impedncia entre as barras 1 e 3

V1 V3 Z13 I13 ,
1,0 V3 0,02 j 0,17 I13 . (7.1)
Por outro lado, sabemos que

S3 V3 I13
*
,
P3
arccos0,9 V3 I13* ,
cos
4,5 MW / 10 MVA
25,84 V3 I13* .
0,9
Isolando a corrente,
0,525,84 (7.2)
I13* .
V3
Substituindo (7.2) em (7.1), vem que
0,5 25,84
1,0 V3 0,02 j 0,17
V3*
0,08559 57,45 (7.3)
1,0 V3 .
V3*
Fazendo V3 V3 3 e reordenando, teremos
0,08559 57,45
1,0 V3 3 .
V3 3
Multiplicando por 1,0 3 ,
0,08559 57,45
1,0 3 V3 .
V3

Separando as partes real e imaginria,

cos 3 V cos57,45 V3
0,08559
3

sen 0,08559 sen57,45
3
V3
ou,
(7.4)
cos 3
0,046051
V3
V3
.
sen 0,072141
3
V3
Elevando ambas as equaes ao quadrado e as somando,

2 2
0,046051 0,072141
cos 3 sen 3
2 2
V3 ,
V3 V3
V3 2 0,046051 V3 2 0,0721412 ,
2
V3 2 0,00212 0,0921V3 2 V3 4 0,00520 ,

V3 4 0,90790 V3 2 0,00732 0 .
A equao acima uma biquadrada que pode ser resolvida diretamente pela frmula de
Bhaskara,
V3 2 0,90790 0,90790 4 0,00732

2
.
2
As razes so,
V3' 0,09023
''
V3 0,94856
A primeira raiz implicaria em uma queda de tenso excessiva entre as barras 1 e 2. As-
sim, escolhemos a segunda raiz.
V3 0,94856 pu
ou,
V3 13,09 kV
Podemos calcular o ngulo de V3 a partir da segunda equao do sistema (7.4)
sen 3
0,072141
3 4,362
0,94856

7.2. Mtodo de Gauss

Dentre os mtodos disponveis para a soluo de problemas de fluxo de potncia, o mtodo de
Gauss um dos mais simples e mais didticos, embora no o mais rpido nem o mais utilizado.
Vamos supor que inicialmente apenas uma equao no linear deva ser resolvida. A pri-
meira etapa do mtodo de Gauss consiste em se escrever a equao na forma x f x . A seguir,
a equao deve ser escrita na forma iterativa. Sendo k-1 o ndice da iterao inicial, teremos
Forma iterativa para

x f x
k
.
k 1
(7.5) utilizao do mtodo de
Gauss.
Finalmente, arbitramos um valor inicial para x e calculamos os valores das iteraes se-
guintes at que o erro x x

k k 1
entre o valor da ltima iterao e o valor da iterao anterior
seja to pequeno quando o desejado.
Em geral desejamos resolver um sistema de n equaes no lineares, que dever ser escrito
como
x1k
f1 x1k 1 , x2k 1 , x3k 1 ,..., xnk 1 (7.6)
k
x2 f x
2
k 1
1 , x2k 1 , x3k 1 ,..., xnk 1
f x
k k 1
x3 3 1 , x2k 1 , x3k 1 ,..., xnk 1

xnk

f n x1k 1 , x2k 1 , x3k 1 ,..., xnk 1
Nota biogrfica3: Carl Friedrich Gauss (1777 1855) foi um matem-

tico, fsico e astrnomo alemo que contribuiu para um grande nmero
de reas, como teoria dos nmeros, estatstica, anlise matemtica,
geometria diferencial, geodsia, geofsica, eletrosttica, astronomia e
tica. Conhecido como o Prncipe dos Matemticos, criana prod-
gio, muitos consideram Gauss o maior gnio da histria e alguns dizem
que seu QI teria girado em torno de 240. A lei de Gauss da distribuio
de erros e a curva normal em forma de sino so hoje conhecidas de todos que trabalham com
estatstica. A lei de Gauss da eletrosttica certamente conhecida de todos os estudantes de
3
DUNNINGTON, G. W. Carl Friedrich Gauss, titan of science. The Mathematical Association of America, 2003.

Engenharia Eltrica. Gauss tambm descobriu a possibilidade de se construir geometrias no-

euclidianas, embora nunca tenha publicado essa descoberta. Tais geometrias libertaram os
matemticos da crena de que os axiomas de Euclides eram todos consistentes e nao-
contraditrios e conduziram, anos mais tarde, teoria da relatividade geral de Albert Einstein,
dentre outras coisas. A partir de 1831 Gauss comeou a trabalhar em colaborao com o fsico
alemo Wilhelm Eduard Weber (1804 1891). Juntos, inventaram o primeiro telgrafo
eletromecnico4, que passou a interligar o observatrio astronmico e o instituto de fsica da
Universidade de Gttingen, na Alemanha. Gauss e Weber, motivados pela descoberta de Oersted
de 1921 (de que uma corrente eltrica produz um campo magntico), passaram a pesquisar se o
inverso no seria possvel, ou seja, se um campo magntico no seria capaz de produzir uma
corrente eltrica. Contudo, eles no foram bem sucedidos, pois no perceberam que, para que
isso acontea, o campo magntico deve ser varivel no tempo. A descoberta do fenmeno da
induo eletromagntica teve de esperar pelos trabalhos do fsico ingls Michael Faraday (1791
1867) e do fsico norte-americano Joseph Henry (1797 1878), em 1831. Gauss tem uma
cratera lunar batizada em sua homenagem.
____
Exemplo 7.2. Resolva o Exerccio 7.1 por meio do mtodo de Gauss. Considere um erro de 10-5.
Soluo. Escrevendo a equao (7.3) na forma iterativa, teremos
0,08559 57,45
V3k 1 1,0 .

V k
3
*
Fazendo V30 1,00 , o valor da primeira iterao ser
0,08559 57,45
V31 1,0 0,95667 4,325 .
1,00
A seguir, usamos V31 0,95667 4,325 para calcular a segunda iterao
0,08559 57,45
V32 1,0 0,94902 4,325 ,
0,95667 4,325
e assim por diante
0,08559 57,45
V33 1,0 0,94863 4,361 .
0,94902 4,325
4
O telgrafo de Gauss-Weber, ao contrrio do telgrafo eletromecnico de Samuel Morse, usava linguagem anal-
gica e no binria, consistindo de uma agulha que se movia distncia sob a influncia de uma fonte de tenso,
localizada distncia e cuja amplitude se variava.

0,08559 57,45
V34 1,0 0,94856 4,361 .
0,94863 4,361
0,08559 57,45
V35 1,0 0,94856 4,362 .
0,94956 4,361
Podemos observar que o erro entre a quarta e a quinta iteraes inferior a 10-5. Assim, o
problema converge com apenas cinco iteraes, mostrando que o mtodo de Gauss funciona de
fato. Contudo, precisamos desenvolver uma formulao que seja capaz de resolver problemas de
n barras e no de apenas uma barra.
Em um sistema de n barras, vale a seguinte equao matricial

I Y V . (7.7)
Para uma barra p qualquer, podemos escrever
I1 Y Y12 Y1 p Y1n V1
11
I 2 Y21 Y22 Y2 p Y2 n V2

.
Ip Yp1 Yp1 Ypp Ypn Vp

I Yn1 Yn 2 Ynp Ynn V
n n
Efetuando a multiplicao e escrevendo a equao para a barra p, obtemos

Ip Yp1V1 Yp 2V2 YppVp YpnVn . (7.8)
Sabemos ainda que

Ip S *p / Vp* . (7.9)
Igualando (7.8) e (7.9)

S *p / Vp* Yp1V1 Yp 2V2 YppVp YpnVn .
Isolando Vp
S * / V * Yp1V1 Yp 2V2 YpnVn

Vp p p .
Ypp
Generalizando e escrevendo na forma iterativa

Forma iterativa do
k 1
Vp
Ypp

1 k k * n k
S p / V p YpqVq .

(7.10)
mtodo de Gauss.
q 1
q p
____
Exemplo 7.3. Resolva o problema do Exemplo 7.1 usando a equao (7.10). Considere um erro
mnimo de 10-5.
Soluo. Devemos determinar V3 , sendo V1 1,00 pu e Sb 10 MVA . A barra 1 assim a
barra de referncia e a barra 3 a de carga.
A matriz admitncia nodal fcil de ser obtida:
j18,3333 j8,3333 0 (7.11)

Y j8,3333 6,8966 j 25,5747 6,8966 j17,2414
0 6,8966 j17,2414 6,8966 j17,2414
Os valores iniciais so:
V10 1,00
V20 1,00
V 0 1,00
3
S20 0
4,5 / 10
S30 arccos(0,9) 0,45 j 0,21794
0,9
A potncia da barra 3 considerada negativa por se tratar de uma barra de carga. A po-
tncia da barra 2 nula para todas as iteraes, por se tratar de uma barra genrica.
De acordo com a equao (7.10), a tenso na barra 2 para a primeira iterao ser:
V21
1

Y22

S20 / V20 Y21V10 Y23V30 .
*
0 j8,3333 1,0 6,8966 j17,2414 1,0

1
V21
6,8966 j 25,5747
6,8966 j 25,5747
V21 1,00
6,8966 j 25,5747

Da mesma forma, a tenso na barra 3 para a primeira iterao ser:
V31
1

Y33

S30 / V30 Y31V10 Y32V20 .
*
V31
0,45 j0,21794 / 1,0 0 1,0 (6,8966 j17,2414) 1,0
6,8966 j17,2414
6,4466 j17,02346
V31 0,98027 1,06
6,8966 j17,2414
Aps a primeira iterao, teremos:
V11 1,00
V 1 1,00
2
V31 0,98027 1,06

S 1 0
2
S31 0,8 j 0,3875

Prosseguindo, a tenso na barra 2 para a segunda iterao ser:
V22
1
Y

S21 / V21 Y21V11 Y23V31 .
22
*

0 j8,3333 1,0 6,8966 j17,2414 0,98027 1,06

1
V22
6,8966 j 25,5747
j8,3333 18,20319 110,74

1
V22
6,8966 j 25,5747
6,44626 j 25,3566
V22 0,98774 0,8272
6,8966 j 25,5747
A tenso na barra 3 para a segunda iterao ser:
V32
1

Y33

S31 / V31 Y31V11 Y32V21
*
(7.12)
V32
0,45 j0,21794 / 0,98027 1,06 0 1,0 (6,8966 j17,2414) 1,0
6,8966 j17,2414
V32
0,45487 j0,230779 6,8966 j17,2414
6,8966 j17,2414
6,44173 j17,01062
V32 0,97953 1,06
6,8966 j17,2414

O erro por enquanto 0,97953 0,98027 7,4 10 4 , ainda longe do erro mnimo
de 105 . O mtodo de Gauss tem convergncia lenta e oscilatria. Erros da ordem de 105 s
comearo a aparecer por volta da 28 iterao e o erro mnimo de 105 s se estabilizar a partir
da 45 iterao, resultando em
V245 0,96970 3,232
V 45 0,94854 4,362
3
A convergncia do mtodo de Gauss ilustrada na Figura 7.2.
Figura 7.2 Convergncia oscilatria do mtodo de Gauss
Pode parecer estranho que um problema que foi resolvido anteriormente com apenas cin-
co iteraes seja resolvido agora com quase dez vezes isso. Contudo, devemos lembrar que a
Equao (7.3) se aplica apenas a problemas de duas ou trs barras, enquanto a Equao 7.9 se
aplica a problemas de n barras (embora a convergncia possa se tornar lenta para muitas barras).
Problemas de fluxo de potncia geralmente so resolvidos por meio de programas de

computador especialmente desenvolvidos ou ento por meio de scrips para MatLab, SciLab ou
aplicativos semelhantes. Entretanto, pouco conhecido que problemas simples, de poucas barras,
podem ser resolvidos por meio do MicroSoft Excel, sem macros ou outro tipo de programao.

7.3. Mtodo de Gauss-Seidel

Como vimos na seo anterior, uma das desvantagens do mtodo de Gauss a convergncia len-
ta e oscilatria. O mtodo de Gauss-Seidel pode melhorar um pouco essa deficincia. Notemos,
por exemplo, na equao (7.12) do Exemplo 7.3, que usamos a tenso V21 para calcular V31 . Con-
tudo, nessa altura dos clculos j dispnhamos de V22 , que uma melhor estimativa de V2 do
que V21 . Podemos, assim, usar os valores das variveis assim que estiverem disponveis. O siste-
ma de equaes (7.6) pode ento ser escrito como
x1k
f1 x1k 1 , x2k 1 , x3k 1 ,..., xnk 1 (7.13)
k
x2 f x , x , x
2
k
1
k 1
2
k 1
3 ,..., xnk 1
f x , x , x
k k 1
x3 3
k
1
k
2 3 ,..., xnk 1

xnk

f n x1k , x2k , x3k , ..., xnk 1
Os mtodos de Gauss e de Gauss-Seidel apresentam, alm das j vistas, as seguintes propri-
edades:
1) Nmero de iteraes para a convergncia funo do nmero de barras.

2) Tempo computacional para a convergncia funo do quadrado do nmero de bar-
ras.
3) No requer inverso de matrizes.
4) Impossibilidade de se utilizar reatncias negativas.
5) Pouca sensibilidade aos valores iniciais.
6) Dificuldade para se encontrar erros de dados e problemas no sistema, pois mesmo
casos no convergentes apresentam resultados coerentes.
Nota biogrfica5: Philipp Ludwig von Seidel (1821 1896) nasceu em Zweibrcken, Alema-
nha, filho de um funcionrio dos correios. Em 1840, Seidel entrou para a Universidade de Ber-
lim, onde foi aluno do matemtico Johann Dirichlet e do astrnomo Johann Franz Encke, mu-
dando-se dois anos depois para a Universidade de Knigsberg, onde foi aluno dos matemticos
Friedrich Bessel, Carl Jacobi e Franz Neumann. Seus principais interesses vieram a ser tica e
anlise matemtica, assuntos de sua tese de doutoramento e de sua dissertao de habilitao,
respectivamente. Seidel desenvolveu, independentemente de outros matemticos como Karl
5
Fonte: O'CONNOR J.J.; ROBERTSON, E.F. Philipp Ludwig von Seidel. Disponvel: http://www-groups.dcs.st-
and.ac.uk/~history/Biographies/Seidel.html .

Weierstrass e George Gabriel Stokes, o conceito analtico de convergncia uniforme, relacionado

a sries de funes e que teve grande importncia na matemtica do final do sculo XIX. Seidel,
assim como Gauss, tem uma cratera lunar batizada em sua homenagem, por causa de seus traba-
lhos em tica e astronomia.
O exemplo a seguir ilustra uma soluo por Gauss-Seidel, assim como apresenta as conside-
raes especiais que devem ser feitas no caso de barramentos de gerao.
____
Exemplo 7.4. Dado o sistema da Figura 7.3, determine V1 e V2 , com um erro de 10-5, sendo
V3 1,020 pu, V1 1,01 pu e P1 30 MW . Considere Sb 100 MVA .
Figura 7.3
Soluo.
(continua...)
7.4. Mtodo de Newton-Raphson

Suponha que desejamos resolver o seguinte sistema de equaes
f1 x1 , x2 y1
, (7.11)
f 2 x1 , x2 y2
onde f1 e f2 so funes, x1 e x2 so variveis e y1 e y2 so constantes.
Sejam x10 e x20 os respectivos valores iniciais de x1 e x2, e x1 e x2 as correes neces-
srias para que as condies (7.11) sejam atendidas. Assim, podemos escrever


f1 x10 x1 , x20 x2 y1

, (7.12)
f 2 x10 x1 , x20 x2 y2
Expandindo as equaes (7.12) em srie de Taylor, teremos
f1 f1
f1 x1 x1 , x2 x2 f1 x1 , x2 x1 x x2 x ... y1
0 0 0 0
1 2
, (7.13)
f x x , x x f x , x x f f
0 0 0 0 2
x2 2
... y2
2 1 1 2 2 2 1 2 1
x1 x2
Desprezando as derivadas de ordem superior primeira e escrevendo o sistema de equa-
es em forma matricial, vem
f1 f1

y1 f1 x10 , x20 x1

0
x2
0 x1
,

0
y2 f 2 x1 , x2 f 2
0 f 2 x2
(7.14)
x1 x2
0 0
onde todas as derivadas devem ser calculadas a partir dos valores iniciais. Podemos escrever
(7.14) abreviadamente como
D J C, (7.15)
onde:
[D] = vetor das diferenas, tambm conhecido como vetor dos mismatches.
[J] = jacobiano de f1 e f2.
[C] = vetor das correes.
Invertendo (7.15) obtemos o vetor [C], que contm as correes desejadas

C J 1 D . (7.16)
Tendo x1 e x2, podemos escrever
x11 x10 x1
1 . (7.17)
x2 x20 x2
Da mesma forma que no caso do mtodo de Gauss-Seidel, o mtodo de Newton-Raphson
iterativo, ou seja, devemos continuar calculando as iteraes at que o vetor das correes se
torne menor do que um erro previamente especificado. Quando tal acontecer, dizemos que o sis-
tema convergiu.

Exemplo 7.5. Resolva o sistema abaixo pelo mtodo de Newton-Raphson, sendo x10 x20 2.
Considere um erro mnimo de 10-10.
2( x1 ) 2 x2 5
. (7.18)
x1 x2 6( x2 ) 2 57
Observao: note que, por exemplo, x12 o valor de x1 na segunda iterao, enquanto
(x1 ) 2 significa x1 elevado ao quadrado.

Soluo.
a) Primeira iterao
O vetor inicial das diferenas

y1 f1 x10 , x20 5 f1 2,2 5
D .

0

y2 f 2 x1 , x2 57 f 2 2,2 29
0 (7.19)
O jacobiano de f1 e f2
f1 f1
4 x 0 8 1
J xf 1 x2
0 0
1
01 . (7.20)
2 f 2 x2 x10 12 x20 x10 2 2 26
x1 0
x2
0
x2 2
0
O jacobiano pode ser invertido facilmente, resultando em

0,126214 0,00485
J 1 . (7.21)
0,00971 0,038835
As correes para a primeira iterao sero
x11 0,126214 0,00485 5 0.77184
1 29 1,174757 . (7.22)

2
x 0,00971 0,038835
Os erros so ainda muito superiores a 10-10. Logo, devemos calcular os novos valores de
x11 e x12 :
x11 x10 x11 2 0,77184 1,228155

1 . (7.23)
x2 x20 x12 2 1,174757 3,174757
b) Segunda iterao
O vetor inicial das diferenas

5 f1 1,228155 ; 3,174757 1,19149
D . (7.24)
57 f 2 1,228155 ; 3,174757 7,3736
O jacobiano de f1 e f2

4 1,3252427
J 5,3009708
1 1

1,3252427 12 2,3980583 2,3980583 30,10194
. (7.25)
2,3980583
Invertendo o jacobiano
0,206958 0,00516
J 1 . (7.26)
0,01671 0,025854
As correes para a primeira iterao sero
x12 0,206958 0,00516 1,19149 0,20778
2 7,3736 0,17073 . (7.27)

2
x 0,01671 0,025854
Os erros j diminuram um pouco, mas ainda so muito maiores do que 10 -10. Os valores
de x12 e x22 para a segunda iterao sero
x12 x11 x12 1,228155 0,20778 1,020372

2 . (7.28)
x2 x12 x22 3,174757 0,17073 3,004029
Prosseguindo com os clculos de maneira semelhante, perceberemos que na quinta itera-
o os erros j sero da ordem de 10-8. O sistema convergir na sexta iterao, com erro zero,
resultando nos seguintes valores:
x16 1,0
6 . (7.29)
x2 3,0
Desejamos agora aplicar o mtodo de Newton-Raphson em problemas de fluxo de pote-
potncia. Considerando um sistema de n barras, onde k e m so duas barras genricas, podemos
escrever
Vk Vk k = tenso na barra k, (7.30)
Vm Vm m = tenso na barra m, (7.31)

Ykm Ykm km Gkm jBkm = elemento da matriz admitncia Y . (7.32)
A potncia injetada na barra k pode ser escrita como
n
Sk Pk jQk Vk Im* Vk E m* Ykm
*
, (7.33)
m1
ou, escrevendo Vk , Vm e Ykm na forma polar

n
Sk Vk Vm* Ykm
*
( k m km ) . (7.34)
m1
Separando as partes real e imaginria, vem

n
V * Y * cos( ),
P
k V k m km k m km
m1
n
(7.35)
Q V
k Vm Ykm sen ( k m km ).
* *
k m1
As equaes (7.35) tambm podem ser escritas como

n
P
k Vk Vm Gkm cos km Bkm sen km ,
m1
n
(7.36)
Q V V G sen B cos ,
k k
m1
m km km km km
onde km k m (logo, sempre teremos kk 0 ).
As derivadas de interesse para a construo do jacobiano so

Pk Pk Qk Qk
, , , , (7.37)
m Vm m Vm
e o jacobiano ter o seguinte aspecto
P1 P1 P1 P1 P1 P1

2 n1 V1 V2 V
1
P2 P2

P2

P2 P2

P2
1 2 n1 V1 V2 V

P Pn1 Pn1 Pn1 Pn1 Pn1
n1
1 2 n1 V1 V2 V
J (7.38)
Q1 Q1 Q1 Q1 Q1 Q1

2 n1 V1 V2 V
1
Q2 Q2

Q2

Q2 Q2

Q2
1 2 n1 V1 V2 V

Q Q Q Q Q Q

1 2 n1 V1 V2 V
A numerao dos barramentos a seguinte:
1) Barramentos de Carga (PQ): 1 a .

2) Barramentos de Gerao (PV): l+ a n 1 .
3) Barramento de referncia (V): n .
O jacobiano pode tambm ser escrito em termos de submatrizes:

H N
J . , (7.39)
M L
onde
P P
H n1n1 , (7.40) Nn1 , (7.41)
V
Q Q
M n1 , (7.42) L . (7.43)
V
As submatrizes tambm podem ser escritas como

Pk
H km VkVm Gkm sen km Bkm cos km , (7.44)
m
Pk
H kk Vk2 Bkk Vk Vm Gkm sen km Bkm cos km , (7.45)
k mk
P
N km k Vk Gkm cos km Bkm sen km , (7.46)
Vm
Pk
N kk Vk Gkk Vm Gkm cos km Bkm sen km , (7.47)
Vk mk
Q
M km k VkVm Gkm cos km Bkm sen km , (7.48)
m
Qk
M kk Vk2Gkk Vk Vm Gkm cos km Bkm sen km , (7.49)
k mk
Q
Lkm k Vk Gkm sen km Bkm cos km , (7.50)
Vm
Qk
Lkk Vk Bkk Vm Gkm sen km Bkm cos km . (7.51)
Vk mk
Algumas caractersticas do jacobino so as seguintes:
1) Sendo n o nmero total de barras do sistema e o nmero de barras de carga, o jaco-

biano ter n 1 linhas e 2(n 1) colunas.
2) Os elementos fora da diagonal principal, correspondentes a barras no conectadas di-
retamente entre si, so nulos. Isso torna o jacobiano altamente esparso.
3) As submatrizes H, M, N e L so esparsas e simtricas.
O algoritmo para solucionar um problema de fluxo de potncia por Newton-Raphson o

seguinte:

1) Montar a matriz-admitncia nodal [ Y ].

2) Arbitrar valores iniciais para as variveis de estado de cada barra k: k0 , Vk0 .
3) Calcular Pk e Qk e verificar a convergncia inicial:
Pk Pk( dado) Pk(calculado) , (7.52)

Qk Qk( dado) Qk(calculado) . (7.53)
onde os sobrendices (e) e (c) significam especificado e calculado, respecti-
vamente.
4) Montar a matriz jacobiana [ J ].
5) Inverter [ J ].
6) Solucionar o sistema
( i 1)
P
(i )
J ,
1
(7.54)

V Q
onde (i) o nmero da iterao.

7) Atualizar as variveis de estado
( i 1) ( i 1)

(i )

. (7.55)

V V V
8) Calcular Pk e Qk e verificar a convergncia.

9) Se o sistema no tiver convergido, voltar ao passo (d).
Nota biogrfica: Joseph Raphson (1648 1715) foi um matemtico ingls que frequentou o
Jesus College, em Cambridge, e se formou em matemtica em 1692. Pouco antes, em 1689, Ra-
phson j havia sido eleito Fellow of the Royal Society, por indicao de Edmund Halley. As da-
tas so incertas, inclusive as de nascimento e morte, pois quase nada se sabe da vida de Raphson.
De qualquer forma, ele ficou conhecido pelo mtodo de Newton-Raphson e por ter traduzido
para o ingls a Arithmetica Universalis de Isaac Newton. Alm disso, Raphson foi um grande
defensor da hiptese de que Newton teria sido o nico inventor do Clculo, em oposio aos
defensores do alemo Gottfried Leibniz. Apesar de tal proximidade de interesses, Raphson e
Newton no eram amigos ou colaboradores, coisa que, na verdade, teria sido difcil a Newton. O
mtodo de Raphson para determinao de razes de equaes no lineares foi apresentado por ele
em seu Analysis Aequationum Universalis, de 1690. Newton desenvolveu um mtodo similar em
seu Mtodo das Fluxes, escrito em 1671, mas publicado somente em 1736. A verso de Ra-

phson, alm de ser superior de Newton, aquela que aparece nos livros atuais e, se no fosse a
importncia de Newton para a matemtica e fsica, seria conhecido apenas como Mtodo de
Raphson.
____
Exemplo 7.6. Monte o jacobiano para o sistema abaixo.
Figura 7.4
Soluo.
O sistema tem n=3 barras. Logo, o jacobiano ter dimenso 22. H apenas uma barra PQ (car-
ga). Logo, 1. Assim, conveniente renumerarmos os barramentos como em azul na figura,
de modo a podermos escrever:
a) Barramentos de Carga (PQ): 1 .
b) Barramentos de Gerao (PV): 1 2 .
c) Barramento de referncia (V): n 3 .
As submatrizes sero
H n1n1 H 22 , Nn1 N21 ,
M n1 M12 , L L11 .
O jacobiano ser
H11 H12 N12
H N H N 22
21 H 22
J . ,
M L M M 22 L22
21

ou, em termos de derivadas

P1 P1 P1
V2
1 2
P2 P2 P2
J 1 2
V2 .

Q Q2 Q2
2
1 2 V2
____
Exemplo 7.7. Resolva o sistema da Figura 7.5, que consiste de uma linha de transmisso repre-
sentada por um modelo , pelo mtodo de Newton-Raphson. Considere P 0,001 pu e
20 0 e demais dados da Tabela 7.1.
Tabela 7.1 Dados para o exemplo 7.7
Barra Tipo P (pu) Q (pu) V (pu)

1 V ? ? 1,0 0,0
2 PV 0,4 ? 1,0 ?
Figura 7.5 Sistema para o Exemplo 7.7. Valores em pu.
A barra 1 de referncia e, logo, conhecemos o mdulo e o ngulo da tenso nela. A bar-

ra 2 de gerao e conhecemos o mdulo da tenso nela, restando calcular o ngulo 2 , que a
nica varivel do problema.
Inicialmente devemos calcular a matriz admitncia:

1 1
Y12 Y21 j1,25 .
z12 j 0,8
1
Y22 Y11 B j 0,02 j1,25 j1,27 .
z12

j1,25
Y j1j1,27
(7.56)
G jB .
,25 j1,27
G
0,0 0,0 (7.57)
.
0,0 0,0
1,27 1,25
B
(7.58)
.
1,25 1,27
Sabendo que a barra 1 de referncia (V), o jacobiano ter apenas um elemento, corres-
pondente barra 2, de gerao (PV), ou seja
J H 22 P2 .
2
De acordo com (7.45), teremos
2
H 22 V B22 V2 Vm G2 m sen 2 m B2 m cos 2 m ,
2
2
m1
ou
H 22 V22 B22 V2 V1 G21sen 21 B21 cos 21 V2 G22sen 22 B22 cos 22 . (7.59)
ou, ainda, sabendo que 22 0 e V1 V2 1,0
H 22 1,25 cos21. (7.60)
A potncia na barra 2 pode ser escrita como

n
P2 V2 Vm G2 m cos 2 m B2 m sen 2 m ,
m1
ou
P2 V2 V1 G21 cos21 B21sen21 V2 G22 cos22 B22sen22 . (7.61)
ou, ainda
P2 1,25 sen21 (7.62)
Um teste inicial de convergncia pode ser feito para 21 2 1 0 0 0 e V2 1,0 .
De (7.52), temos
P2 P2( dado) P2( calc.) 0,4 0 0,4 . (7.63)
Logo, o sistema no converge, pois deveramos ter P2 0,001 pu . Devemos calcular a
primeira iterao, para a qual teremos

I
H 22 1,25 cos 21 1,25 cos(0) 1,25. (7.64)
O novo valor de 2I ser

P2 0,4
2I 0,32 rad.
H 22 1,25 (7.65)
O novo valor de 2 , por sua vez, ser
2I 20 2I 0 0,32 0,32 rad. (7.66)
Agora devemos calcular o novo valor da potncia (observe que os ngulos esto em radi-
anos)
P2I 1,25 sen 21 1,25 sen(0,32) 0,39321 rad. (7.67)
ou,
O erro ser
P2I P2( e) P2I 0,4 (0,39321) 0,00679 . (7.68)
Logo, o sistema ainda no convergiu, sendo necessria uma segunda iterao.

II
H 22 1,25 cos 21 1,25 cos(0,32) 1,18654 . (7.69)
O novo valor de 2II ser
P2I 0,00679
2II 0,00572 rad.
H 22 1,18654 (7.70)
O novo valor de 2 , por sua vez, ser
2II 2I 2II 0,32 0,00572 2II 0,32572 rad (7.71)
O novo valor da potncia ser

P2II 1,25 sen 21 1,25 sen(0,32572 ) 0,39999 pu. (7.72)
ou,
O erro ser
P2II P2( dado) P2II 0,4 (0,39999 ) 0,00001 . (7.73)
Agora o sistema convergiu dentro da margem de erro especificada. Resta somente calcu-
lar P1, Q1 e Q2.
A potncia na barra 1 pode ser escrita como

n
P1 V1 Vm G1m cos1m B1m sen1m ,
m1
ou
P1 V1 V1 G11 cos11 B11sen11 V2 G12 cos12 B12sen12 . (7.74)

ou, ainda, sabendo que V1 V2 1,0 , 11 0 e 12 1 2 0,32572 rad
P1 1,25 sen(0,32572 ) P1 0,39999 pu (7.75)
Notamos que, como no h cargas nem perdas hmicas na linha, as potncias nas barras 1
e 2 so iguais em mdulo. Da mesma forma, de (7.36), temos
2
Q1 V1 Vm G1m sen1m B1m cos1m , (7.76)
m1
ou,
Q1 V1V1 G11sen11 B11 cos11 V2 G12sen12 B12 cos12 , (7.77)
Q1 1,27 cos(0) 1,25 cos(0,32572 ) Q1 2,45428 pu (7.78)
Finalmente
2
Q2 V2 Vm G2 m sen 2 m B2 m cos 2 m , (7.79)
m1
Q2 V2 V1 G12sen12 B12 cos12 V2 G22sen22 B22 cos22 ,

Q2 1,25 cos(0,32572 ) 1,27 cos(0) Q2 2,45428 pu (7.80)
As caractersticas positivas do mtodo de Newton-Raphson aplicado a Sistemas de Po-

tncia so:
1) Converge quase sempre e mais rpido do que o mtodo de Gauss-Seidel.
2) O nmero de iteraes necessrias para a convergncia praticamente independente
do nmero de barras.
3) O tempo computacional para a convergncia diretamente proporcional ao nmero
de barras (e no ao quadrado do nmero de barras).
As caractersticas negativas do mtodo de Newton-Raphson so:

1) Necessidade de armazenar a matriz admitncia.
2) Necessidade de inverter o jacobiano a cada iterao.
____
Exemplo 7.8. Calcule V2 e V3 para trs iteraes pelo mtodo de Newton-Raphson para o siste-
ma a seguir.

Tabela 7.2 Dados para o exemplo 7.8
Barra Tipo V (pu) (rad) P (pu) Q (pu)

1 V(ref.) 1,0 0 ? 0
2 PQ (carga) ? 0 1,2 0
3 PQ (carga) ? 0 -1,5 0
Os elementos da matriz admitncia, em pu, so calculados abaixo

Y11 Y22 10 5 (0,1) 1 (0,2) 1 30 .
Y33 5 5 (0,1)1 (0,1)1 30 .
Y12 Y21 (0,2) 1 5 .
Y Y Y Y (0,1)1 10 .
13 31 23 32
A matriz admitncia pode ser escrita como

G11 G12 G13 30 5 10

Y G21 G22 G23 5 30 10 . (7.81)
G31 G32 G33 10 10 30
Por causa dos valores reais das impedncias, bem como da ausncia de reativos, apenas
P2, P3 e as derivadas em relao a V2 e V3 interessam. Considerando ainda que o sistema tem trs
barras ( n 3 ) e que duas barras so de carga ( 2 ), o jacobiano ter a seguinte forma:
P2 P2
N 23 V2 V3
J
N 22
.
N 33 P3
(7.82)
N 32 P3
V2 V3

De acordo com (7.46), e sabendo que os ngulos so todos nulos, podemos escrever
Pk
N km Vk Gkm .
Vm
Assim, teremos
N 23 V2G23 10V2 ,
N32 V3G32 10V3 .
Da mesma forma, de acordo com (7.47), teremos
3
N kk Vk Gkk VmGkm .
m 1
Assim,
N22 V2 G22 V1G21 V2G22 V3G23 V1G21 V3G23 5V1 10V3 .
N33 V3 G33 V1G31 V2G32 V3G33 V1G31 V2G32 10V1 10V2 .
O jacobiano ser, portanto
5V1 10V3
J
10V2
10V1 10V2
. (7.83)
10V3
Teste inicial de convergncia. Devemos calcular inicialmente P20 e P30 . De acordo com
(7.36), e considerando que Bkm 0 e km 0 , teremos

n
Pk Vk VmGkm .
m1
Alm disso, inicialmente teremos

V10 1,0 (fixo),
V20 1,0 (constante),
V30 1,0 (constante).
Assim,

P20 V20 V10G21 V20G22 V30G23 V20 5V10 30V20 10V30 ,
P20 1,0 5 1,0 30 1,0 10 1,0 15 ,
e

P30 V30 V10G31 V20G32 V30G33 V30 10V10 10V20 30V30 ,
P30 1,0 10 1,0 10 1,0 30 1,0 10 .
Para testar a convergncia, recorremos a (7.52)
P20 1,2 15 13,8 ,
P30 1,5 10 11,5 .

Observando que o sistema no convergiu, devemos passar primeira iterao, calculan-

do o jacobiano e seu inverso
J
15 10
, (7.84)
10 20
0,05
J 1
0,1
. (7.85)
0,05 0,075
Os novos valores das tenses sero
P20
I 0
V2 V2

V V J 1
0 .
3 3 P3
I
V2 1,0 0,1 0,05 13,8
V 1,0 0,05 0,075 11,5
3
I
V2 1,0 0,805 0,195
V 1,0 0,1725 0,8275
3
Calculando as novas potncias,

P2I V2I 5 30V2I 10V3I 0,195 5 30 0,195 10 0,8275 1,4479 ,
e

P3I V3I 10 10V2I 30V3I 0,8275 10 10 0,195 30 0,8275 10,6541 .
Testando a convergncia
P2I 1,2 (1,4479 ) 2,6479 ,
P3I 1,5 10,6541 12,1541 .
A convergncia ainda parece longe. Logo, devemos usar (7.83) para calcular o novo ja-
cobiano da segunda itereo
5V1 10V3
J
10V2
10V1 10V2
,
10V3
5 10 0,8275 10 0,195
J
10 0,8275 10 10 0,195
J
13,275 1,95

8,275 11,95
Invertendo o jacobiano
0,0136842
J
0,0838596
0,0931579
.
0,0580702
Os valores das tenses para a segunda iterao sero

P2I
II I
V2 V2
V V J I .
1
3 3 P3
II
V2 0,195 0,0838596 0,0136842 2,6479
V 0,8275 0,0580702 0,0931579 12,1541
3
II
V2 0,195 0,3884 0,5834
V 0,8275 1,286 0,4585
3
Calculando novamente as potncias

P2II V2II 5 30V2II 10V3II 0,5834 5 30 0,5834 10 0,4585 9,96864 ,
e

P3II V3II 10 10V2II 30V3II 0,4585 10 10 0,5834 30 0,4585 13,5666 .
Testando novamente a convergncia
P2I 1,2 9,9686 8,7686 ,
P3I 1,5 13,5666 15,0666 .
Ainda no convergiu. Logo, devemos repetir o procedimento.
7.5. Mtodo Desacoplado Rpido

Lembrando do mtodo de Newton-Raphson, as derivadas que interessam so do tipo:
P P Q Q
, , , .
V V
Em vrios casos observamos que as variveis PQ e QV so desacopladas, ou seja, inde-
pendentes, e podemos escrever
P P
, (7.86)
V
Q Q
. (7.87)
V
As matrizes M P / V e N Q / podem ento ser desprezadas e o sistema para a

iterao i pode ser escrito como
P
(i ) (i ) (i )
H 0
..... . ..... (7.88)
,
Q 0 L V
ou,

P(i ) H (i ) (i )
. (7.89)
Q(i ) L(i ) V (i )
A separao de (7.88) em duas equaes reflete o desacoplamento desejado. Em uma bar-

ra k qualquer, teremos
Pk Pk( dado) Pk(calc.) , para k de 1 a n-1, (7.90)
Qk Q ( dado)
k Q ( calc.)
k , para k de 1 a . (7.91)
Para acelerar a convergncia, podemos fazer
Pk Pk( dado) Pk( calc.)
, para k de 1 a n-1, (7.92)
Vk Vk
Qk Qk( dado) Qk( calc.)
, para k de 1 a . (7.93)
Vk Vk
O sistema (7.89) pode ser agora escrito como
P (i )
H '
(i ) (i )
V
, (7.94)

(i )

V L' V
Q (i ) (i )

onde
Vm Gkm sen km Bkm cos km ,

H km
H 'km (7.95)
Vk
Vk Bkk Vm Gkm sen km Bkm cos km ,

H kk
H 'kk (7.96)
Vk mk
L
L'km km Gkm sen km Bkm cos km , (7.97)
Vk
V G sen km Bkm cos km .

Lkk 1
L'kk Bkk m km (7.98)
Vk Vk mk
As desigualdade (7.86) e (7.87) so particularmente em sistemas de EAT (Extra-alta Ten-

so, entre 230 kV e 750 kV) e UAT (Ultra Alta Tenso, acima de 750 kV). Em tais sistemas as
seguintes hipteses simpleificadoras podem ser formuladas para o clculo H ' e L' :
a) Se o sistema pouco carregado, ento km pequeno. Logo cos km 1 .
b) A relao Bkm / Gkm elevada, entre 5 e 20. Logo, Bkm Gkm sen km .

c) As reatncias shunt so elevadas, muito maiores do que as reatncias srie. Logo,

BkkVk2 Qk .
d) As tenses Vk e Vm so sempre prximas a 1,0 pu.
Assim, podemos escrever

H 'km Bkm ,
H ' B ,
kk kk
(7.99)
L'km Bkm ,
L'kk Bkk .
Notamos que, com as aproximaes feitas, os elementos dos jacobianos no dependem

mais de tenses e ngulos variveis a cada iterao, dependendo agora somente de parmetros da
matriz admitncia, que so constantes. Em geral se faz agora as seguintes definies:
H ' B' , (7.100)
L' B ' ' . (7.101)
As matrizes B ' e B' ' so o negativo da parte imaginria (susceptncia) da matriz admi-
tncia. Em B ' aparecem as linhas e colunas referentes s barras PQ e PV, enquanto em B' ' apa-
recem apenas as linhas e colunas referentes s barras PQ. Podemos fazer ainda as seguintes
aproximaes adicionais:
a) Na matriz B ' so desprezados os elementos que afetam Q, isto , capacitncias shunt

e transformadores com comutao sob carga. As resistncias tambm so ignoradas
em B ' .
b) Na matriz B' ' so desprezados os elementos que afetam P, tais como transformadores
defasadores.
Sendo assim, sabendo que a reatncia entre as barras k e m xkm, e que o nmero de
barras diretamente ligadas barra k, podemos escrever os elementos de B ' e B' ' da seguinte
forma bastante simplificada:
1
B'km x ,
km

B'kk 1 ,
m1 xkm
(7.102)
m k
B' ' B ,
km km
B' 'kk Bkk .

O exemplo a seguir esclarece a utilizao do mtodo desacoplado rpido.
____
Exemplo 7.9. Formule as equaes para o mtodo desacoplado rpido para o sistema a seguir e
obtenha o ngulo 2 e a tenso V3 para as duas primeiras iteraes.
Tabela 7.3 Dados para o Exemplo 7.9
Barra Tipo V (pu) (rad) P (pu) Q (pu)

1 V(ref.) 1,0 0,0
2 PV (gerao) 1,0 ? 0,4 0,0
3 PQ (carga) ? ? 1,0 0,4
Tabela 7.4 Dados das linhas (em pu) do Exemplo 7.9
Linha r x Bshunt (total)

12 0,01 0,1 1,0
13 0,01 0,1 1,0
23 0,01 0,1 1,0
Continua em breve...por enquanto, s rabiscado em http://bit.ly/JcqmWf.

7.6. Exerccios
xxx

8. ESTABILIDADE ESTTICA E TRANSITRIA

8.1 Estabilidade em regime permanente
Os rotores de mquinas sncronas de p polos, funcionando a uma frequncia f1, giram a uma ve-
locidade sncrona dada, em rpm, por
120 f1
Ns . (8.1)
p
Em regime permanente a velocidade sncrona constante e esse caso que analisare-
mos inicialmente, por ser mais simples. Nesse caso, a potncia fornecida pelo gerador (ou absor-
vida pelo motor) ser tambm constante. Entretanto, podemos entender por regime permanente
tambm aquele no qual a potncia alterada muito lentamente, de modo que alteraes no ngu-
lo de carga , e tambm na frequncia f1, sejam desprezveis. Quando tais alteraes no forem
desprezveis estaremos falando em regime transitrio, que deixaremos para o item 8.2.
A relao (8.1) implica em detalhes construtivos da mquina. Por exemplo, mquinas de

alta velocidade, como o caso de turbogeradores, devero ter um nmero de polos reduzido
(geralmente dois ou quatro). Nesse caso razoavelmente fcil montar as bobinas do enrolamento
de campo, localizado no rotor, em ranhuras sobre uma estrutura cilndrica. A mquina ento
denominada de polos cilndricos ou de polos lisos6. Por outro lado, mquinas de baixa veloci-
dade, como o caso de hidrogeradores, devero ter um nmero elevado de polos (mais do que
dez, podendo chegar a mais de oitenta). Nesse caso mais conveniente enrolar as bobinas de
campo sobre sapatas polares, montadas sobre uma estrutura cilndrica. A presena das sapatas
resulta em salincias e a mquina ento denominada de polos salientes. O fato de a mquina
ser de polos lisos ou salientes afetar o comportamento da mesma, razo pela qual divideremos o
estudo a seguir em duas partes. A Tabela 8.1 resume as caractersticas das mquinas de polos
lisos e de polos salientes.
Tabela 8.1 Diferenas entre mquinas de polos lisos e de polos salientes

Tipo N de polos (p) Velocidade Tipo de gerador
Polos lisos Baixo (2 ou 4) Elevada Turbogerador
Polos salientes Elevado (mais de 10) Baixa Hidrogerador
6
Alguns autores, como JORDO (1980), preferem o termo polos no salientes, por entenderem que a presena das
ranhuras que alojam as bobinas de campo descaracteriza um rotor liso ou perfeitamente cilndrico.

8.2 Mquina de polos lisos em regime permanente

Considere uma mquina sncrona de polos lisos caracterizada por uma reatncia sncrona de eixo
direto xd e por uma resistncia de armadura r1, ambas medidas em ohms por fase. Suponha tam-
bm que, no referencial do estator, o fluxo magntico mtuo entre rotor e estator possa ser es-
crito como
2 (t ) N f 2 sen(et ) , (8.2)
onde N f o nmero de espiras equivalente por fase do estator, 2 o fluxo por polo do rotor,
e 2f1 a frequncia angular eltrica em radianos eltricos (rad-e) por segundo e um an-
gulo de fase. De acordo com a Lei de Faraday, a fora eletromotriz induzida em uma fase do
estator ser
d2 (t )
e f (t ) e N f 2 cos(et ) . (8.3)
dt
Podemos identificar e f (t ) com o seguinte fasor , cuja amplitude foi convertida para volts
eficazes e na qual foi includo ainda um fator de enrolamento kw, resultante das caractersticas
construtivas dos enrolamentos trifsicos
N k
E f e f w 2 . (8.4)
2
Fazendo ainda e N f kw / 2 m (uma constante) e considerando que a f.e.m. em (8.3) es-
t 90 atrasada em relao ao fluxo mtuo em (8.2), podemos escrever

E f jm 2 , (8.5)

onde 2 um fluxo pseudo-fasorial, pois o fluxo em (8.2) s varia no tempo quando no refe-
rencial do estator, mas no quando no referencial do rotor. Logo, trata-se de uma funo diferen-
te de e f (t ) , por exemplo, que varia no tempo quando vista de qualquer referencial.

A linha definida por 2 denominada eixo direto (d) e jaz sobre o eixo de simetria de

um polo norte. E f , estando a 90 de 2 , alinha-se ao longo do eixo em quadratura (q), que
jaz sobre o plano neutro do rotor, conforme ilustrado na Figura 8.1.

Figura 8.1
Definio dos eixos direto (d) e em quadratura (q)

Em situaes de regime podemos imaginar que o fluxo 2 gira no sentido anti-horrio
com velocidade constante N s , juntamente com todo o diagrama fasorial. Se houver um incre-
mente de potncia mecnica na entrada do gerador, o polo norte se adiantar levemente em rela-
o referncia, produzindo tambm um incremento em . Da mesma forma, redues na po-
tncia mecnica sero acompanhadas de redues de , apropriadamente denominado ngulo de
carga (ou de potncia). Os estudos de estabilidade, seja em regime permanente ou transitrio,
tm a finalidade de determinar se os valores iniciais e finais de correspondem a situaes est-
veis, i.e., situaes nas quais o gerador continua a operar em sincronismo com o barramento infi-
nito.
Suponha agora que a tenso interna por fase da mquina seja E f E f , que a tenso
terminal por fase, aqui tomada como referncia, seja V1 V10 e que a corrente de armadura
por fase seja I1 I1 . Podemos ento escrever a seguinte equao fasorial
Equao fasorial para

E f V1 r1 jxd I1 . (8.6) um gerador sncrono de
polos lisos.
O circuito unifilar equivalente para a relao (8.6) mostrado na Figura 8.2. No caso de
mquinas de grande potncia, como sempre ser o nosso caso, a resistncia r1 geralmente des-

prezvel em relao reatncia sncrona xd. A relao (8.6) est escrita para o caso de geradores.
Para o caso de motores, basta multiplicar a corrente por 1.
Figura 8.2
Circuito unifilar equivalente para uma mquina de polos lisos
Antes que o gerador possa produzir energia, ele precisa ser colocado em paralelo com o
barramento infinito. Para isso, necessrio que quatro condies sejam satisfeitas:
1) As sequncias de fases do gerador e do barramento devem ser as mesmas.

2) As frequncias do gerador e do barramento devem iguais.
3) As tenses do gerador e do barramento devem ser iguais.
4) Os ngulos de fase das tenses do gerador e do barramento devem ser iguais.
Sabendo que a tenso do gerador antes do sincronismo E f e que a tenso do barramen-
to V1 , as quatro condies acima podem ser escritas em uma nica equao fasorial.
Condio de paralelismo
E fq V1q . (8.7) entre um gerador e o
barramento infinito.
onde q 1,2,3 o ndice das fases. Podemos verificar que a relao (8.7) satisfaz imediatamente
s condies (a), (c) e (d). A condio (b) satisfeita de maneira implcita, pois a definio de
um fasor decorre da utilizao de um espao onde a frequncia nica e constante.
Supondo que I1 esteja atrasada de um ngulo em relao a V1 e que r1 seja desprez-

vel, podemos desenhar o diagrama fasorial ilustrado na Figura 8.3.
Duas relaes interessantes podem ser obtidas para as potncias ativa e reativa em funo
dos parmetros da mquina ( V1 , E f , xd , ), em vez de em funo dos parmetros da carga
( V1 , I1 , ), como usual. Considere inicialmente, na Figura 3, os tringulos OAB e ABC. Podemos

escrever, respectivamente, que

AB V1sen , (8.8)
AB xd I1sen . (8.9)
Figura 8.3
Diagrama fasorial para um gerador sobreexcitado
Igualando AB em (8.8) e (8.9) e considerando que 90 , teremos
V1sen xd I1sen(90 ) xd I1 cos( ) . (8.10)

Do segmento OBC e do tringulo ABC, podemos agora escrever
BC E f V1 cos , (8.11)
BC xd I1 cos . (8.12)
Igualando BC em (8.11) e (8.12) e considerando novamente que 90 , teremos
E f V1 cos xd I1 cos(90 ) xd I1sen( ) . (8.13)
Multiplicando (8.10) e (8.13) por V1 e desenvolvendo um pouco mais, vem
V12 sen xd V1I1 cos. cos V1I1sen.sen , (8.14)
V1 E f V12 cos xd V1 I1sen. cos V1 I1sen . cos . (8.15)
Substituindo P V1I1 cos e Q V1I1sen em (8.14) e (8.15) e escrevendo em forma

matricial, teremos
1 V1 sen cos sen P

2

xd V1 E f V12 cos sen cos Q
, (8.16)
Invertendo o sistema (8.16), podemos escrever
P 1 cos sen V12 sen

Q x sen ,
cos V1 E f V12 cos
(8.17)
d
Finalmente, desenvolvendo (8.17), teremos as relaes buscadas para P e Q

V1 E f Potncia ativa de uma

P sen . (8.18) mquina de polos lisos,
xd em watts/fase ou pu.
V1 E f V12 Potncia reativa de uma

Q cos . (8.19) mquina de polos lisos,
xd xd em vars/fase ou pu.
Em decorrncia das condies de sincronismo (8.7), um gerador recm-conectado ao bar-

ramento infinito ter V1 E f e Ang (V1 , E f ) 0 . Nessas circunstncias, as relaes (8.18) e
(8.19) indicam que as potncias ativa e reativa inicialmente fornecidas ao barramento sero nu-
las. Para transferir potncia ativa ao barramento devemos aumentar , o que feito por meio do
aumento da potncia mecnica no eixo da mquina. Mantendo-se E f constante, a potncia em
funo do ngulo de carga se comportar conforme mostra a Figura 8.4.
Figura 8.4
Potncia ativa em funo do ngulo de carga para um gerador de polos lisos
A potncia ativa mxima ocorre para ee / 2 rad-e, valor que denominaremos ngulo
de estabilidade esttica, pois a mquina no pode funcionar em regime permanente para ngu-
los acima deste valor e sair de sincronismo. Na regio estvel, supondo excitao constante,
aumentos da potncia mecnica Pm ec no eixo do gerador sero contrabalanados por aumentos da
potncia eltrica resistente oferecida pelo gerador. O ngulo permanecer oscilando levemen-
te em torno do ponto de equilbrio e a velocidade do rotor tambm oscilar levemente em torno
da velocidade sncrona. Acima de / 2 , aumentos da potncia mecnica produziro redu-

es da potncia resistente e a condio de velocidade estvel no poder ser mantida, a no ser

que se atue sobre o controle de velocidade da turbina ou que outra providncia seja tomada.
____
Exemplo 8.1. Seja um gerador trifsico de polos lisos, 10 MVA, 2,2 kV, ligado em estrela. A
reatncia sncrona de eixo direto 1,2 por fase e a resistncia de armadura pode ser despreza-
da. Considerando que a tenso interna seja igual nominal, pede-se: (a) determine a corrente de
armadura correspondente mxima potncia ativa; (b) desenhe o diagrama fasorial correspon-
dente.
Soluo. Para potncia mxima, devemos ter / 2 rad-e. Tomando ainda a tenso terminal
como referncia, teremos, de (8.6)
2.200 2.200 2.200 2.200

90 90
E f V1 3 3 3 3 1.058,5 j1.058,5 A
I1 .
jxd j1,2 j1,2
Convertendo para forma polar
I1 1.496,945 A .
A potncia ativa fornecida ser
V1l E1l 2.200 2.200
P sen sen90 4,03 MW .
xd 1,2
A potncia aparente ser
S 3V1l I1 3 2.200 1.496,9 5,7 MVA .

Note que a potncia acima, inferior ao valor nominal, a mxima potncia que o gerador
pode fornecer na atual condio de excitao. Aumentando-se a excitao o gerador poder che-
gar condio de potncia nominal, desde que se respeite a condio / 2 rad-e. O diagra-
ma fasorial mostrado na Figura 8.5 abaixo.

Figura 8.5
Diagrama fasorial para o Exemplo 8.1
8.3 Curva de capabilidade e curvas V

Alm de conhecermos o limite de estabilidade esttica necessrio conhecermos tambm todos
os limites operacionais da mquina, os quais so representados pela curva de capabilidade
(tambm conhecida como curva de capacidade). Para obter tal curva interessante desenhar-
mos antes o diagrama fasorial das potncias. Comeamos tomando a relao (8.6) e forando o
aparecimento da potncia aparente V1 / xd , multiplicando ambos os lados por V1 / xd e conside-
rando r1 desprezvel
V1 E f V1V1
jV1 I1 . (8.20)
xd xd
Supondo que o fator de potncia seja indutivo e tomando V1 como referncia, a relao
(8.20) poder ser escrita como
2
V1 E f V
1 V1 I1(90 ) . (8.21)
xd xd
O diagrama fasorial correspondente relao (8.21) mostrado na Figura 8.6. Estamos
agora interessados nos limites operacionais da mquina. O crculo ilustrado na Figura 8.6, por
exemplo, representa a potncia aparente mxima da mquina, seja qual for o fator de potncia.
Descartando-se outros componentes do conjunto turbina-gerador ou do conjunto motor-carga, a
potncia operacional da armadura poderia se localizar em qualquer ponto do crculo ilustrado.
Figura 8.6
Diagrama fasorial das potncias

Contudo, existem outros limites a considerar, conforme ilustrado na Figura 8.7. Conside-
re, por exemplo, que a mquina esteja funcionando como gerador. Se o limite da turbina fosse
especificado para coincidir com a potncia ativa mxima do gerador (Pmax), tal coincidncia se
daria apenas em um ponto, ou seja, o de fator de potncia unitrio. Assim, reduzindo um pouco o
limite da turbina que o gerador pode acionar (segmento OC), poderemos ter o conjunto turbina-
gerador funcionando em uma faixa mais extensa de fatores de potncia.
Figura 8.7
Limites operacionais de uma mquina sncrona
Situao semelhante ocorre com o limite de excitao mxima. Se estipulssemos tal li-
mite como coincidindo com Qmax, o campo do gerador estaria sobre-dimensionado, pois o limite
seria atingido somente para fator de potncia zero. Assim, dimensionamos o segmento OA como
o limite mximo de excitao, o que permite que a mquina funcione com excitao mxima em
uma faixa mais ampla de fatores de potncia.
Finalmente, temos o limite de estabilidade. De acordo com a Figura 8.4, a potncia m-

xima ocorre para =90. Contudo, uma margem de segurana deve ser estabelecida, pois o gera-
dor pode sair de sincronismo em =90 e esta situao deve ser evitada. Assim, o ngulo de es-
tabilidade pode ser definido, liustrativamente, em aproximadamente 40, de modo que, na Figura
8.7, o limite de estabilidade corresponde aos segmentos OB, para gerador, e OB, para motor. A
curva de capabilidade completa mostrada na Figura 8.8.

Figura 8.8
Curva de capabilidade de uma mquina sncrona
8.4 Mquina de polos salientes em regime permanente

Nas mquinas de polos lisos o entreferro pode ser considerado aproximadamente constante, o
que resulta em relutncia magntica constante ao longo do entreferro, bem como reatncia de
magnetizao e reatncia sncrona tambm constantes. Em mquinas de polos salientes, contudo,
tal constncia no se verifica e devemos adotar outro modelo para a mquina. E condies de
regime e de equilbrio de fases, os efeitos dos polos salientes podem ser considerados decom-
pondo-se a corrente de armadura I1 em duas componentes: uma componente Id ao longo do eixo
direto e uma componente Iq ao longo do eixo em quadratura, de modo que:
I1 Id Iq . (8.22)
Cada uma das componentes de I1 produzir quedas de tenso diferentes, associadas res-
pectivamente s reatncias sncronas de eixo direto ( xd ) e de eixo em quadratura ( xq ), conforme
mostrado na Fig. 8.9 para o caso de carga indutiva.
Do diagrama obtemos diretamente a equao fasorial do gerador sncrono
Equao fasorial para

E f V1 jxd Id jxq Iq , (8.23) um gerador sncrono de
polos salientes.

onde consideramos, por simplicidade, que r1 0 .
Figura 8.9
Diagrama fasorial para um gerador de polos
salientes alimentando carga indutiva
A potncia ativa igual a uma tenso multiplicada pela componente de corrente em fase
com ela. Por exemplo, sendo V1 a tenso, I1 a corrente e o ngulo entre V1 e I1 , temos a rela-
o convencional P V1I1 cos . J na Fig. 8.9 temos duas tenses ( Vd e Vq , as componentes de
V1 ao longo dos eixos d e q, respectivamente) e duas correntes ( Id e Iq ). Logo, podemos escre-
ver
P Vd I d Vq I q . (8.24)
Da Fig. 8.9, obtemos
Vd V1sen , (8.25)
Vq V1 cos . (8.26)
Conhecendo Vd e Vq , podemos tambm escrever
xq I q V1sen , (8.27)
xd I d E f V1 cos . (8.28)
Substituindo (8.25), (8.26), (8.27) e (8.28) em (8.24), teremos
E V1 cos V sen
P V1sen f V1 cos 1 , (8.29)
xd xq
que pode ser transformada em

V1E f 1 1
P sen V12 sen cos , (8.30)
xd x x
q d
e, finalmente, em
V1E f V12 1 1 Potncia ativa de uma

P sen sen2 . (8.31) mquina de polos salien-
xd 2 xq xd tes, em watts/fase ou pu.
Ao contrrio do que acontece na mquina de polos lisos, cuja potncia ativa varia senoi-
dalmente em funo do ngulo de carga, a potncia da mquina de polos salientes agora uma
soma de duas partes. A primeira parte, conforme mostrado na Fig. 8.10, denominada potncia
de excitao, por ser diretamente proporcional a E f . Tal potncia a mesma daquela da mqui-
na de polos lisos. A segunda parte, que no depende de E f , denominada potncia de relutn-
cia, por depender da diferena entre as relutncias de eixo direto e de eixo em quadratura7.
Figura 8.10
Potncia ativa de um gerador de polos salientes
Podemos usar uma considerao semelhante para deduzir uma expresso para a potncia
reativa. Neste caso, temos uma multiplicao entre as componentes de tenso e de corrente que
se encontram em quadratura, ou
Q Vd I q Vq I d . (8.32)
7
Lembre-se que a relutncia magntica inversamente proporcional reatncia indutiva.

Note o sinal negativo no segundo termo, dado o fato de Vd estar atrasado em relao a
Iq , enquanto Vq est adiantado em relao a Id . Substituindo os valores adequados em (8.32),
teremos
E V1 cos V sen
Q V1 cos f V1sen 1 , (8.33)
xd xq
ou,
V1E f cos2 sen 2

Q cos V1
2
. (8.34)
xd x x
d q
Finalmente, fazendo cos2 1 sen 2 ,
V1E f V2 1 1 Potncia reativa de uma

Q cos 1 V12 sen 2 . (8.35) mquina de polos salien-
xd xd x
d xq tes, em watts/fase ou pu.
Os dois primeiros termos correspondem potncia reativa de excitao e o terceiro termo

corresponde potncia reativa de relutncia.
Na mquina de polos lisos temos xq xd e podemos observar que as expresses (8.31) e
(8.35) redudem-se s expresses (8.18) e (8.19), respectivamente. Em outras palavras, a teoria

para polos salientes mais geral, englobando tambm a teoria para polos lisos.
O clculo das potncias P e Q exige a determinao do ngulo . Na mquina de polos li-
sos isso pode ser feito por meio do clculo da excitao E f , conforme a frmula (8.6). Na m-
quina de polos salientes, por outro lado, E f no pode ser calculada sem o clculo prvio de .
Podemos iniciar com a relao (8.27), aqui reescrita

V1sen xq I q . (8.36)
Para fator de potncia indutivo, a corrente I q pode ser obtida diretamente a partir da
Fig. 8.9
I q I1 cos . (8.37)
Substituindo (8.37) em (8.36), vem
V1sen xq I1 cos . (8.38)

ou,
V1sen xq I1 cos cos sen sen , (8.39)
ou, ainda
sen V1 xq I1sen xq I1 cos cos . (8.40)
Finalmente
xq I1 cos
tg . (8.41)
V1 xq I1sen
A relao (8.41) vale apenas para o caso indutivo. Para o caso capacitivo devemos notar
que a corrente I q deve ser escrita, a partir da Fig. 8.11, como
I q I1 cos . (8.42)
Repetindo o procedimento anterior, teremos
xq I1 cos
tg . (8.43)
V1 xq I1sen
A relao (8.43) foi obtida para carga capacitiva e , mas pode ser rapidamente veri-
ficado que a mesma tambm vlida para carga capacitiva e .
As frmulas (8.41), para carga indutiva, e (8.43), para carga indutiva, podem ser unifica-
das em uma s
xq I1 cos Clculo de para um

tg . (8.44) gerador de polos salien-
V1 xq I1sen tes.
Na frmula (8.44) devemos considerar o sinal adequado para o ngulo, ou seja, <0
para carga indutiva e >0 para carga capacitiva.
O algoritmo mais simples para clculo de E f e das correntes I d e I q , tendo-se a potncia
(ativa ou aparente), tenso, fator de potncia de um gerador e as reatncias xd e xq , o seguinte:
1. Calcule o fasor corrente de armadura, I1 .
2. A partir de I1 , calcule o ngulo , usando a relao (8.44).

3. Calcule Iq I1 cos( ) , onde o ngulo deve ser considerado com o sinal
adequado.
4. Calcule Id I1 Iq .
5. A partir de (8.23), calcule E f . Note que o ngulo de E f deve ser exatamente o
mesmo calculo no item 2.
Figura 8.11
Diagrama fasorial para um gerador de polos salientes alimentando carga capacitiva
____
Exemplo 8.2. Um gerador sncrono de polos salientes, 50 MVA, 2,2 kV, trabalha a plena carga
com fator de potncia 0,8 indutivo. A reatncia de eixo direto 1,0 pu e a reatncia de eixo em
quadratura 0,7 pu. Pede-se: (a) o ngulo de carga ; (b) a tenso interna Ef; (c) a regulao per-
centual de tenso; (d) o diagrama fasorial.
Soluo. Os clculos sero feitos em pu, com Sb=50 MVA e Vb=2,2 kV. A corrente de armadura,
a plena carga
*
I S1 1,036,87 1,0 36,87 pu .
*
1 V
1 1,0
A tangente do ngulo de carga ser
0,7 1,0 cos 36,87

tg 0,3944 .
1,0 0,7 1,0 sen 36,87

O respectivo ngulo de carga ser
atg0,3944 21,52
A corrente de eixo em quadratura ser
Iq 1,0 cos(21,52 (36,87))21,52 0,524121,52 pu .

E a respectiva corrente de eixo direto ser
Id 1,0 36,87 0,524121,52 0,8517 68,48 .

Finalmente, a tenso interna do gerador ser
E f 1,0 j1,0 0,8517 68,48 j 0,7 0,524121,52 ,

ou,
E f 1,0 0,8517 21,52 0,3669 111,52 E f 1,78221,52 pu .
Uma maneira mais simples de calcular E f , desde que no se deseje as correntes I d e I q ,
simplesmente isolar esta varivel em (8.31), resultando na relao a seguir:
V12 1 1 Tenso interna de uma

Ef
xd
P sen2 . (8.45) mquina de polos salien-
V1sen 2 xq xd tes, em volts ou pu.
Na relao acima no podemos nos esquecer que a potncia P dada em watts/fase ou

em pu. Por exemplo, calculando Ef por meio de (8.45), teremos
1,0 (1,0) 2 1 1
Ef 0,8 sen (2 21,52) E f 1,782 pu ,
1,0 sen (21,52) 2 0,7 1,0
que coincide com o resultado calculado anteriormente.
Outro problema importante em mquinas de polos salientes a determinao do ngulo

para potncia ativa mxima, denominado ngulo de estabilidade esttica, ee (veja Fig. 8.10).
Para tanto, reescrevemos a relao (8.30) na forma
P Pa sen Pb sen2 , (8.46)
onde

V1E f
Pa , (8.47)
xd
V12 1 1
Pb . (8.48)
2 xq xd
A seguir, aplicamos a regra da derivada primeira em (8.46)
P
Pa cos ee 2Pb cos2 ee 0 . (8.49)

Com uma pequena transformao trigonomtrica, a relao acima pode ser escrita como
4Pb cos2 ee Pa cos ee 2Pb 0 , (8.50)
que nada mais do que uma equao do segundo grau cuja incgnita cos . Resolvendo
(8.50), teremos
Determinao do ngulo
Pa Pa2 32 Pb2 de estabilidade esttica
cos ee . (8.51)
de uma mquina de po-
8Pb
los salientes.
A Figura 8.xx mostra a variao de ee em funo de Ef, para xd=1,0pu e vrios valores de
xq, de modo a evidenciar o efeito da salincia. Para xq =1,0 pu, por exemplo, o efeito da salincia
nulo (polos lisos), e, para o mesmo valor de Ef, quanto maior a diferena entre xd e xq, menor
ser ee.
Figura 8.XX
Variao do ngulo de estabilidade esttica, para xd =1,0pu e vrios valores de xq.

8.5 Dinmica da mquina sncrona ligada ao barramento infinito

O Barramento Infinito (BI) um sistema eltrico altamente idealizado dotado das seguintes ca-
ractersticas:
1) Tenso e frequncia constantes, independentes de qualquer perturbao externa.
2) Capacidade de absorver e fornecer qualquer potncia, seja ativa ou reativa.
3) Capacidade de absorver qualquer transitrio.
Apesar do grande grau de idealizao, o BI til da mesma forma que um plano sem atri-
to ou um condutor sem resistncia, ou seja, para facilitar a anlise inicial de sistemas mais com-
plexos. Iniciaremos analisando o caso de um nico gerador ligado ao BI e depois passaresmos a
sistemas compostos por dois geradores e finalmente para sistemas multimquinas.
8.6 Equao geral da oscilao

Seja um gerador sncrono cujo rotor cilndrico tem momento de inrcia J (kg.m2), gira com velo-
cidade aproximadamente constante s (rad / s) e est submetido a um torque total Tt=Tmec Te
(Nm), onde Tmec o torque mecnico injetado no eixo do rotor e Te o torque eletromagntico
total, decorrente do amortecimento produzido pelas barras amortecedoras e da potncia eltrica
produzida pelo estator. Sabendo que o ngulo de giro do rotor (rad), medido a partir de um
referencial fixo, podemos ento escrever
d 2
Tt J , (8.52)
dt 2
que nada mais do que a segunda lei de Newton para movimento giratrio.
Queremos escrever as equaes em funo do ngulo de carga m ec , que medido em ra-

dianos mecnicos a partir de um referencial girante. Supondo que em t=0 as referncias fixa e
girante coincidam e que o rotor gire com velocidade s (rad/s), podemos escrever:
(t ) st mec (t ) . (8.53)
Derivando (8.53) duas vezes em relao ao tempo, teremos
d 2 d 2 mec
. (8.54)
dt 2 dt 2
Substituindo (8.54) em (8.52)
d 2 mec
Tt J Tmec Te . (8.55)
dt 2

Multiplicando (8.55) pela velocidade do rotor, R , em rad/s, obtemos uma equao em

funo das potncias
d 2 mec
R J RTmec RTe .
dt 2
Fazendo
M mec R J , Pmec RTmec e Pe RTe
podemos escrever
d 2 mec
M mec Pmec Pe , (8.56)
dt 2
onde Pm ec (MW) a potncia mecnica injetada no eixo da mquina, Pe (MW) a potncia el-
trica total e M mec denominada constante de inrcia da mquina. Nas situaes de interesse
prtico a velocidade do rotor no se desvia muito da velocidade sncrona S 4f / p , pois, caso
isso acontecesse, a mquina perderia o sincronismo. Assim, podemos escrever M mec como
M mec S J , (8.57)
Sabemos tambm que a relao entre m ec , em radianos mecnicos por segundo (rad/s), e
e , em radianos eltricos por segundo (rad-e/s), para uma mquina de p polos,

p
e mec . (8.58)
2
que tambm pode ser escrita como
d e p d mec
. (8.59)
dt 2 dt
Derivando (8.59) mais uma vez, teremos
d 2 e p d 2 mec
. (8.60)
dt 2 2 dt 2
Assim, podemos escrever (8.56) em ngulos eltricos
2M mec d 2 e
Pmec Pe . (8.61)
p dt 2
Podemos definir agora duas outras constantes de inrcia da mquina. Primeiro, uma me-
constante medida em W.s2/rad
2M mec
M ' , (8.62)
p
e, considerando que a potncia base Sb (MVA), podemos tambm escrever uma constante em
pu.s2/rad

M'
M . (8.63)
Sb
Convertendo o resto da equao para pu, podemos escrever
d 2 e
M Pmec ( pu) Pe ( pu) ,
dt 2
ou, ainda, considerando que de agora em diante todas nossas equaes sero escritas em pu e que
todos os ngulos so eltricos
d 2
M Pmec Pe . (8.64)
dt 2
A potncia eletromagntica, por sua vez, pode ser escrita como
d
Pe Pd Ps ( ) , (8.65)
dt
onde Pd a potncia de amortecimento, relacionada ao amortecimento causado pelo ar no en-
treferro e pelas barras amortecedoras, e Ps a potncia sncrona. Para uma mquina de polos
lisos, Ps dada por (8.18) e, para uma mquina de polos salientes, dada por (8.31).
Substituindo (8.65) em (8.64), teremos finalmente a equao diferencial que descreve a

oscilao do ngulo quando de mudanas de potncia
Equao da oscilao
d 2 d
M 2 Pd Ps ( ) Pmec . (8.66) (swing) da mquina
dt dt sncrona.
A equao (8.66) conhecida como equao da oscilao, equao de swing ou equa-

o do balano da mquina sncrona (GRAINGER; STEVENSON, 1994, p. 701).
Alm de M e M, outra constante de inrcia bastante usada em estudos de estabilidade

simplesmente a energia cintica do rotor da mquina dividida pela potncia base, ou seja
Wk
H . (8.67)
Sb
Considerando que Wk medida em MJ e que Sb medida em MVA, a unidade de H ser
o MJ/MVA, ou simplesmente o segundo. fcil obter uma relao entre M e H. Podemos iniciar
escrevendo a energia cintica como
Jmec
2
Wk . (8.68)
2

Lembrando que M mec foi definida como Jmec , teremos
M mecmec
Wk . (8.69)
2
Logo, H ser
M m ecm ec
H , (8.70)
2 Sb
ou, lembrando da definio de M
pSb M mec pMmec
H . (8.71)
2 2 Sb 4
Finalmente, usando (8.57) para expressar a velocidade em radianos eltricos por segundo
pMe 2 Me
H . (8.72)
4 p 2
A constante H escrita na frequncia de referncia do sistema, 0 2f 0 . Logo
H f 0 M . (8.73) Constante de inrcia H.
8.7 Anlise linearizada Mquina de polos lisos

A equao (8.65) no nem um pouco simples. Por um lado, ela no linear, por causa do termo
Ps ( ) . Por outro, ela no homognea, pois o termo direita do sinal de igual no nulo. As-
sim, temos uma equao diferencial de segunda ordem, no linear e no homognea, que, no
caso geral, somente pode ser resolvida por meio de mtodos computacionais (e.g., Runge-Kutta).
Contudo, no caso de pequenas oscilaes (digamos, / 10 / 10 ), podemos considerar
que sen (em radianos). Nesse caso, a potncia sncrona para polos lisos pode ser escrita
como Ps, onde Ps o coeficiente de sincronizao, ou a inclinao da curva Ps na origem.
O resultado a seguinte equao diferencial linearizada
Equao diferencial
d 2 d
M 2 Pd Ps 0 , (8.74) linearizada e homog-
dt dt nea.

na qual Pm ec foi inicialmente igualada a zero (gerador flutuando a vazio), para obtermos uma
equao diferencial homognea, de mais fcil soluo.
A equao (8.74) matematicamente idntica equao do oscilador harmnico amorte-

cido (BRADBURY, 1968, p. 137), cuja soluo geral
(t ) Ae r t Ber t ,
1 2
(8.75)
onde A e B so constantes a determinar e r1 e r2 so as razes da equao caracterstica asso-
ciada a (8.74), que
Mr 2 Pd r Ps 0 . (8.76)
Em analogia soluo do oscilador harmnico amortecido, podemos definir o fator de
amortecimento como
Pd Fator de amortecimento.
, (8.77)
2 MPs
e a frequncia angular natural como
Ps Frequncia angular
n . (8.78) natural.
M
A equao (8.76) pode ento ser reescrita como

r 2 2n r n2 0 , (8.79)
cuja soluo geral
r n n 2 1 . (8.80)
A soluo (8.80) pode ser dividida em quatro casos, analisados a seguir.
Caso I. 0 . Nesse caso a soluo (8.80) se reduz a duas razes imaginrias conjugadas,
r jn , resultando na seguinte soluo para a equao diferencial:
(t ) Ae j t Be j t .
n n
(8.81)
Lembrando que e jnt cos(nt ) jsen (nt ) , a soluo acima pode ser reescrita como
(t ) A cos(nt ) jsen (nt ) B cos(nt ) jsen (nt ),

ou,
(t ) A B cos(nt ) j A B sen(nt ) .
Fazendo A B a0 e j A B b0 , podemos escrever
(t ) a0 cos(nt ) b0 sen(nt ) .
Fazendo ainda a0 A0 cos e b0 A0 sen , podemos finalmente escrever
(t ) A0 cos(nt ) , (8.82)
onde A0 e so constantes a determinar. A relao (8.81) a soluo para o caso do oscilador

harmnico simples.
Caso II. 0 n . Aqui as solues de (8.79) so duas razes complexas conjugadas,
r n j , onde n 1 2 . A soluo da equao diferencial pode agora ser escrita

como:
(t ) e t Ae jt Be jt .
n
(8.83)
Usando os mesmos mtodos do Caso I, podemos escrever
Soluo para oscilador

Aoe t cos(t ) ,
n
(8.84) com amortecimento sub-
crtico.
que uma soluo do tipo oscilador com amortecimento subcrtico. Para 0 , o Caso II se
reduz ao Caso I.
Caso III. n . Aqui as solues de (8.79) so reais e idnticas, r , resultando em uma

soluo criticamente amortecida, sem oscilao.
Aoe t . n
(8.85)
Caso IV. n . Aqui as solues de (8.78) so reais e diferentes, r n n 2 1 , resul-

tando em uma soluo com amortecimento supercrtico, tambm sem oscilao.
n t n 2 1 t 2 1 t

Ao e e

e n
.

ou,

2 1 t
1
Ao 1 e t e
n
n
. (8.86)
e
No se espera solues criticamente ou supercriticamente amortecidas em mquinas sn-
cronas, pois o amortecimento nelas observado no expressivo a esse ponto.
Exemplo 8.3. Um gerador sncrono est flutuando a vazio, com velocidade de oscilao nula,
quando submetido a uma perturbao momentnea o . Considerando que a oscilao rapi-
damente amortecida at voltar situao inicial, determine as constantes e Ao .
Soluo. Em t 0 , devemos ter

(t 0) Ao cos( ) o ,
ou,
o Constante Ao para gera-

Ao . (8.87) dor flutuando a vazio,
cos( ) Casos I e II.
A velocidade de oscilao de d / dt , ou
Aoe t sen(t ) n cos(t ).

n
(8.88)
Velocidade de oscilao
para os Casos I e II.
Sabendo que a velocidade em t 0 , nula,

Ao sen( ) n cos( ) 0 ,
ou,
atg n / . (8.89)
ngulo de fase para os
Casos I e II.
A soluo (8.84) mostra que, para o gerador flutuando a vazio, qualquer perturbao ser
amortecida. Caso isso no ocorra, corremos o risco de que o gerador saia de sincronismo, especi-

almente no caso de potncias mdias muito baixas, quando o ngulo de estabilidade esttica ee
poder ser atingido rapidamente.
Quando Pmec for injetada no eixo do gerador, o ngulo de regime permanente, quando
t , ser . Agora a soluo geral para a equao (8.83) pode ser escrita como
Soluo geral para o

(t ) A0e t cos(t ) ,
n
(8.90) gerador linearizado
amortecido.
que pode ser considerada como a soluo geral linearizada para um gerador flutuando a vazio
( 0 ) ou sumetido a um degrau de potncia mecnica tal que (t ) . A soluo
(8.89) considera que a potncia mecnica subitamente aplicada ao eixo do gerador, em forma
da degrau. Na prtica isso impossvel, por causa da grande inrcia do conjunto turbina-gerador,
mas o que nos interessa aqui a facilidade da modelagem matemtica. A aplicao de potncia
mecnica em forma de rampa, por exemplo, poderia ser modelada como uma grande sequncia
de degraus, na qual o ngulo de uma sequncia entendido como o ngulo 0 da prxima
sequncia. Outros tipos de variao de potncia mecnica so muito difceis, ou mesmo imposs-
veis, de serem resolvidos analiticamente e devero ser objeto de soluo numrica.
Exemplo 8.4. Repita o Exemplo 8.3 para velocidade inicial nula, mas agora com 0 .
Soluo. Em t 0 , devemos ter

(t 0) Ao cos( ) o ,
ou,
o Constante Ao generali-
Ao . (8.91) zada para a oscilao do
cos( ) gerador, Casos I e II.
Os clculos mostram que o ngulo continua sendo dado pela relao (8.89).

Exemplo 8.5. Seja um gerador sncrono de 28 polos, 10 MW, 60 Hz, H=8 MJ/MVA, Ps=20
pu/rad-e, Pd=0,20 pu/rad-e/s, ligado diretamente ao barramento infinito. Escreva a equao da
oscilao do ngulo para: (a) gerador flutuando a vazio, com ngulo inicial de 5 graus eltricos
e ngulo final nulo; (b) carga nominal subitamente aplicada ao eixo do gerador.
Soluo. A equao diferencial a vazio pode ser escrita como
d 2 d
M 2
Pd Ps 0
dt dt
H
A partir dos dados do problema, e considerando que M , teremos
f
d 2 d
0,0424 2 0,20 20 0 (8.92)
dt dt
A frequncia natural de oscilao ser
Ps 20
n 21,708 rad/s. (8.93)
M 0,0424
De acordo com a relao (8.76), o fator de amortecimento ser
Pd 0,20
0,1085 , (8.94)
2 MPs 2 0,0424 20
correspondendo seguinte frequncia angular amortecida
n 1 2 21,71 1 0,1085 2 21,5798 rad/s. (8.95)
O ngulo de fase, de acordo com (8.88), ser
atg n / atg 0,1085 21,7080 / 21,5798 6,23 . (8.96)

A constante Ao ser
o 5
Ao 5,03 . (8.97)
cos( ) cos(6,23)
A resposta para o gerador flutuando a vazio ser, em graus eltricos
(t ) 5,03 e 2,3562t cos(21,5798t 6,23) . (8.98)

A Figura 8.12 ilustra a variao de em relao ao tempo e a Figura 8.13 ilustra a velo-
cidade de em relao ao tempo.

A Figura 8.14 ilustra a evoluo do sistema no chamado espao de fase. Nesse caso, abs-
tramos a varivel tempo e plotamos a velocidade de em funo de . Percebemos que o siste-
ma inicia com =5 e velocidade nula, atingindo a estabilidade no ponto (0, 0).
Figura 8.12
Evoluo do ngulo para o gerador
linearizado flutuando a vazio
Figura 8.13
Velocidade de para o gerador
linearizado flutuando a vazio
Figura 8.14
Diagrama no espao de fase para o
gerador linearizado flutuando a vazio

No caso do gerador que subitamente recebe carga nominal, podemos considerar Pmec=1
pu. A aproximao da anlise linearizada nos permite escrever Pmec Ps , onde o ngulo ex-
presso em radianos eltricos. Assim
Pmec 1,0
0,05 rad, ou 2,87 (8.99)
Ps 20,0
A soluo ser dada agora pela equao (8.89)
(t ) A0e t cos(t ) 2,87 A0e2,3562t cos(21,5798t 6,23) ,
n
(8.100)
onde a constante Ao dada por (8.91)
o 5 2,87
Ao 2,143 (8.101)
cos( ) cos(6,23)
As Figuras 8.15, 8.16 e 8.17 mostram o comportamento do gerador linearizado sob carga.
Percebemos que agora o gerador atinge a estabilidade no ponto (2,87, 0/s). Por causa das limi-
taes do mtodo linearizado, que se restringe a pequenas oscilaes, no sabemos se a potncia
Pmec ser de fato absorvida pelo gerador e se a estabilidade ser atingida em . O mtodo das
reas, visto no item 8.8 a seguir, permitir tal estimativa.
Um atrator um ponto ou conjunto de pontos para o qual evolui um sistema dinmico,

independentemente do ponto de partida. Assim, no caso do gerador a vazio o atrator (0, 0/s),
enquanto no caso do gerador sob carga o atrator (15, 0/s). Em sistemas de potncia todo o
esforo feito para que o atrator seja um nico ponto por gerador. Caso o atrator seja um conjun-
to de pontos, o que pode acontecer se o amortecimento e/ou a potncia sincronizante forem mui-
to baixos, o sistema pode se tornar instvel.
Figura 8.15
Evoluo do ngulo para o
gerador linearizado sob carga

Figura 8.16
Velocidade de para o gerador
linearizado sob carga
Figura 8.17
Diagrama no espao de fase para
o gerador linearizado sob carga
8.8 Mtodo das reas

O mtodo das reas se destina determinao da estabilidade angular de geradores de polos lisos
ou de polos salientes. No um mtodo que permita determinar a evoluo do ngulo de carga,
mas sim um mtodo simplificado que permite determinar o pior caso em determinada situao. O
mtodo das reas se baseia nas seguintes hipteses simplificadoras:
A potncia mecnica se mantm constante durante o perodo de variao angular.
A excitao tambm se mantm constante.
O amortecimento mecnico e as resistncias eltricas so consideradas desprezveis.
Feitas tais simplificaes, a equao da oscilao (8.6) pode ser escrita como
d 2
M Pe Pmec , (8.102)
dt 2

ou,
d 2 d
M 2
M Pmec Pe Pa , (8.103)
dt dt
onde Pa Pmec Pe a potncia acelerante da mquina. Multiplicando (8.26) por
2 2 d
, teremos
M M dt
d P d
2 2 a , (8.104)
dt M dt
ou, ainda
d 2 P d
2 a , (8.105)
dt M dt
Multiplicando (8.104) por dt e integrando em relao a , vem
2 m
2
M 0
Pa d , (8.106)
onde m o ngulo mximo que o rotor atinge durante a oscilao. Finalmente,
2 m

M 0
Pa d , (8.107)
onde a velocidade angular de oscilao do rotor (no confundir com a velocidade angular
sncrona). Considerando que 0 na situao de estabilidade, podemos escrever
m m
Pa d
0
Pmec Pe d 0 , (8.108)
0
ou
m m Equao geral do mto-

0
Pmecd Pe d .
0
(8.109) do das reas.
A integral da potncia em relao ao ngulo de carga igual rea debaixo de determi-

nada poro da curva P( ) . A equao (8.109) representa o mtodo das reas. Em princpio, tal
mtodo pode ser aplicado tanto para mquinas de polos lisos quanto de polos salientes. Iniciare-
mos com o caso dos polos lisos, por facilidade matemtica.

8.8.1 Tomada sbita de carga

Considere um gerador de polos lisos de 650 kW funcionando com ngulo inicial de 10, conec-
tado rede. O problema determinar a mxima potncia mecnica que pode ser aplicada ao eixo
do gerador sem perda de sincronismo.
Conforme mostrado na Fig. 8.18, o rotor poder oscilar at m e atingir a esta-

bilidade no ponto D. A rea embaixo de Pmec ser ento

Pmecd Pmec 0 , (8. 110)
0
que corresponde rea do retngulo ACFGm cujas arestas so arestas Pmec e 0 .
Figura 8.18
rea embaixo de Pmec
A rea embaixo da senoide Pe mostrada na Figura 8.19 e pode ser calculada como
Pmaxsen d Pmax cos 0 Pmax cos 0 cos ,

(8. 111)
0
ou

Pmax sen d Pmax cos 0 cos , (8. 112)
0

Pmec 0 Pmax cos 0 cos . (8. 113)
Alm disso, no ponto de estabilidade D, teremos Pmec Pmaxsen . Logo
Pmax sen 0 Pmax cos 0 cos . (8. 114)

Simplificando e reorganizando um pouco

cos o cos o sen 0 . (8.115)

Equao para tomada
sbita de carga.
Note agora que quando confrontamos as Figuras 8.18 e 8.19, a rea ABDFH comum
aos dois grficos e pode ser eliminada. O problema torna-se assim encontrar o ngulo que
iguala as reas A1 e A2, conforme ilustrado na Figura 8.20. Durante a rea A1 o rotor acelera, pois
a potncia mecnica maior do que a potncia eltrica. Da mesma forma, durante a rea A2 o
rotor desacelera, pois a potncia eltrica torna-se maior do que a potncia mecnica. O rotor
atinge a estabilidade no ponto D, parando de oscilar. A equao (8.109) pode ento ser reescrita
como
m Equao geral do mto-

0
Pa d Pdesd ,

(8.116) do das reas, em funo
de Pa e Pdes.
onde Pdes a potncia desacelerante.
De maneira mais simples, a equao (8.116) pode ser escrita como
Equao geral do mto-

A1 A2 , (8.117) do das reas, em funo
de A1 e A2.
onde A1 e A2 so definidas na Figura 8.20.

Fica como exerccio encontrar a equao (8.115) por meio de (8.116).

Figura 8.19
reas embaixo da senoide Pe
Figura 8.20
Dinmica para tomada sbita de carga
A equao (8.115) no tem soluo analtica, mas pode ser resolvida por tentativas ou por
meio do solver da HP, de alguma planilha numrica ou de um software numrico. Para
0 10 / 18 rad-e, teremos
51,0
Um aspecto interessante desta soluo a independncia dela em relao a Pmax 650 .
Em outras palavras, o resultado acima vale para qualquer mquina de polos lisos funcionando
com 0 10 , independente da potncia nominal, do nmero de polos e da frequncia, respeita-
das as simplificaes do mtodo das reas.
A potncia mecnica na situao de estabilidade ser

Pmec 650 sen (51,0) 505,2 kW . (8. 118)

Outro aspecto interessante da soluo que, apesar da no-linearidade da equao (8.8), a
relao entre 0 e uma reta ascendente o=83,67, caindo a seguir, conforme mostra a Fig.
8.21. No razovel esperar que um gerador de polos lisos esteja operando com o>80 quando
de um aumento subido de carga e, assim, a relao entre 0 e deve ser sempre uma reta as-
cendente.
Figura 8.21
Soluo para tomada sbita de carga
A tomada sbita de carga um caso muito idealizado, pois, nos casos reais, a potncia
mecnica nunca varia como um degrau. Alm disso, a rigor no estamos levando em conta um
dos requisitos do mtodo das reas, que a prpria constncia da potncia mecnica. A seguir
discutimos um caso um pouco mais realista.
8.8.2 ngulo crtico para um gerador e uma linha

A Figura 8.22 ilustra um gerador alimentando um barramento por meio de um transformador e
de uma nica linha. Em dado instante, quando t=t0 e = 0 ocorre uma falta trifsica slida na
linha, que interrompida em ambas as extremidades por disjuntores religadores. Desejamos cal-
cular o ngulo crtico c de excurso do gerador para que os disjuntores religuem a linha sem
que o gerador perca a estabilidade.
Figura 8.22
Circuito simples para clculo do ngulo crtico

A dinmica do problema ilustrada na Figura 8.23. A linha religada em c. Logo, entre

0 e c a potncia eltrica entregue pelo gerador linha nula e, como Pmec>Pe, o rotor acelara.
Entre c e max a linha religada e, como Pmec<Pe, o rotor desacelara. Se A1<A2, o gerador man-
ter a estabilidade quando do religamento. Caso contrrio, a estabilidade no poder ser mantida
e o gerador no poder ser religado.
Figura 8.23
Dinmica para um gerador e uma linha, falta franca
A rea A1 calculada como

A1 Pmec c 0 , (8. 119)
enquanto a rea A2 calculada como
Pmaxsen Pmec d Pmax cos 0 cos c Pmec 0 c ,

max
A2
c
ou,
A2 Pmax cos 0 cos c Pmec 0 c . (8. 120)
Pmec c 0 Pmax cos 0 cos c Pmec 0 c ,
ou,
Pmec 2 0 Pmax cos 0 cos c . (8.121)
Em 0, podemos escrever
Pmec Pmax sen 0 , (8.122)
onde Pmax a potncia mxima antes da falta.
Substituindo (8.122) em (8.120), teremos

cos c 2 0 sen 0 cos 0 . (8.123)

ngulo crtico para uma
linha, falta franca.
A equao (8.123) pode ser resolvida analiticamente: tendo 0, podemos calcular c fa-
cilmente. O resultado ser maior do que o valor real, pois estamos desprezando o amortecimento.
Entretanto, o raciocnio acima s ser vlido caso a falta seja franca, ou seja, caso a transferncia
de potncia entre as barras 2 e 3 seja nula entre 0 e c. Caso isso no acontea, teremos a situa-
o ilustrada na Figura 8.24, na qual r1 definido como:
Pfal
r1 , (8.124)
Pmax
onde Pfal a potncia mxima durante a falta.
A rea A2 continua a mesma calculada pela relao (8.119). A rea A1 , por outro lado,
deve ser calculada como:
c
A1 Pmec r1Pmax sen d c 0 Pmec r1 Pmax cos c cos 0 .
0
Fazendo Pmec Pmax sen 0 , teremos:
A1 Pmax c 0 sen 0 r1 cos c cos 0 . (8.125)

Igualando (8.125) e (8.119), vem
c 0 sen 0 r1 cos c cos 0 cos 0 cos c 0 c sen 0 ,
ou, finalmente
cos c
2 0 sen 0 r1 1 cos 0 ngulo crtico para uma
. (8.126) linha, falta no franca.
1 r1
fcil perceber que (8.126) se reduz a (8.123) quando r1 0 .

Figura 8.24
Dinmica para um gerador e uma linha, falta no franca
Exemplo 8.6. Para o sistema da Figura 8.22, considere que a potncia mecnica injetada no eixo
do gerador 0,5pu, a reatncia sncrona do gerador 0,8pu, a reatncia do transformador 0,2pu
e a reatncia total da linha de transmisso, cuja resistncia desprezvel, 0,4pu. A excitao do
gerador 1,05pu e a tenso no barramento 1,0pu. Em dado momento, uma falta ocorre nas
imediaes da barra 2. Determine o ngulo crtico de extino da falta supondo que: (a) a falta
trifsica e franca; (b) a falta trifsica e tem reatncia de 0,05pu, com resistncia desprezvel.
Soluo. Lorem ipsum lorem ipsum

2 41,81 sen (41,81) 0,0 1cos(41,81)
cos c 180
0,3761 ,
1 0,0
c 67,91 .
Lorem ipsum lorem ipsum

2 41,81 sen (41,81) 0,7778 1cos(41,81)
cos c 180
0,9166 ,
1 0,7778
c 156,43 .
Lorem ipsum lorem ipsum
A Figura 8.25 mostra a variao do ngulo crtico em funo do ngulo inicial para vrios
casos de r1. interessante notar que, quanto maior o valor de r1, ou seja, quanto maior a potncia
mxima durante a falta em relao potncia mxima antes da falta, mais restrita ser a variao

do ngulo inicial. Note tambm que, a partir de o=75, c praticamente igual a o. Logo, no
haver tempo para que o ngulo de carga se estabelea em um novo valor antes de se tornar crti-
co. Para ngulos inciais menores do que 75, c pode atingir maiores valores em relao a o para
maiores impedncias de falta. Para o=50, por exemplo, c=64 para r1=0 e c=115 para r1=0,8.
Assim, faltas francas restrigem mais a excurso de antes que este se torne crtico. Esta e outras
razes tornam necessria a definio de uma margem de segurana para o ngulo de carga.
Figura 8.25
ngulos crticos para vrios casos de r1.
8.8.3 ngulo crtico para um gerador e duas linhas

O circuito da Figura 8.23 do tipo N0, ou seja, no se pode perder nenhum elemento do sis-
tema de transmisso sem que a carga (nesse caso, o BI) perca a alimentao. O circuito da Figura
8.26, por outro lado, do tipo N1: pode-se perder at um elemento do sistema de transmisso
sem que a carga perca a alimentao, embora a potncia deva ser reduzida aps o evento.
Figura 8.26
Circuito duplo para clculo do ngulo crtico
A falta trifsica e ocorre na linha de baixo, podendo ser franca ou no, e extinta pela
abertura ou pela reposio da linha em funcionamento (veremos que a matemtica da situao
indiferente, variando apenas os parmetros).

A Figura 8.27 mostra o circuito reduzida a reatncias, no qual a falta representada tam-
bm por uma reatncia. O gerador representado como uma reatncia atrs de uma tenso de
excitao Ef cujo ngulo de carga em relao ao barramento infinito representado pela barra
3. A falta ocorre no ponto a: se a=0, a falta ocorre na barra 2; se a=1, a falta ocorre na barra 3; se
a=1/2, a falta ocorre na metade da linha de baixo e assim por diante.
Figura 8.27
Circuito reduzido a reatncias para clculo do ngulo crtico
A Figura 8.28 ilustra a dinmica da situao. O parmetro r1 o mesmo dado por (8.124)
e r2 dado por
Ppos
r2 , (8.127)
Pmax
onde Ppos a potncia mxima aps a falta.
No h nada que impea termos r2=1. Nesse caso, a falta seria eliminada por meio da re-
posio da linha de baixo e a situao pr-falta seria idntica situao ps-falta. Contudo, no
podemos ter r2=r1, pois isso significaria que a linha de baixo permanece em falta.
A falta ocorre em o, quando o sistema deixa de funcionar em Pmax sen e passa a funcio-
nar em r1 Pmax sen . A potncia mecnica torna-se superior potncia eltrica, o gerador comea
a acelerar e o ngulo comea a aumentar. Quando o ngulo chega c, os disjuntores do sistema

atuam e a falta extinta. O sistema passa ento a funcionar em r2 Pmax sen , a potncia eltrica
torna-se maior do que a potncia mecnica e o gerador passa a desacelerar. Contudo, por causa
da inrcia do sistema, o ngulo pode atingir o valor mximo max antes de comear a diminuir e
se estabelecer em c. Tal situao estvel, contudo, s ser atingida caso o critrio das reas (8.
117) seja respeitado. Caso contrrio, o sistema perder estabilidade e o gerador dever ser des-
conectado.

Observando a Figura 8.28, a rea A1, correspondente acelerao do rotor, igual a rea
A1 do problema para uma linha e impedncia de falta no nula, cuja demonstrao repetida a
seguir por convenincia.
c
A1 Pmec r1Pmaxsen d c 0 Pmec r1Pmax cos c cos 0 . (8.128)
0
A rea A2, correspondente potncia de desacelerao, ser

max
A2
c
r2 Pmaxsen Pmec d r2 Pmax cos c cos max c max Pmec . (8.129)
Igualando A1 (8.128) e A2 (8.129), vem

c 0 Pmec r1Pmax cos c cos 0 r2 Pmax cos c cos max c max Pmec ,
ou
r1Pmax cos c r2 Pmax cos c r2 Pmax cos max Pmec max 0 Pmec r1Pmax cos 0 ,
ou, ainda
r2 r1 Pmax cos c max 0 Pmec r1Pmax cos 0 r2 Pmax cos max .
Finalmente, fazendo Pmec Pmax sen 0 , vem
cos c
max 0 sen 0 r2 cos max r1 cos 0 ngulo crtico para duas
. (8.130) linhas, com falta em uma
r2 r1 das linhas.
Figura 8.28
Dinmica para um gerador e duas linhas
A Figura 8.29 mostra a variao do ngulo crtico em funo do ngulo inicial para r2=1
e vrios valores de r1. Note que esta figura exatamente igual Figura 8.25, que ilustrava o

comportamento para uma s linha. Assim, se, no caso do sistema com duas linhas, se uma das
linhas entrar em falta e voltar integralmente aps a extino desta, fazendo r2=1, a dinmica do
sistema exatemente mesma do sistema com uma s linha, para o mesmo valor de r1.
Figura 8.29
ngulos crticos para r2=1 e vrios casos de r1.
A Figura 8.30 mostra as mesma curvas anteriores, mas agora para r2=0,85. Note que a
curva para r1=0,9 no mais mostrada, pois no podemos ter r1> r2. Note tambm que a variao
do ngulo crtico agora bem mais restrita. Por exemplo, antes (Figura 8.29), para o=40, o
ngulo crtico era aproximadamente 75 para r1=0,2. Agora, com r2=0,75 e o mesmo valor para
r1, o ngulo crtico aproximadamente 62. Consequentemente, o sistema de proteo deve ser
muito mais rpido na tarefa de extino da falta.
A Figura 8.31 mostra as mesma curvas anteriores, mas agora para r2=0,75, um valor j
bastante factvel, como veremos no exemplo a seguir. Agora, para r1=0,2 e o=40, o ngulo cr-
tico cai para aproximadamente 48.
Note, nas Figuras 8.30 e 8.31, que agora h valores mximos para o ngulo inicial para
manter a estabilidade do sistema. Por exemplo, tanto para r2=0,85 quanto para r2=0,75, se o ge-
rador tiver o>60, o sistema ser levado instabilidade para faltas com quaisquer valores de r1.
Assim, no importa quo rpido for o sistema de proteo, o gerador dever ser desconectado da
rede.

Figura 8.30
ngulos crticos para r2=0,85 e vrios casos de r1.
Figura 8.31
ngulos crticos para r2=0,75 e vrios casos de r1.
Exemplo 8.7. Para o sistema da Figura 8.26, considere que a potncia mecnica injetada no eixo
do gerador 0,8pu, a reatncia sncrona do gerador 0,2pu, a reatncia do transformador 0,2pu
e a reatncia total de cada linha de transmisso, cujas resistncias so desprezveis, 0,4pu. A
tenso de excitao do gerador 1,05pu e a tenso no barramento 1,0pu. Em dado momento,
uma falta ocorre exatamente na metade da linha de baixo e extinta por meio da abertura desta.
Determine o ngulo crtico de extino da falta, supondo que esta trifsica e tem reatncia de
0,02pu, com resistncia desprezvel.
Soluo. Antes da falta o circuito de impedncias pode ser desenhado como na Figura 8.32. A
reatncia antes da falta pode ser facilmente calculada como:

x pre x1'3 xd ' xT xLT 1 // xLT 2 , (8.131)

ou
x pre j 0,2 j 0,2 j 0,4 // j 0,4 j 0,6 pu . (8.132)
Figura 8.32
Circuito de impedncias antes da falta
Depois da falta a linha de transmisso de baixo desligada e o circuito de impedncias

representado na Figura 8.33. A reatncia depois da falta tambm pode ser facilmente calculada:
x pos x1'3 xd ' xT xLT 1 , (8.133)
ou
x pos j 0,2 j 0,2 j 0,4 j 0,8 pu . (8.134)
Figura 8.33
Circuito de impedncias depois da falta
Durante a falta, com a insero da reatncia de falta, o circuito aquele mostrado na Fi-
gura 8.34 e o clculo da reatncia entre as barras 1 e 3 torna-se um pouco mais complicado.
Uma maneira de desenvolver uma expresso geral para x fal x1'3 , correspondente reatncia
entre as barras 1 e 3 durante a falta, transformando o circuito da Figura 8.34 em um circuito

delta entre as barras 1, 3 e o terra.

Figura 8.34
Circuito de impedncias durante a falta
Podemos comear transformando o delta entre as barras 2, 3 e m em um estrela, como in-

dicado na figura. Aplicando a transformao delta-estrela, teremos:
xLT 1axLT 2 ax x
x2 LT 1 LT 2 , (8.135)
xLT 1 axLT 2 (1 a) xLT 2 xLT 1 xLT 2
(1 a) xLT 1 xLT 2
x3 , (8.136)
xLT 1 xLT 2
a(1 a) x 2LT 2
xm . (8.136)
xLT 1 xLT 2
O circuito transformado para estrela representado na Figura 8.35.
Figura 8.35
Circuito da Figura 8.34 transformado para estrela
O circuito deve ser agora transformado de estrela para delta. Apenas a reatncia x1'3 de
interesse, a qual ser igual reatncia durante a falta, x fal :
x1' x3 x3 xn 0 xn 0 x1'
x fal , (8.137)
xn 0
xx
x fal 1' 3 x3 x1' , (8.138)
xn 0

axLT 1 xLT 2 (1 a) xLT 1 xLT 2

xd ' xT
x x LT 2 x LT 1 x LT 2 (1 a) xLT 1 xLT 2 axLT 1 xLT 2 x ' x
x fal LT 1
T, (8.139)
a(1 a) x LT 2 xLT 1 xLT 2
2 d
xf
xLT 1 xLT 2
xLT 1 xLT 2
xd ' xT a (1 a) xLT 1 xLT 2
xLT 1 xLT 2 x x (8.140)
x fal LT 1 LT 2 xd ' xT .
x f ( xLT 1 xLT 2 ) a(1 a) x LT 2
2
xLT 1 xLT 2
Lembrando-nos que xLT 1 xLT 2 /( xLT 1 xLT 2 ) xLT 1 // xLT 2 , e que, de acordo com (8.131),
xd ' xT xLT 1 // xLT 2 x pre , a relao acima pode ser escrita de forma mais compacta. O resultado
suficientemente interessante e geral para ser colocado em destaque:
Reatncia durante a
( x ' x axLT 1 // xLT 2 )(1 a) xLT 1 xLT 2
x fal d T x pre . (8.141) falta para casos seme-
x f ( xLT 1 xLT 2 ) a(1 a) x 2LT 2 lhantes ao do Exemplo
8.7.
No presente exemplo, teremos:
x fal j
0,2 0,2 0,2 a (1 a) 0,4 0,4 j0,6 .
(8.142)
0,02 (0,4 0,4) a(1 a) 0,16
No caso de uma falta na metade da linha de baixo, teremos a=0,5. Assim:
x fal j
0,4 0,2 0,5 0,5 0,16 j0,6 j1,3143 . (8.143)
0,02 0,8 0,5 0,5 0,16
Os prximos passos so calcular as potncias mximas antes, durante e depois da falta:
V1Eg 1,0 1,05
Pmax 1,75 pu . (8.144)
x pre 0,6
V1Eg 1,0 1,05
Pfal 0,7989 pu (8.145)
x fal 1,3143 .
V1Eg 1,0 1,05
Ppos 1,3125 pu (8.145)
x fal 0,8 .
Os parmetros r1 e r2 so calculados pelas relaes (8.124) e (8.127), respectivamente:
Pfal 0,7989
r1 0,4565 , (8.146)
Pmax 1,75
Ppos 1,3125
r2 0,75 . (8.147)
Pmax 1,75

O prximo passo antes de calcularmos o ngulo crtico calcularmos os ngulos o e
max . Analisando a Figura 8.28, podemos escrever:
P ngulo inicial para o

o asen m ec , (8.148) caso de duas linhas.
Pmax
P ngulo mximo para o

max asen mec . (8.149) caso de duas linhas.
r2 Pmax
Aplicando (8.148) e (8.149), teremos, respectivamente:

0,8
o asen 0,4748 rad 27,2 , (8.150)
1,75
0,8
max asen 2,4861 rad 142,44 . (8.151)
0,75 1,75
Finalmente, aplicando a relao (8.130), com os clculos feitos em radianos, vem:
cos c
2,4861 0,4748 sen(0,4748) 0,75 cos(2,4861) 0,4565 cos(0,4748) ,
(8.152)
0,75 0,4565
c 1,8509 rad , (8.153)
c 106,05 (8.154)
A reatncia xfal, dada por (8.141), pode ser plotada em funo de xf, para vrios valores de
a. O resultado mostrado na Figura 8.32. Como pode ser visto, para faltas na barra 2 (a=0), a
reatncia durante a falta xfal varia bastante em funo de xf (quanto mais franca a falta, maior o
valor de xfal e, consequentemente, menor a potncia transferida durante a falta). J quanto mais
perto da barra 3 for a falta, menos dependente de xf torna-se xfal. Para faltas na barra 3, teremos
a=1 e xfal = xpre, indicando que a falta foi absorvida pelo barramento infinito, no tendo sido per-
cebida pelo gerador na barra 1.

Figura 8.32
Variao de xfal em funo de xf, para vrios valores de a
Lorem ipsum...

9. ESTABILIDADE EM SISTEMAS
MULTIMQUINAS
9.1 Introduo
Lorem ipsum...
Em breve...

10. ESTABILIDADE TRANSITRIA NO DOMNIO

DO TEMPO
10.3 Introduo
Em breve...por enquanto, s rabiscado em: http://goo.gl/Kwz0sf.

11. ESTABILIDADE DE TENSO
11.1 Introduo
Lorem ipsum lorem ipsum lorem ipsum lorem ipsum lorem ipsum lorem ipsum lorem ipsum
lorem ipsum lorem ipsum lorem ipsum...
12. OPERAO ECONMICA DE SISTEMAS
12.1 Introduo
Lorem ipsum lorem ipsum lorem ipsum lorem ipsum lorem ipsum lorem ipsum lorem ipsum
lorem ipsum lorem ipsum lorem ipsum...

13. REFERNCIAS
BALDWIN, Thomas. Power system analysis weekly outline, lecture notes and assignments,
2004. Disponvel em: <http://www.eng.fsu.edu/~tbaldwin/eel4213/public/index.html#Outline>.
BARBOSA, F. Maciel. Estabilidade de sistemas elctricos de energia. Apostila, 2007. Dispo-

nvel em: < http://paginas.fe.up.pt/~fmb/RTSEE/Estabilidade%20%20Janeiro%2007.pdf>
BICHELS, Arlei. Sistemas eltricos de potncia mtodos de anlise e soluo. Curitiba:

CEFET-PR, 3 ed., 1998.
BRADBURY, T. C. (Ted Clay). Theoretical mechanics. New York: J. Wiley, c1968. 641 p.
BLACKBURN, J. Lewis. Symmetrical components for power systems engineering. New

York: Marcel Dekker, Inc., 1993, 427 p.
BORGES, Carmen Lucia Tancredo. Anlise de sistemas de potncia. EE-UFRJ, Departamento

de Eletrotcnica, 2005. Disponvel em: <http://www.dee.ufrj.br/lasp/ASP/Apostila%20
Completa.pdf>.
COSTA, Antonio J.A.; SILVA, A.S. Dinmica e controle de sistemas de potncia. Apostila.
Disponvel em: <http://www.labspot.ufsc.br/~simoes/dincont/index.html>
ELETROBRS. Plano nacional de energia eltrica 1993-2015 Plano 2015, 1993.
ELETROBRS. A energia eltrica no Brasil (da primeira lmpada Eletrobrs). Biblioteca do

Exrcito Editora, 1977, 244 p.
ELGERD, Olle Ingemar; VAN DER PUIJE, Patrick D. Electric power enginnering. 2nd ed.
New York: Chapman & Hall, c1998, 449 p.
EL-HAWARY (2008). Introduction to Electrical Power Systems. Wiley-IEEE Press, 2008.
EMAD, Fawzi Philip. ENEE 474 power systems, Department of Electrical and Computer Engi-
neering, University of Maryland, 2004. Disponvel em: <http://www.ece.umd.edu/class/
enee474.F2003/index.htm> .
FITZGERALD, A.E.; KINGSLEY JR., Charles; KUSKO, Alexander. Mquinas eltricas,

McGraw-Hill do Brasil, 1975.
GMEZ-EXPSITO, A.; CONEJO, A.J.; CAIZARES, C. (ed.). Sistemas de energia eltrica

anlise e operao, LTC: Rio de Janeiro, 2011.
GLOVER, J. Duncan; SARMA, Mulukutla S.; OVERBYE, Thomas. Power systems analysis
and design. CL-Engineering, 2007, 768 p.
GROSS, Charles A. Power system analysis. 2nd ed. New York; Chichester: J. Wiley, c1986, 593
p.

JORDO, Rubens Guedes. Mquinas sncronas, Rio de Janeiro: LTC, 2013, 2 ed., 224 p.
KIMBARK, Edward. Power system stability. New York: Wiley-IEEE Press, 1995, 998 p., 3v.
KINDERMANN, Geraldo. Curto-circuito. 4. ed. Florianpolis, SC: UFSC, 2007, 233 p.
KOCHOLIK, Aline; INCOTE, Mrcia C. M. Simulador fasorial para anlise do comporta-

mento do gerador sncrono de polos salientes conectado em barramento infinito operando
em regime permanente. 2009. 139 f. TCC. Departamento Acadmico de Eletrotcnica
DAELT, Universidade Tecnolgica Federal do Paran. 2009. Disponvel em: < http://www.
daelt.ct.utfpr.edu.br/professores/andreal/TCCAlineMarcia.pdf>.
KOTHARI, D.P.; NAGRATH, I.J. Power systems engineering. New Delhi: Tata McGraw-Hill.
2008, 1024p. Disponvel parcialmente em: < http://books.google.com.br/books?id=
yPPYGD88UR0C&printsec=frontcover>.
KUNDUR, Prabha. Power system stability and control. New York, US: McGraw-Hill, 1994.
xxiii, 1176p.
MELLO, F. P. DE; ELETROBRS. Dinmica das mquinas eltricas - I. Santa Maria: Editora
da UFSM: Rio de Janeiro: Eletrobrs, 1979. 223 p. (Curso de engenharia em sistemas eltricos
de potncia.Srie P.T.I.;v.4).
MELLO, F. P. DE; ELETROBRS. Dinmica das mquinas eltricas - II. Santa Maria: Edito-
ra da UFSM: Rio de Janeiro: Eletrobrs, 1979. 103p. (Curso de engenharia em sistemas eltricos
de potncia.Srie P.T.I.;v.5)
MONTICELLI, Alcir. Fluxo de carga em redes de energia eltrica. So Paulo: Edgard Bl-
cher, 1983.
MONTICELLI, Alcir; GARCIA, Ariovaldo. Introduo a sistemas de energia eltrica. Cam-

pinas: Editora da Unicamp, 2003, 251 p.
ROBBA, Ernesto Joo. Introduo a sistemas eltricos de potncia: componentes simtricas.

2. ed. So Paulo: Edgar Blucher, 2000, 467p.
GRAINGER, John J.; STEVENSON, William D. Power system analysis. McGraw-Hill, 1994,
787 p.
WEEDY, B. M. Sistemas eltricos de potncia. So Paulo: Editora da Universidade de So

Paulo: Polgono, 1973. 363p.
ZANETTA Jnior., Luiz Cera. Fundamentos de sistemas eltricos de potncia, Editora Livra-
ria da Fsica, 2005. Disponvel parcialmente em: <http://books.google.com.br/
books?id=Nyep7l_kBFkC&printsec=frontcover> .

Apostila Curto Circuito - Arthur Damaso PDF

Enviado por

Direitos autorais:

Formatos disponíveis

Apostila Curto Circuito - Arthur Damaso PDF

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Apostila Curto Circuito - Arthur Damaso PDF

Enviado por

Direitos autorais:

Formatos disponíveis

SISTEMAS ELTRICOS

Prof. Alvaro Augusto W. de Almeida

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

8.4 MQUINA DE POLOS SALIENTES EM REGIME PERMANENTE ............................................................................ 151

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

Durante toda a existncia do nosso curso de Engenharia, a disciplina de Sistemas Eltri-

A disciplina Sistemas de Potncia 1 da grade atual tem basicamente a mesma ementa da

Longe de pretenderem substituir a literatura existente na rea, as presentes Notas de Aula

Prof. Alvaro Augusto W. de Almeida

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

Glossrio de smbolos gerais

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

A necessidade de alimentar os grandes centros consumidores do Sudeste a partir da ener-

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

O planejamento e operao de sistemas eltricos, do porte do sistema brasileiro ou no,

Fluxo de potncia, tambm conhecido como fluxo de carga, um problema matemti-

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

Fluxo de potncia e curto-circuito formam o contedo bsico de Sistemas de Potncia 1,

Finalmente, estudos de estabilidade tm o objetivo de determinar os limites operacionais

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

4. O SISTEMA POR UNIDADE

Vreal Definio de um valor

Um valor expresso em pu igual a um centsimo do mesmo valor, quando expresso de

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

V % 100 V pu 100 1,0145 101,45% .

As principais vantagens do sistema por unidade so:

1) Os fabricantes de equipamentos tais como geradores, motores e transformadores costu-

2) Equipamentos semelhantes (mesma tenso, mesma potncia, etc.) tm impedncias seme-

3) O uso do fator 3 minimizado nos clculos trifsicos em pu.

4) Como veremos, a impedncia de transformadores, quando expressa em pu, independen-

5) Em pu mais fcil identificar quando os valores de grandezas como tenses e potncias

6) Caso a tenso seja 1 pu, a potncia aparente e a corrente em pu sero numericamente

iguais, por causa do cancelamento do fator 3 , como segue

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

Em princpio, h um grande grau de arbitrariedade na escolha do valor base para determi-

1) Tenso eltrica V (em kV).

2) As bases de tenso, corrente e impedncia transformam-se de acordo com as relaes de

3) Linhas de transmisso e impedncias em srie e em paralelo no afetam as bases de ten-

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

As impedncias-base podem ser obtidas a partir da potncia-base e das tenses-base

Z b3 Zb2 Zb3 952,2

As reatncias do gerador e do transformador podem ser facilmente convertidas para pu

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

A reatncia em pu da linha de transmisso pode ser obtida dividindo-se a reatncia em

4.3. Mudana de base

Sbv = potncia-base velha (equipamento).

Sbn = potncia-base nova (sistema).

Vbv = tenso-base velha (equipamento).

Vbn = tenso-base nova (sistema).

Z = impedncia original do equipamento, em ohms.

Z vpu = impedncia em pu na base velha.

Z npu = impedncia em pu na base nova.

Retomando a definio de pu, podemos agora escrever

Z Zvpu Zbv Z npu Zbn .

Substituindo Z bv Vbv / Sbv e Z bn Vbn / Sbn , teremos

Prof. Alvaro Augusto W. de Almeida UTFPR A Work in Progress

ser obtida a partir de Vb1 , ou seja