Enumerations
enum class	BarDetFType { ELEMENT_CENTER_VALUE , ELEMENT_AVERAGE , NONE }

Functions
BarDetFType	convertStringToDetFBarType (std::string_view const bar_det_f_type_name)
template<int DisplacementDim, int NPOINTS, typename VectorTypeForFbar, typename GradientVectorType, typename DNDX_Type>
VectorTypeForFbar	computeQVector (DNDX_Type const &dNdx, GradientVectorType const &F, bool const)
template<int DisplacementDim, typename GradientVectorType, typename VectorTypeForFbar, typename NodalVectorType, typename ShapeFunction, typename ShapeMatricesType, typename IpData>
std::tuple< double, VectorTypeForFbar >	computeFBarInitialVariablesAverage (std::vector< IpData, Eigen::aligned_allocator< IpData > > const &ip_data, bool const compute_detF0_only, Eigen::VectorXd const &u, NumLib::GenericIntegrationMethod const &integration_method, MeshLib::Element const &element, bool const is_axially_symmetric)
template<int DisplacementDim, typename GradientVectorType, typename GradientMatrixType, typename VectorTypeForFbar, typename ShapeFunction, typename ShapeMatricesType>
std::tuple< double, VectorTypeForFbar >	computeFBarInitialVariablesElementCenter (bool const compute_detF0_only, Eigen::VectorXd const &u, MeshLib::Element const &element, bool const is_axially_symmetric)
template<int DisplacementDim>
MathLib::KelvinVector::KelvinVectorType< DisplacementDim >	identityForF (bool const is_axially_symmetric)
template<int DisplacementDim>
MathLib::KelvinVector::KelvinVectorType< DisplacementDim >	compute2EPlusI (double const alpha_p2, MathLib::KelvinVector::KelvinVectorType< DisplacementDim > const &eps_bar, bool const is_axially_symmetric)

Enumeration Type Documentation

◆ BarDetFType

enum class ProcessLib::NonLinearFbar::BarDetFType

strong

Enumerator
ELEMENT_CENTER_VALUE
ELEMENT_AVERAGE
NONE

Definition at line 23 of file NonLinearFbar.h.

{
    ELEMENT_CENTER_VALUE,
    ELEMENT_AVERAGE,
    NONE
};

Function Documentation

◆ compute2EPlusI()

template<int DisplacementDim>

MathLib::KelvinVector::KelvinVectorType< DisplacementDim > ProcessLib::NonLinearFbar::compute2EPlusI	(	double const	alpha_p2,
		MathLib::KelvinVector::KelvinVectorType< DisplacementDim > const &	eps_bar,
		bool const	is_axially_symmetric )

Definition at line 212 of file NonLinearFbar.h.

{
    return (2.0 * eps_bar +
            identityForF<DisplacementDim>(is_axially_symmetric)) /
           alpha_p2;
}

References identityForF().

Referenced by ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::assembleWithJacobian().

◆ computeFBarInitialVariablesAverage()

template<int DisplacementDim, typename GradientVectorType, typename VectorTypeForFbar, typename NodalVectorType, typename ShapeFunction, typename ShapeMatricesType, typename IpData>

std::tuple< double, VectorTypeForFbar > ProcessLib::NonLinearFbar::computeFBarInitialVariablesAverage	(	std::vector< IpData, Eigen::aligned_allocator< IpData > > const &	ip_data,
		bool const	compute_detF0_only,
		Eigen::VectorXd const &	u,
		NumLib::GenericIntegrationMethod const &	integration_method,
		MeshLib::Element const &	element,
		bool const	is_axially_symmetric )

Definition at line 69 of file NonLinearFbar.h.

