Práctica 5 Máquinas eléctricas I 2023-2024

PRÁCTICA 5: SIMULACIÓN ELECTROMAGNÉTICA DE MÁQUINAS

ELÉCTRICAS

1. MATERIAL NECESARIO.

Software FEMM

2. OBJETO DE LA PRÁCTICA.
El objetivo de la práctica es el estudio de fenómenos electromagnéticos en máquinas
eléctricas mediante el método de elementos finitos, utilizando la herramienta de software libre
FEMM (Finite Element Method Magnetic). Se detallan asimismo ciertas particularidades del
software FEMM necesarias para su uso, relacionándolas con conceptos teóricos del
electromagnetismo. Por tanto, esta práctica sirve como introducción al modelado y diseño de
máquinas eléctricas con el software de elementos finitos FEMM (Finite Element Method Magnetic).

3. EXPLICACIÓN TEÓRICA
El método de elementos finitos (MEF) es un método numérico de resolución de ecuaciones
diferenciales que parte de la discretización de un dominio, dividiéndolo en un número finito de
elementos geométricos que se aproxime al modelo original; en contraste con el análisis diferencial
donde se trabaja con un dominio continuo por lo que los elementos son infinitesimalmente
pequeños. Una vez modelizado la dominio, las ecuaciones diferenciales que afectan al conjunto se
aplican a cada elemento en que se haya dividido, y se resuelven obteniendo un conjunto de
soluciones particularizadas que definen el comportamiento del conjunto.

Figura 1: Proceso de discretización de un dominio continuo.

Desde principios de los sesenta se utilizan métodos numéricos asistidos por ordenador para la
resolución de problemas científicos y tecnológicos. La revolución informática de los últimos años
ha multiplicado la potencia de los ordenadores, existiendo en la actualidad numerosos paquetes de
cálculo numérico comerciales. Uno de los problemas de mayor interés en la ingeniería y en la
ciencia actual es la resolución de ecuaciones diferenciales con derivadas parciales, PDE.
Prácticamente todas las ramas de la física pueden plantearse en forma de PDE, desde problemas
mecánicos, pasando por eléctricos y magnéticos hasta problemas térmicos.

Grado en Ingeniería Eléctrica 1

Práctica 4 Máquinas eléctricas I 2023-2024

Hasta hace pocos años solo era posible la resolución analítica de los problemas de PDE con
simetría suficiente, requiriendo los casos reales simplificaciones drásticas para hacer posible su
resolución. En la actualidad, la potencia de los ordenadores modernos ha permitido un gran
desarrollo de los métodos numéricos de cálculo, por lo que existen diversas posibilidades de
resolución como diferencias finitas (finite differences), elementos finitos (finite elements) o
elementos de frontera (boundary elements).
En la actualidad, por cuestiones de eficiencia numérica y especialmente en problemas con
forma geométrica irregular, no lineales o anisótropos, se prefiere utilizar el método de elementos
finitos. Históricamente se desarrollaron primero programas de elementos finitos para la resolución
de problemas mecánicos, de cálculo de estructuras, por lo que actualmente están más desarrollados
los programas de MEF para la resolución de problemas mecánicos que los dedicados el resto de los
problemas de PDE. En el ámbito de la ingeniería eléctrica interesa especialmente la resolución de
problemas electromagnéticos, es decir, resolver las ecuaciones de Maxwell (un tipo de PDE) para
los casos particulares que interesen en cada situación. Hoy en día existen muchos paquetes
informáticos capaces de abordar el problema con unas u otras limitaciones. La elección más
adecuada entre unos y otros depende de los requerimientos concretos del problema que se pretende
resolver.
Aunque el MEF permite la resolución de problemas muy generales, en la práctica suelen
realizarse programas modulares para la resolución de problemas concretos. De esta manera, en
muchos programas comerciales el problema electromagnético se divide en módulos electrostático,
magnetostático, conducción en baja frecuencia, electromagnético de baja frecuencia y de alta
frecuencia. Esto se debe a la dificultad que plantean muchos casos a la hora de describir un
problema concreto en el marco general.
La reducción del precio de los softwares de elementos finitos comerciales, el acceso a
programas de software libre y la mejora en su facilidad de manejo, reduce el tiempo de
preparación de los modelos y permite que diferentes profesionales puedan utilizar estas técnicas en
el diseño y desarrollo de sus productos. Los programas de elementos finitos están divididos en tres
partes básicas: preprocesador, módulo de resolución y postprocesador.

