Metoda elementów skończonych - Wyjaśnienie: Od siatki do sił w nawodnie

From governing equations to an algebraic system

A loaded bridge is governed by the equations of linear elasticity: equilibrium (∇·σ = 0), the constitutive law relating stress to strain (σ = C:ε), and the strain-displacement relation (ε = ½(∇u + (∇u)ᵀ)). Analytical solutions exist only for the simplest geometries, so FEM instead solves the equivalent weak (variational) form — integrating the equilibrium condition against a virtual displacement rather than solving the differential equation pointwise. This reformulation is what allows a computer to handle arbitrarily complex structures, from a single truss bar to a full bridge deck.

Funkcje kształtu: przybliżanie nieznanego pola

W każdym elemencie nieznane pole przemieszczenia jest przybliżone jako ważona suma funkcji kształtnych Nᵢ(x), gdzie wagi to wartości węzłowe:

u(x) ≈ Σᵢ Nᵢ(x) · uᵢ (w jednym elemencie) Kluczowe właściwości jakiegokolwiek poprawnego zestawu funkcji kształtnych: Nᵢ(xⱼ) = δᵢⱼ (1 w swoim węźle, 0 w każdym innym węźle) Σᵢ Nᵢ = 1 (podział na całość — ruch pędowy jest dokładny) Dla belki z kratownic funkcji kształtnych są po prostu liniowe wzdłuż długości elementu; ogólne elementy ciągłe 2D i 3D używają trójkątnych, kwadratowych, tetraedrycznych lub elementów o kształcie cegły funkcji kształtnych o rosnącym rzędzie wielomianu. Funkcje kształtne wyższego rzędu (p-refinement) poprawiają dokładność bez zmniejszania siatki i mogą osiągnąć wykładniczą zbieżność dla gładkich problemów.

u(x) ≈ Σᵢ Nᵢ(x) · uᵢ   (within one element)

Key properties of any valid shape function set:
  Nᵢ(xⱼ) = δᵢⱼ   (1 at its own node, 0 at every other node)
  Σᵢ Nᵢ = 1      (partition of unity — rigid-body motion is exact)

Element stiffness and assembly

Substituting shape functions into the weak form gives each element’s stiffness matrix K_e = ∫ Bᵀ C B dV, where B contains derivatives of the shape functions and C is the elasticity matrix built from Young’s modulus E and Poisson’s ratio ν. For simple bar and beam elements, this integral has an exact closed form; for general continuum elements it is evaluated numerically via Gauss quadrature. Every element’s K_e is then assembled into the global system K·u = F — for a truss with N nodes and 2 DOF per node, K has size 2N×2N and is sparse, since a node only connects to its neighbouring members.

Warunki brzegowe i rozwiązywanie

Model mostu wymaga dwóch rodzajów warunków brzegowych: Dirichlet (podstawowych) warunków, które nakazują znane przemieszczenia, takie jak podpory zablokowane lub rolki z u = 0, stosowane poprzez modyfikację wierszy i kolumn macierzy sztywności; Neumann (naturalnych) warunków, które nakazują przyłożone siły, które pojawiają się bezpośrednio we wektorze obciążenia F. Dla struktury liniowo sprężystej rozwiązuje się K·u = F raz za pomocą metody bezpośredniej (LU lub rozkład Choleszyńskiego) lub metody iteracyjnej (metoda koniu grzbietowego) — nie jest wymagana iteracja Newtona-Raphson, chyba że pokład ma ulegać dużym odkształceniom lub elementy mogą się plastycznie wyginać.

Zaufanie do wyniku: zbieżność i błąd

Metoda MES produkuje przybliżone rozwiązanie, dlatego dokładność sprawdzana jest poprzez udoskonalanie modelu: dopracowanie (h-refinement) zmniejsza rozmiar elementów (błąd dyslokacji skaluje się w przybliżeniu z h², błąd naprężenia z h dla elementów liniowych); podwyższanie rzędu wielomianu (p-refinement). Oszacowania błędów, takie jak metoda Zienkiewicz-Zhu, porównują wygładzone, odzyskane pole naprężeń z sukowymi, dyskontynującymi rozwiązaniami MES przez granice elementów, oznaczając obszary — jak połączenie w trójkącie przenoszące skoncentrowany obciążenie — gdzie należy dopracować siatkę.

Frequently asked questions

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz the simulation i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.