{
    unsigned const n_integration_points =
        integration_method.getNumberOfPoints();
    // Compute the element volume
    double volume = 0.0;
    for (unsigned ip = 0; ip < n_integration_points; ip++)
    {
        auto const& w = ip_data[ip].integration_weight;
        volume += w;
    }
 
    VectorTypeForFbar averaged_grad_N =
        VectorTypeForFbar::Zero(DisplacementDim * ShapeFunction::NPOINTS);
    NodalVectorType averaged_N_div_r;
    if (is_axially_symmetric)
    {
        averaged_N_div_r = NodalVectorType::Zero(ShapeFunction::NPOINTS);
    }
 
    GradientVectorType F0 =
        MathLib::VectorizedTensor::identity<DisplacementDim>();
 
    for (unsigned i = 0; i < ShapeFunction::NPOINTS; i++)
    {
        Eigen::Vector3d const bar_gradN =
            NumLib::averageGradShapeFunction<DisplacementDim, ShapeFunction,
                                             ShapeMatricesType, IpData>(
                i, element, integration_method, ip_data, is_axially_symmetric) /
            volume;
        averaged_grad_N.template segment<DisplacementDim>(i * DisplacementDim) =
            bar_gradN.template segment<DisplacementDim>(0);
 
        for (int k = 0; k < DisplacementDim; k++)
        {
            F0.template segment<DisplacementDim>(k * DisplacementDim) +=
                u[k * ShapeFunction::NPOINTS + i] *
                bar_gradN.template segment<DisplacementDim>(0);
        }
        if (is_axially_symmetric)
        {
            averaged_N_div_r[i] = bar_gradN[2];
            F0[4] += bar_gradN[2] * u[i];
        }
    }
 
    if (compute_detF0_only)
    {
        return {MathLib::VectorizedTensor::determinant(F0), {}};
    }
 
    VectorTypeForFbar F0InvN =
        VectorTypeForFbar::Zero(DisplacementDim * ShapeFunction::NPOINTS);
 
    Eigen::MatrixXd const F_matrix =
        MathLib::VectorizedTensor::toTensor<DisplacementDim>(F0)
            .template topLeftCorner<DisplacementDim, DisplacementDim>();
    auto const Finv = F_matrix.inverse();
 
    for (unsigned i = 0; i < ShapeFunction::NPOINTS; ++i)
    {
        auto const dNidx = averaged_grad_N.template segment<DisplacementDim>(
            i * DisplacementDim);
        for (int k = 0; k < DisplacementDim; k++)
        {
            F0InvN[k * ShapeFunction::NPOINTS + i] = Finv.col(k).dot(dNidx);
        }
        // TODO: if(is_axially_symmetric)
    }
 
    return {MathLib::VectorizedTensor::determinant(F0), F0InvN};
}

References NumLib::averageGradShapeFunction(), MathLib::VectorizedTensor::determinant(), NumLib::GenericIntegrationMethod::getNumberOfPoints(), MathLib::VectorizedTensor::identity(), and MathLib::VectorizedTensor::toTensor().

Referenced by ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::computeFBarVariables().

◆ computeFBarInitialVariablesElementCenter()

template<int DisplacementDim, typename GradientVectorType, typename GradientMatrixType, typename VectorTypeForFbar, typename ShapeFunction, typename ShapeMatricesType>

std::tuple< double, VectorTypeForFbar > ProcessLib::NonLinearFbar::computeFBarInitialVariablesElementCenter	(	bool const	compute_detF0_only,
		Eigen::VectorXd const &	u,
		MeshLib::Element const &	element,
		bool const	is_axially_symmetric )

Definition at line 149 of file NonLinearFbar.h.

{
    auto const shape_matrices_at_element_center =
        NumLib::initShapeMatricesAtElementCenter<
            ShapeFunction, ShapeMatricesType, DisplacementDim>(
            element, is_axially_symmetric);
 
    auto const N_0 = shape_matrices_at_element_center.N;
    auto const dNdx_0 = shape_matrices_at_element_center.dNdx;
 
    // For the 2D case the 33-component is needed (and the four entries
    // of the non-symmetric matrix); In 3d there are nine entries.
    GradientMatrixType G0(
        DisplacementDim * DisplacementDim + (DisplacementDim == 2 ? 1 : 0),
        ShapeFunction::NPOINTS * DisplacementDim);
 
    auto const x_coord =
        NumLib::interpolateXCoordinate<ShapeFunction, ShapeMatricesType>(
            element, N_0);
    Deformation::computeGMatrix<DisplacementDim, ShapeFunction::NPOINTS>(
        dNdx_0, G0, is_axially_symmetric, N_0, x_coord);
 
    GradientVectorType const grad_u = G0 * u;
    GradientVectorType const F0 =
        grad_u + MathLib::VectorizedTensor::identity<DisplacementDim>();
 
    if (compute_detF0_only)
    {
        return {MathLib::VectorizedTensor::determinant(F0), {}};
    }
 
    VectorTypeForFbar F0InvN =
        computeQVector<DisplacementDim, ShapeFunction::NPOINTS,
                       VectorTypeForFbar, GradientVectorType,
                       typename ShapeMatricesType::GlobalDimNodalMatrixType>(
            dNdx_0, F0, is_axially_symmetric);
 
    return {MathLib::VectorizedTensor::determinant(F0), F0InvN};
}

