Hallo zusammen,
ich habe mir ein APDL-Script geschrieben, mit welchem ich Teilelementvolumen von SOLID 186 Elementen berechnen will.
Generelles Vorgehen ist folgendes:
Element selektieren und prüfen, wieviele Nodes ein Spannungskriterium erfüllen.
Wenn dies nicht bei allen Nodes der Fall ist, auf den jeweiligen Elementkanten den Punkt bestimmen, wo das Kriterium gerade noch erfüllt ist und Berechnung des Volumens des neuen Körpers.
Hierbei tritt nun mein Problem auf. Um die Grenzspannung zu ermitteln müßte ich, wenn ich es genau machen will, eine quadratische Ansatzfunktion f(sigma(x)) erstellen. Hierfür bräuchte ich allerdings zur Bestimmung der Koeffizienten 3 Stützstellen. Diese wären die beiden Nodes an den jeweiligen Ecken und als 3. Stützstelle hatte ich eigentlich gedacht den Spannungswert an der Midside-Node zu nehmen. Auf diesen kann man aber scheinbar in Ansys nicht direkt zugreifen. Versucht habe ich es folgendermaßen:
*get,aknode,ELEM,ak_elem,NODE,j
nodekor(j,1)=NX(aknode)
nodekor(j,2)=NY(aknode)
nodekor(j,3)=NZ(aknode)
*get,nodestr(j),NODE,aknode,S,1
Dies funktioniert allerdings nur für die Nodes 1-8, also die Eckknoten. Für die weiteren Knoten kommt eine Fehlermeldung, dass für die angegebene Node keine Spannungen ausgelesen werden können.
Als Behelfslösung habe ich das ganze jetzt erstmal für eine lineare Ansatzfunktion implementiert, was auch recht brauchbare Ergebnisse liefert, solange die Spannungsgradienten im Element nicht allzu hoch sind. Bei sehr starken Spannungsgradienten allerdings weichen die Ergebnisse jedoch je nach Mesh-Dichte stark von einander ab.
Hat jemamd zufällig eine Idee, wie man an die Spannungsergebnisse der Midside Nodes kommen könnte?
Vielen Dank schonmal,
Gruß,
Kai
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP