Material class for peridynamic correspondence model for finite strain. More...

#include <ComputeFiniteStrainNOSPD.h>

Inheritance diagram for ComputeFiniteStrainNOSPD:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ComputeFiniteStrainNOSPD (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 ()

static MooseEnum	decompositionType ()

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

virtual void	computeQpFhat ()

void	computeQpStrainRotationIncrements (RankTwoTensor &e, RankTwoTensor &r)
	Function to compute strain and rotational increments. More...

void	subtractEigenstrainIncrementFromStrain (RankTwoTensor &strain)
	Function to compute the mechanical strain tensor by subtracting thermal strain from the total strain. 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
std::vector< RankTwoTensor >	_Fhat
	'Incremental' deformation gradient More...

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


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

MaterialProperty< RankTwoTensor > &	_strain_increment

MaterialProperty< RankTwoTensor > &	_rotation_increment


const MaterialProperty< RankTwoTensor > &	_deformation_gradient_old
	Material properties to fetch. More...

const MaterialProperty< RankTwoTensor > &	_mechanical_strain_old

const MaterialProperty< RankTwoTensor > &	_total_strain_old

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


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

enum	DecompMethod { DecompMethod::TaylorExpansion, DecompMethod::EigenSolution }
	Method to decompose into rotation increment and strain increment. More...

const DecompMethod	_decomposition_method

Detailed Description

Material class for peridynamic correspondence model for finite strain.

Definition at line 18 of file ComputeFiniteStrainNOSPD.h.

Member Enumeration Documentation

◆ DecompMethod

enum ComputeFiniteStrainNOSPD::DecompMethod

strongprivate

Method to decompose into rotation increment and strain increment.

Enumerator
TaylorExpansion
EigenSolution

Definition at line 60 of file ComputeFiniteStrainNOSPD.h.

   {
     TaylorExpansion,
     EigenSolution
   };

Constructor & Destructor Documentation

◆ ComputeFiniteStrainNOSPD()

ComputeFiniteStrainNOSPD::ComputeFiniteStrainNOSPD ( const InputParameters & parameters )

Definition at line 37 of file ComputeFiniteStrainNOSPD.C.

   : ComputeStrainBaseNOSPD(parameters),
     _strain_rate(declareProperty<RankTwoTensor>("strain_rate")),
     _strain_increment(declareProperty<RankTwoTensor>("strain_increment")),
     _rotation_increment(declareProperty<RankTwoTensor>("rotation_increment")),
     _deformation_gradient_old(getMaterialPropertyOld<RankTwoTensor>("deformation_gradient")),
     _mechanical_strain_old(getMaterialPropertyOld<RankTwoTensor>("mechanical_strain")),
     _total_strain_old(getMaterialPropertyOld<RankTwoTensor>("total_strain")),
     _eigenstrains_old(_eigenstrain_names.size()),
     _Fhat(2),
     _decomposition_method(getParam<MooseEnum>("decomposition_method").getEnum<DecompMethod>())
 {
   for (unsigned int i = 0; i < _eigenstrains_old.size(); ++i)
     _eigenstrains_old[i] = &getMaterialPropertyOld<RankTwoTensor>(_eigenstrain_names[i]);
 }

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

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

◆ computeQpFhat()

void ComputeFiniteStrainNOSPD::computeQpFhat ( )

protectedvirtual

Reimplemented in ComputePlaneFiniteStrainNOSPD.

Definition at line 97 of file ComputeFiniteStrainNOSPD.C.

Referenced by computeQpStrain().

 {
   // Incremental deformation gradient of current step w.r.t previous step:
   // _Fhat = deformation_gradient * inv(deformation_gradient_old)
   _Fhat[_qp] = _deformation_gradient[_qp] * _deformation_gradient_old[_qp].inverse();
 }

◆ computeQpStrain()

void ComputeFiniteStrainNOSPD::computeQpStrain ( )

overrideprotectedvirtual

Function to compute strain tensors.

Implements ComputeStrainBaseNOSPD.

Definition at line 54 of file ComputeFiniteStrainNOSPD.C.

 {
   computeQpDeformationGradient();
 
   computeQpFhat();
 
   RankTwoTensor total_strain_increment;
 
   // Two ways to calculate these increments: TaylorExpansion(default) or EigenSolution
   computeQpStrainRotationIncrements(total_strain_increment, _rotation_increment[_qp]);
 
   _strain_increment[_qp] = total_strain_increment;
 
   // Remove the eigenstrain increment
   subtractEigenstrainIncrementFromStrain(_strain_increment[_qp]);
 
   if (_dt > 0)
     _strain_rate[_qp] = _strain_increment[_qp] / _dt;
   else
     _strain_rate[_qp].zero();
 
   // Update strain in intermediate configuration
   _mechanical_strain[_qp] = _mechanical_strain_old[_qp] + _strain_increment[_qp];
   _total_strain[_qp] = _total_strain_old[_qp] + total_strain_increment;
 
   // Rotate strain to current configuration
   _mechanical_strain[_qp] =
       _rotation_increment[_qp] * _mechanical_strain[_qp] * _rotation_increment[_qp].transpose();
   _total_strain[_qp] =
       _rotation_increment[_qp] * _total_strain[_qp] * _rotation_increment[_qp].transpose();
 
   // zero out all strain measures for broken bond
   if (_bond_status_var->getElementalValue(_current_elem) < 0.5)
   {
     _strain_rate[_qp].zero();
     _strain_increment[_qp].zero();
     _rotation_increment[_qp].zero();
     _mechanical_strain[_qp].zero();
     _total_strain[_qp].zero();
   }
 }

◆ computeQpStrainRotationIncrements()

void ComputeFiniteStrainNOSPD::computeQpStrainRotationIncrements	(	RankTwoTensor &	e,
		RankTwoTensor &	r
	)

protected

Function to compute strain and rotational increments.

Definition at line 105 of file ComputeFiniteStrainNOSPD.C.

Referenced by computeQpStrain().

 {
   switch (_decomposition_method)
   {
     case DecompMethod::TaylorExpansion:
     {
       // inverse of _Fhat
       RankTwoTensor invFhat(_Fhat[_qp].inverse());
 
       // A = I - _Fhat^-1
       RankTwoTensor A(RankTwoTensor::initIdentity);
       A -= invFhat;
 
       // Cinv - I = A A^T - A - A^T;
       RankTwoTensor Cinv_I = A * A.transpose() - A - A.transpose();
 
       // strain rate D from Taylor expansion, Chat = (-1/2(Chat^-1 - I) + 1/4*(Chat^-1 - I)^2 + ...
       total_strain_increment = -Cinv_I * 0.5 + Cinv_I * Cinv_I * 0.25;
 
       const Real a[3] = {invFhat(1, 2) - invFhat(2, 1),
                          invFhat(2, 0) - invFhat(0, 2),
                          invFhat(0, 1) - invFhat(1, 0)};
 
       Real q = (a[0] * a[0] + a[1] * a[1] + a[2] * a[2]) / 4.0;
       Real trFhatinv_1 = invFhat.trace() - 1.0;
       const Real p = trFhatinv_1 * trFhatinv_1 / 4.0;
 
       // cos theta_a
       const Real C1_squared = p +
                               3.0 * Utility::pow<2>(p) * (1.0 - (p + q)) / Utility::pow<2>(p + q) -
                               2.0 * Utility::pow<3>(p) * (1.0 - (p + q)) / Utility::pow<3>(p + q);
       if (C1_squared <= 0.0)
         mooseException(
             "Cannot take square root of a number less than or equal to zero in the calculation of "
             "C1 for the Rashid approximation for the rotation tensor.");
 
       const Real C1 = std::sqrt(C1_squared);
 
       Real C2;
       if (q > 0.01)
         // (1-cos theta_a)/4q
         C2 = (1.0 - C1) / (4.0 * q);
       else
         // alternate form for small q
         C2 = 0.125 + q * 0.03125 * (Utility::pow<2>(p) - 12.0 * (p - 1.0)) / Utility::pow<2>(p) +
              Utility::pow<2>(q) * (p - 2.0) * (Utility::pow<2>(p) - 10.0 * p + 32.0) /
                  Utility::pow<3>(p) +
              Utility::pow<3>(q) *
                  (1104.0 - 992.0 * p + 376.0 * Utility::pow<2>(p) - 72.0 * Utility::pow<3>(p) +
                   5.0 * Utility::pow<4>(p)) /
                  (512.0 * Utility::pow<4>(p));
 
       const Real C3_test =
           (p * q * (3.0 - q) + Utility::pow<3>(p) + Utility::pow<2>(q)) / Utility::pow<3>(p + q);
 
       if (C3_test <= 0.0)
         mooseException(
             "Cannot take square root of a number less than or equal to zero in the calculation of "
             "C3_test for the Rashid approximation for the rotation tensor.");
 
       const Real C3 = 0.5 * std::sqrt(C3_test); // sin theta_a/(2 sqrt(q))
 
       // Calculate incremental rotation. Note that this value is the transpose of that from Rashid,
       // 93, so we transpose it before storing
       RankTwoTensor R_incr;
       R_incr.addIa(C1);
       for (unsigned int i = 0; i < 3; ++i)
         for (unsigned int j = 0; j < 3; ++j)
           R_incr(i, j) += C2 * a[i] * a[j];
 
       R_incr(0, 1) += C3 * a[2];
       R_incr(0, 2) -= C3 * a[1];
       R_incr(1, 0) -= C3 * a[2];
       R_incr(1, 2) += C3 * a[0];
       R_incr(2, 0) += C3 * a[1];
       R_incr(2, 1) -= C3 * a[0];
 
       rotation_increment = R_incr.transpose();
       break;
     }
     case DecompMethod::EigenSolution:
     {
       FactorizedRankTwoTensor Chat = RankTwoTensor::transposeTimes(_Fhat[_qp]);
       FactorizedRankTwoTensor Uhat = MathUtils::sqrt(Chat);
       rotation_increment = _Fhat[_qp] * Uhat.inverse().get();
       total_strain_increment = MathUtils::log(Uhat).get();
       break;
     }
 
     default:
       mooseError("ComputeFiniteStrainNOSPD Error: Invalid decomposition type! Please specify : "
                  "TaylorExpansion or "
                  "EigenSolution.");
   }
 }

◆ decompositionType()

MooseEnum ComputeFiniteStrainNOSPD::decompositionType ( )

static

Definition at line 17 of file ComputeFiniteStrainNOSPD.C.

Referenced by validParams().

 {
   return MooseEnum("TaylorExpansion EigenSolution", "TaylorExpansion");
 }

◆ 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");
 }

