BOBINADOS
BOBINADOS
BOBINADOS
AUTOR:
DIRECTOR:
CUENCA-ECUADOR
2014
1
Ingeniero FLAVIO QUIZHPI Director de Tesis
2
DECLARACIÓN DE RESPONSABILIDAD
Los conceptos desarrolados, análisis realizados y las conclusiones del presente trabajo,
son de exclusiva responsabilidad de los autores, con trabajo y obtención de los materiales y
herramientas facilitados por el taller de electrcicidad perteneciente a la Universidad Politéc-
nica Salesiana Sede Cuenca, y autorizo la utilización de la misma con fines académicos.
Cuenca, Noviembre del 2014
3
Dedicatoria
A mis queridos padres Mateo y Beatriz que con su esfuerzo, dedicación y apoyo
incondicional siempre me incentivarón a conseguir todos los objetivos que me he propuesto
en la vida, en especial la culminación de la carrera universitaria y del presente trabajo de
tesis
A mi querida hermana Katty que me apoyó en todo momento con el ánimo necesario para
continuar y completar los estudios.
4
Agradecimientos
5
Índice general
2. PRUEBAS DE LABORATORIO 27
2.1. Pruebas de corriente continua. . . . . . . . . . . . . . . . . . . . . . . . . 27
2.2. Pruebas de rotor bloqueado. . . . . . . . . . . . . . . . . . . . . . . . . . 29
2.3. Pruebas de Vacío. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.4. Circuito Equivalente. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.5. Cálculo de la eficiencia del motor de inducción. . . . . . . . . . . . . . . . 36
2.5.1. Comportamiento del motor con la carga del edificio Mario Rizzini . 37
3. ANÁLISIS DE ARMÓNICOS 42
3.1. ANÁLISIS MEDIANTE EL MÉTODO LINEAL. . . . . . . . . . . . . . . 42
6
3.1.1. LÍMITES DE LA DISTORSIÓN ARMÓNICA . . . . . . . . . . . 55
3.2. Análisis mediante un software como el Matlab. . . . . . . . . . . . . . . . 56
3.2.1. Torque vs Tiempo . . . . . . . . . . . . . . . . . . . . . . . . . . 59
3.2.2. Velocidad vs Tiempo . . . . . . . . . . . . . . . . . . . . . . . . . 59
3.2.3. Corriente vs Tiempo . . . . . . . . . . . . . . . . . . . . . . . . . 60
3.2.4. Impedancia vs Tiempo . . . . . . . . . . . . . . . . . . . . . . . . 61
3.3. Análisis técnico de los armónicos del motor de inducción de paso entero y
de paso acortado. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
3.4. Evaluación de las ventajas y desventajas del motor de inducción de paso
acortado. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
3.4.1. Ventajas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
3.4.2. Desventajas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
3.5. Análisis de resultados y aplicaciones factibles para un sistema eléctrico. . . 62
4. Conclusiones y Recomendaciones 63
4.1. Conclusiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
4.2. Recomendaciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
Bibliografía 64
7
Índice de figuras
1.1. Bobina . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.2. Bobinado de paso acortado . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.3. Bobinado trifasico de doble capa . . . . . . . . . . . . . . . . . . . . . . . 18
1.4. Motor de inducción trifásico . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.5. Datos en la placa del motor . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.6. Esquema básico simplificado del devanado trifásico de paso acortado en re-
presentación lineal (a) y circunferencial (b) para un parde polos . . . . . . 23
1.7. Bobinado final . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.1. Fasores de la Io . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.2. Equivalente Thevenin. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.3. Tensión y corriente a lo largo de 24 horas. . . . . . . . . . . . . . . . . . . 37
2.4. Corriente Máxima L1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.5. Corriente Máxima L2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.6. Corriente Máxima L3. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
2.7. Tensión Máxima L1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.8. Tensión Máxima L2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
2.9. Tensión Máxima L3. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.10. Forma de onda L1. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.11. Forma de onda L2. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
2.12. Forma de onda L3. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
8
3.10. (a)Armónicas aguas arriba-(b)Armónicas aguas abajo . . . . . . . . . . . . 52
3.11. Tensiones aguas arriba-aguas abajo L3. . . . . . . . . . . . . . . . . . . . 53
3.12. (a)Armónicas aguas arriba-(b)Armónicas aguas abajo . . . . . . . . . . . . 54
3.13. Limites de armónicos.[1] . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
3.14. Torque vs Tiempo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
3.15. Velocidad vs Tiempo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
3.16. Corriente vs Tiempo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
3.17. Impedancia vs Tiempo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
9
Índice de cuadros
10
Algoritmos
11
INTRODUCCIÓN
La presente tesis se la realiza a partir del análisis de armónicos, cuáles son los problemas
que producen, identificar los armónicos que presentan un mayor número de perturbaciones
en la red eléctrica: a partir de estos planteamientos buscar soluciones que permitan eliminar
o reducir el nivel de afectaciones que provocan como son el funcionamiento incorrecto de
equipos que fuerón diseñados y construidos para operar en condiciones normales, es decir,
con poca distorción armónica, los motores de inducción son sensibles a los armónicos y se
ven sometidos a todas las variaciones de la fuente de potencia, lo que afecta su funciona-
miento y las características propias de trabajo.
Partiendo del trabajo ya realizado por los Ingenieros Eduardo Valle y Gonzalo Contreras[2],
se procede a modificar el paso de las bobinas que se encuentran en el filtro pasivo instala-
do en la red eléctrica de la Universidad Politécnica Salesiana Sede Cuenca, para éste fin se
realiza los cálculos para el diseño de las nuevas bobinas que contendrá el motor, se analiza
cuantos grados eléctricos es lo que se va a recortar las bobinas y una vez obtenidos los resul-
tados de los cálculos, se procede a construir el motor para armónicos con el nuevo bobinado,
se verifica que se encuentre correctamente construído, se realizan las pruebas pertinentes ba-
sadas en las recomendaciones de la IEEE std 4[3], éstas son las pruebas de vacío, pruebas
de corriente contínua y prueba de rotor bloqueado, en cada una de estas pruebas se verifica
diferentes aspectos del funcionamiento del motor, una vez comprobado el mismo, se procede
a instalarlo en la red eléctrica de la Universidad.
El proceso para comprobar el comportamiento del motor consiste en realizar mediciones
constantes en el periodo de mayor carga a la red eléctrica, para este fin en primer lugar
se utiliza el medidor de calidad de energía eléctrica marca FLUKE modelo 1735 POWER
LOGGER ANALYST para determinar cuál es el período de tiempo donde se presenta mayor
carga.
Con los datos obtenidos en las mediciones se puede identificar los lapsos de tiempo con
mayor carga y realizar mediciones para comparar las gráficas obtenidas en las mediciones y
determinar si se obtuvo mejoras o perjuicios al cambiar el tipo de devanado a paso acortado.
Para éste propósito se utiliza el osciloscopio digital Marca: GwINSTEK Modelo: GDS-
1102A-U, y los datos se los obtiene utilizando el software libre FREEWAVE V3.23, las
sondas que se utilizan para medir la corriente estan atenuadas y multiplicadas por un factor
de 100X para obtener la medición en Amperios. Ésta medición se realiza colocando las
sondas de corriente atenuadas simultaneamente en los terminales de la bobina al ingreso y
a la salida de la misma para comparar el comportamiento de las corrientes al mismo tiempo
aguas arriba y aguas abajo del motor.
Para verificar como es el comportamiento de las tensiones se debe tomar mediciónes
entre líneas, para éste propósito se utiliza el osciloscopio digital Marca: GwINSTEK Modelo:
12
GDS-1102A-U, y los datos se los obtiene utilizando el software libre FREEWAVE V3.23,
las sondas que se utilizan para medir la tensión están atenuadas y multiplicadas por un factor
de 500X para obtener la medición en Voltios, como se está trabajando con sistema trifásico
se toma las mediciónes con la combinación de todas las líneas.
Mediante la utilización de simuladores matemáticos, en la presente tésis se utiliza el
Matlab se puede cálcular y gráficar la distorción arm´nica total (THD), que es un parámetro
muy utilizado para el análisis de los armónicos.
El objetivo a verificar es comprobar cuál es el comportamiento de los armónicos cuando
se modifica el paso del bobinado en la construcción de un filtro pasivo para armónicos, con
los datos obtenidos se comprueba si se mejora en la disminución de los armónicos y el
análisis final correspondiente.
13
Capítulo 1
INDUCCIÓN.
14
Figura 1.1: Bobina
15
1.1.1. NÚMERO DE RANURAS.
Es el número total de agujeros o ranuras en el cuerpo del estator donde debemos insertar
el arrollamiento de alambres.
Existe un parámetro denominado el número de grados geométricos, es la separación entre
dos ranuras consecutivas, cuya fórmula1.1 es siempre:
360 ∗ p
α= (1.1)
Nr
p= pares de polos
Nr= número de ranuras
Ésta fórmula se debe a que el número de ranuras Nr están uniformemente distribuidas en los
360 grados que tiene la circunferencia del estator.
Nr
τp = +1 (1.2)
2p
τ p = paso polar
Nr= número de ranuras
2p= número de polos
16
1.1.3.2. PASO ACORTADO.
Se logra un paso acortado1.2 cuando se toma un paso que contenga menos ranuras inter-
medias de las que requiere el paso diametral, por ejemplo un motor de Nr = 30 ranuras y 2p
= 2 polos, el paso diametral será1.2:
30
τ p = Nr
2p + 1 = 2 + 1 = 16 ranuras, es decir, como resultado tenemos 1:16
Un paso acortado en una ranura, si toma el valor de 1:15 ,que es el resultado de restar
una ranura de las 16 originales.
Un paso acortado en dos ranuras si se toma 1:14 ,que es el resultado de restar dos ranuras
de las 16 originales.
En éste caso se tomará un paso que contenga más ranuras intermedias de las que requiere
el paso diametral, por ejemplo un motor de Nr = 30 ranuras y 2p = 6 polos, el paso diametral
será1.2:
30
τ p = Nr
2p + 1 = 6 + 1 = 6 ranuras, es decir, como resultado tenemos 1:6
Un paso alargado en una ranura da al tomar 1:7 ,que es el resultado de aumentar una
ranura a las 6 originales.
También será un paso alargado en dos ranuras si se toma 1:8 ,que es el resultado de
aumentar dos ranura a las 6 originales.
La ecuación 1.3 para el cálculo de las ranuras por polo y fase expresa la cantidad de
ranuras que formarán los polos del estator.
Nr
q= (1.3)
4∗N f
17
1.1.4. BOBINADO TRIFÁSICO DE DOBLE CAPA.
Se define como capa doble1.3 cuando en cada ranura hay dos lados de bobina y están
separadas por un material aislante, en un bobinado como éste, el número total de bobinas es
igual al número de ranuras del núcleo estatórico.
La figura 1.3 demuestra como queda el diagrama del bobinado que se procede a armar
en la presente tesis, y se la puede describir como:
Número de ranuras: 36
Velocidad Nominal en Revoluciones por minuto: 1800
Series:
U - ranura 1
V - ranura 7
W- ranura 13
Se tienen 9 puentes:.
9-18, 11-19, 27-36, 15-24, 17-25, 33-6, 21-30, 23-31, 3-12
Paso concentrico:
6 bobinas 2 elementos
1-8 / 1-10 (activo 1-9)
6 bobinas 1 elemento paso 1-8
18
El motor de inducción trifásico1.4 a utilizar es el que se encuentra funcionando en el
sótano de la Universidad Politécnica Salesiana sede Cuenca.
DATOS
Número de fases 3
Voltaje 220/440V
Velocidad 3450rpm
Frecuencia 60Hz
Corriente 2.2/3.8A
Potencia 3Hp
19
Figura 1.5: Datos en la placa del motor
Para comenzar a diseñar el bobinado del motor de inducción se considera cuáles son las
dimensiones del estator1.2,luego se establece la velocidad que se puede seleccionar depen-
diendo del número de polos, para el presente caso la velocidad nominal es de 1800 rpm que
nos da el equivalente a 4 polos, la característica que debe tener el motor es que las bobinas
sólo contienen una espira, las fórmulas que se utilizan son las que se encuentran en el libro
de Giorgio Crisci[4].
Con lo mencionado anteriormente los parámetros de que disponemos son:
Número de ranuras Nr = 36
Número de polos 2p = 4
Polos p = 4/2 = 2
Número de fases N f = 3
Iniciamos calculando el número de grados eléctricos1.1 que tiene el estator:
360∗p
α= Nr
20
obtenemos p de:
2p = 4
4
p= 2
p=2
entonces se tiene:
360∗2
α= 36
Éste valor de 20° grados eléctricos es un valor de desplazamiento de una ranura a otra,
por lo que representa los grados que se puede recortar el paso de cada bobina para poder
comprobar si es que existe una mejora o no en reducir el número de armónicos motivo de la
presente tesis.
El paso polar1.2 que tiene el filtro originalmente es:
τ p = Nr
2p + 1
τ p = 36
4 +1
τp = 9+1
τ p = 10
La relación es 1:10 que es el paso entero de bobina y que es igual al paso polar.
Como el objetivo es trabajar con paso acortado y de acuerdo a los cálculos anteriores se
reduce el paso en dos ranuras:
Paso entero τ p = 10
Se realiza un paso acortado en dos ranuras, por lo tanto:
10 - 2 = 8
Entonces se toma la relación 1:8
El número de ranuras por polo y fase 1.3 que se tiene es:
Nr
q= 4∗N f
36
q= 4∗3
q = 3rp f
Como es un devanado doble capa, el número de bobinas elementales por par de polos y
fase tiene el valor de:
b=2*q
donde b es el número de bobinas elementales
21
calculando tenemos:
b=2*3=6
Ahora debemos calcular el factor de devanado1.4:
Kw = K p ∗ Kd (1.4)
donde:
Kw=factor de devanado
K p=factor de paso de un devanado1.6
Kd=factor de distribución1.5
Para calcular el factor de distribución Kd se utiliza la siguiente ecuación:
sen q ∗ ( α2 )
Kd = (1.5)
q ∗ sen( α2 )
donde:
q: número de ranuras por polo fase
α: número de grados eléctricos
sen(3∗ 20
2 )
Kd = 3∗sen( 20
2 )
sen(30°)
Kd = 3sen(10°)
Kd = 0,96
ρ
K p = sen (1.6)
2
donde:
ρ: paso polar en grados eléctricos
El paso eléctrico1.7 se lo obtiene de los datos del motor, éste es de 2 polos entonces:
180°
ρ= (1.7)
p
180°
ρ= 2
ρ = 90°
por lo tanto:
22
K p = sen ρ2
K p = sen 90°
2
K p = 0,707
Para éste caso “el factor de paso de un devanado es diferente para cada frecuencia
armónica”[5]
El factor de devanado1.4 es:
K = 0,707 ∗ 0,96
K = 0,679
Donde:
Ncr= 1 Conductor por ranura
Figura 1.6: Esquema básico simplificado del devanado trifásico de paso acortado en representación
lineal (a) y circunferencial (b) para un parde polos
Zf
Ncr = 3 ∗ (1.8)
Nr
de donde se despeja la Z f de 1.8 y se obtiene:
23
Ncr∗Nr
Zf = 3
Z f = 12 conductores
Para conocer con cuantos voltios se va a trabajar en cada fase1.9 se tiene:
Z f ∗C ∗ K p ∗ φint ∗ lg ∗ B
v= (1.9)
p
donde:
p: polos
Z f : número de conductores por fase
C: constante de la frecuencia de 60Hz C = 133
K p: factor de paso de un devanado
φint : diámetro interno del estator[m]
lg: largo geométrico del estator [m]
B: densidad de flujo magnético B = 0.8 [Wb/m2 ]
12∗133∗0,707∗0,098∗0,092∗0,8
v= 2
8,139
v= 2
24
ciones mediante instrumentos de uso general en la mayoría de las pruebas o de uso especifico
cuando se necesita mayor detalle en las pruebas o mediciones.
A continuación se cita el alcance que pretenden estas recomendaciones:
“Alcance.. Está recomendación práctica pretende establecer los principios para el di-
seño de sistemas eléctricos que incluyan cargas lineales y no lineales. Son descritas las
formas de onda de tensión y corriente que pueden existir en todo sistema, y se establecen
los principios de distorsiones de formas de onda para el diseño de sistemas. Es descrita la
interfaz entre fuentes y cargas como el punto de acoplamiento común; y la observación de
los principios de diseño para minimizar la interferencia entre los equipos eléctricos. Está
recomendación práctica está dirigida a limitaciones de estado continuo. Se pueden encon-
trar condiciones transitorias que excedan estas limitaciones. Este documento establece la
calidad de potencia que será proporcionada al punto de acoplamiento común. Este docu-
mento no cubre los efectos de interferencia en radio-frecuencia; sin embargo, incluye la
interferencia electromagnética con los sistemas de comunicación”[1].
En referencia a los armónicos la norma limita la distorsión armónica de tensión total a
5 %, e individual a 3 %, donde nos presenta límites para evitar calentamientos y pérdidas de
vida útil en los motores cuyos límites para los armónicos de secuencia negativa son muy
altos, los niveles de armónicos individuales pares y de secuencia negativa se deben limitar a
menos del 3 %, sobre todo los de menor orden.
La norma nos brinda una medición llamada distorsión armónica total (THD) que es usada
para definir el efecto de los armónicos al voltaje del sistema de potencia que es usada en
sistemas de baja tensión, media tensión, y alta tensión, se la expresa como un porcentaje de
la fundamental
q y usa la siguiente fórmula[1]:
T HD= sumadetodosloscuadradosdelasamplitudesdetodoslosvolta
cuadradodelaamplituddelvolta je f undamental
jesarmónicos
∗ 100 %
q
∑50 2
h=2 V h ∗ 100 %
T HD = V1
Los pasos a seguir para armar las bobinas son:
3. Se va dando forma dentro del estator plástico a los devanados que forman cada fase y
en los extremos de las bobinas colocamos espaguetis que sirven como aislantes y nos
permiten señalar las entradas y las salidas existentes, esto lo repetimos hasta completar
las tres fases.
25
4. Entre las bobinas se debe colocar papel aislante para evitar que se unan las bobinas,
éste proceso se lo repite hasta completar todo el bobinado.
5. Procedemos de la misma forma para completar la otra media fase y se las asegura
realizando un amarre simple en cada bobina,.
6. Se etiquetan las entradas y las salidas como dictan las normas americanas.
8. Por último se introduce el rotor en el estator plástico y se vuelve a armar el motor para
ponerlo en funcionamiento.
26
Capítulo 2
PRUEBAS DE LABORATORIO
Descripción Medida
Voltaje DC(4-3) 1V
Corriente DC(4-3) 95.37A
Voltaje DC(10-9) 1V
Corriente DC(10-9) 92A
Utilizando la Ley de Ohm se calcula la resistencia total 2.1del estator debido a que los
devanados para la prueba son conectados en serie y como éstos estan en doble capa se tiene:
VDC
RT = (2.1)
IDC
1V
R4−3 = 95,37A
27
R4−3 = 0,010485Ω
RT = VIDC
DC
1
R10−9 = 92
R10−9 = 0,010869Ω
RT
Rsdc = (2.2)
3
0,010485
Rs4 − 3dc = 3
RT
Rsdc = 3
0,010869
Rs10 − 9dc = 3
0
Rs = R4−3 k R10−9 (2.3)
0 R4−3 ∗R10−9
Rs = R4−3 +R10−9
0 0,003495∗0,003623
Rs = 0,003495+0,003623
0
Rs = 0,001778924Ω
La resistencia del estator en corriente alterna se ve afectada por el efecto Kelvin por lo
que se aumenta en un 10 %, entonces2.4:
Rs = 1, 1 ∗ R0S (2.4)
Rs = 1, 1 ∗ 0,001778924
Rs = 0,00195681Ω
28
2.2. Pruebas de rotor bloqueado.
El rotor de un motor jaula de ardilla es una bobina simétrica de la barra; por lo tanto,
la impedancia del motor es prácticamente igual para cualquier posición del rotor relativa al
estator. Para las máquinas que tienen un número integral de ranuras por polo, por fase en el
rotor y estator, ésta distancia debe ser igual a dos tercios de un polo para las máquinas trifa-
sicas. Para las máquinas que tienen ranuras para bobinas fraccionarias, la distancia angular
puede ser tanto como un polo completo. El rotor debe ser bloqueado de modo que no pueda
moverse y la diferencia de voltaje será elevará gradualmente hasta obtener una corriente de
valor especificado, en éste caso se asume que el deslizamiento del rotor es s = 1, los valores
del ensayo son presentados en la tabla2.2.
Valores Medidos
V cc = 0,4V
I cc = 180A
f pcc = 0,99
Pc1 112.4W Pc2 107W Pc3 114W
PcT = 333,4W
Qc1 56VAR Qc2 59VAR Qc3 56VAR
QcT = 171VAR
Sc1 125.6VA Sc2 122.18VA Sc3 127VA
ScT = 374,69VA
Para poder calcular la resistencia y las reactancias se utiliza las ecuaciones del reporte
experimental de Ríos [10]de donde2.5:
PcT
Rr = 2
− Rs (2.6)
3 ∗ Icc
333,4W
Rr = 3∗1802
− 0,00195681Ω
Rr = 0,00014732Ω
V cc
Zcc = (2.7)
Icc
29
0,4V
Zcc = 180A = 0,0022222Ω
Xσ e = Xσ r (2.9)
0,000714709
Xσ e = 2
Xσ e = Xσ r = 0,000357354Ω
30
Valores Medidos
Velocidad del eje n = 1124rpm
Vo = 0,5V
Io = 162A
p=4 f = 60Hz
f po = 0,98
Po1 80.7W Po2 79.2W Po3 78W
PoT = 237,9W
Qo1 24VAR Qo2 25VAR Qo3 18.9VAR
QoT = 67,9VAR
So1 82VA So2 82.2VA So3 80VA
SoT = 247,4VA
Con los datos obtenidos en las mediciones se calcula los siguientes parámetros.
Primero se calcula el valor de la velocidad nominal ns 2.10:
120 ∗ f
ns = (2.10)
p
120∗60
ns = 4
ns = 1800rpm
ns − n
s= (2.11)
ns
1800−1124
s= 1800
s = 0,3755
31
Figura 2.1: Fasores de la Io .
φ = cos−1 ( f po ) (2.12)
Ir = Io ∗ sen(φ ) (2.13)
Ir = 162 ∗ sen(0,2003)
Ir = 32,2321A
Ia = Io ∗ cos(φ ) (2.14)
Ia = 162 ∗ cos(0,2003)
Ia = 158,76A
VO
Rm = (2.15)
Ia
0,5
Rm = 158,76 = 0,0031494Ω
VO
Xm = (2.16)
Ir
32
0,5
Xm = 32,2321 = 0,015512Ω
Xer = Xm − Xσ e (2.17)
Xer = 0,0151546Ω
Basados en los cálculos anteriores se determinan los valores de los siguientes parámetros:
Impedancia del estator2.18
Zs = Rs + jXσ e (2.18)
Zs = 0,00195681 + j0,000357354
Zr = jXσ r
Zr = j0,000357354Ω
33
Impedancia de magnetización2.19
Rm ∗ jXm
Zm = (2.19)
Rm + jXm
0,0031494∗ j0,015512
Zm = 0,0031494+ j0,015512
Zm = 0,003024 + j0,0006141
Impedancia de carga2.20
El valor del deslizamientos en el vacío para el caso de la máquina construida es 0.3755,
por lo tanto:
Rr
ZLoad = (2.20)
s
0,00014732
ZLoad = 0,3755 = 0,00039233Ω
Voltaje Thévenin2.21
Vo
Vth = Zm ∗ (2.21)
Zm + Zs
0,5
Vth = 0,003024 + j0,0006141 ∗ (0,003024+ j0,0006141)+(0,00195681+ j0,000357354)
kVth k = 0,304032V
Impedancia Thévenin2.22
Zs ∗ Zm
Zth = + Zr (2.22)
Zs + Zm
(0,00195681+ j0,000357354)∗(0,003024+ j0,0006141)
Zth = (0,00195681+ j0,000357354)+(0,003024+ j0,0006141) + ( j0,000357354)
Rth = 0,001188Ω
Xth = j0,0005838Ω
Vth
Ir = (2.23)
Zth + Zload
34
(0,304023+0,002350 j)
Ir = (0,001188+0,0005838 j)+(0,00039233)
Ir = 169,76164 − 61,225724 j
kIr k = 180,464965A
Potencia en el eje2.24
3 ∗ kVth k2 ∗ Rr ∗ 1−s
s
Pe je = Rr 2
2
(2.24)
Rth + s + Xth
Torque eléctrico2.25
3 ∗ Rr ∗ kIr k2
Te = (2.25)
2∗π ∗ f ∗s
3∗0,00014732∗(180,464965)2
Te = 2∗3,14∗60∗0,3755
Te = 0,101729
Voltaje de magnetización2.26
Vm = 0,0884818 + 0,0366443 j
kVm k = 0,0957696V
Corriente de magnetización2.27
Vm
Im = (2.27)
Zm
desglosando un poco más la fórmula2.27 tenemos:
(Zload + jXσ r)
Im = Zm ∗ Ir
(0,00039233+ j0,000357354)
Im = 0,003024+ j0,0006141 ∗ (169,76164 − 61,225724 j)
35
Im = 30,464362 + 5,931265 jA
kIm k = 31,036386A
Ie = Im + Ir (2.28)
Ie = 200,2260 − 55,294458 j
kIe k = 207,7207A
Potencia Eléctrica2.29
PotenciaE je
E f iciencia = (2.30)
PotenciaEl écrica
Con los datos calculados se encuentra el valor de la eficiencia utilizando la ecuación2.30:
PotenciaE je
E f iciencia = η = PotenciaEl écrica
0,042
η= 0,396
η = 0,106 w 10,6 %
36
2.5.1. Comportamiento del motor con la carga del edificio Mario Riz-
zini .
Con el fin de determinar la cantidad de armónicos que se encuentran en la carga del edifi-
cio Mario Rizzini de la Universidad Politécnica Salesiana sede Cuenca, se utiliza el medidor
de calidad de energía Marca: FLUKE Modelo: 1735 POWER LOGGER ANALYST para
determinar cuál es el período de tiempo donde se presenta la mayor carga. Según las medi-
ciones que se realizo en el período de una semana, se aprecia que la mayor carga se presenta
entre las 17h30 a 18h00, de lunes a viernes, con estos datos se determina los parámetros de
índice de calidad en el sistema2.3.
37
Para la línea L2 se tiene2.5:
38
Figura 2.7: Tensión Máxima L1.
En la línea L22.8:
En la línea L32.9:
39
Figura 2.9: Tensión Máxima L3.
Las formas de onda que se captan en el medidor de calidad de energía son las siguientes:
En la línea L1 se presenta de la siguiente manera2.10:
40
Figura 2.11: Forma de onda L2.
41
Capítulo 3
ANÁLISIS DE ARMÓNICOS
Los armónicos son corrientes, tensiones o las dos al mismo tiempo presentes en un sis-
tema eléctrico, generalmente múltiplo de la frecuencia fundamental. Actualmente con el in-
cremento de cargas no lineales generadas por la electrónica de potencia han ido en aumento
los inconvenientes en las instalaciones eléctricas por la presencia de armónicos ya sean éstos
de corriente ó tensión en el sistema eléctrico. Los problemas que provocan son: el sobreca-
lentamiento de cables, transformadores y motores, altas corrientes en el neutro y fenómenos
de resonancia entre los elementos del circuito, los armónicos provocan un funcionamiento
incorrecto de varios equipos que fueron diseñados y construidos para operar en condiciones
normales es decir con poca distorsión armónica, los motores de inducción son sensibles a los
armónicos y se ven sometidos a todas las variaciones de la fuente de potencia, lo que afecta
su funcionamiento y las características propias de trabajo.
Los problemas que causan los armónicos se los puede evaluar usando la superposición.
Se compara las gráficas obtenidas en las mediciones y se determina si se obtuvo mejoras
o perjuicios al cambiar el tipo de devanado a paso acortado. Para éste propósito se utiliza el
osciloscopio digital Marca: GwINSTEK Modelo: GDS-1102A-U, y los datos se los obtiene
utilizando el software libre FREEWAVE V3.23, las sondas que se utilizan para medir la
corriente están atenuadas a 10mV/A y multiplicadas por un factor de 100X para obtener la
medición en Amperios. Esta medición se realiza colocando las sondas de coriente atenuadas
simultaneamente en los terminales de la bobina al ingreso y a la salida de la misma para
comparar el comportamiento de las corrientes al mismo tiempo aguas arriba y aguas abajo
del motor.
Los resultados son los siguientes:
En el caso de la medición aguas arriba y aguas abajo de la línea L1 se tiene3.1:
42
Figura 3.1: Corrientes aguas arriba-aguas abajo L1.
43
(a)
(b)
Figura 3.2: (a)Armónicas aguas arriba-(b)Armónicas aguas abajo
El valor del THD que se obtiene del matlab para la medición aguas arriba es:
T HD = 19,3140
El valor del THD que se obtiene del matlab para la medición aguas abajo es:
T HD = 19,2572
44
reducción de los armónicos considerable ya que los valores son casi parecidos refiriéndonos
al tercer armónico.
Las gráficas obtenidas para la línea L2 son3.3:
45
(a)
(b)
Figura 3.4: (a)Armónicas aguas arriba-(b)Armónicas aguas abajo
El valor del THD que se obtiene del matlab para la medición aguas arriba es:
T HD = 14,4083
El valor del THD que se obtiene del matlab para la medición aguas abajo es:
T HD = 14,1656
46
reducción de los armónicos considerable ya que los valores son casi parecidos refiriéndonos
al tercer armónico.
En la línea L3 se tiene3.5:
47
(a)
(b)
Figura 3.6: (a)Armónicas aguas arriba-(b)Armónicas aguas abajo
El valor del THD que se obtiene del matlab para la medición aguas arriba es:
T HD = 21,5665
El valor del THD que se obtiene del matlab para la medición aguas abajo es:
T HD = 20,5084
48
reducción de los armónicos considerable ya que los valores son casi parecidos refiriéndonos
al tercer armónico.
Para verificar como es el comportamiento de las tensiones se debe tomar mediciónes
entre lineas, para éste propósito se utiliza el osciloscopio digital Marca: GwINSTEK Modelo:
GDS-1102A-U, y los datos se los obtiene utilizando el software libre FREEWAVE V3.23,
las sondas que se utilizan para medir la tensión están atenuadas a 2mV y multiplicadas por
un factor de 500X para obtener la medición en Voltios. como se está trabajando con sistema
trifásico se toma las mediciónes con la combinación de todas las lineas.
En el caso de las tensiones se tiene los siguientes resultados:
Entre la línea L1-L2:
Como se puede ver en la figura3.7 no existe ninguna diferencia entre voltajes entre la
medición de aguas arriba y aguas abajo ya que los valores que registra el osciloscopio son
identicas entre los dos canales de medición.
La figura 3.8 muestra el valor de los armónicos para la tensión obtenida de las mediciones
en el osciloscopio y cálculadas en el matlab para la L1-L2:
49
(a)
(b)
Figura 3.8: (a)Armónicas aguas arriba-(b)Armónicas aguas abajo
El valor del THD que se obtiene del matlab para la medición aguas arriba es:
T HD = 5,0935
El valor del THD que se obtiene del matlab para la medición aguas abajo es:
T HD = 4,9918
50
reducción de los armónicos considerable ya que los valores son similares refiriéndonos al
quinto armónico.
Entre la línea L1-L3:
Como se puede ver en la figura 3.9 no existe ninguna diferencia entre voltajes entre la
medición de aguas arriba y aguas abajo ya que los valores que registra el osciloscopio son
muy parecidos entre los dos canales de medición.
La figura 3.10 muestra el valor de los armónicos para la tensión obtenida de las medicio-
nes en el osciloscopio y cálculadas en el matlab para la L1-L3:
51
(a)
(b)
Figura 3.10: (a)Armónicas aguas arriba-(b)Armónicas aguas abajo
El valor del THD que se obtiene del matlab para la medición aguas arriba es:
T HD = 4,9030
El valor del THD que se obtiene del matlab para la medición aguas abajo es:
T HD = 4,7408
52
reducción de los armónicos considerable ya que los valores son similares refiriéndonos al
quinto armónico.
Entre la línea L2-L3:
Como se puede ver en la figura 3.11 no existe ninguna diferencia entre voltajes entre la
medición de aguas arriba y aguas abajo ya que los valores que registra el osciloscopio son
identicas entre los dos canales de medición.
La figura 3.12 muestra el valor de los armónicos para la tensión obtenida de las medicio-
nes en el osciloscopio y cálculadas en el matlab para la L2-L3:
53
(a)
(b)
Figura 3.12: (a)Armónicas aguas arriba-(b)Armónicas aguas abajo
El valor del THD que se obtiene del matlab para la medición aguas arriba es:
T HD = 4,6593
El valor del THD que se obtiene del matlab para la medición aguas abajo es:
T HD = 4,5478
54
3.1.1. LÍMITES DE LA DISTORSIÓN ARMÓNICA
Para realizar un análisis de la distorsión armónica la norma IEEE Std 519 de 1992[1]limita
la distorsión armónica de tensión total a 5 %, e individual a 3 %, en ésta se presenta unos lí-
mites para evitar calentamientos y pérdidas de vida útil en los motores, los límites para los
armónicos de secuencia negativa son muy altos; los niveles de armónicos individuales pares
y de secuencia negativa se deben limitar a menos del 3 %, sobre todo los de menor orden, en
la siguiente tabla se detalla entre que rangos deben limitarse los armónicos:
55
3.2. Análisis mediante un software como el Matlab.
Para tener una mejor idea del comportamiento del motor de inducción se utilizan diferen-
tes softwares que ayudan a simular cuál sería el comportamiento del filtro bajo los parámetros
con los cuales se realizó los cálculos, en éste caso el software que se utiliza es el Matlab, para
las líneas de programación de las fórmulas y los programas la referencia es el libro del Dr.
Aller [11] y la tésis de los Ingenieros John Valle y Gonzalo Contreras [2]los cuales permiten
simular el comportamiento de los armónicos del motor de inducción.
También se realiza la simulación del modelo que se utilizó para obtener los valores de
los parámetros del filtro.
56
Algoritmo 3.2 Programación del cálculo de parámetros del motor de inducción. [2]
% Datos
n = 1 4 3 9 ; f = 6 0 ; p = 4 ; w=2 * p i * f ; n s =120 * f / p ; s = ( ns−n ) / n s ;
% C a l c u l o de Re
Vdc = 1 ; I d c = 7 5 ; RT=Vdc / I d c ; Re=RT / 3
% C a l c u l o a r o t o r bloqueado ( s =1)
Pcc = 9 0 0 ; Qcc = 5 0 3 ; Scc = 9 2 4 ; I c c = 1 6 0 ; c c = 3 . 2 ; f p c c = Pcc / Scc ;
% Pcc =3 * I c c ^2 ( Re+Rr )
Rr = ( Pcc / ( 3 * I c c ^2)) − Re ; Zcc=Vcc / I c c ;
% Zcc ^2 = ( Re+Rr ) ^ 2 + ( Xde+Xdr ) ^ 2 ;
% XT=Xde+Xdr ;
XT= s q r t ( Zcc ^2 −(Re+Rr ) ^ 2 ) ; Xde=XT / 2 ; Xdr=Xde ;
% C a l c u l o a v a c i o ( s =0)
Vo = 5 . 5 ; I o = 7 0 ; Po = 2 7 5 . 8 ; Qo = 1 2 9 6 . 2 ; So = 1 2 6 9 . 6 ;
f p o =Po / So ; a= a c o s ( f p o ) ; Ppce =3 * Re * I o ^ 2 ;
Rro = ( Po / ( 3 * I o ^2)) − Re ; Zo=Vo / I o ;
XTo= s q r t ( Zo ^2 −(Re+Rro ) ^ 2 ) ;
Xdeo=XTo / 2 ; Xdro=Xdeo ;
% Po= Ppce +Pmag+ P r o c e
P r o c e = 0 . 0 1 * Po ; Pmag=Po−Ppce−P r o c e ; Iom= I o * s i n ( a ) ; I o r = I o * c o s ( a ) ;
Xm=Vo / Iom ; Rm=Vo / I o r ; Lme=Xm/ w ; Lde=Xde / w ; Ldr =Lde ; L e r =Lme−Lde ;
Ze=Re+ i * Xde ; Zr =Rr+ i * Xdr ; Zm=Rm* i *Xm/ ( Rm+ i *Xm) ;
% Voltaje Equivalente e Impedancia e q u i v a l e n t e
Ve = 5 . 5 ; Zr = i * Xdr ; Zc=Rr / s ; Vth=Zm* Ve / ( Zm+Ze ) ; Vthm= a b s ( Vth ) ;
Z t h = ( Ze *Zm / ( Ze+Zm) ) + Zr ; Rth = r e a l ( Z t h ) ; Xth=imag ( Z t h ) ;
Z t h _ c = Z t h +Zc ; I r =Vth / Z t h _ c ; I r a b s = a b s ( I r ) ;
P e j e = ( 3 * ( Vthm ^ 2 ) * Rr * ((1 − s ) / s ) ) / ( ( ( Rth + ( Rr / s ) ) ^ 2 + ( Xth ^ 2 ) ) )
Te = ( 3 * Rr * ( I r ^ 2 ) ) / ( w* s ) ; Te= a b s ( Te )
Vm= ( Zc+ i * Xdr ) * I r ; Im=Vm/Zm; I e =Im+ I r ; I e a b s = a b s ( I e ) ;
Pe =3 * Ve * I e a b s * 0 . 9 5 ; Pehp=Pe / 7 4 6 ; n= P e j e / Pe ;
% V a l o r e s en P o r Unidad
vb = 5 . 5 ; i b = 7 0 ; s b = 1 1 5 1 . 2 ; zb =vb / ( s q r t ( 3 ) * i b ) ; Lb= zb / w ;
Re_pu=Re / zb ; % R e s i s t e n c i a d e l E s t a t o r en P o r Unidad
% C a l c u l o a r o t o r bloqueado ( s =1)
Scc_pu = Scc / s b ; Pcc_pu = Pcc / s b ; Qcc_pu=Qcc / s b ; I c c _ p u = I c c / i b ;
Vcc_pu=Vcc / vb ; Rr_pu = a b s ( ( Pcc_pu / ( 3 * I c c _ p u ^2)) − Re_pu )
Zcc_pu =Zcc / zb ; XT_pu= s q r t ( Zcc_pu ^2 −(Re_pu+ Rr_pu ) ^ 2 ) ;
% Calculo a vacio
Xde_pu=XT_pu / 2 ; Xdr_pu=Xde_pu ;
Vo_pu=Vo / vb ; I o _ p u = I o / i b ; Po_pu=Po / s b ; Qo_pu=Qo / s b ;
So_pu=So / s b ; Ppce_pu = Ppce / s b ; P r o c e _ p u = 0 . 0 1 * Po_pu ;
Pmag=Po_pu−Ppce_pu−P r o c e _ p u ; Iom_pu= I o _ p u * s i n ( a ) ;
I o r _ p u = I o _ p u * c o s ( a ) ; Xm_pu=Vo_pu / Iom_pu ; Rm_pu=Vo_pu / I o r _ p u ;
L e r _ p u = L e r / Lb ; Lde_pu=Lde / Lb ; Ldr_pu =Lde_pu ;
Los resultados obtenidos son en por unidad los cuales ayudan en el cálculo ya que se
trabaja con cantidades pequeñas y los valores son3.1:
57
Con los resultados obtenidos en la tabla 3.1 se realiza la simulación de las características
que va ha tener el motor de inducción.
dc =−[ s i n ( t h ( k ) ) , s i n ( t h ( k ) + 2 * p i / 3 ) , s i n ( t h ( k ) + 4 * p i / 3 ) ; . . .
s i n ( t h ( k )+4* p i / 3 ) , s i n ( t h ( k ) ) , s i n ( t h ( k )+2* p i / 3 ) ; . . .
s i n ( t h ( k )+2* p i / 3 ) , s i n ( t h ( k )+4* p i / 3 ) , s i n ( t h ( k ) ) ] ;
Te ( k ) = 0 . 5 * L e r * iT ’ * [ z e r o s ( 3 ) , dc ; dc ’ , z e r o s ( 3 ) ] * i T ;
v a e ( k ) = 1 . 4 1 4 2 * Ve * c o s ( T ( k ) ) ;
vbe ( k ) = 1 . 4 1 4 2 * Ve * c o s ( T ( k ) − 2 . 0 9 4 3 9 ) ;
v c e ( k ) = 1 . 4 1 4 2 * Ve * c o s ( T ( k ) − 4 . 1 8 8 7 9 ) ;
VESP ( k ) = s q r t ( 2 / 3 ) * ( v a e ( k ) + vbe ( k ) * exp ( j * 2 * p i / 3 ) + v c e ( k ) * exp ( j * 4 * p i / 3 ) ) ;
IESP ( k ) = s q r t ( 2 / 3 ) * ( i a e ( k ) + i b e ( k ) * exp ( j * 2 * p i / 3 ) + i c e ( k ) * exp ( j * 4 * p i / 3 ) ) ;
Z ( k ) = a b s ( VESP ( k ) / IESP ( k ) ) ;
end
figure (1)
p l o t ( T , Te ) t i t l e ( ’ P a r v s Tiempo ’ )
grid
figure (2)
p l o t ( T , w) t i t l e ( ’ V e l o c i d a d v s Tiempo ( k = 1 ) ’ )
grid
figure (3)
p l o t ( T ( 1 0 : n ) , Z ( 1 0 : n ) ) t i t l e ( ’ IMPEDANCIA DE ENTRADA’ )
grid
figure (4)
p l o t ( T ( 1 0 : n ) , IESP ( 1 0 : n ) ) t i t l e ( ’ CORRIENTE DE ENTRADA’ )
grid
58
3.2.1. Torque vs Tiempo
59
Se muestra en la gráfica3.15 la velocidad que gira el rotor, notamos que acelera hasta que
llega a un punto de equilibrio dependiendo de las especificaciones del motor, la velocidad
esta en revoluciones por minuto.
La corriente se comporta al empezar con valores altos, pero al llegar a una velocidad de
sincronismo tiende a estabilizarse.
60
3.2.4. Impedancia vs Tiempo
3.4.1. Ventajas
Con el motor de inducción de paso acortado no se redujo en mayor medida los efectos
de los armónicos en la red eléctrica de la Universidad.
En las pruebas de laboratorio se vio que el arranque del motor se logro a los 40 A, lo cuál
demuestra que se necesita menos amperaje para logar poner en funcionamiento a éste motor.
61
En la medida de la corriente y la tensión se logro valores menores a los de su predecesor.
El motor no produce vibraciones que produzcan un mal funcionamiento del mismo.
Como se vio en los cálculos y las mediciones se mejoro el torque del motor.
3.4.2. Desventajas
El recalentamiento del motor es un poco elevado cuando se lleva al nivel máximo que
soporta el motor.
El diseño y el armado del devanado de doble capa se lo deber realizar con el mayor
cuidado posible ya que de ello depende el buen funcionamiento del motor.
62
Capítulo 4
Conclusiones y Recomendaciones
4.1. Conclusiones
De las pruebas realizadas se observa que para arrancar el motor se necesitó de 20A
para comenzar a girar, luego el motor alcanza su velocidad nominal. En la corriente
puede alcanzar los 200A que es el valor máximo de corriente que soporta el nuevo
devanado.
La reducción de paso ha permitido tener una mejor forma de onda de la fuerza elec-
tromotriz f.e.m. pero al final no se ve una mayor reducción de los armónicos que era
el principal objetivo de la presente tesis.
Este tipo de filtro es totalmente amable con el medio ambiente ya que no produce
ningún tipo de ruido que perjudique a las personas que se encuentren alrededor del
filtro.
Las bobinas con las que se construyo el motor son más cortas en longitud en relación
a las bobinas de paso entero.
Cuando se calcula el valor de los grados eléctricos se puede tener una perspectiva de
con cuantos grados eléctricos se puede acortar el paso de bobina.
El bobinado trifásico doble capa permite manejar varios niveles de corriente depen-
diendo de la conexión que se realize entre las bobinas.
63
Al establecer el valor del THD se pudo verificar que en las condiciones en que se
desarrolla la tesis, no se ve una disminución de los armónicos aguas arriba.
4.2. Recomendaciones
Como el estator es de nylon hay que tener en cuenta que el motor no se debe recalentar
ya que se podría dañar permanentemente el bobinado del motor y el estator.
Se puede seguir investigando cuáles son las ventajas o desventajas eléctricas y porque
no mecánicas de éste prototipo.
Se puede tratar de reducir el tamaño del filtro para que se realicen aplicaciones al nivel
residencial.
Entre las múltiples aplicaciones que tienen los filtros para armónicos se puede utilizar
éstos en los motores monofásicos y bifásicos para comprobar su funcionamiento.
64
Bibliografía
[6] “IEEE 112 Standard Test Procedure for Polyphase Induction Motors and Generators,”
2004.
[7] “IEEE 115 Standard Test Procedures for Synchronous Machines,” 1995.
[8] “Norma IEEE 95 Ensayos de Alta Tensión Contínua (Máquinas nuevas y usadas ) ,”
2002.
[11] J.M.Aller, “Máquinas Eléctricas Rotativas: Introducción a la Teoría General,” vol. Pri-
mera Edición.
[13] F. Q. F. Breasciani, “Motor de Inducción de una sola espira,” vol. Edición 9, 2013.
65
[14] G. E. P. Solano and M. M. Sánchez, “Análisis de los parámetros de comportamien-
to de un sistema tetrafásico mediante el modelamiento y construcción del sistema de
generación electrónico y su aplicación a un motor de inducción tetrafásico,” p. 134,
2011.
66