Algorytm Lanczosa – zastosowanie w analizie modalnej

Gdy na wolny koniec drewnianej belki, o określonym przekroju poprzecznym, której drugi koniec jest przymocowany do ściany, zostanie nałożone pionowe obciążenie udarowe, belka zacznie drgać w dół w określonym kształcie. Kształt tych drgań zależy od tego, jak szybko i z jaką siłą działające obciążenie zostało przyłożone. Drgania obiektu poddanego obciążeniu dynamicznemu są wynikiem oddziaływania między siłą bezwładności a siłą przywracającą wynikającą ze sztywności obiektu.

Charakterystyka drgań własnych i odwzorowanie układu

Ciało sprężyste – takie jak drewniana belka – może być traktowane jako układ dynamiczny, w którym nieskończona liczba mas punktowych jest połączona z nieskończoną liczbą sprężyn. Takie ciało, nazywane w terminologii technicznej ciągłym, ma nieskończone stopnie swobody. W związku z tym również kształt drgań własnych oraz częstotliwości drgań własnych są nieskończone. Analiza tych drgań własnych obiektu nazywana jest analizą modalną.

Jeśli jednak ciało ciągłe zostanie odwzorowane za pomocą skończonej liczby elementów skończonych, liczba możliwych drgań własnych zostaje ograniczona do liczby stopni swobody siatki, co pozwala uzyskać skończony układ dynamiczny.

Najważniejsze aspekty analizy modalnej

• liczba modów własnych w rzeczywistym obiekcie jest nieskończona,
• modele oparte na elementach skończonych redukują tę liczbę do stopni swobody siatki,
• wyznaczanie wszystkich modów w dużych modelach byłoby czasochłonne i zbędne,
• do poprawnego opisu zachowania konstrukcji zazwyczaj wystarczy kilka pierwszych modów.

W praktyce istnieje potrzeba wyznaczania tylko ograniczonej liczby modów własnych przy użyciu metod numerycznych. W przypadku dużej liczby elementów skończonych pełne obliczenia stają się czasochłonne. Na szczęście wpływ wyższych modów własnych na ogólną odpowiedź dynamiczną struktury maleje wraz z ich rzędem, dlatego w analizach uwzględnia się najczęściej tylko najniższe częstotliwości.

Rola algorytmu Lanczosa

Algorytm Lanczosa to metoda numeryczna umożliwiająca rozwiązanie układu równań macierzowych (n × n) opisujących mody własne przez sprowadzenie go do mniejszego układu (N × N), zawierającego jedynie N najniższych modów własnych. Wybór liczby N powinien być jednak przemyślany i adekwatny do wymagań analizy – zbyt małe N może prowadzić do pominięcia istotnych modów.

Algorytm został opracowany przez Corneliusa Lanczosa (1893–1974), węgierskiego matematyka, asystenta Einsteina i współtwórcę podstaw współczesnej teorii względności.