Trabajo Teórico Experimental
Recepción: 01 Abril 2018
Aprobación: 01 Septiembre 2018
ABSTRACT:
Para estudiar los fenómenos que se originan durante la operación de los transformadores, empleando las herramientas de simulación existentes con resultados de una exactitud aceptable, es necesaria la obtención de datos reales del transformador y su aplicación en un modelo adecuado para cada fenómeno.Entre estos datos, tienen gran relevancia las reactancias del transformador.
En este artículo se presenta y evalúan los métodos de cálculo de la reactancia de dispersión en transformadores de devanados helicoidales concéntricos más generalizados en la literatura especializada y se propone la adecuación del método descrito por Kulkarni-Khaparde para su aplicación a devanados rectangulares. Se realiza el cálculo de la reactancia por los métodos descritos y por la variante propuesta, a un transformador de distribución, comparando los resultados teóricos con los valores obtenidos experimentalmente sobre el mismo transformador, y determinándose el error relativo que arroja cada procedimiento con relación al valor numérico de la medición realizada.
Palabras clave: inductancia, reactancia de dispersión, transformadores.
RESUMEN : In order to study the phenomena that originate during the operation of the transformers, using the existing tools of simulation with results of an acceptable accuracy, is necessary to obtain real data of the transformers and apply these in a suitable model for each phenomenon.Between these data, the reactances of transformers have great relevance.This paper presents and evaluates the two methods more generalized in the specialized literature of calculation of leakage reactance in transformers of helicoidal concentric windings and is proposed the adjusting of a method described by Kulkarni-Khaparde for its application in rectangular windings. It begins with a brief description of these methods and a discussion of them .The calculation of reactance by the described methods and the proposed variant is developed to a distribution transformer, comparing the theoretical results with the values obtained experimentally on the same transformer, and determining the relative error that yields each procedure with respect to the numerical value of the realized measurements.
Key words: inductance, leakage reactance, transformers.
Introducción
Desde hace algunos años, la modelación y simulación de fenómenos transitorios en los transformadores han ganado espacio en el campo de la ingeniería. La exactitud de los resultados obtenidos mediante un proceso de simulación depende del grado de detalle que se tenga en el modelo empleado y de la exactitud de los datos aportados, donde el valor de la reactancia de dispersión de los transformadores, tiene gran relevancia.
La reactancia de dispersión del transformador puede determinarse experimentalmente, mediante mediciones cuando existen las condiciones y es posible realizar el ensayo. También puede obtenerse de forma analítica con resultados razonablemente buenos. Algunos autores han documentado métodos teóricos de obtención de las reactancias de dispersión estableciendo una altura equivalente para cada bobina del transformador, como el resultado de la suma de la longitud axial de cada bobina más dos veces su espesor [1- 3]. Las expresiones de reactancias de dispersión que aportan estas publicaciones no pueden ser aplicadas directamente a transformadores con devanados de sección rectángular, dado que los autores realizan sus desarrollos sobre transformadores de devanados de secciones circulares.
El método descrito por Kulkarni-Khaparde realiza su desarrollo también sobre transformadores de devanados de secciones circulares, pero establece una altura equivalente común para todas las bobinas del devanado como el cociente de la altura real del devanado entre un factor menor que la unidad llamado Factor de Rogowski (KR) [4; 5]. Las expresiones de reactancias de dispersión que aporta tampoco pueden ser aplicadas directamente a transformadores con devanados de sección rectángular.
Ambos autores plantean expresiones de aplicación en transformador de devanados helicoidales cuyos perímetros estén inscritos en un círculo, sin embargo, la gran mayoría de los transformadores de mediana potencia y de distribución tienen devanados cuyos perímetros están inscritos en un rectángulo, por lo que es necesario realizar modificaciones a estos métodos para poder aplicarlos a estos últimos.
Materiales y métodos
En el campo magnético de dispersión de un transformador tiene notables diferencias entre la distribución del flujo de dispersión que existe en la dirección axial de la ventana del núcleo y la que existe en la dirección perpendicular a esta. Cuando las bobinas de alta tensión (AT) y baja tensión (BT) tienen la misma altura, los amper-vueltas, se encuentran distribuidos uniformemente en la altura de los devanados, entonces el flujo de dispersión es predominantemente axial, excepto en los extremos, donde se deforma buscando el camino más corto para cerrar a través del yugo o la columna del núcleo [1 - 5]. Ver figuras 1a y 1b.

