Evolutionary Topology Optimization

TOPOLOGY OPTIMIZATION / RESEARCH ATLAS

Optimisation topologique tenant compte
du champ magnétique et des contraintes mécaniques

Explorer la répartition de l’acier et des vides dans un moteur en tenant compte du champ magnétique et des contraintes mécaniques.
Découvrez les solutions de Pareto finales obtenues par la méthode NGnet on/off et le rôle de leur géométrie.

01 / PARETO EXPLORER

Résultats pour un moteur dimensionné pour un véhicule hybride

Maximiser le couple moyen · Minimiser l’ondulation du couple
Contrainte : contrainte maximale de von Mises ≤ 300 MPa
Cliquez sur un point pour afficher sa géométrie.

Chargement de la dernière génération…
● Solutions de Pareto de la dernière génération (212 points)Solutions admissibles uniquement

Placez le focus sur un point puis appuyez sur Entrée. Tous les individus sont aussi accessibles dans la liste ci-dessous.

Individu —
■ Acier■ Aimant□ Vide
Couple moyen—Nm
Ondulation du couple—Nm
Contrainte maximale—MPa

Les distributions d’induction magnétique et de contrainte ainsi que les courbes de couple de la solution représentative figurent dans la comparaison ci-dessous.

OPTIMIZATION HISTORY

Voir l’évolution de l’optimisation

Cette vidéo de 24.5 secondes présente le front de Pareto et les géométries représentatives sur 200 générations. Une solution proche de 110 Nm est choisie à chaque génération ; il ne s’agit pas d’une lignée parent–descendant.

Les textes de la vidéo sont en japonais.

Lire ou télécharger la vidéo (MP4)

02 / METHOD

Répartition des matériaux par la méthode NGnet on/off

Au lieu d’activer ou de désactiver directement de nombreux éléments,
les poids des fonctions de base gaussiennes servent de variables de conception.

NGnet / NORMALIZED GAUSSIAN NETWORK
φ(x)=∑i=1Nwigi(x)∑i=1Ngi(x)
gi(x)=exp(−‖x−ci‖22σ2)
φ(x) ≥ 0 Acier ONφ(x) < 0 Vide OFF

100 fonctions de base gaussiennes ; paramètre de covariance 2×10⁻⁶ m². La répartition des matériaux varie dans la zone de conception, en préservant les régions fixes, notamment les aimants.

  1. 01

    Champ de conception lisse

    Les poids gaussiens déterminent la position, la forme et la connectivité des cavités. La répétition du domaine fondamental symétrique reconstruit le moteur entier.

  2. 02

    Évaluer le champ et les contraintes

    L’analyse magnétique non linéaire fournit le couple à chaque position du rotor. Des éléments P2 évaluent la contrainte maximale de von Mises due à la force centrifuge.

  3. 03

    Sélection et mise à jour par NSGA-II

    Les solutions concilient la maximisation du couple moyen et la minimisation de son ondulation crête à crête. Les violations des limites de contrainte et des autres exigences sont pénalisées.

Fondement de la méthode : Sato et al., Optimisation topologique d’un moteur IPM avec NGnet (2015). Les valeurs et géométries présentées ici proviennent des calculs de ce projet.

BOUNDARY RESOLUTION

Précision des interfaces entre matériaux par remaillage

La discrétisation géométrique influence les contraintes locales. La fonction NGnet est réévaluée aux nœuds ajoutés pour mieux localiser les interfaces, et pas seulement pour augmenter le nombre d’éléments.

01Raffiner le maillage de fond près des contours
02Réévaluer NGnet aux nouveaux nœuds

La présence de matériau est évaluée aux nœuds ajoutés pour localiser l’interface.

03Adapter le maillage aux interfaces avec MMG

Rayon de congé : 0,5 mm ; 16 segments par arc de 90°. Les interfaces liées et séparées sont préservées.

04Inspecter les petites arêtes avant l’analyse mécanique

Seules les petites arêtes libres admissibles près des extrémités liées sont rabattues, avec vérification de la connectivité et de la qualité des éléments.

Références sur le remaillage et l’évaluation des contraintes

Takenouchi, Hiruma et leurs collègues présentent une optimisation de moteur IPM avec des contraintes de résistance, de connectivité et de qualité du maillage. Des maillages conformes aux interfaces assurent la précision des analyses magnétiques et élastiques. Ce travail est cité pour son lien avec le remaillage et l’évaluation des contraintes.

K. Takenouchi, S. Hiruma, T. Mifune and T. Matsuo (2024), Topology and parameter optimization of IPM motor considering mechanical strength by stress and connection constraints, COMPEL, 43(4), 839–851.
DOI: 10.1108/COMPEL-11-2023-0551 · Dépôt de l’Université de Hokkaido

