Civil Engineering Reference
In-Depth Information
2.2.3 Konsistente Linearisierung
Newton-Raphson-Verfahren
Ziel: Iterative Losung eines nichtlinearen Gleichungssystems der Form R ( u )= 0 .
Taylorreihenentwicklung:
R ( u )= R ( u )+ D R ( u ) · Δ u + ...
(2.31)
u = u u :Losungsvektor des nichtlinearen Gleichungssystems
u : Zuvor ermittelte Naherungslosung
Abbruch nach linearem Glied D R ( u ) · Δ u mit
D R ( u ) = grad R = d d x
= d d u
= K T
(2.32)
R = I P : Fehlkraftvektor
I :
Vektor der inneren Krafte
P :
Vektor der außeren Krafte
u k ( i , j : FHG)
k )= ∂R i
∂u j
K T :
Tangentiale Steifigkeitsmatrix mit K T ,ij ( u
Lineares Gleichungssystem:
k )
k )
K T ( u
·
Δ u =
R ( u
(2.33)
Iterationsvorschrift ( k = Iterationszahler):
k +1 = u
k u
u
(2.34)
Abbruchbedingungen:
Inkrement Δ u (Losung des LGS) hinreichend klein (Verschiebungskonvergenz)
Residuum R hinreichend klein (Kraftkonvergenz)
Energienorm“ R · Δ u hinreichend klein (kann auch entfallen)
Uberprufung der quadratischen Konvergenz (in der Nahe der Losung):
( R · Δ u ) k +1
( R · Δ u ) k 2
( R · Δ u ) k +2
( R · Δ u ) k +1
(2.35)
Beispiel :Aus zwei nichtlinearen Gleichungen bestehendes Gleichungssystem R = 0 bzw.
( R 1 ,R 2 )=(0 , 0) mit den Unbekannten (Freiheitsgraden) u 1 und u 2 :
Losung u =( u 1 , u 2 )lasst sich geometrisch als gemeinsamer Schnittpunkt der beiden
Flachen R 1 = R 1 ( u 1 ,u 2 ) und R 2 = R 2 ( u 1 ,u 2 ) mit u 1 - u 2 -Ebene interpretieren.
Zur iterativen Losung werden zunachst zwei Tangentialebenen an die durch die
Naherungslosung u
k =( u 1 ,u 2 )des k -ten Iterationsschrittes definierten Punkte
R 1 ( u 1 ,u 2 ) und R 2 ( u 1 ,u 2 ) gelegt.
Schnittpunkt dieser beiden Ebenen mit der u 1 - u 2 -Ebene liefert die verbesserte
Naherungslosung u
k +1 =( u k + 1 ,u k + 2 ) (Startwert eines erneuten Iterationsschrittes).
Search WWH ::




Custom Search