IEEE INTERCON 2002 Optimizacin de Diseo de Transformadores

Area de Ingeniera Elctrica Mquinas Elctricas: Optimizacin en el Diseo de Mquinas Elctricas Ttulo: Optimizacin de Diseo de Transformadores Autor : Manuel Carranza Arvalo CIP 4152 Universidad Nacional de Ingeniera Direccin: Telfono: Los Abedules 215, La Molina 435 9127 1. Introduccin El Transformador de Potencia, como dispositivo electromagntico destinado a la conversin de los factores de la potencia elctrica, consiste, esencialmente en el enlace fsico de un circuito magntico con un circuito elctrico. La potencia pasante est determinada por el producto de las reas de conduccin del flujo magntico y la corriente elctrica. Es posible, por tanto, obtener soluciones tcnicas en una variedad de combinaciones de cantidad de materiales bsicos y solicitaciones especficas. Esto es: Diferentes alternativas de pesos y prdidas, tanto para el fierro magntico, como para el cobre. Los pesos de los materiales inciden directamente en el costo de adquisicin y las prdidas, en los gastos operativos del transformador La alternativa ms econmica para el usuario ser aquella que arroje un mnimo para el costo de adquisicin, ms el valor actualizado de los costos de las prdidas durante la explotacin. Este trabajo se refiere a los transformadores trifsico de columnas cilndricas. Presenta un enfoque relativamente simple y completo para evaluar los costos de la parte activa as como aquellos de las prdidas en servicio. Se indica una solucin analtica mediante el mtodo de las derivadas, y un procedimiento de clculo numrico para determinar las dimensiones bsicas de la solucin: Dimetro de columna, Ancho y Altura de ventana. Este procedimiento ha sido plasmado en un programa ejecutable en PC. Tambin se presenta un grfico de sensibilidad donde se aprecia la importancia de ubicarse en el dimetro adecuado, quedando cierta holgura para escoger el ancho de ventana y altura de la columna. Indice 1 2 3 4 5 6 7 8 9 Introduccin Simbologa Pesos Costos Procedimiento Matemtico Procedimiento numrico Programa de aplicacin Conclusiones Bibliografa

Pag

1/7

E- mail [email protected] Resumen En el Transformador Elctrico la potencia pasante est determinada por el producto de las reas de fierro y cobre. Es posible, por tanto, obtener soluciones tcnicas en una variedad de combinaciones de cantidad de materiales bsicos y solicitaciones especficas. Los pesos de los materiales inciden directamente en el costo de adquisicin y las prdidas, en los gastos operativos del transformador Este trabajo se refiere a los transformadores trifsico de columnas cilndricas Presenta el clculo de los pesos de fierro y cobre a partir de las dimensiones bsicas del ncleo: Dimetro, Ancho y Altura de ventana, as como los factores de plenitud de las reas respectivas. Para la potencia, se considera, adicionalmente las densidades de flujo magntico y de la corriente elctrica. Para la evaluacin econmica se considera los costos de materiales, prdidas especfica y coeficientes de capitalizacin de la prdidas. Se presenta la solucin analtica mediante el mtodo de las derivadas y un procedimiento numrico convergente, para determinar el Dimetro, Altura y Ancho de ventana ptimos. En forma grfica se muestra el efecto de las desviaciones. Se incluye tambin un programa de aplicacin que permite resolver los casos de costo de material mnimo, o de costo de operacin mnimo. ________________________________________ ODT02 / Manuel Carranza A. / 03 Mayo 2002

IEEE INTERCON 2002 Optimizacin de Diseo de Transformadores

Wf 2. D H A Ff Fc Vf Vc U I F B J Sf Sc P E Sc f c Gf Gc Cf Cc qf qc Simbologa. Dimetro del Ncleo Altura de Ventana Ancho de Ventana Factor de Plenitud del Fierro Factor de Plenitud del Cobre Volumen del Fierro Volumen del Cobre Tensin de Fase Intensidad de Fase Frecuencia Densidad de Flujo Densidad de Corriente Seccin del Ncleo Seccin del Arrollamiento Potencia Solicitacin especfica Seccin del Arrollamiento Peso especfico del Fierro Peso especfico del Cobre Peso del Fierro Peso del Cobre Precio del Fierro Precio del Cobre Prdidas especficas en el Fierro Prdidas especficas en el Cobre m m m m3 m3 V A Hz T A / m2 m
2

