ComputeDamageStress computes the stress for a damaged elasticity material. More...

#include <ComputeDamageStress.h>

Inheritance diagram for ComputeDamageStress:

Public Member Functions
	ComputeDamageStress (const InputParameters &parameters)

void	initialSetup () override

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	computeQpStress () override
	Compute the stress and store it in the _stress material property for the current quadrature point. More...

virtual void	initQpStatefulProperties () override

virtual void	computeQpProperties () override

bool	hasGuaranteedMaterialProperty (const MaterialPropertyName &prop, Guarantee guarantee)

Protected Attributes
MaterialProperty< Real > &	_matl_timestep_limit
	Property that stores the time step limit. More...

DamageBase *	_damage_model
	Pointer to the damage model. More...

const std::string	_elasticity_tensor_name
	Name of the elasticity tensor material property. More...

const MaterialProperty< RankFourTensor > &	_elasticity_tensor
	Elasticity tensor material property. More...

const MaterialProperty< RankTwoTensor > &	_strain_increment
	Strain increment material property. More...

const MaterialProperty< RankTwoTensor > &	_rotation_increment
	Rotation increment material property. More...

const MaterialProperty< RankTwoTensor > &	_stress_old
	Old state of the stress tensor material property. More...

const MaterialProperty< RankTwoTensor > &	_elastic_strain_old
	The old elastic strain is used to calculate the old stress in the case of variable elasticity tensors. More...

const std::string	_base_name
	Base name prepended to all material property names to allow for multi-material systems. More...

const MaterialProperty< RankTwoTensor > &	_mechanical_strain
	Mechanical strain material property. More...

MaterialProperty< RankTwoTensor > &	_stress
	Stress material property. More...

MaterialProperty< RankTwoTensor > &	_elastic_strain
	Elastic strain material property. More...

const MaterialProperty< RankTwoTensor > &	_extra_stress
	Extra stress tensor. More...

std::vector< const Function * >	_initial_stress_fcn
	initial stress components More...

MaterialProperty< RankFourTensor > &	_Jacobian_mult
	derivative of stress w.r.t. strain (_dstress_dstrain) More...

Private Attributes
const InputParameters &	_gc_params
	Parameters of the object with this interface. More...

FEProblemBase *const	_gc_feproblem
	Reference to the FEProblemBase class. More...

BlockRestrictable *const	_gc_block_restrict
	Access block restrictions of the object with this interface. More...

Detailed Description

ComputeDamageStress computes the stress for a damaged elasticity material.

This model must be used in conjunction with a damage model (derived from DamageBase)

Definition at line 24 of file ComputeDamageStress.h.

Constructor & Destructor Documentation

◆ ComputeDamageStress()

ComputeDamageStress::ComputeDamageStress ( const InputParameters & parameters )

Definition at line 27 of file ComputeDamageStress.C.

   : ComputeFiniteStrainElasticStress(parameters),
     _matl_timestep_limit(declareProperty<Real>("matl_timestep_limit")),
     _damage_model(nullptr)
 {
 }

Member Function Documentation

◆ computeQpProperties()

void ComputeStressBase::computeQpProperties ( )

overrideprotectedvirtualinherited

Definition at line 50 of file ComputeStressBase.C.

 {
   computeQpStress();
  
   // Add in extra stress
   _stress[_qp] += _extra_stress[_qp];
 }

◆ computeQpStress()

void ComputeDamageStress::computeQpStress ( )

overrideprotectedvirtual

Compute the stress and store it in the _stress material property for the current quadrature point.

Reimplemented from ComputeFiniteStrainElasticStress.

Definition at line 47 of file ComputeDamageStress.C.

 {
   ComputeFiniteStrainElasticStress::computeQpStress();
  
   _damage_model->setQp(_qp);
   _damage_model->updateDamage();
   _damage_model->updateStressForDamage(_stress[_qp]);
   _damage_model->finiteStrainRotation(_rotation_increment[_qp]);
   _damage_model->updateJacobianMultForDamage(_Jacobian_mult[_qp]);
  
   _matl_timestep_limit[_qp] = _damage_model->computeTimeStepLimit();
 }