El modelo del flujo de dispersión típico que se muestra en la figura 1a, puede ser reemplazado por líneas de flujo paralelas, todas de una longitud igual a una altura equivalente como se muestra en la figura 2a.

En la figura:
A.T.Devanado de alta tensiónB.T.Devanado de baja tensiónδ1Espacio que ocupa radialmente el conductor en A.Tδ2Espacio que ocupa radialmente el conductor en B.TδoEspacio entre el conductor de A.T. y B.T. También llamado entre cara activaLILado más interno del devanadoLELado más externo del devanadohdAltura real del devanadoAvAltura de la ventana del núcleoMétodo de Juan Corrales Martín [1].
Plantea que las alturas equivalentes primaria y secundaria para las trayectorias del flujo de dispersión son las siguientes. Vea ecuaciones (1-6).
Donde:
heq1Altura equivalente del devanado primarioheq2Altura equivalente del devanado secundariohd1Altura real del devanado primariohd2Altura real del devanado secundarioX1Reactancia de dispersión del devanado primarioX2Reactancia de dispersión del devanado secundarioFFrecuencia del sistema eléctricoμoPermeabilidad del vacío; 4π . 10-9 H/cmN1Número de vueltas del devanado primarioN2Número de vueltas del devanado secundarioLm1Perímetro medio del devanado primarioLm2Perímetro medio del devanado secundarioLmoPerímetro medio de la entre cara activaMétodo de Kulkarni-Khaparde [4, 5].
Plantea que una altura equivalente para todos los devanados como se muestra en las ecuaciones (7) y (8).
Donde:
KRFactor de Rogowski.
El factor de corrección de Rogowski es utilizado para tener en cuenta el flujo de dispersión en la parte superior e inferior de los devanados y el efecto del núcleo magnético [6;7].
Los autores consideran un trayecto cerrado de flujo a una distancia x desde el lado interior (LI) del devanado, estableciendo la densidad de flujo como se muestra en la ecuación (9).
Para llegar a la ecuación de la reactancia, se obtiene una expresión general de las concatenaciones de flujo de un tubo de flujo con una altura equivalente heq y una pared de grosor δ, como se aprecia en la figura 3.

En la figura 3; NI es el valor de los amper-vueltas nominales del devanado, a(NI) y b(NI) son las fracciones de amper-vueltas correspondientes a LI y LE respectivamente.
Para un punto a una distancia x desde LI la f.m.m. se puede plantear según la ecuación (10).
Sustituyendo la ecuación (9) en la ecuación (8) y sacando factor común NI; el valor eficaz de la densidad de flujo a una distancia x de LI queda según la ecuación (11).
Las concatenaciones de flujo de un tubo de flujo incremental de ancho de x emplazado en x, son las expresadas en la ecuación (12).
Donde:
ϕxFlujo que concatena con las vueltas NxNxCantidad de vueltas en el punto x. Se puede obtener de acuerdo con la ecuación (13)
Si se tiene en cuenta lo que se expresa en la ecuación (14):
Donde:
AÁrea de la pared del tubo de flujo, m2Entonces la ecuación (11), se puede transformar en ecuación (15).
Donde dA es el área incremental del tubo de flujo.
Sustituyendo las ecuaciones (11) y (13) en (15), quedan las ecuaciones (16) y (17).
En devanados circulares concéntricos, como se muestra en la figura 4, el área incremental del tubo de flujo se obtiene según la ecuación (18).
Sustituyendo la ecuación (18) en (17) queda la ecuación (19).

