Abstract— In this paper, it is presented a methodology for hree-phase distribution transformer modeling, considering several types of transformer configuration, to be used in algorithms of power flow in three-phase radial distribution networks. The paper provides a detailed discussion about the models and the results from an implementation of the power flow algorithm. The results, taken from three different networks, are presented for several transformer configurations and for voltage regulators as well. Keywords— Distribution transformer modeling, Nodal admittance matrices, Three-Phase power flow. I. NOMENCLATURA bus trafoY Matriz de admitância nodal do transformador; , pp ssY Y Submatrizes da matriz de admitância nodal bus trafoY do transformador para o lado primário e secundário respectivamente; , ps spY Y Submatrizes da matriz de admitância nodal bus trafoY do transformador para o lado primário- secundário e o lado secundário-primário respectivamente; , , YI YII YIII Submatrizes derivadas da bus trafoY ; ty Admitância de dispersão do transformador em p.u.; , , a b c p p pI I I Injeção de corrente trifásica no lado primário; , , a b c s s sI I I Injeção de corrente trifásica no lado secundário do transformador; , , a b c p p pV V V Tensão trifásica no lado primário do transformador; , , a b c s s sV V V Tensão trifásica no lado secundário do transformador; José. L. Choque é aluno de pós graduação na Unesp de Ilha Solteira, jluischoque@yahoo.es. Dario. E. Rodas é professor da Universidad tecnologica de Pereira – Colômbia, drodas@utp.edu.co. Antonio. Padilha-Feltrin é professor da UNESP campus de Ilha Solteira, padilha@dee.feis.unesp.br. M Conjunto de segmentos de linha conectados para barra j ; jn Segmentos de linha que pertencem a M ; , , a b c j j jI I I Injeção de correntes trifásicas da carga na barra j ; , , a b c jn jn jnS S S Potência trifásica na linha jn ; , , a b c j j jV V V Tensão trifásica na barra j ; 1 2 0, , p p pV V V Componentes de seqüências positiva, negativa e zero da tensão trifásica do lado primário respectivamente; , , a b c i i iY Y Y Admitância shunt na barra i ; , , a b c Fases no sistema trifásico. II. INTRODUÇÃO ários métodos de cálculo de fluxo de potência trifásico em redes radiais de distribuição surgiram nos últimos anos. Atualmente, os engenheiros de distribuidoras de energia têm disponíveis mais e melhores conjuntos de dados, de componentes da rede e de comportamento das cargas, para utilizarem em programas computacionais de cálculo de fluxo de potência. Como conseqüência, mais e melhores modelos de componentes são frequentemente necessários para serem implementados nos algoritmos de fluxo de potência e em outros. Este artigo apresenta e discute em detalhes uma implementação de modelos de transformadores trifásico de distribuição, para uso em programas computacionais baseados no método backward / forward sweep [1], [2], [3] e [4]. Particularmente em [2] é dada ênfase à modelagem de linhas considerando-se redes a quatro fios, agora o objetivo é a consideração dos transformadores. Dentre os métodos dedicados a redes radiais (ou com poucas malhas, que são solucionados em [3]) de distribuição, este método tem se destacado pela boa característica de convergência, pouco uso de memória computacional e simplicidade de compreensão e implementação. Neste tipo de método, e em muitos outros também, surgem problemas para representação correta dos vários tipos de conexões de transformadores de distribuição. Alguns trabalhos anteriores propuseram soluções para resolver o problema, dentre eles destacam-se as implementações [5] e [6]. Assim sugerem-se estas referências V Distribution Transformer Modeling for Application in Three-Phase Power Flow Algorithm J. L.Choque , D. Rodas, A.Padilha-Feltrin, Senior Member 196 IEEE LATIN AMERICA TRANSACTIONS, VOL. 7, NO. 2, JUNE 2009 para encontrar uma boa revisão sobre os temas fluxo de potência trifásico e de modelagem de transformadores de rede, e também sobre as dificuldades de modelagem devido a problemas como singularidade de matrizes. Cabe destacar que algumas publicações clássicas na área, dentre elas [7], [8] e [9], não detalham como implementar modelos diferentes do Yg-Yg em programas de fluxo de potência trifásico. Neste trabalho segue-se a linha de modelagem discutida em [5] e [6], procurando-se detalhar e discutir as implementações de forma a dar uma contribuição adicional na implementação de modelos de transformadores trifásicos. III. MODELAGEM DE TRANSFORMADORES DE DISTRIBUIÇÃO Os transformadores trifásicos podem ser representados partindo-se da matriz admitância primitiva. Considerando esta primitiva pode-se obter a matriz admitância nodal para as diversas conexões, conforme descrito em [7], [8] e [9]. Por facilidade de representação das diferentes conexões, pode-se colocar a matriz admitância nodal na forma de submatrizes (3x3) conforme a equação (1). pp ps bus trafo sp ss Y Y Y Y Y ⎡ ⎤ = ⎢ ⎥ ⎣ ⎦ (1) O arranjo das submatrizes de busY do transformador permite representar as diferentes conexões dos transformadores de distribuição. Nas Tabelas 1 e 2 são mostradas as submatrizes correspondentes aos tipos mais comuns de conexões, tanto para transformadores abaixadores (Tabela 1) quanto transformadores elevadores (Tabela 2). Estas tabelas podem ser comparadas com as apresentadas em [1], [6], [10], [11] e [12]. Para ilustrar como são obtidas estas tabelas, no Apêndice A mostra-se o caso da obtenção da bus trafoY para as conexões Yg -Yg e D - Yg para transformador elevador, incluindo taps de ambos os lados. TABELA I SUBMATRIZES PARA CONEXÕES DO TRANSFORMADOR ABAIXADOR Primário Secundário Ypp Yss Yps Ysp Yg Yg YI YI -YI -YI Yg Y YII YII -YII -YII Yg Δ YI YII YIII YIIIT Y Yg YII YII -YII -YII Y Y YII YII -YII -YII Y Δ YII YII YIII YIIIT Δ Yg YII YI YIII YIIIT Δ Y YII YII YIII YIIIT Δ Δ YII YII -YII -YII Admitância PrópriaConexão Admitância Mútua TABELA II SUBMATRIZES PARA CONEXÕES DO TRANSFORMADOR ELEVADOR Primário Secundário Ypp Yss Yps Ysp Yg Yg YI YI -YI -YI Yg Y YII YII -YII -YII Yg Δ YI YII YIIIT YIII Y Yg YII YII -YII -YII Y Y YII YII -YII -YII Y Δ YII YII YIIIT YIII Δ Yg YII YI YIIIT YIII Δ Y YII YII YIIIT YIII Δ Δ YII YII -YII -YII Admitância MútuaConexão Admitância Própria Em que: t 1 0 0 YI 0 1 0 .y 0 0 1 ⎡ ⎤ ⎢ ⎥= ⎢ ⎥ ⎢ ⎥⎣ ⎦ (2) t 2 1 1 1YII 1 2 1 .y 3 1 1 2 − −⎡ ⎤ ⎢ ⎥= − −⎢ ⎥ ⎢ ⎥− −⎣ ⎦ (3) t 1 1 0 1YIII 0 1 1 .y 3 1 0 1 −⎡ ⎤ ⎢ ⎥= −⎢ ⎥ ⎢ ⎥−⎣ ⎦ (4) Podem-se relacionar correntes e tensões primárias e secundárias do transformador (para fins de cálculo de fluxo de potência) com uso da matriz da equação (1): ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ ⎥ ⎦ ⎤ ⎢ ⎣ ⎡ =⎥ ⎦ ⎤ ⎢ ⎣ ⎡ s p sssp pspp s p V V . yy yy I I (5) IV. REPRESENTAÇÃO DOS TRANSFORMADORES NO ALGORITMO DE FLUXO DE POTÊNCIA BACKWARD / FORWARD SWEEP As tradicionais implementações do algoritmo Backward / Forward [3], [5] e [13] para o fluxo de potência necessitam basicamente, obter: Etapa Backward - Cálculo das injeções de correntes em todas as barras devido a cargas diretamente conectadas, assim como, cálculo dos fluxos de correntes em todas as linhas, usando-se os valores de tensão estimados (primeira iteração) ou calculados (demais iterações). Nesta etapa os cálculos são realizados partindo-se dos extremos do alimentador e seguindo-se para a subestação. Etapa Forward - Cálculo das tensões em todos os nós, utilizando-se os valores de fluxo de correntes obtidos na etapa backward. Nesta etapa parte-se da subestação (nó fonte) e segue-se para os nós extremos. Na Fig. 1 destaca-se uma seção de linha de uma rede radial, na qual foi incluído o transformador de distribuição. A partir CHOQUE et al.: DISTRIBUTION TRANSFORMER MODELING 197 desta figura serão detalhados os cálculos necessários nas etapas Backward e Forward (entre duas barras conectadas por um transformador). a jlS b jlS c jlS a b c nó jnó i jViV pp ps sp ss Y Y Y Y ⎡ ⎤ ⎢ ⎥ ⎣ ⎦ a b c nó fonte nó l Transformador a sI b sI c sI a pI b pI c pI nó ka jkS b jkS c jkS a b c Carga b kI c kIa kI a iY b iY c iY Figura 1. Seção de linha trifásica com o transformador embutido. A. Etapa Backward 1) Injeção de Corrente na Barra Terminal do Secundário Se o secundário do transformador está conectado diretamente na barra “ j ”, como mostrado na Fig. 1, pode-se obter a corrente no secundário a partir da diferença entre as correntes que entram e que saem do nó “ j ” ( jI representa a contribuição de cargas, se houver, conectadas diretamente na barra). Assim, tem-se: ( ) ( ) ( ) * * * / / / a a jn j b b s j jn j jn M c c jn j S V I I S V S V ∈ ⎡ ⎤ ⎢ ⎥ ⎢ ⎥= − ⎢ ⎥ ⎢ ⎥ ⎢ ⎥⎣ ⎦ ∑ (6) 2) Injeção de Corrente na Barra Terminal do Primário No procedimento backward sweep sV e sI são conhecidos, enquanto que pV e pI são os cálculados. A corrente no terminal primário não pode ser obtida de forma similar àquela que foi utilizada para a corrente no secundário, porém deve-se usar a equação (5), ou mais especificamente: p pp p ps sI Y V Y V= + (7) Mas a tensão do lado primário pV a ser utilizada não pode ser estimada, e sim calculada de forma a levar para o lado primário o efeito da injeção de corrente existente no secundário. Assim, embora tradicionalmente não existam cálculos de tensão na etapa Backward, agora eles são necessários para ser possível considerar corretamente os modelos dos transformadores. Portanto, conhece-se a tensão sV (estimada para a primeira iteração ou calculada para as seguintes iterações na etapa Forward) e a corrente do secundário sI (calculada no passo anterior), enquanto que a tensão pV e a corrente pI do primário devem ser calculadas. A corrente do primário pI pode ser obtida a partir da equação (7), desde que antes seja calculada a tensão pV . Esta tensão deve ser obtida a partir da equação (5), ou seja: ( )1 p sp s ss sV Y I Y V−= − (8) É importante notar que a submatriz spY é singular para todas as conexões dos transformadores, exceto para Yg-Yg. Isto é observado melhor através das equações (2), (3) e (4), onde se pode ver que ambas submatrizes YII e YIII são singulares, e YI é não singular. Em essência, a singularidade das matrizes dessas configurações de transformadores aparece devido a falta de um ponto de referência de tensão para ambos os lados do transformador. Para poder representar todos os tipos de conexões de transformadores no algoritmo Backward / Forward é preciso evitar o problema da singularidade da spY na etapa Backward. a) Resolvendo o problema de Singularidade na Etapa Backward A fim de resolver o problema de tal singularidade, faz-se a transformação da tensão do primário pV em suas componentes simétricas. ( )1 2 0 p p pV V V+= + (9) Em que ( )1 2 pV + representa as componentes de seqüência positiva e negativa, e 0 pV a componente de seqüência zero. Assim pode-se re-escrever a equação (8) como: ( )( )1 2 0 sp p p s ss sY V V I Y V+ + = − (10) O produto de spY por 0 pV é sempre zero para todas as configurações nas quais spY é representada por YII ou YIII (ou ainda tYIII ), portanto: 0. 0pYII V = (11) 0. 0pYIII V = e 0 . 0t pYIII V = (12) Então, a equação (10), exceto para a conexão Yg – Yg, pode ser reduzida a: ( )1 2 sp p s ss sY V I Y V+ = − (13) 198 IEEE LATIN AMERICA TRANSACTIONS, VOL. 7, NO. 2, JUNE 2009 A equação (13) indica que a componente de seqüência zero de pV não afeta o cálculo na etapa Backward para transformadores com uma matriz spY singular. A análise acima mostra que a equação (13) pode ser usada para calcular o vetor pV e seus componentes de seqüência não zero ( )1 2 pV + . Então, desde que ( )1 2 sV + não contem a componente de seqüência zero, pode-se escrever: [ ] ( )1 21 1 1 . 0pV + = (14) As equações (13) e (14) podem ser combinadas, resultando em: ( )1 2r r r sp p s ss sY V I Y V+ = − (15) Em que: r spY é obtida pela substituição da última linha de spY com [ ]1 1 1 , enquanto que r sI e r ssY são iguais a sI e ssY (com as últimas linhas iguais a zero), respectivamente. Com uso da equação (15), calculam-se as componentes de seqüência não zero das tensões das duas fases do lado primário, ou seja: ( )-1r r r r p sp s ss sV Y I Y V= − (16) 3) Cálculo das Injeções de Potência: A injeção de potência sobre o lado primário pode ser calculada como em [5] e [11] por: * p p pS V I= (17) Em que: p pp p ps sI Y V Y V= + (18) Note-se que a pI acima, pode ser obtida da equação (5), a qual pode ser função das tensões de ambos os lados do transformador. B. Etapa Forward 1) Cálculo da Tensão do Secundário: O cálculo da tensão no lado secundário do transformador é realizado com uso da equação (19), deduzida de (5): ( )-1 - s ss s sp pV Y I Y V= (19) Note-se que em (19) pode haver problemas, tendo em vista que ssY pode ser uma matriz singular, dependendo da conexão do transformador. A submatriz ssY tem inversa somente para as conexões Yg-Yg e D-Yg. Para a conexão Yg- Yg, as equações (7), (8) e (19) podem ser diretamente usadas para o cálculo das etapas backward / forward. Para a configuração D-Yg somente a (19) pode ser usada na etapa forward. Para todos os demais tipos de conexões aparece a singularidade de matrizes, devido a falta de um ponto de referência de tensão em um ou em ambos os lados do transformador. a) Resolvendo o Problema de Singularidade na Etapa Forward Para evitar o problema de singularidade de matrizes, as componentes de seqüência não zero das tensões podem ser determinadas. Assim, pode-se re-escrever (19) como: - ss s s sp pY V I Y V= (20) Sabendo-se que ( )1 2 sV + representa as componentes de seqüência positiva e negativa, isto é: ( )1 2 0 s s sV V V+ = − (21) Onde o vetor 0 sV é tensão de seqüência zero sobre o lado secundário. Desta forma, substituindo-se a equação (21) em (20), obtém-se: ( )( )1 2 0 ss s s s sp pY V V I Y V+ + = − (22) O produto de ssY por 0 sV é sempre zero para todas as configurações de transformadores, exceto para as configurações Yg-Yg e D-Yg. Isto acontece porque ssY é representado por YII em todas as outras configurações, portanto, de (11) e (12), pode-se obter: 0. 0sYII V = (23) Então a equação (22), exceto para as conexões Yg - Yg e D - Yg , pode ser reduzida a: ( )1 2 ss s s sp pY V I Y V+ = − (24) A equação (24) indica que a componente de seqüência zero de sV não afeta o cálculo na etapa forward para transformadores com uma matriz ssY singular. A análise acima mostra que a equação (19) pode ser usada para calcular ambos: o vetor sV , e suas componentes de seqüência não zero ( )1 2 sV + . Então, desde que ( )1 2 sV + não contem a componente de CHOQUE et al.: DISTRIBUTION TRANSFORMER MODELING 199 seqüência zero, pode-se escrever: [ ] ( )1 21 1 1 . 0sV + = (25) As equações (24) e (25), podem ser combinadas como: ( )1 2w w w ss s s sp pY V I Y V+ = − (26) Em que: ( )1 2 sV + é a componente de seqüência não zero de sV ; a submatriz w ssY é a mesma submatriz ssY , exceto que a última linha é substituída com [ ]1 1 1 ; o vetor w sI e a submatriz w spY são obtidas pelo substituição dos elementos na última linha de sI e spY por zero, respectivamente. Uma vez que as componentes de seqüência não zero de pV e sV são calculados, a componente de seqüência zero é agregada a eles, para formar a tensão linha – neutro, tal que o procedimento backward / forward possa continuar. V. TESTES E RESULTADOS Nesta seção, apresentam-se alguns testes e os resultados obtidos com três sistemas de distribuição. A. Sistema de 2 barras Neste caso, tem-se um sistema de duas barras [5] e [10], como mostrado na Fig. 2, que são ligadas com um transformador trifásico abaixador em conexão delta – estrela aterrada (D-Yg). A carga no secundário do transformador é de 400+j300 kVA. Assume-se que o transformador está operando em condições nominais, portanto os taps do primário e do secundário são iguais à unidade. A tensão na barra 1 é assumida como sendo 1.0 p.u., a carga é desbalanceada: 50% para a fase a; 30% para a fase b; e 20% para a fase c. Os dados nominais do transformador são: tensões de 13,8 kV– 208 kV; potência de 1000 kVA; Z=6%; e X/R=5. Carga 1 2 Barra Fonte D - Yg Figura 2. Sistema exemplo de dois barras Na Tabela III apresenta-se o perfil de tensão na barra de carga, estando a carga conectada em estrela. A tolerância para a convergência foi de 0,01. As magnitudes das tensões |V| estão em p.u., e os ângulos “θ” de cada fase estão em graus. Os valores mostrados na Tabela 3 são os de cada iteração até a convergência. TABELA III PERFIL DE TENSÃO DO SISTEMA DE DUAS BARRAS |V| θ |V| θ |V| θ 0 1,0 0,0 1,0 -120,0 1,0 120,0 1 0,9572 28,27 0,9819 -91,51 0,9867 148,27 2 0,9648 27,95 0,9825 -90,91 0,9934 149,11 3 0,9645 28,23 0,9794 -91,06 0,9861 149,30 Iterac. No. Fase a A Fase b Fase c B. Sistema IEEE 4 barras Na Fig. 3 mostra-se o sistema IEEE- 4 barras [14]. O principal propósito de usar este alimentador é para proporcionar diversos testes com várias conexões de transformadores trifásicos abaixadores, bem como elevadores. Os casos testados para este alimentador incluem as seguintes conexões de transformadores: (Estrela aterrado – Estrela aterrado (Yg-Yg); Estrela aterrado – Delta (Yg-D); Estrela – Delta (Y-D); Delta – Estrela aterrado (D-Yg); Delta – Delta (D-D); Estrela aberto – Delta aberto (Yo-Do)). A carga pode ser balanceada bem como desbalanceada. 1 2 3 4 Carga Barra Fonte Figura 3. Sistema IEEE-4 barras Na Tabela 4 são apresentados os valores de tensão trifásica para cada fase (a, b e c) na barra de carga (barra 4), para as várias configurações do transformador elevador entre as barras (2-3), sendo a carga balanceada, conectada em estrela e modelada como potência constante. A tensão é elevada de 12,47 kV (barra 2) para 24,9 kV (barra 3). A tolerância utilizada foi de 10-3, estando os taps tanto no primário como no secundário iguais a unidade, e aplicando-se uma tensão inicial de 1,0 p.u. na barra de referência. TABELA IV PERFIL DE TENSÃO DO SISTEMA IEEE-4 BARRAS |V| θ |V| θ |V| θ Yg - Yg 3 0,9489 -3,4711 0,9524 -123,4870 0,9510 116,4404 Yg - D 4 0,9469 -33,2231 0,9495 -153,1654 0,9492 86,6907 Y - D 4 0,9576 26,7543 0,9592 -93,3034 0,9578 146,6640 D - Yg 4 0,9541 26,3908 0,9567 -93,6639 0,9545 146,2547 D - D 3 0,9493 -3,4688 0,9522 -123,4660 0,9510 116,4071 Yo - Do 4 0,9751 -30,0690 0,9887 -150,1749 0,9804 89,2181 Conexão Trafo Iterações No. Fase a Fase b Fase c C. Sistema IEEE 34 barras Na Fig.4 mostra-se o diagrama unifilar do sistema IEEE- 34 barras, como descrito como em [13]. Na figura observa-se um transformador abaixador no trecho 832-888. O trecho 888- 890 é um trecho longo e responsável por 25% do carregamento do sistema. 200 IEEE LATIN AMERICA TRANSACTIONS, VOL. 7, NO. 2, JUNE 2009 800 806 808 812 814 810 802 850 818 824 826 816 820 822 828 830 854 856 852 832 888 890 838 862 840 836860834 842 844 846 848 864 858 Figura 4. Sistema IEEE-34 barras. Na Fig. 5 é apresentado o perfil de tensão do sistema IEEE- 34 barras, para a configuração do transformador trifásico abaixador em estrela aterrado – estrela aterrado (Yg-Yg), que abaixa a tensão de 24,9 kV para 4,16 kV. A barra terminal (890) apresenta um péssimo perfil de tensão, como pode ser visto na Fig. 5, com as cargas do sistema desbalanceadas em cada fase (a, b e c). Os resultados mostrados neste gráfico foram obtidos sem consideração dos dois reguladores de tensão trifásicos. Na subestação (barra 800) a tensão inicial foi de 1,05 p.u. (barra de referência). A tolerância usada para a convergência foi de 10-4. O número de iterações obtidas para alcançar a convergência foi igual a 7. Observa-se também na Fig.5, que a curva desenhada foi somente para as barras (800, 802, 806, 808, 812, 814, 850, 816, 824, 828, 830, 854, 852, 832, 888 e 890), apenas para ilustrar as tensões nessas barras do ramal principal. Perfil de Tensão 0,75 0,85 0,95 1,05 800 802 806 808 812 814 850 816 824 828 830 854 852 832 888 890 Barras Te ns ão (p u) fase A fase B fase C Figura 5. Perfil de Tensão do sistema sem os reguladores. Na Fig. 6 é apresentado o perfil de tensão do sistema IEEE- 34 barras, com a configuração do transformador trifásico abaixador em estrela aterrado – estrela aterrado (Yg-Yg), com as mesmas condições anteriores. Neste caso aplicaram-se os dois reguladores de tensão encontrados entre as barras (814 – 850), e entre as barras (852-832). O número de iterações obtidas para alcançar a convergência foi de 18. Neste caso pode-se verificar que a tensão na barra (890) melhorou devido aos reguladores, atingindo um nível de tensão razoável. Perfil de Tensão 0,75 0,85 0,95 1,05 800 802 806 808 812 814 850 816 824 828 830 854 852 832 888 890 Barras Te ns ão (p u) fase A reg2 fase B reg2 fase C reg2 Figura 6. Perfil de Tensão do sistema com a atuação dos dois reguladores. VI. CONCLUSÕES Neste trabalho foi proposta uma metodologia para modelar matematicamente os transformadores trifásicos de distribuição, considerando suas diversas configurações como abaixador ou elevador, para aplicação em fluxo de potência trifásico para sistemas de distribuição radial. A metodologia proposta apresenta uma formulação em valores em p.u., para representar o transformador trifásico e suas diferentes configurações encontradas dentro de um sistema de energia. Apresentaram-se detalhadamente as tabelas para as conexões do transformador abaixador bem como elevador, fazendo uso das três submatrizes YI , YII e YIII (derivadas da Ybus trafo) para montar a conexão desejada, nas quais pode-se levar em conta a defasagem angular e a incidência de taps. Um dos problemas da aplicação dos modelos de transformadores em fluxo de potência é a singularidade de matrizes. Neste trabalho foi mostrado, com detalhes, quando isto ocorre e uma forma eficiente de resolver este problema. Finalmente foram apresentados resultados de cálculo de fluxo de potência, que foram obtidos considerando a modelagem apresentada. APÊNDICE A: Transformador Estrela aterrado – Estrela aterrado. A Fig.A1 mostra o modelo do circuito acoplado do transformador estrela aterrado – estrela aterrada. 2R 2R 2R B 2L 2L 2L 4i 1R 1L1i a 1R 1L3i b 1R 1L5i c 2i 6i A C M M M Figura A.1. Modelo do circuito acoplado do transformador estrela aterrado – estrela aterrada (Yg-Yg). CHOQUE et al.: DISTRIBUTION TRANSFORMER MODELING 201 A matriz bus trafoY⎡ ⎤⎣ ⎦ do transformador mostrada na Tabela 1, para este conexão é detalhada a seguir, onde as quatro submatrizes são não singulares. Os taps do lado primário do transformador “α” e os taps do lado secundário “β” podem ser diferentes da unidade. ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ − − − − − − = 2 tt 2 tt 2 tt t 2 t t 2 t t 2 t trafobus β y00 αβ y00 0 β y00 αβ y0 00 β y00 αβ y αβ y00 α y00 0 αβ y00 α y0 00 αβ y00 α y Y APÊNDICE B: Transformador Delta – Estrela aterrado. A Fig.B1 mostra o modelo do circuito acoplado do transformador delta – estrela aterrada. M M M 1i 3i 5i 1R 1R 1R 2R 2R 2R a b c B A C 1L 1L 1L 2L 2L 2L 4i 2i 6i Figura B1. Modelo do circuito acoplado do transformador delta – estrela aterrada (D-Yg). A matriz bus trafoY⎡ ⎤⎣ ⎦ do transformador elevador mostrada na Tabela 2, para este conexão é detalhada a seguir, onde as três submatrizes ,pp psY Y e spY são singulares, e a submatriz ssY é não singular. Os taps do lado primário do transformador “α” e os taps do lado secundário “β” podem ser diferentes da unidade. ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎢ ⎣ ⎡ − − − −−− −−− −−− = 2 ttt 2 ttt 2 ttt tt 2 t 2 t 2 t tt 2 t 2 t 2 t tt 2 t 2 t 2 t trafobus β y00 αβ3 y0 αβ3 y 0 β y0 αβ3 y αβ3 y0 00 β y0 αβ3 y αβ3 y αβ3 y αβ3 y0 3α 2y 3α y 3α y 0 αβ3 y αβ3 y 3α y 3α 2y 3α y αβ3 y0 αβ3 y 3α y 3α y 3α 2y Y REFERÊNCIAS [1] C. S. Cheng, and D. Shirmohammadi, “A three-phase power flow method for real-time distribution system analysis”, IEEE Transaction on Power Systems, vol. 10, no. 2, pp. 671-679, May 1995. [2] R. M. Ciric, A. Padilha-Feltrin, and L. F. Ochoa, “Power flow in four- wire distribution networks-general approach”, IEEE Transaction on Power Systems, vol. 18, no. 4, pp. 1283-1290, Nov. 2003. [3] D. Shirmohammadi, H. W. Hong, A. Semlyen, and G. X. Luo, “A compensation-based power flow method for weakly meshed distribution and transmission networks”, IEEE Transaction on Power Systems, vol. 3, no. 2, pp. 753-762, May. 1988. [4] B. Alencar Souza, H. D. de Mac~edo Braz, J. M. Cavalcante de Albuquerque, J. G. Gerlach Gutteres, “Radial Distribution Systems Power Flow with Distributed Generation: Modified Power Summation Method”, IEEE LATIN AMERICA TRANSACTIONS, Vol. 4, No. 3, pp. 192-197, May 2006. [5] Z. Wang, F. Chen, and J. Li, “Implementing transformer nodal admittance matrices into backward / forward sweep-based power flow analysis for unbalanced radial distribution systems”, IEEE Power Transaction on Power Systems, vol. 19, no. 4, pp. 1831-1836, Nov. 2004. [6] Peng Xiao, and David C. Yu, “A unified three-phase transformer model for distribution load flow calculations”, IEEE Transaction on Power Systems, vol. 21, no. 1, pp. 153-159, February 2006. [7] M. S. Chen and W. E. Dillon, “Power system modeling”, IEEE Proc., pp. 901-915, July 1974. [8] J. Arrillaga, C. P. Arnold, and B. J. Harker, “Computer modeling of electrical power systems”, IEEE PES. 2nd ed. New York: Wiley, 2001. [9] W. H. Kersting, “Distribution system modeling and analysis”, Boca Raton, Florida, CRC Pres, 2002. [10] T. H. Chen, M. S. Chen, T. Inoue, P. Kotas, and E. A. Chebli, “Three- phase cogenerator and transformer models for distribution system analysis”, IEEE Transaction on Power Delivery, vol. 6, no. 4, pp. 167- 1681, Oct. 1991. [11] D. E. Rodas, “Modelagem de transformadores para sistemas de distribuição”, ETELEC- UNESP, Ilha Solteira, Aug. 2005. [12] A. Tan, W. H. E. Liu, and D. Shirmohammadi, “Transformer and load modeling in short circuit analysis for distribution systems”, IEEE Transaction on Power Systems, vol 12, no. 3, pp. 1315- 1322, Aug. 1997. [13] G. X. Luo and A. Semlyen, “Efficient load flow for large weakly meshed networks”, IEEE Transactions on Power Systems, vol. 5, no. 4, pp. 1309-1316, Nov. 1990. [14] W. H. Kersting, "Radial distribution test feeders," IEEE Transactions on Power Systems, vol. 6, no. 3, pp. 975-985, Aug. 1991. José L. Choque graduado da Universidade Nacional de São Agustín de Arequipa (U.N.S.A.) – Arequipa – Peru, em 2001. Atualmente ele é um estudante de M.Sc. em Sistemas Elétricos de Potência na UNESP – Ilha Solteira, SP, Brazil. Seu principal interesse é analise de sistemas de distribuição. Dario E. Rodas Rendón Engenheiro Eletricista da Universidade Tecnológica de Pereira (UTP), Colômbia, 1975, Especialista em sistemas de transmissão e distribuição de energia elétrica na Universidade dos Andes, 1989. Atualmente ele é professor na UTP. Sua área de interesse é planejamento em sistemas de energia elétrico, maquinas elétrica. Antonio Padilha Feltrin (SM’06) cursou engenharia elétrica na EFEI,obteve M.Sc. e Ph.D. da UNICAMP, Brazil. Atualmente, é professor titular da UNESP – Ilha Solteira, SP. De 1995 até 1997, foi prof, visitante no departamento ECE da Universidade de Wisconsin – Madison. Sua área de interesse é em analise e controle de sistemas de potência. 202 IEEE LATIN AMERICA TRANSACTIONS, VOL. 7, NO. 2, JUNE 2009