Gleichungslösermethoden bei nichtlinearen Berechnungen

Fachbeitrag zum Thema Statik und Anwendung von Dlubal Software

  • Knowledge Base

Fachbeitrag

In diesem Beitrag wird der Gleichungslöser bei einer nichtlinearen Berechnung mit einer Newton-Raphson-Iteration erläutert.

Einführung

Die Finite-Elemente-Methode kommt immer dann zum Einsatz, wenn sich mechanische Problemstellungen analytisch nicht lösen lassen. Häufig werden hierbei auch nichtlineare Effekte wie Ausfall bei Druck (geometrische Nichtlinearität), Plastifizierungen (materielle Nichtlinearität) und Kontakt- beziehungsweise kinematische Freiheitsgrade berücksichtigt. Diese Effekte, insbesondere bei nichtlinearen Materialmodellen, lassen sich über eine iterative Berechnungsmethode berücksichtigen.

FEM-Formulierung

Grundlegende Schritte einer FEM-Formulierung (weitere Informationen finden sich in [1]):

  1. Schwache Form des Gleichgewichts

    Schwache Form

    BgradδuσdVnichtlinearer Anteil = Bδu · ρbdV + Bδu · todA

    δu Virtuelle (Test-)Verschiebung
    t0 Anfangslastfaktor
    σdV Innere Kräfte
    ρbdV Volumenkräfte
    B integriertes Gebiet
  2. Umformung in Voigt Notation mit Tensor 4. Stufe

    Schwache Form in Voigt Notation

    BδεεdV = BδuTbdVVolumenkräfte + BδuT · t0dAOberflächenkräfte

    C Steifigkeitsmatrix
    δε Variation des Verzerrungszustands
    B integriertes Gebiet
    ε Dehnung
    u Verformung

    Mit dieser Notation wird im Weiteren die Näherungslösung des FE-Ansatzes bei nichtlinearen Materialien gelöst.
  3. Das Verschiebungsfeld wird dazu mit den Ansatzfunktionen multipliziert.

    Verschiebungsfeld über Ansatzfunktion

    u(x,t) = H(x)u^(t)

    u(x,t) Verschiebung über Zeit (Lastinkrement)
    H Formfunktion
    û Knotenverschiebung
  4. Ableitung der Verschiebung in die schwache Form einsetzen. Über numerische Integration werden Knotenverschiebung und im Post Processing die Spannungen und Verzerrungen über das Stoffgesetz berechnet.

Ablauf einer Newton-Raphson-Iteration

Aufgrund des nichtlinearen Materialverhaltens ändert sich die Materialmatrix C in obiger Gleichung 2 mit jedem Dehnschritt. Das Standard-Rechenverfahren zur Lösung dieses Problems ist die sogenannte Newton-Raphson-Iteration. Mit ihr wird die Funktion in einem Ausgangspunkt linearisiert. In der Iteration wird immer die Steifigkeitsmatrix C des Vorschritts verwendet. In dem linearisierten Iterationsschritt wird eine Tangente an die Nullstelle der Funktion gelegt.

Die zum Ablaufschema im obigen Bild gehörenden Gleichungen:

  1. Last in Lastschritte zerlegen.

    Last in Lastschritte

    fextt + t = fextt + f

    fext Äußere Last
    t Zeit/Lastschritt
    Δf Neues Lastinkrement
  2. Prädiktorschritt

    Prädiktorschritt

    K0t0φ = fextt+t - fint0t

    K Steifigkeitsmatrix des vorherigen Zeitschritts
    t+Δtfext


    Äußere Kräfte um eine weitere Laststufe erhöht

    0tfint

    Innere Kräfte des vorherigen Zeitschritts
    ϕ Verzerrung
  3. Im Punkt 3 Iteration des Ablaufschemas wird die Gesamtverzerrung abzüglich der plastischen Verzerrung berechnet (Korrektorschritt).

Das Ziel der iterativen Berechnung ist immer, dass die Summe der Lasten Null ergibt. Numerisch ist dies jedoch nicht möglich. Daher wird eine Abbruchrate ε definiert, bei der die Berechnung als ausreichend genau abgebrochen wird.

Abbruchrate

R = fextt + t - fintnt + t < ε

R Abbruchrate
fext äußere Kräfte
fint innere Kräfte
ε Epsilon Abbruchrate
t Zeitschritt

Im Programm kann die Abbruchrate unter den Berechnungsparametern gesteuert werden.

Im folgenden Bild wird der Ablauf einer Newton-Raphson-Iteration dargestellt. In der ersten Iteration

1. Iteration

Rt + t - F0t + t

R Abbruchrate
t Zeitschritt
F Kraft

wird die Abbruchrate R beziehungsweise ε nicht erreicht. Auch in der zweiten Iteration (rot) wird die Toleranzschranke nicht erreicht. Erst in der dritten Iteration ist der Abstand der Tangentensteifigkeit so klein, dass Konvergenz erreicht wird.

Wie bereits erwähnt, wird die Verformung während der Iteration immer weiter aufsummiert.

Zusammenfassung

Die Newton-Raphson-Iteration hat die Konsistenz- beziehungsweise Konvergenzordnung 2. Die Anzahl der "richtigen" Stellen der Iteration verdoppelt sich bei jedem Schritt. Demzufolge konvergiert eine Newton-Raphson-Iteration quadratisch und die Genauigkeit wird mit jeder Iteration größer, wenn das Verfahren konvergiert. Wenn das Verfahren jedoch nicht konvergiert, geht der Fehler ins Unendliche und die Berechnung wird abgebrochen.