En la figura:
D1Diámetro medio del devanado primarioD2Diámetro medio del devanado secundarioDmDiámetro medio de todo el devanado o la entre cara activaPor lo tanto, el flujo total del tubo se obtiene de la siguiente integral. Vea ecuación (20).
Integrando y resolviendo, la ecuación (20) queda de la siguiente manera. Vea ecuación (21)
El último término dentro de los corchetes puede ser despreciado sin que esto arroje un error apreciable [4, 5]. Vea ecuación (22).
El término [LI + 3δ/2] puede ser considerado igual al diámetro medio (Dm) del tubo de flujo (esto es admisible, principalmente, para arrollamientos que tienen valores de diámetro comparativamente más grande que los valores de sus espesores radiales) [4; 5]. Vea ecuación (23).
Se considera. Vea ecuación (24).
Que corresponde al área del Diagrama de Amper - Vueltas para el caso del devanado de forma circular. La inductancia de dispersión de un transformador con n tubos circulares de flujo se puede obtener como se muestra en la ecuación (25). [4; 5]:
Y la correspondiente expresión para la reactancia de cortociruito Xcc es como se muestra en la ecuación (26). [8-10]:
Las constantes a y b de la ecuación (24) tienen los valores de 0 y 1 para BT, 1 y 1 para la entre cara activa y 1 y 0 para AT. Entonces, utilizando la ecuación (26) se obtiene el área del Diagrama de Amper - Vueltas para devanados cilíndricos como se muestra en la ecuación (27).
Modificación del método de Kulkarni-Khaparde para devanados helicoidales con forma rectangular.
Cuando el núcleo y el devanado tienen sección rectangular el Diagrama de Amper - Vueltas no es el mismo dado que el área incremental del tubo de flujo y los perímetros de los devanados tienen otra forma de obtención.
Haciendo el análisis para un trayecto cerrado de flujo a una distancia x desde el lado interior de un devanado de forma rectángular como se muestra en la figura 5, se tiene:

En la figura:
AfeLongitud del lado menor de la columna del núcleo magnético rectangularProLongitud del lado mayor de la columna del núcleo magnético rectangularLI1Longitud del lado interior menor del devanadoLI2Longitud del lado interior mayor del devanadoLmLongitud media de la espira del devanado a una distancia xHolaHolgura entre el núcleo y el devanado por su lado de mayor longitudHolpHolgura entre el núcleo y el devanado por su lado de menor longitudEn este caso, el área incremental del tubo rectangular de flujo se obtiene según la ecuación (28).
Donde:
Dado que (2LI1+2LI2) es el perímetro interno del tubo rectangular, la ecuación (28) se escribe como se muestra en la ecuación (29)
Donde:
Pintperímetro interno del tubo rectangularSustituyendo la ecuación (29) en la ecuación (17), queda la ecuación (30).
Por lo tanto, el total de las concatenaciones de flujo de un tubo de forma rectangular está dado por la ecuación (31).
Resolviendo la integral, aplicando la regla de integración por partes, queda de la siguiente manera. Vea ecuación (32).
Sumando y restando a2 en el último término dentro del paréntesis de la ecuación (32), y sacando factor común, queda de la siguiente manera. Vea ecuación (33).
Sustituyendo la ecuación (33), en la ecuación (31), queda la ecuación (34):
El Diagrama de Amper - Vueltas para el caso de devanado de forma rectangular, sería como se muestra en la ecuación (35).
El ATDR para el devanado primario, se tiene la ecuación (36).
El ATDR para la entre cara activa del devanado, es la ecuación (37).
Y el ATDR para el devanado secundario, se tiene la ecuación (38).
El ATDR para el cálculo de la inductancia tanto del primario como secundario, es la suma del área del devanado primario ecuación (36), y/o devanado secundario ecuación (38), más la mitad del área de la entre cara activa ecuación (37), por lo tanto las reactancias de dispersión primaria y secundaria tendrán la siguiente forma. Vea las ecuaciones (39) y (40).
De la reactancia de dispersión del devanado primario y secundario, se puede calcular la reactancia de cortocircuito, por la ecuación (41).
Resultados y discusión
Cálculo de la reactancia por los métodos descritos
Se muestran los resultados del cálculo de la reactancia de dispersión de un transformador de distribución de dos devanados cuya reactancia de cortocircuito medida según [11] es 39,618 Ω. Los parámetros físicos del transformador, se muestra en la tabla 1.

