next up previous contents
: 実装 : 三角形 : 自然座標   目次

領域積分

要素 $Ie$ における 関数 $ {}^{t} f$ について、


$\displaystyle {}^{t} f = f ( \! ( \{ \xi \} ) \! )$     (5.32)

その領域積分は、


    $\displaystyle \int_{S^{(Ie)}}
{}^{t} f
\mathrm{d} S$  
  $\textstyle =$ $\displaystyle \int_0^1 \int_0^{1-\eta}
{}^{t} f {}^{t} J^S
\mathrm{d} \xi \mathrm{d} \eta$  
  $\textstyle =$ $\displaystyle \sum_{(Ip)}
w_{(Ip)}
f ( \! ( \{ \xi \} _{(Ip)} ) \! )
J^S ( \! ( \{ \xi \} _{(Ip)} ) \! )$ (5.33)

ここで、 面積積分の変換子 $ {}^{t} J^S$ は、 Jacobianのdeterminantの $1/2$ 、すなわち、三角形の面積となり、


$\displaystyle {}^{t} J^S$ $\textstyle =$ $\displaystyle J^S ( \! ( \{ \xi \} ) \! )$  
  $\textstyle =$ $\displaystyle 1/2 \mathrm{det} \; \left[ \frac{ \partial {}^{t} \{ x \} }{ \partial \{ \xi \} } \right] a$ (5.34)

$ \{ \xi \} _{(Ip)}$ は、三角形の積分点 $Ip$ の自然座標、 $w_{(Ip)}$ は、その重み係数である。

各方向の数値積分の座標と重みについては、 「数値積分」「ガウス積分」「2次元三角形」を参照。



Hiroshi KAWAI 平成15年4月19日