◆ subtractEigenstrainIncrementFromStrain()

void ComputeFiniteStrainNOSPD::subtractEigenstrainIncrementFromStrain ( RankTwoTensor & strain )

protected

Function to compute the mechanical strain tensor by subtracting thermal strain from the total strain.

Definition at line 203 of file ComputeFiniteStrainNOSPD.C.

Referenced by computeQpStrain().

 {
   for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
   {
     strain -= (*_eigenstrains[i])[_qp];
     strain += (*_eigenstrains_old[i])[_qp];
   }
 }

◆ validParams()

InputParameters ComputeFiniteStrainNOSPD::validParams ( )

static

Definition at line 23 of file ComputeFiniteStrainNOSPD.C.

Referenced by ComputePlaneFiniteStrainNOSPD::validParams().

 {
   InputParameters params = ComputeStrainBaseNOSPD::validParams();
   params.addClassDescription(
       "Class for computing nodal quantities for residual and jacobian calculation "
       "for peridynamic correspondence models under finite strain assumptions");
 
   params.addParam<MooseEnum>("decomposition_method",
                              ComputeFiniteStrainNOSPD::decompositionType(),
                              "Methods to calculate the strain and rotation increments");
 
   return params;
 }

Member Data Documentation

◆ _Cijkl

const MaterialProperty<RankFourTensor>& ComputeStrainBaseNOSPD::_Cijkl

