OGS
NumLib::ShapePrism6 Class Reference

Detailed Description

Shape function for a 6-nodes prism element in natural coordinates

Definition at line 15 of file ShapePrism6.h.

#include <ShapePrism6.h>

Public Types

using MeshElement = MeshLib::Prism

Static Public Member Functions

template<class T_X, class T_N>
static void computeShapeFunction (const T_X &x, T_N &N)
template<class T_X, class T_N>
static void computeGradShapeFunction (const T_X &x, T_N &dN)

Static Public Attributes

static constexpr std::array reference_element_centre
static const unsigned DIM = MeshElement::dimension
static const unsigned NPOINTS = MeshElement::n_all_nodes
static constexpr int ORDER = 1

Member Typedef Documentation

◆ MeshElement

using NumLib::ShapePrism6::MeshElement = MeshLib::Prism

Definition at line 39 of file ShapePrism6.h.

Member Function Documentation

◆ computeGradShapeFunction()

template<class T_X, class T_N>
void NumLib::ShapePrism6::computeGradShapeFunction ( const T_X & x,
T_N & dN )
static

Evaluate derivatives of the shape function at the given point

Parameters
[in]xnatural coordinates (r,s,t)
[out]dNa matrix of the derivatives

Definition at line 21 of file ShapePrism6-impl.h.

22{
23 const double L1 = x[0];
24 const double L2 = x[1];
25 const double t = x[2];
26 // dN/dL1
27 dN[0] = -0.5 * (1.0 - t);
28 dN[1] = -dN[0];
29 dN[2] = 0.0;
30 dN[3] = -0.5 * (1.0 + t);
31 dN[4] = -dN[3];
32 dN[5] = 0.0;
33 // dN/dL2
34 dN[6] = dN[0];
35 dN[7] = 0.0;
36 dN[8] = -dN[0];
37 dN[9] = dN[3];
38 dN[10] = 0.0;
39 dN[11] = -dN[3];
40 // dN/dt
41 dN[12] = -0.5 * (1.0 - L1 - L2);
42 dN[13] = -0.5 * L1;
43 dN[14] = -0.5 * L2;
44 dN[15] = -dN[12];
45 dN[16] = -dN[13];
46 dN[17] = -dN[14];
47}

◆ computeShapeFunction()

template<class T_X, class T_N>
void NumLib::ShapePrism6::computeShapeFunction ( const T_X & x,
T_N & N )
static

Evaluate the shape function at the given point

Parameters
[in]xnatural coordinates (r,s,t)
[out]Na vector of calculated shape functions

Definition at line 7 of file ShapePrism6-impl.h.

8{
9 const double L1 = x[0];
10 const double L2 = x[1];
11 const double t = x[2];
12 N[0] = 0.5 * (1.0 - L1 - L2) * (1.0 - t);
13 N[1] = 0.5 * L1 * (1.0 - t);
14 N[2] = 0.5 * L2 * (1.0 - t);
15 N[3] = 0.5 * (1.0 - L1 - L2) * (1.0 + t);
16 N[4] = 0.5 * L1 * (1.0 + t);
17 N[5] = 0.5 * L2 * (1.0 + t);
18}

References NumLib::N.

Member Data Documentation

◆ DIM

const unsigned NumLib::ShapePrism6::DIM = MeshElement::dimension
static

Definition at line 40 of file ShapePrism6.h.

◆ NPOINTS

const unsigned NumLib::ShapePrism6::NPOINTS = MeshElement::n_all_nodes
static

Definition at line 41 of file ShapePrism6.h.

◆ ORDER

int NumLib::ShapePrism6::ORDER = 1
staticconstexpr

Definition at line 42 of file ShapePrism6.h.

◆ reference_element_centre

std::array NumLib::ShapePrism6::reference_element_centre
staticconstexpr
Initial value:
= {1. / 3., 1. / 3.,
0.0}

Definition at line 36 of file ShapePrism6.h.

36 {1. / 3., 1. / 3.,
37 0.0};
The documentation for this class was generated from the following files: