ComputeFiniteStrainElasticStress computes the stress following elasticity theory for finite strains. More...

#include <ComputeFiniteStrainElasticStress.h>

Inheritance diagram for ComputeFiniteStrainElasticStress:

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

ComputeFiniteStrainElasticStress computes the stress following elasticity theory for finite strains.

Definition at line 24 of file ComputeFiniteStrainElasticStress.h.

Constructor & Destructor Documentation

◆ ComputeFiniteStrainElasticStress()

ComputeFiniteStrainElasticStress::ComputeFiniteStrainElasticStress ( const InputParameters & parameters )

Definition at line 24 of file ComputeFiniteStrainElasticStress.C.

   : ComputeStressBase(parameters),
     GuaranteeConsumer(this),
     _elasticity_tensor_name(_base_name + "elasticity_tensor"),
     _elasticity_tensor(getMaterialPropertyByName<RankFourTensor>(_elasticity_tensor_name)),
     _strain_increment(getMaterialPropertyByName<RankTwoTensor>(_base_name + "strain_increment")),
     _rotation_increment(
         getMaterialPropertyByName<RankTwoTensor>(_base_name + "rotation_increment")),
     _stress_old(getMaterialPropertyOldByName<RankTwoTensor>(_base_name + "stress")),
     _elastic_strain_old(getMaterialPropertyOldByName<RankTwoTensor>(_base_name + "elastic_strain"))
 {
 }

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 ComputeFiniteStrainElasticStress::computeQpStress ( )

overrideprotectedvirtual

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

Implements ComputeStressBase.

Reimplemented in ComputeMultipleInelasticStress, ComputeMultipleInelasticCosseratStress, ComputeSmearedCrackingStress, and ComputeDamageStress.

Definition at line 47 of file ComputeFiniteStrainElasticStress.C.

 {
   // Calculate the stress in the intermediate configuration
   RankTwoTensor intermediate_stress;
  
   intermediate_stress =
       _elasticity_tensor[_qp] * (_elastic_strain_old[_qp] + _strain_increment[_qp]);
  
   // Rotate the stress state to the current configuration
   _stress[_qp] =
       _rotation_increment[_qp] * intermediate_stress * _rotation_increment[_qp].transpose();
  
   // Assign value for elastic strain, which is equal to the mechanical strain
   _elastic_strain[_qp] = _mechanical_strain[_qp];
  
   // Compute dstress_dstrain
   _Jacobian_mult[_qp] = _elasticity_tensor[_qp]; // This is NOT the exact jacobian
 }

Referenced by ComputeDamageStress::computeQpStress().

◆ 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 initialSetup(), ComputeSmearedCrackingStress::initialSetup(), ComputeLinearElasticPFFractureStress::initialSetup(), CriticalTimeStep::initialSetup(), and ComputeMultipleInelasticStress::initialSetup().

◆ initialSetup()

void ComputeFiniteStrainElasticStress::initialSetup ( )

override

Definition at line 39 of file ComputeFiniteStrainElasticStress.C.

 {
   if (!hasGuaranteedMaterialProperty(_elasticity_tensor_name, Guarantee::ISOTROPIC))
     mooseError("ComputeFiniteStrainElasticStress can only be used with elasticity tensor materials "
                "that guarantee isotropic tensors.");
 }

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

static

Definition at line 17 of file ComputeFiniteStrainElasticStress.C.

 {
   InputParameters params = ComputeStressBase::validParams();
   params.addClassDescription("Compute stress using elasticity for finite strains");
   return params;
 }