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

registerMooseObject("PeridynamicsApp", NodalRankTwoScalarPD);

InputParameters
NodalRankTwoScalarPD::validParams()
{
  InputParameters params = NodalRankTwoUserObjectBasePD::validParams();
  params.addClassDescription(
      "Class for calculating scalar quantities 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.addRequiredParam<MooseEnum>(
      "scalar_type", RankTwoScalarTools::scalarOptions(), "Type of scalar output");
  params.addParam<Point>(
      "point1",
      Point(0, 0, 0),
      "Start point for axis used to calculate some cylindrical material tensor quantities");
  params.addParam<Point>("point2",
                         Point(0, 1, 0),
                         "End point for axis used to calculate some material tensor quantities");
  params.addParam<Point>("direction", Point(0, 0, 1), "Direction vector");

  return params;
}

NodalRankTwoScalarPD::NodalRankTwoScalarPD(const InputParameters & parameters)
  : NodalRankTwoUserObjectBasePD(parameters),
    _scalar_type(getParam<MooseEnum>("scalar_type")),
    _point1(parameters.get<Point>("point1")),
    _point2(parameters.get<Point>("point2")),
    _input_direction(parameters.get<Point>("direction") / parameters.get<Point>("direction").norm())
{
}

void
NodalRankTwoScalarPD::gatherWeightedValue(unsigned int id, dof_id_type dof, Real dg_vol_frac)
{
  Point curr_point = *_current_elem->node_ptr(id);
  Real scalar_value = RankTwoScalarTools::getQuantity(
      _tensor[id], _scalar_type, _point1, _point2, curr_point, _input_direction);

  _aux.solution().add(dof, scalar_value * dg_vol_frac);
}