---

Guten Tag,

Ich benötige aus einer Wärestromrechnung die flächennormalen Wärmestromdichte. Bisher waren die Flächen immer entlang einer Hauptachse ausgerichtet, so dass ich einfach HFL1/2/3 verwenden konnte. Nun habe ich eine gekrümmte Fläche (zur Vorstellung: ähnlich einer Meereswelle). Für mich ist diese Berechnung nur ein kleiner Zwischenschritt in einer Versuchsauswertung, ich verfüge kaum über Erfahrung mit Abaqus.

Prinzipiell ist mein Ansatz nun:
1. Eine Liste der Knoten/Elemente der Oberfläche zu erzeugen
2. Die HFL-Werte aus dem odb lesen
3. Die Oberflächennormale zu berechnen
4. HFL1/2/3 und Flächennormale miteinander multiplizieren.

Wie setze ich das am besten mit einem Abaqus Skript um?

Ich habe verschiedene Wege versucht um an die HFL Werte zu kommen, z.B.:
----------
surf_nset = odb.rootAssembly.instances['MESSPLATTE-FK-OHNE-TE-BOH-1'].nodeSets['ENDWALL']
surf_eset = odb.rootAssembly.instances['MESSPLATTE-FK-OHNE-TE-BOH-1'].elementSets['ENDWALL']

hfle = last_frame.fieldOutputs['HFL'].getSubset(region=surf_eset)
hfln = last_frame.fieldOutputs['HFL'].getSubset(region=surf_nset)
-------

len(hfle.values) ist allerdings deutlich länger als die Elementanzahl in meinem Set (4e5 statt 6e4 die bei der Erzeugung in cae angezeigt werden).
len(hfln.values) ist 0.

Wie komme ich an die HFL Werte an der Oberfläche? Und wie komme ich dann an die Normale an diesen Stellen?

Vielen Dank im Voraus

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

---

Du könntest Folgendes probieren: Die fragliche Oberfläche mit geeigneten 2D-Elemente überziehen (wichtig: gemeinsame Knoten mit dem 3D-Netz) und nach Berechnung die Auswertung auf die Flächenelemente beschränken. Deren Normalenrichtung ist immer die 3-Richtung. Damit könntest Du Dir sämtliches Skripting ersparen (wenns funktioniert)

------------------
Pam

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

---

Vielen Dank,

ich habe vom Support eine ähnliche Antwort erhalten. Um an den normalen Wärmestrom zu gelangen ist eine Kontaktbedingung notwendig. Es sollte möglich sein Schalenelemente zu erzeugen und zwischen diesen Dummyzellen und der betroffenen Oberfläche einen Kontakt zu erzeugen.

Ich konnte die Lösung bisher nicht testen. Wenn ich dazu komme, werde ich mich nochmals mit dem Ergebnis melden.
Als kurzfristige Lösung berechne ich mir die Normale der Wand und die Größe auf Basis der HFL1/2/3 in Matlab.

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