Metoda mnożników Lagrange’a – podstawy

W wielu zjawiskach fizycznych zachowanie obiektu podlega określonym ograniczeniom, które muszą być spełnione. Mogą one dotyczyć geometrii, kinematyki lub właściwości materiałowych. W rozwiązaniach analitycznych ograniczenia te można narzucić dokładnie, jednak w metodach numerycznych, takich jak metoda elementów skończonych (MES), ich ścisłe spełnienie bywa trudne.

W praktyce często stosuje się metodę kar, która narzuca ograniczenia w sposób przybliżony. Gdy jednak wymagane jest dokładne spełnienie ograniczeń, stosuje się metodę mnożników Lagrange’a.

Zasada działania metody mnożników Lagrange’a

Zasada działania metody mnożników Lagrange’a polega na wprowadzeniu do problemu dodatkowych niewiadomych – mnożników Lagrange’a – które odpowiadają za egzekwowanie nałożonych ograniczeń. Każde ograniczenie skutkuje pojawieniem się dodatkowej niewiadomej w układzie równań.

Z punktu widzenia fizycznego metoda ta:

• wprowadza wirtualne siły związane z ograniczeniami,
• dodaje iloczyn ograniczenia i odpowiadającego mu mnożnika do energii wewnętrznej układu,
• zapewnia dokładne spełnienie warunków narzuconych na model.

Konsekwencją takiego podejścia jest zwiększenie rozmiaru macierzy globalnej, co prowadzi do większych wymagań obliczeniowych.

Porównanie metody mnożników Lagrange’a z metodą kar

Metoda mnożników Lagrange’a i metoda kar różnią się zarówno pod względem dokładności, jak i efektywności obliczeniowej:

• metoda mnożników Lagrange’a
  • zapewnia dokładne spełnienie ograniczeń,
  • zwiększa liczbę niewiadomych o liczbę ograniczeń,
  • prowadzi do dłuższego czasu obliczeń,
• metoda kar (penalty method)
  • narzuca ograniczenia w sposób przybliżony,
  • nie zwiększa liczby niewiadomych,
  • jest bardziej wydajna obliczeniowo.

Z tego względu metoda mnożników Lagrange’a bywa mniej praktyczna w dużych modelach zawierających wiele ograniczeń.

Zastosowanie w analizach MES

Dobrym przykładem zastosowania metody mnożników Lagrange’a jest analiza materiałów nieściśliwych. Aby spełnić warunek nieściśliwości, konieczne jest nałożenie ograniczenia na:

• każdy element skończony,
• każdy węzeł,
• lub każdy punkt całkowania numerycznego.

Prowadzi to do znacznego wzrostu liczby niewiadomych, co może istotnie obniżyć wydajność obliczeń. W takich przypadkach częściej stosuje się metodę kar lub formulacje mieszane, które pozwalają uzyskać kompromis między dokładnością a efektywnością.