Liebe Gemeinschaft,
nach langer Sucherei versuche ich nun hier Hilfe zu bekommen. Ich habe eine VUHARD-Subroutine geschrieben. Diese lässt sich auch kompilieren, weswegen ich schon mal davon ausgehen kann, dass kein Syntaxfehler vorliegt.
In der VUHARD muss ich neben der Fließkurve auch die Ableitungen nach dem Umformgrad, der Dehnrate und der Temperatur hinterlegen. Die Ableitungen habe ich mit Matlab berechnet und auch mit Matlab den Fortrancode generiert. Ich kann also auch erstmal davon ausgehen, dass die Ableitungen richtig sind und ich keinen Tippfehler habe.
Wenn ich nun die Subroutine einbinde und die Simulation starten lasse, dann bekomme ich sofort (im 3 von xtausend incrementen) die Fehlermeldung: 280 von 280 highly distorted elements. Darauf hin habe ich mein Einheitensystem geprüft, die Materialdaten überprüft und bin eigentlich zu dem Entschluss gekommen, dass alles stimmt.
Folgende Anomalität konnte ich jedoch noch beobachten: Kommentiere ich die Ableitung nach dem Umformgrad aus, läuft die Simulation durch. Ob die Ergebnise plausibel sind, kann ich noch nicht sagen. Aber Abbrechen tut es schon mal nicht.
Meinem Verständnis nach würden distorted elements bedeuten, dass die Verformung viel zu groß ist. Unabahängig davon, dass nach meiner Definition der Randbedingungen die Verformung nicht zu groß sein dürfte und auch beim Auskommentieren der Ableitung nach dem Umformgrad keine distorted elements auftreten, habe ich eine Begrenzung für den Umformgrad eingeführt: Dieser darf nicht null sein und muss min. phi = 0.05 betragen. Aber auch das funktioniert nicht.
Hat irgendeiner einer Idee, wie ich diese Subroutine vernünftig debuggen kann? Ich vermute dass sich irgendwo in INF oder NAN versteckt, das ich nicht sehe oder das aus mir noch unbekannter weise entsteht...
Danke.
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP