Detailed Description

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>
class ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >

The normal traction boundary condition is a special type of Neumann boundary condition where the given value is applied in the direction of the element's normal vector \(\mathbf{n}\):

\[ \bar{t} := \sigma \mathbf{n} = p \mathbf{n}, \]

where \(p\) is the value on the boundary given by the parameter tag.

Definition at line 34 of file NormalTractionBoundaryCondition.h.

#include <NormalTractionBoundaryCondition.h>

Inheritance diagram for ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >:

[legend]

Collaboration diagram for ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >:

[legend]

Public Member Functions
	NormalTractionBoundaryCondition (unsigned const integration_order, unsigned const shapefunction_order, MeshLib::Mesh const &bulk_mesh, NumLib::LocalToGlobalIndexMap const &dof_table_bulk, int const variable_id, MeshLib::Mesh const &bc_mesh, ParameterLib::Parameter< double > const &pressure)
void	applyNaturalBC (const double t, std::vector< GlobalVector * > const &x, int const process_id, GlobalMatrix K, GlobalVector &b, GlobalMatrix Jac) override
Public Member Functions inherited from ProcessLib::BoundaryCondition
virtual void	getEssentialBCValues (const double, GlobalVector const &, NumLib::IndexValueVector< GlobalIndexType > &) const
	Writes the values of essential BCs to `bc_values`.
virtual void	preTimestep (const double, std::vector< GlobalVector * > const &, int const)
virtual void	postTimestep (const double, std::vector< GlobalVector * > const &, int const)
virtual	~BoundaryCondition ()=default

Private Attributes
MeshLib::Mesh const &	_bc_mesh
std::vector< MeshLib::Node * >	_nodes_subset
std::unique_ptr< MeshLib::MeshSubset const >	_mesh_subset_all_nodes
std::unique_ptr< NumLib::LocalToGlobalIndexMap >	_dof_table_boundary
unsigned const	_integration_order
	Integration order for integration over the lower-dimensional elements.
std::vector< std::unique_ptr< NormalTractionBoundaryConditionLocalAssemblerInterface > >	_local_assemblers
	Local assemblers for each element of number of _elements.
ParameterLib::Parameter< double > const &	_pressure
std::vector< Eigen::Vector3d >	_element_normals
	Normal vectors for each element in the boundary condition mesh.

Constructor & Destructor Documentation

◆ NormalTractionBoundaryCondition()

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::NormalTractionBoundaryCondition	(	unsigned const	integration_order,
		unsigned const	shapefunction_order,
		MeshLib::Mesh const &	bulk_mesh,
		NumLib::LocalToGlobalIndexMap const &	dof_table_bulk,
		int const	variable_id,
		MeshLib::Mesh const &	bc_mesh,
		ParameterLib::Parameter< double > const &	pressure )

Create a boundary condition process from given config, DOF-table, and a mesh subset for a given variable and its component. A local DOF-table, a subset of the given one, is constructed.

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

    : _bc_mesh(bc_mesh),
      _integration_order(integration_order),
      _pressure(pressure)
{
    // Create component ids vector for the current variable.
    auto const& number_of_components =
        dof_table_bulk.getNumberOfVariableComponents(variable_id);
    std::vector<int> component_ids(number_of_components);
    std::iota(std::begin(component_ids), std::end(component_ids), 0);
 
    // BC mesh subset creation
    std::vector<MeshLib::Node*> const bc_nodes = _bc_mesh.getNodes();
    DBUG("Found {:d} nodes for Natural BCs for the variable {:d}",
         bc_nodes.size(), variable_id);
 
    MeshLib::MeshSubset bc_mesh_subset(_bc_mesh, bc_nodes);
 
    // Compute normal vectors for each element in the boundary condition mesh.
    auto const* const bulk_element_ids = MeshLib::bulkElementIDs(_bc_mesh);
    assert(bulk_element_ids != nullptr);
    auto const& elements = _bc_mesh.getElements();
    _element_normals =
        elements |
        ranges::views::transform(
            [&](const MeshLib::Element* e_ptr)
            {
                // Compute the corresponding bulk element for normal orientation
                auto const* bulk_element =
                    bulk_mesh.getElement((*bulk_element_ids)[e_ptr->getID()]);
                assert(bulk_element != nullptr);
 
                return computeElementNormal<GlobalDim>(*e_ptr, *bulk_element);
            }) |
        ranges::to<std::vector<Eigen::Vector3d>>();
 
    // Create local DOF table from the BC mesh subset for the given variable and
    // component ids.
    _dof_table_boundary = dof_table_bulk.deriveBoundaryConstrainedMap(
        variable_id, component_ids, std::move(bc_mesh_subset));
 
    BoundaryConditionAndSourceTerm::detail::createLocalAssemblers<
        GlobalDim, LocalAssemblerImplementation>(
        *_dof_table_boundary, shapefunction_order, _bc_mesh.getElements(),
        _local_assemblers, NumLib::IntegrationOrder{integration_order},
        _bc_mesh.isAxiallySymmetric(), _pressure, _element_normals);
}

