Base material class for correspondence material model. More...

#include <ComputeStrainBaseNOSPD.h>

Inheritance diagram for ComputeStrainBaseNOSPD:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ComputeStrainBaseNOSPD (const InputParameters &parameters)

virtual void	initQpStatefulProperties () override

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

void	validateDerivativeMaterialPropertyBase (const std::string &base)

const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	computeProperties () override

virtual void	computeBondStretch () override

virtual void	computeQpDeformationGradient ()
	Function to compute deformation gradient for peridynamic correspondence model.

virtual void	computeConventionalQpDeformationGradient ()
	Function to compute conventional nonlocal deformation gradient.

virtual void	computeBondHorizonQpDeformationGradient ()
	Function to compute bond-associated horizon based deformation gradient.

virtual void	computeQpStrain ()=0
	Function to compute strain tensors.

Protected Attributes
const MooseEnum	_stabilization
	Option of stabilization scheme for correspondence material model: FORCE, BOND_HORIZON_I or BOND_HORIZON_II.

const bool	_plane_strain
	Plane strain problem or not, this is only used for mechanical stretch calculation.


const MaterialProperty< RankFourTensor > &	_Cijkl
	Material properties to fetch.

std::vector< MaterialPropertyName >	_eigenstrain_names

std::vector< const MaterialProperty< RankTwoTensor > * >	_eigenstrains


MaterialProperty< RankTwoTensor > &	_shape2
	Material properties to store.

MaterialProperty< RankTwoTensor > &	_deformation_gradient

MaterialProperty< RankTwoTensor > &	_ddgraddu

MaterialProperty< RankTwoTensor > &	_ddgraddv

MaterialProperty< RankTwoTensor > &	_ddgraddw

MaterialProperty< RankTwoTensor > &	_total_strain

MaterialProperty< RankTwoTensor > &	_mechanical_strain

MaterialProperty< Real > &	_multi

MaterialProperty< Real > &	_sf_coeff
	fictitious force coefficient for force stabilized model

Private Member Functions
bool	haveMaterialProperty (const std::string &prop_name)

std::vector< VariableName >	buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)

void	validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)

bool	isNotObjectVariable (const VariableName &name)

Private Attributes
FEProblemBase &	_dmi_fe_problem

Detailed Description

Base material class for correspondence material model.

Definition at line 19 of file ComputeStrainBaseNOSPD.h.

Constructor & Destructor Documentation

◆ ComputeStrainBaseNOSPD()

ComputeStrainBaseNOSPD::ComputeStrainBaseNOSPD ( const InputParameters & parameters )

Definition at line 37 of file ComputeStrainBaseNOSPD.C.

  : DerivativeMaterialInterface<MechanicsMaterialBasePD>(parameters),
    _stabilization(getParam<MooseEnum>("stabilization")),
    _plane_strain(getParam<bool>("plane_strain")),
    _Cijkl(getMaterialProperty<RankFourTensor>("elasticity_tensor")),
    _eigenstrain_names(getParam<std::vector<MaterialPropertyName>>("eigenstrain_names")),
    _eigenstrains(_eigenstrain_names.size()),
    _shape2(declareProperty<RankTwoTensor>("rank_two_shape_tensor")),
    _deformation_gradient(declareProperty<RankTwoTensor>("deformation_gradient")),
    _ddgraddu(declareProperty<RankTwoTensor>("ddeformation_gradient_du")),
    _ddgraddv(declareProperty<RankTwoTensor>("ddeformation_gradient_dv")),
    _ddgraddw(declareProperty<RankTwoTensor>("ddeformation_gradient_dw")),
    _total_strain(declareProperty<RankTwoTensor>("total_strain")),
    _mechanical_strain(declareProperty<RankTwoTensor>("mechanical_strain")),
    _multi(declareProperty<Real>("multi")),
    _sf_coeff(declareProperty<Real>("stabilization_force_coeff"))
{
  for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
  {
    _eigenstrain_names[i] = _eigenstrain_names[i];
    _eigenstrains[i] = &getMaterialProperty<RankTwoTensor>(_eigenstrain_names[i]);
  }
}

Member Function Documentation

◆ computeBondHorizonQpDeformationGradient()

