---

Hallo,

ich automatisiere im Augenblick die Analyse eines simplen Hohlzylinders unter Torsion.
Zu diesem Zweck transformiere ich meine Spannungswerte in ein zylindrisches Koordinatensystem. Leider weichen die angezeigten Spannungswerte in CAE von denen durch Python ausgelesenen Werten voneinander ab.

Bevor ich das Problem genauer schildere, würde ich gerne wissen, ob meine folgende Koordinatentransformation in Python äquivalent zur der in CAE durch "Options-Results-Transformation... etc." ist.

myOdb = openOdb('Job-1_py.odb')
myDatumCsys = myOdb.rootAssembly.datumCsyses[name_of_datum_csys]
myStep = myOdb.steps[name_of_step]
myFrame = myStep.frames[-1]
myNode_A = myOdb.rootAssembly.nodeSets[name_of_set_A] #das ist der Knoten, dessen Spannungen mich interessieren
myFieldOutput_S  = myFrame.fieldOutputs['S']
myFieldOutput_S_transformed  = myFieldOutput_S.getTransformedField(datumCsys=myDatumCsys)
mySubset_from_Set_A = myFieldOutput_S.getSubset(region=myNode_A,position=ELEMENT_NODAL)       

Die Spannungswerte finden sich dann in  mySubset_from_Set_A.values.data[0-5]

Wenn ich den Knoten in CAE probe, bekomme ich S23 = 17.2586250 .
Python liefert aber S23 = 17.2585487365723 .
Die Werte unterscheiden sich also aber der 4. Nachkommastelle. Woran könnte das liegen?

       
       

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

---

Vielleicht greift du mit Python einen ungemittelten Wert an dem Knoten ab und mit CAE den gemittelten.

Da es aber um die 4. Nachkommastelle geht, frage ich mich gerade, ob es wirklich Sinn macht sich über sowas Gedanken zu machen...

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

---

Danke für deine Antwort Mustaine.

Ich schaue mir sowohl in CAE, als auch in Python explizit ungemittelte Werte an.

Die Abweichung ist, in der Tat, nicht sehr groß. Aber ich muss so gut wie möglich an meine analytische Lösungen herankommen. Solange ich nicht weiß, woher der Fehler in der Koordinatentransformation mit Python kommt, kann ich sie nicht benutzen. Das ist nicht wissenschaftlich und in der Praxis gefährlich.

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

---

Pfff, dann hinterfrage ertmal woher die Formel für die analytische Lösung kommt und wieviele Ungenauigkeiten du in der FEM bis zu deinem Spannungsergebnis hast. Danach dürfte dir die 4. Nachkommastelle herzlich egal sein.

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

---

Es geht nicht um die größe des Fehlers, sondern um das Auftreten des Fehlers an sich. Der Fehler bedeutet, dass die Koordinatentransformation über CAE/Viewer nicht equivalent zu meiner Transformation in Python ist. Wie soll ich jetzt entscheiden, welches Ergebnis ich nehme?

Ich verstehe, dass der Fehler sehr klein ist, aber ich brauche die Genauigkeit.

Weiß niemand sonst Rat?

[Diese Nachricht wurde von maxmeister am 10. Mrz. 2016 editiert.]

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

---

Update:

Also ich habe es tatsächlich geschafft Kontakt zum Abaqus Support aufzunehmen und dieser hat mir den kleinen Unterschied in den Transformationen bestätigt. Anscheinend gibt es einen Unterschied in den Berechnungswegen bzw. in der Numerik wenn Python benutzt wird. In der Praxis ist der unterschied jedoch vernachlässigbar, zumindest im allgemeinen.

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