Rank two tensor based failure ctriteria to update the bond status for non-ordinary state-based model. More...

#include <RankTwoBasedFailureCriteriaNOSPD.h>

Inheritance diagram for RankTwoBasedFailureCriteriaNOSPD:

Public Member Functions
	RankTwoBasedFailureCriteriaNOSPD (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeFailureCriterionValue () override
	Function to compute the value of a given failure criterion. More...

virtual Real	computeValue () override

Protected Attributes
MooseEnum	_failure_criterion
	MooseEnum used to control which failure criterion to use. More...

const MaterialProperty< RankTwoTensor > *	_tensor
	Material property containing the rank two tensor. More...

MooseVariable *	_bond_status_var
	Bond_status variable. More...

const VariableValue &	_critical_val
	Critical AuxVariable. More...

PeridynamicsMesh &	_pdmesh
	Reference to peridynamic mesh object. More...

const unsigned int	_dim
	Problem dimension. More...

Detailed Description

Rank two tensor based failure ctriteria to update the bond status for non-ordinary state-based model.

Definition at line 24 of file RankTwoBasedFailureCriteriaNOSPD.h.

Constructor & Destructor Documentation

◆ RankTwoBasedFailureCriteriaNOSPD()

RankTwoBasedFailureCriteriaNOSPD::RankTwoBasedFailureCriteriaNOSPD ( const InputParameters & parameters )

Definition at line 34 of file RankTwoBasedFailureCriteriaNOSPD.C.

   : BondStatusBasePD(parameters),
     _failure_criterion(getParam<MooseEnum>("failure_criterion")),
     _tensor(nullptr)
 {
   if (hasMaterialProperty<RankTwoTensor>("rank_two_tensor"))
     _tensor = &getMaterialProperty<RankTwoTensor>("rank_two_tensor");
   else
     mooseError("Error in RankTwoBasedFailureCriteriaNOSPD! Required rank two tensor is not "
                "available for current peridynamics model!");
 }

Member Function Documentation

◆ computeFailureCriterionValue()

Real RankTwoBasedFailureCriteriaNOSPD::computeFailureCriterionValue ( )

overrideprotectedvirtual

Function to compute the value of a given failure criterion.

Positive value indicates the critical value is exceeded. Otherwise, the critical value is not reached yet

Returns: failure criterion value

Implements BondStatusBasePD.

Definition at line 48 of file RankTwoBasedFailureCriteriaNOSPD.C.

 {
   Real val = 0.0;
   RankTwoTensor avg_tensor = 0.5 * ((*_tensor)[0] + (*_tensor)[1]);
   Point dirc;
  
   switch (_failure_criterion)
   {
     case 0:
       val = RankTwoScalarTools::axialStress(
           avg_tensor, *_current_elem->node_ptr(0), *_current_elem->node_ptr(1), dirc);
       break;
     case 1:
       val = RankTwoScalarTools::maxPrincipal(avg_tensor, dirc);
       break;
     case 2:
       val = 4.0 / 3.0 * RankTwoScalarTools::maxShear(avg_tensor);
       break;
     case 3:
       val = 2.0 * RankTwoScalarTools::maxShear(avg_tensor);
       break;
     case 4:
       val = RankTwoScalarTools::effectiveStrain(avg_tensor);
       break;
     case 5:
       val = RankTwoScalarTools::vonMisesStress(avg_tensor);
       break;
     default:
       mooseError("Unsupported rank two tensor-based failure criterion. Choose from: Axial "
                  "MaxPrincipal TrescaStrain TrescaStress VonMisesStrain VonMisesStress");
   }
  
   return val - _critical_val[0];
 }

◆ computeValue()

Real BondStatusBasePD::computeValue ( )

overrideprotectedvirtualinherited

Definition at line 35 of file BondStatusBasePD.C.

 {
   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;
 }