---

Hallo liebe Gemeinde,

ich möchte eine Instance in ABAQUS gerne um alle drei Koordinatenachsen gleichzeitig drehen. Ziel ist es, die Orientierung mit maximale vMises Spannung zu identifizieren, welche später die versagenskritische Stelle darstellt. Alles soll später über ein Python Script ablaufen, was zur Klärung jedoch zunächst von untergeordneter Bedeutung ist.
Mein Vorgehen bisher war, mittels

x = range(0, 360, 90)
y = range(0, 360, 90)
z = range(0, 360, 90)
...
for i in x:
  for j in y:
    for k in z:
      ...
      mdb.models['Model'].rootAssembly.rotate(angle=float(i), axisDirection=(1.0, 0.0, 0.0), axisPoint=(0.0, 0.0, 0.0), instanceList=('Model-1', ))  # Drehung x-Achse
      mdb.models['Model'].rootAssembly.rotate(angle=float(j), axisDirection=(0.0, 1.0, 0.0), axisPoint=(0.0, 0.0, 0.0), instanceList=('Model-1', ))  # Drehung y-Achse
      mdb.models['Model'].rootAssembly.rotate(angle=float(k), axisDirection=(0.0, 0.0, 1.0), axisPoint=(0.0, 0.0, 0.0), instanceList=('Model-1', ))  # Drehung z-Achse
      ...

meine Instance um die Achsen zu drehen, wobei die Winkel zumeist alle unterschiedlich sind. Problem hierbei ist, dass die Instance also erst um die x-Achse, von dieser Orientierung aus um y-Achse und hiervon schließlich um die z-Achse gedreht wird; nicht also um alle Achsen gleichzeitig mit dem gewünschten Winkeln.
Gibt es in ABAQUS eine Möglichkeit, die Rotation um alle Achsen gleichzeitig durchzuführen?
Ein Ansatz von mir war, die Instance um die x-Achse zu drehen und über die inverse Rotationsmatrix(x) mir die axisDirection von y auf den ursprünglichen Zustand umzurechnen. Dies hatte bislang leider keinen Erfolg, ebenso wie der Ansatz, über die Gesamt-Rotationsmatrix eine für die drei Winkel gültige axisDirection zu berechnen. Schwierigkeit bei letzterem war, dass ich aus den drei unterschiedlichen Winkel einen einzigen Winkel für obige Definition berechnen müsste. Falls es in ABAQUS keine Funktion gibt, sieht jemand mit Python eine Möglichkeit zur Bestimmung der für obige Definition benötigten Parameter?

Gruß Marcel   

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hat sich erledigt. Die Methode über Rückrechnung der axisDirection funktioniert. Nur darf man nicht die inverse Rotationsmatrix nehmen, sondern die klassische Rotationsmatrix.
Sollte jemand irgend ein ähnliches Problem haben, kann er sich auf diese Weise behelfen!

[Diese Nachricht wurde von MWicke90 am 20. Jan. 2016 editiert.]

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP