Clculos de Curto-Circuito Trifsico

Alexandre de Oliveira Goulart e Gabriela Torllone de Carvalho Ferreira

Engenharia Eltrica 1o. Perodo
Professor: Rodney Josu Biezuner
Disciplina: Geometria Analtica e lgebra Linear

1. Introduo
Com o desenvolvimento da eletrnica, os circuitos eltricos, tal como as aplicaes para os mesmos, comearam a
ficar cada vez maiores e, conseqentemente, mais complexos. Desta forma, todos problemas relacionados sistemas
eltricos, em meados da dcada de 50, eram resolvidos mo ou com o uso de um analisador de redes. Este instrumento
era um dispositivo para fazer analogias eltricas, representando o circuito eltrico mas em uma escala bastante reduzida.
Uma opo para a investigao de circuitos eltricos seria atravs de mtodos matriciais que, embora existissem,
eram inviveis em virtude da falta de recursos computacionais realmente eficientes. Isso restringiu as investigaes de
circuitos a pequenos sistemas de equaes, envolvendo somente matrizes muito pequenas. Porm, com os avanos nas
ferramentas computacionais (entres elas est o computador digital), um mtodo que at ento se restringia apenas a um
instrumento de pesquisa pode se revelar eficiente para a anlise e soluo de sistemas eltricos de maior porte: os mtodos
matriciais.

2. O problema e seu modelamento

Quando o computador comeou a ser usado em grande escala, os clculos de fluxo de carga que antes eram feitos
mo puderam ser realizados mais facilmente, porm os primeiros programas de computador escritos para anlise de
circuitos apenas automatizaram estes mtodos manuais.
Isso subtilizava a capacidade do computador. O primeiro programa que realmente teve sucesso foi o desenvolvido
por Ward e Hale, que criaram um algoritmo que descreve o circuito eltrico atravs de um mtodo interativo de Newton
modificado. Aps isso, os programas que surgiram implementaram o algoritmo de Gauss-Seidel.
Todos esses algoritmos utilizam os dados do circuito sob a forma matricial, onde cada matriz guarda caractersticas
como impedncia, capacitncia e tenso.
Tudo isso foi desenvolvido para a anlise de problemas de fluxo de carga. Alm desse problema, circuitos sob
condies de curto circuito seguiram um desenvolvimento semelhante.

3. Clculos de Curto-Circuito Trifsico

Pode-se considerar um circuito eltrico, em condies de curto-circuito, como um sistema alimentado por vrias
fontes (geradoras) com uma nica carga ligada ao sistema no n sujeito ao curto-circuito. As correntes normais de carga dos
consumidores so usualmente ignoradas, pois so pequenas em comparao corrente de curto-circuito. Esta simplificao
no afeta a preciso do estudo de curto-circuito. Isto equivalente anlise estrutural de uma ponte apoiada em vrios
pilares e sujeita a uma s carga concentrada, desprezando-se os pesos dos vrios elementos da prpria estrutura. O restante
deste texto se concentra na anlise dos problemas de circuitos, mas deve-se lembrar que as tcnicas aqui desenvolvidas
aplicam-se igualmente s estruturas.
A anlise completa de curto-circuito de um sistema pode ser feita por simples operaes aritmticas, assim que as
tenses dos ns tenham sido determinadas para uma condio de falta particular. Coombe e Lewis determinaram as tenses
no n em estudos de curto-circuito pelo procedimento interativo de Gauss-Seidel desenvolvido por Ward e Hale para a
soluo dos problemas normais de fluxo de carga [1,2]. O mtodo pouco conveniente para estudos de curto-circuito, pois
cada condio de falta exige uma soluo interativa. Uma anlise completa de um sistema com mil ns (barras) pode exigir
cerca de vinte e trinta mil condies de falta, cada uma das quais teria que ser resolvida por um procedimento interativo
completamente diferente.
1

A anlise de curto-circuito de sistemas eltricos muito grandes feita da maneira mais eficiente, usando-se a matriz
Z-barra [3].

3.1 Descrio da matriz Z-barra

A matriz Z-barra contm as impedncias no ponto de cada n com relao a um n de referncia escolhido
arbitrariamente. A impedncia no ponto de um n a impedncia equivalente entre ele e a referncia. A matriz Z-barra
contm tambm a impedncia de transferncia entre cada barra do sistema e cada outra barra, com relao ao n de
referncia. As impedncias de transferncia so determinadas calculando-se as tenses que existiriam em cada uma das
outras barras do sistema, com relao referncia, quando uma barra em particular recebe uma injeo de corrente de uma
unidade (veja a Fig. 3.1).
A equao matricial que relaciona a matriz Z-barra s correntes injetadas nos ns e s tenses dos ns, a equao
3.1 a seguir:
ZI = E
Reconheceu-se logo que, se a matriz Z-barra usando como barra de referncia escolhida a barra comum atrs das
reatncias transitrias do gerador estivesse disponvel, a anlise completa do circuito poderia se facilmente feita com pouca
computao adicional. Lembrado-se que um circuito em condies de falta poderia ser considerado como tendo somente
uma corrente nodal, a equao matricial pode ser escrita como na equao 3.2 abaixo:

Z11 Z12 Z13 ... Z1k ... Z1n

Z21 Z22 Z23 ... Z2k ... Z2n
...................................................
...................................................
Zk1 Zk2 Zk3 ... Zkk ... Zkn
...................................................
Zn1 Zn2 Zn3 ... Znk ... Znn

0
0
..
..
Ik
..
0

E1
E2
..
..
Ek
..
En

onde o circuito est sujeito a uma nica injeo de corrente Ik aplicada no n k, que a barra sujeita condio de falta.
Obviamente, a coluna k permite a determinao do perfil de tenses do circuito quando ocorre um curto-circuito no n k,
desde que Ik tenha sido, ou possa ser determinada. O elemento Zkk a impedncia no ponto da barra k. Os elementos fora
da diagonal Zik so impedncias de transferncia entre as outras barras e a barra k.

Fig 3.1. Impedncia no ponto: Ek - r = IZkk ; Ik = 1,0; Zkk = Ek - r.

Impedncia de transferncia: Ik=1,0; Eir = IZik; Zik = Eir
2

Nos clculos de curto-circuito usual admitir que todos os geradores ligados ao circuito esto operando com uma
tenso de 1,0 p.u. atrs de suas reatncias internas. Este ponto comum atrs das reatncias dos geradores usado como
referncia. Portanto, pode-se considerar o circuito como sendo alimentado por uma nica fonte comum (veja a Fig. 3.2).
Fig. 3.2 Representao simplificada do sistema para estudo de curto-circuito.

O uso de 1,0 p.u. de tenso atrs da reatncia interna da mquina pode ser justificado pela aplicao do teorema de
Helmholtz-Thvenin. A tenso de circuito-aberto antes da falta no ponto de falta de aproximadamente 1,0 p.u. e pode,
portanto, ser tomada como sendo este valor. A impedncia de curto-circuito do sistema a partir do ponto de falta
determinada pela impedncia dos elementos do circuito (incluindo as impedncias internas das mquinas), com as fontes,
que so nulas, e a corrente devido a uma fonte de tenso superposta que reduz a tenso no ponto de falta para zero. Esta
tenso , obviamente, igual ao negativo da tenso antes da falta, sendo portanto, igual a 1,0 p.u. da terra barra comum
atrs das reatncias transitrias.
Quando qualquer n curto-circuitado, fica ligado terra. A tenso total , portanto, aplicada entre o n de
referncia e o n sujeito condio de falta. Por exemplo, para uma falta no n 6, o diagrama pode ser desenhado como
mostra a Fig. 3.3.
Como os elementos da matriz Z da equao 3.2 so as impedncias no ponto (elementos da diagonal) e as
impedncias de transferncia (elementos fora da diagonal), com relao barra de referncia as tenses da Eq. 3.2 so todas
medidas em relao ao n de referncia atrs das reatncias transitrias do gerador. O n de referncia est, portanto, a um
potencial zero em relao a si prprio, mas est a tenso total em relao terra no sistema real.

Fig. 3.3 O n 6 est em condio de curto-circuito.

A tenso obtida da Eq. 3.2 para a barra sob condio de falta a tenso total dos geradores em relao referncias.
No sistema real, a barra curto-circuitada est com um potencial nulo em relao terra. Esta diferena de tenses, depende
do ponto de referncia tomado, no deve causar nenhuma dificuldade, mas deve ser levada em considerao ao se expressar
os resultados de um clculo:
3

E pg = 1 E m
p
onde Ep a tenso da barra p com a terra como referencia, como seria medida no sistema real; Ep a tenso obtida do
clculo matricial da Eq. 3.2 e medida em relao barra de referncia atrs das reatncias transitrias do gerador.
Quando a barra k est sob condies de curto-circuito, a restrio de que a tenso total seja aplicada barra k pode
ser satisfeita injetado-se a corrente Ik, que determinada pela Eq. 3.3:

Ik =

1
Z kk

A corrente total de falta para qualquer barra , portanto, obtida tomando-se o recproco do elemento diagonal
correspondente da matriz Z. As tenses que aparecem nas outras barras do sistema, quando a barra k est em condies de
falta, dependem das impedncias de transferncia dadas pelos elementos fora da diagonal da coluna k da matriz Z-barra. Por
exemplo, a tenso com relao referncia na barra p para um curto-circuito na barra k seria dada pela Eq. 3.4.

E p = Z pk

1
Z kk

A corrente que vai da barra p barra q pela linha p-q, cuja impedncia Zlinha pq, dada pela Eq. 3.5:

I pq =

Z qk Z pk 1
Z linhap q Z kk

A corrente de falta total para uma falta em qualquer barra obtida pela Eq. 3.3 e a corrente em qualquer linha para um
curto-circuito em uma barra particular obtida pela Eq. 3.5. A anlise completa do sistema obtida por estas operaes
aritmticas simples, logo que a matriz Z tenha sido obtida.

3.2 O algoritmo de construo da matriz Z

