Diese Arbeit beschäftigt sich mit der Simulation und Formoptimierung von elektromechanischen Energiewandlern unter Unsicherheit. Genauer wird eine Asynchronmaschine betrachtet, deren elektromagnetischen Felder durch die magnetoquasistatische Approximation der Maxwell Gleichungen beschrieben werden können, welche mit Netzwerkgleichungen für den Käfigläufer und für den anregenden Dreiphasenstrom gekoppelt werden. Komplettiert wird das Zustandssystem durch eine Bewegungsgleichung, welche durch das elektromagnetische Drehmoment angeregt wird. Dies führt auf ein System von partiell differential-algebraischen Gleichungen. Zur numerischen Lösung der Zustandsgleichung wird ein Finite Elemente Ansatz mit einem Zeitschrittverfahren verwendet. Wir betrachten Unsicherheiten in Material und Geometrie der Maschine und verwenden einen Worst-Case-Ansatz, um diesen Unsicherheiten zu begegnen. Dies führt auf ein zweistufiges Optimierungsproblem. Da diese Probleme numerisch schwierig zu lösen sind, verwenden wir Approximationen bis zur zweiten Ordnung als Ersatzmodell. Insbesondere verwenden wir Taylormodelle in Kombination mit einer adaptiven Strategie zur Verbesserung der Approximationsgüte und ableitungsfreie Interpolationsmodelle, die ebenfalls iterativ verbessert werden können. Die Diskretisierung von partiellen Differentialgleichungen führt auf Systeme mit vielen Freiheitsgraden. Zusätzlich erhöht die Betrachtung von Unsicherheiten in der Optimierung den Berechnungsaufwand. Um unsere Berechnungen zu beschleunigen, verwenden wir Techniken zur Modellreduktion.