NodalRankTwoComponentPD.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 "NodalRankTwoComponentPD.h"
#include "RankTwoTensor.h"
#include "AuxiliarySystem.h"

registerMooseObject("PeridynamicsApp", NodalRankTwoComponentPD);

InputParameters
NodalRankTwoComponentPD::validParams()
{
  InputParameters params = NodalRankTwoUserObjectBasePD::validParams();
  params.addClassDescription(
      "Class for calculating components of nodal rank-two stress and strain tensors "
      "from material properties (stress and strain) for edge elements (i.e., "
      "bonds) connected at that node. NOTE: This UserObject only applies to the NOSPD model.");

  params.addRequiredRangeCheckedParam<unsigned int>(
      "index_i",
      "index_i >= 0 & index_i <= 2",
      "The index i of ij for the tensor to output (0, 1, 2)");
  params.addRequiredRangeCheckedParam<unsigned int>(
      "index_j",
      "index_j >= 0 & index_j <= 2",
      "The index j of ij for the tensor to output (0, 1, 2)");

  return params;
}

NodalRankTwoComponentPD::NodalRankTwoComponentPD(const InputParameters & parameters)
  : NodalRankTwoUserObjectBasePD(parameters),
    _i(getParam<unsigned int>("index_i")),
    _j(getParam<unsigned int>("index_j"))
{
}

void
NodalRankTwoComponentPD::gatherWeightedValue(unsigned int id, dof_id_type dof, Real dg_vol_frac)
{
  _aux.solution().add(dof, _tensor[id](_i, _j) * dg_vol_frac);
}