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

広告


・内挿関数N3の微分

4面体1次要素 の内挿関数N3の微分値を示します。



x微分

\[ \begin{aligned} \frac{\partial N_3}{\partial x} &= \frac{\partial L_3}{\partial x} = \frac{\partial(V_3/V)}{\partial x} = \frac{1}{V}\frac{\partial V_3}{\partial x} \\[6pt] &= \frac{1}{V}\frac{\partial}{\partial x}\frac{1}{6} \Big( x_2y z_4-x_2z y_4-y_2x z_4+y_2z x_4+z_2x y_4-z_2y x_4\\[6pt] &\qquad -x_1y z_4+x_1z y_4+x_1y_2z_4-x_1y_2z-x_1z_2y_4+x_1z_2y\\[6pt] &\qquad +y_1x z_4-y_1z x_4-y_1x_2z_4+y_1x_2z+y_1z_2x_4-y_1z_2x\\[6pt] &\qquad -z_1x y_4+z_1y x_4+z_1x_2y_4-z_1x_2y-z_1y_2x_4+z_1y_2x \Big) \\[6pt] &= \frac{1}{6V}\Big( -y_2z_4+z_2y_4+y_1z_4-y_1z_2-z_1y_4+z_1y_2 \Big) \end{aligned} \]


y微分

\[ \begin{aligned} \frac{\partial N_3}{\partial y} &= \frac{\partial L_3}{\partial y} = \frac{\partial(V_3/V)}{\partial y} = \frac{1}{V}\frac{\partial V_3}{\partial y} \\[6pt] &= \frac{1}{V}\frac{\partial}{\partial y}\frac{1}{6} \Big( x_2y z_4-x_2z y_4-y_2x z_4+y_2z x_4+z_2x y_4-z_2y x_4\\[6pt] &\qquad -x_1y z_4+x_1z y_4+x_1y_2z_4-x_1y_2z-x_1z_2y_4+x_1z_2y\\[6pt] &\qquad +y_1x z_4-y_1z x_4-y_1x_2z_4+y_1x_2z+y_1z_2x_4-y_1z_2x\\[6pt] &\qquad -z_1x y_4+z_1y x_4+z_1x_2y_4-z_1x_2y-z_1y_2x_4+z_1y_2x \Big) \\[6pt] &= \frac{1}{6V}\Big( x_2z_4-z_2x_4-x_1z_4+x_1z_2+z_1x_4-z_1x_2 \Big) \end{aligned} \]


z微分

\[ \begin{aligned} \frac{\partial N_3}{\partial z} &= \frac{\partial L_3}{\partial z} = \frac{\partial(V_3/V)}{\partial z} = \frac{1}{V}\frac{\partial V_3}{\partial z} \\[6pt] &= \frac{1}{V}\frac{\partial}{\partial z}\frac{1}{6} \Big( x_2y z_4-x_2z y_4-y_2x z_4+y_2z x_4+z_2x y_4-z_2y x_4\\[6pt] &\qquad -x_1y z_4+x_1z y_4+x_1y_2z_4-x_1y_2z-x_1z_2y_4+x_1z_2y\\[6pt] &\qquad +y_1x z_4-y_1z x_4-y_1x_2z_4+y_1x_2z+y_1z_2x_4-y_1z_2x\\[6pt] &\qquad -z_1x y_4+z_1y x_4+z_1x_2y_4-z_1x_2y-z_1y_2x_4+z_1y_2x \Big) \\[6pt] &= \frac{1}{6V}\Big( -x_2y_4+y_2x_4+x_1y_4-x_1y_2-y_1x_4+y_1x_2 \Big) \end{aligned} \]

prev | up | next

カウンタ

(2011.3.15~)