ProcessLib::NormalTractionBoundaryCondition Namespace Reference

Classes
struct	IntegrationPointData
class	NormalTractionBoundaryCondition
class	NormalTractionBoundaryConditionLocalAssembler
class	NormalTractionBoundaryConditionLocalAssemblerInterface

Functions
template<int GlobalDim>
Eigen::Vector3d	computeElementNormal (const MeshLib::Element &element, const MeshLib::Element &bulk_element)
template<int GlobalDim>
std::unique_ptr< NormalTractionBoundaryCondition< GlobalDim, NormalTractionBoundaryConditionLocalAssembler > >	createNormalTractionBoundaryCondition (BaseLib::ConfigTree const &config, MeshLib::Mesh const &bc_mesh, MeshLib::Mesh const &bulk_mesh, NumLib::LocalToGlobalIndexMap const &dof_table, int const variable_id, unsigned const integration_order, unsigned const shapefunction_order, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters)
template<int GlobalDim>
std::unique_ptr< NormalTractionBoundaryCondition< GlobalDim, NormalTractionBoundaryConditionLocalAssembler > >	createNormalTractionBoundaryCondition (BaseLib::ConfigTree const &config, MeshLib::Mesh const &bc_mesh, NumLib::LocalToGlobalIndexMap const &dof_table, int const variable_id, unsigned const integration_order, unsigned const shapefunction_order, std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &parameters)

Function Documentation

computeElementNormal()

template<int GlobalDim>

Eigen::Vector3d ProcessLib::NormalTractionBoundaryCondition::computeElementNormal	(	const MeshLib::Element &	element,
		const MeshLib::Element &	bulk_element )

Computes the normal vector for a boundary element.

Parameters

element	The boundary element.
bulk_element	The corresponding bulk element.

Template Parameters

GlobalDim

Global dimension of the problem.

Returns

The computed normal vector.

Definition at line 34 of file NormalTractionBoundaryCondition-impl.h.

{
    if (element.getGeomType() == MeshLib::MeshElemType::LINE)
    {
        auto const n = MeshLib::FaceRule::getSurfaceNormal(bulk_element);
        Eigen::Vector3d const e = element.getNode(1)->asEigenVector3d() -
                                  element.getNode(0)->asEigenVector3d();
        Eigen::Vector3d const normal = e.cross(n).normalized();
 
        // Compute center of the bulk element to correctly orient the
        // normal.
        auto const c =
            MeshLib::getCenterOfGravity(bulk_element).asEigenVector3d();
        // Flip n if the normal points toward c:
        // <normal, c - n0> > 0.
        if (normal.dot(c - element.getNode(0)->asEigenVector3d()) > 0)
        {
            return -normal;
        }
        return normal;
    }
 
    Eigen::Vector3d normal = MeshLib::FaceRule::getSurfaceNormal(element);
    normal.tail<3 - GlobalDim>().setZero();
    return normal.normalized();
}

References MathLib::Point3d::asEigenVector3d(), MeshLib::getCenterOfGravity(), MeshLib::Element::getGeomType(), MeshLib::Element::getNode(), MeshLib::FaceRule::getSurfaceNormal(), and MeshLib::LINE.

createNormalTractionBoundaryCondition() [1/2]

template<int GlobalDim>

std::unique_ptr< NormalTractionBoundaryCondition< GlobalDim, NormalTractionBoundaryConditionLocalAssembler > > ProcessLib::NormalTractionBoundaryCondition::createNormalTractionBoundaryCondition	(	BaseLib::ConfigTree const &	config,
		MeshLib::Mesh const &	bc_mesh,
		MeshLib::Mesh const &	bulk_mesh,
		NumLib::LocalToGlobalIndexMap const &	dof_table,
		int const	variable_id,
		unsigned const	integration_order,
		unsigned const	shapefunction_order,
		std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &	parameters )

Input File Parameter

prj__process_variables__process_variable__boundary_conditions__boundary_condition__type

Input File Parameter

prj__process_variables__process_variable__boundary_conditions__boundary_condition__NormalTraction__parameter

Definition at line 133 of file NormalTractionBoundaryCondition-impl.h.

{
    DBUG("Constructing NormalTractionBoundaryCondition from config.");
    config.checkConfigParameter("type", "NormalTraction");
 
    auto const parameter_name =
        config.getConfigParameter<std::string>("parameter");
    DBUG("Using parameter {:s}", parameter_name);
 
    auto const& pressure = ParameterLib::findParameter<double>(
        parameter_name, parameters, 1, &bc_mesh);
    return std::make_unique<NormalTractionBoundaryCondition<
        GlobalDim, NormalTractionBoundaryConditionLocalAssembler>>(
        integration_order, shapefunction_order, bulk_mesh, dof_table,
        variable_id, bc_mesh, pressure);
}

References BaseLib::ConfigTree::checkConfigParameter(), DBUG(), ParameterLib::findParameter(), and BaseLib::ConfigTree::getConfigParameter().

Referenced by ProcessLib::createBoundaryCondition().

createNormalTractionBoundaryCondition() [2/2]

template<int GlobalDim>

std::unique_ptr< NormalTractionBoundaryCondition< GlobalDim, NormalTractionBoundaryConditionLocalAssembler > > ProcessLib::NormalTractionBoundaryCondition::createNormalTractionBoundaryCondition	(	BaseLib::ConfigTree const &	config,
		MeshLib::Mesh const &	bc_mesh,
		NumLib::LocalToGlobalIndexMap const &	dof_table,
		int const	variable_id,
		unsigned const	integration_order,
		unsigned const	shapefunction_order,
		std::vector< std::unique_ptr< ParameterLib::ParameterBase > > const &	parameters )