Calcular a matriz de impedncias no ponto e de transferncia para um sistema eltrico de transmisso completo seria
totalmente impossvel. Entretanto, possvel, de maneira bastante simples, modificar a matriz Z de um sistema para a
adio de uma nica linha. Desta maneira, o sistema pode ser construdo comeando-se com uma nica linha de
transmisso, adicionando-se uma linha de cada vez, modificando-se a matriz para cada linha adicionada, e construindo o
sistema desejado e a matriz que corresponde ao sistema [3].

Preparao de dados
Desenha-se um diagrama do sistema. Os pontos de juno, onde duas ou mais impedncias de linhas de transmisso
de transformadores ou de geradores so ligadas, recebem um nico nmero de barra (n). O nmero zero reservado para a
barra de referncia. Em estudos de curto-circuito, a barra de referncia escolhida como sendo o ponto comum atrs de
todas as reatncias dos geradores.
Os dados so preparados descrevendo-se cada elemento do sistema de transmisso pelas duas barras nos terminais da
linha e sua impedncia em uma base p.u. comum. Estes dados do sequenciados por um algoritmo, partindo de uma ordem
aleatria para uma seqncias, de maneira que, medida que cada linha selecionada da lista de dados para processamento,
possvel lig-lo ao sistema previamente construdo. A primeira linha na lista deve ser uma linha que liga a referncia a
alguma barra do sistema, de maneira a fornecer uma trajetria at a referncia para a corrente injetada em qualquer n do
sistema que est sendo construdo.
Cada linha selecionada da lista deve pertencer a um dos trs tipos:
123-

Uma linha da referncia para uma nova barra.

Uma linha radial de uma barra existente para uma nova barra.
Uma linha entre duas barras j includas no sistema (linha de fechamento de lao).

So necessrias trs diferentes rotinas para modificar a matriz para a adio de uma linha ao sistema, dependendo do
tipo de linha a se adicionada.

Uma linha de referncia para uma nova barra

Uma linha da referncia para uma nova barra do sistema identificada ao se verificar que uma das barras a barra de
referncia e a outra barra no est includa no sistema j construdo.
A corrente injetada na nova barra k, que est ligada por uma linha radial referncia, no produzir tenso nas outras
barras no sistema ( veja a Fig. 3.4).
Uma injeo de corrente em qualquer barra do sistema que j tenha sido construdo no produzir nenhuma tenso
nova na barra k . Todos os elementos fora da diagonal da nova fila e da nova coluna so, portanto, nulos.

Z ik = Z ki = 0
(3.6)

ik

A impedncia no ponto da nova barra a impedncia da nova linha que est sendo adicionada. O elemento da
diagonal no novo eixo da matriz correspondente barra k dado pela Eq. 3.7.

Z kk = Z linhao k
Para a adio de uma nova linha radial da referncia a uma nova barra, aumente a matriz de uma fila e uma coluna de zeros.
O elemento na diagonal deste novo eixo a impedncia da nova linha que est sendo adicionada. O nmero de barra k
adicionado lista de barras que compe o sistema.

Fig 3.4. Adio de uma nova linha da referncia nova barra.

Adio de uma nova linha radial a uma nova barra

Uma linha radial ligando uma barra do sistema a uma nova barra identificada ao observar que nenhuma das barras
que descrevem a linha a referncia, e somente uma barra que descreve a linha est includa no sistema. Uma injeo de
uma unidade de corrente na barra q produz tenses em todas as outras barras do sistema idnticas s tenses que seriam
produzidas se a corrente fosse injetada na barra p (veja a Fig. 3.5).

Z qk = Z pk
kq
Z kq = Z kp
(3.8)
A impedncia no ponto da barra q igual a impedncia no ponto da barra p mais a impedncia da linha que est
sendo adicionada (veja a Fig. 3.5).
5

Z qq = Z pp + Z linhapq
(3.9)
Um novo eixo adicionado matriz correspondente nova barra q. Os elementos fora da diagonal da nova fila e
coluna so iguais aos elementos da fila e coluna da barra p do sistema existente. O elemento da diagonal obtido da Eq. 3.9.
A barra q adicionada lista de barras do sistema.

Fig. 3.5. Adio de uma linha radial de uma barra do sistema a uma nova barra.

Adio de uma linha de fechamento de lao

Uma linha de fechamento de lao identificada pelo fato de que ambas as barras que descrevem a linha esto
includas na lista de barras do sistema construdo. A adio de uma linha de fechamento de lao no produz um novo n no
sistema, que forneceria um novo eixo para a matriz. H, entretanto, uma nova descrio. A resposta do sistema introduo
de uma unidade de corrente no lao criado pela adio desta nova linha uma restrio que deve ser satisfeita.
A injeo de uma unidade de corrente na barra p faz com que apaream tenses em cada barra dos sistemas idnticos
aos elementos da coluna p da matriz Z (veja a Fig. 3.6). A injeo de uma corrente de 1,0 na barra q produz tenses nas
barras do sistema iguais aos elementos da coluna q da matriz Z, mas de sinal oposto. Uma corrente de lao igual unidade
pode ser considerada como um conjunto de correntes de Ip= 1,0 e Iq = 1,0 agindo simultaneamente. As tenses que
aparecem nas barras do sistema so iguais diferena das colunas correspondentes s barras p e q, como dado na Eq. 3.10.