Pag
Valorizacin de Prdidas Cobre Costo del Material Costo de las Prdidas Costo total Parmetro de relacin de costos Coeficiente de costo Variable auxiliar (dimetro)

2/7
$/W $ $ $ -

Valorizacin de Prdidas Fierro $ / W

Wc Km Kw Kt T Kp X Y Z

Variable auxiliar (ancho ventana) Variable auxiliar (costo) -

Variable auxiliar = X / Y Variable auxiliar = = X 4 Ff2 / T

3.

Pesos

m2 VA VA / m m2 kg / m3 kg / m3 kg kg $ / kg $ / kg W / kg W / kg 3.2
4

Se utiliza las dimensiones fundamentales: Dimetro de columna, Ancho y Altura de ventana, como se muestra en la Fig1
3.1 Peso del Ncleo La Seccin del Ncleo es: Sf = / 4 D2 Ff (1)

donde Ff es el factor de plenitud del fierro Ff = Seccin efectiva/ Seccin bruta Este factor de plenitud del fierro flucta entre 0.8 y 0.9, dependiendo del nmero de escalones de la seccin inscrita en el crculo (9.2) Peso de Ncleo: Gf = Sf( 3H + 4A + 6D ) f (2) donde la expresin :3H + 4A +6C es una aproximacin suficiente para la longitud media del ncleo(9.2). Peso del Arrollamiento:

________________________________________ ODT02 / Manuel Carranza A. / 03 Mayo 2002

IEEE INTERCON 2002 Optimizacin de Diseo de Transformadores

Para el clculo del peso por columna la seccin del cobre es: Sc
=

Pag
E = ( 32 / 4)FJB

3/7
(7)

A H Fc

(3)

La Potencia queda expresada como el producto de la Solicitacin Especfica, por la Seccin del ncleo, por la seccin del Arrollamiento: P = E Sf Sc (8)

Donde Fc es el factor de plenitud del cobre Fc = Seccin efectiva/ Seccin bruta Este factor de plenitud del cobre flucta entre 0.2 y 0.6, dependiendo de la tensin y tipo del arrollamiento (9.2) Segn el teorema del Guldin para cuerpos de revolucin, el volumen de cobre, por columna es: V c = Sc ( D + A / 2 ) (4) 5 El peso de cobre para las tres columnas es: Gc = 3Sc ( D + A / 2 ) c 4 Frmula de la Potencia Segn la ley de Faraday, para un flujo sinusoidal se tiene la siguiente tensin inducida: Tensin de fase U = 2 F N B Sf Haciendo uso de la densidad de corriente como parmetro y tomando en cuenta que a cada columna le corresponde la mitad de la ventana y que para la definicin de la potencia del transformador se toma slo uno de los dos arrollamientos (primario o secundario) se tiene: Intensidad de fase: I = J Sc / N Para la potencia de las tres fases se tiene: Potencia P=3UI (6) (5)

La potencia del transformador queda expresada por el producto de dos superficies: Seccin del ncleo y seccin del arrollamiento El coeficiente numrico correspondiente representa las densidades de Flujo magntico y corriente elctrica (solicitaciones). Procedimiento Matemtico Para el costo de la parte activa considera los costos de fierro cobre como costo por kg de material trabajado, a nivel de precio de venta del producto final de modo que: El costo de la Parte Activa es: Km = C f Gf + C c Gc (9)

Para considerar tambin las perdidas en el proceso de optimizacin, se aplica primero el valor de las prdidas especficas ( W/kg) que dependen de las densidades de flujo y de corriente, para el fierro y el cobre respectivamente. Despus se multiplica las prdidas por sus Coeficientes de Capitalizacin para poner sus costos a nivel de costo de adquisicin (9.5). Se tiene as: Costo de las Prdidas: Kw= Wf qf Gf + Wc qc G c (10)

Resulta el costo de material y prdidas: Costo Total: Kt = (Cf + Wf qf )Gf + (Cc+ Wc qc) G c (11) Para iniciar el anlisis matemtico reemplazamos en las frmulas de los pesos,

P = (32 /4)FJBSf Sc

Se define como parmetro de Solicitacin Especfica : ________________________________________ ODT02 / Manuel Carranza A. / 03 Mayo 2002

IEEE INTERCON 2002 Optimizacin de Diseo de Transformadores