References ProcessLib::Deformation::computeGMatrix(), computeQVector(), MathLib::VectorizedTensor::determinant(), MathLib::VectorizedTensor::identity(), NumLib::initShapeMatricesAtElementCenter(), and NumLib::interpolateXCoordinate().

Referenced by ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::computeFBarVariables().

◆ computeQVector()

template<int DisplacementDim, int NPOINTS, typename VectorTypeForFbar, typename GradientVectorType, typename DNDX_Type>

VectorTypeForFbar ProcessLib::NonLinearFbar::computeQVector	(	DNDX_Type const &	dNdx,
		GradientVectorType const &	F,
		bool const	)

Definition at line 35 of file NonLinearFbar.h.

{
    // Note: MathLib::VectorizedTensor::toTensor returns a (3 x 3) matrix by
    //       Eigen::Matrix3D. The size of the return matrix is wrong, which
    //       should be (DisplacementDim x DisplacementDim).
    // Besides, Eigen::Matrix member inverse does not work with the mapped
    // matrix, e.g. that by using topLeftCorner. Therefore a hard copy is
    // needed for the following inverse, i.e. declaring F0inv as type
    // Eigen::MatrixXd instead of type auto.
    //
    Eigen::MatrixXd const F_matrix =
        MathLib::VectorizedTensor::toTensor<DisplacementDim>(F)
            .template topLeftCorner<DisplacementDim, DisplacementDim>();
    auto const Finv = F_matrix.inverse();
 
    VectorTypeForFbar FInvN =
        VectorTypeForFbar::Zero(DisplacementDim * NPOINTS);
    for (int i = 0; i < NPOINTS; ++i)
    {
        auto const dNidx = dNdx.col(i);
        for (int k = 0; k < DisplacementDim; k++)
        {
            FInvN[k * NPOINTS + i] = Finv.col(k).dot(dNidx);
        }
    }
 
    return FInvN;
}

References MathLib::VectorizedTensor::toTensor().

Referenced by ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::assembleWithJacobian(), and computeFBarInitialVariablesElementCenter().

◆ convertStringToDetFBarType()

BarDetFType ProcessLib::NonLinearFbar::convertStringToDetFBarType ( std::string_view const bar_det_f_type_name )

Definition at line 11 of file NonLinearFbar.cpp.

{
    if (boost::iequals(bar_det_f_type_name, "element_center_value"))
    {
        INFO(
            "Use F bar method with the element center value of F for the "
            "det(F) modification.");
        return BarDetFType::ELEMENT_CENTER_VALUE;
    }
    if (boost::iequals(bar_det_f_type_name, "element_average"))
    {
        INFO(
            "Use F bar method with the element average of F for the det(F) "
            "modification.");
        return BarDetFType::ELEMENT_AVERAGE;
    }
 
    return BarDetFType::NONE;
}

References ELEMENT_AVERAGE, ELEMENT_CENTER_VALUE, INFO(), and NONE.

Referenced by ProcessLib::LargeDeformation::createLargeDeformationProcess().

◆ identityForF()

template<int DisplacementDim>

MathLib::KelvinVector::KelvinVectorType< DisplacementDim > ProcessLib::NonLinearFbar::identityForF ( bool const is_axially_symmetric )

Definition at line 192 of file NonLinearFbar.h.

{
    if (DisplacementDim == 3 || is_axially_symmetric)
    {
        return MathLib::KelvinVector::Invariants<
            MathLib::KelvinVector::kelvin_vector_dimensions(
                DisplacementDim)>::identity2;
    }
 
    auto identity2 = MathLib::KelvinVector::Invariants<
        MathLib::KelvinVector::kelvin_vector_dimensions(
            DisplacementDim)>::identity2;
 
    identity2[2] = 0.0;
 
    return identity2;
}

References MathLib::KelvinVector::kelvin_vector_dimensions().

Referenced by compute2EPlusI(), and ProcessLib::LargeDeformation::LargeDeformationLocalAssembler< ShapeFunction, DisplacementDim >::computeOutputStrainData().

Enumerations

Functions