BondStatusBasePD.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "BondStatusBasePD.h"
#include "PeridynamicsMesh.h"
 
template <>
InputParameters
validParams<BondStatusBasePD>()
{
  InputParameters params = validParams<AuxKernelBasePD>();
  params.addClassDescription("Base class for different failure criteria to update the bond status");
 
  params.addRequiredCoupledVar("critical_variable", "Name of critical AuxVariable");
 
  params.set<ExecFlagEnum>("execute_on", true) = EXEC_TIMESTEP_END;
 
  return params;
}
 
BondStatusBasePD::BondStatusBasePD(const InputParameters & parameters)
  : AuxKernelBasePD(parameters),
    _bond_status_var(&_subproblem.getStandardVariable(_tid, "bond_status")),
    _critical_val(coupledValue("critical_variable"))
{
}
 
Real
BondStatusBasePD::computeValue()
{
  Real val = 0.0;
 
  if (_bond_status_var->getElementalValue(_current_elem) > 0.5) // unbroken bond
  {
    Real failure_criterion_val = computeFailureCriterionValue();
    if (failure_criterion_val < 0.0) // unmet failure criterion
      val = 1.0;                     // bond is still unbroken
  }
 
  return val;
}