Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
1
UNIVERSIDAD NACIONAL DE INGENIERIA FACULTAD DE INGENIERIA ELÉCTRICA Y ELECTRÓNICA INSTITUTO DE INVESTIGACIÓN
PROYECTO DE INVESTIGACIÓN
“ESTUDIO DE LAS MÁQUINAS ELÉCTRICAS, SIMULACIÓN DIGITAL”
JEFE DEL PROYECTO: MSc. Ing. FÈLIX VÌCTOR CÁCERES CÁRDENAS INTEGRANTES: Ing. MODESTO TOMÁS PALMA GARCÌA Sr. DAVID FLORES RODRIGUEZ Sr. ERIC WASHINGTON TORRES PADILLA LIMA – PERU 2002 “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
2 INDICE GENERAL CAPITULO 1. CIRCUITOS MAGNETICOS 1.1 Propiedades magnéticas de la materia 1.2 Características de magnetización 1.2.1 Representación matemática 1.3 Circuito análogo resistivo 1.4 Ecuaciones de Maxwell 1.4.1 Ley circuital de ampere 1.4.2 Configuración de laminas 1.5 Excitación de circuitos magnéticos 1.5.1 Excitación con corriente continua 1.5.2 Excitación con corriente alterna 1.6 Reactor con núcleo de hierro CAPITULO 2. TRANSFORMADOR NÚCLEO DE HIERRO 2.1 Introducción 2.2 Transformador ideal 2.2.1 Polaridad 2.2.2 Impedancias y variables reflejadas 2.3 Transformador núcleo de hierro 2.4 Determinación de parámetros ( condiciones nominales ) 2.4.1Prueba de cortocircuito 2.4.2 Prueba de vacío 2.4.3 Prueba bajo carga 2.4.3 Pruebas especiales 2.5 Eficiencia 2.5.1 Eficiencia durante todo el día 2.6 Regulación 2.6.1 Regulación de tensión bajo carga 2.7 Calentamiento en las máquinas eléctricas 2.7.1 Balance Termodinámico 2.8 Balance por unidad 2.9 Autotransformadores 2.9.1 Circuito equivalente 2.9.2 Ventajas y desventajas CAPITULO 3. TRANSFORMADORES EN SISTEMAS TRIFÁSICOS 3.1 Introducción. 3.2 Especificaciones técnicas. 3.3 Conexionado de transformadores trifásicos. 3.4 Tipos y grupos de conexión. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
3 3.4.1Tipo de conexión. 3.4.2 Grupos de conexión. 3.5 Banco de transformadores trifásicos. 3.6 Pruebas a banco de transformadores trifásicos. 3.6.1 Prueba de cortocircuito. 3.6.2 Prueba de vacío 3.7 Puesta en paralelo de transformadores. 3.8 Reparto de carga. 3.9 Sistemas eléctricos desbalanceados 3.10 Transformadores de tres arrollamientos. 3.10.1 Circuito equivalente. CAPITULO 4. CONVERSIÓN DE ENERGÍA ELECTROMECÁNICA 4.1 Introducción. 4.2 Ley Universal de Faraday. 4.3 Sistemas Eléctricos y mecánicos. Variables y funciones de estado. 4.3.1 Análisis para un elemento inductivo. 4.3.2 Análisis para un elemento inductivo. 4.4 Forma restringida de la ecuación de Lagrange. 4.5 Determinación de la fuerza mecánica a través del Principio de los trabajos virtuales. 4.6 Sistemas de excitaciones. 4.6.1 Excitación simple. 4.6.2 Excitación doble. CAPITULO 5. CONCEPTOS BÁSICOS PARA EL ESTUDIO DE LAS MÀQUINAS ELÉCTRICAS Y TRANSFORMACIONES. 5.1 Campo magnético generado por devanados estatóricos. 5.2 Campo magnético generado por devanados rotóricos. 5.3 La máquina primitiva. 5.4 Parámetros de inducción estacionaria de la máquina d-q. 5.5 Tensiones inducidas por rotación. 5.6 Ecuaciones de equilibrio de la máquina d-q. 5.7 Vectores espaciales. 5.8 Ecuación de equilibrio mecánico. CAPITULO 6. LA MÀQUINA DE INDUCCIÓN. 6.1 Planteamiento matemático de las ecuaciones eléctricas. 6.1.1 Determinación del flujo común de los devanados estatóricos y rotóricos. 6.1.2 Determinación del flujo de dispersión. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
4 6.1.3 Consideración conjunta del flujo común y del flujo de dispersión. 6.2 Interpretación física del concepto de vector espacial. 6.3 Segundo Método para la determinación de las ecuaciones eléctricas de la máquina. 6.4 Significado físico de las magnitudes Ls, Lr y Lm. 6.5 Deducción de la ecuación mecánica. Ecuación del torque. 6.5.1 Primer Método: Interacción campo magnético – corriente. 6.5.2 Segundo Método: Energía de campo. 6.6 Sistema general de ecuaciones electromecánicas. Solución del Sistema para el caso de régimen estático. 6.7 Sistema general de ecuaciones del régimen dinámico incluyendo el vector espacial de corriente de magnetización. 6.8 Sistema general de ecuaciones eléctricas considerando el efecto de la saturación magnética. 6.8.1 El efecto de la saturación. CAPITULO 7. ANÁLISIS TRANSITORIO DE LA MÀQUINA SINCRONA. 7.1 Inductancias de la máquina síncrona de polos salientes. 7.2 La transformación de Park. 7.3 Cortocircuito trifásico del generador síncrono en vacío. 7.4 Simulación: cortocircuito trifásico de la máquina síncrona. 7.4.1 Programa en Matlab. 7.4.2 Resultados de la simulación. 7.4.3 Comentarios. 7.5 Modelos Simplificados de la máquina síncrona para l análisis de transitorios. 7.5.1 Determinación de las constantes de tiempos τ`` , τ` de la máquina síncrona. 7.6 Estabilidad dinámica de la máquina síncrona. 7.6.1 Método de igualdad de área. CAPITULO 8. PROBLEMAS Bibliografía
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
5
CAPITULO 1 CIRCUITOS MAGNETICOS 1.1 PROPIEDADES MAGNETICAS DE LA MATERIA Toda materia está constituida fundamentalmente de átomos, y cada átomo consiste en electrones en movimiento. Estos circuitos de electrones, cada uno de los cuales esta confinado a un solo átomo son los que se define como corrientes atómicas. Parecería que se tuviera dos clases de corriente: 1) Una corriente verdadera, que consiste en transporte de carga; esto es el movimiento de electrones libres o de iones cargados, y 2) Corrientes atómicas, que son corrientes puras que circulan sin dar origen a transporte de carga. Sin embargo ambas clases de corrientes pueden producir campos magnéticos. En las máquinas eléctricas tales como: transformadores, generadores y motores se utilizan una gran variedad de materiales magnéticos en cuanto a su tipo y formas. Todos los materiales ferromagnéticos utilizados en las máquinas eléctricas se caracterizan por poseer una permeabilidad relativa elevada ( µr ) y una relación no lineal entre la densidad magnética ( B ) y la intensidad magnética ( H ). PRINCIPALES SISTEMAS DE UNIDADES Tenemos los siguientes Sistemas de Unidades. 1. CGS ues, cuya cuarta magnitud es la permitividad (∈). 2. CGS uem, cuya cuarta magnitud es la permeabilidad ( µ ). 3. Sistemas mixtos de Guass (utilizado por los físicos). (coule ) 2 4. MKSA Giorgi no racionalizado µ0=10-7; ∈=1,11*10-10 N.m 2 5. MKSA Giorgi racionalizado (utilizado por los electricistas) µ0=4π∗10-7 Vs/Am ó H/m, ∈0=8.85∗10-12 As/Vm ó F/m. Los dos sistemas racionalizados y no racionalizado van a distinguirse únicamente por la relación del ángulo sólido. Sistema MKSA no racionalizado, unidad de ángulo sólido = estereorradián. Sistema MKSA racionalizado, unidad ángulo sólido espacio total spat. Spat = 4π estereorradián. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
6
Magnitud 1) Angulo Sólido.
Unidad Estereorradiá n
2) Trabajo Joule y energía.
3) Potencia
Símbolo
Descripción
Sr
Ángulo sólido que subtiende una superficie de cualquier esfera con centro en su vértice. Equivale a una superficie de 14 de la superficie total de la misma.
J
Trabajo producido por una fuerza de 1,0 Newton al trasladar su punto de aplicación en la misma dirección y sentido una distancia de 1,0 m. Trabajo realizado por un sistema de potencia constante, de un watt en un segundo. Energía adquirida por un electrón acelerado por una diferencia de potencial de 1,0 voltio. Potencia constante de un sistema que desarrolla uniformemente un trabajo de 1,0 Joule durante un intervalo de tiempo de 1,0 seg. Intensidad de una corriente eléctrica invariable que mantenida en dos conductores paralelos, rectilíneos de longitud infinita y sección despreciable situados a una distancia de 1,0 m uno del otro en el vacío origina una fuerza de 2 x 10-7N. Sobre cada trecho de 1,0 m de longitud de cada uno de los dos conductores. Carga eléctrica que atraviesa durante un segundo una sección transversal de cualquier conductor recorrido por una corriente eléctrica de intensidad constante igual a 1,0 Amperio. Diferencia de potencial eléctrico existente entre dos puntos de un conductor recorrido por una corriente eléctrica de intensidad constante de 1,0 Amperio cuando la potencial disipada entre esos dos puntos es de 1,0 Watt.
Watt-Seg.
Ws
Electrón Volt.
EV
Watt
W
4) Intensidad Ampere de Corriente Eléctrica
A
5) Carga Eléctrica
Coulomb
C
6) Tensión Eléctrica
Voltio
V
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
7 1.2 CARACTERISTICAS DE MAGNETIZACIÓN Cada material ferromagnético como por ejemplo: acero, fierro silicoso, grano orientado u otros; presentan su propia característica de magnetización, que son reflejados mediante la representación de curvas o de formulaciones matemáticas ( modelos). Los materiales ferromagnéticos se caracterizan por poseer uno ó varios de los atributos siguientes : a) Pueden magnetizarse mucho más fácilmente que los demás materiales ( p.e. material diamagnético ). Esta característica viene indicada por que los materiales ferromagnéticos tienen una alta permeabilidad relativa ( µr ). b) Se imanan con una facilidad muy diferente según sea el valor del campo magnético. Este atributo lleva a una relación no lineal entre los módulos de la densidad magnética ( B ) y la intensidad magnética ( H ). c) Un incremento del campo magnético les origina una variación del flujo diferente de la variación que origina una disminución igual de campo magnético. Este atributo indica que las relaciones que expresan la inducción magnética y la permeabilidad µ como funciones del campo magnético, no son lineales ni uniformes. d) Conservan la imanación cuando se suprime el campo magnético. e) Tienden a oponerse a la inversión del sentido de la imanación una vez imantados. f) Tienen una densidad magnética intrínseca máxima Bmáx elevada.
Fig.1 Curva de Característica de magnetización de un material ferromagnético y las zonas definidas
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
8 1.2.1 REPRESENTACION MATEMATICA Existe una variedad de modelos matemáticos que son formulados para poder representar la curva de la característica de magnetización de los materiales ferromagnéticos, siendo algunas de estos los siguientes: A. Análisis por Zonas, esta forma de representar la curva de la característica de magnetización de un material ferromagnético está relacionado con la definición de tres zonas típicas como son: Zona I lineal Zona II Codo de saturación Zona III Saturación
: : :
y = c0 + c1x y = c0 + c1x + c2x2 y = c0 + c1x + c2x2 + c3x3
En las relaciones funcionales anteriores, están representadas la densidad magnética y la intensidad magnética como las variables x e y respectivamente. Las constantes co, c1 y c2 son propias de cada material ferromagnético. B. Aproximación Lineal, este método consiste en linealizar la función de la curva de característica de magnetización, para ello se aplica en cada punto el concepto de la derivada. C. Formulación Polinomial, este método consiste en expresar una relación polinomial de las variables magnéticas B y H, donde se distingue la No linealidad entre ellas, siendo expresadas por la relación siguiente: ó
H = a0 + a1B + a2B3 + a3B5 + ... 3
5
i = b0 + b1λ+ b2λ + b3λ + ...
……….. (1.1) ………. (1.2)
Donde: H: Intensidad magnética B: Densidad magnética i : Corriente eléctrica λ: Acoplo inductivo a0, a1, ... an , b0, b1, ... bn : Constantes, obtenidas del ensayo de laboratorio. D. Ecuación de Froelich, esta relación funcional entre la densidad (B) e intensidad (H) magnética es conveniente en la zona del codo de saturación, está expresada por:
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
9 B=
a 0H a1H + a 2
V=
;
c 0i c 1i + c 2
........... ( 1.3 ) 1.3 CIRCUITO ANÁLOGO RESISTIVO Las máquinas eléctricas, como los transformadores, motores y generadores están constituidos por estructuras ferromagnéticas conformadas por un número de láminas en el primario(estator) y secundario(rotor); asimismo provistas por un conjunto de bobinas que contienen un número determinado de espiras. Para su análisis del comportamiento estacionario y transitorio se requiere realizar la formulación de las máquinas eléctricas en modelos circuitales, que se encuentran constituidas por elementos puros como son: inductancias, resistencias y capacitancias. Los circuitos análogos resistivos dependen de la configuración del circuito magnético, y para su formulación se requiere tomar en cuenta los conceptos siguientes: : A. Continuidad de las líneas de flujo, nos indica que, en una superficie cerrada el flujo magnético entrante es igual al flujo magnético saliente de la superficie. La expresión matemática es la siguiente: Σφ=0 ........... ( 1.4 ) B. La Ley circuital de Ampere, esta expresada por:
∫ H - dl = I
Ejemplo Ilustrativo:
............. (1.5 )
wx
A2
A1
φ1
A
B
φ2
C
φ I λ1 λ1
F
N A
E
λ 2 D
Fig. 1.2 Circuito Magnético de tres columnas “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
10 Rn =
lm um A m
lm: longitud media por donde recorre el flujo magnético
;
um: Permeabilidad magnética Rm: Reluctancia magnética Am: Área del material magnético por donde pasa el flujo
Rm = 1/Pm
;
Pm: Permeancia magnética Rm: Reluctancia magnética
R eq = R qm + R qg
φ n ( Rqm + Rqg ) = Ni
........... ( 1.6 )
Circuito eléctrico Análogo:
φ1
φ2 R
φ
φ
R1
R2 NI
R
R1/2
NI
Fig.1.3 Circuito Análogo- resistivo. Las reluctancias en el lado 1 y 2 son iguales, es decir: R1 = R2 Del circuito equivalente de la fig.1.3, se expresa la relación siguiente: λ λ R NI = φ R + 1 → NI = φ + 1 2 µA 2 µA1
NI =
; Si se tiene: A = 2 A1
φ [λ + λ1 ] → NI = BA [λ + λ1 ] → NI = B [λ + λ1 ] µA µA µ ........... (1.7 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
11 Si: NI es conocido Æ del grafico B-H, se calcula B y µ. B
µ
H Fig. 1.4 Curva de la Característica de magnetización NI = H (λ + λ1 ) → H =
NI (λ + λ1 )
........... (1.8 ) ELECTRICOS MAGNETICOS
Ve I Re J G
Vm φm Rm um P
CONSIDERACIONES:
a) Se considera que el flujo magnético es uniforme en una sección transversal del núcleo. b) Las longitudes del circuito magnético son consideradas como longitudes medias (lm). c) De la ecuación continuidad de las líneas de flujos, tenemos: ∑ φ i = 0 i n
d) La ecuación de la ley circuital de Ampere se considera: NI = ∑ Hi λi i
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
12 De la figura del ejemplo ilustrativo, podemos expresar: Para la trayectoria A-B-E-F-A, tenemos: λ1*H1+λ*H.=Ni Para la trayectoria A-C-D-F-A, tenemos: λ1H1+λ2H2=Ni donde:
...........( 1.9 ) .......... ( 1.10 )
Ni : FMM fuerza magnetomotriz Vmi = λiHi diferencia de potencial magnético en la rama i.
para la rama 1, tenemos: Vm1=λ1H1; Vm1=λ1B1/µ1; Vm1=λ1φ1/A1µ1
.…….( 1.11 )
λ λ Vm1 = 1 φ1 , se define: R1 = 1 Reluctancia de la rama 1. µ 1A 1 µ 1A 1 Vm1 = R1φ1
Similarmente, para las otras ramas: Vm 2 = R 2 φ 2
;
Vm3 = R 3 φ 3
........... (1.12 )
1.4 ECUACIONES DE MAXWELL
Es un conjunto de Ecuaciones que representan una generalización de algunas observaciones experimentales, dichas ecuaciones son:
rotH = J +
rotE = −
∂D ∂t
∂B ∂t
…........( 1.13 )
...........( 1.14 )
divD = P
............( 1.15 )
divB = 0
...........( 1.16 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
13 La ecuación ( 1.13 ), representa una extensión de la Ley de Ampere. La ecuación ( 1.14 ), es la forma diferencial de la Ley de Inducción Electromagnética de Faraday. La ecuación ( 1.15 ), es la Ley de Gauss que a su vez se deduce de la Ley de Coulomb. La ecuación ( 1.16 ), representa la continuidad de las líneas de flujos, el hecho de que los polos magnéticos simples no se han observado. Las ecuaciones de Maxwell representan expresiones matemáticas de algunos resultados experimentales, siendo evidente que no pueden demostrarse, sin embargo la aplicabilidad a cualquier situación puede ser verificada. 1.4.1 LEY CIRCUITAL DE AMPERE
Para el análisis en bajas frecuencias, de la ecuación ( 1.13 ) se tiene que D =0, luego: rotH = J además: B = µ0H; siendo µ0: Permeabilidad del vacío o aire. luego: rotB= µ0J
........... ( 1.17 )
Aplicando el Teorema de Stokes, tenemos:
∫ rotB οnds = ∫ B οdλ S
........... ( 1.18 )
C
Luego, utilizando la ecuación ( 1.17 ) para el rotacional de B, se tiene:
∫ B οdλ = µ ∫ J οnds C
0
............ (1.19 )
S
Para una espira, tenemos:
∫ B ο dλ = µ I
Para “N” espiras, tenemos:
∫ B ο dλ = µ NI
C
C
0
0
Luego, podemos expresar:
∫ H ο dλ = NI C
“Estudio de las Máquinas Eléctricas. Simulación Digital”
............ (1.20 )
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
14 1.4.2 CONFIGURACION DE LAMINAS
Los circuitos magnéticos están conformados por un conjunto de láminas, cuyos espesores pueden ser de 0,30 mm. ó 0,50 mm; asimismo presentan diversas configuraciones geométricas. La característica de que los circuitos magnéticos sean laminados es para reducir las pérdidas en el núcleo magnético. Las configuraciones de las láminas pueden ser de tipo C-I, E-I, E-E ú otras. Los materiales Magnéticos (ferromagnéticos) pueden ser :Acero-Silicoso, fierro fundido, fierro-silicoso, grano orientado y otros. 1.5
EXCITACION DE CIRCUITOS MAGNETICOS
1.5.1 EXCITACIÓN CON CORRIENTE CONTINUA
Los circuitos magnéticos pueden presentar diversas configuraciones geométricas y ser excitados con fuentes de corriente continua, que para conocer sus variables eléctricas y magnéticas se requiere la utilización de métodos de solución, como son los siguientes: MÈTODO DE SOLUCION
Se presentan dos métodos de solución, como son: el método gráfico y el método iterativo. A) METODO GRÀFICO
Para la solución de los circuitos magnéticos mediante el método gráfico, se deberán tomar en cuenta los procedimientos siguientes: 1. Tener en cuenta la curva de la característica de magnetización del material ferromagnético en cuestión, es decir la curva B –H, que en la practica es proporcionado por el fabricante ó se determina de los ensayos en laboratorio; por ejemplo en el caso de un Transformador monofásico de núcleo de hierro, dicha curva B –H es obtenida realizando la prueba de vació. 2. De acuerdo a la configuración geométrica del circuito magnético, se formula la ecuación de equilibrio electromagnético, obteniéndose la ecuación de una recta con pendiente negativa. 3. La curva característica B – H, se debe llevar a una linealización debido a que la ecuación obtenida en el punto 2 es lineal. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
15 4. Relacionando el punto 2 y punto 3, determinamos el punto de intercepto de ambas características, obtenemos el punto de operación requerido.
B
BfB
Bt
P H Ht
HfH
Fig. 1.5 Determinación del punto de operación Ejemplo ilustrativo: Dado el circuito magnético siguiente:
I
Lf
Lg N
φ Fig. 1.6 Circuito magnético con entrehierro Donde: λf λg Af Ag
: Longitud media del material ferromagnético. : Longitud del entrehierro. : Área del material : Área del entrehierro.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
16 Bf Bg Hf
: Densidad de flujo en el material (Wb/m2). : Densidad de flujo en el entrehierro. : Intensidad de campo magnético en el material (Amp/m).
Aplicando la Ley Circuital de Ampere, tenemos: H f λ f + Hg λg = NI Hg =
.........( 1.21 )
Bg
.........( 1.22 )
µ0
por continuidad de flujo: φf = φg Æ BfAf = BgAg Æ B g =
Af Bf Ag
en la ecuación(1), tenemos: Hf λf +
A f λg µ0 A g
B f = NI
µ0 A g µ A λ NI − 0 g f Bf = A λ A λ f g f g
........... ( 1.23 ) H f
...........( 1.24 )
La ecuación ( 1.24 ) es una recta con pendiente negativa. De la ecuación ( 1.24 ), tenemos: Si: Bf = 0 Æ H fH =
NI ( intercepto con H ) λf
µ0 A g NI ( intercepto con B ) Si: Hf = 0 Æ B fB = A λ f g
B) MÉTODO ITERATIVO “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
17
El método Iterativo consiste en realizar iteraciones sucesivas hasta obtener un % Er ≤ 5% Donde: Er =
NI( dato ) − NI( calculado ) x100 NI( dato )
...........( 1.25 )
Para la solución de los circuitos magnéticos por el método Iterativo, se deberá considerar el procedimiento siguiente: 1. Se considera la permeabilidad magnética del material mucho mayor que la del vació o aire, es decir: um >> u0. 2. De acuerdo a la ley de Ampere, se tiene lo siguiente: NI = φnReq φn = BnAn Req = Reqm + Reqg Bn An =
NI NI ; Bn = An Rmeq + R geq Req
[
………… ( 1.26 )
]
Si: um << u0
Rm , g =
l m, g u m,c Am, g
; Bn =
NI An R geq
...........( 1.27 )
[ ]
3. Con el valor de Bn obtenido en la relación ( 1.27 ), obtenemos um y Hm de la característica B versus H; es decir: Bn(1), um(1), Hm(1) 4. Con um(1) obtenido reemplazamos dicho valor en la ecuación ( 1.26 ) y obtenemos Bn(2), um(2), Hm(2). 5. Luego se realizan las iteraciones sucesivamente hasta conseguir un porcentaje de error menor al 5 %;es decir: % Er ≤ 5%. NI( dato ) − NI( calculado ) Siendo: Er = x100 NI( dato )
Donde:
∑ H .l = NI ⇒NI
( calculado )
= H mlm +
………… ( 1.28 )
λ g An B n
u 0 Ag
………… ( 1.29 )
CALCULO DE INDUCTANCIA: “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
18 I
λ
Fig. 1.7 Elemento inductivo λ: Acoplo inductivo. λ = Li
dφ d → υ = (N φ ) dt dt di d υ=L y si : L = cte → υ = (Li) dt dt Nφ Nφ= Li Æ L = i υ = N
además: L =
……….. ( 1.30 ) …… ….. ( 1.31 )
NBn An NµH n An NµNiAn →L= →L= i i iλ
…………. ( 1. 32 )
B N2 N2 →L = , donde: µ = n del gráfico B -H L= λ R Hn µA f
EFECTO DE BORDE:
a
Lg
b
Fig. 1.8 Influencia del entrehierro Siendo: “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
19 a: ancho del paquete de láminas. b: altura del paquete de láminas. n: numero de láminas. t: espesor de láminas. fa: factor de apilamiento. N: Número de espiras lg: longitud del entrehierro Área del entrehierro: Donde : Área del material:
Ag = (a + l g )(b g + l g )
……….. (1.33 )
nt fa
……….. (1.34 )
Am = (nt ).a
……….. (1.35 )
bg . =
Donde :
bef.= nt.
Ademas :
Am = bef.a = (fa.bg) a
……….. (1.36 )
1.5.2 EXCITACIÓN CON CORRIENTE ALTERNA
Los circuitos magnéticos al ser excitados con una fuente de tensión alterna, presentan pérdidas en el núcleo magnético y pérdidas en el conductor eléctrico ( generalmente se utiliza el cobre ), por lo que dichas pérdidas deben ser cuantificadas. Las pérdidas en el núcleo están conformadas por las pérdidas debidas a las corrientes parásitas( Eddy) y las pérdidas por Histéresis. Pérdidasnucleo = Phistéresis + Pfoucoult
............ ( 1.37 )
Pérdidas en el Cobre( efecto Joule): I2R
............ ( 1.38 )
Las pérdidas en el cobre están expresadas en función de la resistencia y la corriente que pasa por el conductor (Joule). Siendo :
R = ρ L/A
............ ( 1.39 )
Donde: “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
20 ρ: Resistividad del material. L: Longitud del conductor. A: Área de la sección transversal del conductor.
Formas de medir Resistencias:
Método Tensión – Corriente (V, I ) Se aplica una tensión continua entre los bordes extremos del conductor y se registra los valores de las corrientes que son tabuladas. El área de la sección transversal es menor que el área geométrica del conductor. Método del Ohmimetro Con el instrumento llamado ohmimetro, se coloca las puntas en los extremos del conductor y se mide la resistencia en corriente continua. Para ambos casos se deberá realizar la corrección del valor de la resistencia medida en corriente continua, siendo afectada por un factor de corrección K , es decir: RAC > RDC ………… ( 1.40 ) RAC = K RDC K>1 1 < K < 1,5 Para un caso particular, se tiene: k = 1,25 Pérdida en el Núcleo (Pnúcleo):
Las pérdidas en el núcleo del material ferromagnético son debida a la Histérisis y las corrientes parásitas, es decir: Pnúcleo = Phistéresis + Pfoucoult Phistéresis:
Phist = Kn f Bmaxn
………… ( 1.41 )
Donde: Kn: Cte. de histéresis Bmax: Densidad máxima magnética f: frecuencia “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
21 n: Exponente de Steimetz n< 2 ; función de las líneas de flujos magnéticos. Pérdidas por Foucoult (Corrientes parásitas) Pfouc = Kf f2 B2max
………… ( 1.42 )
(πt ) 2 6ρ
………… ( 1.43 )
Donde: Kf =
ρ = % concentración de Silicio. t = espesor de la lamina. f = frecuencia. Bmax = Densidad máxima magnética. Phistéresis Pfoucoult
Pnúcleo W
Pcobre
Æ Pnúcleo =
W
- RI2
......... ( 1.44 )
lectura Vatímetro
Además, se tiene:
Eeficaz = 4,44 f N Bmax Aneta
........... ( 1.45 )
E ef 4( f .f )fN.A neta
........... ( 1.46 )
Luego: B max =
Phisteresis = K n f
E ef
n
[ 4( f .f )f .N.A n ]
n
1 Kn x n−1 xE neq Phisteresis = n [ 4,4( f .f )N.A n ] f
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
22
Phisteresis =
K 1E neq
........... ( 1.47 )
f n−1
Pfoucoult = K f f
2
E eq
2
[4(f.f )NA n ]
2
xf
2
→ Pfoucoult = K 2E eq
2
........... ( 1.48 )
Determinación de las Pérdidas en el núcleo, cuando el Flujo magnético es constante( φ = constante )
Se deben expresar las pérdidas por histérisis y corrientes parásitas, cuando el flujo magnético se mantiene constante, es decir: Phist, foucoult Æ φ =cte = Bmax An Para una onda cualesquiera se tiene la expresión de la tensión eficaz inducida expresada por: Eef = 4( ff ) Nφ max E eq = 4( ff ) fNφ max → f 2 Pfoucoult = K f f 2 Bmax = bf
2
........... ( 1.49 )
n Phisteresis = K n fBmax = af
........... ( 1.50 ) Pnucleo = Phist + Pfouc = af + bf 2
........... ( 1.51 )
Pnuc = a + bf f
Para determinar los valores de las constantes a y b, se deberán realizar por lo menos dos ensayos a tensión y frecuencias diferentes, obteniéndose los resultados: p f 2 − p n 2 f12 2 a = n1 2 Pn1 = af 1 + bf1 f1 f 2 ( f 2 − f1 ) ........... ( 1.52 ) Pn 2 = af 2 + bf 2
2
b=
p n1 f 2 − p n 2 f 1 f1 f 2 ( f1 − f 2 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
23 Luego, para un tercer valor de tensión y frecuencia se tiene las pérdidas en el núcleo expresada por: Pn3 = af 3 + bf 3
2
........... ( 1.53 )
Siendo, las pérdidas por histérisis y Foucoult : Phistéresis
= af3
........... ( 1.54 )
Pfoucoult
= bf32
........... ( 1.55 )
1.6 Reactor con Núcleo de Hierro
El reactor de núcleo de hierro es un dispositivo electromagnético estático conformado por una estructura ferromagnética ( laminas ) y un paquete eléctrico (bobinas) tiene diversas aplicaciones en las instalaciones eléctricas industriales y en los sistemas eléctricos de potencia. Para el circuito equivalente, se plantea la ecuación de equilibrio electromagnético siguiente: o
v( t ) = λ( t ) = Ri( t ) +
∂ψ (t) ∂t
...........( 1.56 )
Siendo expresada el flujo magnético por: M
ψ K ( t ) = ψ KK + ∑ ψ iK
............( 1.57 )
i =1 i ≠K
El flujo magnético para una bobina está expresado por: ψ 1(t ) = ψ {11
= L 1i1
flujo propio
Luego, reemplazando en la expresión ( 1.56 ), se obtiene: v ( t ) = Ri( t ) +
∂ (Li) ∂t
“Estudio de las Máquinas Eléctricas. Simulación Digital”
............( 1.58 )
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
24 Desarrollando en derivada parciales la ecuación ( 1.58 ), obtenemos:
v( t ) = Ri( t ) + L
di dL dx +i . dt dt dx
v( t ) = Ri( t ) +
di L {dt
o
+ix
dL dx
...........( 1.59 )
Transformación
o
ix
dL : Tensión producida por movimiento ( traslación o rotación ) dx
Para el caso del reactor de núcleo de hierro que no experimenta movimiento de la bobina, obtenemos:
v( t ) = Ri( t ) + L
di o ; λ = v(t ) dt
...........( 1.60 )
ο
q : corriente que absorbe el reactor.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
25
CAPITULO 2
TRANSFORMADORES NÚCLEO DE HIERRO 2.1 Introducción. Forma Constructiva, Utilización, tipos. Los transformadores de núcleo de hierro son máquinas eléctricas que están conformados por estructuras magnéticas laminadas(de diversas configuración geométrica) y paquetes eléctricos conformados por dos o más devanados estacionarios que se encuentran acoplados magnéticamente. Tienen diversos usos como son: en los sistemas eléctricos de potencia, en las redes de distribución primarias y secundarias, en las subestaciones de las plantas industriales, en edificaciones de múltiple usos y centros de transformación. El principal uso de los transformadores son para el cambio de la magnitud de la tensión alterna. Por convención el devanado AC de entrada es usualmente referido como el devanado primario, y los otros devanados de salida son referido como devanados secundario o terciario. Los transformadores de núcleo de hierro son las máquinas eléctricas que presentan menores valores de pérdidas ( pérdidas en el núcleo, pérdidas en el cobre ) en relación a las otras máquinas eléctricas, como son: los generadores y motores eléctricos; por lo que los transformadores tienen un valor de alta eficiencia( pueden llegar hasta el 99 %). Trabajan a baja frecuencia entre 25 a 400 Hz, para aplicaciones de alta frecuencia se tienen núcleo de ferrita o núcleo con entrehierro para contrarrestar la excesiva pérdidas en el núcleo. Las pérdidas por corrientes de Eddy o corrientes parásitas en el núcleo de hierro pueden reducirse construyendo el núcleo laminado; para una frecuencia de 60 Hz la lamina del núcleo de hierro puede ser 0,35 mm de espesor. Los transformadores de potencia según su utilización, pueden ser: a. Transformador de unidad, conectado a la salida de un generador para elevar la tensión a los niveles de transmisión. b. Transformador de Subestación, conectado en el otro extremo de la línea de transmisión para reducir la tensión a los niveles de Distribución. “Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
26 c. Transformador de Distribución, conectados en las redes de distribución primaria para reducir la tensión a los niveles de utilización. Los devanados del primario y secundario están físicamente enrollados uno sobre el otro, el devanado de baja tensión está situado en la parte interna; esta disposición cumple dos objetivos: 1. Resulta menor flujo de dispersión en relación a disponer los dos devanados en el núcleo separados. 2. Simplifica el problema del aislamiento del devanado de alta tensión desde el núcleo. Los transformadores se construyen de dos tipos: 1. Tipo Núcleo, consta de una pieza rectangular de acero laminado con los devanados enrollados sobre dos de los lados. 2. Tipo acorazado, consta de un núcleo de acero laminado de tres columnas, cuyos devanados están enrollados en la columna central. 2.2 Transformador Ideal El Transformador ideal es un dispositivo que consta de un devanado de entrada y otro de salida, las relaciones entre las variables de tensión de entrada y tensión de salida, y entre las corrientes de entrada y de salida, están dadas por la relación entre el número de espiras de cada devanado del transformador. Un transformador ideal, tiene las características siguientes: 1. El flujo magnético está confinado en el núcleo, esto quiere decir que no existe flujo de dispersión. Las pérdidas en el núcleo son consideradas despreciables; es decir no deben tener pérdidas por Histérisis ni corrientes parásitas. 2. Las resistencias de los devanados primarios y secundarios son despreciables, por lo que las pérdidas en el cobre son nulas. 3. La permeabilidad magnética del núcleo ( µm ) es considerada infinita. Al ser conectado el devanado primario a una fuente de tensión ν1 variable en el tiempo, en el núcleo magnético se establece un flujo φm. Una tensión e1 es inducido en el devanado y será igual a la tensión aplicada si la resistencia del devanado es despreciable( transformador ideal ), es decir se cumple: ν1 = e1 = N1 dφm dt
....................( 2.1)
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
27 El flujo en el núcleo enlaza al devanado secundario y induce una tensión secundaria e2, siendo igual a la tensión en los terminales ν2, es decir se cumple: ν2 = e2 = N2 dφm dt
.....................(2.2)
Relacionando las ecuaciones ( 2.1 ) y ( 2.2 ), obtenemos: ν1 = N1 = a ν2 N2
.....................( 2.3 )
Donde ¨a¨ es denominada la relación de transformación. La ecuación ( 2.3 ) nos muestra que las tensiones en los devanados de un transformador ideal son directamente proporcional al número de espiras de los devanados. Asimismo, la fuerza magnetomotriz (FMM) neta que se establece en el transformador ideal es cero, es decir : N1 i1 – N2 i2 = FMM neta = 0
...........................( 2.4 )
N1 i1 = N2 i2
........................... ( 2.5 )
N2 = i1 = 1 N1 i2 a
..........................( 2.6 )
Relacionando la ecuación (2.3) y (2.6), establecemos: ν1i1 = ν2i2
.
.................... ( 2.7)
La ecuación ( 2.7 ) nos indica que la potencia instantánea de entrada del transformador es igual a la potencia instantánea de salida. Esto es debido a que todas las pérdidas son consideradas despreciables en un transformador ideal 2.2.1 Polaridad Los transformadores utilizan la convención de puntos que aparecen en cada extremo de los devanados, éstos indican la polaridad de la tensión y la corriente en el lado del transformador, la relación es la siguiente: “Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
28 a. Si la tensión primaria es positivo en el extremo del devanado marcado con un punto, respecto al extremo que no tiene marca; la tensión secundaria será positivo también en el extremo marcado con punto. b. Si la corriente primaria fluye hacia dentro del devanado primario por el extremo marcado con punto, la corriente secundaria fluirá hacia fuera del devanado secundario por el extremo marcado con punto. 2.2.2 Impedancias y Variables reflejadas Para el análisis operativo de los transformadores conectados bajo carga o en vacío, es de importancia de establecer el circuito equivalente del transformador referido al lado primario o secundario; para lo cual es recomendable realizar un tratamiento de las impedancias y variables reflejadas. Valores Reflejados: V2´= aV2
V1´=
V1 a
I1´= aI1 I I 2´= 1 a
- Variables Reflejadas. - Impedancias Reflejadas. ;
Tensión secundaria reflejada al lado primario.
;
Tensión primaria reflejada al lado secundario.
;
Corrientes primaria y secundaria reflejadas.
V1I1 = V2I2 Z1I12 = Z 2 I 22 →
...........( 2.08 ) Z1 = a2 Z2
..........( 2.09 )
Z1 = a2Z2 = Z2´
...........( 2.10 )
Z2´ = a2Z2
.......... ( 2.11 )
Z1 = R1 + jX1
.......... ( 2.12 )
Z2 = R2 + jX2
.......... ( 2.13 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
29 R2’+jX2’ = a2(R2 + jX2)
............( 2.14 )
1 ' R2 a2 1 X '2 = a 2 X 2 → X 2 = 2 X '2 a R '2 = a 2 R 2 → R 2 =
............( 2.15 )
2.3 Transformador núcleo de hierro Para el análisis del comportamiento dinámico del transformador se plantean las ecuaciones de equilibrio electromagnético y para el análisis en régimen estacionario se requiere el circuito equivalente. A continuación se planteamos las ecuaciones de equilibrio electromagnético:
V1 ( t ) = R1i1 ( t ) +
dψ 1 ( t ) dt
V2 ( t ) = R 2i 2 ( t ) +
dψ 2 ( t ) dt
............ ( 2.16 ) ..........( 2.17 )
Se denota: n
ψ K = ψ KK + ∑ψiK i =1 i≠K
ψ 1 = ψ 11 (t ) + ψ 12 (t )
............. (2.18 )
ψ 2 = ψ 22 (t ) + ψ 21 (t )
.............(2.19 )
Reemplazando las ecuaciones ( 2.18 ) y ( 2.19 ) en las ecuaciones (2.16) y (2.17 ), obtenemos las expresiones siguientes:
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
30
V1 (t ) = R1i1 (t ) +
dψ 11 dψ 12 (t ) + dt dt
V2 (t ) = R2i2 (t ) +
dψ 22 dψ 21 (t ) + ∂t dt
ψ 11 = L1i1 ψ 22 = L2i2 ψ 12 = M 12i2 ψ 21 = M 21i1 ............( 2.20)
Luego: V1 ( t ) = R1i1 ( t ) + L 1
di1 ( t ) di ( t ) + M12 2 dt dt
V2 ( t ) = R 2i 2 ( t ) + L 2
di 2 ( t ) di ( t ) + M21 1 dt dt ............( 2.21 )
La ecuación ( 2.21 ) esta expresada en términos de régimen transitorio o función del tiempo. El análisis en régimen estacionario, las ecuaciones en función del tiempo o ecuaciones diferenciales lineales son transformadas en ecuaciones de variable compleja que, a continuación son desarrollados:
V = V∠0º Análisis en régimen estacionario: d p = = jw dt ........ ( 2.22 ) Luego: V 1 = R1 I 1 + L1 jw I 1 + M 12 jw I 2 + jwM 12 I 1 − jwM 12 I 1 V 2 = R2 I 2 + L2 jw I 2 + M 21 jw I 1
........ ( 2.23 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
31 V 1 = R1 I 1 + ( L1 − M 12 ) jw I 1 + jwM 12 ( I 1 + I 2 ) V 1 = R1 I 1 + jwl1λ I 1 + jX m ( I 1 + I 2 )
........ ( 2.24 ) Análogamente: V 2 = R2 I 2 + jwl 2 λ I 2 + jX m ( I 1 + I 2 )
........ ( 2.25 ) Luego, los parámetros del circuito equivalente están expresados por:
Req = R1 + R21 xeq = x11 + X 21λ X eq R eq
θ eq = Tg −1
........ ( 2.26 ) Circuito Equivalente
I
I
1
E
E
2
2
1
I
I
M
p
Fig. 2.1 Circuito equivalente exacto de un transformador de núcleo de hierro “Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
32 2.4 Determinación de parámetros (condiciones nominales) Los parámetros de una máquina eléctrica son obtenidos a partir de la realización de pruebas de laboratorio, siendo estas: - Prueba de cortocircuito. - Prueba de vacío. Los parámetros de las máquinas eléctricas son obtenidos a condiciones nominales ( C.N. ), es decir en la prueba de cortocircuito se realiza a corriente nominal y en la prueba de vacío se aplica la tensión nominal 2.4.1 Prueba de cortocircuito (C.C): Con esta prueba se determinan los parámetros R1, R2, x1l, x2l.
Fig. 2.2 Circuito para la prueba de cortocircuito Procedimiento: 1. Se recomienda cortocircuitar el lado de baja tensión e instalar un amperímetro. 2. Se suministra la energía o corriente alterna por el lado de alta tensión (A.T), hasta conseguir que la corriente obtenida en el lado de Baja Tensión (B.T) sea la corriente nominal (INBT). 3. El valor de la tensión aplicada en el lado primario, para obtener la corriente secundaria equivalente a la nominal (INBT) es llamada Tensión de cortocircuito (Vcc) V Ucc = cc x100 Ucc: % tensión de C.C VN 4. La tensión de C.C porcentual esta comprendida en: Transformadores: 3% < Ucc < 15% Motores: 25% - 40%
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
33 Circuito Equivalente del Transformador: para determinar el circuito se requiere obtener los parámetros del transformador, es decir:
Z eq =
Vcc I cc
;
Z eq = Req cc + jX eq cc ........... ( 2.27 )
Cálculo de Reqcc: W
= I2Reqcc
R eq cc =
W I 12
R eq cc = R1 + R 12 ............ ( 2.28 ) Se cumple: R1 =
R12 =
R eq cc 2 R eq cc 2
R21 = a 2 R2 1 Req cc R21 ; R2 = a 2 2 R2 = 2 a ............ (2.29 ) análogamente obtenemos: X1 =
X 21 =
X eq cc 2 X eq cc 2
;
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
34 X 21 = a 2 X 2 X2 =
X 21 a2
X2 =
1 X eq cc a2 2
.......... (2.30 )
2.4.2 Prueba de vacío (C. A.) Se determinan los parámetros Rp (resistencia de perdidas) y Xm ( Reactancia de magnetización) .
Fig. 2.3 Circuito para la prueba de vacío Procedimiento: 1. Se recomienda mantener a circuito abierto el lado de A.T. 2. Se suministra la energía o corriente alterna por el lado de B.T, hasta conseguir, que la tensión aplicada en el lado de B.T, sea la tensión nominal (VNB.T). 3. El valor de la corriente en lado secundario para obtener la tensión secundaria equivalente a la nominal (VNB.T) es llamada corriente de vacío (Ic-a). I ; i0: % de corriente en vacío. i 0 = 0 x100 IN 4. La corriente de vacío porcentual esta comprendida en: Transformadores núcleo de hierro: 2% < i0 < 10% Para Transformadores de núcleo de aire: i0 < 80%IN “Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
35 Determinación de Rp: = Ph + R2 I22
W
Ph =
W
………. ( 2.31 )
- R2 I22 Æ Ph = RP IP2 ,
E = RP IP Æ I P =
RP =
RP =
Ph IP2
E P → RP = h2 xRP2 RP E
E2 Ph
........... (2.32 )
Además:
V 2 = E + I2 Z 2
;
Z 2 = R 2 + jX 2λ
E = V 2 − I2 Z 2 ........... (2.33 ) cos θ 2 =
W W → θ 2 = cos −1 → θ 2 = θ c.a VI VI
θ cc ≠ θ c.a sc
Luego:
Xm =
E ;Im = (I02 − IP2 )1/ 2 Im
............ (2.34 )
Vcci Vcc 2 = I ccN I cc 2 V Vcc = cc 2 xI ccN I cc 2 ........... (2.35 ) “Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
36 2.4.3 Prueba bajo carga Los Transformadores de núcleo de Hierro son diseñado y construidos para que trabajen a diferentes condiciones de carga, es decir las cargas pueden ser : resistiva, inductiva y capacitiva; según el tipo de carga los transformadores presentan unos indicadores operativos , como son : eficiencia y regulación. FACTOR DE POTENCIA
CARGA
Resistiva (R) : horno Predominante inductiva (XL): motor inducción Predominante capacitiva (XC): banco de capacitores.
- Potencia - Factor de Potencia
- Resistivo: cosφ =0.0 -Inductivo,capacitivo:cosφ <±1.0
Se realizan las pruebas bajo carga con ciertas consideraciones de las potencias y el factor de potencia, siendo estas según mostradas en el cuadro siguiente: . 100 KVA 25% 50% 75% 100% 110% - 120% 0,5 0,8 0,9 0,5 0,8 0,9 0,5 0,8 0,9 0,5 0,8 0,9 0,5 0,8 0,9 cosφ
2.4.4 Pruebas especiales Estas pruebas son realizadas a transformadores de alta tensión (p.ej. 220 kV ) y niveles de potencia de cientos de MVA; siendo estas pruebas:
• • •
Prueba a frecuencia industrial (50, 60 Hz) Pruebas impulsionales ( alta frecuencias ). Simulaciones. Pruebas de descargas parciales.
BiL: Nivel básico de aislamiento. Por ejemplo, para realizar ensayos a transformadores de potencia de 220 kV, la Norma IEC recomienda tener un generador de 500 KV.
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
37 Estas pruebas especiales se realizan para transformadores mayores de 100 MVA a niveles de tensión de 220KV. Se debe de tener en cuenta el BiL (nivel básico de aislamiento): interno e externo 2.5 Eficiencia ( η ) En las máquinas eléctricas, la eficiencia es de importancia medirla por lo que nos indica la relación de potencia o energía que entrega la máquina eléctrica a la carga con respecto a la potencia o energía que absorbe de la fuente, dicho indicador operativo se expresa en porcentaje, siendo de la manera siguiente:
η
=
P salida P entrada
x 100 .......... (2.36 )
En toda máquina eléctrica real existe una potencia o energía de pérdida, por lo que está expresada por : Pentrada = Psalida + Ppérdida
............ ( 2.37 )
Luego, reemplazando en la expresión de la eficiencia, tenemos:
n=
Psalida x100 Psalida + Pperdida ........... (2.38 )
La máquina eléctrica con mayor eficiencia creada por el hombre es el auto transformador , que alcanza eficiencia hasta el 99%. eficiencia = f (Pot, cosφ) El Organismo Latinoamericano de Desarrollo de la Energía (OLADE), para Transformadores de distribución recomienda que, la eficiencia máxima se obtiene cuando el transformador se encuentra operando entre el 80 % al 120 % de su potencia nominal. nmáx e < 80% - 120% > PN
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
38 Las pérdidas de potencia o de energía de una máquina eléctrica depende del material a ser utilizado; cuya expresión esta dada por: 5
E perd . = EcuTD + EFeTD = ∑α i2 xPcuN xH i + PFe xH i i =1
5
ETD = ∑α i2 xS N xH i x cos φ Li ; i =1
E TD E +E perd TD
η=
x100
........... (2.39 ) Psalida = VL IL cosφL Pperdida =
Pnucleo +
Pcu(bob)
PFe F(V)
I22R eq f(I)
Reflejado al primario, tenemos: Req = R1 + R 12 R=λ
L A
Por lo que, reemplazando las relaciones respectivas obtenemos: η=
VLIL cos φ L VLIL cos φ L + PFe + I22R eq
........... (2.40 ) La ecuación (2.40) nos representa la relación de la eficiencia para cualquier condición de carga, por lo que tomaría valores como se indican a continuación:
• • •
Eficiencia máxima ( ηmaxima. ) Eficiencia a plena carga ( ηp.c. ) Eficiencia a un porcentaje de plena carga ( η%pc ) 50% p.c. - sub carga 110% p.c. - sobre carga.
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
39 El factor de Carga ( α ), expresa la relación existente de la corriente que absorbe la carga entre la corriente nominal del transformador, expresada por:
α=
I2 ; I 2 = αI 2 N I2N
...........( 2.41 )
Reemplazando la ecuación ( 2.40 ) en ( 2.41 ), obtenemos:
η=
αVLI2 N cos φL αVLI2 N cos φL + PFe + α 2 (I22 NR eq ) 14 2 43 PcuN
............ (2.42 )
η=
αVL I 2 N cos φ L ; αVL I 2 N cos φ L + PFe + α 2 Pcu
∂n =0 ∂α ........... ( 2.43 )
Para obtener la eficiencia máxima, se deriva la relación de eficiencia con respecto al factor de carga (
∂n ∂α
= 0 ) o igualando las pérdidas en el fierro a las del cobre
a una condición de carga cualesquiera ( PfeN = α2 PcuN ), obteniéndoselo siguiente:
α=
PFeN PcuN
............ ( 2.42 )
2.5.1 Eficiencia durante todo el día (ηTD ) Para poder determinar la eficiencia durante todo el día ( ηTD ) en un transformador de núcleo de hierro se debe conocer y/o analizar el diagrama de carga diario y ver como ha ido desarrollando la tensión de carga ( VL ) y su factor de potencia ( cos φL ). El diagrama de carga diario deberá ser potencia versus tiempo.
ηTD =
PsalidaTD PsalidaTD + PfeTD
ηTD =
Energia TD EnergiaTD + Energia PérdidaTD ........... (2.43 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
40 Energía = Potencia x tiempo x cosφ En el cuadro adjunto, se muestra la variación de la potencia de la carga y de su factor de potencia durante todo el día, siendo: %PC Pot. Horas CosφL
0% 0 6 0
50% 75% 100% 110% 50KVA 75KVA 100KVA 110KVA 6 6 3 3 0,6 0,8 0,9 1,0
2.6 Regulación (%r) La regulación es un indicador operativo que expresa la relación existente entre la variación de la tensión de salida en la carga con respecto a la tensión sin carga; usualmente se indica en porcentaje, siendo:
%r =
V entrada − V salida V salida
x100
.......... (2.44 ) Asimismo, el porcentaje de regulación se puede expresar de la manera siguiente:
1) %r =
Vvacio − Vc arg a Vc arg a
2
x100 ............. (2.45 )
2) %r =
Vvacio − Vc arg a Vvacio
x100
............. ( 2.46 )
3) %r = u cca cos φL + u ccr senφL +
(u cca senφL − u ccr cos φL ) 2 200 ............ ( 2.47 )
En la ecuación (2.47 ) se tiene las relaciones siguientes:
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
41
2 2 1/ 2 ucc = ucca + juccr ; ucc = (ucca + uccr ) ;
2 1/ 2 uccr = (ucc2 − ucca )
ucc =
Vcc x100% VN
Z cc =
Vcc I cc
........... (2.48 ) u cca
u cca
IccN Rcc IN = x100 x VN IN 64 P7cuN48 I 2 ccN Rcc 10 x100 x = VN I N 10
u cca =
PcuN (watt ) 10PN (KVA)
........... (2.49 )
u θ cc = Tg −1 ccr u cca donde: ucc : Tensión de cortocircuito porcentual equivalente ucca: Tensión de cortocircuito porcentual de la componente activa. uccr: Tensión de cortocircuito porcentual de la componente reactiva. Para una carga inductiva , se tiene la relación siguiente:
V12 = (V21 ) 2 + ( Z eq I 21 ) 2 + 2V21Z eq I 21 cos(θ T − φ L ) 1 4 2 4 3 V12 = (V21 + Z eq I 21 ) 2 V1 = V21 + Z eq I L → θ T − θ L = 0 → θ T = θ L ............. (2.50 ) “Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
42 Por lo que, para regulación máxima se tiene: %Rmax Æ φL = θT Para carga inductiva. Para porcentaje de regulación mínima se tiene: %Rmin Æ φL ≠ θT
2.6.1 Regulación de tensión bajo carga Los transformadores de núcleo de hierro utilizados para los sistemas de distribución y de potencia, el devanado primario es construido con tomas variables para obtener la regulación de tensión variando el número de espiras de dicho devanado. En la práctica, las cargas a ser servidas se encuentran distantes del centro de generación o de las subestaciones transformadoras, por lo que es necesario regular la tensión del primario de los transformadores; comercialmente se proporciona las tomas siguientes: a)
( ±2 x 2,5%)xV1
b)
( ±5 x 2,0%)xV1
c)
( ±10x 1,0%)xV1
.......... (2.51)
Los transformadores presentan una impedancia equivalente, que al paso de la corriente eléctrica cuando se conecta la carga en el lado secundario se produce una caída de tensión, ocasionando una variación de la tensión en la carga por lo que es necesario variar la tensión de entrada del primario para garantizar que la tensión en la carga se mantenga constante, esta situación es dada cuando la solicitud de potencia de la carga es variable. Para el caso de un transformador que es representado como modelo eléctrico por una impedancia equivalente, al ser conectado en el lado secundario una carga a través de un conductor, la tensión en la carga se encuentra expresada por la ecuación fasorial siguiente: V A = Z T + ZC I 1 + V L V A = (Z T + Z C )I j + V L ........... (2.52 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
43 2.7 Calentamiento en las máquinas eléctricas En las máquinas eléctricas, una de las variables de importancia a ser medida y controlada es la temperatura, por lo que se requiere utilizar técnicas adecuadas para determinar los valores de temperatura en tiempo real en las partes constituyentes de una máquina eléctrica , como son : el núcleo magnético y las bobinas primarias (estator) y secundarias( rotor ). Existe un elemento utilizado para la medición de la temperatura en las máquinas eléctricas, conocido como: RTDS ( Detector Resistivo de Temperatura ). 2.7.1 Balance termodinámico Aplicando la primera ley de la Termodinámica, tenemos:
P: Potencia Pdt = Cdt + K(θ − θa)dt C: Capacidad K: Disipacion θ: Temperatura. Luego, obtenemos:
θ(t) = θa + θf (1 - e-t/T) ........... ( 2.53 )
Por ejemplo , para algunas máquinas rotativas , se tiene T: 5 periodo En la práctica existen diversos modelos termodinámicos para realizar un análisis riguroso de la temperatura en las máquinas eléctricas, donde se modelan las partes del rotor, estator, núcleo y la relación entre ellos y el medio ambiente. VALORES POR UNIDAD (P.U.)
Los valores por unidad son de mucha utilidad en los procesos de cálculo en el campo de la ingeniería eléctrica, en particular en las máquinas eléctricas, por lo que se obtienen resultado numérico con menor porcentaje de error. La expresión del valor por unidad esta dado por la relación existente entre el valor real y el valor base de la variable o parámetro a ser evaluado, siendo dada por: V V p.u = Re al VBase ........... (2.54 ) “Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
44
VReal : Valor medido. VBase : Valor de referencia.
Z p.u =
Z Base
Z Re al . Z Base 2 VBase = SBase
.......... (2.55 )
; n p.u .sistema = X p.u .propia
Z Bpropia Z Bsistema .......... (2.56 )
2.9 Autotransformadores
Son máquinas eléctricas que presentan un solo bobinado en cada fase del autotransformador, existiendo una conexión física entre la entrada y la salida de la potencia o energía de cada bobinado. Los Autotransformadores tienen diversas aplicaciones como son: en los Sistemas eléctricos de potencia y en el arranque de motores de inducción. Para el arranque de motores de inducción 3φ mediante la utilización por autotransformador 3φ, la norma VDE 0530 recomienda que se aplique una tensión del 50% del valor de la tensión nominal, es decir: VAplicada ≈ 50% Vn Lo indicado por la Norma VDE 0530 es referencial por que, en algunos casos para motores de inducción trifásicos de gran potencia ( p.ej. 1 200 kW ), la tensión aplicada para el arranque sería del 75 % de la tensión nominal, esto es debido a que en el arranque el motor tiene que ser capaz de vencer su gran masa inercial. 2.9.1 Circuito equivalente
El autotransformador de núcleo de hierro, de acuerdo al análisis de las ecuaciones de equilibrio electromagnético, presenta los valores de su parámetros dadas por:
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
45 Z eq autotransf . = Req autotransf . + jX eq autotransf . Req autotransf . = R1 + R2 (a − 1) 2 X eq autotransf . = X 1λ + X 2 λ (a − 1) 2 Z eq autotransf . < Z eq transf .
............ (2.57 ) 2.9.2 Ventajas y Desventajas con respecto al transformador
1.- El auto transformador presenta menor impedancia equivalente que el transformador 2.- El auto transformador presenta una mayor regulación de tensión que el trabajo. 3.- El auto transformador presenta una mayor capacidad de corriente de corto circuito. 4.-El auto transformador presenta una mayor eficiencia que el transformador (> 99%). Ejemplo: Se tiene un transformador monofásico de 1,0 KVA, 220/110, 60Hz. Se requiere hacerlo trabajar como un autotransformador. Calcular la máxima potencia a entregar sin sobrecargarse. I ba =
1000 = 9,09amp. 110
I bc =
1000 = 4,54amp. 220
I c arg a = I ba + I bc = 13,63amp. KVA =
330 x13,63 = 2,9986 KVA ≈ 3KVA 1000
“Estudio de las Máquinas Eléctricas. Simulación Digital”. Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
46
CAPITULO 3 TRANSFORMADORES EN SISTEMAS TRIFÁSICOS 3.1 Introducción Los transformadores trifásicos son aquellas máquinas eléctricas que tienen múltiples aplicaciones en el sector eléctrico, siendo utilizados en: Centrales Eléctricas, centros de transformación, subestaciones de distribución, plantas industriales, electrificación urbana y rural, centros comerciales y otros. Debido a su importancia se requiere realizar un análisis del comportamiento de operación bajo carga cuando se encuentran conectados en paralelo dos o más unidades, asimismo de especificar las características técnicas en cuanto a su tipo y grupo de conexión para realizar una correcta puesta en paralelo e indicar las condiciones de las mismas. Los transformadores trifásicos pueden ser de dos arrollamientos ( primario y secundario ) y de tres arrollamientos
( primario, secundario y terciario ). A
diferencia de los transformadores monofásicos, presentan otras características sus arrollamientos ,como son: a.- Transformadores trifásicos de dos arrollamientos.- están denotados por: Nivel de tensión
: V1/V2
Potencia
: kVA , MVA
Impedancia porcentual de cortocircuito : %Xpu Tipo y grupo de conexión
: Ynd11
b.- Transformadores de tres arrollamientos.- están denotados por: Nivel de tensión
: V1 / V2 / V3
Potencia
: S1 / S2 / S3
Impedancia porcentual de cortocircuito : %X1 / %X2 / %X3 “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
47 Tipo y grupo de conexión
: Ynynd5
A continuación, se indican una relación de características que debe tenerse en cuenta para la correcta especificación técnica de los transformadores trifásicos, que serán utilizados en las diversas instalaciones eléctricas 3.2 Especificaciones técnicas Para la solicitud de adquisición de los transformadores trifásicos se deberán tener en cuenta los criterios de operación para lo cual se indicaran la potencia, nivel de tensión, condiciones de operación, altitud, nivel de aislamiento entre otros, es decir las características técnicas siguientes: 1.- Potencia: kVA, MVA 2.- Nivel de tensión: kV 3.- Frecuencia : Hz. 4.- Pérdidas: W, kW - núcleo - cobre 5.-Tipo de conexión: Estrella- delta Delta-estrella. 6.- Grupo de conexión: 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,12 7.- Porcentaje de tensión de cortocircuito: %µ 8.- Tipos de refrigeración: ONAN, ONAF. 9.- Nivel de altura de trabajo: msnm. 10.- Características del neutro:
- aislado. - directamente a tierra. - a través de impedancia.
11.- Nivel básico de aislamiento ( BIL): kV - Interno. - Externo. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
48 12.- Regulador de tensión bajo carga: ( +-5x2,0 % )V1 ( +-10x1,0 % )V1 13.- Sistema de Protección: - Tierra (Puesta a Tierra). - Relé diferencial
3.3 Conexionado de transformadores trifásicos Los transformadores trifásicos para ser conectados en paralelo, uno de los requisitos que se deben tomar en cuenta es el tipo y grupo de conexionado, por lo que se debe verificar su correcta conexión para evitar que se presenten diferencias de potenciales entre diferentes puntos de una misma fase . 3.4 Tipos y grupos de conexión Dentro de las características de los transformadores trifásicos se encuentran el tipo y grupo de conexión que son de importancia para la operación de dos o más transformadores conectados en paralelo, siendo descrito de la manera siguiente: 3.4.1 Tipo de conexión.- Los tipos de conexión pueden ser estrella, delta, zigzag, delta abierto o la combinación de cualesquiera de los mencionados, en el siguiente cuadro se establecen los diversos tipos de conexionado: TIPOS - Estrella – estrella: - Estrella – delta: - Delta – estrella: - Delta – delta: - Estrella – zigzag: - Delta – zigzag: - Delta - abierto
AT BT (mayúsculas) (minúsculas) Y y Y d D y D d Y z D z
Cuadro 3.1 Tipos de conexionado “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
49 Dependiendo de la utilización de los transformadores trifásicos son requeridos los tipos de conexionado por ejemplo: en los centro de generación se tiene el tipo delta- estrella, en los centros de transformación del tipo estrella- delta, para electrificación rural del tipo delta abierto.
Fig. 3.2 tipos de conexionado a) Estrella-estrella, b) Delta-delta 3.4.2 Grupos de conexión .- De acuerdo a las recomendaciones de la Norma IEC, se tienen en cuenta los grupos de conexión siguientes: Grupos
Índice horario
Grupo I Grupo II Grupo III Grupo IV
0, 4, 8, 12 2, 6, 10 1, 5, 9 3, 7, 11
Como podemos observar, los grupos I y II son aquellos que tienen los índices pares y los grupos III y IV son aquellos que cuentan con los índices impares. En la “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
50 práctica para realizar la puesta en paralelo de dos o más transformadores de un mismo grupo o de diferentes grupos se deberán realizar corrimientos o permutaciones de cada fase, según sea el caso. Cada índice representa el desfase vectorial existente entre la tensión primaria con respecto a la tensión secundaria, por lo que podemos denotar dos características que son : a. Permutación Se deben realizar entre los transformadores que presentan índices horarios de diferentes grupos, por ejemplo: - Grupo III de índice 1 se permuta ( 180º ) con el grupo IV de índice 11 - Grupo III de índice 5 se permuta ( 180º ) con el grupo IV de índice 7 b. Corrimiento Se deben realizar entre los transformadores que presentan diferentes índices horarios de un mismo grupo, por ejemplo: - Grupo III de índice 1 se realiza corrimiento ( 120º ) con el de índice 5. - Grupo III de índice 11 se realiza corrimiento ( 120º ) con el de índice 7. En el caso que se presentara realizar permutación y corrimiento a un transformador trifásico o banco de transformador trifásico conformado por transformadores monofásicos, como regla práctica se deberá tomar en cuenta lo siguiente: 1.- Primero, se debe realizar corrimiento solamente en el lado secundario.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
51 2.- Luego se debe realizar permutación en el lado primario así como en el lado secundario. Regla:
1º Corrimiento 2º Permutación
: 120º : 180º
Por ejemplo, dentro los transformadores utilizados comercialmente de los grupos III y IV, tenemos: Grupo III
:
1
5
Grupo IV
:
7
11
A continuación presentamos algunos casos: A.- Corrimiento: 1 Æ 5
}corrimiento
7 Æ 11 B.- Permutación 1 Æ 11
} permutacion
7 Æ 5 C.- Corrimiento y permutación 1 Æ 7
;
5 Æ
11
1°) 1 Æ 5 : Corrimiento 2°) 5 Æ 7 : Permutación
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
52 3.5 Banco de transformadores trifásicos Los bancos de transformadores trifásicos son utilizados en los sistemas eléctricos de potencia como son en los centrales eléctricas, en los centros de transformación y subestaciones de distribución. Pueden ser utilizados como una unidad trifásica o un banco de transformadores trifásicos conformado por unidades monofásicas. 3.5.1
Casos
que
se
presentan
en
el
Conexionado
de
banco
de
transformadores A continuación se plantean los posibles casos en que se presentan para los conexionados de los transformadores o bancos trifásicos conformados por transformadores monofásicos, siendo dados por: a.- Se tiene como dato el conexionado, es decir el tipo y grupo lo que se solicita es el conexionado físico de las bobinas, ejemplo Dy1. b.- Se tiene el conexionado físico de las bobinas del transformador o banco trifásico de transformadores monofásico y se solicita el grupo de conexión a que pertenece.
Fig. 3.3 Banco trifásico tipo D-y conformado por transformadores monofásicos “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
53 3.- Se tiene el tipo de conexión y a través de ensayos de laboratorio, se solicita el grupo de conexión, esta parte es conocida como la determinación experimental del índice de conexión. 3.6 Pruebas a banco de transformadores trifásicos Para obtener los parámetros eléctricos de los bancos de transformadores trifásicos es necesario realizar las pruebas en laboratorio siendo las pruebas de cortocircuito y de vacío ( circuito abierto ), estos parámetros son utilizados para establecer el circuito equivalente por fase de los bancos de transformadores y utilizados en el análisis de la operación del transformador conectado bajo carga, ya sea como una unidad trifásica o conectado en paralelo con otras unidades trifásicas. 3.6.1 Prueba de corto circuito Esta prueba se realiza para obtener los parámetros R1, R2, X1 e X2 del banco trifásico, se recomienda aplicar tensión por el lado de alta tensión y cortocircuitar el lado de baja tensión, según mostrado en la figura 3.4.
Fig. 3.4 Esquema eléctrico para la prueba de cortocircuito
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
54 3.6.2 prueba de vacío Esta prueba se realiza para obtener los parámetros Rp, e Xm del banco trifásico, se recomienda aplicar tensión por el lado de baja tensión y mantener a circuito abierto el lado de baja tensión, según como se muestra en la figura 3.5.
Fig. 3.5 Esquema eléctrico para la prueba de circuito abierto 3.7 Puesta en paralelo de transformadores En
las
diversas
Instalaciones
Electromecánicas,
existen
transformadores
operando como una unidad trifásica o transformadores conectados en paralelo, dependiendo de las características de operación de las cargas. Para la operación de dos o más unidades de transformadores conectados en paralelo, se deberán tener en cuenta las condiciones siguientes: 1. Igualdad de Potencia (KVA) 2. Igualdad de Frecuencia (Hz) 3. Igualdad de Tensión de Cortocircuito Porcentual 4. Igualdad de relación de Transformación 5. Igualdad de Grupo de conexión.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
55
En la práctica es poco probable que se cumpla con los cincos puntos planteados por lo que, se deberán establecer medidas correctivas para garantizar un buen conexionado para la puesta en paralelo de dos o más transformadores U
V
W
U ρ A
ρ B
ρ C
u
v
ρ b
ρ a
ρ A
ρ c w
t
w
Banco
r
u
Yd1 : Grupo III (Referencia)
ρ a
ρ C
ρ B
ρ b
W
ρ c
1
V v s
U
V
W
U ρ A
ρ C
ρ B u
ρ c
v
ρ b
ρ a
w
ρ C
v
ρ b
W
ρ c
U
ρ A
ρ a
V
ρ B
ρ C
u
v
ρ b
Banco
ρ a
Yd5 : Grupo III
ρ B
V u
W u
ρ A
w
ρ c
ρ a w
ρ C
5
1 1ρ A
ρ c
Banco ρ b ρ B
Yd11 : Grupo IV v
w
Fig. 3.6 Conexionado físico de las bobinas de los bancos de transformadores
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
56
R S T
Yd1
U u
V v
W w
Yd5
U u
V v
W w
Yd11
U u
V v
W w
r s t
Fig. 3.7 Conexionado de la puesta en paralelo de 3 bancos de transformadores 3.8 Reparto de carga Cuando dos o más transformadores son conectados en paralelo, para su análisis de la contribución de potencia a la carga es necesario establecer las condiciones de cada unidad transformadora; es decir su potencia, impedancia porcentual de cortocircuito y su relación de transformación.
Fig. 3.8 Conexionado en paralelo de dos transformadores “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
57 Para determinar la contribución o aporte de potencia de cada unidad transformadora hacia la carga, es importante establecer el circuito equivalente por fase, según mostrado en la fig. 3.7 donde se debe indicar todas las variables eléctricas y los parámetros por fase. En el caso de los transformadores trifásicos se deberá tener en cuenta que las relaciones de transformaciones serán por fase .
Fig.3.7 Circuito equivalente por fase del conexionado en paralelo Con respecto a la relación de transformación, se presentan dos casos, siendo: CASO 1. Cuando los transformadores presentan relación de transformación iguales, es decir: a1 = a2 .....= an Del análisis circuital se demuestra que:
SK =
Y K S NK n
∑Y
i
xS L
S Ni
i =1
........... ( 3.1 ) Donde: Y K : Admitancia en p.u del transformador k-esimo. S NK : Potencia nominal del transformador k-esimo. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
58
S L : Potencia aparente de las cargas. CASO 2.- Cuando los transformadores presentan relación de transformación diferentes, es decir: a1 , a2 , ….. , an. son diferentes.. En la práctica esta situación planteada se presenta con mayor frecuencia, por lo que amerita su análisis de reparto de carga resolviéndose aplicando las ecuaciones de mallas. Siendo las siguientes: Malla `1´ :
Malla `2´ :
V1φ
( )( )
…....... ( 3.2 )
V 1φ = Z 2φ I 2φ + V Lφ a 2φ
…....... ( 3.3 )
V 1φ = Z nφI nφ+ V L φ a nφ
...........( 3.4 )
a1φ
= Z 1φ I 1φ + V Lφ
. . . Malla `n´ :
En el nodo N se cumple: I1φ + I2φ + ….. + Inφ = ILφ
.......... ( 3.5 )
El conjunto de ecuaciones formuladas, se plantean en forma matricial de la manera siguiente: [ V - V L ] φ = [Z] [I] φ a
.......... ( 3.6 )
Asimismo, la expresión de la potencia sería: S
Kφ
=V 1φI
Kiφ φ
S3φK = 3SφK
“Estudio de las Máquinas Eléctricas. Simulación Digital”
.......... ( 3.7 ) .......... ( 3.8 )
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
59 Como un caso particular desarrollamos las ecuaciones de equilibrio para determinar la potencia que entregarían cada transformador a la carga, siendo las siguientes: Malla `1´ :
V 1φ = Z 1φI 1φ +V a2 φ
Malla `2´ :
V 1φ = Z 2φI 2φ+V a 2φ
Lφ
Lφ
……….. ( 3.9 )
…..........( 3. 10 )
Ecuación de nodo: I1φ + I2φ = ILφ Cuando dos transformadores
……..... ( 3. 11 ) de diferentes relaciones de transformación son
conectados en paralelo, se presenta la corriente circulante por lo que la corriente que contribuye cada unidad trifásica tiene dos componentes que son : la corriente de carga ( I carga ) y la corriente circulante ( I circulante ), es decir:
I1φ = Icarga + Icirculante
……..... ( 3. 12 )
Reemplazando la ecuación ( 3.11 ) en las ecuaciones ( 3.9 ) y ( 3.10 ) y reagrupando términos obtenemos las corrientes I1φ e I 2 φ , es decir: 1 1 − V1φ a a 2φ Z2φ 1φ I1φ = xI Lφ + ( Z1φ + Z2φ ) Z1φ + Z2φ
I 2φ
1 1 − V a 1φ a Z1φ 1φ 2φ = I Lφ − ( Z1φ + Z2φ ) (Z1φ + Z2φ ) ........... ( 3.13 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
60 Para determinar la expresión de la tensión primaria V 1φ , la expresión de la corriente I1φ es reemplazada en la ecuación (1) y reagrupando términos, tenemos: 1
1
V 1φ =
1
1
Z 1φ Z 2φ I Lφ + ( Z 1φ + Z 2φ )V Lφ Z 11φ Z 12φ + a2φ a1φ ............. ( 3.14 )
Es importante tener en cuenta en los transformadores que, la relación de transformación por fase con respecto a la relación de transformación trifásica depende del tipo de conexionado, es decir: Para el caso de tipo de conexión Y-d, tenemos:
a 3φ =
Vf VL1 3Vf 1 1 = = 3 VL2 Vf 2 Vf 2 14 2 43 a1φ
............. ( 3.15 )
a 3φ = 3∂1φ → a 1φ =
1 ∂3φ 3 ............. ( 3.16 )
Para el caso de tipo de conexión D-y, tenemos:
a 3φ =
VL1 = VL2
Vf 1 Vf1 = 1 x 3Vf 2 Vf 2 3 123 a 1φ
a 3φ =
1 a1φ → a1φ = 3∂ 3φ 3 ............. ( 3.17 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
61 3.9 Sistemas eléctricos desbalanceados En los sistemas eléctricos de potencia y en las máquinas el eléctricas, se presentan diversos tipos de contingencias o eventos de características eléctricas y/o mecánicas, ocasionando transitorios en el sistema o máquina eléctrica; por lo que amerita su análisis de los sistemas desbalanceados. Uno de los métodos de análisis de sistemas desbalanceados es mediante la componentes simétricas, planteándose circuitos equivalentes para los diversos componentes del sistema eléctrico, como son: generadores, líneas de transmisión, motores, cargas, transformadores de potencia y otros. Una variable vectorial puede ser expresada en función de tres componentes, es decir: V = f (V1 ,V2 ,V0 )
Donde: V1: Componente de secuencia positiva ( directa ). V2: Componente de secuencia negativa ( inversa ). V0: Componente de secuencia cero ( homopolar ). La configuración del neutro en un sistema o máquina eléctrica, puede ser : ·
Neutro aislado.
·
Conectado directamente a tierra.
·
Conectado a tierra mediante una impedancia ( Zn )
3.5 Transformadores de tres arrollamientos Son aquellos transformadores que presentan en su estructura 3 arrollamiento, siendo estas: Arrollamiento primario (A.T.) Arrollamiento secundario ( M.T) Arrollamiento terciario ( B.T) “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
62 En su forma constructiva el devanado terciario es instalado cercano al núcleo magnético, el devanado secundario a continuación del primario y por último el devanado terciario. Son utilizados generalmente en los sistemas eléctricos de potencia (SEP). Las especificaciones técnicas son similares a los transformadores trifásicos, con el adicional que se deben especificar las características del tercer arrollamiento, es decir: V1 / V2 / V3
;
220 / 60 / 10KV
S1 / S2 / S3
; 30 / 10 / 10 MVA
%X1 / %X2 / %X3
; 3,5% / 4,0% / 4,5%
3.10.1 Circuito equivalente Para establecer el circuito equivalente del transformador de tres arrollamiento, previamente es conveniente determinar las reactancias porcentuales de cada lado, esto es debido a que los fabricantes proporcionan las reactancias correspondiente vista entre dos lados del transformador, para lo cual se plantean las ecuaciones siguientes: XPS = XP + XS
…..…..( 3.18 )
XPT = XP + XT
…..…..( 3.19 )
XST = XS + XT
……....( 3.20 )
Resolviendo, obtenemos las reactancias siguientes: XP =
X PS + X PT − X ST 2
XS =
X PS + X ST − X PT 2
XT =
X PT + X ST − X PS 2
“Estudio de las Máquinas Eléctricas. Simulación Digital”
............ ( 3.21 ) Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
63
CAPÍTULO IV
4. CONVERSIÖN DE ENERGIA ELECTROMECÁNICA 4.1 Introducción El análisis de la Conversión de energía electromecánica es de importancia en el estudio de las máquinas eléctricas, por lo que se establecen las ecuaciones de equilibrio electromecánicos de la máquina eléctrica lográndose comprender la transformación de la energía eléctrica en energía mecánica y viceversa. Dentro las características de las máquinas eléctricas, podemos establecer lo siguiente: A.- Generadores eléctricos, de acuerdo a su conexionado y tipo de generación de tensión, tenemos: - Generadores AC ( Alternadores ) : Síncronos y Asíncronos - Generadores DC. B.- Motores eléctricos, de acuerdo a su conexionado y tipo de fuente de alimentación de tensión, tenemos: - Motores AC
A) Asíncronos : -Anillos rozantes: Rotor Bobinado. -Jaula de ardilla B) Síncronos: compensadores
- Motores DC
- Shunt. - Serie. - Excitación Independiente. - Coumpound.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
64 4.2 LEY UNIVERSAL DE FARADAY La Ley Universal de Faraday es aplicable para el estudio de los generadores eléctricos, siendo expresada por: ei = etransformacion + erotacion ei =
∂B
∫ ∂t xdS + ∫ v{ x B.dλ S
6 rx w
............ ( 4.1 )
B= B0 cosα B B= Bmaxcos(wt−α)
........... (4.2 )
fresult = fcampo µ fmecan. ( B variable en el tiempo )
............ ( 4.3 )
Cuando B es constante la frecuencia inducida es la misma que la máquina prima. 4.3 Sistemas eléctricos y mecánicos Para el análisis de la Conversión de energía electromecánica, se requiere tomar en cuenta : los elementos puros, variables de estado, funciones estado energía y coenergía siendo expresados en el cuadro siguiente: SISTEMA Eléctrico
Mecánico
ELEMENTO Inductancia Capacitancia Conductancia B.1 Traslación Masa Elástico B.2 Rotación Inercia Elástico
SIMBOLOGIA VARIABLES DE ESTADO L λ , q C λ , q G
FUNCIONES DE ESTADO Wm, W´m Wm, W´m
M K
P, x P, x
T , T´ V, V´
J Ko
l, l,
T , T´ V, V´
θ θ
Cuadro 4.1 representación de las variables y funciones de estado de energía y Coenergía. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
65 4.3.1 Análisis para un elemento inductivo En las máquinas eléctricas rotativas tales como generadores y motores tienen arrollamientos o bobinas en el estator y rotor , por lo que existe un predominio numérico del valor del parámetro inductivo; siendo desarrollado como variable de estado. El área encima de la curva representa la energía magnética: Wm El área debajo de la curva representa la coenergía magnética: Wm´. (Wm´): La coenergia magnética no tiene un significado físico pero si matemático. o
o
dWm = q dλ → Wm = ∫ q dλ o
...............(4.4)
o
dWm1 = λd q → Wm1 = ∫ λd q o
Wm + Wm1 = λ q ........... (4.5 )
Fig. 4.1 Representación grafica de la energía y coenergía para un elemento inductivo. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
66 Elemento inductivo:
dλ
L=
o
o
→ dλ = Ld q
............ (4.6 )
dq * L = cte ⇒ Lineal (Wm = Wm1 ) Wm (λ , x) = ∫ q dλ → ∫ q Ld q o * L = f (q ) ⇒ No Lineal o
o
o
............ (4.7 ) o
o
Wm1 (q, x) = ∫ λd q ........... (4.8 ) Si L = cte. W m ( λ, x ) =
1 o2 Lq 2
........... (4.9 )
4.3.2 Análisis para un elemento de elasticidad y masa Ecuaciones de equilibrio mecánico en función de estado de energía. Elemento elasticidad, tenemos: La
V (Energia Potencial ) =
o
∫ p a( X
1 a
)dX a1
0
∂V ( Xa) = pao ( Xa) Fuerza sobre el resorte. ∂Xa
Para el elemento masa: o
o
∫ dT = ∫ x dp → T = ∫ x dp o
o
1 1 ∫ dT = ∫ pd x → T = ∫ pd x
0
Xb
T1 =
∫
θ
............ (4.10 )
o1 o1 pb1 X b d X b
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
67
∂T 1 o
∂ Xb
d ∂T 1 o = pb Fuerza sobre la masa. dt ∂X b
o
= pb ( X b ) →
............ (4.11 )
Por el equilibrio de fuerzas: f(t): Fuerza sobre el resorte, Fuerza sobre la masa.
f (t ) =
∂ V ( Xa ) ∂ Xa
d + dt
∂T 1 o ∂ X b
........... ( 4.12 )
Puerta: Es aquel punto por el cual ingresa o sale la energía.
4.4 Forma restringida de la ecuación de Lagrange
Para el elemento de elasticidad: o1
(
)
V ( X α , X β ,.....) = ∑ ∫ p K X α1 , X β1 ,..... dX K1
........... ( 4.13 )
Para el elemento de masa: o
o
o
T 1 ( X a , X b ,.....) = ∑ K
XK
∫ 0
o o o1 p K1 X a , X b ,..... d X K
........... ( 4.14 )
Por lo tanto la ecuación general, resulta: d ∂T 1 dt 0 ∂ X K
o o o ∂V ( X 1, X 2 ,....., X n ) = QK X 1, X 2 ,....., X n + ∂X K ........... ( 4.15 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
68 Como ejemplo de aplicación, analizaremos para un caso donde un motor eléctrico es conectado a una carga, se establecen las ecuaciones de equilibrio. QK: Fuerza sobre el nodo ´K´. Coenergia Cinética del Sistema: 1 o o o o o T 1 θ 1,θ 2 = ∑ ∫ λ g θ 1,θ 2 d θ i
........... ( 4.16 )
Energía Potencial (elemento elástico) θ3
V (θ3 ) = ∫ λ13 (θ31 )dθ31
........... ( 4.17 )
0
Considerando elementos lineales: o
λ1 = J 1 θ 1
o
λ2 = J 2 θ 2
;
.......... ( 4.18 )
Reemplazando en la expresión para obtener T´, se tiene: T 1 (θ11,θ21 ) =
1 o2 1 o2 J θ 1 + J2 θ 2 2 1 2
Para el elemento elástico, se tiene: 1 2 V (θ3 ) = θ 2Kθ 3
......... ( 4.19 )
......... ( 4.20 )
Al hacer las uniones mecánicas tenemos: o
o
o
o
o
o
θ1 = θ a θa =θb
........ ( 4.21 ) o
o
o
θ 3 = θ a −θ b = θ 1−θ 2 En las expresiones de T´, V se obtiene: o
o
T 1 (θ a ,θ b ) =
1 o2 1 o2 J θ a + J2 θ b 2 1 2
........ ( 4.22 ) V (θa ,θb ) =
1 2 θa − θb ) ( 2Kθ
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
69
Fuerzas externas: Se analiza para un motor eléctrico conectado a una carga, siendo expresadas por: Motor: T1 = θa
........ ( 4.23 )
Carga: T2 = -θb
........ ( 4.24 )
Luego reemplazando en la ecuación restringida de Lagrange, tenemos: Para la coordenada θa: o0
J1 θ a +
1 (θ − θb ) = T1 (Torque electromagnetico) Kθ a
........ ( 4.25 )
Para la coordenada θb: o0
J2 θ b +
1 (θ − θb ) = −T2 ( torque de la carga) . Kθ a
........ ( 4.26 )
4.5 Determinación de la fuerza mecánica a través del principio de los trabajos virtuales.
La primera Ley de la Termodinámica, establece: o o o o o ∂Wn q λ dt + p x dt = ∆We + ∆Wm → {p = ∂x fuerza
o
λ=
T{
= M ( Momento)
Torque Electromagnetico
............ ( 4.27 ) T = [i][G][i] t Torque Electromagnetico T = Kφd i T = K i s i r sen γ 1 .......... ( 4.28 ) “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
70
además:
dWm ( λ , x ) = dWe (q, x ) =
o ∂ Wm ∂ Wm ∂ Wm . dλ + . dx = q dλ + . dx ∂λ ∂x ∂x
o ∂ We ∂ We ∂ We . dq + . dx = λ dq + . dx ∂q ∂x ∂x
......... (4.29 )
........( 4.30 )
Reemplazando ( 4.29 ) y ( 4.30 ) en ( 4.27 ), obtenemos: o
o
o
λ dq + p dx = λ dq +
o ∂ We ∂ Wm . dx + q dλ + . dx ∂ x ∂ x
o ∂ W e o dλ ∂ W m p = +q + { x ∂ ∂ x dx fuerza
........ (4.31 )
........ (4.32 )
La ecuación (4.32 ) representa la expresión general de la fuerza mecánica de origen electromagnética.
Fig. 4.2 Máquina simple reluctancia con una fase y dos polos en el estator y rotor
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
71
Fig. 4.3 Posición alineada
Fig. 4.4.Posición desalineada
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
72
Fig. 4.5 Variación de la inductancia y torque con la posición del rotor.
Fig. 4.6 Variación de la inductancia, corriente, flujo concatenado y fem con la posición del rotor:
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
73
Fig. 4.7 Cálculo del flujo distribuido por el análisis de elemento finito Para el caso de la influencia de elementos de almacenamiento de energía electromagnética, tenemos: o
p=
o ∂ Wm ∂ W m1 o ( λ, x ) → p = − (q, x ) ∂ x ∂ x
........... (4.33 ) o
W m + W m1 = λ q
..........( 4.34 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
74 4.6 Sistemas de excitaciones
Las máquinas eléctricas son sometidas a excitaciones del tipo mecánica y eléctrica, se denomina puertas al punto por donde sale o ingresa la energía. Se tienen las excitaciones del tipo: simple, doble y múltiple. 4.6.1 Excitación Simple
δ
i 0 o
N
λ
Wr
Fig. 4.8 Máquina eléctrica reluctante X ο
q K ο
λ
N
3,0 Cm 3,0 Cm
Fig. 4.9 Dispositivo electromecánico de excitación Simple.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
75 A continuación se desarrolla la expresión de la fuerza, para lo cual se establecen las ecuaciones siguientes: o
f = p = ∂ Wm ( λ,x)
........... ( 4.35 )
Tθ = Mθ = ∂ Wm ( λ,θ )
............ (4.36 )
o
o
o
dW m1 = λ d q → W m1 (q, x ) = ∫ λ d q
............ (4.37 )
De la configuración del circuito magnético de la figura 4.9, tenemos: o
∑ Hλ = Ni = N q
........... (4.38 )
o
H mλm + 2 H g λg = N q
Por continuidad de líneas de flujos magnéticos, se tiene: φm = φ g Bm Am = Bg Ag = u0 H g Ag A Hg = n u A 0 g
...........( 4.39 )
βm
Reemplazando la ecuación (4.39 ) en la ecuación ( 4.38 ), obtenemos: o An Bm = N q H m λ m + 2λ g ...........( 4.40 ) u 0 Ag
βm = u m H m → H m =
βm um
λ = Nφ = NBm Am → βm =
λ NAn
o λm An βm = N q + 2λ g u 0 Ag um
o
βm =
Nq λm An + 2λ g u 0 Ag um
“Estudio de las Máquinas Eléctricas. Simulación Digital”
........... ( 4.41 ) Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
76
luego: o
λ = NAn
Nq λm An + 2λ g u 0 Ag um
........... ( 4.42 ) o
λ=
An N 2 q λm An + 2λ g u 0 Ag um
→
λ o
q
=
An N 2 λm An + 2λ g u 0 Ag um
=L
.......... ( 4.43 )
Fig. 4.10 Característica λ – qº para diferentes longitudes de entrehierro Para un sistema lineal, se tiene: 1 o 2 dL ∂ Wm p= q = 2 dX ∂ X o
........... (4.44 )
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
77
Luego: o 1 o d An N 2 p= q 2 dX λ m 2 XAn u + u A 0 g m
............. ( 4.45 )
Para un sistema no lineal , se tiene la relación siguiente: H m = C0 β
m
+ C1 β
3 m
........... ( 4.46 )
Para determinar la fuerza, se deberá obtener la energía o coenergía del sistema y posteriormente derivar con respecto a la posición espacial.
Fig. 4.11 Lugar geométrico de operación del sistema a corriente constante
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
78
Fig. 4.12 Lugar geométrico de operación del sistema a flujo constante 4.6.2 Excitación doble
Fig. 4.13 Dispositivo electromecánico de excitación Doble “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
79
De acuerdo a la figura 4.13 se plantean las ecuaciones de equilibro electromagnético, siendo estas: Para el devanado estatorico “1”: o
o
λ 1 = V1 (t ) = R1 q1 (t ) +
dλ , λ = Nφ = ψ dt
............( 4.47 ) o
o
λ 1 = V1 (t ) = R1 q1 (t ) +
dψ 1 (t) dt
............ (4.48 ) Además:
ψ 1 = ψ 11 + ψ 12 o
o
ψ 1 = L1 q1 + M q 2 o o dλ1 d o = V1 (t ) = R1 q1 (t ) + L1 q1 + M q 2 dt dt
o
o
o dλ1 dq o dL d q2 o dM = V1 (t ) = R1 q1 (t ) + L1 1 + q1 wr 1 + M + q 2 wr dt dt dθ dt dθ .............. ( 4.49 )
Para el devanado rotórico “2”: o
o
o dλ 2 d q2 o dL d q1 o dM = V2 (t ) = R2 q 2 (t ) + L2 + q 2 wr 2 + M + q1 wr dt dt dθ dt dθ
............ (4.50 ) o
dW m = q dλ o
W m = ∫ q dλ
;
T=
∂ W dθ
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
80
Para la determinación del acoplo inductivo, de las ecuaciones anteriores, se tiene: o
o
o
λ 1 (t ) = ∫ R1 q 1 dt + ∫ L1d q 1 + ∫ q 1 W r
o o dL1 dM . dt + ∫ Md q 2 + ∫ q 2 W r dt dθ dθ
......... (4.51 ) La
característica de la inductancia puede ser: a. Constante b. K1 + K2 cos 2θ o
o
c. a + b q 1 + c q 12 Para el análisis de las máquinas eléctricas se presentan varios casos, siendo estos: Caso I.o
λ 1= VDC1
motor DC Excitacion Independ .
M (θ ) = M cosθ
o
λ 2 = VDC 2
L1 = K1
L2 = K 2
Caso II.o
o
q1 = q 2
,
Motor DC
M(θ ) = K 3 - K q cos2θ
VDC
Conexion seria
Caso III.o
λ 2= 0
,
o
λ 1(t) = λ
max
coswt ;
M (θ ) = a + b cos2θ
Motor de induccion
L1 , L 2 constante
Para el análisis de la velocidad, se tiene: − cons tan te wr = − Variable Modulo Sentido “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
81 Caso IV.O
λ
2
O
λ
1
= VDC − Tension = − Frecuencia
V = f ( w, β ,θ ) V , f f = f (w ) r { Frec .
Donde la velocidad puede ser: − wr > ws wr r − w < ws
; ws =
w P
; w = 2Π f
− Cons tan te wr − Variable Caso V.O
λ 2= λ
max
cos wt
− wr > ws w r − w < ws r
V = f ( w, β ,θ ) V, f r f resul tan te = f ( f , f c ) r f r = f c ± f mec
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
CAPITULO V EL CAMPO MAGNETICO EN EL ENTREHIERRO
El objetivo del presente capítulo es plantear las ideas básicas necesarias para abordar en profundidad el estudio de las máquinas eléctricas de corriente alterna. Una de la funciones principales de los devanados es la creación de un campo magnético cuya distribución espacial en el entrehierro sea sinusoidal (o muy próxima a la forma senoidal). Los devanados cumplen esta misión mediante una construcción y distribución adecuada de las bobinas a lo largo de la periferia del estator y rotor, la alimentación de dichos devanados con corriente polifásica contribuyen a conseguir la aproximación a la forma de onda senoidal. La polifasidad ejerce dos funciones principales como son: la creación de un campo magnético giratorio y la supresión adicional de armónicos de la onda de campo.
5.1 Campo magnético en una máquina eléctrica ideal. Se establece una máquina eléctrica cuyo estator lo constituye un cilindro perfecto hueco a base de chapas magnéticas en cuyo interior concéntricamente y con análoga constitución se encuentra alojado el rotor. Se plantean las hipótesis siguientes: 1.- La permeabilidad magnética del hierro( μm ) es muy grande comparada con la del vacío ( μo ). Se desprecian las pérdidas en el hierro. 2.- La dimensión del entrehierro( δ ) , que es simétrico y constante es muy pequeña comparada con el diámetro del estator( Ds ) y rotor( Dr ) . 3.- Los conductores están distribuidos en forma simétrica sobre la superficie cilíndrica exterior del rotor e interior del estator, dispuestos paralelamente al eje de la máquina y sus dimensiones radiales se consideran despreciables (conductores puntuales). 4.- En un primer análisis no se considera los fenómenos de una eventual saturación magnética. No se toma en cuenta la dispersión del campo de los bordes laterales de la máquina. En la Fig. 5.1, se muestra una sección transversal de la máquina eléctrica ideal donde se ha dispuesto arbitrariamente que existen seis conductores en el estator y rotor.
a’ b
c
c’
b’ a
Fig. 5.1 Distribución
Los conductores del estator pueden estar por ejemplo, conectados en serie entre sí y exteriormente a una sola fase de tensión continua o alterna, o bien estar cada dos conductores opuestos diametralmente unidos entre sí, resultando tres fases, o bien adoptar otra forma de interconexión diferentes de las anteriores. Similar consideraciones se pueden hacer para los conductores del rotor.
5.1.1 Fuerza magnetomotriz creada por una espira simple de paso diametral. Supongamos que no existe devanado en el rotor (o bien que se encuentre abierto) y, que el devanado del estator está conformado por una espira simple diametral aa' según fig. 5.2 por la cual circula una corriente continua I.
A
D
α M
a´
B Q
C R
P
S
N
T
a V
Fig. 5.2
Debido al supuesto que la permeabilidad en el hierro( μm ) es infinita, las líneas del campo magnético en el entrehierro son radiales y el campo magnético en el hierro será nulo (es decir, toda la tensión magnética cae en el entrehierro). La primera ley de Maxwell establece:
r r r ∂D rotH = J + ∂t
.......(5.1)
Para corriente continua ó para bajas frecuencias, tenemos: r ∂D ≈0 ∂t Luego, la ecuación (5.1), podemos escribirlo: r r r r rotH * dA = J * dA
…..(5.2)
Aplicando el teorema de Stokes, tenemos:
r r H * dl = ∑ I
…..(5.3)
ΣI representa la corriente encerrada por la superficie a la cual se ha extendido la integral. La ecuación (5.3) es conocida también como el teorema de Ampere. r r H ∫ . dl = ∫ H AB * dl + ∫ H DC * dl =0 ABCD
AB
CD
......(5.4)
Debido a las dimensiones muy pequeñas del entrehierro, se considera constante el valor del campo magnético en la dirección radial del mismo siendo δ la anchura del entrehierro, de la ecuación (5.4) podemos establecer: HAB * δ = -HCD * δ = HCD * δ HAB = HDC Por consiguiente, la distribución del campo magnético en el entrehierro es uniforme, siendo por simetría de valor opuesto en la zona N ( Flujo entrante) y S (flujo saliente). N H N = - HN
S En la Fig. 5.3, se representa mediante sus líneas de fuerzas correspondientes a la g. 5.3 distribución del campo magnético en laFimáquina creado por una espira de paso Distribución
diametral. Se observa que las líneas de campo en el entrehierro muestran una distribución uniforme a lo largo de él, característica que traduce gráficamente la constancia del campo en el entrehierro. Aplicando la ecuación (5.3) a las trayectorias cerradas MNPQ y RSTV de la Fig. 5.2, tenemos: I f N = HN * δ = 2 f S = HS * δ = −
I 2
......(5.5) El valor H * δ en un punto de coordenada angular (fig. 5.2) recibe el nombre de tensión magnética ó fuerza magnetomotriz (f.m.m.) de entrehierro en dicho punto, y esta representada por fα.. Si desarrollamos el entrehierro de la máquina en forma lineal y se representa en función del ángulo α la f.m.m. de entrehierro, obtenemos la figura 5.4 f
I/2
a
a’
0
π
α 2π
- I/2
Fig. 5.4
Si en lugar de una espira elemental, consideramos una bobina elemental con Ns espiras supuestas idealmente de dimensiones infinitesimales, tendríamos:
fN =
NS * I 2
;
fS = −
NS * I 2
....(5.6)
5.1.2 Fuerza magnetomotriz de una bobina de paso recortado. Se considera una espira elemental según mostrada en la Fig. 5.5, cuyos lados no se hallan situados sobre el mismo diámetro.
α
Aplicando el mismo razonamiento que del apartado anterior (5.1.1), se demuestra que el campo en el entrehierro a lo largo del arco aNa' posee en todas las partes el mismo valor. Simultáneamente ocurre para el arco aSa'. La tercera ley de Maxwell establece: r Div B = 0
......(5.7)
La ecuación (5.7) expresa que no existen cargas magnéticas libres, y por lo tanto el flujo magnético es conservativo. Por consiguiente, los módulos de las fuerzas magnetomotrices norte y sur tendrán que cumplir la relación:
f N arco aNa ′ = fS arco aSa ′
..... (5.8)
La relación (5.8), traduce matemáticamente el hecho físico de entrantes y salientes en el estator (fig. 5.5) son iguales.
que
los
flujos
Asimismo, aplicando la ecuación (5.3) a dos superficies que contengan a cada una un conductor de la espira, obtenemos: fN + fS = I Si en lugar de una espira, se considera una bobina elemental con Ne espiras, tenemos: fN + fS = Ne * I
f
a
a’
α
En la fig. 5.6 se muestra la representación de la f.m.m. a lo largo del entrehierro. Como podemos observar en la fig. 5.5 se han dibujado más próximas las líneas de fuerza del campo en la zona norte, lo que significa que la densidad de líneas de fuerza (campo de f.m.m. de entrehierro) sean mayores en dicha zona. Consideremos ahora una segunda espira, análoga a la anterior pero desplazada un ángulo 2β ; según mostrada en la fig. 5.7.
a
a’ 2β
b
b’
Fig. 5.7
Las formas de ondas de la f.m.m. de entrehierro debidas a cada espira así como su resultante aparecen en la fig. 5.8. Podemos observar que la onda resultante tiene la tendencia ha aproximarse a la forma teórica senoidal buscada (posee menor contenido de armónicos que las ondas representadas en las fig. 5.4 y fig. 5.6.
f
α
5.1.3 Fuerza magnetomotriz de una bobina múltiple. Existe la posibilidad de repartir las NS espiras de la bobina elemental en una serie de 'q' bobinas elementales distribuidas a lo largo de la periferia conectada en serie entre sí y formada cada una de ellas a su vez por Ne espiras, constituyendo una bobina múltiple. En la figura 5.9 está representada una distribución de bobinas múltiples formadas por 5 bobinas elementales diametrales.
a’
e γe
b’
d
c’
c
d’
b e’
a
Fig. 5.9
La onda resultante del f.m.m. de entrehierro puede obtenerse directamente superponiendo las ondas correspondientes a cada bobina elemental, y se muestra en la Fig. 5.10.. Es apreciable la aproximación a la forma senoidal que se obtiene mediante este procedimiento. Se puede lograr una mayor aproximación, distribuyendo el devanado en un número superior de ranuras y empleando bobinas de paso acortado adecuado (asimismo, la
polifasidad contribuye a eliminar nuevos armónicos). En el caso general de una máquina de 'p' pares de polos con 'q' bobinas elementales en serie por par de polos y con Ne espiras cada una, el valor de la f.m.m. máxima es: fm = q*
Ne * I NS * I = 2 2* p
......(5.9)
Siendo Ns el número total de espiras conectadas en serie.
qNeI 2 a
b
c
d
e
a’
b’
c’
d’
e’
NeI 2
Fig. 5.10
5.1.4 Análisis de armónico. Las curvas obtenidas relativas a la distribución espacial de la f.m.m. en el entrehierro, representa a una cierta escala y con los supuestos admitidos en el estudio de la máquina ideal, las curvas de campo y de inducción en el entrehierro. Considerando la ecuación 5.5, tenemos:
Bα = μ0 * Hα = μ0 *
fα
δ
.....(5.10)
Siendo μ0 yδ (entrehierro) constantes, queda demostrada la semejanza de las curvas.
H
ONDA CUADRADA FUNDAMENTAL
3 er ARMONICO 5 to ARMONICO
En la figura 5.11 se muestra la curva de distribución de campo en el entrehierro creada por una única bobina de paso diametral de Ns espiras, siendo el valor máximo, Hm : Hm =
fm
δ
I * NS 2* p*δ
=
.....(5.11)
Si la intensidad I varia siguiendo una función del tiempo I = I(t), la configuración de la curva ser la misma variando únicamente la altura del rectángulo al ritmo de las variaciones de I, siendo: H m (t) =
Ns * I( t ) 2* p*δ
.....(5.12)
El desarrollo en serie de Fourier de una onda rectangular, tomando como origen del ángulos el origen del rectángulo, tenemos: 4 .....(5.13) Hα = * H m ( t )* [ sen α + 1 / 3 sen 3α + 1 / 5 sen 5α + ...] Π La expresión de la onda fundamental esta dada por : H1 =
4
Π
* H m ( t )* sen α =
4
Π
*
Ns * I ( t )* sen α 2* p* δ
El valor máximo (en el espacio, para un instante determinado) es: H m1 ( t ) =
4
Π
* Hm ( t )
.....(5.15)
......(5.14)
El armónico de orden 'h' está expresado por: Hh ( t ) =
4
Π
* H m ( t )*
1 * sen( hα ) h
.....(5.16)
Su valor máximo es: H mh ( t ) =
1 4 1 * * H m ( t ) = * H m1 ( t ) h Π h
.....(5.17)
Es conveniente mencionar: a) No existen Armónicos pares. b) La amplitud de los armónicos es inversamente proporcional al orden del armónico.
5.2 Factores de Paso y de Distribución. 5.2.1 Factor de Paso Si consideramos el caso de dos bobinas de paso acortado conectados en serie, se obtiene la curva mostrada en la figura 5.8. El desarrollo en serie de Fourier de dicha onda rectangular, eligiendo como origen de ángulos el punto intermedio entre los rectángulos esta expresado por: ∞
H(α ,β ) = ∑ h= 1
4
Π
*
1 * H m * cos( hβ )* sen( hα ) h
Desarrollando, tenemos: H(α ,β ) =
1 ⎡ ⎤ H m * ⎢ sen α * cos β + * sen 3α * cos 3β + .....⎥ Π 3 ⎣ ⎦ 4
..... (5.18)
El ángulo 2β es el acortamiento de paso expresado en magnitudes angulares eléctricas. Generalmente el acortamiento se da por una función (1/3,1/4,..., en general 1/k) del paso polar, como el paso expresado en ángulo eléctrico corresponde a Π radianes, tenemos: 2β =
Π k
;
β=
Π 2* k
Luego, la ecuación (5.18), podemos expresarla: Hα =
Π 1 3Π 1 5Π ⎛ ⎞ H m ⎜ sen α cos + sen 3α .cos + sen 5α cos + .....⎟ .....(5.19) ⎝ ⎠ Π 2k 3 2k 5 2k 4
Se define el factor de paso ( kp ) como la relación que existe entre el valor máximo de la onda fundamental o armónico producido por la bobina de paso acortado, y el que correspondería si las bobinas fueran de paso diametral. De las ecuaciones (5.13) y (5.19), se deduce para la componente fundamental: K p1 = cos β = cos
Π 2k
.....(5.20)
Y para el armónico de orden 'h', tenemos: ⎡h* Π ⎤ K ph = cos( hβ ) = cos⎢ ⎣ 2 * k ⎥⎦
.....(5.21)
La ecuación (5.21) expresa que cualquier armónico 'h' puede eliminarse de la curva de inducción ó de la curva de campo acortando el paso en una fracción h=k. Podemos observar que el acortamiento de paso que suprime por completo un armónico, no influye demasiado a la componente fundamental, aunque el armónico suprimido sea de orden bajo. Por ejemplo: Para : h = k = 5 , resulta: kp1 = 0.951 h = k = 7 , resulta: kp1 = 0.975
5.2.2 Factor de Distribución Para este caso consideramos una bobina múltiple formada por varias bobinas elementales distribuidas sobre la periferia (según fig. 5.9). Se define el factor de distribución (Kd) como la relación entre el módulo del vector resultante y la suma aritmética de los módulos de los vectores componentes. En otra expresión, como la relación entre la amplitud de la onda obtenida con un devanado distribuido y aquella que se obtendrá si todas las bobinas elementales estuvieran concentradas formando una única bobina diametral. Sea 'Q' el número total de ranuras por polo, si la máquina tiene 2p polos (p: número de pares de polos) el número total de ranuras Nr será: Nr = 2 * p * Q El ángulo geométrico τg formado por dos ranuras, está expresado por:
τg =
2*Π 2*p*Q
Y el ángulo eléctrico τe, es:
τe = p* τg =
Π Q
.....(5.22)
En el caso más general de un devanado de 'm' fases, se da como dato el valor 'q' que representa el número de ranura por polo y fase. Podemos establecer: Q=q*m
....(5.23)
Dividiendo una circunferencia en 2Q partes iguales (según fig. 5.12). Cada lado del polígono inscrito AB , BC , etc. representa en magnitud y fase el vector representativo de la onda fundamental de campo debido a cada una de las bobinas elementales alojadas en la ranuras 1, 2, ...., q correspondiente a una bobina múltiple de una fase. E D C B
A
qγe
r
γe
Fig. 5.12
El vector de campo resultante será: AE = AB + BC + CD + DE El factor de Distribución de acuerdo a su definición es:
Kd =
2 * r * sen( q * τ e / 2 ) AE = q * AB q * 2 * r * sen( τ e / 2 )
Kd =
sen( q * τ e / 2 ) q * sen( τ e / 2 ) ...........(5.24)
En el caso ideal, de que el número de bobinas elementales fuera muy grande podrá admitirse teóricamente que están distribuidas ocupando todo el arco. Luego el valor límite de Kd, es: ( K d )limite =
( K d )lim . =
cuerda 2 * r * sen( q * τ e / 2 ) = arco r * (q * τe )
sen( q * τ e / 2 ) (q * τe / 2)
.......(5.25) El factor de distribución Kdh correspondiente a un armónico cualquiera de orden h, se expresa:
K dh =
sen( q * h * τ e / 2 ) q * sen( h * τ e / 2 )
.....(5.26)
Para el caso de bobinados trifásicos, se tiene Q = 3 * q. De la fórmula (5.22) se tiene:
τe =
Π 3* q
Podemos considerar que la distribución del devanado actúa como un filtro pasa banda que disminuye la amplitud de todos los armónicos, mientras que el acortamiento de paso equivale a un filtro selectivo que elimina un armónico particular.
5.2.3 Factor de Devanado. Devanado Eléctrico Equivalente En el caso general de un devanado con bobinas distribuidas y de paso acortado, el valor de la onda fundamental ó de la onda de cualquier armónico, se calcula a partir del valor de la onda que se obtiene con una bobina única de paso diametral que tenga la totalidad de las espiras concentradas, multiplicando por un factor igual al producto Kp * Kd, para la onda fundamental y Kph * Kdh para el armónico de orden h. A este factor se le denomina factor de devanado K,Kh respectivamente. K = Kp * Kd Kh = kph * Kdh Por consiguiente, a efectos de creación de ondas de f.m.m. o de campo, una bobina de paso y distribución cualquiera con un número real q * Ne de espiras en serie por par de polos y fase, puede sustituirse por una única bobina concentrada de paso diametral y con un número efectivo de espiras por par de polos y fase Nef, variable para cada armónico, siendo: Nef = (q * Ne) * K
Nefh = (q * Ne) * Kh
....(5.27)
Luego, de las ecuaciones (5.12),(5.17) y (5.27) se deduce que el valor máximo del campo en el entrehierro correspondiente al armónico de orden h, producido por un devanado distribuido y de paso recortado, para un instante cualquiera es: H mh ( t ) =
1 4 N S * Kh * I ( t ) * * h Π 2* p*δ
Si se considera la componente fundamental, su valor máximo es: H m1 ( t ) =
4
*
Π
NS * K * I(t ) 2* p*δ
Como podemos observar, el valor máximo de esta onda así como el de cualquier armónico no es constante, variando en el tiempo de acuerdo con las variaciones de I(t). Para el caso práctico de que la corriente a través de la fase sea senoidal, tenemos: I(t) = √2 * Ief * sen (w1 * t) Reemplazando I(t) en la expresión de Hm1(t), tenemos: H m1 ( t ) =
H m1 ( t ) =
4
Π* 2
*
N S * K * I ef p*δ
0 .9 * K * N S * I ef p*δ
* sen( w1 * t )
* sen( w1 t ) ......(5.28)
La ecuación (5.28) proporciona el valor de la amplitud en todo instante de la onda fundamental del campo magnético en el entrehierro, creada por una fase alimentada por corriente alterna sinusoidal. 5.3 Onda de Campo en el entrehierro senoidal pura. Devanado Eléctrico equivalente ideal. Se busca la expresión matemática del devanado eléctrico equivalente a aquel que genera un campo senoidal puro en el entrehierro, introduciéndose para ello el concepto de capa de espiras que es definida como la densidad de conductores por unidad de longitud periférica por los que circula la corriente en sentido positivo. Se demuestra, que si la capa de espiras presenta una distribución senoidal tal como aparece en la figura 5.13, se consigue una onda de f.m.m. de entrehierro también senoidal pura.
M
Q
N
P
α
Debido a que todas las espiras están recorridas por la misma corriente, la densidad lineal de corriente a lo largo de la periferia (llamada también capa de corriente) seguir también una ley senoidal representada por : Aα = Am * senα Siendo Aα el valor de la densidad lineal de corriente en el punto de coordenada angular α y, Am el valor máximo. Por simetría, el campo magnético en los puntos α y Π-α tiene el mismo módulo pero de sentido opuesto. Aplicándose el teorema de Ampere a la superficie cerrada MNPQ (fig. 5.13), podemos expresar:
r
r
Π −α
∫ H * dl = ∫α A
m
* r * sen α * dα
Resolviendo, tenemos : − Hα * MN − Hα * PQ = −2 * Hα * δ = 2 Am * r * cos α Hα * δ = − Am * r * cos α
Π⎞ ⎛ fα = Am * r * sen⎜ α − ⎟ ⎝ 2⎠ ....(5.29) La ecuación (5.29) expresa que, la curva de f.m.m. de entrehierro es también senoidal pura y que se encuentra retrasada 90° en el espacio con respecto a la onda de la capa de corriente(según fig.5.14). Así mismo podemos observar en la fig. 5.13, que la densidad de las líneas de campo varia sinusoidalmente a lo largo del entrehierro, presentando su valor máximo positivo con un retraso de Π/2 radianes con respecto al punto de máxima densidad de conductores positivos.
f (α) = fm.sen (α - π/2)
Por consiguiente, todo devanado formado por un conjunto de espiras conectadas en serie construidas y distribuidas de modo cualquiera, pero tal que produzca una onda f.m.m. de entrehierro senoidal pura puede ser sustituido a efectos de f.m.m.‚ inducción por otro cuya capa de espiras supuestas diametrales siga una ley de distribución sinusoidal. En la práctica por razones constructivas, no se realiza físicamente un devanado como en la fig. 5.13 sino que, en general todas las ranuras contienen igual número de conductores. La forma senoidal del campo se consigue por otros procedimientos, habiéndose desarrollado los más importantes para máquinas de entrehierro constantes en los capítulos anteriores. Con un sistema polifásico de infinitas fases (en la práctica muchas fases) formada cada una de ellas por una sola espira diametral se conseguiría una curva de f.m.m. senoidal pura. Las ventajas de trabajar con capas de espiras en lugar de capas de corrientes se lograrán al efectuar el análisis dinámico y estático de las máquinas eléctricas de corriente alterna polifásicas mediante vectores espaciales.
5.4 Campos Magnéticos Giratorios 5.4.1 Obtención mediante un elemento móvil Debido a un solo devanado formado por un conjunto de espiras en serie, que se agrupan formando diversas bobinas elementales puede lograrse un campo en el entrehierro cuya configuración espacial se asemeja a una senoide, según desarrollado en el apartado 5.1 si la corriente que alimenta el devanado es constante, el campo no cambia; si se alimenta con una corriente que varía en el tiempo siguiendo una ley determinada (se suponen variaciones cuya rapidez esté comprendida en el campo de las frecuencias industriales), el módulo del campo variará al mismo ritmo que la corriente, pero su configuración espacial (senoidal) permanece inalterable. Es decir, resulta un campo de módulo variable en el tiempo pero de eje fijo en el espacio. En la máquinas de corriente alterna de campo giratorio se desea obtener lo contrario; es decir un campo de distribución senoidal de módulo constante en el tiempo, pero de eje variable en el tiempo (campo giratorio).
Una primera forma de conseguirlo consiste en desplazar físicamente mediante un acoplamiento mecánico exterior, el estator, por cuyo devanado se hace circular una corriente continua. Este procedimiento se emplea en los generadores síncronos, con la única diferencia de hacer girar exteriormente el rotor en lugar del estator. Supongamos el devanado del estator de una máquina constituido por un conjunto de espiras conectadas en serie, alimentadas con corriente continua y distribuidas de tal manera que dan origen a un campo senoidal en el entrehierro de 'p' pares de polos. Manteniendo la corriente constante, se hace girar al estator con una velocidad angular Wmec. constante en sentido positivo; al cabo de un tiempo t, la onda senoidal de campo en el entrehierro que inicialmente ocupaba la posición '1' en la figura 5.15 estará representada por la curva de posición '2'.
H ‘2’
‘1’
0
β
P=2
2π rad. geom.
En un punto 'p' cualquiera, el valor del campo en el instante t=0, es: Fig. 5.15
H(0,α) = Hm * cosα En un instante t, posterior, el ángulo eléctrico recorrido por la onda de campo en su giro será p * wmec * t. Luego, el valor del campo en el mismo punto 'p' para el instante considerado es: H(t,α) = Hm * cos(α - p * wmec * t) = Hm * cos(p * wmec t - α) H(t,α) = Hm * cos(w1 * t - α)
.....(5.30)
Siendo: w1 = p * wmec, velocidad angular eléctrica. Si el giro del estator es de sentido contrario, la expresión del campo es: H(t,α) = Hm * cos(w1 * t + α)
.....(5.31)
Generalizando, podemos establecer: H(t, α ) = Hm cos(w1t + α)
....(5.32)
La ecuación (5.32), constituye la expresión general de un campo giratorio de velocidad angular eléctrica constante w1 (correspondiente a una velocidad mecánica w1/p) y de modulo constante en el tiempo.
5.4.2 Obtención mediante procedimientos estáticos. Estudio de la Componente fundamental. Considerando un devanado trifásico colocado en el estator y alimentado por un sistema de corrientes trifásicas senoidales. En la figura 5.16 cada fase esta representada por una bobina concentrada de paso diametral; en estas condiciones cada fase originará un campo magnético en el entrehierro de configuración espacial pero de amplitud variable en el tiempo, el campo real en todo instante es una resultante originada debido a la interacción de los tres campos en cada momento.
R´
T
S
S´
T´
α
En un primer análisis consideramos las componentes fundamentales, con la suposición que los armónicos de cada fase resultanRdespreciables, y que éstas generan campos senoidales puros. Fig. 5.16 Se elige como origen de tiempo el instante en que la intensidad por la fase 'R' alcanza su valor máximo. En un instante 't' posterior, los valores máximos del campo para cada fase son: HRm(t) = Hm * cos(w1 * t) HSm(t) = Hm * cos(W1 * t – 2Π/3) HTm(t) = Hm * cos(w1 * t - 4Π/3) Hm corresponde al 'maximun maximorum' es decir, al valor máximo posible que puede alcanzar el campo creado por una fase y que se presenta justamente cuando la corriente de dicha fase pasa por el máximo. Las expresiones HRm, HSm y HTm representan los valores máximos instantáneos (cambiantes con el tiempo) y que, en el espacio se presentan en aquellos puntos fijos del entrehierro situados en los ejes de la fase correspondiente. Como dicho puntos tienen una desfase entre si de 2Π/3 radianes eléctricos en el espacio, el valor del campo magnético creado por cada una de las fases en un punto del entrehierro distante
en ángulo 'α' del origen de espacios (que en la fig. 2.16, para simplificar las expresiones matemáticas se ha hecho coincidir con el eje de la bobina de la fase 'R') en el instante 't' es: HR(t, α)= HRm * cos α = Hm * cos(w1 * t) * cosα HS(t, α)= HSm * cos(α - 2Π/3) = Hm * cos(w1t - 2Π/3) * cos(α - 2Π/3) HS(t, α)= HTm * cos(α - 4Π/3) = Hm * cos(w1t - 4Π/3) * cos(α - 4Π/3) .....(5.33) Puesto que el campo en los tres casos es radial, el campo resultante H(t, α) ser también radial y de valor igual a la suma aritmética de HR(t, α), HS(t, α) y HT(t, α); es decir: H(t, α) = HR(t, α) + HS(t, α) + HT(t, α) Reemplazando la ecuación (2.33), tenemos: ⎡ cos( w1 * t ) cos α + cos( w1 * t − 2Π / 3 )* cos( α − 2Π / 3 ) + ⎤ H(t,α ) = H m ⎢ ⎥ ⎣ cos( w1 * t − 4Π / 3 )* cos( α − 4Π / 3 ) ⎦
Aplicando identidades trigonométricas, la expresión H(t, α) resulta: H(t,α ) =
Hm 2
⎡ 3 * cos( w1 * t − α ) + cos( w1 * t + α ) + cos( w1 * t + α − 2Π / 3 ) + ⎤ ⎢ cos( w * t + α − 4Π / 3 ) ⎥ 1 ⎣ ⎦
Como podemos observar, los tres últimos sumandos del corchete representan tres vectores de módulo igual a la unidad y desfasados 2Π/3 rad., cuya suma es nula. Luego podemos expresar: H(t,α ) =
3 * H m * cos( w1 * t − α ) 2
.....(5.34)
La ecuación (5.34) coincide con la (5.30) y representa la expresión matemática de un campo giratorio de módulo constante y de velocidad angular igual a w1 radianes eléctricas por segundo. Los efectos magnéticos del sistema trifásico analizado se traducen en la creación de un campo magnético en el entrehierro cuya configuración espacial es una senoide, dicho campo gira físicamente en tono al eje axial de la máquina manteniendo constante su forma, con una velocidad angular geométrica de w1/p radianes por segundo. Desde el punto de vista del campo magnético nos ubicamos exactamente en la misma situación
que la analiza en el apartado 5.4.1, pero sin necesidad de efectuar el giro mecánico del rotor o del estator. En el caso general de un sistema de 'm' fases, operando de la manera análoga se obtiene la expresión siguiente: H(t,α ) =
m H * cos( w1 * t m α ) 2 m
.....(5.35)
5.4.3 Obtención mediante sistema polifásico. Análisis de los armónicos. Realmente, la onda de campo ó de f.m.m. de entrehierro producida por cada fase no es una senoide pura, la descomposición en serie de Fourier de dicha onda conduce a una componente fundamental y una serie de ondas armónicas senoidales. De la misma manera que lo analizado para la onda fundamental, los armónicos de un mismo orden correspondiente a las tres fases y cuya resultante no sea nula, se combinarán entre sí originando campos armónicos giratorios. Para el análisis, se consideran como origen de escalas y espacios los mismos que los señalados en el figura 5.16, siendo 'α' un ángulo eléctrico referido a la componente fundamental. Como origen de tiempo se fija el instante en que la intensidad por la fase 'R' representa su valor máximo. En un instante 't' cualquiera el valor del campo correspondiente al armónico de orden 'h' creado por la fase 'R' en su eje (que coinciden con el origen de coordenadas) es: HRh(t) = Hmh * cos(w1 * t - 2Π/3)
....(5.36)
siendo:
Hm h Hmh, representa el valor máximo correspondiente a dicho armónico. H mh =
Similarmente, el valor del campo del armónico de orden 'h' creado por las fases S y T en los ejes de sus bobinas correspondientes, para un instante 't' cualquiera, están expresados por: HSh(t) = Hmh * cos(w1 * t - 2Π/3)
....(5.37)
HTh(t) = Hmh * cos(w1 * t - 4Π/3)
.....(5.38)
Además, el punto de entrehierro dado por la coordenada angular 'α' corresponde a un ángulo 'hα' referido al armónico h. Por lo tanto el valor del campo en dicho punto correspondiente al armónico de orden h, debido a la fase 'R' es: HRh(t, α) = Hmh * cos(w1 * t) * cos(h * α)
.....(5.39)
El campo debido a las fases S y T en dicho punto y para el mismo armónico son: HSh(t, α) = Hmh * cos(w1 * t - 2Π/3) * cos[(α - 2Π/3) * h]
....(5.40)
HTh(t, α) = Hmh * cos(w1 * t - 4Π/3) * cos[(α - 4Π/3) * h
.....(5.41)
Debido a que el campo en los tres casos es radial, el campo resultante será también radial y de valor igual a la suma aritmética de los componentes de cada fase. Sumando las ecuaciones (5.39), (5.40) y (5.41), tenemos: 2Π ⎞ 2Π ⎞ ⎤ ⎡ ⎛ ⎛ ⎢ cos( w1 t ) + cos( αh ) + cos⎜⎝ w1 t − 3 ⎟⎠ cos⎜⎝ α − 3 ⎟⎠ h + ⎥ ⎥ H RST,h (t,α ) = H m ⎢ ⎢ ⎛ ⎥ 4Π ⎞ 4Π ⎞ ⎛ ⎢ cos⎜⎝ w1 t − 3 ⎟⎠ cos⎜⎝ α − 3 ⎟⎠ h ⎥ ⎣ ⎦ Efectuando operaciones y simplificando, obtenemos las siguientes expresiones: 1.- Para h=3, ó múltiplo de 3; es decir h = 3 * K, tenemos: HRST, 3K (t,α) = 0 .....(5.42) 2.- Para h = 5, en general para todos los armónicos de orden h = 6 * K - 1, tenemos: HRST, (6K - 1) (t, α) = 3/2 * Hm(6K - 1) * cos[w1 * t + (6 * K - 1) * α] ....(2.43) 3.- Para h = 7, en general para todos los armónicos de orden h = 6 * K + 1, tenemos: HRST, (6K - 1) (t, α) = 3/2 * Hm(6K + 1) * cos[w1 * t - (6 * K + 1) * α] ....(2.44) De las ecuaciones (5.42),(5.43) y (5.44), podemos establecer las siguientes conclusiones: a.- Los campos armónicos de la fases individuales de orden (6K + 1) de eje fijo en el espacio y de módulo variable senoidalmente en el tiempo, se combinan dando lugar a un campo magnético giratorio armónico de orden (6K + 1) que gira en sentido contrario que el campo fundamentalmente y con una velocidad geométrica igual a ns/(6K + 1). El valor del módulo del campo es: H m ( 6 K +1) =
1 * Hm (6 * K + 1)
Siendo Hm, el módulo del campo giratorio fundamental.
b.- Los campos armónicos de las fases individuales de orden (6K - 1), de eje fijo en el espacio y módulo variable senoidalmente en el tiempo de combinan dando lugar a un campo magnético giratorio armónico de orden (6K - 1) que gira en el mismo sentido que en el campo fundamental y con una velocidad geométrica igual a ns / (6K - 1). El valor del módulo de dicho campo es: H m ( 6 K −1) =
1 * Hm (6 * K − 1)
c.- Los campos armónicos de las fases individuales de orden 3 * K, quedan eliminados. La conclusión 'c' permite comprender como un sistema trifásico (en general, polifásico) suprime adicionalmente armónicos no eliminados por las fases individuales, acercando la forma de onda del campo aún más a la senoide, por esta razón los devanados de las máquinas reales tienden a un acortamiento reducido, y al armónico restante de orden más bajo es el 11, cuya amplitud es muy pequeña.
5.5 Campo Magnético en el entrehierro creado por el devanado del rotor. Los análisis de los capítulos anteriores, donde se consideró a los conductores localizados en el estator y el devanado del rotor en circuito abierto, son aplicables para el caso opuesto (circuito de estator abierto e corriente circulante por las bobinas del rotor, es decir cualquiera que sea la disposición de los conductores o bobinas del rotor, cuando por ellos circula una corriente, se produce en el entrehierro de la máquina el mismo campo que el producido por un devanado idéntico situado en el estator y atravesado por la misma corriente, puesto que las dimensiones radiales de los conductores se han supuesto despreciables, así mismo las dimensiones del entrehierro frente al diámetro del rotor.
5.6 Campo magnético en el entrehierro creado por la acción conjunta de los devanados del rotor y estator en el caso de una máquina eléctrica ideal. En el caso teórico ideal, los devanados del rotor y estator producen individualmente dos ondas del campo de entrehierro perfectamente senoidales, con independencia de que las corrientes que atraviesan ambos devanados sean diferentes. En el caso teórico de devanados cuya capas de espiras sea senoidal (caso ideal equivalente al cual se aproximan los devanados reales) las formas de los campos de entrehierro de rotor y estator no dependen del valor instantáneo de la corriente que atraviesa los devanados, sino que presenta siempre una configuración senoidal; lo que si dependen de las corrientes es el valor efectivo de dichas ondas de campo. Además, como sucede siempre en la práctica el número de pares de polos de las ondas de campo correspondientes a los bobinados de rotor y estator es lo mismo, dichas ondas senoidales se suman produciendo en todo instante una onda de campo senoidal, y al suponer la ausencia de saturación, la onda de inducción de entrehierro de la máquina ideal presenta también en estas condiciones una configuración senoidal.
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
172
CAPITULO VIII MAQUINAS ELECTRICAS I EXAMEN FINAL Prob. 1: Se tiene un contactor según mostrado en la figura adjunta, cuyo núcleo ferromagnético esta conformado por 60 laminas de espesor t = 0,5 mm c/u. Cuando las secciones del núcleo son mantenidas juntas el entrehierro inherente es de 1,0 mm. La bobina consta de 250 espiras y de una resistencia de 7,5 ohmios medida en DC. Considerar: Um >> U0, RAC = 1,2 RDC. X ο
q K ο
λ
N
3,0 Cm 3,0 Cm
1.1 Si la bobina es excesitada con una fuente de tensión continua, tal que, para mantener las dos secciones juntas desarrolla una fuerza de 333,25Kg.f. 1.2a) La tensión aplicada. 1.1b) La corriente que absorbe el contactor. 1.1c) La energía magnética.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
173 1.2 Si al contactor aplicamos una tensión alterna a una frecuencia de 60 Hz, tal que, para mantener las dos secciones del núcleo juntas desarrolla una fuerza de 69.51 Kg.f (prom.) 1.2a) La energía magnética. 1.2b) La tensión aplicada. Prob. 2: Una maquina de inducción 3φ, según mostrado en la figura adjunta, de rotor bobinado, conectada en estrella, 2 polos, 3kW, 220V, 60Hz, cuando esta parada tiene una tensión de terminal de línea para el rotor en circuito abierto, que es dos veces la tensión aplicada entre líneas al estator. Con el circuito del rotor abierto, el eje de éste gira a 1600 rpm, en el sentido del campo giratorio. Se aplican al estator la tensión y frecuencia nominales.
Anillos Pac deslizantes
Fuente 3φ
+ VDC Pdc
n P2
Maquina
Maquina
de Inducc.
DC
Frecuencia Variable
2.1a) Hallar las expresiones convenientes para las tensiones que, en circuito abierto, aparecen en cada par de los anillos de deslizamiento. 2.1b) A qué velocidad tendría que ser impulsado el rotor, para suministrar una frecuencia de 80 Hz? 2.1c) Si el rotor se mueve a 1800 rpm. Hallar la magnitud y frecuencia de la tensión inducida en él, para los dos casos, es decir, en el sentido del campo y contraria al campo. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
174 Prob. 3: Un generador síncrono 3φ, 60Hz, conexión estrella tiene 8 polos, 96 ranuras en el estator y un acortamiento de paso 9/12. La densidad de flujo en el entrehierro muestra que el tercer y quinto armónicos están presentes, y tienen amplitudes del 30% y 15% de la onda fundamental.
Calcular la relación
proporcional entre la tensión de línea y la tensión de fase. Considerar: E α Bmax KW, E L = 3 ( E12 + E 52 + .....)
1 2
Donde: KW: Factor de devanado. E f = ( E12 + E 32 + E 52 + .....)
“Estudio de las Máquinas Eléctricas. Simulación Digital”
1 2
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
175
EXAMEN FINAL 2 MAQUINAS ELÉCTRICAS I Prob. 1: Se tiene un transductor rotativo de excitación simple, cuya bobina está conformada por 450 espiras y una permeancia característica (P), según como se muestra en a figura adjunta.
P
Π 2 K − ( Ni ) Π , 0 < θ ≤ 2 Tei = K Π + ( Ni )2 , <θ ≤ Π Π 2
0,05
Π
Π
2
3Π 2
2Π
0
1a. Deducir el torque instantáneo. 1b. Calcular el torque promedio si la bobina es excitada con una corriente i(t) = 2 x 15 sen (377t + 15°) Te AV
2 KN 2Imax = − sen(2δ ) , TeAV= -612,25 N.m. ΠW
Prob. 2: Se ha formado un banco trifásico con tres transformadores monofásicos de relación de transformación 380/220 Voltios, el tipo y grupo de conexión del banco trifásico obtenido es Yd1. Se reúne el borde U de A.T. con el borde µ de “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
176
BiT, y se alimental el lado de A.T. con una fuente de tensión trifásica de valor eficaz 380 voltios. YWu = 380 Voltios. YWυ = YWW = 277, 37 V Calcular las tensiones entre los bornes W - w, W - υ, W - u. Prob. 3: Un generador síncrono trifásico, conexión estrella, 6 polos, tiene una
densidad de flujo espacial expresado por B(e) = 1,00 sen θ + 0,45 sen ∂e + 0,25 sen 5θ. El estator tiene 36 ranuras y un devanado trifásico balanceado de doble capa, el paso de bobina es 120°.
En términos de la tensión eficaz fundamental E1.
Determinar: Ef 1 = EL 3
3a. La relación entre la tensión de fase (ELN) con respecto a la de línea (Eu). Prob. 4: Un rotor gira a una velocidad angular Wr rad/seg. de manera que la
inductancia de los devanados excitadores del dispositivo electromagnético que se muestra en la figura varia en forma sinusoidal entre un máximo de 2L y un mínimo de cero.
δ
i 0 N
o
λ
Wr
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
177
4a. Determinar el torque promedio para W = Wr. TAV =
1 2 LIef sen 2δ 2
4b. Determinar el voltaje en el devanado excitador cuya resistencia de la bobina es despreciable. V(t) = -Wr L Im sen Wr t (1 + cos 2δ) - 2Wr L Im cos Wr t sen 2θ
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
178
EXAMEN PARCIAL 1 Prob. 1. Se desea diseñar con núcleo de hierro para ser utilizado con una lámpara
fluorescente, cuya característica es: 220V, 60Hz, 9W; el cual debe tener una reactancia de 565.5 ohmios.
La plancha utilizada es de material H-23, cuyo
espesor de lamina es 0.5 mm.
las dimensiones son mostradas en la figura
adjunta. Considerar: factor de apilamiento igual a 0.98, Bmax = 1.15T, la bobina tiene 300 espiras. 6a
5a lg
a
a
a
2a
1a) Calcular el numero de laminas utilizado. 1b) Si practicamos un entrehierro en la columna central (λg), de tal modo que se consiga un valor de la inductancia que sea igual al 90% del valor inicial. Calcular dicha longitud de entrehierro. 1c) considerando una característica no lineal de la inductancia, donde se cumple: i = C0λ +C1λ3. Utilizando el circuito equivalente del reactor y aplicando una tensión v ( t ) = 2 x 220sen377t . Plantear la ecuación diferencial de la corriente de magnetización. 1d) Si al circuito magnético aplicamos un tensión de amplitud Y = 220 voltios y de forma onda cuadrada, a una frecuencia de 60Hz. Despreciar la resistencia de la bobina. 1d.1) Calcular la densidad magnética máxima. Graficar la tensión aplicada y el φ máximo. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
179
1d.2) Si se desea obtener una densidad magnética máxima de 1.25 T. Calcular el voltaje máximo aplicado. Prob. 2. Un transformador 1φ, 10 KVA, (10± 2 X 2.5%) / 0.23, 60 Hz, fue sometido a ensayos de laboratorio obteniéndose los resultados siguientes:
Prueba de vacío: 230V, 1.30A, 150W; Prueba de cortocircuito: 600V, 1.50A, 815W. 2.1) Si el transformador alimenta una carga de 15 KVA, f.d.p. 0.81 capacitivo y la tensión aplicada es de valor nominal. 2.1a) Calcular la tensión en la carga (despreciar la corriente del circuito de Rp y Xm). 2.1b) Calcular el % de regulación y eficiencia. 2.1c) Si se mantiene la tensión y corriente de carga constante, y esto se varia el f.d.p. de la carga. Determinar el ángulo de la carga para obtener regulación cero. 2.2) Al transformador aplicamos una tensión primaria nominal y conectamos una carga de 11.50 KVA con f.d.p. 0.866 en atraso. Si deseamos mantener la tensión en la carga 230 Voltios. Determinar la posición del tap. 2.3) si el transformador durante todo el día presenta el ciclo de carga siguiente: %P.C. 25 50 75 100 110 f.d.p. (ind) 0.95 0.98 0.85 0.90 0.99 horas 6 6 6 3 3 2.3a) Calcular las perdidas totales. 2.3b) Calcular la eficiencia.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
180
EXAMEN PARCIAL 2 Prob. 1: Se desea diseñar un reactor con núcleo de hierro, cuya plancha a ser
utilizada es de material H23 de espesor 0.50 mm cada lámina; cuyas características técnicas son 110V, 60Hz, 10W, reactancia del reactor 7180 ohmios, 600 vueltas, factor de apilamiento 0.98, densidad máxima 1,25 Tesla. La configuración geométrica del núcleo es la siguiente: tres columnas, dos ventanas, el ancho de la columna central es igual al doble del ancho de la columna del extremo, el ancho de la ventana es igual al ancho de la columna extrema, la altura de la ventana es igual al triple del ancho de la ventana y la altura del núcleo es cinco veces el ancho de la columna del extremo. 1a) Calcular el peso del material magnético a ser utilizado. 1b) Si practicamos un entrehierro en la columna central (lg) de tal forma que se consiga un valor de la inductancia que sea igual al 80% del valor inicial. Calcular dicha longitud de entrehierro. Prob. 2: Un transformador monofásico de 100 kVA, 60Hz, (4800 +-5x1%)/230
Voltios, tiene un circuito equivalente exacto cuyos parámetros son: R1 = 0.805 Ohmios, X1 = 2.034 Ohmios, Xm = 5.487 Ohmios. R2 = 0.003 Ohmios, X2 = 0.007 Ohmios, Rp = 86.02 Ohmios El transformador es conectado a diferentes tipos de cargas, según lo indicado en las preguntas 2a), 2b) y 2c) respectivamente. 2a) Si se conecta una carga en el lado de baja tensión que, consume una corriente 1.10 ln a una tensión Nominal y el factor de potencia 0.866 en atraso. Calcular: 2a1) La tensión que es necesario aplicar al lado de alta tensión. 2a2) La eficiencia del transformador.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
181
2b) Si el transformador alimenta a una carga que consume una corriente de 1,25 ln y con un factor de potencia, tal que, se obtenga un valor del porcentaje de regulación máxima. Calcular: 2b1) La tensión en la carga si, la tensión primaria aplicada es igual a la tensión nominal. 2b2) La posición del tap si, la tensión de alimentación del primario es 4 796 voltios y, se desea mantener en la carga la tensión nominal. 2c) Si el transformador alimenta una carga que consume una corriente de 1,15 ln y con una tensión primaria de 4944 Voltios. Calcular la tensión y el factor potencia en la carga para conseguir un porcentaje de regulación nula (utilizar circuito equivalente aproximado). 2d) Si el transformador conectado a cargas variables durante todo el día presenta el ciclo de carga siguiente: %P.C.
20
40
60
100
110
f.d.p. (induct.)
0.92
0.94
0.95
0.91
0.97
Tiempo (Hr.)
6
6
5
4
3
Calcular la eficiencia del transformador.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
182
EXAMEN PARCIAL 3 Prob. 1: Se desea diseñar un reactor con núcleo de hierro, cuya plancha a ser
utilizada es de material H23 de espesor 0,50 mm cada lamina; cuyas características técnicas son 220V, 60Hz, 10W, reactancia del reactor 7180 ohmios, 600 vueltas, factor de apilamiento 0,98, densidad máxima 1,30 Tesla. La configuración geométrica del núcleo es la siguiente: tres columnas, dos ventanas, el ancho de la columna central es igual al doble del ancho de la columna del extremo, el ancho de la ventana es igual al ancho de la columna extrema, la altura de la ventana es igual al triple del ancho de la ventana y la altura del núcleo es cinco veces el ancho de la columna del extremo. 1a) Calcular el peso del material magnético a ser utilizado. 1b) Si practicamos un entrehierro en la columna central (lg) de tal forma que se consiga un valor de la inductancia que sea igual al 90% del valor inicial. Calcular dicha longitud de entrehierro. Prob. 2: Un transformador monofásico de 100 kVA, 60Hz, (4800 +-5x1%)/230
Voltios, tiene un circuito equivalente exacto cuyos parámetros son: R1 = 0.805 Ohmios, X1 = 2.034 Ohmios , Xm = 5.487 Ohmios R2 = 0.003 Ohmios, X2 = 0.007 Ohmios , Rp = 86.02 Ohmios El transformador es conectado a diferentes tipos de cargas, según lo indicado en las preguntas 2a), 2b) y 2c) respectivamente. 2a) Si se conecta una carga en el lado de baja tensión que, consume una corriente 1.10 ln a una tensión Nominal y factor de potencia 0.866 en atraso. Calcular: 2a1) El porcentaje de variación de la tensión aplicada con respeto a la tensión Nominal. 2a2) La eficiencia del transformador.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
183
2b) Si el transformador alimenta a una carga que consume una corriente de 1,25 ln y con un factor de potencia tal que se obtenga un valor del porcentaje de regulación máxima. Calcular: 2b1) La tensión en la carga si, la tensión primaria aplicada es igual a la tensión nominal. 2b2) La posición del TAP si, la tensión de alimentación del primario es 4796 voltios y, se desea mantener en la carga la tensión nominal. 2c) Si el transformador alimenta una carga que consume una corriente de 1,15 ln y con una tensión primaria de 4944 Voltios.
Calcular la tensión y el factor de
potencia en la carga para conseguir un porcentaje de regulación nula (utilizar circuito equivalente aproximado). 2d) Si el transformador conectado a cargas variables durante todo el día presenta el ciclo de carga siguiente: %P.C.
20
40
60
100
110
f.d.p. (induct.)
0.92
0.94
0.95
0.91
0.97
Tiempo (Hr.)
6
6
5
4
3
Calcular la eficiencia del transformador
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
184
SEGUNDA PRACTICA CALIFICADA Prob. 1: Un transformador monofásico de 100KVA, 60Hz, (3800 ± 5 x 2.0%) / 230
Voltios, fue sometido a los ensayos de vacío y cortocircuito a condiciones nominales, ambas lecturas se tomaron con alimentación por el lado de baja tensión, de las que se obtuvieron: Prueba de vacío: 230V. 18A, 3312W. Prueba de cortocircuito: 12V. La eficiencia máxima se obtiene para una carga del 86% de plena carga. R1 = 3.22 Ohm. R2 = 0.0118 Ω RP = 15.97 Ω Xm = 24.30 Ω X1 = 1.95 Ω X2 = 0.0071 Ω 1.a) Determinar los parámetros del circuito equivalente exacto. 1.b) Si se alimenta una carga que consume una potencia 110 KVA a un f.d.p. igual a 0.866 capacitivo a una tensión de 230 Voltios. Calcular la tensión del primario. V4 = 3913∠2.79º Voltios. 1.c) Utilizando el circuito equivalente aproximado. Si se alimenta al primario con una tensión de 3800 Voltios y se conecta al lado de baja tensión una carga que consume una potencia de 120 KVA y f.d.p. tal que el valor del porcentaje de regulación sea máxima. Calcular la tensión en la carga. V1 = 215 Voltios. “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
185
1.d) Si el transformador alimenta la misma carga de la pregunta 1c), siendo la tensión de alimentación al primario de 3700 voltios, y se desea mantener la tensión en la carga a 230 Voltios. Determinar la posición del tap. posición: 10, (-4 x 2.0%) VNP. Prob. 2: Se dispone de tres transformadores monofásicos de 1,0 KVA, 60 Hz,
110/440. Se desea alimentar una carga trifásica resistiva en 380 voltios desde una red de alimentación trifásica. 2.1) Dibujar el esquema de conexiones. Smax = 6.0 KVA 2.2) La potencia máxima que se puede suministrar a la carga, sin sobrecargar los trafos 4φS. 2.3) Calcular la tensión de la fuente de suministro. Vsum = 190 Voltios. Prob. 3: Marcar Verdadero (V) ó Falso (F) y justificar en el cuadernillo.
3.1) La regulación de tensión del transformador no depende necesariamente del tipo de carga para un mismo factor de potencia.
()
3.2) Para un transformador de 500 KVA donde las perdidas en el fierro representan el 62% de las perdidas en el cobre nominal la eficiencia máxima se obtiene cuando alimenta una carga de 394 KVA
()
3.3) Si un transformador trabaja a tensión y frecuencia nominal, luego si se alimenta con la misma tensión eficaz pero una frecuencia igual al doble de la nominal, entonces la densidad de flujo máxima y las perdidas por histéresis se duplican. “Estudio de las Máquinas Eléctricas. Simulación Digital”
() Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
186
SEGUNDO CONTROL DE MAQUINAS ELÉCTRICAS 1.
Se tiene un transformador monofásico de 10 kVa, 10/0.23 KV, 60 hz, se han realizado la prueba de cortocircuito con 550V, 1.20 Amp, y 600 Watts; y la prueba de vacío son 230V, 1.739Amp, y 200 Watts.
Determinar: 1.a) El circuito equivalente exacto referido al lado de T A.T. 1.b) El diagrama fasorial cuando se alimenta una carga de potencia nominal y f.d.p. = 0.866 (capacitivo). 2. Si el transformador del problema 1, alimenta una carga cuya potencia es el 10% más que la potencia del transformador y f.d.p. = 0.809 (capacitivo) manteniendo la tensión de entrada en 10 kV. Calcular: 2.a) La tensión en la carga (despreciar la corriente de vacío.)
224 V.
2.b) La regulación para la carga mencionada. %r = 2.678%. 3. Para el transformador utilizado, si se alimenta una carga nominal con f.d.p. 0.707 inductivo. Calcular: 3.a) La eficiencia a plena carga. np.c. = 91.97%. 3.b) La corriente de carga para obtener la eficiencia máxima. I = 30.121 Amp. 4. Poner verdadero (V) o Falso (F) las sentencias siguientes justificando su respuesta en el cuadernillo.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
187
4.a) Cuando un transformador este operando a media carga y a tensión nominal, pasa a operar a plena carga y a tensión nominal; entonces las perdidas en el fierro aumentan y las perdidas en el cobre se mantienen constantes.
( )
4.b) La tensión en vacío de un transformador con núcleo de hierro es siempre mayor que la tensión a plena carga.
( )
4.c) La corriente en vacío de un transformador de potencia es siempre menor que la de un transformador de audio.
( )
4.d) Para un transformador de audiofrecuencia, la ganancia máxima es obtenida cuando trabaja a altas frecuencias.
( )
4.e) La eficiencia de un transformador bajo las mismas condiciones de tensión y corriente es mayor con carga resistiva que con respecto a una carga inductiva. ( ) 5. La toma variable “k” del autotransformador de la fig. mostrada se encuentra en el punto medio del devanado, siendo las corrientes I1, I2, I3. Cómo varían estas corrientes cuando “k” se desplaza hacia arriba.
I1 V
I3
K` I2
“Estudio de las Máquinas Eléctricas. Simulación Digital”
RL
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
188
TERCERA PRACTICA CALIFICADA PROB. 1: Se ha formado un banco trifásico con tres transformadores monofásicos
de relación de transformación 220/110 Voltios, el tipo y grupo de conexión del banco trifásico obtenido es Yd1. Si se une los bornes U de A.T. con el borne u de B.T., y se alimenta el lado de .T. con una fuente de tensión trifásica de valor eficaz 220 Voltios. Calcular las tensiones entre los bornes W-w, W-u, W-v. PROB. 2: Se tiene dos transformadores monofásicos, cuyas características
técnicas son las siguientes: Potencia (KVA)
relación de transformación
Impedancia de Cortocircuito
64
10000/230 Voltios
Z 1 = 125∠85º
50
10000/220 Voltios
Z 2 = 140∠82º
2.a) Si ambos transformadores se conectan en paralelo para alimentar a una carga que consume una potencia de 120 KVA y f.d.p. 0.85 inductivo a una tensión de carga de 230 Voltios. Calcular: 2.a.1) La potencia que entrega cada transformador. 2.a.2) El porcentaje de variación de la tensión aplicada con respecto a la tensión nominal primaria. 2.b) Si súbitamente se desconecta la carga y manteniendo la tensión en la carga a 230 Voltios. Calcular el porcentaje de variación de la tensión aplicada.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
189 PROB. 3: Se tiene tres bancos de transformadores trifásicos según mostrados en
la figura: U
V
W
u
v
w
BANCO N°
U
u
V
v
W
w
BANCO N°
U
V
W
u
v
w
BANCO N°
3.a) Determinar los grupos de conexión de cada banco. 3.b) Realizar las conexiones físicas de los bancos trifásicos para la puesta en paralelo.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
190 PROBLEMA DE APLICACIÓN:
Se ha diseñado un reactor de las siguientes características: 220 V, 60Hz. -Numero de espiras - Numero de placas - Espesor de placa - Ancho de ventana - Ancho de núcleo - Altura de ventana - Factor de apilamiento - Densidad del fierro
: N = 40.5 : n = 80 : t = 0.5 mm : 1” : 2” : 3” : 0.98 : 7.65 gr/cc.
a) Si se aplica a la bobina de excitación un voltaje senoidal de 200 volt. (eficaz) y 60 Hz, la perdida por histéresis es de 40 W/ocl y la debida a las corrientes parásitas de 20 W/ocl. La densidad de flujo es 0.93 Wb/m2. determinar las perdidas por histéresis y corrientes parásitas, cuando se aplica una tensión de: υ(t) = 250 senWt + 71.5 sen 3Wt. Asumir n= 1.6 b) Si el reactor hubiera sido ensamblado con 51 planchas del mismo material y dimensiones, pero de 0.70 mm. de espesor (manteniendo inalterable el numero de vueltas de la bobina). Luego al ser operado a los mismos valores de tensión y frecuencia del caso a) con voltaje senoidal (200V, 60Hz). Determinar las perdidas respectivas.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
191 Prob. 1: Un circuito magnético compuesto en el que varia la sección transversal se muestra en la figura. La porción de fierro es de grado H = 23 de 0,55 mm de espesor con f.a. = 0,9 los datos son:
N = 1000 espiras. L1 = 4 L2 = 40 cm. S1 = 2S2 = 10 cm2 Lg = 2 mm Sg = 10,1 cm2 φd = 0,01 mwb Calcúlese la corriente necesaria para establecer en el entrehierro una densidad De flujo de 0,6 Tesla (βg = 0,6t). Prob. 2: El núcleo del circuito magnético que se muestra en la figura es de aleación de hierro y Si/grado H-23, formado por 90 laminas de 0,5 mm de espesor y un factor de aislamiento de 0,96 (f.a), la bobina tiene 200 espiras. Determinar: a) El flujo magnético para un entrehierro de: n = 90 laminas (0,5 mm) fa = 0,96 N = 200 Prob. 3: El material empleado en la estructura de la figura es de Acero-Si con espesor de laminas 0,5 mm; se sabe que el fa =0,92. además I1= I2 = 5A, N1 = 1000 espiras; N2 = 924 esp. Determinar el flujo en el entrehierro, sabiendo además que A = C, B = D = E. (todo en milímetros) Prob. 4: El núcleo mostrado en la figura es de material ferromagnético. Grado H23 (0,5 mm de espesor), fa = 0,9, 110 V, 60 Hz. Presenta una reactancia inductiva XL = 137 con unas perdidas en el núcleo no mayor de PT = 1w y lg = 2 mm, PFe = Kg 7750 3 (densidad del material ferromagnético). Se sabe que PFe = 1,8 m Watt An ; β2max = 0,8 . K A g g
Determinar: n = ? (# de laminas) N = ? (# de espiras de la bobina). “Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002
Universidad Nacional de Ingeniería Facultad de Ingeniería Eléctrica y Electrónica Instituto de Investigación
192
Prob. 5: Las perdidas en el núcleo de un reactor de Hierro son 548 W; de los cuales 402 W son perdidas por histéresis esto cuando el voltaje aplicado en 240 Volt. Y la frecuencia 30 Hz.
a) Determinar las perdidas en el núcleo cuando el voltaje y la frecuencia se duplican. b) Determinar la frecuencia si el voltaje es 220v, la f = 60Hz. Prob. 6: Con el núcleo mostrado en la figura se diseñara un reactor para las siguientes condiciones: Eléctricas: V=180v, f=60Hz, XL=150Ω, Bmax=1,02T Núcleo: Laminas: E-I, 0,5mm Grado: H-23 Densidad: 7750 Kg/m3 Dimensiones: 2a = 25,4 mm b = 24,74 mm fa = 0,95 Perdidas = 2,15β2max(watt/Kg).
a) b) c) d)
Calcular N (# de vueltas de la bobina) Calcular lg (entrehierro necesario) Calcular IP (I de perdida), Im (magn.), f de P que presenta la bobina. Para un incremento del 50% de la longitud del entrehierro ( g) cual será el valor de la inductancia y el f de P.
“Estudio de las Máquinas Eléctricas. Simulación Digital”
Junio 2002