---

Hi,

Das ist mein erster Beitrag hier und gleich eine (sicher nicht ganz neue) Frage. Arbeite derzeit mit Abaqus an einigen Homogenisierungen für meine Abschlussarbeit und bin auch fast fertig. Zur Untersuchung der Spannungsverhältnisse an meinem Schalenmodell möchte ich die Mises Spannungen an bestimmten Knoten (Sets sind dort vorhanden) auslesen.
Mittels ProbeValue ist das im CAE kein Problem. Leider finde ich trotz Manual nicht den korrekten Befehl das mittels Python zu erledigen und weder rpy-File noch mit Macro wird der ProbeValue aufgezeichnet.

Zahlreiche Ansätze bringen mich nicht recht zum Ziel.

o1 = session.openOdb(name=odb_name)
o1.steps['Step1'].frames[-1].fieldOutputs['S'].values[0].mises

liefert mir einen Zahlenwert. Leider habe ich keine Ahnung zu welchem Punkt dieser gehört.

Einen bestimmten Punkt/Set wähle ich mit:
o1.steps['Step1'].frames[-1].fieldOutputs['S'].getSubset(region=o1.rootAssembly.nodeSets['K-4M']).values[0].mises
IndexError: Sequence index out of range

Generell ist mir die Bedeutung des IndexErrors klar, habe aber leider keinen Ansatz auf welchen Index er sich bezieht.

Vermutlich ist es nur ein kleiner Verständnisfehler, leider komme ich aber nicht weiter.
Die Arbeit ist quasi fertig, es fehlen nur noch die Spannungen.

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

---

Hi,
ich glaube dein erster Beitrag liefert dir zum ersten Element und zum ersten Integrationspunkt die entsprechende Mises-Spannung.
Über

o1.steps['Step1'].frames[-1].fieldOutputs['S'].values[0].elementLabel

kannst du dir die entsprechende Elementnummerierung ausgeben lassen.
Jetzt musst du nur noch herausfinden welches elementLabel das Element in deinem Set hat.
Das sollte über die getSubset funktion funktionieren.
Dann kannst du über die entsprechende value-Nummerierung auf die gewünschte Mises-Spannung zu greifen.

P.S.: Verwende doch mal die Funktion prettyPrint(...). Da kann man sehr schön die vorhandenen instanzen der Python-Objekte sehen.

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

---

Hi, erstmal Danke für den Ansatz!

Der erste Befehl liefert wirklich die Misesspannung am ersten Integrationspunkt des ersten Elements. Allerdings habe ich 4 Integrationspunkte (S4 Elemente) und 5 Sectionpoints. Wahllos den ersten Punkt abzugreifen erscheint mir etwas willkürlich. Wie komme ich denn an die anderen Werte des Elements? Und genaugenommen suche ich ja auch nach der Spannung in einem Knoten (meine Sets bestehen nur aus genau einem Knoten) und nicht einem Element.

Der prettyPrint ist schonmal ein guter Hinweis zum Weitersuchen.
Leider liefert er mir
prettyPrint(o1.steps['Step1'].frames[-1].fieldOutputs['S'].getSubset(region=o1.rootAssembly.nodeSets['K-4M']).values)
[]
Bedeutet das, dass für Set K-4M keine Daten vorliegen? Mit dem ProbeValue erhalte ich eine Misesspannung am Knoten. Irgendwo muss die doch hinterlegt sein...

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

---

Hi
hast du den prettyPrint mal nachgeschaut?
du verwendest den vielleicht nicht ganz korrekt.
Besser wäre

prettyPrint(o1.steps['Step1'].frames[-1].fieldOutputs['S'].getSubset(region=o1.rootAssembly.nodeSets['K-4M']).values,3)

Jetzt kann man bis zu 3 instanzen tief schauen. Bei einem großen Model wird die ausgabe allerdings sehr länglich. Von daher fange doch am besten so an:

prettyPrint(o1.steps['Step1'].frames[-1].fieldOutputs['S'].getSubset(region=o1.rootAssembly.nodeSets['K-4M']).values[0],3)

Dann kannst du dir die spezifischen Werte des aktuellen Integrationspunktes heraussuchen und ein Kritierum finden mit dem du dann den Integrationspunkt deines Sets identifizieren kannst.

Weiterhin glaube ich stehen nur Werte wie Mises-Spannung usw. für die Integrationspunkte der Elemente zur Verfügung.
Wenn du die auf einem Knoten brauchst musst du interpolieren. Sprich alle Elemente die an dem Knoten deines Sets angrenzen heraussuchen. Dann alle Mises-Spannungen der Integrationspunkte bestimmen und über die Formfunktionen (ABAQUS-Manual) auf den Knotenpunkt interpolieren.
Vielleicht hilft ja auch erst mal nur den Mittelwert der Integtationspunktwerte weiter.
So wäre jedenfalls mein Ansatz. Vielleicht irre ich mich ja auch aber ich denke erst der Viewer macht es möglich den Knoten Werte wie Spannungen usw. zuzuweisen und das auch nur im Modus "Banded".

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

---

Hi,

danke für die Infos. Ich sehe mittlerweile das Problem. Da ich sowieso nicht an exakten Werten, sondern bloß an dem generellen Verhalten der Spannungen auf Geometrieänderungen interessiert bin, werde ich mit den Elementen rund um meine Ecken/Kanten arbeiten und deren Spannungen an den Integrationspunkten untersuchen.

Viele Grüße,
Daniel

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