---

Hallo liebe CAD-User,

könnt ihr mir bei folgender Problemstellung helfen?

Wie implementiere ich eine Creep Subroutine?

Ich habe ein *Creep, law=USER in die Materialbeschreibung eingefügt und beim Starten des Jobs ein user=creep.f hinzugefügt. Die Subroutine wird auch ohne Fehlermeldung compiliert, aber die Rechnung läuft dann sehr schnell durch und auch im Ergebnis ist keine plastische Verformung zu erkennen. Da ich leider im Internet und im Handbuch nur die beiden bereits verwendeten Einstellungen gefunden habe, möchte ich euch einmal fragen, ob man noch etwas Spezielles beachten muss?

Danke schon einmal für eure Hilfe!

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

---

Rechnest Du *VISCO oder *STATIC ?

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

---

Graham_Walles_KriechenV2.txt

Ich rechne mit Coupled-Temp-Displacement.

Wenn ich die bereits mitgelieferten Creep_Gleichungen (zB Time-Hardening) aktiviere, rechnet Abaqus einwandfrei und es kommt zu einer plastischen Verformung.

Was mir noch aufgefallen ist, dass Abaqus die max Temperaturänderung pro Increment ignoriert und mit viel zu wenigen Incrementen die Temperaturschwankungen durchrechnet. Aber die Temperatur wird im .odb richtig angezeigt.

Ich habe einmal meine Creep-Subroutine angefügt. Es wäre eine super Hilfe, wenn ich wüsste, dass diese nicht die Fehlerursache ist!

Vielen Dank für die Hilfe!

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

---

Starte doch erst einmal ganz einfach indem Du den DECRA-Vektor mit konstanten Werten besetzt anstatt ihn zu berechnen. So kannst Du testen, ob kompilieren und linken funktioniert hat.

Die Koeffizienten Deines Kriechgesetzes kommen mir arg anders vor als die die ich üblicherweise verwende. Ist aber nur so ein Gefühl. Ist Dein Einheitensystem konsistent?

Die üblichen Kriechgesetze können Ärger bei t=0 machen, da hier die Dehnrate unendlich ausfällt. K.A., ob das zum Problem beiträgt.

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

---

Hallo Goldstein,

ich habe jetzt mal die Subroutine einfach mit dem Inhalt Decra(1)=0.00001 gefüllt. Leider entsteht auch hier kein Kriechen. Je nachdem wie groß ich den Wert wähle bekommt Abaqus ein Problem mit der Konvergenz der Kraft. Daher musste ich bei größeren Decra-Werten die Incrementstartgröße verkleinern. Aber im Ergebnis ist kein Kriechen erfolt..

In der Log-Datei wird nur die fehlende Sprachvorgabe bemängelt:

Begin Compiling Abaqus/Standard User Subroutines
Thu 14 Feb 2013 07:56:13 AM CET
perl: warning: Setting locale failed.
perl: warning: Please check that your locale settings:
LANGUAGE = (unset),
LC_ALL = (unset),
LANG = " "
    are supported and installed on your system.
perl: warning: Falling back to the standard locale ("C").
End Compiling Abaqus/Standard User Subroutines
Thu 14 Feb 2013 07:56:13 AM CET
Begin Linking Abaqus/Standard User Subroutines
Thu 14 Feb 2013 07:56:13 AM CET
perl: warning: Setting locale failed.
perl: warning: Please check that your locale settings:
LANGUAGE = (unset),
LC_ALL = (unset),
LANG = " "
    are supported and installed on your system.
perl: warning: Falling back to the standard locale ("C").
End Linking Abaqus/Standard User Subroutines
Thu 14 Feb 2013 07:56:14 AM CET
Begin Analysis Input File Processor
Thu 14 Feb 2013 07:56:14 AM CET
Run pre.exe
Thu 14 Feb 2013 07:56:17 AM CET
End Analysis Input File Processor
Begin Abaqus/Standard Analysis
Thu 14 Feb 2013 07:56:17 AM CET
Run standard.exe
Thu 14 Feb 2013 08:32:19 AM CET
End Abaqus/Standard Analysis
Abaqus JOB Creep_V10 COMPLETED

