111

Scientia et Technica Ao XII, No 30, Mayo de 2006 UTP. ISSN 0122-1701

OPERACIN DE TAPS EN TRANSFORMADORES DE DISTRIBUCIN PARA

REDUCCIN DE PRDIDAS
RESUMEN
El presente artculo muestra una metodologa que busca disminuir las prdidas
de tipo tcnico en los ramales secundarios de los sistemas de distribucin en
sistemas trifsicos equilibrados. Su aplicacin en sistemas desbalanceados
conserva la misma lgica.

ALEXANDER MOLINA C
Ingeniero Electricista
Profesor
Universidad Tecnolgica de Pereira
almo@utp.edu.co
DARIO ELIECER RODAS
Ingeniero Electricista
Profesor Titular
Universidad Tecnolgica de Pereira
drodas@utp.edu.co

Se considera el transformador con cambiador de derivaciones y distintos tipos de

carga. Para la asignacin del Tap ptimo se usan reglas difusas que consiguen
minimizar las perdidas sin violar los lmites de tensin. Se presentan resultados
sobre un sistema de prueba donde se mejoran los valores de prdidas cuyo costo
se reduce al pago de personal que efecte las modificaciones del Tap.

DIEGO GONZALEZ
Ingeniero Electricista
Profesor
Universidad Tecnolgica de Pereira
diegogo@utp.edu.co

PALABRAS CLAVES: sistemas de distribucin secundarios, prdidas.

ABSTRACT
This paper shows a methodology for reduce technical losses on secundary
systems. The Tap is included in the trafo model. In the Tap assignation Fuzzy
Logic is used to get losses disminution whitout limits violation. A test system
shows the advantage of the methodology in a better value of energy losses which
is reached with low cost.
KEYWORDS: Distribution systems, technical losess.
1. INTRODUCCIN
En Colombia los sistemas de distribucin son
responsables del mayor porcentaje de las prdidas
tcnicas. Las metodologas empleadas de manera
tradicional dan soluciones parciales en niveles de alta
tensin descuidando el problema en baja tensin [1],[2].

ser reformulado abandonando su solucin mediante

tcnicas de optimizacin y a cambio se emplee reglas
difusas.
En la seccin 4 se presenta la metodologa con su
respectivo flujograma y en 5 se exponen de algunos
resultados particulares.

En los niveles bajos de voltaje se plantea tan solo el

cambio de acometidas para el control de prdidas y el
Tap pocas veces se ha considerado como variable para
control de ellas. Tan solo [6] analiza esta posibilidad y su
aplicacin se efecta a niveles de alta tensin y en
sistemas automatizados, caso contrario a las posibilidades
regionales.

2. MODELAMIETO DEL PROBLEMA

En las ecuaciones 1 y 2 se presentan la funcin objetivo y
las restricciones del problema respectivamente. Se inicia
con la formulacin clsica como un problema de
optimizacin para luego justificar el uso de un modelo
desacoplado.
2.1 Funcin objetivo

Dado que el Tap (relacin de espiras) controla los niveles

de tensin en el secundario y estos niveles de tensin
determinan corrientes; es posible establecer su relacin
con los estados de prdidas. Es as como mediante la
articulacin de modelos adecuados de la carga y del
transformador de distribucin y de la carga se consigue
modelar el comportamiento general de las
En el numeral 2 se presenta el modelamiento general del
problema y algunas consideraciones a tener en cuenta
para facilitar su tratamiento.
En el numeral 3 se analiza el sistema secundario, nivel de
baja tensin y se obtienen algunas expresiones e
informaciones claves para el estudio del problema. Aqu
toma gran importancia el hecho que el problema pueda
Fecha de Recepcin: 31 Enero de 2006
Fecha de Aceptacin: 7 Abril de 2006

min k

hTm Ril ( I (Vl )il _ t )

T =T
l =1 i =1
w

Tn

vl

(1)

+ O( Mov)
Esta incluye el costo de las prdidas y el costo de mover
los Taps en O(movimientos). La suma se realiza para los
Tn perodos, en los w transformadores que alimentan
ramales secundarios y para los vl tramos de conductor de
cada ramal secundario.
k es la constante $/kWh. Se considera invariable para
todos los intervalos de la curva de carga de un mismo
ao. hTm es el ancho en horas de cada intervalo Tm que
compone la curva de carga anual discretizada.