• Preprocesador
En esta primera etapa se define el problema imponiendo unas condiciones. Es aquí donde se
quitarán grados de libertad posibilitando la obtención de un sistema compatible determinado,
pasando a un sistema de n ecuaciones con n incógnitas y de solución única. Las acciones a definir
son:
 Dominio: espacio físico, su definición geométrica.

 Condiciones iniciales y/o de frontera: valores iniciales y de contorno

 Otros: propiedades físicas, entorno y otros parámetros El último paso es el
mallado. Se aproxima el dominio con una red de elementos sencillos (por ejemplo,
triángulos) que cubra la totalidad del dominio, formando lo que se llama la malla.
Es decir, se pasa de una formulación continua del problema a una formulación
discreta. Cuanto menor sea el tamaño de los elementos más cercana estará la
solución obtenida por MEF a la solución exacta. Sin embargo, por muy pequeños
que sean los elementos, siempre se introduce cierto error al mallar.

2 Departamento de Ingeniería Eléctrica

Práctica 5 Máquinas eléctricas I 2023-2024

• Módulo de resolución
Una vez conformado el dominio en pequeños elementos, se formulan las ecuaciones para
cada elemento y se combinan para obtener la solución en todo el cuerpo. Existen multitud de
métodos y literatura científica sobre los procesos para estas soluciones, pero el rasgo común se
basa en aproximar la función solución mediante muchas funciones sencillas definidas a trozos.
Los métodos más utilizados se fundamentan en:
 Método de Galerkin
 Principio de Trabajos virtuales

• Postprocesador
Una vez obtenida la solución del problema, se visualizan los resultados a través de una
interfaz gráfica, lo que permite al diseñador tomar decisiones acerca de las solicitaciones críticas e
intervenir en el diseño si fuese necesario. Para interpretar los resultados habitualmente es necesario
representarlos gráficamente de diversas formas (mapas equipotenciales, representación de los
vectores de campo eléctrico o magnético, cálculo de integrales de superficie, p.e. para calcular el
flujo magnético, integrales de línea, para calcular diferencias de potencial o fuerzas
magnetomotrices, etc.)
En el documento titulado “Anexos-Teoría-electromagnética y fundamentos de MEF” se pueden
encontrar descripciones más detalladas sobre la formulación matemática de MEF y de las
ecuaciones de Maxwell involucradas en los diferentes análisis.

4. MANUAL DE SOFTWARE FEMM

Antes de realizar la práctica se van a exponer las opciones que presenta el programa en sus 3
paquetes, a saber, preprocesador, módulo de resolución y postprocesador.

4.1 Preprocesador

Al abrir el programa, lo primero que nos va a aparecer es el tipo de problema que queremos abordar.
Tenemos 4 opciones

• Magnetics Problem

• Electrostatics Problem

• Heat Flow Problem

• Current Flow Problem

Elegiremos Magnetics problem en esta práctica.

Grado en Ingeniería Eléctrica 3

Práctica 4 Máquinas eléctricas I 2023-2024

Figura 2: Elección del tipo de problema a resolver.

Después, podremos seleccionar en el menú Problem,

Figura 3: Barra de menú y botonera del preprocesador

los detalles que definen el tipo de problema:

• Problem Type: Aquí definimos el tipo de simetría del problema:

planar o axisimétrico. Elegiremos Planar
• Length Units: Seleccionamos las unidades sobre las que vamos a
referir todas las cotas. Elegiremos Millimeters
• Frequency (Hz): es la frecuencia de operación del sistema, si se
indica 0 Hz, será un problema magnetostático; para otros valores será
magnetodinámico. Según los ejercicios elegiremos un valor u otro.
• Depth: aquí seleccionamos la profundidad, únicamente si tenemos
simetría planar. Según los ejercicios elegiremos un valor u otro.
• Solver Precision: indicamos la precisión en la resolución numérica,
la mínima es 1x10-8 . Lo dejamos como está.
• Min Angle/Max Angle: Ángulo mínimo/máximo de los triángulos
que componen el mallado de la herramienta TRIANGLE. Lo
dejamos como está.
• Smart mesh on, Succ. Approx
Figura 4: Detalles de la definición del problema