Fehlerursachen sind zum Beispiel eine zu starke Steigung der Last Verformungs-Kurve und zu flache Neigung der Kurve im plastischen Bereich. Wenn die Last-Verformungs-Kurve im obigen Bild im zweiten Iterationsschritt einen zu starken Knick aufwiese, würde die Materialtangente und damit die Steifigkeitsmatrix die Steigung des elastischen Bereichs nicht korrekt abbilden. Die Steigung der Nullstelle wäre in diesem Fall für den plastischen Bereich falsch. Dies ist einer der Gründe, weshalb eine Erhöhung der Lastschritte mit einer verbesserten Konvergenz einhergeht.

Autor

Dipl.-Ing. (FH) Bastian Kuhn, M.Sc.

Dipl.-Ing. (FH) Bastian Kuhn, M.Sc.

Product Engineering & Customer Support

Herr Kuhn ist mit der Entwicklung und Qualitätssicherung im Bereich Holzbau betraut und im Kundensupport aktiv.

Schlüsselwörter

Konvergenz Lastschritt Zeitschritt

Literatur

[1]   Nackenhorst, U.: Vorlesungsskript Numerische Mechanik. Hannover: Institut für Baumechanik und Numerische Mechanik, Gottfried Wilhelm Leibniz Universität, 2013

Links

Schreiben Sie einen Kommentar...

Schreiben Sie einen Kommentar...

  • Aufrufe 1654x
  • Aktualisiert 2. Juni 2021

Kontakt

Kontakt zu Dlubal

Haben Sie Fragen oder brauchen Sie einen Rat? Kontaktieren Sie uns über unseren kostenlosen E-Mail-, Chat- bzw. Forum-Support oder nutzen Sie die häufig gestellten Fragen (FAQs) rund um die Uhr.

+49 9673 9203 0

info@dlubal.com

Online Training | German

RFEM | Grundlagen

Online-Schulung 6. Juli 2021 8:30 - 12:30 CEST

Online-Schulung | Englisch

RFEM | Grundlagen

Online-Schulung 13. Juli 2021 9:00 - 13:00 CEST

Online Training | German

Eurocode 2 | Stahlbetontragwerke nach DIN EN 1992-1-1

Online-Schulung 27. Juli 2021 8:30 - 12:30 CEST

Online-Schulung | Englisch

Eurocode 2 | Betonbauten nach DIN EN 1992-1-1

Online-Schulung 29. Juli 2021 8:30 - 12:30 CEST

Online-Schulung | Englisch

RFEM | Strukturdynamik und Erdbebenbemessung nach EC 8

Online-Schulung 11. August 2021 8:30 - 12:30 CEST

Online-Schulung | Englisch

RFEM für Studenten | USA

Online-Schulung 11. August 2021 13:00 - 16:00 EDT

Online-Schulung | Englisch

Eurocode 3 | Stahlbau nach DIN EN 1993-1-1

Online-Schulung 25. August 2021 8:30 - 12:30 CEST

Online Training | German

Eurocode 3 | Stahltragwerke nach DIN EN 1993-1-1

Online-Schulung 9. September 2021 8:30 - 12:30 CEST

Online Training | German

RSTAB | Grundlagen

Online-Schulung 16. September 2021 8:30 - 12:30 CEST

Online Training | German

Eurocode 5 | Holztragwerke nach DIN EN 1995-1-1

Online-Schulung 21. September 2021 8:30 - 12:30 CEST

Online-Schulung | Englisch

Eurocode 5 | Holzbau nach DIN EN 1995-1-1

Online-Schulung 23. September 2021 8:30 - 12:30 CEST

Online Training | German

RFEM | Baudynamik und Erdbebenbemessung nach EC 8

Online-Schulung 28. September 2021 8:30 - 12:30 CEST

Glass Design mit Dlubal Software

Glass Design mit Dlubal Software

Webinar 8. Juni 2021 14:00 - 14:45 CEST

Glasbemessung mit Dlubal-Software

Glasbemessung mit Dlubal-Software

Webinar 1. Juni 2021 14:00 - 14:45 CEST

Blast Time History Analysis in RFEM

Blast Time History Analysis in RFEM

Webinar 13. Mai 2021 14:00 - 15:00 EDT

Holzbau | Teil 2: Bemessung

Holzbalken und Flächenstrukturen | Teil 2: Bemessung

Webinar 11. Mai 2021 14:00 - 15:00 CEST

Tragwerke aus Holz | Teil 2: Bemessung

Stab- und Flächentragwerke aus Holz | Teil 2: Bemessung

Webinar 27. April 2021 14:00 - 14:45 CEST

Platten- und Schalenbeulen mit Dlubal-Software

Platten- und Schalenbeulen mit Dlubal-Software

Webinar 30. März 2021 14:00 - 14:45 CEST

Plattenbeulen und Schalenbeulen mit Dlubal-Software

Plattenbeulen und Schalenbeulen mit Dlubal-Software

Webinar 23. März 2021 14:00 - 14:45 BST

CSA S16:19 Stahlbemessung in RFEM

CSA S16: 19 Stahlbemessung in RFEM

Webinar 10. März 2021 14:00 - 15:00 EDT

Die häufigsten Anwenderfehler mit RFEM und RSTAB

Die häufigsten Anwenderfehler mit RFEM und RSTAB

Webinar 4. Februar 2021 14:00 - 15:00 BST

Die häufigsten Anwenderfehler mit RFEM und RSTAB

Die häufigsten Anwenderfehler mit RFEM und RSTAB

Webinar 28. Januar 2021 14:00 - 15:00 BST

}
RFEM
RFEM 5.xx

Basisprogramm

Das FEM-Programm RFEM ermöglicht die schnelle und einfache Modellierung und Berechnung von Tragkonstruktionen mit Stab-, Platten-, Scheiben-, Faltwerk-, Schalen- und Volumen-Elementen aus verschiedenen Materialien.

Erstlizenzpreis
3.540,00 USD