ZeroBondStatusUserObjectPD.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 "ZeroBondStatusUserObjectPD.h"
#include "AuxiliarySystem.h"
 
registerMooseObject("PeridynamicsTestApp", ZeroBondStatusUserObjectPD);
 
InputParameters
ZeroBondStatusUserObjectPD::validParams()
{
  InputParameters params = GeneralUserObjectBasePD::validParams();
  params.addClassDescription("Class to assign value of zero to auxilary variable bond status for "
                             "user specified bonds list");
 
  params.addRequiredParam<std::vector<unsigned int>>(
      "bond_ids_list", "List of bond IDs to be assigned with value 0");
 
  params.set<ExecFlagEnum>("execute_on", true) = EXEC_INITIAL;
 
  return params;
}
 
ZeroBondStatusUserObjectPD::ZeroBondStatusUserObjectPD(const InputParameters & parameters)
  : GeneralUserObjectBasePD(parameters), _list(getParam<std::vector<unsigned int>>("bond_ids_list"))
{
}
 
void
ZeroBondStatusUserObjectPD::initialize()
{
  _aux.solution().close();
}
 
void
ZeroBondStatusUserObjectPD::execute()
{
  for (const auto & bid : _list)
  {
    const Elem * current_bond = _pdmesh.elemPtr(bid);
    if (current_bond->processor_id() == _tid)
    {
      dof_id_type dof = current_bond->dof_number(_aux.number(), _bond_status_var->number(), 0);
      _aux.solution().set(dof, 0); // set the status of bonds in the given list to 0
    }
  }
}
 
void
ZeroBondStatusUserObjectPD::finalize()
{
  _aux.solution().close();
}