FEM in der Elektromagnetik: Funktionsweise, Vernetzung und Genauigkeit

Wer zum ersten Mal ANSYS Maxwell, COMSOL Multiphysics oder CST Studio Suite öffnet, stößt unweigerlich auf dieselbe Frage: Was passiert hinter der Benutzeroberfläche, wenn die Software ein Magnetfeld berechnet? Die Finite-Elemente-Methode (FEM) ist das rechnerische Fundament dieser Werkzeuge. Wer sie versteht, erkennt, warum ein zu grobes Netz im Luftspalt das berechnete Drehmoment verfälscht und warum falsch gesetzte Randbedingungen die Lösung grundlegend verändern.

Warum elektromagnetische Felder numerische Methoden brauchen

Elektromagnetische Felder variieren kontinuierlich im Raum. Sie folgen den Maxwell-Gleichungen, einem System partieller Differentialgleichungen, das elektrisches Feld, magnetisches Feld, Ladungsdichte und Stromdichte verknüpft. Analytische Lösungen existieren nur für einfachste Geometrien: die unendlich lange Spule, die ebene Welle im freien Raum, den idealen Plattenkondensator.

Sobald Sie einen realen Elektromotor mit Nuten, Wicklungen und nichtlinearem Eisenkern betrachten, versagen diese Formeln. Dasselbe gilt für eine Patch-Antenne auf einem Substrat mit Massefläche oder einen Hohlraumresonator in der MRT-Technik. Die Geometrie ist zu komplex, die Materialübergänge zu abrupt, die Randbedingungen zu unregelmäßig. FEM löst dieses Problem, indem sie den kontinuierlichen Raum in handhabbare Teilbereiche zerlegt.

Das Grundprinzip: Vom kontinuierlichen Feld zum diskreten Netz

FEM unterteilt das Simulationsgebiet in kleine Elemente. In zwei Dimensionen sind das typischerweise Dreiecke, in drei Dimensionen Tetraeder oder Hexaeder. Innerhalb jedes Elements wird das Feld durch einfache mathematische Funktionen – die sogenannten Ansatzfunktionen – approximiert. Diese Funktionen sind an den Knoten des Netzes miteinander verknüpft.

Ein anschauliches Beispiel: die Temperaturverteilung auf einer unregelmäßig geformten Metallplatte. Teilt man die Platte in kleine Dreiecke auf und nimmt in jedem Dreieck einen linearen Temperaturverlauf an, ist das Ergebnis an Stellen mit starken Gradienten ungenau – mit hinreichend vielen kleinen Dreiecken nähert man sich der Realität jedoch beliebig weit an. Die FEM in der Elektromagnetik folgt demselben Prinzip, mit dem Unterschied, dass das Feld kein Skalar, sondern ein Vektor ist, was die Wahl der Ansatzfunktionen erheblich anspruchsvoller macht.

Das globale Feldverhalten entsteht durch das Zusammensetzen aller Elementlösungen zu einem konsistenten Gesamtsystem. Jedes Element liefert einen Beitrag zur Systemmatrix, und die Lösung des resultierenden Gleichungssystems ergibt die Feldwerte an allen Knoten des Netzes.

Von Maxwell zu Matrizen: Die schwache Formulierung

Die Maxwell-Gleichungen in ihrer starken Form verlangen punktweise Gültigkeit im gesamten Simulationsgebiet. Das ist numerisch kaum handhabbar, weil es exakte Differenzierbarkeit an jedem Punkt voraussetzt – eine Bedingung, die diskrete Netze nicht erfüllen können.

Die schwache Formulierung, die der FEM zugrunde liegt, umgeht dieses Problem. Sie fordert keine punktweise, sondern nur eine integrale Übereinstimmung zwischen der approximierten Lösung und den Feldgleichungen. Technisch geschieht das durch Multiplikation der Differentialgleichungen mit Testfunktionen und anschließende Integration über das gesamte Gebiet – ein Verfahren, das nach dem russischen Mathematiker Boris Galerkin als Galerkin-Methode bezeichnet wird.

Das Ergebnis ist ein lineares Gleichungssystem der Form K · u = f, wobei K die Steifigkeitsmatrix, u der Vektor der unbekannten Feldwerte und f der Lastvektor ist. Die Lösung dieses Systems liefert die elektrischen und magnetischen Feldwerte an allen Netzknoten. Bei nichtlinearen Materialien – etwa dem sättigenden Eisenkern eines Transformators – muss dieser Prozess iterativ wiederholt werden.

Ansatzfunktionen: Wie das Feld zwischen den Knoten lebt