protectedinherited

Material properties to fetch.

Definition at line 59 of file ComputeStrainBaseNOSPD.h.

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

◆ _ddgraddu

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_ddgraddu

protectedinherited

Definition at line 68 of file ComputeStrainBaseNOSPD.h.

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

◆ _ddgraddv

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_ddgraddv

protectedinherited

Definition at line 69 of file ComputeStrainBaseNOSPD.h.

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

◆ _ddgraddw

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_ddgraddw

protectedinherited

Definition at line 70 of file ComputeStrainBaseNOSPD.h.

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

◆ _decomposition_method

const DecompMethod ComputeFiniteStrainNOSPD::_decomposition_method

private

Definition at line 65 of file ComputeFiniteStrainNOSPD.h.

Referenced by computeQpStrainRotationIncrements().

◆ _deformation_gradient

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_deformation_gradient

protectedinherited

Definition at line 66 of file ComputeStrainBaseNOSPD.h.

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

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainNOSPD::_deformation_gradient_old

protected

Material properties to fetch.

Definition at line 49 of file ComputeFiniteStrainNOSPD.h.

Referenced by ComputePlaneFiniteStrainNOSPD::computeQpFhat(), and computeQpFhat().

◆ _eigenstrain_names

std::vector<MaterialPropertyName> ComputeStrainBaseNOSPD::_eigenstrain_names

protectedinherited

Definition at line 60 of file ComputeStrainBaseNOSPD.h.

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

◆ _eigenstrains

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

protectedinherited

Definition at line 61 of file ComputeStrainBaseNOSPD.h.

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

◆ _eigenstrains_old

std::vector<const MaterialProperty<RankTwoTensor> *> ComputeFiniteStrainNOSPD::_eigenstrains_old