4 Departamento de Ingeniería Eléctrica

Práctica 5 Máquinas eléctricas I 2023-2024

La forma de dibujar con FEMM se basa en colocar puntos distribuidos en la rejilla según
sus coordenadas y después realizar segmentos entre ellos, bien líneas rectas, bien curvas.

Figura 5: Botonera para la definición de geometría y etiquetas de materiales

También existe la posibilidad de importar el diseño desde archivos en formato .dxf (dibujos de
AutoCAD, Catia, etc.) desde la pestaña File > Import DXF. Hay que tener en cuenta que al
dibujar una geometría es necesario poner todos los puntos, ya que FEMM no reconoce los
“atajos” de edición que tiene AutoCAD, por lo que es recomendable dibujar todos y cada
uno de los puntos de la geometría. Cuidado con darle al ratón para definir puntos. No puede
haber puntos sueltos o FEMM no compilará
Para poder recrear una geometría mediante el dibujo, es necesario introducir los puntos, a
través de la tecla TAB,

Figura 6: Tabulador para poder introducir las coordenadas de la geometría

con el botón

Figura 7: Botón que ha de estar activado para la implementación de geometría

Grado en Ingeniería Eléctrica 5

Práctica 4 Máquinas eléctricas I 2023-2024

activado. Al primer punto de todos le asignaremos el origen de coordenadas (x,y)=(0,0), el

cual estará asociado a 1 esquina de la geometría por comodidad. Si nos equivocamos y creamos
por error 1 punto, siempre podemos ir al menú Edit/Undo. O seleccionar ese punto con botón
izquierdo del ratón y darle a suprimir.
Una vez tengamos los puntos necesarios, iremos a unirlos con el botón de la línea activado.

Figura 8: Botón que ha de estar activado para la unión rectilínea de puntos

Si quisiéramos trazar arcos, tendríamos que ir al símbolo correspondiente,

Figura 9: Botón que ha de estar activado para la unión curvilínea de puntos

Y clicaríamos un punto y luego el otro entre los cuales queremos formar el arco. Hecho esto,
aparecería una ventana para indicar los grados del arco:

Figura 10: Ventana para especificar el ángulo de arco de segmento entre 2 puntos

Terminada la definición de la geometría, nos quedaría nuestra pieza recreada en el plano 2D:

6 Departamento de Ingeniería Eléctrica

Práctica 5 Máquinas eléctricas I 2023-2024

Figura 11: Geometría de la pieza a estudiar finalizada

Ahora tendríamos que asignar materiales a cada una de las celdas de nuestro problema. Para
ello, pondremos etiquetas de materiales, allí donde vaya a haber contorno cerrado. Utilizaremos el
símbolo verde de la barra de botones:

Figura 11: Botón que ha de estar activado para poner etiquetas de materiales

Y el resultado nos quedaría así:

Figura 12: Geometría de la pieza a estudiar finalizada y con etiquetas de materiales sin asignar

Para este caso concreto, tendríamos 3 etiquetas distintas de materiales:

• Hierro (para el núcleo)

• Cobre (para la bobina)

• Aire (para el entrehierro y todo el espacio que rodea la pieza)

Dando clic izquierdo a cada etiqueta podremos asignar el material concreto. Pero antes, hay que
importarlos a nuestra librería de materiales. Para ello vamos al menú Properties/Materials Library
y nos sale esto:

Grado en Ingeniería Eléctrica 7

Práctica 4 Máquinas eléctricas I 2023-2024

Figura 13: Librería de materiales y carpeta de materiales del modelo

Habrá que arrastrar los materiales necesarios de nuestro problema a la carpeta Model
Materials. Sin embargo, también podemos definir un nuevo material dando clic izquierdo a la
columna de la derecha. Y nos aparecerá el siguiente cuadro:

Figura 14: Ventana de definición de nuevo material

Este cuadro lo podremos modificar también en materiales que hayamos seleccionado de la