void ComputeStrainBaseNOSPD::computeBondHorizonQpDeformationGradient ( )

protectedvirtual

Function to compute bond-associated horizon based deformation gradient.

Definition at line 163 of file ComputeStrainBaseNOSPD.C.

{
  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(_qp));
  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(_qp));
 
  dof_id_type nb_index =
      std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - _qp)) -
      neighbors.begin();
  std::vector<dof_id_type> dg_neighbors =
      _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(_qp), nb_index);
  Real dg_sub_vol = _pdmesh.getHorizonSubsetVolume(_current_elem->node_id(_qp), nb_index);
  Real dg_sub_vol_sum = _pdmesh.getHorizonSubsetVolumeSum(_current_elem->node_id(_qp));
 
  // calculate the shape tensor and prepare the deformation gradient tensor
  Real vol_nb, weight_nb;
  RealGradient origin_vec_nb, current_vec_nb;
 
  for (unsigned int nb = 0; nb < dg_neighbors.size(); ++nb)
    if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb]])) > 0.5)
    {
      vol_nb = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb]]);
      origin_vec_nb = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb]]) -
                      _pdmesh.getNodeCoord(_current_elem->node_id(_qp));
 
      for (unsigned int k = 0; k < _dim; ++k)
        current_vec_nb(k) =
            origin_vec_nb(k) +
            _disp_var[k]->getNodalValue(*_pdmesh.nodePtr(neighbors[dg_neighbors[nb]])) -
            _disp_var[k]->getNodalValue(*_current_elem->node_ptr(_qp));
 
      weight_nb = _horizon_radius[_qp] / origin_vec_nb.norm();
      for (unsigned int k = 0; k < _dim; ++k)
      {
        for (unsigned int l = 0; l < _dim; ++l)
        {
          _shape2[_qp](k, l) += weight_nb * origin_vec_nb(k) * origin_vec_nb(l) * vol_nb;
          _deformation_gradient[_qp](k, l) +=
              weight_nb * current_vec_nb(k) * origin_vec_nb(l) * vol_nb;
        }
        // calculate derivatives of deformation_gradient w.r.t displacements of node i
        _ddgraddu[_qp](0, k) += -weight_nb * origin_vec_nb(k) * vol_nb;
        _ddgraddv[_qp](1, k) += -weight_nb * origin_vec_nb(k) * vol_nb;
        if (_dim == 3)
          _ddgraddw[_qp](2, k) += -weight_nb * origin_vec_nb(k) * vol_nb;
      }
    }
  // finalize the deformation gradient and its derivatives
  if (_shape2[_qp].det() == 0.)
    computeConventionalQpDeformationGradient(); // this will affect the corresponding jacobian
                                                // calculation
  else
  {
    _deformation_gradient[_qp] *= _shape2[_qp].inverse();
    _ddgraddu[_qp] *= _shape2[_qp].inverse();
    _ddgraddv[_qp] *= _shape2[_qp].inverse();
    _ddgraddw[_qp] *= _shape2[_qp].inverse();
  }
 
  // force state multiplier
  if (_stabilization == "BOND_HORIZON_I")
    _multi[_qp] = _horizon_radius[_qp] / _origin_vec.norm() * _node_vol[0] * _node_vol[1] *
                  dg_sub_vol / _horizon_vol[_qp];
  else if (_stabilization == "BOND_HORIZON_II")
    _multi[_qp] = _node_vol[_qp] * dg_sub_vol / dg_sub_vol_sum;
}

Referenced by computeQpDeformationGradient().

◆ computeBondStretch()

void ComputeStrainBaseNOSPD::computeBondStretch ( )

overrideprotectedvirtual

Definition at line 241 of file ComputeStrainBaseNOSPD.C.

{
  _total_stretch[0] = _current_len / _origin_vec.norm() - 1.0;
  _total_stretch[1] = _total_stretch[0];
  _mechanical_stretch[0] = _total_stretch[0];
 
  Real factor = 1.0;
  if (_plane_strain)
    factor = 1.0 + ElasticityTensorTools::getIsotropicPoissonsRatio(_Cijkl[0]);
 
  for (auto es : _eigenstrains)
    _mechanical_stretch[0] -= 0.5 * factor * ((*es)[0](2, 2) + (*es)[1](2, 2));
 
  _mechanical_stretch[1] = _mechanical_stretch[0];
}

