RBMPresetOldValuePD.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#include "RBMPresetOldValuePD.h"
 
registerMooseObject("PeridynamicsApp", RBMPresetOldValuePD);
 
InputParameters
RBMPresetOldValuePD::validParams()
{
  InputParameters params = DirichletBCBase::validParams();
  params.addClassDescription("Class to apply a preset BC to nodes with rigid body motion (RBM).");
 
  // Forcefully preset the BC
  params.set<bool>("preset") = true;
  params.suppressParameter<bool>("preset");
 
  return params;
}
 
RBMPresetOldValuePD::RBMPresetOldValuePD(const InputParameters & parameters)
  : DirichletBCBase(parameters),
    _pdmesh(dynamic_cast<PeridynamicsMesh &>(_mesh)),
    _u_old(_var.dofValuesOld()),
    _bond_status_var(&_subproblem.getStandardVariable(_tid, "bond_status"))
{
}
 
Real
RBMPresetOldValuePD::computeQpValue()
{
  return _u_old[_qp];
}
 
bool
RBMPresetOldValuePD::shouldApply() const
{
  // check whether the number of active bonds is less than number of problem dimension
  unsigned int active_bonds = 0;
  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_node->id());
 
  for (unsigned int nb = 0; nb < bonds.size(); ++nb)
    if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[nb])) > 0.5)
      active_bonds++;
 
  return (active_bonds < _pdmesh.dimension() + 1);
}