Funkcje kształtu MES



Czym są funkcje kształtu?

Funkcje kształtu danego elementu opisują rozkład parametrów w elemencie skończonym.

Co to oznacza w praktyce?

Jeżeli policzymy korzystając z MES np. wartość przemieszczeń w węzłąch trójkątnego elementu skończonego to nie znaczy, że znamy przemieszczenie dowolnego punktu tego elementu. Musimy więc dokonać interpolacji naszych wyników.

Funkcje kształtu muszą spełniać następujące warunki:

1. W i-tym węźle (czyli tym względem którego w danym momencie liczymy) musi być równa jeden

2. W pozostałych węzłąch musi być równa zero

3. Suma funkcji kształtu ZAWSZE wynosi jeden


Obliczenie funkcji kształtu dla elementu trójkątnego

Wprost z powyższych wynika w jaki sposób możemy wyznaczyć funkcje kształtu dowolnego elementu

Zapisujemy funkcje dla poszczególnych węzłów i rozwiązujemy powstałe układy równań

wezeł 1

wezeł 2

𝑁𝑗(𝑥,𝑦)=𝐴𝑥+𝐵𝑦+𝐶𝑁2(0,2)=𝐴0+𝐵2+𝐶=0𝑁2(1,0)=𝐴1+𝐵0+𝐶=1𝑁2(2,2)=𝐴2+𝐵2+𝐶=0𝐴=0𝐵=12𝐶=1𝑁2(𝑥,𝑦)=12𝑦+1

wezeł 3

𝑁𝑘(𝑥,𝑦)=𝐴𝑥+𝐵𝑦+𝐶𝑁3(0,2)=𝐴0+𝐵2+𝐶=0𝑁3(1,0)=𝐴1+𝐵0+𝐶=0𝑁3(2,2)=𝐴2+𝐵2+𝐶=1𝐴=12𝐵=14𝐶=12𝑁3(𝑥,𝑦)=12𝑥+14𝑦12

W praktyce jeżeli obliczamy funkcje kształu korzystając z komputera nie z kartki/długopisu/kalkluatora dużo więcej sensu ma równoważny zapis:

𝑥1𝑦1𝑥2𝑦2𝑥3𝑦3=021022𝑎1𝑎2𝑎3𝑏1𝑏2𝑏3𝑐1𝑐2𝑐3=𝑥1𝑦11𝑥2𝑦21𝑥3𝑦3111201214121412112𝑁1(𝑥,𝑦)=𝑎1𝑥+𝑏1𝑦+𝑐1𝑦4+(12𝑥2)𝑁2(𝑥,𝑦)=𝑎2𝑥+𝑏2𝑦+𝑐2𝑦2+1𝑁3(𝑥,𝑦)=𝑎3𝑥+𝑏3𝑦+𝑐3𝑦4+(𝑥212)

Jak widać wyniki są jednoznaczne i w obu przypadkach suma funkcji kształtu to jeden:

𝑁1(𝑥,𝑦)+𝑁2(𝑥,𝑦)+𝑁3(𝑥,𝑦)=1𝑦4+12𝑥2𝑦2+1+𝑦4+𝑥212=1

Jednoznaczność wyników wynika oczywiście z tego że drugie podejście nie jest niczym innym jak przedstawieniem macierzowym równań z podejścia pierwszego:

𝑁1=𝑎1𝑥+𝑏1𝑦+𝑐1 𝑥1𝑦11𝑥2𝑦21𝑥3𝑦31𝑎1𝑏1𝑐1=100=𝑎1𝑏1𝑐1𝑥1𝑦11𝑥2𝑦21𝑥3𝑦311100 𝑁2=𝑎2𝑥+𝑏2𝑦+𝑐2 𝑥1𝑦11𝑥2𝑦21𝑥3𝑦31𝑎2𝑏2𝑐2=010=𝑎2𝑏2𝑐2𝑥1𝑦11𝑥2𝑦21𝑥3𝑦311010 𝑁3=𝑎3𝑥+𝑏3𝑦+𝑐3 𝑥1𝑦11𝑥2𝑦21𝑥3𝑦31𝑎3𝑏3𝑐3=001=𝑎3𝑏3𝑐3𝑥1𝑦11𝑥2𝑦21𝑥3𝑦311001 𝑎1𝑎2𝑎3𝑏1𝑏2𝑏3𝑐1𝑐2𝑐3=𝑥1𝑦11𝑥2𝑦21𝑥3𝑦311
Łukasz Cichowicz
GRZEGORZ
MAZUR
Korepetytor
200 PLN
60 MIN
Zobacz mój profil na
Edupanda Logo