Multi-Material Topology Optimization with Continuous Magnetization Direction for Permanent Magnet Synchronous Reluctance Motors

Resumen Técnico: Optimización Topológica Multi-material con Dirección de Magnetización Continua para Motores PMSynRM

1. Planteamiento del Problema

Los motores de reluctancia síncrona asistidos por imanes permanentes (PMSynRM) ofrecen un par motor promedio significativamente mayor que los motores de reluctancia síncrona (SynRM) tradicionales. Sin embargo, el diseño óptimo de estos motores representa un desafío complejo de optimización topológica multi-material, donde se debe distribuir de manera óptima tres materiales dentro del rotor:

1. Material ferromagnético (hierro).
2. Aire.
3. Imanes permanentes.

Limitaciones de enfoques anteriores:

• La mayoría de las optimizaciones existentes fijan la dirección de magnetización o la limitan a un conjunto discreto de ángulos.
• Cuando se consideran direcciones continuas, el diseño final suele requerir un filtrado posterior para cumplir con restricciones de fabricación, lo que puede degradar el rendimiento.
• Los métodos paramétricos tradicionales dependen de ingenieros experimentados y modelos analíticos complejos, limitando la exploración de nuevas geometrías.

El objetivo de este estudio es desarrollar un esquema de optimización que permita una dirección de magnetización continua sin depender de la intuición del usuario ni de parametrizaciones geométricas complejas, maximizando el par motor promedio bajo restricciones de volumen.

2. Metodología

El enfoque propuesto se basa en la optimización topológica basada en densidad y se estructura en los siguientes pilares:

A. Formulación del Problema y Variables de Diseño

• Se define un dominio computacional que incluye el estator, el rotor, las bobinas y los imanes.
• Se introducen tres variables de densidad en el rotor:
  • $\rho_\nu$: Densidad para el material ferromagnético.
  • $\rho_{Mx}, \rho_{My}$: Dos variables de densidad para representar las componentes de la magnetización, permitiendo una dirección continua.
• Se utiliza un método de 4 puntos estáticos para aproximar el par promedio, evaluando el rotor en las posiciones angulares $\theta \in \{0, \pi/12, \pi/6, \pi/4\}$. Esto reduce drásticamente el tiempo de cálculo en comparación con evaluar 500 posiciones.

B. Esquemas de Interpolación de Materiales
El artículo propone y compara varios esquemas para interpolar las propiedades entre los materiales:

1. SIMP (Polinomial): El método estándar, pero con problemas de simetría.
2. Esquema de D. Lukàš: Utiliza una función arctan para penalizar materiales intermedios sin favorecer simetrías no deseadas.
3. Nueva Propuesta (Inspirada en la Derivada Topológica): Se introduce un nuevo esquema de interpolación cúbica de Hermite. Este método asegura que la derivada de la función de interpolación coincida con la derivada topológica del problema en los límites de densidad 0 y 1, mejorando la convergencia y la calidad del diseño final.

C. Transformación del Vector de Magnetización
Para manejar la dirección continua sin ambigüedades angulares (problemas en el cálculo de gradientes con coordenadas polares) ni dependencia de la magnitud máxima según la dirección (problemas con coordenadas cartesianas puras), se propone una transformación de cuadrado a disco:

• Las variables $\rho_{Mx}$ y $\rho_{My}$ (en $[0,1]$) se mapean a un disco unitario utilizando una malla elíptica.
• Esto garantiza que la magnitud máxima de magnetización sea constante independientemente del ángulo, resolviendo problemas de escalado en la optimización.

D. Tratamiento de la Rotación y Restricciones

• Mortaring de tipo Nitsche: Se aplica para permitir la rotación libre de la malla del rotor respecto a la del estator sin necesidad de remallado, mejorando la precisión y eficiencia.
• Restricciones de Volumen: Se imponen límites superiores en el volumen de hierro y de imanes mediante el marco de Lagrangiano Aumentado.
• Filtrado y Proyección: Se utiliza filtrado Helmholtz (PDE) para evitar patrones de ajedrez y métodos de proyección (función tangente hiperbólica) para obtener fronteras definidas y eliminar materiales "grises".

