Magnetfeldberechnung in Elektromotoren: FEM-Simulation von der Geometrie zum Ergebnis

Ein Permanentmagnet-Synchronmotor (PMSM) mit unerwarteter Drehmomentwelligkeit ist ein Entwurfsproblem, das analytische Näherungsformeln nicht lösen können. Genau hier setzt die Finite-Elemente-Methode (FEM) an: Sie berechnet das Magnetfeld im gesamten Motorquerschnitt numerisch, berücksichtigt nichtlineare Materialeigenschaften und bildet den realen Betrieb präzise ab. Dieser Artikel führt Sie Schritt für Schritt durch den vollständigen FEM-Workflow – von der Geometriemodellierung bis zur Auswertung der Flussdichteverteilung.

Grenzen analytischer Methoden

Klassische analytische Methoden beschreiben den Elektromotor als idealisiertes System: gleichmäßige Luftspaltgeometrie, lineare Materialien, vernachlässigte Streufelder. Für einfache Auslegungsschritte ist das ausreichend. Sobald Sie aber Nutgeometrien anpassen, Permanentmagnete in Segmente aufteilen oder Sättigungseffekte im Elektroblech berücksichtigen wollen, stoßen Näherungsformeln an strukturelle Grenzen.

Die Drehmomentwelligkeit eines PMSM entsteht aus dem Zusammenspiel von Rastmoment, Oberwellen im Luftspaltfeld und der genauen Nutgeometrie. Analytisch lässt sich dieser Effekt allenfalls grob abschätzen. Die FEM-Simulation berechnet das Magnetfeld direkt im gesamten Querschnitt, erfasst lokale Sättigungszonen im Zahnkopf und liefert den Drehmomentverlauf über eine vollständige Rotorumdrehung – mit einer Vorhersagegenauigkeit, die den Modellierungsaufwand rechtfertigt.

Physikalische Grundlagen der Magnetfeldberechnung

Maxwellsche Gleichungen und Vektorpotenzial

Die FEM-Magnetfeldberechnung basiert auf den Maxwellschen Gleichungen, konkret auf dem Durchflutungsgesetz und dem Induktionsgesetz. In der stationären Formulierung gilt: rot H = J, wobei H die magnetische Feldstärke und J die Stromdichte bezeichnet. Statt H und B direkt zu lösen, führen FEM-Codes das magnetische Vektorpotenzial A ein, definiert durch B = rot A. Diese Formulierung verringert die Anzahl der Unbekannten und vereinfacht die Implementierung von Randbedingungen.

Nichtlineare B-H-Kennlinien als Herausforderung

Elektroblech verhält sich nicht linear. Die Permeabilität μ sinkt mit steigender Flussdichte, sobald das Material in die Sättigung geht. FEM-Solver müssen die B-H-Kennlinie iterativ auswerten – bei jedem Berechnungsschritt passt sich die lokale Permeabilität an den aktuellen Feldzustand an. Das Newton-Raphson-Verfahren ist hier das Standardverfahren. Wer eine falsche oder zu grob abgetastete B-H-Kurve eingibt, erhält Ergebnisse, die im Sättigungsbereich systematisch falsch liegen. Die Qualität der Materialdaten bestimmt unmittelbar die Qualität der Ergebnisse.

Vom CAD-Modell zum FEM-Netz: Geometrie und Vernetzung

Periodische Symmetrie nutzen

Ein vollständiger Motorquerschnitt mit zwölf Nuten und acht Polen muss nicht vollständig modelliert werden. Bei einem vierpoligen Motor mit gleichmäßiger Nutverteilung reicht ein Polpaar als Modellausschnitt – Nutzen Sie diese periodische Symmetrie.. Das reduziert die Knotenanzahl um den Faktor vier und verkürzt die Rechenzeit entsprechend. Antiperiodische Randbedingungen an den Schnittflächen stellen sicher, dass das Feld am Rand korrekt fortgesetzt wird.

Vernetzungsstrategie für genaue Ergebnisse

Die Netzqualität entscheidet über die Genauigkeit der Drehmomentwelligkeit. Der Luftspalt ist die kritischste Zone: Hier wechselt das Feld von Stator zu Rotor, und das Gradientenprofil ist steil. Als Faustregel gilt mindestens drei bis fünf Elementlagen über die Luftspaltbreite. Zahnköpfe und Magnetkanten erfordern ebenfalls ein feineres Netz, da dort lokale Sättigungseffekte auftreten. Im Joch hingegen, wo die Feldgradienten deutlich kleiner sind, genügen gröbere Elemente.