Hast du noch eine Idee woran es liegen kann? Ich habe bereits versucht mit Expliziter Rechnung voran zu kommen, aber leider auch ohne Erfolg?

Ich schätze, dass irgendwo ein falscher Haken sitzt, weiß aber nicht was alles voreingestellt werden muss um die Subroutine Creep richtig zu verlinken.

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

---

Geht abaqus denn nun durch die subroutine oder nicht ? Schreib doch mal was in eine Datei.

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

---

Wenn ich ehrlich bin, kann ich leider Fortran nicht gut genug um etwas in eine externe Datei zu schreiben.
Kannst du mir vielleicht mit einer einfachen Befehlszeile aushelfen?

Vielen Dank!

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

---

      write(*,*) "CREEP"

einbauen. Ausgabe sollte ins log-file erfolgen

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

---

Hallo Goldstein,

also die Ausgabe funktioniert. Im Logfile wurden innerhalb weniger Sekunden sehr viele Einträge geschrieben, dh die Subroutine wird richtig aufgerufen und läuft durch.

Woran könnte es denn noch liegen, dass es zu keinem Kriechen kommt?

Nachdem was ich gelesen habe, reicht die Rückgabe von Decra(1) und Decra(5)-Werten. Müssen sonst noch Werte in der Routine berechnet werden?

Vielen Dank für die Hilfe!

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

---

Ferndiagnose schwierig. Lade mal das Rechenmodell hoch (*inp, *for  ...)

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

---

Kriechplatten_Test.inp.txt Graham_Walles_Kriechen_write.f.txt

Das stimmt natürlich.

Ich habe das Inputdeck und die Subroutine angehängt.

Um den Subroutinenprozess zu überprüfen habe ich auch andere Versuchsversionen ausprobiert. Aber diese hier ist eigentlich die gewünschte Endversion. Die Parameter stimmen definitiv noch nicht, da ich noch auf die Kriechkurven warten muss. Ich habe daher ab 400°C relativ hohe Werte eingetragen um auf sicheres Kriechen zu kommen.

Die Routine läuft auch soweit durch, da sobald eine Temperatur von über 400°C in der Simulation auftritt der Write-Befehl ausgeführt wird.

Vielen Dank für die Hilfe!

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

---

Also beim zehnten Blick ist mir selbst nun ein Fehler aufgefallen. Ich habe vergessen die Variable sig von qtild abhängig zu machen. Werde mal sehen ob es nun funktioniert.

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

---

Aktuell hast du in dieser Alu-Platte ein Lastniveau von 2-3 MPa. Um Kriechen auszulösen ist das viel zu wenig.

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

---

Zitat:

---

Original erstellt von Goldstein:
Die üblichen Kriechgesetze können Ärger bei t=0 machen, da hier die Dehnrate unendlich ausfällt. K.A., ob das zum Problem beiträgt.

---

Hallo Goldstein,

ich bin ganz neu hier angemeldet (und auch überhaupt das erste mal bei einem Forum) deswegen bitte ich um etwas Nachsicht... 

Ich wollte diesen Thread nochmal in Hoffnung auf Hilfe aufgreifen, da ich hier in einer deiner Antworten den einzigen Bezug zu meinem aktuellen Problem gefunden habe.

