ProcessLib::LinearBMatrix::detail Namespace Reference

Functions
template<int NPOINTS, typename DNDX_Type , typename BMatrixType >
void	fillBMatrix2DCartesianPart (DNDX_Type const &dNdx, BMatrixType &B)
template<int DisplacementDim, int NPOINTS, typename BBarMatrixType , typename BMatrixType >
void	applyBbar (BBarMatrixType const &B_bar, BMatrixType &B, const bool is_axially_symmetric)

Function Documentation

applyBbar()

template<int DisplacementDim, int NPOINTS, typename BBarMatrixType , typename BMatrixType >

void ProcessLib::LinearBMatrix::detail::applyBbar	(	BBarMatrixType const &	B_bar,
		BMatrixType &	B,
		const bool	is_axially_symmetric )

Definition at line 141 of file LinearBMatrix.h.

{
    if constexpr (DisplacementDim == 3)
    {
        for (int i = 0; i < NPOINTS; ++i)
        {
            auto const B_bar_i = B_bar.col(i);
 
            // The following loop is based on the following facts:
            // B1 (dN/dx) is the 1st element of the 1st column of B,
            // B2 (dN/dy) is the 2nd element of the 2nd column of B,
            // B3 (dN/dz) is the 3rd element of the 3rd column of B.
            for (int k = 0; k < 3; k++)
            {
                // k-th column of B matrix at node i
                auto B_i_k = B.col(k * NPOINTS + i);
 
                B_i_k.template segment<3>(0) -=
                    Eigen::Vector3d::Constant((B_i_k[k] - B_bar_i[k]) / 3.0);
            }
        }
 
        return;
    }
 
    // 2D or axisymmetry
    for (int i = 0; i < NPOINTS; ++i)
    {
        auto B_i_0 = B.col(i);
        auto B_i_1 = B.col(NPOINTS + i);
 
        auto const B_bar_i = B_bar.col(i);
 
        if (is_axially_symmetric)
        {
            double const b0_dil_pertubation =
                (B_i_0[0] - B_bar_i[0] + B_i_0[2] - B_bar_i[2]);
            B_i_0.template segment<3>(0) -=
                Eigen::Vector3d::Constant((b0_dil_pertubation) / 3.);
            B_i_1.template segment<3>(0) -=
                Eigen::Vector3d::Constant((B_i_1[1] - B_bar_i[1]) / 3.);
            continue;
        }
 
        // Plane strain
        B_i_0.template segment<2>(0) -=
            Eigen::Vector2d::Constant((B_i_0[0] - B_bar_i[0]) / 2.);
        B_i_1.template segment<2>(0) -=
            Eigen::Vector2d::Constant((B_i_1[1] - B_bar_i[1]) / 2.);
    }
}

fillBMatrix2DCartesianPart()

template<int NPOINTS, typename DNDX_Type , typename BMatrixType >

void ProcessLib::LinearBMatrix::detail::fillBMatrix2DCartesianPart	(	DNDX_Type const &	dNdx,
		BMatrixType &	B )

Definition at line 28 of file LinearBMatrix.h.

{
    for (int i = 0; i < NPOINTS; ++i)
    {
        B(1, NPOINTS + i) = dNdx(1, i);
        B(3, i) = dNdx(1, i) / std::sqrt(2);
        B(3, NPOINTS + i) = dNdx(0, i) / std::sqrt(2);
        B(0, i) = dNdx(0, i);
    }
}