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`.

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, NumLib::IndexValueVector< typename >::ids, and NumLib::IndexValueVector< typename >::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 (g_idx < 0)
        {
            continue;
        }
 
        if ((*x[process_id])[g_idx] <= irreversibleDamage)
        {
            _bc_values.ids.emplace_back(g_idx);
            _bc_values.values.emplace_back(0.0);
        }
    }
}