librería. Lo primero de todo es darle un nombre coherente y después, según la especificación, bien
implementar su curva B-H con zonas de saturación o tratarlo como material lineal con una
pendiente de esa curva concreta que podamos especificar. Si decidimos tratarlo como un material
no lineal, podemos asignarle Nonlinear B-H Curve y clicar en Edit B-H curve. Así podremos
meter los datos tabulados de la curva B-H.

8 Departamento de Ingeniería Eléctrica

Práctica 5 Máquinas eléctricas I 2023-2024

Figura 15: Implementar/Editar curva B-H associada a un material

E incluso visualizarla con el botón Plot B-H curve. Por lo demás, deberemos decidir si
incluimos un valor de conductividad o si configuramos el material como laminado, con el espesor
de chapa y factor de relleno. Una vez hallamos definido todos los materiales en nuestra carpeta,
podremos ir a la pieza a asignar en cada etiqueta el tipo de material:

Figura 16: Geometría de la pieza a estudiar finalizada y con etiquetas de materiales asignadas

A la hora de asignar etiqueta de material a las secciones de la bobina, tendremos que tener en
cuenta el sentido de circulación de la corriente eléctrica. Antes, habrá que haber definido un
circuito con el valor de la corriente. Para ello vamos a Properties/ Circuits y le damos a Add new
property. Le damos un nombre coherente, por ejemplo “I” o “Coil” y le asignamos el valor
numérico en Amperes. Cuando se trata de bobinas, tendremos que pulsar la opción “series”,
indicando que tendrá una serie de espiras en serie.
Grado en Ingeniería Eléctrica 9
Práctica 4 Máquinas eléctricas I 2023-2024

Hecho esto, vamos a las etiquetas de material correspondientes a las secciones de las bobinas
y con clic izquierdo pulsamos y nos sale el siguiente cuadro:

Figura 17: Propiedades de bobinas

En block type tenemos que seleccionar el nombre del material que hemos definido en nuestra
carpeta. (Aunque en esta captura y en las siguientes pone un tipo de material concreto para los
conductores de la bobina, no quiere decir que sea el que tenéis que escoger, porque su elección es
una tarea de la práctica. Deberéis elegir el material de acuerdo a unos cálculos que se pedirán en la
hoja de resultados para no superar límites físicos).
En “In Circuit” tenemos que seleccionar el nombre del circuito que hemos definido. Y en
“Number of Turns”, el número de espiras, siendo de sentido positivo si la corriente sale. Como en
esta sección hemos definido +50, en la otra sección de la bobina tendremos que obligatoriamente
definir el opuesto, -50.
Ya solo nos queda definir las condiciones de contorno. El problema magnético se considera
un problema de frontera abierta, y por tanto, el campo magnético creado por nuestro problema iría
menguando conforme nos fuéramos alejando, pero nunca sería completamente cero. Es por ello
que necesitamos definir 1 frontera donde trunquemos el campo magnético y lo forcemos a que sea
igual a cero. Si no lo hiciéramos, al PC le demandaríamos unas capacidades infinitas de cálculo y
no podría resolver por pedirle infinitos nudos. Por tanto, hay que fijar una frontera donde nuestro
campo se anule. Elegimos una frontera circular que envuelva nuestra pieza. Para ello creamos 2
puntos y entre ellos formamos 2 arcos a +180º y a -180º. Para ello, con la tecla del arco activada,
primero pulsamos 1 punto y después el otro, marcando después 180º. Para hacer el arco contrario,
seguimos el mismo proceso pero en orden invertido y después, marcamos180º. Esto es lo que nos
quedaría:

10 Departamento de Ingeniería Eléctrica

Práctica 5 Máquinas eléctricas I 2023-2024

Figura 18: Definición de frontera geométrica

Ahora, vamos a definir la condición de contorno asociada a nuestra frontera. Para ello, vamos
al menú Properties/Boundary y clicamos en Add property. Nos sale por defecto lo siguiente:

Figura 19: Definición de condición de frontera

Con BC type, lo que seleccionamos es el tipo de condición de frontera que queremos incluir
en nuestro problema. Con “Prescribed A” tenemos la condición de campo nulo allí donde la
asignemos.
Seleccionamos ahora ambos arcos con la tecla del arco activada y clic izquierdo, y le damos a
la barra espaciadora. Nos sale la siguiente ventana