Die Wahl der Ansatzfunktionen bestimmt, wie genau das Feld innerhalb eines Elements approximiert wird. Lineare Ansatzfunktionen sind einfach und recheneffizient, aber bei starken Feldgradienten ungenau. Quadratische oder kubische Ansatzfunktionen erhöhen die Genauigkeit erheblich, ohne zwingend mehr Elemente zu erfordern.

Warum Kantenelemente in der Elektromagnetik die richtige Wahl sind

Für Vektorfelder in der Elektromagnetik haben sich sogenannte Nédélec-Elemente, auch Kantenelemente genannt, als besonders geeignet erwiesen. Der Grund ist physikalischer Natur: An Materialübergängen muss die Tangentialkomponente des elektrischen Feldes stetig sein, die Normalkomponente der magnetischen Flussdichte dagegen ebenfalls. Herkömmliche nodale Elemente erzwingen Stetigkeit an Knoten – was zu unphysikalischen Scheinlösungen führen kann, den sogenannten „spurious modes“.

Kantenelemente ordnen die Freiheitsgrade den Kanten statt den Knoten zu. Sie erfüllen die Tangentialstetigkeitsbedingungen automatisch und unterdrücken Scheinlösungen von Haus aus. Softwarewerkzeuge wie COMSOL Multiphysics und ANSYS Maxwell verwenden Kantenelemente standardmäßig für Hochfrequenzprobleme.

Vernetzung: Wo Genauigkeit entsteht oder verloren geht

Die Vernetzungsstrategie ist der Bereich, in dem erfahrene FEM-Anwender die größten Qualitätsunterschiede erzielen. Felder mit starken Gradienten – etwa an Materialübergängen, Kanten oder im Luftspalt eines Elektromotors – erfordern lokal verfeinerte Netze. Ein zu grobes Netz in diesen Bereichen führt zu Fehlern, die sich nicht durch eine feinere Auflösung an anderer Stelle kompensieren lassen.

Elementqualität und ihre Auswirkung auf die Lösung

Schlechte Elementformen destabilisieren die numerische Lösung. Sehr spitze Winkel in Dreiecken oder stark verzerrte Tetraeder erhöhen die Konditionszahl der Systemmatrix und können zu numerischen Instabilitäten führen. Erfahrungsgemäß sollte der minimale Innenwinkel eines Dreieckelements 20° nicht unterschreiten, und das Seitenverhältnis von Tetraederelementen sollte möglichst nahe bei 1 liegen.

Bei Hochfrequenzproblemen kommt eine weitere Anforderung hinzu: Die Elementgröße muss klein genug sein, um die Wellenlänge aufzulösen. Als Richtwert gelten mindestens 6–10 Elemente pro Wellenlänge im jeweiligen Material. Bei 2,4 GHz in Luft bedeutet das eine maximale Elementgröße von etwa 1,25 cm – in einem Substrat mit hoher Permittivität sinkt dieser Wert entsprechend.

Adaptive Netzverfeinerung als Qualitätssicherung

Moderne FEM-Software bietet adaptive Netzverfeinerung: Der Solver berechnet zunächst eine Lösung auf einem groben Netz, schätzt dann den lokalen Fehler und verfeinert das Netz gezielt dort, wo der Fehler groß ist. Dieser Prozess wiederholt sich, bis ein Konvergenzkriterium erfüllt ist. Das Verfahren ist recheneffizienter als eine gleichmäßige Verfeinerung des gesamten Netzes und führt zu einer Fehlerordnung, die mit der Polynomordnung der Ansatzfunktionen skaliert.

Wer eigene Simulationssetups aufbauen möchte, kann mit dem kostenlosen Open-Source-Tool FEMM (Finite Element Method Magnetics) Magnetostatik-Probleme berechnen und dabei den Einfluss der Netzauflösung auf das Ergebnis direkt beobachten.

Randbedingungen in der elektromagnetischen FEM

Falsch gesetzte Randbedingungen gehören zu den häufigsten Fehlerquellen in elektromagnetischen FEM-Simulationen. Perfekt elektrisch leitende Ränder (PEC) erzwingen das Verschwinden der Tangentialkomponente des elektrischen Feldes – eine physikalisch korrekte Bedingung für ideale Metallwände, aber eine Fehlerquelle, wenn sie versehentlich auf offene Gebiete angewendet wird.

Für Strahlungsprobleme – etwa bei der Antennensimulation – muss das Simulationsgebiet nach außen hin begrenzt werden, ohne dass Reflexionen entstehen. Hier kommen absorbierende Randbedingungen zum Einsatz: die Perfectly Matched Layer (PML) oder die Absorbing Boundary Condition (ABC). Eine PML ist eine künstliche Schicht am Rand des Rechengebiets, die einlaufende Wellen dämpft, ohne sie zu reflektieren. Sie bildet einen offenen Strahlungsraum nach, ohne dass das Rechengebiet unendlich groß sein muss.

