Lösungsverfahren und Algorithmus zur Laststeuerung
Ein standardmäßiger vollständiger Newton-Raphson-Algorithmus (NR) wird verwendet, um die Lösung eines nichtlinearen FEM-Problems zu finden.
Im Allgemeinen konvergiert der NR-Algorithmus oft nicht, wenn die volle Last in einem einzigen Schritt aufgebracht wird. Ein üblicher Ansatz, der auch hier verwendet wird, besteht darin, die Last sequentiell in mehreren Inkrementen aufzubringen und das Ergebnis des vorherigen Lastinkrements zu verwenden, um die Newton-Lösung eines nachfolgenden Inkrements zu starten. Zu diesem Zweck wurde ein Lastkontrollalgorithmus zusätzlich zum Newton-Raphson implementiert. Für den Fall, dass die NR-Iterationen nicht konvergieren, wird das aktuelle Lastinkrement auf die Hälfte seines Wertes reduziert und die NR-Iterationen werden erneut versucht.
Ein zweiter Zweck des Lastkontrollalgorithmus besteht darin, die kritische Last zu finden, die bestimmten "Stoppkriterien" entspricht - insbesondere der maximalen Dehnung im Beton, dem maximalen Schlupf in den Verbundelementen, der maximalen Verschiebung in den Verankerungselementen und der maximalen Dehnung in den Bewehrungsstäben. Die kritische Last wird mit Hilfe der Halbierungsmethode ermittelt. Wird das Stoppkriterium irgendwo im Modell überschritten, werden die Ergebnisse des letzten Lastinkrements verworfen und ein neues Inkrement mit der Hälfte der Größe des vorherigen berechnet. Dieser Vorgang wird so lange wiederholt, bis die kritische Last mit einer bestimmten Fehlertoleranz gefunden ist.
Für Beton wurde das Abbruchkriterium standardmäßig auf eine Druckdehnung von 5 % (d. h. etwa eine Größenordnung größer als die tatsächliche Versagensdehnung von Beton) und eine Zugdehnung von 7 % an den Integrationspunkten der Schalenelemente festgelegt. Im Zug wurde der Wert so gewählt, dass die Grenzdehnung in der Bewehrung, die ohne Berücksichtigung der Zugversteifung in der Regel bei 5 % liegt, zuerst erreicht wird. Im Druckbereich wurde der Wert aus mehreren Alternativen so gewählt, dass er groß genug ist, um die Auswirkungen der Quetschung in den Ergebnissen sichtbar werden zu lassen, aber klein genug, um nicht zu viele Probleme mit der numerischen Stabilität zu verursachen.
\[ \textsf{\textit{\footnotesize{Abb. 21\qquad Konstitutive Beziehung von Verbund- und Verankerungselementen, die für den Nachweis der Verankerungslänge verwendet werden:}}}\]
\[ \textsf{\textit{\footnotesize{((a) Scherspannungs-Schlupfverhalten eines Verbundelements; (b) Kraft-Weg-Verhalten eines Verankerungselements}}}\]
Für die Bewehrung wird das Anhaltekriterium in Form von Spannungen definiert. Da die Spannungen am Riss modelliert werden, entspricht das Kriterium auf Zug der Zugfestigkeit der Bewehrung unter Berücksichtigung des Sicherheitsbeiwerts. Derselbe Wert wird für das Kriterium in Druckrichtung verwendet.
Das Anhaltekriterium in Verbund- und Verankerungselementen ist α·δumax, wobei δumaxder maximale Schlupf ist, der in den Code-Nachweisen verwendet wird, und α = 10.