next up previous contents
: 実装 : 二次元シンプレックス要素 : 実装   目次

位置 $ \{ x \} $

二次元シンプレックス要素では、 位置(全体座標) $ \{ x \} $ は二次元ベクトルである。


$\displaystyle \{ x \}
=
{
\left \langle \begin{array}{cc}
x & y
\end{array} \right \rangle
} ^ { T }$     (1.10)

指定の面積座標 $L0, L1, L2$ に対応する 位置(全体座標) $ \{ x \} $ は、 以下のようになる。


$\displaystyle \{ x \}
=
\{ x \} ( \! ( L0, L1, L2 ) \! )
=
L0 \{ x0 \} + L1 \{ x1 \} + L2 \{ x2 \}$     (1.11)

また、 指定の位置(全体座標) $ \{ x \} $ に対応する面積座標 $L0, L1, L2$ は、 以下のように求まる。


$\displaystyle L0$ $\textstyle =$ $\displaystyle L0 ( \! ( \{ x \} ) \! ) = a0 + b0 x + c0 y$  
$\displaystyle L1$ $\textstyle =$ $\displaystyle L1 ( \! ( \{ x \} ) \! ) = a1 + b1 x + c1 y$  
$\displaystyle L2$ $\textstyle =$ $\displaystyle L2 ( \! ( \{ x \} ) \! ) = a2 + b2 x + c2 y$ (1.12)

ここで、 $a0$, $b1$, $c2$ などは以下のように定義される。


$\displaystyle a0$ $\textstyle =$ $\displaystyle \frac{1}{2 S} (x1 y2 - x2 y1)$  
$\displaystyle a1$ $\textstyle =$ $\displaystyle \frac{1}{2 S} (x2 y0 - x0 y2)$  
$\displaystyle a2$ $\textstyle =$ $\displaystyle \frac{1}{2 S} (x0 y1 - x1 y0)$ (1.13)


$\displaystyle b0$ $\textstyle =$ $\displaystyle \frac{1}{2 S} (y1 - y2)$  
$\displaystyle b1$ $\textstyle =$ $\displaystyle \frac{1}{2 S} (y2 - y0)$  
$\displaystyle b2$ $\textstyle =$ $\displaystyle \frac{1}{2 S} (y0 - y1)$ (1.14)


$\displaystyle c0$ $\textstyle =$ $\displaystyle \frac{1}{2 S} (x2 - x1)$  
$\displaystyle c1$ $\textstyle =$ $\displaystyle \frac{1}{2 S} (x0 - x2)$  
$\displaystyle c2$ $\textstyle =$ $\displaystyle \frac{1}{2 S} (x1 - x0)$ (1.15)

さらに、 面積座標の勾配、すなわち、位置(全体座標) $ \{ x \} $ による微分は、


$\displaystyle \frac{ \partial L0 }{ \partial \{ x \} }$ $\textstyle =$ $\displaystyle {
\left \langle \begin{array}{cc}
b0 & c0
\end{array} \right \rangle
} ^ { T }$  
$\displaystyle \frac{ \partial L1 }{ \partial \{ x \} }$ $\textstyle =$ $\displaystyle {
\left \langle \begin{array}{cc}
b1 & c1
\end{array} \right \rangle
} ^ { T }$  
$\displaystyle \frac{ \partial L2 }{ \partial \{ x \} }$ $\textstyle =$ $\displaystyle {
\left \langle \begin{array}{cc}
b2 & c2
\end{array} \right \rangle
} ^ { T }$ (1.16)




next up previous contents
: 実装 : 二次元シンプレックス要素 : 実装   目次
Hiroshi KAWAI 平成15年8月11日