Z11 Z12 ... Z1p ... Z1q ... Z1n

................... Z2p ... Z2q ... Z2n
................... Z3p ... Z3q ............
........................................................
................... Znp ... Znq ... Znn

0
0
Ip = 1,0
Iq = -1,0
0

Z1p Z1q
Z2p Z2q
....
....
Znp - Znq

(3.10)

Fig. 3.6. Simulao de uma linha de fechamento de lao por duas correntes injetadas.
A tenso que deve ser introduzida no lao para causar a circulao de uma unidade de corrente no lao criado pela
adio de uma nova linha pode ser calculada por 3.11 como fica evidenciado na Fig. 3.7.

E lao = I lao ( Zpp Zpq + Zqq Zpq + Z linhapq )

(3.11)
A impedncia no ponto do lao, Z laolao determinada fazendo-se I lao = 1,0 em 3.11.

Z laolao = Z pp + Z qq 2 Z pq + Z linhapq
(3.12)

Fig. 3.7. Determinao da impedncia no ponto do lao.

Um eixo de lao adicionado matriz Z da Eq. 3.10 na qual

Z i lao = Z ip Z iq

i lao

(3.13)
Z laoi = Z pi Z qi

i lao

(3.14)

O elemento da diagonal obtido de 3.12. o eixo do lao eliminado da matriz por uma reduo de Kron [4], pela
aplicao da Eq. 3.15. Todos os elementos que no pertencem fila ou coluna do lao so modificados. A fila e a coluna
do lao so eliminadas. A lista das barras do sistema permanece inalterada.
A reduo de Kron muito simples se um nico eixo estiver sendo eliminado. Isto ser mostrado posteriormente na
ilustrao do algoritmo de formao.

Reduo de Kron de uma matriz

Considere a matriz da Fig. 3.8. a matriz pode ser modificada para refletir as alteraes no

Z1
Z3

Z2
Z4

Fig. 3.8. Matriz particionada para ser reduzida.

sistema quando os eixos correspondentes s filas de Z 3 e s colunas de Z 2 so eliminados. Esta matriz modificada Z 1
pode ser considerada como uma reduo de circuito de Kron ou como uma eliminao algbrica. A reduo de Kron dada
por

Z 1 `= Z 1 Z 2 Z 41 Z 3
(3.15)
A validade da reduo 3.15 pode ser provada considerando-se a equao matricial 3.16.

A1 A2
A3 A4

X1
X2

B1
B2
(3.16)

na qual A1, A2, A3 e A4 podem ser consideradas matrizes ou coeficientes simples e X1, X2, B1 e B2 so vetores ou
variveis simples, respectivamente.
A Eq. 3.16 em forma expandida

A1 X 1 + A2 X 2 = B1
(3.17)
A3 X 1 + A4 X 2 = B2
(3.18)
Reescrevendo 3.18 vem

A4 X 2 = B2 A3 X 1
(3.19)
1

A pr-multiplicao por A4 d
8

X 2 = A41 B2 A41 A3 X 1
(3.20)
A substituio de X 2 em 3.17 d

A1 X 1 + A2 ( A41 B2 A41 A3 X 1 ) = B1
(3.21)
Rearranjando e reunindo os termos vem

( A1 A2 A41 A3 ) X 1 = B1 A2 A41 B2
(3.22)
V-se que o coeficiente da incgnita X 1 o resultado de uma reduo de Kron em A1 , como dado por 3.15, quando
a incgnita X 2 eliminada.

4. Um exemplo ilustrativo do algoritmo de construo

Considere o circuito mostrado na Fig. 3.9 e os seguintes dados:

Fig. 3.9.
A matriz do circuito construda montando-se o sistema linha por linha e modificando-se a matriz para refletir a
alterao nas impedncias equivalentes do circuito com a adio da linha. A lista de linhas j foi reordenada a partir de uma
lista aleatria para uma seqncia tal que seja possvel ligar cada linha ao sistema, quando ela for selecionada da lista para
processamento.
Adio da primeira linha. A primeira linha deve sempre ser uma linha ligada referncia. Assim, a linha 0-1 pode
ser a primeira linha a ser processada. Neste ponto, nenhum circuito, nem matriz e nem elementos na lista de barras do
sistema para determinar o tipo de linha e o algoritmo a ser usado. V-se que a linha est ligada referncia e a barra 1. A
barra 1 comparada com a lista de barras do sistema. Nesta altura, ainda no h nenhuma barra no sistema. A linha ,
portanto, um caso degenerado da adio de uma linha a partir da referncia at uma nova barra (veja a Fig. 3.4).
impossvel adicionar uma fila e uma coluna de zeros matriz, pois ainda no h nenhuma matriz neste ponto. O
elemento da diagonal do novo eixo a impedncia da linha que est sendo adicionada. A nova barra adicionada `a lista de
barras. Aps adicionada a primeira linha, temos:

Matriz

1
1 [0,01]

lista de barras 1

O diagrama do sistema mostrado a seguir.