In dem Zitierten Beitrag schreibst du, dass einige Subroutines Probleme beim Start t=0 haben, da hier das Kriechdehninkrement unendlich hoch wird. Vor dieser Schwierigkeit stehe ich nun auch:
Ich versuche grad die in Abaqus mit *Creep, Law=Strain hinterlegete Subroutine "nach zuprogrammieren" und stoße da auf etwas Schwierigkeiten am Start. Bei mir erfolgt der Abbruch wenn die Vergleichskriechdehnungen EC(1) bzw. EC(2), die zur Berechnung des Kriechdehninkrements benötigt werden zu klein sind, was auf implizite Art ja Zeitpunkten von t->0 entspricht wenn ich es richtig verstehe. Diese kleinen Vergleichskriechdehnungen führen dann zu sehr großen Spannungen und letztendlich dazu, dass die Lösungen divergieren.
Meine erste Idee war nun bei der Berechnung des Kriechdehninkrements für das erste Zeitinkrement eine kleine Vergleichskriechdehnung vorzugeben um einen von Null verschiedenen Wert zu haben. Die Berechnung funktionierte jetzt auch für das erste Zeitinkrement ganz gut, da die Kriechdehnungen dann beim zweiten Zeitschritt aber immer noch sehr gering waren führte dies hier zu unrealsitisch hohen Spannungen (QTILD) und so zum Abbruch. In meiner aktuellen Version die auch realistische Ergebnisse liefert und auch im Vergleich zur hinterlegten Law=Strain subroutine ganz gut aussieht (Differenzen von ca. 0.01%) sage ich, dass zur Berechnung des Kirechdehninkrements ein bestimmter Wert anstatt der Vergleichskriechdehnungen benutzt werden soll, solange diese unterhalb des bestimmten Wertes liegen(Dieser Wert ist so groß, dass es nicht zu unrealistisch hohen Kriechdehnungen kommt). Sobald die Vergleichskriechdehnungen diesen Wert überschreiten werden sie zur Berechnung benutzt. Da die größe dieses Wertes in Abhängigkeit von den wirkenden Spannungen stark variiert habe ich hier noch einen Bezug zu QTILD eingebaut. Die Subroutine trifft jetzt auch für Stark unterschiedliche Spannungen eigentlich sehr gut den Verlauf der
hinterlegten Subroutine und bleibt auch bei spontanen Spannungsänderungen stabil.

Mich würde jetzt halt nur sehr interessieren wie bei diesem Problem üblicherweise vorgegangen wird, ob du noch eine bessere Lösungsidee hättest oder wie dieses Problem in Abaqus gelöst ist um halt genau den Verlauf der in Abaqus hinterlegten LAW=STRAIN Subroutine zu treffen (falls das möglich ist).

Ich freue mich über jede Idee und Anregung natürlich auch von jedem anderen USER und entschuldigt, dass der Text jetzt so arg lang geworden ist!
Vielen Dank schon mal und viele Grüße
NICW

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

---

Verwende testweise mal ein sehr kleines Anfangszeitinkrement (z.B. 1E-25 oder noch weniger). I.d.R. hilft das

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

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

---

Werde ich Morgen versuchen! Vielen Dank für die Antwort 

...aber eigentlich müsste es doch auch anders gehen oder? Weil ich suche ja eigentlich den Unterschied zwischen meiner Subroutine zur Berechnung des Kriechens und der in Abaqus hinterlegten. Und mit der hinterlegten Subroutine läuft es ja auch mit größerem Start-Zeit-Inkrement...

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

---

Hallo, tut mir leid, dass ich bisher keine Rückmeldung gegeben habe! Ich habe es letzte Woche nicht mehr geschafft...

Habe es aber grad mit dem sehr kleinen Zeitinkrement versucht, aber an dem Verlauf der durch meine Subroutine berechneten Kriechdehnungen ändert eine Variation des Zeitinkrements nichts! Das heißt im natürlich auch, dass ich den Verlauf der in Abaqus hinterlegten Subroutine auch immernoch nicht genau treffe.

Ich müsste also nochmal um Ideen oder Hilfe bitten, um die Frage zu klären wie man in der Programmierung das Startproblem der Kriechgleichungen (unendlich hohes Kriechdehninkrement am Start) vernünftig lösen kann.

Vielen Dank und viele Grüße
NICW

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