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

registerMooseObject("PeridynamicsApp", ComputeSmallStrainVariableHorizonMaterialOSPD);

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

  return params;
}

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

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

  // _b = _bij * _horizon_i + _bji * _horizon_j
  _b = _dim * _dim * (_bulk_modulus / 2.0 - _a) * (1.0 / _horizon_vol[0] + 1.0 / _horizon_vol[1]);

  // _d_i = _di * _horizon_i = _dim / _nvsum_i
  _d[0] = _dim / _horizon_vol[0];

  // _d_j = _dj * _horizon_j = _dim / _nvsum_j
  _d[1] = _dim / _horizon_vol[1];
}