Referenced by computeProperties().

◆ computeConventionalQpDeformationGradient()

void ComputeStrainBaseNOSPD::computeConventionalQpDeformationGradient ( )

protectedvirtual

Function to compute conventional nonlocal deformation gradient.

Definition at line 113 of file ComputeStrainBaseNOSPD.C.

{
  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(_qp));
  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(_qp));
 
  // calculate the shape tensor and prepare the deformation gradient tensor
  Real vol_nb, weight_nb;
  RealGradient origin_vec_nb, current_vec_nb;
 
  for (unsigned int nb = 0; nb < neighbors.size(); ++nb)
    if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[nb])) > 0.5)
    {
      vol_nb = _pdmesh.getNodeVolume(neighbors[nb]);
      origin_vec_nb =
          _pdmesh.getNodeCoord(neighbors[nb]) - _pdmesh.getNodeCoord(_current_elem->node_id(_qp));
 
      for (unsigned int k = 0; k < _dim; ++k)
        current_vec_nb(k) = origin_vec_nb(k) +
                            _disp_var[k]->getNodalValue(*_pdmesh.nodePtr(neighbors[nb])) -
                            _disp_var[k]->getNodalValue(*_current_elem->node_ptr(_qp));
 
      weight_nb = _horizon_radius[_qp] / origin_vec_nb.norm();
      for (unsigned int k = 0; k < _dim; ++k)
      {
        for (unsigned int l = 0; l < _dim; ++l)
        {
          _shape2[_qp](k, l) += weight_nb * origin_vec_nb(k) * origin_vec_nb(l) * vol_nb;
          _deformation_gradient[_qp](k, l) +=
              weight_nb * current_vec_nb(k) * origin_vec_nb(l) * vol_nb;
        }
        // calculate derivatives of deformation_gradient w.r.t displacements of node nb
        _ddgraddu[_qp](0, k) += -weight_nb * origin_vec_nb(k) * vol_nb;
        _ddgraddv[_qp](1, k) += -weight_nb * origin_vec_nb(k) * vol_nb;
        if (_dim == 3)
          _ddgraddw[_qp](2, k) += -weight_nb * origin_vec_nb(k) * vol_nb;
      }
    }
 
  // finalize the deformation gradient tensor
  if (_shape2[_qp].det() == 0.)
    mooseError("Singular shape tensor is detected in ComputeStrainBaseNOSPD! Use "
               "SingularShapeTensorEliminatorUserObjectPD to avoid singular shape tensor!");
 
  _deformation_gradient[_qp] *= _shape2[_qp].inverse();
  _ddgraddu[_qp] *= _shape2[_qp].inverse();
  _ddgraddv[_qp] *= _shape2[_qp].inverse();
  _ddgraddw[_qp] *= _shape2[_qp].inverse();
}

Referenced by computeBondHorizonQpDeformationGradient(), and computeQpDeformationGradient().

◆ computeProperties()

void ComputeStrainBaseNOSPD::computeProperties ( )

overrideprotectedvirtual

Definition at line 230 of file ComputeStrainBaseNOSPD.C.

{
  setupMeshRelatedData();     // function from base class
  computeBondCurrentLength(); // current length of a bond from base class
  computeBondStretch();       // stretch of a bond
 
  for (_qp = 0; _qp < _nnodes; ++_qp)
    computeQpStrain();
}

◆ computeQpDeformationGradient()

void ComputeStrainBaseNOSPD::computeQpDeformationGradient ( )

protectedvirtual

Function to compute deformation gradient for peridynamic correspondence model.

Definition at line 75 of file ComputeStrainBaseNOSPD.C.

