next up previous contents
: 実装 : 二次元アイソパラメトリック要素 : 実装   目次

ベクトル未知数の勾配 $ \frac{ \partial {}^{t} \{ u \} }{ \partial \{ x \} } $

自然座標 $ \{ \xi \} $ における ベクトル未知数の勾配 $ \frac{ \partial {}^{t} \{ u \} }{ \partial \{ x \} } $ は、


$\displaystyle \frac{ \partial {}^{t} \{ u \} }{ \partial \{ x \} }$ $\textstyle =$ $\displaystyle \frac{ \partial \{ u \} }{ \partial \{ x \} } ( \! ( \{ \xi \} ) \! )$  
  $\textstyle =$ $\displaystyle \sum_{In} \sum_{Id}
{}^{t} U_{(In Id)} [ N^x ] _{(In Id)}$ (4.23)

ここで、 $ [ N^x ] _{(In Id)}$ は 節点 $In$ 自由度 $Id$ での補間関数勾配テンソルであり、 要素タイプごとにそれぞれ 形状関数 $N_{(In)}$ やその他の付加情報などを用いて定められる。

もし、 ベクトル未知数の自由度が各軸方向の2成分により表現される場合、 $Id = 0 \sim 1$ であるから、 補間関数勾配テンソル $ [ N^x ] _{(In Id)}$ は、


$\displaystyle [ N^x ] _{(In Id)}$ $\textstyle =$ $\displaystyle [ N^x ] _{(In Id)} ( \! ( \{ \xi \} ) \! )$  
  $\textstyle =$ $\displaystyle \frac{ \partial N }{ \partial x } _{(In)} \{ e \} _{Id} \otimes \...
..._x
+ \frac{ \partial N }{ \partial y } _{(In)} \{ e \} _{Id} \otimes \{ e \} _y$ (4.24)





Hiroshi KAWAI 平成15年8月11日