ComputeStrainIncrementBasedStress computes stress considering list of inelastic strain increments. More...

#include <ComputeStrainIncrementBasedStress.h>

Inheritance diagram for ComputeStrainIncrementBasedStress:

Public Member Functions
	ComputeStrainIncrementBasedStress (const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

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

virtual void	computeQpJacobian ()

virtual void	initQpStatefulProperties () override

virtual void	computeQpProperties () override

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 > &	_stress_old
	Old state of the stress tensor material property. More...

std::vector< MaterialPropertyName >	_inelastic_strain_names
	Names of the inelastic strain material properties for all inelastic models. More...

unsigned int	_num_inelastic_strain_models
	Number of inelastic models. More...

const std::string	_base_name
	Base name prepended to all material property names to allow for multi-material systems. 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...


const MaterialProperty< RankTwoTensor > &	_mechanical_strain
	Current and old states of the mechanical strain material property. More...

const MaterialProperty< RankTwoTensor > &	_mechanical_strain_old


std::vector< const MaterialProperty< RankTwoTensor > * >	_inelastic_strains
	Vectors of current and old states of the inelastic strain material properties. More...

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

Detailed Description

ComputeStrainIncrementBasedStress computes stress considering list of inelastic strain increments.

Definition at line 22 of file ComputeStrainIncrementBasedStress.h.

Constructor & Destructor Documentation

◆ ComputeStrainIncrementBasedStress()

ComputeStrainIncrementBasedStress::ComputeStrainIncrementBasedStress ( const InputParameters & parameters )

Definition at line 27 of file ComputeStrainIncrementBasedStress.C.

   : ComputeStressBase(parameters),
     _elasticity_tensor_name(_base_name + "elasticity_tensor"),
     _elasticity_tensor(getMaterialPropertyByName<RankFourTensor>(_elasticity_tensor_name)),
     _stress_old(getMaterialPropertyOld<RankTwoTensor>(_base_name + "stress")),
     _mechanical_strain(getMaterialPropertyByName<RankTwoTensor>(_base_name + "mechanical_strain")),
     _mechanical_strain_old(
         getMaterialPropertyOldByName<RankTwoTensor>(_base_name + "mechanical_strain")),
     _inelastic_strain_names(getParam<std::vector<MaterialPropertyName>>("inelastic_strain_names"))
 {
   _num_inelastic_strain_models = _inelastic_strain_names.size();
  
   if (_num_inelastic_strain_models > 0)
   {
     _inelastic_strains.resize(_num_inelastic_strain_models);
     _inelastic_strains_old.resize(_num_inelastic_strain_models);
  
     for (unsigned int i = 0; i < _num_inelastic_strain_models; ++i)
     {
       _inelastic_strains[i] = &getMaterialProperty<RankTwoTensor>(_inelastic_strain_names[i]);
       _inelastic_strains_old[i] =
           &getMaterialPropertyOld<RankTwoTensor>(_inelastic_strain_names[i]);
     }
   }
 }

Member Function Documentation

◆ computeQpJacobian()

void ComputeStrainIncrementBasedStress::computeQpJacobian ( )

protectedvirtual

Definition at line 68 of file ComputeStrainIncrementBasedStress.C.

 {
   _Jacobian_mult[_qp] = _elasticity_tensor[_qp];
 }

Referenced by computeQpStress().

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

protectedvirtual

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

Implements ComputeStressBase.

Definition at line 55 of file ComputeStrainIncrementBasedStress.C.

 {
   RankTwoTensor elastic_strain_increment = (_mechanical_strain[_qp] - _mechanical_strain_old[_qp]);
  
   for (unsigned int i = 0; i < _num_inelastic_strain_models; ++i)
     elastic_strain_increment -= (*_inelastic_strains[i])[_qp] - (*_inelastic_strains_old[i])[_qp];
  
   _stress[_qp] = _stress_old[_qp] + _elasticity_tensor[_qp] * elastic_strain_increment;
  
   computeQpJacobian();
 }

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

static

Definition at line 17 of file ComputeStrainIncrementBasedStress.C.

 {
   InputParameters params = ComputeStressBase::validParams();
   params.addClassDescription("Compute stress after subtracting inelastic strain increments");
   params.addParam<std::vector<MaterialPropertyName>>("inelastic_strain_names",
                                                      "Names of inelastic strain properties");
  
   return params;
 }