流体力学から数値計算まで

広告


・4面体要素の体積

4面体要素の体積を示します。

以降の式で V は正の要素体積を表します。節点順序によって体積行列式が負になる場合は、絶対値を用いるか、正の向きになるよう節点順序を統一します。

\[ \begin{aligned} V&=\frac{1}{6}\det \begin{pmatrix} 1&x_1&y_1&z_1\\ 1&x_2&y_2&z_2\\ 1&x_3&y_3&z_3\\ 1&x_4&y_4&z_4 \end{pmatrix}\\[6pt] &=\frac{1}{6}\Big( x_2y_3z_4-x_2z_3y_4-y_2x_3z_4+y_2z_3x_4+z_2x_3y_4-z_2y_3x_4\\[6pt] &\qquad-x_1y_3z_4+x_1z_3y_4+x_1y_2z_4-x_1y_2z_3-x_1z_2y_4+x_1z_2y_3\\[6pt] &\qquad+y_1x_3z_4-y_1z_3x_4-y_1x_2z_4+y_1x_2z_3+y_1z_2x_4-y_1z_2x_3\\[6pt] &\qquad-z_1x_3y_4+z_1y_3x_4+z_1x_2y_4-z_1x_2y_3-z_1y_2x_4+z_1y_2x_3 \Big) \end{aligned} \]

ここで、

\[ \begin{aligned} &a_{11}=1,\ a_{12}=x_1,\ a_{13}=y_1,\ a_{14}=z_1\\[6pt] &a_{21}=1,\ a_{22}=x_2,\ a_{23}=y_2,\ a_{24}=z_2\\[6pt] &a_{31}=1,\ a_{32}=x_3,\ a_{33}=y_3,\ a_{34}=z_3\\[6pt] &a_{41}=1,\ a_{42}=x_4,\ a_{43}=y_4,\ a_{44}=z_4 \end{aligned} \]

同様に

\[ \begin{aligned} V_1&=\frac{1}{6}\det \begin{pmatrix} 1&x&y&z\\ 1&x_2&y_2&z_2\\ 1&x_3&y_3&z_3\\ 1&x_4&y_4&z_4 \end{pmatrix}\\[6pt] &=\frac{1}{6}\Big( x_2y_3z_4-x_2z_3y_4-y_2x_3z_4+y_2z_3x_4+z_2x_3y_4-z_2y_3x_4\\[6pt] &\qquad-xy_3z_4+xz_3y_4+xy_2z_4-xy_2z_3-xz_2y_4+xz_2y_3\\[6pt] &\qquad+yx_3z_4-yz_3x_4-yx_2z_4+yx_2z_3+yz_2x_4-yz_2x_3\\[6pt] &\qquad-zx_3y_4+zy_3x_4+zx_2y_4-zx_2y_3-zy_2x_4+zy_2x_3 \Big) \end{aligned} \]

\[ \begin{aligned} V_2&=\frac{1}{6}\det \begin{pmatrix} 1&x_1&y_1&z_1\\ 1&x&y&z\\ 1&x_3&y_3&z_3\\ 1&x_4&y_4&z_4 \end{pmatrix}\\[6pt] &=\frac{1}{6}\Big( xy_3z_4-xz_3y_4-yx_3z_4+yz_3x_4+zx_3y_4-zy_3x_4\\[6pt] &\qquad-x_1y_3z_4+x_1z_3y_4+x_1yz_4-x_1yz_3-x_1zy_4+x_1zy_3\\[6pt] &\qquad+y_1x_3z_4-y_1z_3x_4-y_1xz_4+y_1xz_3+y_1zx_4-y_1zx_3\\[6pt] &\qquad-z_1x_3y_4+z_1y_3x_4+z_1xy_4-z_1xy_3-z_1yx_4+z_1yx_3 \Big) \end{aligned} \]

\[ \begin{aligned} V_3&=\frac{1}{6}\det \begin{pmatrix} 1&x_1&y_1&z_1\\ 1&x_2&y_2&z_2\\ 1&x&y&z\\ 1&x_4&y_4&z_4 \end{pmatrix}\\[6pt] &=\frac{1}{6}\Big( x_2yz_4-x_2zy_4-y_2xz_4+y_2zx_4+z_2xy_4-z_2yx_4\\[6pt] &\qquad-x_1yz_4+x_1zy_4+x_1y_2z_4-x_1y_2z-x_1z_2y_4+x_1z_2y\\[6pt] &\qquad+y_1xz_4-y_1zx_4-y_1x_2z_4+y_1x_2z+y_1z_2x_4-y_1z_2x\\[6pt] &\qquad-z_1xy_4+z_1yx_4+z_1x_2y_4-z_1x_2y-z_1y_2x_4+z_1y_2x \Big) \end{aligned} \]

\[ \begin{aligned} V_4&=\frac{1}{6}\det \begin{pmatrix} 1&x_1&y_1&z_1\\ 1&x_2&y_2&z_2\\ 1&x_3&y_3&z_3\\ 1&x&y&z \end{pmatrix}\\[6pt] &=\frac{1}{6}\Big( x_2y_3z-x_2z_3y-y_2x_3z+y_2z_3x+z_2x_3y-z_2y_3x\\[6pt] &\qquad-x_1y_3z+x_1z_3y+x_1y_2z-x_1y_2z_3-x_1z_2y+x_1z_2y_3\\[6pt] &\qquad+y_1x_3z-y_1z_3x-y_1x_2z+y_1x_2z_3+y_1z_2x-y_1z_2x_3\\[6pt] &\qquad-z_1x_3y+z_1y_3x+z_1x_2y-z_1x_2y_3-z_1y_2x+z_1y_2x_3 \Big) \end{aligned} \]

prev | up | next

カウンタ

(2011.3.15~)