◆ hasGuaranteedMaterialProperty()

bool GuaranteeConsumer::hasGuaranteedMaterialProperty	(	const MaterialPropertyName &	prop,
		Guarantee	guarantee
	)

protectedinherited

Definition at line 28 of file GuaranteeConsumer.C.

 {
   if (!_gc_feproblem->startedInitialSetup())
     mooseError("hasGuaranteedMaterialProperty() needs to be called in initialSetup()");
  
   // Reference to MaterialWarehouse for testing and retrieving block ids
   const auto & warehouse = _gc_feproblem->getMaterialWarehouse();
  
   // Complete set of ids that this object is active
   const auto & ids = (_gc_block_restrict && _gc_block_restrict->blockRestricted())
                          ? _gc_block_restrict->blockIDs()
                          : _gc_feproblem->mesh().meshSubdomains();
  
   // Loop over each id for this object
   for (const auto & id : ids)
   {
     // If block materials exist, look if any issue the required guarantee
     if (warehouse.hasActiveBlockObjects(id))
     {
       const std::vector<std::shared_ptr<MaterialBase>> & mats = warehouse.getActiveBlockObjects(id);
       for (const auto & mat : mats)
       {
         const auto & mat_props = mat->getSuppliedItems();
         if (mat_props.count(prop_name))
         {
           auto guarantee_mat = dynamic_cast<GuaranteeProvider *>(mat.get());
           if (guarantee_mat && !guarantee_mat->hasGuarantee(prop_name, guarantee))
           {
             // we found at least one material on the set of block we operate on
             // that does _not_ provide the requested guarantee
             return false;
           }
         }
       }
     }
   }
  
   return true;
 }

Referenced by ComputeFiniteStrainElasticStress::initialSetup(), ComputeSmearedCrackingStress::initialSetup(), ComputeLinearElasticPFFractureStress::initialSetup(), CriticalTimeStep::initialSetup(), and ComputeMultipleInelasticStress::initialSetup().

◆ initialSetup()

void ComputeDamageStress::initialSetup ( )

override

Definition at line 35 of file ComputeDamageStress.C.

 {
   MaterialName damage_model_name = getParam<MaterialName>("damage_model");
   DamageBase * dmb = dynamic_cast<DamageBase *>(&getMaterialByName(damage_model_name));
   if (dmb)
     _damage_model = dmb;
   else
     paramError("damage_model",
                "Damage Model " + damage_model_name + " is not compatible with ComputeDamageStress");
 }

◆ initQpStatefulProperties()

void ComputeStressBase::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Reimplemented in ComputeMultipleInelasticStress, ComputeMultipleInelasticCosseratStress, ComputeLinearViscoelasticStress, ComputeSmearedCrackingStress, ComputePFFractureStressBase, FiniteStrainHyperElasticViscoPlastic, FiniteStrainUObasedCP, FiniteStrainCrystalPlasticity, ComputeMultiPlasticityStress, FiniteStrainPlasticMaterial, and ComputeCosseratStressBase.

Definition at line 43 of file ComputeStressBase.C.

 {
   _elastic_strain[_qp].zero();
   _stress[_qp].zero();
 }

Referenced by ComputeCosseratStressBase::initQpStatefulProperties(), FiniteStrainPlasticMaterial::initQpStatefulProperties(), ComputeMultiPlasticityStress::initQpStatefulProperties(), and ComputeMultipleInelasticStress::initQpStatefulProperties().

◆ validParams()

InputParameters ComputeDamageStress::validParams ( )

static

Definition at line 18 of file ComputeDamageStress.C.

 {
   InputParameters params = ComputeFiniteStrainElasticStress::validParams();
   params.addClassDescription(
       "Compute stress for damaged elastic materials in conjunction with a damage model.");
   params.addRequiredParam<MaterialName>("damage_model", "Name of the damage model");
   return params;
 }