las secciones de fierro y cobre en funcin de las variables fundamentales A,D, y H: El peso de fierro, ecc. (2) se convierte en: Gf = Sf ( 3H + 4A + 6D ) f Gf =(/4)D2 Ff ( 3H + 4A + 6D)) f (12) El peso de Cobre, ecc. (5) se convierte en: Gc = 3 Sc ( D + A / 2 ) c Gc = 3( H A /2 )Fc ( D + A / 2 ) c ( 13) 5.1 Reduccin de variables Independientes La expresin del costo (9): Km = C f Gf + C c Gc 3P f Cf Km = --------------- + AD2Ff fCf + EAFc 3 + ---------- D3 Ff f Cf + 2 6Pc Cc 3PAc Cc + ------------- + -----------------E D2 Ff E D Ff (18) 5.2 Cambio de Variables Se define como mdulo de longitud la expresin: M = ( P / E ) 1/4 (15) (19) Gc = Arrollamiento (13): Gc = 3(H A /2 )Fc ( D + A / 2 ) c 6P 3PA ---------- + -------------ED2 Ff EDFf (17)

Pag

4/7

Reemplazando en la ecc.(8) las secciones de fierro y cobre , se tiene:

P = E Sf Sc P= E(/ 4)D2 Ff AH Fc (14)

De aqu se puede despejar H, para continuar el anlisis con slo dos variables. P H=--------------------------E( / 4) D2 Ff A Fc

Con este valor de H los pesos se convierten en: Ncleo (12):

Gf =(/4)D Ff (3H + 4A+ 6D ) ) f 3Pf Gf = ---------- + AD2Ff f + EAFc 3f + ---------- D3 Ff f 2 (16) ________________________________________ ODT02 / Manuel Carranza A. / 03 Mayo 2002
2

Adems se define las variables auxiliares X e Y, tales que: D = X M A = Y M (20) (21)

Los trminos de la funcin de costo (18) se transforman, respectivamente en: 3 f Cf 3P f Cf ----------------- = ----- ----------- M 3 EA Fc Y Fc (22)

IEEE INTERCON 2002 Optimizacin de Diseo de Transformadores

AD2 fCfFf = X2 Y f CcFf M 3 (23 3 3 ----- D 3 fCfFf = ----- X 3 f CcFf M 3 2 2 (24) 6 c Cc 6P c Cc --------------- = ------ ---------- M 3 E D Ff X Ff (25) 3Y cCc 3PAcCc ------------------ = ------ ---------- M 3 X2 Ff E D2 Ff

Pag

5/7

Definiendo un nuevo parmetro: T , como Relacin de Costo Volumtricos: T = c Cc -----------f Cf

(29)

la ecuacin (28) se convierte en : 3 Km =M3fCf ------- + X2YFf + YFc 3YT 3 6T c + ------X3 Ff + --------- + ---------- 2 XFf X2 Ff (30)

(26)

Se tiene as, para la expresin del costo: 3f Cf Km=M3------------ + X2YfCfFf + Y Fc 3 + ---- --X3 f Cf Ff + 2 6cCc 3YcCc + ----------- + -- ------------- X Ff X2 Ff (27) o tambin: 3 Km =M3fCf ------ + X2YFf + YFc 3 + --------X3 Ff + 2 3YcCc 6cCc + --------------- + ----------------- XFffCf X2 FffCf (28) ________________________________________ ODT02 / Manuel Carranza A. / 03 Mayo 2002

5.2

Optimizacin del costo Para encontrar el costo mnimo, la funcin objetivo es: 3 3 Z =------ + X2YFf + -----X3 Ff + YFc 2 6T + -- ------- + XFf 3Y T -----------X2 Ff

(31)

Las condiciones para el valor mnimo se obtienen de las derivadas parciales: dZ 9 ----- = 2XYFf + ------- X2Ff dX 2 6T 6YT - -------- - -----------X2 Ff X3 Ff

=0 (32)

dZ

3T

IEEE INTERCON 2002 Optimizacin de Diseo de Transformadores

----- = - --------- + X Ff dY Y2 Fc
2

Pag

6/7

+ ------- = 0 X2 Ff (33)

dZ 3 3T ----- = - ------- + X2Ff + ------- = 0 dY Y2 Fc X2 Ff

Procedimiento numrico De la ecuacin. (32) se tiene, multiplicando por 2X3 Ff : 4X 4YFf 2 + 9X 5Ff 2 - 12X T - 12 YT
4 2

3 X2 Ff - ------------- + X 4 Ff 2 + 3T Y2 Fc

=0

De donde se obtiene una nueva expresin para : Fc( 3T + X 4 Ff2 ) 1 / 2 X ---- = = ------------------------------ Y 3 Ff = 0 (38)

= 0

X ( 9 X Ff

- 12T )

- Y ( 12T - 4X 4Ff 2 )

12T - 4 X 4 Ff 2 X ---- = ------------------------------ = 9X4 Ff 2 - 12T Y (34) El rango de aplicacin de X , se encuentra de las condiciones para hacer nulo el numerador, respectivamente el denominador de esta expresin. Para = 0 : X4 = ( 3 T ) / ( Ff 2 ) (35) Para = : X4 = ( 4T ) / ( 3 Ff 2 ) (36) Haciendo uso de como dato, en la ecuacin (34), se obtiene para X 4 : X 4 =12T(1 + )/( ( 4 + 9 ) ( Ff 2 ) ) (37) De la ecuacin. (33) se tiene, multiplicando por X 2 Ff :

6.1

Clculo numrico El problema numrico consiste en determinar X e Y partiendo de un valor dado para T, como variable independiente. Un procedimiento convergente aplicable para el clculo, consiste en: A) B) C) D) E) Asumir para X4 su valor mximo segn (33) 35 Calcular con la frmula (36) 38 Calcular X4 con la frmula (37 ) Calcular la diferencia entre los dos ltimos valores de X4 Si esta diferencia no es aceptable, volver a al paso B)

