Software: FEM - Tutorial - Magnetfeld - Wirbelfeld-Ansatz
Am Beispiel einer 2D-Magnetfeld-Berechnung soll verdeutlicht werden, was sich hinter der Begriffs- und Formelwelt verbirgt, mit welcher man bei der FEM-Simulation konfrontiert wird.
Gewünschte Simulationsergebnisse
In einem abgegrenzten Gebiet der x,y-Ebene sollen die magnetische Induktion B(x,y) und die elektrische Feldstärke E(x,y) als grundlegende Größen berechnet werden:
- Die magnetische Induktion B(x,y) ist ein Vektor in der x,y-Ebene
- Die elektrische Feldstärke E(x,y) ist ein Vektor in z-Richtung (Senkrecht zur Induktion!)
- Alle anderen elektrischen und magnetischen Größen lassen sich über Abhängigkeiten aus B und E berechnen.
Maxwellsche Gleichungen
Der Wirbelfeld-Ansatz berücksichtigt die Verkopplung zwischen dem elektrischen und dem magnetischen Feld über die ersten beiden Maxwellschen Gleichungen sowie die Quelleneigenschaften von magn. und elektrischem Feld:
- Durchflutungsgesetz (rot H = J + δD / δt)
Das Umlaufintegral der magn. Feldstärke (MMK) ist gleich der Summe der von diesem Umlauf umfassten Ströme, dies verdeutlicht die Integralform (Gesetz von Stokes). - Induktionsgesetz (rot E = -δB / δt)
Das Umlaufintegral der elektrischen Feldstärke ist gleich der von diesem Umlauf umfassten Flussänderungsgeschwindigkeit (EMK) - Quellenfreiheit des Magnetfeldes (div B=0)
Der in eine Hülle "hineinströmende" Fluss ist gleich dem aus dieser Hülle "herausströmenden" Fluss. - Elektrische Ladungsquelle (div D=ρ)
Die Differenz aus der in eine Hülle hineinströmenden Ladung und der herausströmenden Ladung entspricht der Ladungsänderung des umhüllten Volumens.
Magnetisches Vektorpotential
Das Berechnungsproblem lässt sich numerisch vereinfachen, wenn man es auf die Berechnung des so genannten magnetischen Vektorpotentials reduziert:
- Für beliebige Wirbelfelder (die immer quellenfrei sind!) lässt sich ein "allgemeines" Vektorpotential A definieren (div rot A=0).
- Für magnetische Felder gilt wegen der Quellenfreiheit div B=0.
- Daraus schlussfolgerte Maxwell B=rot A, womit der Zusammenhang zwischen "konkreter" Induktion und "abstraktem" Vektorpotential hergestellt ist.
- Da A als zu berechnende Unbekannte in unserem 2D-Beispiel nur eine Komponente in z-Richtung besitzt, reduziert sich der Rechenaufwand für das Gleichungssystem gewaltig.
- A(x,y) kann man sich als Lösungsfläche über der betrachteten x,y-Ebene vorstellen, wenn man die Länge z des Vektors A(x,y) als Höhenwert über dem Punkt (x,y) aufträgt.
- Die Neigung dieser Fläche entspricht in jedem Punkt (x,y) der dort vorhandenen Induktion B(x,y):
- Das ergibt sich aus der vollständigen Koordinaten-Darstellung des Differential-Operators rot. Da A(x,y,z) nur ein Vektor in z-Richtung ist, sind die partiellen Ableitungen von Ax und Ay in alle 3 Koordinatenrichtungen gleich Null:
===>>> Hier geht es bald weiter!
Script vom vorigen Jahr siehe: http://www.ifte.de/lehre/cae/fem/06_magnet/wirbelfeld.html
