19namespace LinearBMatrix
23template <
int NPOINTS,
typename DNDX_Type,
typename BMatrixType>
26 for (
int i = 0; i < NPOINTS; ++i)
28 B(1, NPOINTS + i) = dNdx(1, i);
29 B(3, i) = dNdx(1, i) / std::sqrt(2);
30 B(3, NPOINTS + i) = dNdx(0, i) / std::sqrt(2);
37template <
int DisplacementDim,
45 const bool is_axially_symmetric)
47 static_assert(0 < DisplacementDim && DisplacementDim <= 3,
48 "LinearBMatrix::computeBMatrix: DisplacementDim must be in "
51 BMatrixType B = BMatrixType::Zero(
53 NPOINTS * DisplacementDim);
55 switch (DisplacementDim)
58 for (
int i = 0; i < NPOINTS; ++i)
60 B(2, 2 * NPOINTS + i) = dNdx(2, i);
61 B(4, NPOINTS + i) = dNdx(2, i) / std::sqrt(2);
62 B(4, 2 * NPOINTS + i) = dNdx(1, i) / std::sqrt(2);
63 B(5, i) = dNdx(2, i) / std::sqrt(2);
64 B(5, 2 * NPOINTS + i) = dNdx(0, i) / std::sqrt(2);
66 detail::fillBMatrix2DCartesianPart<NPOINTS>(dNdx, B);
69 detail::fillBMatrix2DCartesianPart<NPOINTS>(dNdx, B);
70 if (is_axially_symmetric)
72 for (
int i = 0; i < NPOINTS; ++i)
74 B(2, i) = N[i] / radius;
constexpr int kelvin_vector_dimensions(int const displacement_dim)
Kelvin vector dimensions for given displacement dimension.
void fillBMatrix2DCartesianPart(DNDX_Type const &dNdx, BMatrixType &B)
BMatrixType computeBMatrix(DNDX_Type const &dNdx, N_Type const &N, const double radius, const bool is_axially_symmetric)
Fills a B-matrix based on given shape function dN/dx values.