Scientia et Technica Ao XII, No 30, Mayo de 2006. UTP

112
R

l
i ,

( I ( V l ) il _ T )

son los valores de Resistencia y

corrientes, respectivamente, asociados a la lnea i-sima,
del ramal l-simo en cada intervalo Tm de la curva de
carga
2.2 Restricciones
vl

Snodoil _ Tm = Dnodoil _ Tm
i =1

Vmin < Vil < Vmax

Vslack _ l =

V prim _ l
Tapl

(2)

Moverl [0,1]
Tapl min Tapl Tapl max

Tapl = {Tapl min , " ,1 paso, 1 ,1+paso," , Tapl max }

vl

Snodoil _Tm = Dnodoil _Tm --

Considera primera y segunda

i =1

ley de Kirchhoff. El cumplimiento de estas restricciones

las garantiza el flujo de carga. La suma de las potencias
que llegan al nodo i-simo del ramal l-simo es igual a
la demanda de dicho nodo.
Vmin < Vil < Vmax

-- Tiene en cuenta el cumplimiento de las

restricciones de voltaje para cualquier nodo.

Vslack _ l =

V prim _ l

-- Aqu se relacionan los voltajes en

Tapl
primario y secundario con el valor de Tap.
Moverl [0,1]

-- Considera la posibilidad de mover o no

el Tap del transformador l-simo. 0 no se mueve; 1, si
permite mover. Esto para incluir transformadores en los
que no se puedan cambiar las posiciones.
Tapl = {Tapl min , " , 0.975 , 1 , 1.025 ," , Tapl max } --Considera

el

movimiento de los Taps dentro de un rango especificado,

al igual que se obliga a que ellos tomen un valor discreto
determinado por el paso que se establezca por norma
tcnica [10].

de Taps y el nmero de transformadores donde es posible

moverlo, es decir, (#Posiciones_Tap) ^(w)
Al momento de efectuar los flujos de carga en los niveles
secundarios se observa que los cambios en las prdidas
no superan el 3% de las iniciales, de manera que los
voltajes aguas arriba no sufren mayor modificacin. En
general no se presentan cambios substanciales en el
estado del sistema al mover el Tap y esto puede
verificarse analizando el Jacobiano de los sistemas de
distribucin [4].
Con base en lo anterior se propone desacoplar el ramal
secundario del primario de manera que la funcin
objetivo se puede presentar como la suma de los mnimos
de cada subproblema. Esto porque el estado en cada
ramal depende slo del voltaje que tenga en su lado
primario y este valor es entregado por el flujo de carga
como condicin inicial.
min k

Tn

vl

l
l
Ri ( I (Vl )i _ t )
hTm
T =T
l =1 i =1
m

w
vl
hTm min Ril I (Vl )il _ t

T =T
l =1
i =1
m 1
Tn

(3)

Con las anteriores consideraciones y teniendo en cuenta

la ecuacin 3, se plantea que el mnimo de la suma de
las prdidas en todos los ramales es igual a la suma de las
mnimas prdidas en cada uno de ellos, permitiendo
estudiar cada problema de manera separada.
3. ANLISIS DEL SISTEMA SECUNDARIO
3.1 Modelo del transformador
El modelo del transformador se presenta a continuacin
para el caso monofsico. Figura 1. Las variaciones en las
espiras del secundario son nulas. Este modelo del
transformador en pu se toma de [7] y permite integrarlo
al flujo de carga. Las variaciones en el Tap no tienen
mayor influencia en las prdidas del mismo ni en el flujo
de reactivos [8].
j
ZTH
+

VTH

Vj
-

2.3 Caracterstica del modelo

En los sistemas de distribucin el nmero de
transformadores que alimentan ramales secundarios es
grande, con las consiguientes consecuencias para el
espacio de soluciones. Adicionalmente la variable
principal de este problema, el Tap, es discreta.
Los trminos de la funcin objetivo se componen de las
prdidas en cada uno de los ramales secundarios. Dichas
prdidas dependen del valor de los voltajes en secundario
de los transformadores y stos a su vez dependen del
voltaje entregado por el flujo en el lado de alta tensin.
El espacio de soluciones tiene crecimiento No-Polinomial
determinado por las posibles combinaciones del nmero

Figura 1. Modelo equivalente del transformador, con tap, usado

como slack del ramal secundario.

VTH

ZTH

V prim

p.u.

Z t p.u.

Tap

Tabla 1. Valores sin variacin de espiras en baja.

3.2 Modelo de las cargas

V
V
Si = Pn i +j Qn i
Vn
Vn

(4)

Scientia et Technica Ao XII, No 30, Mayo de 2006. U.T.P

113

El modelo utilizado por la mayora de analistas en las

empresas distribuidoras es el de modelo de carga tipo
potencia constante. En esta ocasin los modelos que
incluyen cargas tipo impedancia constante, potencia
constante y corriente constante se incluyen de manera
que se vean los cambios que se presentan en diferentes
estados de voltaje. Los subndices i y n indican estados
generales y nominales respectivamente. Se opt por el
modelo exponencial, ecuacin 3, y una generalizacin
mayor se hizo pues se consider, para los casos de
prueba, variacin en igual sentido de la potencia activa y
reactiva, es decir, se asumen cargas con valores =
[1].

El posible espacio de soluciones de un sistema de

distribucin con dos transformadores, alimentando cada
uno a un sistema secundario, se puede representar por la
figura 5. En cada uno de los ejes horizontales se
encuentra la variacin del Tap de cada ramal y en el eje
vertical se presentan los valores de prdidas.

3.3 Modelo de las lneas

Las lneas se modelan de acuerdo a lo planteado en [5].
3.4 Modelo matemtico de las prdidas
Se estudia el modelo del sistema que se presenta en la
figura 2 donde se denomina PL al valor de las prdidas
del sistema. El subndice n indica valore nominal y el
subndice i indica otro estado.

Figura 3. Valores de prdidas normalizados para diferentes

valores de

Prdidas

Figura 2. Sistema de dos nodos.

En la ecuacin 8 se determina el comportamiento del

valor de prdidas para distintos valores de . En un valor
de = 1 las prdidas permanecen constantes pues el
exponente de Vi se hace cero. Con un valor de > 1 si Vi
crece, crece PL. En el caso en el que < 1 el exponente
de Vi se hace negativo, lo que hace que el
comportamiento de las prdidas tengan relacin inversa
con el voltaje.
PL = RI i2
Vi
2
Sn
V
Si
PL = R = R n
Vi

Vi

Vi

2 Vn
PL = R * S n
Vi 2
PL =

(5)

Tap 1

Tap 1

(6)

(7)

Ahora, suponiendo una cada despreciable entre V y Vi,

es decir; V Vi se tiene:

R * Pn2 2( 1)
Vi
Vn2

Tap 2

Tap 2

Aunque se aprecia que los cambios en las prdidas

obedecen a la formulacin planteada en la ecuacin 8, se
realizan ensayos sobre un sistema de prueba para
confirmar lo planteado.

2
n

PL =

Tap 3
Tap 3

Figura 4. Espacio de soluciones para sistema de 2 ramales.

R * S 2( 1)
V
Vn2

Ramal 1

Ramal 2

(8)

La anterior expresin muestra la relacin existente entre

el voltaje que alimenta los ramales, el modelo de carga y
las prdidas. La figura 3 ilustra relaciones entre prdidas
y voltajes para distintos valores de .

3.5 Sistema de prueba

Se estudia un caso que valida la ecuacin 8 y permite
apuntar elementos para la obtencin de la poltica de
manejo de Taps. Se genera un sistema de prueba
secundario que es alimentado por una seal de voltaje de
120 Voltios. Este mismo voltaje se establece como
voltaje nominal. El sistema de prueba tiene una carga de
41 kW distribuida de manera uniforme en sus nodos y
valores de impedancia de lneas de Z=1,017+j0,097125
(Ohm/km). Los tramos que unen los nodos tienen una
longitud de 5 metros. Para la impedancia de cortocircuito
del transformador se uso
un valor tpico de
Zt=0,0138+j0.02785 referidas al lado de baja tensin. La
topologa se presenta en la figura 5.

Scientia et Technica Ao XII, No 30, Mayo de 2006. UTP

114
Ntese en la tabla 2 que para el caso en el que el valor del
exponente es de 0.5, se mudan las prdidas desde un
valor de 1,1812 kW hasta 0,9894. En el caso del valor
del exponente 1.5 se encuentran prdidas de 0,94706 kW
en 114 V y a 126 V se encuentran prdidas de valor de
1,0483 kW. Cuando el exponente es igual a 1 las prdidas
permanecen constantes. Esto ya haba sido visto desde la
aparicin del trmino ( 1) en la ecuacin 7 y 8.

3.6.3. Carga tipo potencia constante o exponente

menor que 1
Para las cargas modeladas como potencia constante o en
general con valor de exponente menor que 1, las prdidas
mnimas se presentan cuando el valor de voltaje se ubica
en su valor mximo. Al igual que para los dems tipos de
carga, los valores de voltaje estn restringidos.

Figura 5. Sistema de prueba de 42 nodos.

3.6.4. Ubicacin de tap

Para la aplicacin de la poltica de manejo de Tap se
construye un sistema de inferencia difuso que represente
las reglas lingsticas mencionadas en los numerales
3.6.1 a 3.6.3. Esta funcin tiene como entradas el tipo de
carga que domina el ramal secundario, figura 6 y el valor
del voltaje en el lado primario del transformador, figura
7. El tipo de carga de cada ramal se puede determinar a
travs de una ponderacin de los tipos de carga que a l
pertenezcan.

Voltaje

114

117

120

123

126

1,2260

1,1594

1,0982

1,0118

0,9708

0,5

1,1812

1,1209

1,0623

1,0351

0,9894

1,0289

1,0289

1,0289

1,0289

1,0289

1,5

0,9470

0,9724

0,9977

1,0230

1,0483

0,8740

0,9206

0,9684

1,0175

1,0677

Tabla 2. Prdidas en kW en el sistema de 32 nodos.

Figura 6. Funcin de Inferencia de la entrada del modelo de

carga (exponente).

3.6 Poltica operativa de taps

Ayudados de los resultados del sistema de prueba, tabla 2
y de la ecuacin 8 se puede afirmar que para cada ramal
es posible encontrar un valor de Tap que minimice las
prdidas y a su vez cumpla con los lmites de voltaje. De
aqu se desprende que la reduccin de prdidas se puede
obtener identificando el tipo de carga en cada ramal y
asignando un valor de Tap adecuado que permita cumplir
restricciones de voltaje.
3.6.1 Carga tipo corriente constante, exponente igual
a1
En las cargas modeladas con exponente 1, el valor de las
prdidas permanece constante a pesar de los cambios en
el nivel de tensin. Se concluye que en las cargas de tipo
corriente constante se puede variar el voltaje, el Tap, de
manera que se cumpla con las restricciones de voltaje y
en ningn caso tales movimientos afectarn las prdidas
[1],[5].
3.6.2. Carga tipo impedancia constante o exponente
mayor que 1
Se puede afirmar que para cargas con valor de exponente
mayor que 1, las prdidas mnimas en ramales
secundarios se presentan cuando el valor del voltaje se
ubica en su valor mnimo, es decir, cuando el Tap est en
su mximo. Dicho valor mnimo de voltaje est
restringido por los lmites de regulacin.

Figura 7. Funcin de inferencia para la entrada de voltaje.

En el sistema difuso que se menciona, el tipo de carga

permite establecer en que direccin se debe mover el
Tap, es decir, si en direccin del aumento de voltaje o al
contrario. Aparte de dicha direccin, es necesario
establecer qu tan grande debe ser dicho avance. Lo
anterior porque no se puede asignar el mismo Tap a dos
transformadores que, aunque estn alimentando el mismo
tipo de carga, presenten voltajes primarios diferentes. Es
decir; cada cierto de intervalo de voltajes en el primario y
cada uno de los tipos de carga requiere una posicin de
Tap diferente que se asigna con el sistema de inferencia
difuso.
Dos ejemplos de reglas del sistema pueden ser las
siguientes:
1.

If(exponente>1)and(Voltaje= bajo)then

t = 1

Scientia et Technica Ao XII, No 30, Mayo de 2006. U.T.P

2.

If(exponente>1)and(Voltaje=alto)then

115

t = 1.025

Ntese que las anteriores reglas, an teniendo exponentes

de la misma caracterstica, el hecho que el voltaje en
primario sea diferente cambia el valor de Tap a asignar.
Mientras la regla 1 mantiene el voltaje bajo, asignando un
valor de 1 al Tap, en la regla 2 dado que el voltaje se
encuentra en valor alto y que la carga es de tipo
impedancia constante, asigna un Tap de 1.025, para
disminuir su voltaje y de esta manera bajar las prdidas.
4. METODOLOGA
Con los anteriores elementos presentados y ayudados del
anlisis elctrico del problema, se abandona la idea de
solucionar el modelo a travs de una metodologa de
optimizacin y se busca la solucin a travs de reglas.
Dado que es posible obtener para cada ramal secundario
un estado de prdidas mnimo basado en una posicin de
Tap se replantea la estrategia

Tap en un valor tal que en los ramales se presenten las

mnimas prdidas.
La metodologa supone conocimiento pleno de la
informacin que se requiere de cada uno de los ramales
secundarios. Su flujograma se presenta en la figura 9.
Entrada de datos

Determinacin del estado inicial del

sistema. Flujo de carga

Anlisis de ramales secundarios.

Restricciones.
Prdidas iniciales.

Implementacin de la poltica de
manejo de taps.

no

Primera
Iteracin?

si

En primera instancia se determina el estado del sistema a

nivel de alta tensin. En esta primera etapa se concentran
las cargas, junto con el valor de prdidas iniciales, en los
nodos primarios de los transformadores. Estas
condiciones permiten aplicar las reglas generales de
asignacin de Tap de 3.6.4. Una vez se cuenta con los
voltajes en los nodos primarios y ayudado de la
informacin que describe el tipo de carga que cada uno
de los transformadores alimenta, se ubica la nueva
variable del transformador. Con la nueva posicin de Tap
se ejecutan los flujos de carga en los ramales secundarios
para verificar el cumplimiento de las restricciones. En
caso de violarse alguna, entra en funcionamiento una
nueva etapa de correccin que acta como se muestra.

CORRIGE TAP
ASIGNA TAP

Implementa reglas.

Evita restricciones
Implementa reglas.

Anlisis de ramales secundarios.

no

Restricciones.
Prdidas iniciales.
Costos

cumple
restricciones?

si
Salida de Datos
Costos, Tap, Prdidas, Ahorros

FIN

t t paso _ tap
(Vmnimo_ramal l < Vmnimo) then =
t t + paso _ tap
If (Vmximo_ramal l > Vmximo) then =

If

En la figura 8 se presentan el estado inicial en 1. En el

proceso de asignacin de tap es posible que se viole la
restriccin ubicando el voltaje en el punto 2. Al corregir
el Tap, el voltaje regresa a la regin factible. A diferencia
del Flujo de Potencia ptimo [9], que fija el voltaje en el
lmite aprovechando la continuidad de las variables, esta
correccin lleva al cumplimiento pero no lo fija en la
frontera.
Correccin
3

Asigna Tap

Violacin de
restriccin

Lmite

Figura 8. Correccin de las restricciones

El criterio de parada del algoritmo es el cumplimiento de

las restricciones dado que las polticas siempre ubican el

Figura 9. Flujograma de la metodologa

5. PRUEBAS Y RESULTADOS
El sistema IEEE de 33 nodos modificado se muestra a
continuacin. Los ramales secundarios con su respectivo
transformador se conectan en cada uno de los nodos del
sistema. El tipo de carga usado para esta prueba, en todos
los nodos, tiene exponente menor que 1. En otros casos
estudiados, que no se referencian aqu, se logr reducir
los niveles de prdidas.
Para el caso de prueba se usa variaciones en el Tap de
2.5% dado que son transformadores usados a niveles de
13200 v. Para transformadores en niveles de tensin
superiores las variaciones son generalmente de 1.25%.
La tabla 3 muestra cmo se cambia el nivel de prdidas
para un perodo de 1 ao, con curva de carga discretizada
en 3 bloques y valor de 75$kWh. En la tabla 4 se nota
que los valores de Tap en los transformadores toma
valores bajos. Esto es coherente con la necesidad de la
carga de tipo potencia constante, dado que un valor bajo

Scientia et Technica Ao XII, No 30, Mayo de 2006. UTP

116
de Tap hace que los valores de voltaje sean altos con la
consiguiente reduccin de corriente por las lneas.
El ejercicio se repiti para el mismo sistema adicionando
restricciones de movimiento en los nodos 2, 5, 8, 10, 1318, 22, 25, 28 y 33. En esta ocasin los ahorros fueron
menos favorables.

Generalmente los cambiadores de derivacin pueden

moverse sin problemas an en transformadores antiguos.
Por otro lado los transformadores de distribucin actuales
incluyen la variacin del Tap, lo que no implica costos
adicionales a la empresa que opera la red.
6. CONCLUSIONES
La metodologa planteada mostr ser capaz de reducir las
prdidas tcnicas en los ramales de prueba secundarios
estudiados. La restriccin de movimiento de Tap que se
puede presentar en transformadores defectuosos evita que
las prdidas disminuyan su valor .
El movimiento de Tap en transformadores de distribucin
se debe manejar de manera cuidadosa y planeada y no
permitir su movimiento a discrecin de los miembros de
cuadrillas en los sistemas de distribucin reales.

Figura 12. Sistemas de 33 nodos

Prdidas Ini

Prdidas fin

Ahorro%

Ahorro

735426,143 kWh

587635,6133 kWh

20,09

10345337,07 $

Tabla 3. Ahorros en el sistema

Nod

Tap Ini

Tap Fi

Mov

Nod

Tap Ini

Tap Fi

Mov

1,05

0,9

18

0,95

0,95

1,05

0,9

19

0,95

0,9

0,9

0,9

20

0,975

0,9

0,95

0,9

21

1,025

0,9

0,925

0,9

22

1,025

0,9

1,025

0,9

23

0,9

1,025

0,9

24

0,9

1,05

0,9

25

1,05

0,9

10

0,9

26

1,025

0,9

11

0,9

27

1,05

0,9

12

1,05

0,9

28

1,025

0,9

13

0,9

29

0,9

0,9

14

0,95

30

0,9

0,9

31

1,025

0,9

15

1,05

0,975

16

0,95

32

0,9

17

0,9

0,95

33

0,9

Tabla 4. Resultados del empleo de la metodologa.

Prdidas Ini

Prdidas fin

735426,143 kWh

656388,4686 kWh

Ahorro
%
10,75

Ahorro
dinero
5532637,206 $

Tabla 5. Ahorros con restriccin en algunos trafos.

Se implement esta metodologa junto con otras que se

aplican para corregir prdidas a nivel primario como
inyeccin de reactivos, mostrando ser compatibles y
complementarias ya que los ahorros hechos en alta
tensin se pueden sumar con los hechos en baja tensin.
Un elemento adicional es que se corrigen las violaciones
de los voltajes en los secundarios de los transformadores.

La metodologa se implement suponiendo sistemas

balanceados. A nivel de 120 V los sistemas tienen un alto
desbalance pero la lgica que propone la metodologa
sigue siendo vlida porque a pesar de las diferencias de
voltaje de las fases, la tendencia para mover el Tap es
igual para ellas. Por otro lado, no existe Tap para cada
fase.
7. BIBLIOGRAFA
[1] IEEE Task Force on Load Representation.
Bibliography on load models for power flow and
dynamic
performance
simulation,
IEEE
Transactions on Power Systems, Vol. 10, No. 1.
February 1995.
[2] Garcs A., Reconfiguracin de Alimentadores
primarios usando una tcnica combinatorial bajo un
modelamiento trifsico, Trabajo de grado Magster
en Ingeniera elctrica. UTP. PIE. 2005.
[3] Turan G..Electric Power Distribution system.
Editorial Mc-Graw Hill. 1986.
[4] Grainger, J., Stevenson W. Anlisis de Sistemas de
Potencia. Editorial Mc-Graw Hill. 1996
[5] Rodas D. Modelamiento de cargas en sistemas de
distribucin . Scientia et Tcnica. N 7 Mayo de
1998.
[6] Augugliaro A, et al. Voltaje regulation and power
losses minimization in automated distribution
networks by an evolutionary multiobjective
approach, IEEE transactions on systems.Vol 19.
N3. August 2004.
[7] Rodas D, et al. Transformador monofsico de
distribucin y de potencia, modelamiento con
TAPS, Scientia et Tcnica. N 30. 2005.
[8] Molina A., Poltica operativa de transformadores de
distribucin, Trabajo de grado Magster en
Ingeniera elctrica. UTP. PIE. 2005.
[9] Tinney W., Dommel H. Optimal power flor
solutions. IEEE trans. Power App. Sist. Vol Pass87. Oct. 1968.
[10] NTC 737 Transformadores Especificaciones de
devanados y sus derivaciones. NTC 2003.