Sterowanie iteracją w analizie nieliniowej

Sterowanie iteracją w analizie nieliniowej

Powodzenie analizy nieliniowej uzależnione jest od uzyskania zbieżności iteracji w każdym z przyrostów. W przypadku pewnych modeli, charakteryzujących się nagłym wzrostem wrażliwości na deformacje w ściśle określonym momencie narastania sił, zbieżność iteracji staje się trudno osiągalna, niezależnie od przyjętego kryterium zbieżności i wielkości przyrostów. 

W takich – wydawałoby się – beznadziejnych przypadkach zastosowanie jednej z dwóch opcji zawartych w sekcji Sterowanie iteracją okna dialogowego nieliniowej analizy statycznej: Automatyczne ograniczenie przyrostu odkształcenia i Ograniczenie przyrostu przemieszczenia może być niezbędne, by przekroczyć newralgiczny punkt analizy, a tym samym – by sfinalizować obliczenia. Niniejszy artykuł przybliża działanie tych funkcjonalności.


1. Automatyczne ograniczenie przyrostu odkształcenia

Pierwszą z metod sterowania iteracją stosuje się przede wszystkim, gdy w modelu występują elementy o skokowo zmiennej sztywności – wykresy obrazujące zależność naprężenie-odkształcenie tych elementów mają charakter – jak na poniższym rysunku – np. linii łamanej.


Wraz z narastaniem obciążenia, w tego typu elementach w którejś z iteracji może dojść do przekroczenia wartości granicznej (dla danego segmentu linii wykresu) odkształcenia elementu – niebieski punkt na kolejnym obrazku. W rezultacie siła wewnętrzna staje się większa, niż wymagana do osiągnięcia równowagi sił przy λ*P. 


W odpowiedzi na tę sytuację program w kolejnej iteracji wyznaczyć musi przemieszczenie 𝛿U wykorzystując bieżącą sztywność styczną, aby zmniejszyć siłę wewnętrzną. Przemieszczenie to może okazać się jednak na tyle duże, że punkt stanu przejdzie na jeszcze inny segment charakterystyki.


W kolejnej iteracji program ponownie podejmuje próbę osiągnięcia równowagi sił, a punkt stanu przesuwa się na inny segment. W efekcie powstaje pętla, przerywana dopiero przez wyczerpanie liczby iteracji przewidzianej dla danego przyrostu. Dalsze obliczenia muszą zostać zarzucone, a zadanie nieliniowe pozostaje nierozwiązane.


Włączenie funkcjonalności automatycznego ograniczenia dla przyrostów odkształcenia pozwoli programowi rozpoznać moment przejścia przez punkt zmiany nachylenia linii funkcji naprężenie-odkształcenie.


W takim przypadku przyrost przemieszczenia jest ograniczany do wartości 𝛿Uc, mniejszej niż obliczana standardowo wartość 𝛿U.

Dzięki temu punkt stanu może zostać przeniesiony na właściwy segment wykresu, co umożliwia wyznaczenie równowagi sił. W ten sposób udało się uniknąć zawiązania pętli kolejnych iteracji. Obliczenia mogą być efektywnie kontynuowane.​

2. Ograniczenie przyrostu przemieszczenia

Drugi rodzajem sterowania iteracją to Ograniczenie przyrostów przemieszczenia. Polega ona na przyjęciu ograniczenia dla maksymalnego przemieszczenia 𝛿U. W efekcie wektor 𝛿U jest zmniejszany do wartości 𝛿Uc, tak aby jego największa składowa równała się ustalonej wartości granicznej. Dzięki temu ograniczeniu możliwe jest przerwanie nieskończonego cyklu naprzemiennych iteracji. Należy jednak pamiętać, że im mniejsza ustawiona wartość ograniczenia przemieszczenia, tym bardziej czasochłonne obliczenia.

Praktycznym przykładem efektywnego zastosowania tej opcji jest prosty model z imperfekcją przechyłową i stężeniem z prętów kratowych o nieliniowej charakterystyce "tylko rozciąganie". W pierwszej iteracji w obu prętach kratowych, z powodu obciążenia pionowego, wystąpią ujemne (ściskające) siły osiowe. Tym samym, oba pręty kratowe zostaną de facto wyłączone z pracy. W kolejnej iteracji pozbawiona stężenia konstrukcja dozna nadmiernych deformacji – brak będzie możliwości uzyskania zbieżności, a analiza zakończy się niepowodzeniem.

Rzeczywista konstrukcja wykazałaby jednak większą sztywność, osiągając stan równowagi z lekkim przechyleniem. Ów stan równowagi będzie mógł zostać znaleziony w drodze analizy nieliniowej po uwzględnieniu opcji ograniczenia przyrostów przemieszczenia. Odpowiednio dobierając wielkość maksymalnego przyrostu, możliwe będzie wyznaczenie wartości odkształceń na tyle małych, by w modelu numerycznym wyłonił się układ z jednym z dwóch prętów kratowych zdolnym do efektywnej pracy (rozciąganym). 

Na poniższej grafice przedstawiono dwie ścieżki rozwoju deformacji układu: bez ograniczenia przyrostów przemieszczenia (ścieżka górna) i z ograniczeniem przyrostów przemieszczenia (ścieżka dolna).


    • Related Articles

    • Obliczanie ugięć elementu drewnianego w AxisVM

      1. Wprowadzenie Kwestia wyznaczenia przemieszczenia (najczęściej w postaci ugięcia) w elemencie drewnianym np. krokwi czy płatwi, wbrew pozorom nie jest zadaniem banalnym. Na końcowe przemieszczenie elementu składa się przemieszczenie chwilowe ...
    • Brak sztywności w węźle

      Jeżeli komunikat taki zostanie wyświetlony, to uruchomiona analiza układu nie będzie mogła zostać przeprowadzona, ze względu na błędne zamodelowanie warunków brzegowych i powstanie układu geometrycznie zmiennego. Jedną z najczęstszych przyczyn jest ...
    • Umieszczenie plików licencji w folderze

      Wraz z licencją AxisVM dostarczany jest plik o rozszerzeniu key. Plik ten jest niezbędny, aby uruchomić program główny AxisVM, ponieważ zapisana jest w nim wykupiona konfiguracja. W przypadku rozbudowy AxisVM o kolejne moduły lub przejście na nową ...
    • Jakie moduły mam dostępne w swojej konfiguracji?

      Aby sprawdzić jakie moduły są dostępne w danej konfiguracji programu Kliknij w menu główne Pomoc > O programie ... Wyświetlone zostanie okienko, w którym zakupione (odblokowane) moduły oznaczone są czarnym tekstem. Pozostałe moduły oraz normy ...
    • Czy do sprawdzenia przebicia w stopie fundamentowej konieczny jest moduł RC3 (Sprawdzanie przebicia w płytach)?

      Sprawdzenie przebicia w stopie fundamentowej wykonywane jest niezależnie w ramach modułu RC4 Wymiarowanie fundamentów bezpośrednich. Oznacza to, że nie jest konieczne posiadanie modułu RC3 Sprawdzanie przebicia od słupów i ścian, który weryfikuje ...