VIVES,
J. C., ALVARADO, C. A., OLORTEGUI, J. A. Modelado, simulación e
implementación de módulo de levitación magnética para
investigación rotodinámica en el pregrado. Julio - diciembre
2020. Vol. 7 / Nº 2, pp. 128-141 - ISSN: 2313-1926
MODELADO, SIMULACIÓN E IMPLEMENTACIÓN DE MÓDULO DE LEVITACIÓN MAGNÉTICA PARA INVESTIGACIÓN ROTODINÁMICA EN EL PREGRADO
MODELING, SIMULATION AND IMPLEMENTATION OF A MAGNETIC LEVITATION MODULE FOR ROTODYNAMIC RESEARCH IN THE UNDERGRADUATE
Juan Carlos Vives Garnique1a
Carlos Alexis Alvarado Silva1b
Jorge Antonio Olortegui Yume1c
Fecha de recepción : 15 de Noviembre del 2020
Fecha de aprobación : 02 de Diciembre del 2020
DOI : https://doi.org/
Resumen
En este trabajo se desarrolló una interfaz gráfica de usuario (GUI) para simulación de un banco de pruebas de levitación magnética de utilidad en el diseño e instrucción virtual en el pregrado. Esta herramienta virtual basa su lógica en el diseño paramétrico y modelos matemáticos de electromagnetismo y resistencia de materiales. Esta GUI permite al usuario tomar decisiones de diseño a medida que la información de los valores de las variables relevantes se ingresa mostrando su influencia en curvas paramétricas comparativas. Como herramienta de instrucción la GUI puede utilizarse en los cursos de Electromagnetismo, Maquinas Eléctricas y Diseño Electromecánico. Mediante el uso de la GUI se logró desarrollar el diseño preliminar del propio módulo de levitación magnética el cual se encuentra en su estadio inicial de fabricación.
Este módulo permite manipular las variables de entrada y variables geométricas para calcular los datos de la corriente en polo, carga ejercida por el polo, flujo de corriente; dimensionamiento de los elementos del banco de ensayos, además de calcular las fuerzas electromagnéticas necesarias para levitar el eje de la volante. Finalmente, el dimensionamiento del eje y rodamiento magnético activo se seleccionan en base a los valores disponibles proporcionados por el fabricante.
Palabras claves: Rodamiento magnético activo, fenómeno de levitación, diseño paramétrico.
Abstract
In this research, a programming of a graphical user interface (GUI) was developed in the MATLAB environment to be used as a virtual teaching technological tool to design a rotor test bench with magnetic bearings. The programming of this virtual tool was based on the parametric design grounded on mathematical models of electromagnetism and resistance of materials, so that step by step, the user recognizes the influence of the information as it is added and determines through comparative curves the most appropriate data, improving learning in the subject through the analysis of results. This module allows to manipulate the input variables and geometric variables to calculate the data of the current in pole, load exerted by the pole, current flow; sizing of the elements of the test bench, in addition to calculating the electromagnetic forces necessary to levitate the axis of the steering wheel. Finally, the sizing of the shaft and active magnetic bearing are selected based on the available values provided by the manufacturer.
Keywords: Magnetic bearings, levitation phenomenon, parametric design, virtual module.
Los laboratorios experimentales forman parte integral de cualquier curso de ingeniería para el aprendizaje significativo de los estudiantes en todos los niveles; los educadores expertos en estas áreas han determinado que existen grandes beneficios para aprender ejecutando actividades en estos ambientes. (Ibrahim, 2010).
En la actualidad, el e-learning y los laboratorios virtuales han ganado una popularidad sustancial en la educación científica. En medio de los cierres de COVID-19, la enseñanza presencial y los cursos de laboratorio regulares han sido suspendidos en varios países del mundo. En este escenario, las clases virtuales y los recursos en línea podrían servir de manera más efectiva como una posible forma alternativa de aprender ciencia desde casa. (Ray, S., Srivastava, S., 2020)
MATLAB presenta una extensión de programación que puede ser usada como una herramienta de interfaz gráfica de usuario, o GUI, que permite al profesor de ingeniería producir softwares pedagógicos combinando la capacidad informática y gráfica de MATLAB mediante interfaces que brinden un entorno amigable de aprendizaje para los estudiantes. (R. DeMoyer , E. E. Mitchell, 1999)
Por otro lado, se sabe que las máquinas rotativas cubren una amplia gama de instalaciones críticas y constituyen la columna vertebral de numerosas industrias, desde turbinas a gas utilizadas en la producción de electricidad hasta turbo-maquinarias utilizadas para generar energía en la industria aeroespacial. Es vital que estas máquinas funcionen de forma segura a lo largo del tiempo y en diferentes condiciones operativas, para garantizar productividad y prevenir cualquier falla catastrófica, que conduciría a reparaciones extremadamente costosas y también poner en peligro la vida del personal de operaciones. (F. Al-Badour, 2011).
Uno de los mayores desafíos de las máquinas rotativas consisten en la eliminación de la fricción, lo que genera aumento de calor y vibración en todos los componentes, haciendo que las máquinas no alcancen toda su capacidad de trabajar con velocidades de rotacion y cargas elevadas, un segundo desafío, está relacionado con la lubricacion, puesto que desecharlo puede afectar directamente al medio ambiente. Todas estas desventajas pueden reducirse con el uso de los cojinetes magnéticos, los que al utilizar las fuerzas electromagnéticas para su levitación eliminan el rozamiento entre componentes haciendo que el uso de lubricantes sea innecesario. (Sayed Mahmoud, 2012)
Las pesquisas en sistemas de levitación magnética en nuestro país han ido tomando importancia, ya que en la actualidad a nivel internacional se vienen teniendo resultados favorables sobre este desarrollo tecnológico, razón por la cual este tipo de investigaciones contribuirá con el conocimiento y sustento para continuar investigando su aplicabilidad por parte de los grupos de investigación, empresas y sector educación en el progreso del país y su tecnología. Sin embargo, la conceptualización de estas teorías puede tornarse retador durante el aprendizaje de los estudiantes (Soto & Villegas, 2006)
Por lo tanto, dada la coyuntura social y de salud actual, se presenta una propuesta educativa virtual a través de una GUI (Graphical User Interface) para el estudio de caso en máquinas rotatorias que pueda ser considerada como metodología de enseñanza. Esta GUI está implementada para la enseñanza de sistemas rotodinámicos que se puede utilizar como complemento educativo en las disciplinas de Dinámica, Electromagnetismo y diseño de máquinas dando lugar a la posibilidad de estudiar casos reales propuestos evaluados por el docente.
La investigación es de tipo pre-experimental, ya que se pretende proponer el diseño paramétrico que debe efectuarse concibiendo una apropiada metodología de diseño. Para ello se plantean cuatro grupos de variables: de análisis de diseño, de solución, variables intervinientes y de conceptualización del problema siguiendo la literatura admitida para diseño mecánico; con ello se logró concebir el diseño conceptual. (Figura 1)
Figura 1. Metodología a seguir, para el desarrollo de la GUI.
Para establecer un punto de partida al diseño de la GUI y considerando que el uso de esta herramienta es dual, diseño e instrucción de estudiantes de pregrado, se realizaron dos actividades iniciales asociadas. Primero, dado que el diseño de sistemas de levitación magnética aun es de interés comercial no existe normativa definida para su fabricación, por ello, se obtuvieron recomendaciones de investigadores el área (Sikanen, 2014; Pastor & Rubio, 2015; Gaurav et al, 2016), lo que permitio identificar los modelos matematicos, componentes, características y técnicas relacionadas. En segundo lugar, se realizaron entrevistas a los docentes de la Escuela de Ingeniería Mecánica Eléctrica en conexión con los cursos de Electromagnetismo, Maquinas Eléctricas y Diseño Electromecánico. Estas entrevistas estuvieron orientadas a obtener de información de campo respecto de las necesidades del currículo a ser cubiertas por una herramienta virtual como la GUI (Tabla 1). Considerando los datos de la Tabla No1 y las recomendaciones de diseño se elaboró la tabla de especificaciones de ingeniería (Tabla 2)
Tabla 1.
Cuadro resumen de entrevistas con las Necesidades de la Institución.
|
Preguntas |
Consolidado Respuestas |
ASPECTO DE INFLUENCIA ACADÉMICA E INSTITUCIONAL |
|
Fractura de piezas del rodamiento, corrosión, desgaste, fatiga en los elementos, deformación rotacional |
|
Montaje incorrecto, insuficiente lubricación, excesiva vibración, sobrecarga, desbalance. |
|
|
Fractura, fatiga, velocidad crítica, desbalance de piezas. |
|
|
Mejorará la aplicación de la teoría de vibraciones y generará mejoras con nuevos estudios de investigación. |
|
ASPECTO DE DISEÑO DE LA MÁQUINA |
|
Ancho 1 m, largo 1.5m, altura 0.45 m. |
|
Costo, facilidad de mantenimiento, sistema de transmisión, control, resistencia. |
|
|
Revoluciones, nivel de vibración, fuerza electromagnética. |
|
ASPECTO ECONÓMICO |
|
Entre S/. 20,000.00 a S/. 50,000.00 nuevos soles. |
Tabla 2.
Especificaciones de ingeniería.
Subfunción |
Características Ingenieriles |
Unidad |
Limites |
|
|||
Dimensiones adecuadas para un laboratorio |
Largo, Ancho, Altura |
|
|
||||
Vencer la Inercia de partes Móviles |
Potencia de Arranque |
Hp |
|
||||
Controlar la estabilidad de la levitación
|
Número de polos |
Adimensional |
3, 4 ,6 ,8 |
||||
Mantener Espaciamiento de seguridad |
Holgura entre el rodamiento y eje. |
mm |
|
||||
Rotar a velocidad coherente con equipos industriales |
Velocidad máxima del eje. |
rev/min |
|
||||
Sostener la masa giratoria |
Máxima capacidad de rodamiento magnético en el rotor para mantenerlo fijo en su posición. |
N
|
130 |
Fuente: Elaboración Propia
Teniendo en cuenta que se tiene como objetivo implementar físicamente un módulo de levitación magnética que incluya un banco de ensayo de rotodinámica tradicional, i.e., rodamientos de bolas regulares, y pueda fácilmente adaptarse a un banco de ensayos con rodamientos magnéticos (levitación magnética), la GUI debe partir de un dimensionamiento preliminar como punto de partida. Tomando como base las Tablas 1 y 2 se dimensiono el banco el ensayo según la Figura 2.
Figura 2. Dimensiones generales de la máquina
Fuente: Elaboración Propia
El modelo matemático para el desarrollo de la GUI se basó en el trabajo de Matsuda & Kanemitsu (2007). De este ultimo trabajo se adapto la Figura No 3 la que se asocia a la ecuacion principal para el diseño de un rodamiento magnético activo viene dada por:
(1)
Donde:
Figura 3. Esquema del rodamiento magnético radial activo
Fuente: (Matsuda & Kanemitsu, 2007)
Tomando en cuenta que la cantidad de vueltas revertidas para los polos, entonces, ; para las dimensiones de los polos tendremos: ; la holgura para todos los casos será: Obteniendo finalmente la expresión para queda expresada como
(2)
En la Figura 4 se puede apreciar la GUI terminada mostrando los parámetros geométricos (restricciones o parámetros de definición de problema), las variables de entrada, las variables de diseño (Variables de solución de diseño) y los parámetros de solución.
Figura 4. Interfaz gráfica de usuario del diseño paramétrico
Fuente: Elaboración Propia
La figura 5 muestra las gráficas comparativas de las variables de solución como fuerza electromotriz, número de vueltas de la bobina y corriente en el polo para seleccionar la mejor opción que soportarán los rodamientos magnéticos a las solicitaciones de cargas externas.
Figura 5. Diagrama de fuerzas electromagnéticas generadas según el número de bobinas.
Fuente: elaboración propia
El diagrama de la figura 6 muestra el análisis convencional mediante el diagrama de cuerpo libre de las cargas de la volante, eje y peso del mango de laminación para un caso específico con el que se validarán los datos de la GUI.
Figura 6. Diagrama de cuerpo libre del sistema
Fuente: Elaboración Propia
La determinación de las fuerzas magnéticas necesarias del rodamiento para generar la levitación se realiza una vez establecida las cargas del eje principal. (figura 7).
Figura 7. Reacción del rodamiento magnético
Fuente: Elaboración Propia
Los resultados de los parámetros geométricos y las variables (entrada, diseño y solución) (ver tabla 3, 4, 5 y 6), con similitudes de investigaciones como la de Zapata Sánchez (2010) y Zhang & Zhu (2017):
Tabla 3.
Variables Dimensionales
Variables |
Valores |
Unidades |
Número de polos del rodamiento magnético. |
4 |
adim |
Angulo del polo estator |
25 |
Radian |
Velocidad de rotación |
3600 |
Rpm |
Longitud del magneto |
0.05 |
M |
Longitud del rotor |
0.05 |
M |
Longitud de la volante |
0.01 |
M |
Longitud del eje |
0.8 |
M |
Diámetro exterior de la volante |
0.2 |
M |
Diámetro interior de la volante |
0.015 |
M |
Límite de corriente del sistema |
2 |
A |
Relación de hierro |
0.5 |
Adim |
Densidad de corriente |
3.6 |
/ |
Radio del polo |
0.0086 |
M |
Permeabilidad del aire |
|
T*m/A |
Fuente: Elaboración Propia
Tabla 4.
Variables de ingreso
Variables de ingreso |
||
Variables |
Valores |
Unidad |
Densidad del eje y volante |
7850 |
|
Gravedad |
9.81 |
|
Densidad del rotor |
7650 |
|
Fuente: elaboración Propia
Tabla 5.
Variables de diseño
Variables de diseño |
||
Variables |
Valores |
Unidad |
Factor de espaciamiento |
0.278 |
adim |
Diámetro del rotor óptimo |
0.0172 |
m |
Masa del rotor |
0.021 |
kg |
Masa de la volante |
2.452 |
kg |
Masa del eje |
1.1098 |
kg |
Carga total sobre el eje |
3.6054 |
kg |
Área transversal del alambre |
|
m2 |
Área del polo estator |
0.000187 |
m2 |
Área transversal bobina |
0.00067 |
m2 |
Ancho del polo estator |
0.0039 |
m |
Diámetro exterior del estator |
0.038 |
m |
Radio interior del estator |
0.015 |
m |
Ángulo de la bobina |
50 |
radian |
Fuente: Elaboración Propia
Tabla 6.
Variables de solución
Variables |
Valor |
Unidad |
||
Número de vueltas |
600 |
adim |
||
Corriente en polo 01 |
1.248 |
amper |
||
Corriente en polo 02 |
1.248 |
amper |
||
Corriente en polo 03 |
1.006 |
amper |
||
Corriente en polo 04 |
0.0896 |
amper |
||
Carga ejercida por el Polo 02 |
264 |
newton |
||
Carga ejercida por el Polo 03 |
171.69 |
newton |
||
Carga ejercida por el Polo 04 |
136. 1 |
adim |
Fuente: Elaboración Propia
La figura 8 muestra el resultado del comportamiento de las fuerzas generadas por flujo eléctrico y la intensidad de corriente para cuatro tipos de configuración de bobina (bobinas de 100, 300, 450 y 600 vueltas); las intensidades mostradas son resultado de los polos 03 y 04.
Figura 8. Fuerzas generadas en los polos 03 y 04 generadas por el flujo eléctrico y la intensidad de corriente eléctrica para cuatro configuraciones de bobinas
Fuente: elaboración propia
La figura 9 muestra como resultado una fuerza máxima de 9.239 Newton generada por la bobina que de 100 vueltas al paso del flujo eléctrico.
Figura 9. Resultado de la fuerza máxima generada para una bobina con 100 vueltas de alambre de cobre calibre AWG17.
Fuente: elaboración propia
La figura 10 muestra como resultado una fuerza máxima de 83.1508 Newton generada por la bobina que de 300 vueltas al paso del flujo eléctrico.
Figura 10. Resultado de la fuerza máxima generada para una bobina con 300 vueltas de alambre de cobre calibre AWG17.
Fuente: elaboración propia
La figura 11 muestra como resultado una fuerza máxima de 187.0892 Newton generada por la bobina que de 450 vueltas al paso del flujo eléctrico.
Figura 11. Resultado de la fuerza máxima generada para una bobina con 450 vueltas de alambre de cobre calibre AWG17.
Fuente: elaboración propia
La figura 12 muestra como resultado una fuerza máxima de 332.6 Newton generada por la bobina que de 450 vueltas al paso del flujo eléctrico. Esta configuración es la más aceptable para soportar la fuerza máxima de 264 Newton que deberá soportar el rodamiento electromagnético como elemento soporte.
Figura 12. Resultado de la fuerza máxima generada para una bobina con 600 vueltas de alambre de cobre calibre AWG17
Fuente: elaboración propia
La elección de los componentes para el banco de ensayos, se obtuvieron de acuerdo a los resultados alcanzados con el software MATLAB a través de la GUI. Con este análisis se dedujo la potencia de 1 Hp y 3600 rpm del motor idóneo para la máquina. Mediante el análisis de falla por fatiga se comprobó que el diámetro del eje de transmisión es el adecuado para soporta cargas solicitadas. El análisis de resonancia realizado se usó para determinar la frecuencia natural en los diferentes modos se obtuvo mediante el uso del software SolidWorks aplicando el método de los elementos finitos comprobando la factibilidad técnica de los parámetros seleccionados a través de la GUI (ver tabla 6)
Tabla 7.
Análisis por Elementos Finitos (FEA) para obtener la respuesta de frecuencia natural.
Modo |
Velocidad (rad/seg) |
Frecuencia (Hertz) |
Periodo (seg) |
Deformación (mm) |
1 |
413.07 |
65.742 |
0.015211 |
0.185 |
2 |
852.65 |
135.7 |
0.007369 |
0.381 |
3 |
1040.2 |
165.55 |
0.0060406 |
0.464 |
4 |
1638.1 |
260.72 |
0.0038356 |
0.730 |
5 |
3070.9 |
488.75 |
0.002046 |
1.368 |
Fuente: elaboración propia
3.1. Elaboración de protocolo de laboratorio piloto para el curso de Diseño de Sistemas Electromecánicos
En la Figura 13 y 14 se muestra un extracto de una guía de laboratorio (piloto) para el curso de Diseño de Sistemas Electromecánicos. Esta guía se ve enriquecida por el uso de la GUI desarrollada en este trabajo y se adapta muy fácilmente a varias metodologías de diseño, teniendo en cuenta todas ellas incluyen el diseño paramétrico en alguna de sus etapas.
Figura 13: Extracto de Guía de Laboratorio Piloto para el curso de Diseño de Sistemas Electromecánicos
Fuente: Elaboración propia.
Figura 14 : Resultados de la GUI asociados a la Guía de Laboratorio Piloto para el curso de Diseño de Sistemas Electromecánicos.
Fuente: Elaboración propia.
Atendiendo a las necesidades académicas la GUI propuesta para realizar pruebas con diferentes rotores que giran a una velocidad máxima de 3600 rpm con una masa prueba permitirán la implementación de las aulas virtuales de los cursos como electromagnetismo, circuitos eléctricos, máquinas eléctricas y diseño de máquinas podrán cubrir sus necesidades académicas al utilizar este banco de ensayos para el estudio del comportamiento eléctrico, electromagnético y mecánico de los rotores.
La elaboración de la GUI concluyo satisfactoriamente dado que los gráficos y resultados obtenidos concuerdan razonablemente con valores de bancos de ensayos en investigación.
En uso para el diseño de bancos de ensayo la GUI elaborada brindo valores de diseño que permiten realizar una selección de componentes para compra con miras a la implementación del banco de ensayos
La utilización de la GUI facilito el desarrollo de un protocolo de laboratorio piloto para el curso de Diseño de Sistemas Electromecánicos, el cual, servirá en el futuro inmediato para el enriquecimiento de la didáctica en dicho curso.
Con el diseño paramétrico programado en la GUI y análisis frecuencial de FEA, se comprobó la integridad estructural del eje de transmisión con un análisis de falla por carga estática y fatiga, validando los resultados que se obtienen con la simulación de la GUI.
Ibrahim, D. (2010). Engineering simulation with MATLAB: improving teaching and learning effectiveness. ELSEVIER.
Ray, S., Srivastava, S. (2020). Virtualization of science education: a lesson from the COVID-19 pandemic. J Proteins Proteom, 11, 77–80. doi:10.1007/s42485-020-00038-7
R. DeMoyer , E. E. Mitchell. (1999). Use of the MATLAB graphical user interface development environment for some control system applications. IEEE.
F. Al-Badour, M. S. (2011). Vibration analysis of rotating machinery using time–frequency. ELSEVIER.
Sayed Mahmoud, I. E. (2012). Estimación de la posición y Alternativa de Control de Sistemas de Cojinetes Magnéticos. Gijón: Universidad de Oviedo.
Soto, A., & Villegas, D. (2006). Diseño de una Plataforma de Adquisición de Datos y Control para un Banco de Pruebas de Cojinetes Magnéticos Activos. Venezuela: Universidad Simón Bolivar.
Sikanen, E. (2014). Mechanical Design and analysis of modular active magnetic bearing testrig. Finlandia: Lappeenranta University of Technology
Pastor Duran, L., & Rubio Garcia, H. (2015). Prototipo de un rodamiento magnético activo.
Gaurav, K., Madhurjya Dev , C., Sivaramakrishnan , N., & Karuna, K. (2016). Design and Analysis of a Radial Active Magnetic Bearing for Vibration Control. 810-816.
Matsuda, K., & Kanemitsu, Y. (2007). En Optimal Number of Stator Poles for Compact Active Radial Magnetic Bearings.
Zapata Sánchez, P. M. (2010). Diseño e implementación de electroimán y actuador electrónico aplicado a un sistema de levitación magnética. Lima: Creative Commons.
Zhang, W., & Zhu, H. (2017). Radial Magnetic Bearings: An Overview. Physics, 1-14.
1a1 Universidad Señor de Sipán, Pimentel – Chiclayo – Perú
a Magister en Ingeniería Mecánica, jvives@crece.uss.edu.pe
1b Magister en Ingeniería Mecánica, calvarados@crece.uss.edu.pe
1c Doctor en Ingeniería Mecánica, jolortegui@crece.uss.edu.pe