Ein häufiger Anfängerfehler ist ein zu grobes Netz im Luftspalt bei gleichzeitig feinem Netz im Joch. Das Ergebnis ist ein Drehmomentverlauf mit numerischen Artefakten, der die echte Welligkeit überlagert. Prüfen Sie die Netzqualität daher anhand einer Konvergenzstudie: Verfeinern Sie das Netz schrittweise und beobachten Sie, ab welcher Netzdichte sich die Ergebnisse nicht mehr wesentlich ändern.

Materialien und Randbedingungen richtig definieren

Materialdaten: Bibliothek oder eigene Messkurven?

ANSYS Maxwell und COMSOL Multiphysics liefern ANSYS Maxwell und COMSOL Multiphysics liefern Materialbibliotheken mit B-H-Kurven für gängige Elektrobleche und Permanentmagnetwerkstoffe. mit B-H-Kurven für gängige Elektrobleche und Permanentmagnetwerkstoffe. Für Standardmaterialien wie M270-35A oder M400-50A reichen diese Daten in der Regel aus. Sobald Sie ein spezifisches Blechpaket eines bestimmten Herstellers verwenden oder Alterungseffekte berücksichtigen müssen, sind eigene Messkurven aus dem Epstein-Rahmen-Verfahren vorzuziehen.

Für Permanentmagnete bestimmen zwei Parameter das Simulationsergebnis Für Permanentmagnete legen zwei Parameter das Simulationsergebnis fest: die Remanenzflussdichte B_r und die Koerzitivfeldstärke H_c. B_r und die Koerzitivfeldstärke H_c. Beide Werte hängen von der Betriebstemperatur ab – ein NdFeB-Magnet verliert bei einer Temperaturerhöhung um 100 °C typischerweise rund 10–15 % seiner Remanenz. Wer das in der Simulation ignoriert, überschätzt das Drehmoment im Hochtemperaturbetrieb systematisch.

Randbedingungen physikalisch verstehen

Eine Dirichlet-Randbedingung legt den Wert des Vektorpotenzials A am Modellrand fest, typischerweise A = 0 an einem äußeren Eisenring. Das entspricht physikalisch einer magnetischen Abschirmung: Kein Fluss tritt durch diesen Rand aus. Eine Neumann-Randbedingung setzt die Normalkomponente von B zu null – das Feld verläuft tangential zum Rand. Periodische und antiperiodische Randbedingungen verknüpfen gegenüberliegende Schnittflächen und ermöglichen die Nutzung der Motorsymmetrie.

2D- versus 3D-Simulation: Aufwand und Genauigkeit abwägen

Die 2D-FEM-Simulation deckt den Großteil der relevanten Magnetfeldphänomene in Elektromotoren ab. Sie setzt voraus, dass das Feld im Wesentlichen in der Querschnittsebene liegt – eine Annahme, die für gerade Nuten und axial homogene Geometrien gut erfüllt ist. Je nach Netzdichte und Simulationstyp sind Rechenzeiten von wenigen Minuten bis zu einer Stunde typisch.

Drei Situationen erzwingen den Schritt zur 3D-Simulation: Schrägung der Nuten oder Magnete, Stirnstreuung an den Wicklungsköpfen und axiale Asymmetrien im Magnetkreis. Ein geschrägter PMSM lässt sich in 2D nicht korrekt abbilden, weil die axiale Versetzung der Magnetpole das Rastmoment grundlegend verändert. 3D-Modelle benötigen typischerweise zehn- bis hundertmal mehr Rechenzeit als vergleichbare 2D-Modelle. Für frühe Entwurfsphasen empfiehlt sich daher die 2D-Simulation; 3D-Analysen sollten gezielt für die Feinoptimierung eingesetzt werden.

Ergebnisse interpretieren: Flussdichte, Feldlinien und Drehmoment

Magnetische Flussdichteverteilung lesen

Die farbcodierte Darstellung der magnetischen Flussdichte B ist das wichtigste Ergebnisbild einer FEM-Magnetfeldsimulation. Bereiche mit Flussdichten oberhalb von 1,8–2,0 T zeigen Sättigung im Elektroblech an – erkennbar an der Farbskala, die ins Rot wechselt. Sättigungszonen in Zahnköpfen erhöhen die Reluktanz lokal und beeinflussen Drehmomentverlauf sowie Verluste. Zeigen sich solche Zonen im Modell, lohnt eine Geometrieanpassung des Zahns oder eine Überprüfung der Flussdichteauslegung.

Drehmomentverlauf und Welligkeit auswerten

Das Drehmoment wird in FEM-Codes über zwei Methoden berechnet: den Maxwell-Spannungstensor entlang einer Integrationsfläche im Luftspalt oder das Prinzip der virtuellen Arbeit. Beide Methoden liefern bei korrekter Vernetzung vergleichbare Ergebnisse. Der Drehmomentverlauf über eine Rotorumdrehung zeigt die Welligkeit direkt. Eine Drehmomentwelligkeit von unter 5 % gilt für viele Antriebsanwendungen als akzeptabel; für Servoantriebe in der Robotik sind Werte unter 1 % gefordert. Überschreiten FEM-Ergebnisse diesen Zielwert, bilden sie den Ausgangspunkt für geometrische Optimierungen – etwa Nutschrägung oder Magnetformgebung.