Dirichlet-Randbedingungen legen den Feldwert an einem Rand fest (etwa das elektrische Potenzial auf einer Elektrode). Neumann-Randbedingungen schreiben die Normalableitung des Feldes vor – was bei symmetrischen Problemen zur Reduktion des Rechengebiets genutzt werden kann.

FEM im Vergleich: Wann ist es die richtige Wahl?

FEM ist nicht für jedes elektromagnetische Problem die beste Methode. Ein Vergleich mit den wichtigsten Alternativen hilft bei der Entscheidung:

Kriterium	FEM	FDTD	MoM
Geometrieflexibilität	Sehr hoch	Mittel	Hoch
Inhomogene Materialien	Sehr gut	Gut	Eingeschränkt
Breitbandige Analyse	Eingeschränkt	Sehr gut	Mittel
Offene Strahlungsprobleme	Gut (mit PML)	Gut (mit PML)	Sehr gut

FDTD arbeitet im Zeitbereich und eignet sich für breitbandige Probleme sowie für transiente Feldverläufe. Gekrümmte Geometrien bereiten FDTD jedoch Schwierigkeiten, weil das kartesische Gitter Treppenstufen erzeugt. Die Momentenmethode (MoM) ist für Strahlungsprobleme im freien Raum sehr kompakt, weil sie nur Oberflächenströme diskretisiert – bei inhomogenen Volumenmaterialien verliert sie diesen Vorteil. FEM ist bei komplexen Geometrien, anisotropen Materialien und frequenzbereichsbasierten Analysen klar im Vorteil.

FEM in der elektromagnetischen Entwicklungspraxis

Bei der Elektromotorauslegung berechnet FEM Wirbelstromverluste und Sättigungseffekte im Eisenkern unter realen Betriebsbedingungen. Der Luftspalt zwischen Rotor und Stator – oft nur wenige Zehntelmillimeter breit – erfordert eine besonders feine Vernetzung, weil dort die Feldgradienten am stärksten sind. Ein zu grobes Netz im Luftspalt führt zu einer systematischen Unterschätzung der Reluktanz und verfälscht das berechnete Drehmoment.

Im Antennenentwurf erlaubt FEM die vollständige Simulation von Feldverteilung, Eingangsimpedanz und Strahlungsdiagramm, bevor der erste Prototyp gefertigt wird. Bei 5G-Antennen im Millimeterwellenbereich, wo die Wellenlänge unter einem Zentimeter liegt, ist die Vernetzungsdichte besonders kritisch.

Hochfrequenzresonatoren in der Medizintechnik – etwa die Spulen von MRT-Geräten – profitieren von FEM-Analysen der Feldhomogenität und der spezifischen Absorptionsrate (SAR) im Gewebe. Das Standardwerk von Jin, „The Finite Element Method in Electromagnetics“, sowie das Buch von Volakis, Chatterjee und Kempel bieten für diese Anwendungen eine solide mathematische Grundlage.

Häufig gestellte Fragen zur FEM in der Elektromagnetik

Wie funktioniert FEM bei der Berechnung von Magnetfeldern?

FEM unterteilt das Simulationsgebiet in kleine Elemente, approximiert das Magnetfeld in jedem Element durch Ansatzfunktionen und löst das resultierende lineare Gleichungssystem. Das Ergebnis sind die Feldwerte an allen Netzknoten, aus denen Größen wie Flussdichte, Induktivität oder Wirbelstromverluste abgeleitet werden.

Warum ist die Vernetzung bei FEM so wichtig?

Die Netzauflösung bestimmt direkt, wie genau das Feld approximiert wird. In Bereichen mit starken Gradienten – Materialübergänge, Kanten, Luftspalte – führt ein zu grobes Netz zu systematischen Fehlern, die sich nicht durch feinere Netze an anderer Stelle ausgleichen lassen.

Wann ist FEM besser als FDTD?

FEM ist vorzuziehen, wenn die Geometrie komplex und gekrümmt ist, wenn inhomogene oder anisotrope Materialien vorliegen oder wenn eine frequenzbereichsbasierte Analyse genügt. FDTD ist effizienter für breitbandige Zeitbereichsprobleme mit einfacherer Geometrie.

Wie fein muss das Netz bei 2,4 GHz sein?

Als Richtwert gelten 6–10 Elemente pro Wellenlänge. Bei 2,4 GHz in Luft entspricht das einer maximalen Elementgröße von etwa 1–2 cm. In Materialien mit hoher Permittivität verkürzt sich die effektive Wellenlänge, und die Elementgröße muss entsprechend reduziert werden.