6.2

Convergencia del mtodo Si para la ecc. (34) se toma como variable auxiliar: = X 4 Ff2 / T (39)

________________________________________ ODT02 / Manuel Carranza A. / 03 Mayo 2002

IEEE INTERCON 2002 Optimizacin de Diseo de Transformadores

La ecc(34) se convierte en: 8 = ( 3 - ) / ( 2.25 * - 3 ) (40) La ecc(38) se convierte en: = (Fc *T* ( 1 + / 3 ) / Ff ) 1 / 2 (41) La ecc (40) representa una curva hiperblica donde el rango positivo de se da para valores de entre 1.333 y 3 La Ecc(41) representa una parbola donde el rango positivo de se da para valores de entre 3 y 3. La Solucin buscada est en la interseccin de estas curvas, como se muestra en la Fig 2 En la fig 3 se muestra como, el procedimiento empleado converge rpidamente hacia la solucin. 7 Programa de aplicacin El procedimiento expuesto ha sido utilizado en la elaboracin de un programa de clculo ejecutable en PC, denominado: TRAFOTEC. Se adjunta copia de la pantalla de trabajo. Con este programa se ha verificado tambin la convergencia del procedimiento empleado. 8.1 Conclusiones Este procedimiento de anlisis conduce a una ubicacin directa en la problemtica econmica del diseo, incluyendo los costos de las prdidas, durante el servicio. Este mtodo permite obtener resultados con a datos actuales, reemplazando as con ventaja la informacin referencial de los textos clsicos, contenida en tablas y curvas empricas, correspondientes a tecnologas y costos de la poca en que fueron editados. . En el grfico de sensibilidad se aprecia la importancia de ubicarse en el dimetro adecuado, quedando cierta holgura para elegir el ancho y altura de ventana, como sucede generalmente en la prctica para acomodarse a los requerimientos de disipacin de calor y tensin de corto circuito. El procedimiento numrico presentado puede servir de inspiracin para la solucin de otros casos de similar complejidad. Bibliografa Richter, Rudolf Die Transformatoren Verlag Birkhaeuser, Basel 1954 Sshaefer, Wilhelm: Transformatoren : Sammlung Goeschen Band 952. Berlin 1957 Corrales Martin, Juan: Teora, Clculo y Construccin de Transformadores Editorial Labor S.A. , Mjico 1960 Feinberg, Raphael :Modern Power Transformer Practice The McMillan Press Ltd. London 1980 Carranza, Manuel: Evaluacin de Prdidas de Transformadores Pag. Web de ABB Per. Lima 2001 : http:// www.abb.com/pe.

Pag

7/7

8.2

8.3.1

8.4

9 9.1 9.2 9.3

7.1

Anlisis de Sensibilidad En la figura 4 se muestra el efecto de las desviaciones de las dimensiones del ncleo, para el caso de tener: T = 1 ; Ff = 0.9 ; Fc = 0.6$ Se puede apreciar la importancia de ubicarse en el dimetro adecuado, quedando cierta holgura para elegir el ancho y altura de ventana.

9.4

9.5

________________________________________ ODT02 / Manuel Carranza A. / 03 Mayo 2002