Softwarewerkzeuge im Überblick

Kriterium	ANSYS Maxwell	COMSOL Multiphysics	FEMM (Open Source)
Stärke	Transiente Analyse rotierender Maschinen	Thermisch-elektromagnetische Kopplung	Einfache 2D-Magnetfeldprobleme
Lizenzmodell	Kommerziell	Kommerziell	Kostenlos
3D-Unterstützung	Vollständig	Vollständig	Nein
Einstiegshürde	Mittel bis hoch	Mittel	Niedrig

ANSYS Maxwell bietet spezialisierte Lösungstypen für rotierende Maschinen: transiente Simulationen mit bewegtem Rotor und Frequenzbereichsanalysen für Wirbelstromverluste. COMSOL Multiphysics punktet bei Aufgaben, die elektromagnetische und thermische Physik direkt koppeln – etwa wenn Kupfer- und Eisenverluste gleichzeitig als Wärmequellen in der Wärmeströmungssimulation erscheinen sollen. FEMM ist als kostenlose Lösung für stationäre 2D-Probleme an Hochschulen weit verbreitet und eignet sich gut als Einstieg in die FEM-Magnetfeldsimulation.

Von der Simulation zur Motoroptimierung

FEM-Ergebnisse bilden nicht den Abschluss des Entwurfsprozesses, sondern seinen Ausgangspunkt. Parametrische Studien – etwa die systematische Variation der Nutöffnungsbreite oder des Magnetwinkels – lassen sich mit Skriptschnittstellen in ANSYS Maxwell und COMSOL automatisieren. Jede Geometrievariante erzeugt einen neuen Drehmomentverlauf; die Auswertung über alle Varianten zeigt, welche geometrische Änderung die Welligkeit am stärksten reduziert.

Die Kopplung mit thermischer Analyse bewertet Kupfer- und Eisenverluste gemeinsam und zeigt, wo Kühlmaßnahmen priorisiert werden müssen. Surrogatmodelle, die auf FEM-Ergebnissen trainiert werden, beschleunigen diesen Entwurfsprozess weiter: Statt jede Geometrievariante erneut zu simulieren, interpoliert das Modell zwischen bekannten Punkten im Entwurfsraum. Wer den FEM-basierten Motorentwurf effizienter gestalten will, sollte diese Kombination aus numerischer Simulation und datengetriebener Optimierung frühzeitig in den Entwicklungsprozess einbinden.

Häufig gestellte Fragen zur FEM-Magnetfeldsimulation

Wie wird ein Elektromotor für eine FEM-Magnetfeldsimulation geometrisch modelliert?

Die Geometrie wird entweder direkt in der FEM-Software gezeichnet oder als CAD-Datei importiert. Anschließend werden Symmetriebereiche identifiziert und das Modell auf ein Polpaar reduziert. Luftspalt, Nuten, Magnete und Joch erhalten separate Materialdefinitionen.

Welche physikalischen Gleichungen liegen der FEM-Magnetfeldberechnung zugrunde?

Die Grundlage bilden die Maxwellschen Gleichungen, insbesondere das Durchflutungsgesetz. FEM-Codes lösen diese Gleichungen in der Vektorpotenzialformulierung, was die Implementierung von Randbedingungen vereinfacht und die Anzahl der Unbekannten verringert.

Wie fein muss das Netz im Luftspalt eines Elektromotors sein?

Als Richtwert gelten mindestens drei bis fünf Elementlagen über die Luftspaltbreite. Eine Konvergenzstudie – schrittweise Netzverfeinerung mit Beobachtung des Drehmomentwerts – zeigt, ab welcher Netzdichte die Ergebnisse stabil bleiben.

Welche Software eignet sich am besten für die FEM-Simulation von Elektromotoren?

ANSYS Maxwell ist die erste Wahl für transiente Analysen rotierender Maschinen mit hohem Genauigkeitsanspruch. COMSOL Multiphysics bietet Vorteile bei gekoppelten thermisch-elektromagnetischen Simulationen. FEMM eignet sich als kostenloser Einstieg für stationäre 2D-Probleme.

Wie interpretiert man die Ergebnisse einer FEM-Simulation?

Die farbcodierte Flussdichteverteilung zeigt Sättigungszonen im Elektroblech. Der Drehmomentverlauf über eine Rotorumdrehung liefert die Welligkeit als Qualitätsmerkmal. Feldlinienbilder veranschaulichen den magnetischen Flussweg und helfen, Streufelder zu identifizieren.