One Dimensional Finite Element Method
Parabolic Element-2

もう1つ形状関数の条件がありました。それは、下式です。この条件は、要素内の全てのxについて満足しなくてはなりません。

N1(x) + N2(x) + N3(x) = 1

これらの条件を満足する形状関数を探し出すのは大変ですが、1次元の2次要素程度だと、Lagrangeの補間式で作ることが出来ます。下図に形状関数を示しますので、貴方はそれらが条件を満たしているかどうかをチェックして下さい。

上の形状関数を一筆で書くと、下の様になる。式中の変な記号は、ギリシャのπの 大文字です。掛け算をするという意味です。参考までに、Σは足し算です。

そして、u(x)は下の様に書けます。

以外と簡単に形状関数できてしまいましたが、有限要素法のプログラムを完成するまでは、まだまだ沢山の作業が残っています。

BACK NEXT
Menu View Helm wrm Lin Element Rmrk Vari Para Non-L Wire