protected

Definition at line 52 of file ComputeFiniteStrainNOSPD.h.

Referenced by ComputeFiniteStrainNOSPD(), and subtractEigenstrainIncrementFromStrain().

◆ _Fhat

std::vector<RankTwoTensor> ComputeFiniteStrainNOSPD::_Fhat

protected

'Incremental' deformation gradient

Definition at line 56 of file ComputeFiniteStrainNOSPD.h.

Referenced by ComputePlaneFiniteStrainNOSPD::computeQpFhat(), computeQpFhat(), and computeQpStrainRotationIncrements().

◆ _mechanical_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_mechanical_strain

protectedinherited

Definition at line 73 of file ComputeStrainBaseNOSPD.h.

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

◆ _mechanical_strain_old

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainNOSPD::_mechanical_strain_old

protected

Definition at line 50 of file ComputeFiniteStrainNOSPD.h.

Referenced by computeQpStrain().

◆ _multi

MaterialProperty<Real>& ComputeStrainBaseNOSPD::_multi

protectedinherited

Definition at line 75 of file ComputeStrainBaseNOSPD.h.

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

◆ _plane_strain

const bool ComputeStrainBaseNOSPD::_plane_strain

protectedinherited

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

Definition at line 56 of file ComputeStrainBaseNOSPD.h.

Referenced by ComputeStrainBaseNOSPD::computeBondStretch().

◆ _rotation_increment

MaterialProperty<RankTwoTensor>& ComputeFiniteStrainNOSPD::_rotation_increment

protected

Definition at line 45 of file ComputeFiniteStrainNOSPD.h.

Referenced by computeQpStrain().

◆ _sf_coeff

MaterialProperty<Real>& ComputeStrainBaseNOSPD::_sf_coeff

protectedinherited

fictitious force coefficient for force stabilized model

Definition at line 78 of file ComputeStrainBaseNOSPD.h.

Referenced by ComputeStrainBaseNOSPD::computeQpDeformationGradient().

◆ _shape2

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_shape2

protectedinherited

Material properties to store.

Definition at line 65 of file ComputeStrainBaseNOSPD.h.

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

◆ _stabilization

const MooseEnum ComputeStrainBaseNOSPD::_stabilization

protectedinherited

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 ComputeStrainBaseNOSPD::computeBondHorizonQpDeformationGradient(), and ComputeStrainBaseNOSPD::computeQpDeformationGradient().

◆ _strain_increment

MaterialProperty<RankTwoTensor>& ComputeFiniteStrainNOSPD::_strain_increment

protected

Definition at line 44 of file ComputeFiniteStrainNOSPD.h.

Referenced by computeQpStrain().

◆ _strain_rate

MaterialProperty<RankTwoTensor>& ComputeFiniteStrainNOSPD::_strain_rate

protected

Material properties to store.

Definition at line 43 of file ComputeFiniteStrainNOSPD.h.

Referenced by computeQpStrain().

◆ _total_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBaseNOSPD::_total_strain

protectedinherited

Definition at line 72 of file ComputeStrainBaseNOSPD.h.

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

◆ _total_strain_old

const MaterialProperty<RankTwoTensor>& ComputeFiniteStrainNOSPD::_total_strain_old

protected

Definition at line 51 of file ComputeFiniteStrainNOSPD.h.

Referenced by computeQpStrain().

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

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

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ DecompMethod

Constructor & Destructor Documentation

◆ ComputeFiniteStrainNOSPD()

Member Function Documentation

◆ computeBondHorizonQpDeformationGradient()

◆ computeBondStretch()

◆ computeConventionalQpDeformationGradient()

◆ computeProperties()

◆ computeQpDeformationGradient()

◆ computeQpFhat()

◆ computeQpStrain()

◆ computeQpStrainRotationIncrements()

◆ decompositionType()

◆ initQpStatefulProperties()

◆ subtractEigenstrainIncrementFromStrain()

◆ validParams()

Member Data Documentation

◆ _Cijkl

◆ _ddgraddu

◆ _ddgraddv

◆ _ddgraddw

◆ _decomposition_method

◆ _deformation_gradient

◆ _deformation_gradient_old

◆ _eigenstrain_names

◆ _eigenstrains

◆ _eigenstrains_old

◆ _Fhat

◆ _mechanical_strain

◆ _mechanical_strain_old

◆ _multi

◆ _plane_strain

◆ _rotation_increment

◆ _sf_coeff

◆ _shape2

◆ _stabilization

◆ _strain_increment

◆ _strain_rate

◆ _total_strain

◆ _total_strain_old