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

registerMooseObject("PeridynamicsApp", ComputeSmallStrainConstantHorizonMaterialOSPD);

InputParameters
ComputeSmallStrainConstantHorizonMaterialOSPD::validParams()
{
  InputParameters params = ComputeSmallStrainMaterialBaseOSPD::validParams();
  params.addClassDescription(
      "Class for computing peridynamic micro elastic moduli for ordinary state-based model "
      "using regular uniform mesh");

  return params;
}

ComputeSmallStrainConstantHorizonMaterialOSPD::ComputeSmallStrainConstantHorizonMaterialOSPD(
    const InputParameters & parameters)
  : ComputeSmallStrainMaterialBaseOSPD(parameters)
{
}

void
ComputeSmallStrainConstantHorizonMaterialOSPD::computePeridynamicsParams()
{
  _a = 0.5 * (_bulk_modulus - (8.0 - _dim) / 3.0 * _shear_modulus);

  // _b = 2 * _b * _horizon_(i/j) * _origin_vec.norm() //_origin_vec.norm() will be cancelled out in
  // parent material model
  _b = _origin_vec.norm() * (3.0 * _dim + 6.0) * _shear_modulus / M_PI /
       std::pow(_horizon_radius[0], _dim + 1);

  // _d_i = _di * _horizon_(i/j)
  _d[0] = (_dim / 4.0 + 1.5) / M_PI / std::pow(_horizon_radius[0], _dim);
  _d[1] = _d[0];
}