References NormalTractionBoundaryCondition(), _bc_mesh, _dof_table_boundary, _element_normals, _integration_order, _local_assemblers, _pressure, MeshLib::bulkElementIDs(), ProcessLib::NormalTractionBoundaryCondition::computeElementNormal(), ProcessLib::BoundaryConditionAndSourceTerm::detail::createLocalAssemblers(), DBUG(), and MeshLib::Element::getID().

Referenced by NormalTractionBoundaryCondition(), and applyNaturalBC().

Member Function Documentation

◆ applyNaturalBC()

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

void ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::applyNaturalBC	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		int const	process_id,
		GlobalMatrix *	K,
		GlobalVector &	b,
		GlobalMatrix *	Jac )

overridevirtual

Calls local assemblers which calculate their contributions to the global matrix and the right-hand-side.

Reimplemented from ProcessLib::BoundaryCondition.

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

{
    GlobalExecutor::executeMemberOnDereferenced(
        &NormalTractionBoundaryConditionLocalAssemblerInterface::assemble,
        _local_assemblers, *_dof_table_boundary, t, x, K, b, Jac);
}

References NormalTractionBoundaryCondition(), _dof_table_boundary, _local_assemblers, applyNaturalBC(), ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryConditionLocalAssemblerInterface::assemble(), and NumLib::SerialExecutor::executeMemberOnDereferenced().

Referenced by applyNaturalBC().

Member Data Documentation

◆ _bc_mesh

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

MeshLib::Mesh const& ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::_bc_mesh

private

Definition at line 54 of file NormalTractionBoundaryCondition.h.

Referenced by NormalTractionBoundaryCondition().

◆ _dof_table_boundary

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

std::unique_ptr<NumLib::LocalToGlobalIndexMap> ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::_dof_table_boundary

private

Local dof table, a subset of the global one restricted to the participating number of _elements of the boundary condition.

Definition at line 64 of file NormalTractionBoundaryCondition.h.

Referenced by NormalTractionBoundaryCondition(), and applyNaturalBC().

◆ _element_normals

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

std::vector<Eigen::Vector3d> ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::_element_normals

private

Normal vectors for each element in the boundary condition mesh.

Definition at line 77 of file NormalTractionBoundaryCondition.h.

Referenced by NormalTractionBoundaryCondition().

◆ _integration_order

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

unsigned const ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::_integration_order

private

Integration order for integration over the lower-dimensional elements.

Definition at line 67 of file NormalTractionBoundaryCondition.h.

Referenced by NormalTractionBoundaryCondition().

◆ _local_assemblers

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

std::vector< std::unique_ptr<NormalTractionBoundaryConditionLocalAssemblerInterface> > ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::_local_assemblers

private

Local assemblers for each element of number of _elements.

Definition at line 72 of file NormalTractionBoundaryCondition.h.

Referenced by NormalTractionBoundaryCondition(), and applyNaturalBC().

◆ _mesh_subset_all_nodes

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

std::unique_ptr<MeshLib::MeshSubset const> ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::_mesh_subset_all_nodes

private

Definition at line 60 of file NormalTractionBoundaryCondition.h.

◆ _nodes_subset

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

std::vector<MeshLib::Node*> ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::_nodes_subset

private

Intersection of boundary nodes and bulk mesh subset for the variable_id/component_id pair.

Definition at line 58 of file NormalTractionBoundaryCondition.h.

◆ _pressure

template<int GlobalDim, template< typename, int > class LocalAssemblerImplementation>

ParameterLib::Parameter<double> const& ProcessLib::NormalTractionBoundaryCondition::NormalTractionBoundaryCondition< GlobalDim, LocalAssemblerImplementation >::_pressure

private

Definition at line 74 of file NormalTractionBoundaryCondition.h.

Referenced by NormalTractionBoundaryCondition().

The documentation for this class was generated from the following files:

ProcessLib/BoundaryConditionAndSourceTerm/NormalTractionBoundaryCondition.h
ProcessLib/BoundaryConditionAndSourceTerm/NormalTractionBoundaryCondition-impl.h

Detailed Description

Public Member Functions

Private Attributes

Constructor & Destructor Documentation

◆ NormalTractionBoundaryCondition()

Member Function Documentation

◆ applyNaturalBC()

Member Data Documentation

◆ _bc_mesh

◆ _dof_table_boundary

◆ _element_normals

◆ _integration_order

◆ _local_assemblers

◆ _mesh_subset_all_nodes

◆ _nodes_subset

◆ _pressure