Detailed Description

Definition at line 20 of file PhaseFieldIrreversibleDamageOracleBoundaryCondition.h.

#include <PhaseFieldIrreversibleDamageOracleBoundaryCondition.h>

Inheritance diagram for ProcessLib::PhaseFieldIrreversibleDamageOracleBoundaryCondition:

Collaboration diagram for ProcessLib::PhaseFieldIrreversibleDamageOracleBoundaryCondition:

Public Member Functions
	PhaseFieldIrreversibleDamageOracleBoundaryCondition (NumLib::LocalToGlobalIndexMap const &dof_table, MeshLib::Mesh const &mesh, int const variable_id, int const component_id)

void	getEssentialBCValues (const double t, const GlobalVector &x, NumLib::IndexValueVector< GlobalIndexType > &bc_values) const override
	Writes the values of essential BCs to `bc_values`. More...

void	preTimestep (const double t, std::vector< GlobalVector * > const &x, int const process_id) override

Public Member Functions inherited from ProcessLib::BoundaryCondition
virtual void	applyNaturalBC (const double, std::vector< GlobalVector * > const &, int const, GlobalMatrix &, GlobalVector &, GlobalMatrix *)

virtual void	postTimestep (const double, std::vector< GlobalVector * > const &, int const)

virtual	~BoundaryCondition ()=default

Private Attributes
NumLib::LocalToGlobalIndexMap const &	_dof_table

MeshLib::Mesh const &	_mesh

int const	_variable_id

int const	_component_id

NumLib::IndexValueVector< GlobalIndexType >	_bc_values

Constructor & Destructor Documentation

◆ PhaseFieldIrreversibleDamageOracleBoundaryCondition()

ProcessLib::PhaseFieldIrreversibleDamageOracleBoundaryCondition::PhaseFieldIrreversibleDamageOracleBoundaryCondition	(	NumLib::LocalToGlobalIndexMap const &	dof_table,
		MeshLib::Mesh const &	mesh,
		int const	variable_id,
		int const	component_id
	)

inline

Definition at line 24 of file PhaseFieldIrreversibleDamageOracleBoundaryCondition.h.

         : _dof_table(dof_table),
           _mesh(mesh),
           _variable_id(variable_id),
           _component_id(component_id)
     {
         if (variable_id >= static_cast<int>(dof_table.getNumberOfVariables()) ||
             component_id >=
                 dof_table.getNumberOfVariableComponents(variable_id))
         {
             OGS_FATAL(
                 "Variable id or component id too high. Actual values: ({:d}, "
                 "{:d}), "
                 "maximum values: ({:d}, {:d}).",
                 variable_id, component_id, dof_table.getNumberOfVariables(),
                 dof_table.getNumberOfVariableComponents(variable_id));
         }
     }

References NumLib::LocalToGlobalIndexMap::getNumberOfVariableComponents(), NumLib::LocalToGlobalIndexMap::getNumberOfVariables(), and OGS_FATAL.

Member Function Documentation

◆ getEssentialBCValues()

void ProcessLib::PhaseFieldIrreversibleDamageOracleBoundaryCondition::getEssentialBCValues	(	const double	,
		const GlobalVector &	,
		NumLib::IndexValueVector< GlobalIndexType > &
	)		const

overridevirtual

Writes the values of essential BCs to bc_values.

Reimplemented from ProcessLib::BoundaryCondition.

Definition at line 21 of file PhaseFieldIrreversibleDamageOracleBoundaryCondition.cpp.

 {
     bc_values.ids.clear();
     bc_values.values.clear();
  
     // convert mesh node ids to global index for the given component
     bc_values.ids.reserve(bc_values.ids.size() + _bc_values.ids.size());
     bc_values.values.reserve(bc_values.values.size() +
                              _bc_values.values.size());
  
     std::copy(_bc_values.ids.begin(), _bc_values.ids.end(),
               std::back_inserter(bc_values.ids));
     std::copy(_bc_values.values.begin(), _bc_values.values.end(),
               std::back_inserter(bc_values.values));
 }

References _bc_values, MathLib::LinAlg::copy(), NumLib::IndexValueVector< IndexType >::ids, and NumLib::IndexValueVector< IndexType >::values.

◆ preTimestep()

void ProcessLib::PhaseFieldIrreversibleDamageOracleBoundaryCondition::preTimestep	(	const double	t,
		std::vector< GlobalVector * > const &	x,
		int const	process_id
	)

overridevirtual

Reimplemented from ProcessLib::BoundaryCondition.

Definition at line 40 of file PhaseFieldIrreversibleDamageOracleBoundaryCondition.cpp.

 {
     // phase-field variable is considered irreversible if it loses more than 95%
     // of the stiffness, which is a widely used threshold.
     double irreversibleDamage = 0.05;
  
     _bc_values.ids.clear();
     _bc_values.values.clear();
  
     auto const mesh_id = _mesh.getID();
     auto const& nodes = _mesh.getNodes();
     for (auto const* n : nodes)
     {
         std::size_t node_id = n->getID();
         MeshLib::Location l(mesh_id, MeshLib::MeshItemType::Node, node_id);
         const auto g_idx =
             _dof_table.getGlobalIndex(l, _variable_id, _component_id);
  
         if ((*x[process_id])[node_id] <= irreversibleDamage)
         {
             _bc_values.ids.emplace_back(g_idx);
             _bc_values.values.emplace_back(0.0);
         }
     }
 }