Wartość węzłowa

Wartość węzłowa to wynik analizy MES przypisany bezpośrednio do węzła siatki — np. przemieszczenie, temperatura czy potencjał. Stanowi ona podstawową wielkość obliczaną w metodzie elementów skończonych i punkt wyjścia do wyznaczania pozostałych pól, takich jak naprężenia, odkształcenia czy strumienie cieplne.

Czym jest wartość węzłowa?

W metodzie elementów skończonych rozwiązanie zadania fizycznego sprowadza się do rozwiązania równania macierzowego, którego niewiadomymi są wartości pewnej wielkości w węzłach siatki. Po jego rozwiązaniu w każdym węźle uzyskujemy konkretną wartość — i to właśnie ona nazywana jest wartością węzłową.

W zależności od rodzaju analizy wartości węzłowe mogą reprezentować różne wielkości fizyczne:

• przemieszczenia węzłowe — w analizach mechanicznych,
• temperatury węzłowe — w analizach termicznych,
• ciśnienia węzłowe — w zagadnieniach przepływowych,
• potencjały elektryczne — w analizach elektromagnetycznych,
• siły węzłowe — jako reakcje lub wartości równoważne obciążenia.

Skąd biorą się wartości węzłowe?

W MES poszukiwane pole (np. pole przemieszczeń) przybliża się za pomocą sumy funkcji kształtu pomnożonych przez nieznane współczynniki. Te współczynniki są wyznaczane z równania macierzowego — i to właśnie one stają się wartościami węzłowymi. Pole rozwiązania w dowolnym punkcie elementu otrzymuje się przez interpolację wartości węzłowych za pomocą funkcji kształtu.

Funkcje kształtu i własność interpolacyjna

Większość klasycznych funkcji kształtu w MES ma szczególną własność — w przypisanym sobie węźle przyjmują wartość 1, a we wszystkich pozostałych węzłach wartość 0. Jest to tzw. własność interpolacyjna (delta Kroneckera). Dzięki niej:

• współczynnik przy danej funkcji kształtu jest równy wartości pola w odpowiadającym węźle,
• wynik równania macierzowego można bezpośrednio interpretować jako wartość węzłową,
• nie jest potrzebna dodatkowa procedura wyznaczania wartości w węzłach.

Jeśli zastosowane funkcje kształtu nie mają tej własności (np. funkcje hierarchiczne, funkcje B-spline w analizie izogeometrycznej czy funkcje w metodach bezsiatkowych), wartość w węźle trzeba dodatkowo obliczyć — wstawiając współrzędne węzła do wyrażenia opisującego pole rozwiązania.

Wartości pierwotne a wielkości pochodne

W MES warto rozróżnić dwa rodzaje wyników:

• wartości pierwotne — bezpośrednio wyznaczane z układu równań (np. przemieszczenia, temperatury, ciśnienia),
• wielkości pochodne — obliczane na podstawie pochodnych pola pierwotnego (np. naprężenia, odkształcenia, strumień ciepła).

Pierwsze są wyznaczane dokładnie w węzłach. Drugie — ze względu na sposób ich obliczania — są zwykle wyznaczane w innych punktach elementu.

Punkty Gaussa a węzły

Naprężenia, odkształcenia i inne wielkości pochodne oblicza się najczęściej w punktach całkowania numerycznego, zwanych punktami Gaussa. Wynika to z kilku faktów:

• pochodne funkcji kształtu są w tych punktach najdokładniejsze (zjawisko superzbieżności),
• punkty Gaussa są naturalnym miejscem obliczeń przy całkowaniu macierzy sztywności,
• bezpośrednie obliczanie naprężeń w węzłach prowadzi często do większego błędu lokalnego.

Z tego powodu w wynikach MES często można spotkać wartości naprężeń podane zarówno w punktach całkowania (raw), jak i ekstrapolowane do węzłów.

Ekstrapolacja na węzły i uśrednianie

Na potrzeby wizualizacji i interpretacji wyniki z punktów Gaussa są przenoszone do węzłów. W węzłach współdzielonych przez wiele elementów wartości pochodne z różnych elementów zazwyczaj się różnią. Standardowo postępuje się wtedy następująco:

• ekstrapolacja wartości z punktów Gaussa do węzłów dla każdego elementu osobno,
• uśrednianie wartości w węzłach współdzielonych przez kilka elementów,
• wyświetlanie wyniku jako pola ciągłego (averaged) lub nieciągłego (unaveraged).

Wynik nieuśredniony pozwala dostrzec lokalne nieciągłości, które mogą świadczyć o niewystarczającym zagęszczeniu siatki lub o granicach między materiałami.

Praktyczne wskazówki dla interpretacji wartości węzłowych

Podczas interpretacji wartości węzłowych warto pamiętać o kilku zasadach:

• przemieszczenia i temperatury są wiarygodne wprost z węzłów,
• naprężenia i odkształcenia najlepiej oceniać przy włączonej wizualizacji nieuśrednionej,
• duże skoki naprężeń między sąsiednimi elementami sygnalizują potrzebę zagęszczenia siatki,
• na styku różnych materiałów uśrednianie naprężeń węzłowych może maskować rzeczywiste skoki,
• siły reakcji odczytuje się bezpośrednio w węzłach z zadanymi warunkami brzegowymi.