Grado en Ingeniería Eléctrica 11

Práctica 4 Máquinas eléctricas I 2023-2024

Figura 20: Asignación de condición de frontera a frontera geométrica

Le damos a OK y ya hemos impuesto esa condición de contorno. Con lo cual, hemos definido
la frontera de nuestro dominio con el círculo que hemos definido y FEMM no calculará más allá
de él.
Ya solo nos queda mallar para ir al módulo de resolución. Para ello vamos al botón:

Figura 21: Botón para mallar nuestro modelo

Y nos malla nuestro dominio, indicando el número de nodos empleados:

Figura 22: Modelo mallado y número de nodos empleados

Nuestro problema está preparado para ser resuelto.

12 Departamento de Ingeniería Eléctrica

Práctica 5 Máquinas eléctricas I 2023-2024

4.2 Módulo de resolución

Activamos el botón

Figura 23: Botón para resolver el problema

para poder resolver. Y nos sale la pantalla de ejecución de proceso, que tardará más o menos
en función de la complejidad de nuestro problema.
4.3 Postprocesador
Una vez ha completado podemos ver los resultados en el postprocesador con el botón:

Figura 24: Botón para visualizar resultados de simulación

• Valores de punto.

Si seleccionamos un punto cualquiera en nuestro dibujo, el programa nos permite

visualizar datos, como el flujo en ese punto, la intensidad de campo etc…

Figura 25: Botón para visualizar valores de punto

• Valores de línea e integrales.

Si lo que queremos es observar las propiedades de una línea del dibujo

deberemos seleccionarla para que el sistema nos represente una gráfica con la
variación de la magnitud que deseemos:

Grado en Ingeniería Eléctrica 13

Práctica 4 Máquinas eléctricas I 2023-2024

Figura 26: Botones para visualizar valores de línea e integrales

• Valores de flujo
A nivel visual e informativo este apartado es el más importante pues nos
da una imagen muy clara de cómo se reparten los valores de las diferentes
magnitudes. Podremos visualizar las líneas de flujo magnético, elegir su
número, su rango máximo etc.

Figura 27: Botones para visualizar valores de flujo

14 Departamento de Ingeniería Eléctrica

Práctica 5 Máquinas eléctricas I 2023-2024

• Valores de circuito

Con el botón

Figura 28: Botón para visualizar valores de circuito

activado, podemos comprobar resultados en los circuitos eléctricos definidos, tales

como:
 Caída de tensión
 Potencia activa consumida
 Potencia reactiva consumida

 Potencia aparente consumida

 Impedancia asociada
 Inductancia asociada

5. REALIZACIÓN DE LA PRÁCTICA
En esta práctica se van a resolver 3 casos prácticos de electromagnetismo y máquinas
eléctricas con el software FEMM, los cuales de describen a continuación. Para favorecer el
aprovechamiento de las 3 horas de práctica se facilitarán en la sesión las medidas geométricas de
cada ejercicio, pudiendo diferir para diferentes parejas de estudiantes. En esta sección se
describirán los ejercicios a nivel conceptual para preparar al alumno de cara a su resolución.
5.1 Estudio magnetostático de un electroimán accionado por una bobina excitada en DC

Los electroimanes son dispositivos electromagnéticos que están presentes en multitud de

aplicaciones dentro de la ingeniería eléctrica, desde relés, pasando por cerraduras y sistemas de
seguridad electromagnéticos hasta escáneres de resonancia magnética. Implementaremos una de
las topologías más simples que existen y abordaremos el problema como un caso
magnetoestático, es decir, con frecuencia nula, y donde todas las magnitudes estarán en DC.
Tendremos, por tanto, un electroimán accionado por una bobina excitada en DC, el cual está
compuesto por 2 piezas, una fija (la que lleva la bobina) y otra móvil separadas por un gap, de
espesor g. La bobina, al ser alimentada con una corriente IDC va a crear un campo magnético DC
que permeará por el circuito magnético de la figura y producirá en la pieza móvil una fuerza de
atracción magnética.

Figura 29: Geometría de electroimán empleada

Grado en Ingeniería Eléctrica 15
Práctica 4 Máquinas eléctricas I 2023-2024