A matriz mostra que uma unidade de corrente injetada na barra 1 faria aparecer nesta uma tenso de 0,01 medida em
relao barra de referncia.
Adio da segunda linha. A prxima linha (0-2) selecionada da lista de dados para processamento. O exame dos
nmeros de barra 0 e 2 e a comparao com a lista de barras do sistema mostram que esta linha vai da barra de referncia
(0) a uma nova barra 2. tambm um linha do primeiro tipo. Aumente a matriz com uma fila e uma coluna de zeros. O
elemento da diagonal do novo eixo a impedncia da nova linha. Adicione a barra 2 lista de barras.

Matriz

1
2

0,01
0

0
0,015

Lista de barras 1,2

O diagrama do sistema mostrado a seguir.

10

A matriz mostra-nos que a injeo de uma unidade de corrente na barra 1, saindo pela referncia, faz aparecer uma
tenso nula na barra 2. a injeo de uma unidade de corrente na barra 2 produzir uma tenso de 0,015 na barra 2 e uma
tenso nula na barra 1.
Adio da terceira linha. A prxima linha (1-2) selecionada para processamento. O exame dos nmeros da barra
mostra que esta no uma linha de referncia. A comparao dos nmeros das barras da linha com os nmeros das barras
do sistema, mostra que esta uma linha do terceiro tipo (fechamento de lao).
A matriz aumentada de uma fila e uma coluna de lao tomando-se as diferenas das filas (e colunas)
correspondentes s barras 1 e 2, Eq. 3.13. o elemento da diagonal obtido pela Eq. 3.12.

1
Matriz

1
2
lao

0,01
0
0,01

lao

0
0,01
0,015 -0,015
-0,015 0,109

Lista de barras 1,2

O diagrama do sistema mostrado abaixo

Fig. 3.10
A matriz reduzida pela aplicao da reduo de Kron usando a Eq. 3.15 como j salientamos, a reduo de Kron
1

grandemente simplificada quando nela est envolvido somente um eixo. Neste caso Z 4 torna-se 1/Z e no necessria a
inverso de uma matriz. A modificao dos elementos no pertencentes fila e coluna do lao pode ser executada mais
facilmente elemento por elemento que pela aplicao da Eq. 3.15 como uma equao matricial.
Pode-se verificar facilmente que a modificao de um elemento, Z ij , simplesmente:
Zij =Zij Zi-lao

___1___
Zlao-lao

Zlao-j

(3.23)

11

lao

Zij

Zi-lao

lao

Zlao-j

Zlao-lao

Fig. 3.11. elementos matriciais usados na modificao do elemento Zij pela reduo de Kron
A aplicao da Eq. 3.23 usada para modificar todos os elementos no pertencentes ao eixo do lao na matriz. O
eixo do lao ento eliminado.
Z11 = Z11 Z1-lao

___1___
Zlao-lao

Zlao-1 =

= 0,01 (0,01)(0,01)
0,109
= 0,01 0,00091743 = 0,00908257
Z12 = Z12 Z1-lao

___1___
Zlao-lao

Zlao-2 =

= 0 (0,01)(-0,015)
0,109
= 0 + 0,00137615 = 0,00137615
Z22 = 0,015 (-0,015)(-0,015)
0,109
= 0,015 0,00206421
=0,0129359

12

A matriz modificada

1
Matriz

1
2

0,00908257
0,00137615

0,00137615
0,01293579

Lista de barras 1,2

A Reduo da Matriz Uma Converso Tringulo-Estrela. A reduo da matriz pode ser vista como uma
reduo de circuito tipo tringulo-estrela. O circuito da Fig. 3.10 pode ser convertido na estrela equivalente pela converso
tringulo-estrela normal.

Z1 = __(0,01)(0,015)___ = (0,01)(0,015) = 0,00137615

0,01+0,015+0,084
0,109
Z2 = (0,01)(0,084) = 0,00770642
0,109
Z3 = (0,015)(0,084) = 0,01155964
0,109
Se o circuito da Fig. 3.12 excitado por uma unidade de corrente na barra 1, a impedncia no ponto Z 11 a soma de

Z1 e Z 2 .
Z11 = 0,00137615 + 0,0077062 = 0,00908257
referncia
Z1
Z2

Z3
2

