---

Hallo zusammen !

Ich habe ein Skript um nachträglich die Odb zu bearbeiten um die Triaxialität und den Lodewinkel einzufügen.
Jedoch bekomme ich einen Fehler in der letzten Zeile:

Code:

---

triax = allFields[Stress].getScalarField(invariant=PRESS)/allFields[Stress].getScalarField(invariant=MISES)       
xi = power(allFields[Stress].getScalarField(invariant=INV3)/allFields[Stress].getScalarField(invariant=MISES),3)  
                                                                                     
lode = 1-2/pi*acos(xi)

---

damit xi im Definitionsbereich zwischen -1 und 1 liegt habe ich folgende if Bedingung eingebaut.

Code:

---

if xi < -1:
  xi = -1
elif xi > 1:
  xi = 1

---

diese if Funktion scheint hier nicht zu Funktionieren.

Ich hoffe es kann mir jemand sagen, mit welchem Code ich xi im Definitionsbereich der Arcus-Cosinus Funktion halten kann.
Vielen Dank im Vorraus

[Diese Nachricht wurde von Harr am 02. Okt. 2013 editiert.]

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

---

edit_odb_Start_Holger_neu2.txt

Das Skript findet Ihr im Anhang.
Vielleicht kann mir dann jemand behilflich sein.

in Zeile 136 kommt der Fehler. (Error: Value out of range)
Also hat es wohl nichts mit dem Definitionsbereich zu tun...

Vielen Dank

[Diese Nachricht wurde von Harr am 04. Okt. 2013 editiert.]

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

---

Die Operation
(FieldOutput / FieldOutput) ist nicht vorgesehen,
siehe 9.8.2 Valid mathematical operations (Scripting User's Manual)
Der erste Fehler liegt bereits in Zeile 134 vor.

in Zeile 136 kommt der Fehler. (Error: Value out of range)
[/QUOTE]

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

---

Vielen Dank für den Hinweis adaptiveMesh!!
Ich werde dem sofort nachgehen und mich dann zurückmelden!

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

---

Also beim ersten Umwandeln der Odbs hat es Funktioniert, und die Werte sind auch i.O. was die Triaxialität angeht. (PRESS/MISES) -- (FieldOutput/FieldOutput)

Code:

---

triax = allFields[Stress].getScalarField(invariant=PRESS)/allFields[Stress].getScalarField(invariant=MISES)

---

Die Originalen Odbs habe ich jetzt aber nicht mehr und muss mit den bereits Umgewandelten weiter machen.

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

---

Also wenn ich die Zeile mit der Arcuscosinus-Funktion wie folgt ergänze funktioniert es
Wobei die komplette If-Funktion mit der Variable xi weg bleibt

Code:

---

lode = 1-2/pi*acos(xi*pi/180)

---

Ich bin verwirrt.. Abaqus rechnet in Grad und mit Python muss ich mit Radiant rechnen Oo

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

---

9.8.1:
Arguments to the trigonometric functions must be in radians.

Table 9-3:
all / float
acos(all)

Wie gesagt: Deine Division ist so für FieldOutput objects nicht vorgesehen. Acos(FieldOutput objects) geht.
Ärgerlich ist aber die unvollständige abaqus-Doku. Nicht einmaml die __doc__ Function ist für die Methoden sinnvoll belegt worden...

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