# Linear Brick Shape Functions

Shape functions define the geometry of an element in its computational (intrinsic) domain. As was seen in Finite Element Formulation, physical coordinates are transformed into simpler computational intrinsic coordinates so that integration of values is numerically more efficient.

Where, $r\equiv \xi$ , $s\equiv \eta$ , and $t\equiv \zeta$ .

The shape functions of an 8 node brick element, shown in Figure 1, are given by:

${\Phi }_{1}=\frac{1}{8}\left(1-\xi \right)\left(1-\eta \right)\left(1-\zeta \right)$
${\Phi }_{2}=\frac{1}{8}\left(1-\xi \right)\left(1-\eta \right)\left(1+\zeta \right)$
${\Phi }_{3}=\frac{1}{8}\left(1+\xi \right)\left(1-\eta \right)\left(1+\zeta \right)$
${\Phi }_{4}=\frac{1}{8}\left(1+\xi \right)\left(1-\eta \right)\left(1-\zeta \right)$
${\Phi }_{5}=\frac{1}{8}\left(1-\xi \right)\left(1+\eta \right)\left(1-\zeta \right)$
${\Phi }_{6}=\frac{1}{8}\left(1-\xi \right)\left(1+\eta \right)\left(1+\zeta \right)$
${\Phi }_{7}=\frac{1}{8}\left(1+\xi \right)\left(1+\eta \right)\left(1+\zeta \right)$
${\Phi }_{8}=\frac{1}{8}\left(1+\xi \right)\left(1+\eta \right)\left(1-\zeta \right)$

The element velocity field is related by:

${v}_{i}=\sum _{I=1}^{8}{\Phi }_{I}.{v}_{iI}$

Where, ${v}_{iI}$ are the nodal velocities.