Les valeurs de validation ci-dessous et les résultats d’optimisation proviennent de ce projet, et non de l’article cité.

03 / REFERENCE TOPOLOGY

Géométrie de référence et conditions d’analyse

La référence est un moteur à aimants permanents intérieurs à 8 pôles et 48 encoches. Les aimants en V sont fixes ; la répartition de l’acier du rotor est optimisée.

La géométrie de référence à droite comporte des barrières de flux aux deux extrémités des aimants. Issue du CAO d’origine, elle a été analysée à 80 Arms et pour une longueur axiale de 50 mm.

Aimants
Disposition en V / Induction rémanente 1,25 T
Longueur active
50 mm
Analyse magnétique
Courant de phase 80 Arms · Phase 30° / 0–15°, 17 points (premier exclu)
Analyse mécanique
12 000 tr/min / Déformation plane / Éléments P2
Interfaces des aimants
Seule la grande face plaquée contre l’acier par la force centrifuge est liée ; les autres faces sont séparées et libres.

Rotor reconstruit par symétrie. Gris : acier de conception ; rouge : aimants.

04 / REFERENCE COMPARISON

Comparaison des géométries optimisée et de référence

Une solution représentative de la dernière génération (couple moyen 110.247 Nm, ondulation 4.264 Nm) a été sélectionnée
et réévaluée avec la géométrie de référence sur 64 intervalles.

La géométrie de référence s’appuie sur le maillage du modèle de moteur public de l’IEEJ et conserve les barrières de flux aux extrémités des aimants du CAO d’origine.

Comparaison sur un pôle des géométries de référence et optimisée
Réévaluation à 80 Arms, longueur axiale 50 mm, 64 intervalles
GéométrieCouple moyen
Nm
Ondulation
Nm
Contrainte maximale
MPa
Référence (avec barrières de flux)109.07017.368327.916
Géométrie optimisée110.2634.526298.944

La contrainte maximale est la contrainte de von Mises projetée aux nœuds. Les valeurs supérieures à 300 MPa sont en rouge.

Comparer l’induction magnétique, les contraintes et les courbes de couple
Induction magnétique des géométries de référence et optimisée
Induction magnétique à 80 Arms et phase 30°.
Contrainte de von Mises projetée aux nœuds des géométries de référence et optimisée
12 000 tr/min · Déformations planes · Éléments P2. Seule la longue face extérieure de chaque aimant est liée ; les autres faces sont séparées. « Max » correspond à la contrainte maximale projetée aux nœuds du tableau. La référence utilise les angles du CAO d’origine ; la géométrie optimisée comporte des congés d’usinage.
Courbes de couple à 64 points des géométries de référence et optimisée
Le couple moyen est la moyenne arithmétique ; l’ondulation est le maximum moins le minimum.

Résultats comparatifs CSV ↓ Courbes de couple à 64 points CSV ↓

Comparaison selon le courant et l’angle d’avance

Comparaison du couple de détente

Comparaison du couple de détente
Comparaison des courbes de couple à courant nul (0 Arms).

Comparaison couple–angle d’avance du courant

Comparaison couple–angle d’avance du courant
Comparaison du couple moyen à courants identiques pour des angles d’avance de 0° à 90°.

Ondulation du couple aux points de couple maximal

Ondulation du couple aux points de couple maximal
Pour chaque géométrie et courant, le point de couple maximal est choisi parmi les angles d’avance évalués. Le taux d’ondulation est le couple crête à crête divisé par le couple moyen, multiplié par 100.

Couple de détente CSV ↓ Courbes couple–avance CSV ↓ Ondulation au couple maximal CSV ↓

05 / INDUSTRY–ACADEMIA COLLABORATION

Demandes de collaboration industrie–université

Nous accueillons les demandes de collaboration sur l’optimisation topologique des moteurs et la conception tenant compte du champ magnétique et des contraintes mécaniques.

Coût indicatif du calcul dans le cloud

Dans cette optimisation, chaque individu nécessite 17 étapes d’analyse magnétique, un remaillage et une étape d’analyse mécanique. Les coûts estimés pour 200 générations sur un service cloud tel qu’AWS figurent ci-dessous.

Coûts de calcul indicatifs selon les conditions d’analyse
Individus par générationGénérationsCoût estimé du calcul
640 individus200 générationsEnviron 30 000 JPY
1 024 individus200 générationsEnviron 50 000 JPY

Ces coûts de calcul sont indicatifs. Les coûts réels dépendent des conditions d’analyse, de la configuration du cloud et des tarifs du service.