El valor de heq se calcula por medio de la ecuación (6), en la cual el factor de Rogowski está dado por la ecuación (7).
KR=0,8819
heq=17,9059cm
Cálculo de las reactancias por el método de Kulkarni-Khaparde
Para la aplicación del método, dado que la sección de la columna del núcleo del transformador en cuestión es rectangular, el diámetro interior (LI) del devanado de BT, se tiene que obtener como el valor promedio entre los valores LI1 y LI2 del tubo rectangular de flujo. Estos son la suma de los lados de la columna del núcleo más las holguras, ver figura 5 y ecuación (42):
Los diámetros medios de los devanados de BT y AT, y del espacio entre BT y AT se calcula a partir de las ecuaciones (43), (44) y (45).
El valor de se calcula a partir de la ecuación (26):
Aplicando la ecuación (25), para calcular la reactancia de cortocircuito para devanados de forma circular, se tiene:
Xcc = 34,055 [Ω]
El error relativo que se obtiene con la aplicación de este método se obtiene por la ecuación (46):
Cálculo de las reactancias con la modificación del método de Kulkarni-Khaparde
Para la aplicación del método de Kulkarni-Khaparde modificado para devanados rectangulares, las reactancias de dispersión y de cortocircuito se calculan mediante las ecuaciones (39), (40) y (41):
El error relativo que se obtiene con la aplicación de este método se obtiene por la ecuación (46):
Conclusiones
Referencias
[1] CORRALES MARTÍN, Juan. Cálculo industrial de máquinas eléctricas. ISBN: 8426704387. Idioma: SPA. P. imprenta: Barcelona, Marcombo, 1982. 2 v. : 23 cm. Edición ; 1. ed.
[2] DEL VECCHIO, Robert M., et al. Transformer Design Principles With Applications 3e: With Applications to Core-Form Power Transformers. CRC press, 2017. ISBN 9781498787536 - CAT# K30258. Edición; 3. ed
[3] MCLYMAN, Colonel Wm T. Transformer and inductor design handbook. CRC press, 2011. ISBN 9781439836873 - CAT# K11667. Edición; 4. Ed
[4] KULKARNI, Shrikrishna V.; KHAPARDE, S. A. Transformer engineering design and practice. Editorial Marcel Deckkerinc. .New York. EE.UU, 2004. ISBN 0824756533, 9780824756536.
[5] KULKARNI, Satish V.; KHAPARDE, S. A. Transformer engineering: design, technology, and diagnostics. CRC Press, 2012. ISBN 9781439853771 - CAT# K12483. Edición; 2. Ed
[6] ARIVANMUDHAN, M.; SANTHI, S. et al. FEM based reactance computation for transformers winding deformation detection. International Journal of Computer Applications (0975-8887). National Conference Potential Research Avenues and future Opportunities in Electrical and Instrumentation Engineering . 2014. [Consultado octubre 2017]. Disponible en: Disponible en: https://www.ijcaonline.org/proceedings/eteiac/number1/17326-1401
[7] MEMIJE, D. Diseño del núcleo y bobinas para un transformador por medio de un simulador. Tesis de grado de ingeniería eléctrica. Instituto Politécnico Nacional. México. 2014.[Consultado junio 2017]. Disponible en: Disponible en: https://es.scribd.com/document/342677611/DISENO-DEL-NUCLEO-Y-BOBINAS-PARA-UN-TRANSFORMADOR-POR-MEDIO-DE-UN-SIMULADOR-pdf
[8] OLIVEIRA, Luís MR; CARDOSO, AJ Marques. Modeling and characterization of leakage inductances for transformer winding fault studies. 2016. p.423-430. Disponible en: https://link.springer.com/chapter/10.1007/978-3-642-37291-9_45
[9] AVEDAÑO, A.; MORK,, V.A.; HOIDALEN, H.K.. Transformer internal fault modeling in ATP. Paper submitted to the International Conference on Power Systems Transient (IPST2011) in Delph. Netherlands. 2011. [Consultado octubre 2017]. Disponible en: Disponible en: http://www.ipstconf.org/papers/Proc_IPST2011/11IPST017.pdf
[10] CHACMANA, S. Diseño de un transformador multitap para pozos petroleros. Tesis de grado de ingeniería eléctrica. Universidad Nacional de Ingeniería. Lima. Perú. 2013.[Consultado septiembre 2017]. Disponible en: Disponible en: http://cybertesis.uni.edu.pe/handle/uni/8435
[11] RAMOS GUARDARRAMA, et al. Ensayos de pérdidas en vacío y con carga en transformadores mediante la adquisición de datos. Ingeniería Energética, 2016, vol. 37, n 1, p. 73-80. ISSN 1815-5901.[Consultado octubre 2017]. Disponible en: Disponible en: http://scielo.sld.cu/scielo.php?script=sci_arttext&pid=S1815-59012016000100009
Notas de autor
*Autor para correspondencia: Cristian Ivan Torres Suquilanda. E-mail:cristivantorres@gmail.com