• Univ.Prof. Dr. Joachim Schöberl
• Univ.Prof. Dr. Oszkár Bíró

• Haik Jan Silm, Dipl.-Ing.
• Richard, Prüller BSc (Student)

Fonds zur Förderung der wissenschaftlichen Forschung FWF, öffnet eine externe URL in einem neuen Fenster
Grant number: P 27028
Funding periode: from September 1, 2014 to April 30, 2018

Die Simulation von Wirbelströmen in geschichteten Eisenkernen mit der Methode der finiten Elemente (FEM) ist im Entwurf von elektrischen Maschinen und Transformatoren von großem Interesse. Der Eisenkern wird aus ferromagnetischen kornorientierten Blechen hergestellt. Die Materialeigenschaften sind anisotrop und weisen eine magnetische Hysterese auf. Die Größenordnungen sind sehr verschieden, sie variieren von den Gesamtabmessungen eines Eisenkerns von bis zu einigen Metern bis zur Dicke von einzelnen Blechen in einem Bereich von 0,2-0,3mm. Das Modellieren jedes einzelnen Blechs ist verständlicherweise keine geeignete Lösung. Viele finite Elemente (FEe) müssen in einem solchen Modell verwendet werden, was zu extrem großen nichtlinearen Gleichungssystemen führt. Daher ist die genaue Simulation der Wirbelströme und der Eisenverluste in geschichteten ferromagnetischen Kernen mit angemessener Rechnerleistung keinesfalls zufriedenstellend gelöst. Sie stellt eine der großen Herausforderungen in der Berechnung elektromagnetischer Felder dar. 

Geschichtete Eisenkerne bilden eine periodische Mikrostruktur und sind daher prädestiniert für die FEM mit Homogenisierung. Simulationen mit FEM und Homogenisierung zeigen eine Grenzschicht, ganz ähnlich wie jene in den entsprechend stark vereinfachten Modellen von solchen Kernen mit anisotropen Materialeigenschaften. Eine genaue Approximation der Grenzschicht ist wesentlich für die exakte Bestimmung der Eisenverluste, jedoch sind viele FE-Schichten dafür erforderlich, welche die gesamte Anzahl der FEe im Model beträchtlich erhöhen. Die periodische Natur der Bleche wird durch step lap joints und Lüftungsschlitze unterbrochen oder durch Schrägung gestört, was zu komplizierten Geometrien führt, welche für sich genommen schon aufwendig in der FE-Modellierung sind. 

Eine genaue Approximation mit der FEM mit Standardpolynomen im Fall von Gleichungen mit stark oszillierenden Koeffizienten, z.B. bei Materialien mit einer Mikrostruktur, und von Problemen mit einer Grenzschicht, benötigt ein sehr feines FE-Gitter. Wir werden daher neue Mehrskalen-Finite-Elemente-Methoden (MSFEM) entwickeln, um die Mikrostruktur effizient aufzulösen. Dabei wird die Standardpolynombasis mit speziellen Funktionen vergrößert, welche a priori Information im Ansatzraum berücksichtigt um feine FE Gitter zu vermeiden. Die MSFEM wird mit dem multiharmonischen Ansatz kombiniert, um den Rechenaufwand zu verringern. Ansätze werden für die Grenzschicht und die geometrischen Probleme entworfen und in die MSFEM integriert. Die Hysterese wird durch ein Preisach- Modell berücksichtigt werden. Schnelle numerische Integrationsmethoden werden für eine effiziente MSFEM entwickelt, welche die Genauigkeit der Approximation nicht beeinflussen. 

Alle MSFEM Ansätze werden in den open source hp-FEM code Netgen/NGSolve integriert. Ein Benchmark mit Messdaten und der Supercomputer VSC, um sehr große Referenzlösungen berechnen zu können, werden eine optimale Entwicklung der neuen MSFEM-Ansätze sicherstellen. Das Ziel ist es, sehr genaue numerische Lösungen zu erstellen, welche minimale Computerkapazitäten erfordern und problemlos auf einem Personal Computer laufen.

Finite element package ngsolve, öffnet eine externe URL in einem neuen Fenster for electromagnetic problems.