E. Post-procesamiento: Agrupamiento K-Means
Dado que la optimización produce direcciones de magnetización continuas (difíciles de fabricar), se aplica un algoritmo de agrupamiento K-Means en la etapa final.

• Este algoritmo agrupa los elementos del rotor en $k$ clusters basándose en su posición y ángulo de magnetización.
• Esto homogeneiza la dirección de los imanes dentro de cada cluster, cumpliendo con las restricciones de fabricación (imanes con direcciones discretas) sin sacrificar significativamente el rendimiento.

3. Resultados Clave

Los autores validaron su estrategia mediante tres conjuntos de diseños utilizando el marco de trabajo NGSolve:

1. Optimización Hierro-Aire (SynRM):

   • Se obtuvieron diseños coherentes con la literatura para motores de reluctancia con bobinado distribuido.
   • El par promedio aumentó significativamente al incrementar la fracción de volumen de hierro (hasta 1.11 N·m con 40% de hierro).

2. Optimización Multi-material (Hierro-Aire-Imanes):

   • Se demostró que los imanes se distribuyen naturalmente en las barreras de aire del rotor.
   • Comparación de rendimiento:
     • Con 30% de volumen de imanes, el par optimizado fue de 1.90 N·m.
     • Tras el post-procesamiento con K-Means ($k=5$), el par se mantuvo en 1.87 N·m, mostrando una pérdida mínima (<2%) al imponer restricciones de fabricación.
   • Se comparó un diseño con fase de corriente $\phi = 5\pi/6$ (máximo par) y otro con $\phi = 3\pi/32$, obteniendo estructuras de rotor óptimas diferentes para cada caso.

3. Precisión del Método de 4 Puntos:

   • La aproximación del par promedio usando solo 4 posiciones estáticas mostró un error inferior al 0.1% en comparación con el cálculo de 500 puntos, validando la eficiencia computacional del método.

4. Contribuciones Principales

1. Dirección de Magnetización Continua: Propuesta de un esquema de densidad que permite optimizar la dirección de los imanes de forma continua durante el proceso, evitando sesgos iniciales.
2. Nueva Función de Interpolación: Desarrollo de un esquema de interpolación inspirado en la derivada topológica, superior a los métodos SIMP y Lukàš para este contexto no lineal.
3. Transformación Geométrica: Uso de una mapeo cuadrado-disco para manejar las variables de densidad de la magnetización, eliminando ambigüedades angulares.
4. Eficiencia Computacional: Combinación del método de Nitsche para la rotación y la aproximación de 4 puntos estáticos, reduciendo drásticamente el tiempo de cálculo.
5. Estrategia de Fabricabilidad: Integración de un algoritmo K-Means en el post-procesamiento para convertir diseños teóricos continuos en diseños manufacturables con imanes de dirección discreta.

5. Significado e Impacto

Este trabajo representa un avance significativo en el diseño automático de máquinas eléctricas. Al permitir la optimización simultánea de la topología del rotor y la dirección de magnetización continua, el método descubre configuraciones que los ingenieros humanos o los métodos paramétricos no podrían imaginar.

La capacidad de generar diseños de alto rendimiento (maximizando el par) y luego filtrarlos para cumplir con restricciones de fabricación reales (mediante K-Means) cierra la brecha entre la teoría de optimización y la aplicación industrial. Además, la reducción del tiempo de cálculo mediante la aproximación de 4 puntos y el método Nitsche hace viable la aplicación de estos algoritmos complejos en entornos de ingeniería práctica.

En conclusión, el estudio demuestra que la optimización topológica multi-material con magnetización continua es una herramienta poderosa para el desarrollo de la próxima generación de motores PMSynRM más eficientes y compactos.