Material class for 2D correspondence material model for small strain: plane strain, generalized plane strain, weak plane stress. More...

#include <ComputePlaneSmallStrainNOSPD.h>

Inheritance diagram for ComputePlaneSmallStrainNOSPD:

[legend]

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ComputePlaneSmallStrainNOSPD (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	computeQpTotalStrain () override
	Function to compute the total strain tensor for small strain case. More...

Real	computeOutOfPlaneStrain ()
	Function to compute out-of-plane component of strain tensor for generalized plane strain and weak plane stress. More...

virtual void	computeQpStrain () override
	Function to compute strain tensors. More...

virtual void	computeProperties () override

virtual void	computeBondStretch () override

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

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

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

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

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


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

std::vector< MaterialPropertyName >	_eigenstrain_names

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


MaterialProperty< RankTwoTensor > &	_shape2
	Material properties to store. More...

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 More...

Private Attributes

const bool	_scalar_out_of_plane_strain_coupled
	Scalar out-of-plane strain for generalized plane strain. More...

const VariableValue &	_scalar_out_of_plane_strain


const bool	_out_of_plane_strain_coupled
	Out-of-plane strain for weak plane stress. More...

const VariableValue &	_out_of_plane_strain

Detailed Description

Material class for 2D correspondence material model for small strain: plane strain, generalized plane strain, weak plane stress.

Definition at line 19 of file ComputePlaneSmallStrainNOSPD.h.

Constructor & Destructor Documentation

◆ ComputePlaneSmallStrainNOSPD()

ComputePlaneSmallStrainNOSPD::ComputePlaneSmallStrainNOSPD ( const InputParameters & parameters )

Definition at line 30 of file ComputePlaneSmallStrainNOSPD.C.

   : ComputeSmallStrainNOSPD(parameters),
     _scalar_out_of_plane_strain_coupled(isCoupledScalar("scalar_out_of_plane_strain")),
     _scalar_out_of_plane_strain(_scalar_out_of_plane_strain_coupled
                                     ? coupledScalarValue("scalar_out_of_plane_strain")
                                     : _zero),
     _out_of_plane_strain_coupled(isCoupled("out_of_plane_strain")),
     _out_of_plane_strain(_out_of_plane_strain_coupled ? coupledValue("out_of_plane_strain") : _zero)
 {
 }

Member Function Documentation

◆ computeBondHorizonQpDeformationGradient()

void ComputeStrainBaseNOSPD::computeBondHorizonQpDeformationGradient ( )

protectedvirtualinherited

Function to compute bond-associated horizon based deformation gradient.

Definition at line 163 of file ComputeStrainBaseNOSPD.C.

Referenced by ComputeStrainBaseNOSPD::computeQpDeformationGradient().

 {
   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;
 }

◆ computeBondStretch()

void ComputeStrainBaseNOSPD::computeBondStretch ( )

overrideprotectedvirtualinherited

Definition at line 241 of file ComputeStrainBaseNOSPD.C.

Referenced by ComputeStrainBaseNOSPD::computeProperties().

 {
   _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];
 }

◆ computeConventionalQpDeformationGradient()

void ComputeStrainBaseNOSPD::computeConventionalQpDeformationGradient ( )

protectedvirtualinherited

Function to compute conventional nonlocal deformation gradient.

Definition at line 113 of file ComputeStrainBaseNOSPD.C.

Referenced by ComputeStrainBaseNOSPD::computeBondHorizonQpDeformationGradient(), and ComputeStrainBaseNOSPD::computeQpDeformationGradient().

 {
   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();
 }

◆ computeOutOfPlaneStrain()

Real ComputePlaneSmallStrainNOSPD::computeOutOfPlaneStrain ( )

protected

Function to compute out-of-plane component of strain tensor for generalized plane strain and weak plane stress.

Returns: The value of out-of-plane strain

Definition at line 51 of file ComputePlaneSmallStrainNOSPD.C.

Referenced by computeQpTotalStrain().

 {
   if (_scalar_out_of_plane_strain_coupled)
     return _scalar_out_of_plane_strain[0];
   else
     return _out_of_plane_strain[_qp];
 }

◆ computeProperties()

void ComputeStrainBaseNOSPD::computeProperties ( )

overrideprotectedvirtualinherited

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 ( )

protectedvirtualinherited

Function to compute deformation gradient for peridynamic correspondence model.

Definition at line 75 of file ComputeStrainBaseNOSPD.C.

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

 {
   _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();
   }
 }

◆ computeQpStrain()

void ComputeSmallStrainNOSPD::computeQpStrain ( )

overrideprotectedvirtualinherited

Function to compute strain tensors.

Implements ComputeStrainBaseNOSPD.

Definition at line 31 of file ComputeSmallStrainNOSPD.C.

 {
   computeQpDeformationGradient();
 
   computeQpTotalStrain();
 
   _mechanical_strain[_qp] = _total_strain[_qp];
 
   // Subtract Eigen strains
   for (auto es : _eigenstrains)
     _mechanical_strain[_qp] -= (*es)[_qp];
 
   // zero out all strain measures for broken bond
   if (_bond_status_var->getElementalValue(_current_elem) < 0.5)
   {
     _mechanical_strain[_qp].zero();
     _total_strain[_qp].zero();
   }
 }

◆ computeQpTotalStrain()

void ComputePlaneSmallStrainNOSPD::computeQpTotalStrain ( )

overrideprotectedvirtual

Function to compute the total strain tensor for small strain case.

Reimplemented from ComputeSmallStrainNOSPD.

Definition at line 42 of file ComputePlaneSmallStrainNOSPD.C.

 {
   // the green-lagrange strain tensor
   _total_strain[_qp] = 0.5 * (_deformation_gradient[_qp].transpose() + _deformation_gradient[_qp]) -
                        RankTwoTensor(RankTwoTensor::initIdentity);
   _total_strain[_qp](2, 2) = computeOutOfPlaneStrain();
 }

◆ initQpStatefulProperties()

void ComputeStrainBaseNOSPD::initQpStatefulProperties ( )

overridevirtualinherited

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 ComputePlaneSmallStrainNOSPD::validParams ( )

static

Definition at line 15 of file ComputePlaneSmallStrainNOSPD.C.

 {
   InputParameters params = ComputeSmallStrainNOSPD::validParams();
   params.addClassDescription(
       "Class for computing nodal quantities for residual and jacobian calculation "
       "for peridynamic correspondence models under planar small strain assumptions");
 
   params.addCoupledVar("scalar_out_of_plane_strain",
                        "Scalar out-of-plane strain variable for generalized plane strain");
   params.addCoupledVar("out_of_plane_strain",
                        "Nonlinear out-of-plane strain variable for plane stress condition");
 
   return params;
 }