Fig. 3.12.
A tenso da barra 2 igual `a queda de tenso na impedncia Z 1 = 0,00137615 . Isto a impedncia de
transferncia Z 12 .
A impedncia no ponto da barra 2 ento obtida injetando uma unidade de corrente na barra 2.

Z 22 = Z 1 + Z 3 = 0,00137615 + 0,01155964 = 0,01293579

Estes valores concordam com a matriz obtida pela reduo de Kron. Assim, a reduo pode ser encarada como uma
reduo de circuito ou como uma manipulao algbrica.
Adio da quarta linha. Continua-se o algoritmo de construo, selecionando-se a prxima linha 0-3 para
processamento. Esta linha identificada como uma linha do primeiro tipo. A matriz aumentada de uma fila e uma coluna
de zeros e a diagonal do novo eixo a impedncia da linha, 0,005. a nova barra adicionada lista de barras.

13

1
1
2
3

Matriz

0,00908257
0,00137615
0

0,00137615
0,01293579
0

0
0
0,005

Lista de barras 1,2,3

Adio da quinta linha. A prxima linha na lista de dados (2-3)_ uma linha de fechamento de lao, pois ambas as
barras esto na lista de barras do sistema. O eixo do lao a diferena das colunas correspondentes s barras 2 e3. o
elemento da diagonal obtido da Eq. 3.12.

Zlao-lao = Z22 + Z33 2Z23 + Zlinha 2-3

=0,01293578 + 0,005 (2)(0) + 0,122
=0,13993578
A matriz aumentada :

1
2
3
lao

lao

0,00908257
0,00137615
0
0,00137615

0,00137615
0,01293579
0
0,01293579

0
0
0,005
-0,005

0,00137615
0,01293579
-0,005
0,13993579

A reduo da matriz pela aplicao da Eq. 3.23 d

Matriz

1
2
3

0,00906904
0,00124893
0,00004917

0,00124893
0,01173999
0,00046220

0,00004917
0,00046220
0,00482135

O diagrama do sistema mostrado a seguir.

Adio da sexta linha. A linha 2-4 identificada como do segundo tipo, uma linha que vai da barra existente 2 at a
nova barra 4. um novo eixo adicionado matriz. O elemento da diagonal determinado pelo uso da Eq. 3.9.

Z 44 = Z 22 + Z linha 2 4 = 0,011739999 + 0,084 = 0,095739999

14

Os elementos fora da diagonal so obtidos da Eq. 3.8 na qual q=4 e p=2. a fila correspondente barra 4 idntica
fila da barra 2.

1
1
2
3
4

0,00906904
0,00124893
0,00004917
0,00124893

0,00124893
0,01173999
0,00046220
0,01173999

0,00004917
0,00046220
0,00482135
0,00046220

0,00137615
0,01173999
0,00046220
0,09573999

O diagrama do sistema dado a seguir.

Adio da stima linha. Adio da linha 3-5 d:

Z 55 = Z 33 + Z linha 35 = 0,00482135 + 0,037 = 0,04182135

A coluna 5 uma duplicata da coluna 3.

1
1
2
3
4
5

0,00906904
0,00124893
0,00004917
0,00124893
0,00004917

0,00124893
0,01173999
0,00046220
0,01173999
0,00046220

0,00004917
0,00046220
0,00482135
0,00046220
0,00482135

0,00137615
0,01173999
0,00046220
0,09573999
0,00046220

0,00004971
0,00046220
0,00482135
0,00046220
0,04182135
listas de Barras 1,2,3,4,5

Adio da oitava linha. A adio ao sistema da linha 1-6, que uma linha ligando uma barra existente 1 a uma nova
barra 6, executada como indicado para a linha 3-5 na etapa 7 (linha do tipo 2).
Adio da nona linha. A linha 6-7 tambm uma linha do tipo 2, ligando uma barra existente 6 a uma nova barra 7.
o processo ilustrado na etapa 7. aps a adio destas duas linhas a matriz

1
2
3
4
5
6
7

0,00906904
0,00124893
0,00004917
0,00124893
0,00004917
0,00906904
0,00906904

0,00124893
0,01173999
0,00046220
0,01173999
0,00046220
0,00124893
0,00124893

0,00004917
0,00046220
0,00482135
0,00046220
0,00482135
0,00004917
0,00004917

0,00137615
0,01173999
0,00046220
0,09573999
0,00046220
0,00124893
0,00124893
15

0,00004971
0,00046220
0,00482135
0,00046220
0,00482135
0,00004917
0,00004917

0,00906904
0,00124893
0,00004917
0,00124893
0,00004917
0,13506904
0,13506904

0,00906904
0,00124893
0,00004917
0,00124893
0,00004917
0,13506904
0,13506904

Adio da dcima linha. A linha 4-7 uma linha de fechamento de lao, pois ambas as barras j constam do
sistema. A fila e a coluna do lao so obtidas tomando-se a diferena das filas 4 e 7 e das colunas 4 e 7. o elemento da
diagonal deste novo eixo obtido da Eq. 3.12.

Z laolao = Z 44 + Z 77 2Z 47 + Z linha 47 =
= 0,09573999 + 0,30306904 2(0,00124893) + 0,084
O eixo do lao eliminado como ilustrado nas etapas 3 e 5.
Adio da dcima-primeira linha. A linha 5-8 uma linha do tipo 2 e sua adio feita da maneira ilustrada na
etapa 6.
Adio da ltima linha. A linha 7-8 uma linha de fechamento de lao. O mtodo j foi ilustrado na etapa 3. V-se
facilmente que, independente da complexidade do sistema, ela pode ser obtida por este mtodo simples de adicionar uma
linha de cada vez.
A matriz completa do sistema-exemplo dada para fins de referncia.

1
2
3
4
5
6
7
8

0,00906904
0,00124893
0,00004917
0,00124893
0,00004917
0,00906904
0,00906904
0,00102639

0,00124893
0,01173999
0,00046220
0,01173999
0,00046220
0,00124893
0,00124893
0,00220239

0,00004917
0,00046220
0,00482135
0,00046220
0,00482135
0,00004917
0,00004917
0,00375267

0,00137615
0,01173999
0,00046220
0,09573999
0,00046220
0,00124893
0,00124893
0,01464439

0,00004971
0,00046220
0,00482135
0,00046220
0,00482135
0,00004917
0,00004917
0,03199261

0,00906904
0,00124893
0,00004917
0,00124893
0,00004917
0,13506904
0,13506904
0,01219496

7
0,00906904
0,00124893
0,00004917
0,00124893
0,00004917
0,13506904
0,13506904
0,02708638

8
0,00102639
0,00220239
0,00375267
0,01464439
0,03199261
0,01219496
0,02708638
0,06023255
16

Anlise de faltas em um sistema

Uma anlise completa do sistema possvel to logo a matriz Z do sistema tenha sido completada. Usando os valores
que esto na matriz do sistema exemplo para fins de ilustrao, considere uma condio de falta no n 3. o elemento da
matriz Z 33 = 0,00475959 mostra que, se uma tenso daquele valor for aplicada entre a barra 3 e a referncia, uma
corrente total de 1,0 ir percorrer o circuito. A tenso total do gerador produzir uma corrente que pode ser determinada
considerando-se

I = E
I
E
(3.24)
Onde I=1,0 quando E=0,00475959. deseja-se determinar Iquando E=1,0.

I'=

1,0
= 210,10 p.u.
0,00475959

Este o resultado que obtido considerando-se a matriz 3.2 e a Eq.3.3.

O valor total da falta pode ser obtido tanto em ampres, p.u., ou MVA dividindo os ampres-base, a unidade ou base
MVA p.u. dos dados da linha pelos elementos correspondentes da diagonal da matriz (veja a Eq. 3.5).
A contribuio de uma linha falta calculada usando-se a Eq. 3.5. a contribuio da barra 2 falta

I 23 =

Z 33 Z 32 1,0
= (0,00475959 0,00055371) /(0,122)(0,00475959) = 7,2 p.u
Z linha 23 Z 33

Nota: a base aqui foi 1,0 e a corrente est em p.u. Poderia tambm ser usados ampres-base ou MVA.

Tenses nas barras durante uma condio de falta

A tenso na barra 2 quando injetada uma unidade de corrente na barra 3 Z 23 = 0,00055371 . A tenso na barra
3 neste instante 0,00475959. a diferena de tenso Z 33 Z 32 = 0,00475959 0,00055371 = 0,00420588 . Esta
diferena ocorreria se a corrente fosse 1,0 mas a corrente de falta 1

Z 33

conforme determinado. A diferena de tenso

entre as barras 2 e 3 , portanto,

Z 33 Z 32 0,00420588
=
= 0,8835 p.u.
Z 33
0,00475959
mas a tenso da barra 3 em condies de falta zero. Portanto, a tenso na barra 2 = 0,8835.

Abertura de uma linha durante um estudo

Uma linha de um sistema pode ser aberta ou removida adicionando-se uma linha em paralelo linha existente. A
impedncia da nova linha a ser adicionada o negativo da impedncia da linha original. So usadas as equaes de
fechamento de lao (Eqs. 3.12 e 3.13) e a eliminao do lao por uma reduo de Kron (Eq. 3.15).

17

No decorrer de um estudo completo de faltas freqentemente desejvel abrir todas as linhas ligadas barra em
condio de falta, uma de cada vez, e obter a nova falta total e a contribuio das linhas restantes. No desejvel modificar
a matriz do sistema total, pois uma grande quantidade de clculos desnecessrios seria feita em elementos que no so
necessrios anlise. Mais ainda, tambm indesejvel (por causa do arredondamento) remover uma linha, adiciona-la de
novo, remover outra, adiciona-la de novo, e assim por diante. Os erros se acumulariam nos elementos da matriz Z devido
modificao repetida da matriz.
O melhor mtodo consiste em extrair da matriz total uma pequena matriz que inclua as impedncias no ponto e de
transferncia da barra que estar em falta e de suas barra vizinha prximas.
Por exemplo, se a barra 3 deve sofrer uma falta, extrai-se a pequena matriz.

3
2
5

0,00475959
0,00055371
0,00425613

0,00055371
0,01134623
0,00137805

0,00425613
0,00137805
0,03662437

Para abrir a linha 3-5 adiciona-se uma linha ( de fechamento de lao) cuja impedncia o negativo da impedncia da
linha original 0,037. o elemento da diagonal do eixo do lao

Z laolao = Z 33 + Z 55 2 Z 35 + Z linha35 = 0,00475959 + 0,03662437 2(0,00425613) 0,037

= 0,00412830
A coluna do lao obtida subtraindo-se a coluna 5 da coluna 3.

3
3
2
5
lao

0,00475959
0,00055371
0,00425613
0,00050346

0,00055371
0,01134623
0,00137805
-0,00082434

lao

0,00425613
0,00137805
0,03662437
-0,03236824

0,00050346
-0,00082432
-0,03236824
-0,00412830

Para maior eficincia, a reduo de Kron usada para modificar somente a fila 3, pois uma falta na barra 3 pode ser
completamente analisada, usando-se somente esses valores.
O vetor fila modificado que indica a abertura da linha 3-5

3 [0,00482099 0,00045317 0,00030870]

(3.25)
A nova falta total

base
1
=
= 207,43 p.u.
Z 33
0,00482099
A contribuio da barra 2 de 7,4 p.u. a corrente da barra 5 para a barra 3 atravs da linha X=0,037 tem o mesmo
valor que a corrente na linha X=-0,037 mas com sinal contrrio. A contribuio lquida de 5 para 3 , portanto, zero.

18

O vetor fila, que d os valores para a falta na barra 3 com a linha 3-5 aberta, fornece uma excelente oportunidade
para a determinao do valor da corrente de falta que ocorreria se a falta fosse retirada da barra 3 e colocada no lado da
linha correspondente ao disjuntor na barra 5, na linha 3-5 (veja a Fig. 3.13).
O vetor fila d as impedncias no ponto e de transferncia da barra 3 com a linha aberta. A impedncia no ponto da
barra 9 obtida da Eq. 3.9. Como a linha 3-9 radial a partir da barra 3,

Z qq = Z 33 + Z linha 35 = 0,00482099 + 0,037 = 0,04182099

A falta total no ponto 9 = base/0,04182099 =- 23,91 p.u.

Fig. 3.13 Falta transferida para a extremidade aberta da linha.

Disjuntores de acoplamento de barras

Alguns sistemas eltricos tm disjuntores de acoplamento de barras que podem ser abertos ou fechados por motivos
de operao. Por exemplo, uma estao geradora pode Ter dois geradores ligados a um sistema de barras como mostrado
pela Fig. 3.14.
Se qualquer dos geradores estiver fora de servio para a manuteno, ou por alguma outra razo de operao, pode
ser necessrio fechar o disjuntor de acoplamento de barras. A matriz pode ser modificada para representar a condio do
disjuntos de acoplamento de barras fechado adicionando-se uma linha de impedncia nula entre as barras P e Q. Este um
mtodo perfeitamente satisfatrio de modelar o sistema para o disjuntor de acoplamento de barras fechado, mas seria
impossvel abrir novamente o disjuntor mais adiante no estudo. Isso se torna evidente considerando-se a adio de uma
linha de impedncia 0,0 numa tentativa de abrir a linha de impedncia 0,0 correspondente ao disjuntos fechado. Quando o
disjuntor fechado, as duas colunas (P e Q) se tornam idnticas em valor, pois existe uma impedncia nula ligada entre as
barras. O elemento da diagonal para a adio do elemento de impedncia 0,0 numa tentativa de abrir o disjuntor, torna-se

Zlao-lao = Zpp + Zqq 2Z pq + (-0,0)

mas

Zpp = Zqq =Zpq

19

quando as barras so ligadas e Zlao-lao =0.

A substituio deste valor na Eq. 3.23 no permitida. Portanto o disjuntor de acoplamento de barras no pode ser
aberto se esta representao for usada.

Disjuntor de acoplamento de barras que pode ser aberto

Para se ligar duas barras com uma impedncia de acoplamento de barras que possa ser aberta, introduz-se uma barra
fictcia entre as barras a serem ligadas e adiciona-se uma linha de P a T com uma impedncia Z1 e uma linha de T e Q com
uma impedncia de Z1 (veja a Fig. 3.15).
A impedncia da barra P para a barra Q nula, mas agora o disjuntor pode ser aberto adicionando-se uma linha da
barra T barra Q com impedncia +Z1 . Isto remove um lado do disjuntor de acoplamento de barras.
Z1

Barra P

Barra T

-Z1

Barra Q

Fig. 3.15. Simulao de um disjuntor de acoplamento de barras fechado.

4. Limitao do mtodo da matriz Z

Como a matriz Z cheia, ou seja, existe uma impedncia de transferncia para cada barra em relao a cada outra
barra, as exigncias de memria para seu armazenamento so muito grandes. Devido a relao X para R dos sistemas de alta
tenso ser elevada (de ordem de 10, por exemplo) a matriz Z pode ser calculada usando-se somente os valores de X sem
afetar seriamente a preciso do estudo. Mais ainda, a matriz simtrica e pode, portanto, ser armazenada na forma
triangular superior. Um sistema de 100 barras necessitar de uma matriz de 5.050 elementos N(N+1)/2. Um sistema de 200
barras tem 20.100 elementos. At cerca de 1958 os estudos de falta nos grandes sistemas eltricos eram feitos em modelos
analgicos com menos de 100 barras. A introduo dos programas de computadores digitais aumentou muito rapidamente o
nmero de barras usadas na representao. Tornou-se necessrio imaginar mtodos para aumentar o tamanho do programa.

5. Referncias Bibliogrficas
1- J. R. Ward e H. W. Hale, Digital solution of power flow problems, Trans. AIEE, Vol. 75.
2- R. A . Frazer, W. J. Duncan, e A R. Collar, Elementary Matrices, Cambridge Press, 1995.
3- Homer E. Brown, Grandes Sistemas Eltricos, LTC/EFEI.

20