ComputeSmallStrainNOSPD.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 "ComputeSmallStrainNOSPD.h"

registerMooseObject("PeridynamicsApp", ComputeSmallStrainNOSPD);

InputParameters
ComputeSmallStrainNOSPD::validParams()
{
  InputParameters params = ComputeStrainBaseNOSPD::validParams();
  params.addClassDescription(
      "Class for computing nodal quantities for the residual and Jacobian calculation "
      "for the peridynamic correspondence models under small strain assumptions");

  return params;
}

ComputeSmallStrainNOSPD::ComputeSmallStrainNOSPD(const InputParameters & parameters)
  : ComputeStrainBaseNOSPD(parameters)
{
}

void
ComputeSmallStrainNOSPD::computeQpStrain()
{
  computeQpDeformationGradient();

  computeQpTotalStrain();

  _mechanical_strain[_qp] = _total_strain[_qp];

  // Subtract Eigen strains
  for (auto es : _eigenstrains)
    _mechanical_strain[_qp] -= (*es)[_qp];

  // zero out all strain measures for broken bond
  if (_bond_status_var->getElementalValue(_current_elem) < 0.5)
  {
    _mechanical_strain[_qp].zero();
    _total_strain[_qp].zero();
  }
}

void
ComputeSmallStrainNOSPD::computeQpTotalStrain()
{
  // the small strain tensor
  _total_strain[_qp] = 0.5 * (_deformation_gradient[_qp].transpose() + _deformation_gradient[_qp]) -
                       RankTwoTensor(RankTwoTensor::initIdentity);
}