Emplearemos este primer caso para “aprender” a manejar el software con una geometría
sencilla y a partir de cálculos analíticos podremos comprobar con el software los resultados.
Obtendremos los valores de campo B(T) en diferentes partes de la geometría y el valor de flujo, así
como la fuerza que ejerce la pieza fija sobre la móvil debido a la alimentación de la bobina con
corriente DC, y cuánta corriente es necesaria para asegurar una inducción magnética B en el gap
concreta. Se recomienda al estudiante que, antes de la práctica, trate de resolver de forma
analítica el circuito magnético y exprese las magnitudes como reluctancia de pieza fija, de pieza
móvil, reluctancia de gap, fuerza magnetomotriz, valor esperado de flujo etc…. en función de las
variables de diseño IDC, N, g…etc

5.2 Estudio magnetodinámico de un transformador monofásico de baja potencia

Implementaremos una topología de transformador monofásico de 2 columnas y abordaremos
el problema como un caso magnetodinámico de 50 Hz, es decir, con frecuencia distinta de
cero, y donde las corrientes y tensiones estarán en AC y serán sus valores de pico. Los
resultados de corrientes y tensiones también serán valores de pico. Para elegir un diámetro
coherente del conductor de las bobinas consideraremos tablas comerciales de selección de cables
normalizados AWG, según la corriente máxima que estos pueden soportar.

Figura 30: Geometría de transformador monofásico empleada

Cambiaremos las condiciones de corriente en las bobinas primarias y secundarias para poder
plantear los ensayos de vacío y cortocircuito. A partir de las simulaciones en los 2 ensayos
obtendremos resultados que nos permitirán calcular los parámetros del circuito equivalente
aproximado.
5.3 Estudio magnetodinámico de un motor asíncrono trifásico
Partiremos de una topología de máquina asíncrona trifásica con rotor de jaula de ardilla y
abordaremos el problema con la ayuda de varios casos magnetodinámicos, donde todas las
magnitudes de tensión y corriente estarán en AC y serán sus valores de pico. Esta vez
partiremos de una geometría ya implementada y nos dedicaremos a profundizar en las
simulaciones.
Se realizarán varios casos magnetodinámicos puesto que en una máquina asíncrona las
magnitudes del estator pulsan a una frecuencia diferente a la frecuencia de las magnitudes del
rotor, y la relación entre ambas está en el deslizamiento. Por tanto, no podremos analizar lo que
sucede en el estator y rotor de forma simultánea, pero sí de forma separada.

16 Departamento de Ingeniería Eléctrica

Práctica 5 Máquinas eléctricas I 2023-2024

Figura 31: Geometría de motor assíncrono empleada

De esta forma, plantearemos un ensayo de motor asíncrono al vacío donde sólo miraremos lo que
ocurre en el estator e inicialmente a la frecuencia estatórica de f1 =50 Hz. Después, analizaremos los
resultados de la máquina a otras frecuencias estatóricas f1, 40 y 60 Hz, respectivamente, y veremos
cómo afecta el cambio de esta frecuencia a los resultados estatóricos: caída de tensión en 1 fase
estatórica, pérdidas por efecto Joule en devanado estatórico, pérdidas magnéticas por corrientes
parásitas en núcleo estatórico, campo en el entrehierro.
Posteriormente, plantearemos otros casos magnetodinámicos diferentes, correspondientes a
ensayos del motor en carga, donde podamos implementar las diferentes frecuencias a la que pulsan las
magnitudes del rotor, f2, 0.1, 0.3, 0.6 y 1 Hz, correspondientes a diferentes regímenes de deslizamiento
s, cumpliendo f2=s*f1 y considerando una frecuencia estatórica fija, f1 =50 Hz. Veremos cómo afecta
el cambio de esta frecuencia a los resultados rotóricos: pérdidas por efecto Joule en el rotor, pérdidas
magnéticas por corrientes parásitas en núcleo rotórico, campo en el entrehierro, par motor,
deslizamiento, velocidad en el eje, potencia útil. Así mismo, también extrapolaremos algún resultado
estatórico como la caída de tensión en 1 fase estatórica Veremos el efecto de invertir la secuencia de
fases del estator en el signo del par.

Grado en Ingeniería Eléctrica 17