{
  _shape2[_qp].zero();
  _deformation_gradient[_qp].zero();
  _ddgraddu[_qp].zero();
  _ddgraddv[_qp].zero();
  _ddgraddw[_qp].zero();
  _multi[_qp] = 0.0;
  _sf_coeff[_qp] = 0.0;
 
  if (_dim == 2)
    _shape2[_qp](2, 2) = _deformation_gradient[_qp](2, 2) = 1.0;
 
  if (_bond_status_var->getElementalValue(_current_elem) > 0.5)
  {
    if (_stabilization == "FORCE")
    {
      computeConventionalQpDeformationGradient();
 
      _sf_coeff[_qp] = ElasticityTensorTools::getIsotropicYoungsModulus(_Cijkl[_qp]) *
                       _pdmesh.getNodeAverageSpacing(_current_elem->node_id(_qp)) *
                       _horizon_radius[_qp] / _origin_vec.norm();
      _multi[_qp] = _horizon_radius[_qp] / _origin_vec.norm() * _node_vol[0] * _node_vol[1];
    }
    else if (_stabilization == "BOND_HORIZON_I" || _stabilization == "BOND_HORIZON_II")
      computeBondHorizonQpDeformationGradient();
    else
      paramError("stabilization",
                 "Unknown stabilization scheme for peridynamic correspondence model!");
  }
  else
  {
    _shape2[_qp].setToIdentity();
    _deformation_gradient[_qp].setToIdentity();
  }
}

Referenced by ComputeFiniteStrainNOSPD::computeQpStrain(), and ComputeSmallStrainNOSPD::computeQpStrain().

◆ computeQpStrain()

virtual void ComputeStrainBaseNOSPD::computeQpStrain ( )

protectedpure virtual

Function to compute strain tensors.

Implemented in ComputeFiniteStrainNOSPD, and ComputeSmallStrainNOSPD.

Referenced by computeProperties().

◆ initQpStatefulProperties()

void ComputeStrainBaseNOSPD::initQpStatefulProperties ( )

overridevirtual

Definition at line 62 of file ComputeStrainBaseNOSPD.C.

{
  _mechanical_strain[_qp].zero();
  _total_strain[_qp].zero();
  _deformation_gradient[_qp].setToIdentity();
 
  if (_qrule->n_points() < 2) // Gauss_Lobatto: order = 2n-3 (n is number of qp)
    mooseError(
        "NOSPD models require Gauss_Lobatto rule and a minimum of 2 quadrature points, i.e., "
        "order of FIRST");
}

◆ validParams()

InputParameters ComputeStrainBaseNOSPD::validParams ( )

static

Definition at line 16 of file ComputeStrainBaseNOSPD.C.

{
  InputParameters params = MechanicsMaterialBasePD::validParams();
  params.addClassDescription(
      "Base material strain class for the peridynamic correspondence models");
 
  MooseEnum stabilization_option("FORCE BOND_HORIZON_I BOND_HORIZON_II");
  params.addRequiredParam<MooseEnum>(
      "stabilization",
      stabilization_option,
      "Stabilization techniques used for the peridynamic correspondence model");
  params.addParam<bool>("plane_strain",
                        false,
                        "Plane strain problem or not, this will affect the mechanical stretch "
                        "calculation for problem with eigenstrains");
  params.addParam<std::vector<MaterialPropertyName>>(
      "eigenstrain_names", {}, "List of eigenstrains to be applied in this strain calculation");
 
  return params;
}

Referenced by ComputeFiniteStrainNOSPD::validParams(), and ComputeSmallStrainNOSPD::validParams().

Member Data Documentation

◆ _Cijkl

const MaterialProperty<RankFourTensor>& ComputeStrainBaseNOSPD::_Cijkl

protected

Material properties to fetch.

Definition at line 59 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondStretch(), and computeQpDeformationGradient().

◆ _ddgraddu

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_ddgraddu

protected

Definition at line 68 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondHorizonQpDeformationGradient(), computeConventionalQpDeformationGradient(), and computeQpDeformationGradient().

◆ _ddgraddv

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_ddgraddv

protected

Definition at line 69 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondHorizonQpDeformationGradient(), computeConventionalQpDeformationGradient(), and computeQpDeformationGradient().

◆ _ddgraddw

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_ddgraddw

protected

Definition at line 70 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondHorizonQpDeformationGradient(), computeConventionalQpDeformationGradient(), and computeQpDeformationGradient().

◆ _deformation_gradient

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_deformation_gradient

protected

