Zmodyfikowana metoda Newtona-Raphsona

W wielu zagadnieniach inżynierskich równanie macierzowe Ax = b ma charakter nieliniowy, ponieważ macierz A lub wektor b zależą od niewiadomej x. W takich przypadkach rozwiązanie nie może zostać wyznaczone jednorazowo i konieczne jest zastosowanie metod iteracyjnych, takich jak metoda Newtona-Raphsona.

Klasyczna metoda Newtona-Raphsona polega na:

• przyjęciu początkowego przybliżenia rozwiązania,
• obliczeniu macierzy układu dla tej wartości,
• wyznaczeniu poprawionej wartości rozwiązania,
• ponownym przeliczeniu macierzy w każdej kolejnej iteracji.

Chociaż metoda ta charakteryzuje się bardzo dobrą zbieżnością, jej wadą jest wysoki koszt obliczeniowy wynikający z konieczności aktualizacji macierzy w każdej iteracji, szczególnie w dużych układach równań.

Idea zmodyfikowanej metody Newtona-Raphsona

W celu ograniczenia czasu obliczeń stosuje się zmodyfikowaną metodę Newtona-Raphsona. W podejściu tym:

• macierz styczna (np. macierz sztywności w analizach mechanicznych) obliczana jest tylko raz — zazwyczaj w pierwszej iteracji,
• w kolejnych iteracjach wykorzystuje się tę samą macierz bez jej aktualizacji,
• aktualizowana jest jedynie wartość wektora niewiadomych.

Takie rozwiązanie znacząco zmniejsza koszt obliczeń pojedynczej iteracji.

Zalety i ograniczenia

Zmodyfikowana metoda Newtona-Raphsona posiada zarówno zalety, jak i ograniczenia:

Zalety:

• krótszy czas pojedynczej iteracji,
• mniejsze wymagania obliczeniowe dla dużych modeli,
• prostsza implementacja w dużych układach równań.

Ograniczenia:

• wolniejsza zbieżność w porównaniu z pełną metodą Newtona-Raphsona,
• większa liczba iteracji potrzebna do osiągnięcia wymaganej dokładności,
• mniejsza skuteczność w silnie nieliniowych problemach.

Zastosowanie w analizach MES

Zmodyfikowana metoda Newtona-Raphsona jest często stosowana w analizach nieliniowych metodą elementów skończonych (MES), szczególnie w sytuacjach, gdy:

• liczba stopni swobody modelu jest bardzo duża,
• pełna aktualizacja macierzy w każdej iteracji byłaby zbyt kosztowna,
• nieliniowość problemu nie jest bardzo silna.

Pozwala to uzyskać kompromis pomiędzy dokładnością obliczeń a czasem obliczeniowym.