26 GMatrixType& g_matrix,
27 const bool is_axially_symmetric,
31 static_assert(0 < DisplacementDim && DisplacementDim <= 3,
32 "LinearGMatrix::computeGMatrix: DisplacementDim must be in "
37 switch (DisplacementDim)
44 for (
int d = 0; d < DisplacementDim; ++d)
46 for (
int i = 0; i < NPOINTS; ++i)
48 g_matrix(d + 0 * DisplacementDim, i + 0 * NPOINTS) =
50 g_matrix(d + 1 * DisplacementDim, i + 1 * NPOINTS) =
52 g_matrix(d + 2 * DisplacementDim, i + 2 * NPOINTS) =
62 for (
int d = 0; d < DisplacementDim; ++d)
64 for (
int i = 0; i < NPOINTS; ++i)
66 g_matrix(d, i) = dNdx(d, i);
67 g_matrix(d + DisplacementDim, i + NPOINTS) = dNdx(d, i);
70 if (is_axially_symmetric)
72 for (
int i = 0; i < NPOINTS; ++i)
74 g_matrix(4, i) = N[i] / radius;