Definition at line 66 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondHorizonQpDeformationGradient(), computeConventionalQpDeformationGradient(), computeQpDeformationGradient(), ComputeFiniteStrainNOSPD::computeQpFhat(), ComputePlaneFiniteStrainNOSPD::computeQpFhat(), ComputeSmallStrainNOSPD::computeQpTotalStrain(), ComputePlaneSmallStrainNOSPD::computeQpTotalStrain(), and initQpStatefulProperties().

◆ _eigenstrain_names

std::vector<MaterialPropertyName> ComputeStrainBaseNOSPD::_eigenstrain_names

protected

Definition at line 60 of file ComputeStrainBaseNOSPD.h.

Referenced by ComputeFiniteStrainNOSPD::ComputeFiniteStrainNOSPD(), and ComputeStrainBaseNOSPD().

◆ _eigenstrains

std::vector<const MaterialProperty<RankTwoTensor> *> ComputeStrainBaseNOSPD::_eigenstrains

protected

Definition at line 61 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondStretch(), ComputeSmallStrainNOSPD::computeQpStrain(), ComputeStrainBaseNOSPD(), and ComputeFiniteStrainNOSPD::subtractEigenstrainIncrementFromStrain().

◆ _mechanical_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_mechanical_strain

protected

Definition at line 73 of file ComputeStrainBaseNOSPD.h.

Referenced by ComputeFiniteStrainNOSPD::computeQpStrain(), ComputeSmallStrainNOSPD::computeQpStrain(), and initQpStatefulProperties().

◆ _multi

MaterialProperty<Real>& ComputeStrainBaseNOSPD::_multi

protected

Definition at line 75 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondHorizonQpDeformationGradient(), and computeQpDeformationGradient().

◆ _plane_strain

const bool ComputeStrainBaseNOSPD::_plane_strain

protected

Plane strain problem or not, this is only used for mechanical stretch calculation.

Definition at line 56 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondStretch().

◆ _sf_coeff

MaterialProperty<Real>& ComputeStrainBaseNOSPD::_sf_coeff

protected

fictitious force coefficient for force stabilized model

Definition at line 78 of file ComputeStrainBaseNOSPD.h.

Referenced by computeQpDeformationGradient().

◆ _shape2

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_shape2

protected

Material properties to store.

Definition at line 65 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondHorizonQpDeformationGradient(), computeConventionalQpDeformationGradient(), and computeQpDeformationGradient().

◆ _stabilization

const MooseEnum ComputeStrainBaseNOSPD::_stabilization

protected

Option of stabilization scheme for correspondence material model: FORCE, BOND_HORIZON_I or BOND_HORIZON_II.

Definition at line 53 of file ComputeStrainBaseNOSPD.h.

Referenced by computeBondHorizonQpDeformationGradient(), and computeQpDeformationGradient().

◆ _total_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_total_strain

protected

Definition at line 72 of file ComputeStrainBaseNOSPD.h.

Referenced by ComputeFiniteStrainNOSPD::computeQpStrain(), ComputeSmallStrainNOSPD::computeQpStrain(), ComputeSmallStrainNOSPD::computeQpTotalStrain(), ComputePlaneSmallStrainNOSPD::computeQpTotalStrain(), and initQpStatefulProperties().

The documentation for this class was generated from the following files:

peridynamics/include/materials/ComputeStrainBaseNOSPD.h
peridynamics/src/materials/ComputeStrainBaseNOSPD.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ComputeStrainBaseNOSPD()

Member Function Documentation

◆ computeBondHorizonQpDeformationGradient()

◆ computeBondStretch()

◆ computeConventionalQpDeformationGradient()

◆ computeProperties()

◆ computeQpDeformationGradient()

◆ computeQpStrain()

◆ initQpStatefulProperties()

◆ validParams()

Member Data Documentation

◆ _Cijkl

◆ _ddgraddu

◆ _ddgraddv

◆ _ddgraddw

◆ _deformation_gradient

◆ _eigenstrain_names

◆ _eigenstrains

◆ _mechanical_strain

◆ _multi

◆ _plane_strain

◆ _sf_coeff

◆ _shape2

◆ _stabilization

◆ _total_strain