Scalar damage model computed as the combination of multiple damage models. More...

#include <CombinedScalarDamage.h>

Inheritance diagram for CombinedScalarDamage:

Public Member Functions
	CombinedScalarDamage (const InputParameters &parameters)

void	initialSetup () override

virtual void	initQpStatefulProperties () override

virtual void	updateDamage () override
	Update the internal variable(s) that evolve the damage. More...

virtual void	updateStressForDamage (RankTwoTensor &stress_new) override
	Update the current stress tensor for effects of damage. More...

virtual void	updateJacobianMultForDamage (RankFourTensor &jacobian_mult) override
	Update the material constitutive matrix. More...

virtual void	computeUndamagedOldStress (RankTwoTensor &stress_old) override

virtual Real	computeTimeStepLimit () override
	Compute the limiting value of the time step for this material. More...

const Real &	getQpDamageIndex (unsigned int qp)

const std::string	getDamageIndexName () const

virtual void	finiteStrainRotation (const RankTwoTensor &rotation_increment)
	Perform any necessary rotation of internal variables for finite strain. More...

void	setQp (unsigned int qp)
	Sets the value of the member variable _qp for use in inheriting classes. More...


void	resetQpProperties () final
	Retained as empty methods to avoid a warning from Material.C in framework. These methods are unused in all inheriting classes and should not be overwritten. More...

void	resetProperties () final

Static Public Member Functions
static InputParameters	validParams ()

Protected Types
enum	CombinationType { CombinationType::Maximum, CombinationType::Product }

Protected Member Functions
virtual void	updateQpDamageIndex () override
	Update the damage index at the current qpoint. More...

Protected Attributes
const CombinationType	_combination_type
	Type of expansion. More...

std::vector< MaterialName >	_damage_models_names

std::vector< ScalarDamageBase * >	_damage_models

const std::string	_damage_index_name
	Name of the material property where the damage index is stored. More...

const bool	_use_old_damage
	If true, use the damage index from the old state (rather than the current state) More...

const Real &	_residual_stiffness_fraction
	Residual fraction of stiffness used for material that is fully damaged. More...

const Real &	_maximum_damage_increment
	Maximum damage increment allowed for the time step. More...

const std::string	_base_name
	Base name optionally used as prefix to material tensor names. More...


MaterialProperty< Real > &	_damage_index
	Material property that provides the damage index. More...

const MaterialProperty< Real > &	_damage_index_old

const MaterialProperty< Real > &	_damage_index_older

Detailed Description

Scalar damage model computed as the combination of multiple damage models.

Definition at line 23 of file CombinedScalarDamage.h.

Member Enumeration Documentation

◆ CombinationType

enum CombinedScalarDamage::CombinationType

strongprotected

Enumerator
Maximum
Product

Definition at line 35 of file CombinedScalarDamage.h.

   {
     Maximum,
     Product
   };

Constructor & Destructor Documentation

◆ CombinedScalarDamage()

CombinedScalarDamage::CombinedScalarDamage ( const InputParameters & parameters )

Definition at line 35 of file CombinedScalarDamage.C.

   : ScalarDamageBase(parameters),
     _combination_type(getParam<MooseEnum>("combination_type").getEnum<CombinationType>()),
     _damage_models_names(getParam<std::vector<MaterialName>>("damage_models"))
 {
 }

Member Function Documentation

◆ computeTimeStepLimit()

Real ScalarDamageBase::computeTimeStepLimit ( )

overridevirtualinherited

Compute the limiting value of the time step for this material.

Returns: Limiting time step

Reimplemented from DamageBase.

Definition at line 100 of file ScalarDamageBase.C.

 {
   Real current_damage_increment = (_damage_index[_qp] - _damage_index_old[_qp]);
   if (MooseUtils::absoluteFuzzyEqual(current_damage_increment, 0.0))
     return std::numeric_limits<Real>::max();
  
   return _dt * _maximum_damage_increment / current_damage_increment;
 }

◆ computeUndamagedOldStress()

void ScalarDamageBase::computeUndamagedOldStress ( RankTwoTensor & stress_old )

overridevirtualinherited

Implements DamageBase.

Definition at line 83 of file ScalarDamageBase.C.

 {
   Real damage_index_old =
       std::max((1.0 - (_use_old_damage ? _damage_index_older[_qp] : _damage_index_old[_qp])), 0.0);
  
   if (damage_index_old > 0.0)
     stress_old /= damage_index_old;
 }

◆ finiteStrainRotation()

void DamageBase::finiteStrainRotation ( const RankTwoTensor & rotation_increment )

virtualinherited

Perform any necessary rotation of internal variables for finite strain.

Parameters

rotation_increment The finite-strain rotation increment

Definition at line 57 of file DamageBase.C.

58 {

59 }

Referenced by ComputeDamageStress::computeQpStress(), and ComputeMultipleInelasticStress::computeQpStress().

◆ getDamageIndexName()

const std::string ScalarDamageBase::getDamageIndexName ( ) const

inlineinherited

Definition at line 44 of file ScalarDamageBase.h.

44 { return _damage_index_name; }

◆ getQpDamageIndex()

const Real & ScalarDamageBase::getQpDamageIndex ( unsigned int qp )

inherited

Definition at line 60 of file ScalarDamageBase.C.

 {
   setQp(qp);
   updateQpDamageIndex();
   return _damage_index[_qp];
 }

Referenced by updateQpDamageIndex().

◆ initialSetup()

void CombinedScalarDamage::initialSetup ( )

override

Definition at line 43 of file CombinedScalarDamage.C.

 {
   for (unsigned int i = 0; i < _damage_models_names.size(); ++i)
   {
     ScalarDamageBase * model =
         dynamic_cast<ScalarDamageBase *>(&getMaterialByName(_damage_models_names[i]));
     if (model)
       _damage_models.push_back(model);
     else
       paramError("damage_model",
                  "Damage Model " + _damage_models_names[i] +
                      " is not compatible with CombinedScalarDamage");
   }
 }

◆ initQpStatefulProperties()

void ScalarDamageBase::initQpStatefulProperties ( )

overridevirtualinherited

Definition at line 54 of file ScalarDamageBase.C.

 {
   _damage_index[_qp] = 0.0;
 }

◆ resetProperties()

void DamageBase::resetProperties ( )

inlinefinalinherited

Definition at line 72 of file DamageBase.h.

72 {}

◆ resetQpProperties()

void DamageBase::resetQpProperties ( )

inlinefinalinherited

Retained as empty methods to avoid a warning from Material.C in framework. These methods are unused in all inheriting classes and should not be overwritten.

Definition at line 71 of file DamageBase.h.

71 {}

◆ setQp()

void DamageBase::setQp ( unsigned int qp )

inherited

Sets the value of the member variable _qp for use in inheriting classes.

Definition at line 40 of file DamageBase.C.

 {
   _qp = qp;
 }

Referenced by ComputeDamageStress::computeQpStress(), ComputeMultipleInelasticStress::computeQpStress(), and ScalarDamageBase::getQpDamageIndex().

◆ updateDamage()

void ScalarDamageBase::updateDamage ( )

overridevirtualinherited

Update the internal variable(s) that evolve the damage.

Reimplemented from DamageBase.

Definition at line 68 of file ScalarDamageBase.C.

 {
   updateQpDamageIndex();
 }

◆ updateJacobianMultForDamage()

void ScalarDamageBase::updateJacobianMultForDamage ( RankFourTensor & jacobian_mult )

overridevirtualinherited

Update the material constitutive matrix.

Parameters

jacobian_mult Material constitutive matrix to be modified for effects of damage

Implements DamageBase.

Definition at line 93 of file ScalarDamageBase.C.

 {
   jacobian_mult *= std::max((1.0 - (_use_old_damage ? _damage_index_old[_qp] : _damage_index[_qp])),
                             _residual_stiffness_fraction);
 }

◆ updateQpDamageIndex()

void CombinedScalarDamage::updateQpDamageIndex ( )

overrideprotectedvirtual

Update the damage index at the current qpoint.

Implements ScalarDamageBase.

Definition at line 59 of file CombinedScalarDamage.C.

 {
   switch (_combination_type)
   {
     case CombinationType::Maximum:
       _damage_index[_qp] = _damage_index_old[_qp];
       for (unsigned int i = 0; i < _damage_models.size(); ++i)
         _damage_index[_qp] = std::max(_damage_index[_qp], _damage_models[i]->getQpDamageIndex(_qp));
       break;
     case CombinationType::Product:
       _damage_index[_qp] = 1.0;
       for (unsigned int i = 0; i < _damage_models.size(); ++i)
         _damage_index[_qp] *= 1.0 - _damage_models[i]->getQpDamageIndex(_qp);
       _damage_index[_qp] = 1.0 - _damage_index[_qp];
       break;
   }
  
   _damage_index[_qp] =
       std::max(_damage_index_old[_qp], std::max(0.0, std::min(1.0, _damage_index[_qp])));
 }

◆ updateStressForDamage()

void ScalarDamageBase::updateStressForDamage ( RankTwoTensor & stress_new )

overridevirtualinherited

Update the current stress tensor for effects of damage.

Parameters

stress_new Undamaged stress to be modified by the damage model

Implements DamageBase.

Definition at line 74 of file ScalarDamageBase.C.

 {
   // Avoid multiplying by a small negative number, which could occur if damage_index
   // is slightly greater than 1.0
   stress_new *=
       std::max((1.0 - (_use_old_damage ? _damage_index_old[_qp] : _damage_index[_qp])), 0.0);
 }

◆ validParams()

InputParameters CombinedScalarDamage::validParams ( )

static

Definition at line 17 of file CombinedScalarDamage.C.

 {
   InputParameters params = ScalarDamageBase::validParams();
  
   params.addClassDescription(
       "Scalar damage model which is computed as a function of multiple scalar damage models");
  
   params.addRequiredParam<std::vector<MaterialName>>("damage_models",
                                                      "Name of the damage models used to compute "
                                                      "the damage index");
  
   MooseEnum combination_type("Maximum Product", "Maximum");
   params.addParam<MooseEnum>(
       "combination_type", combination_type, "How the damage models are combined");
  
   return params;
 }

Member Data Documentation

◆ _base_name

const std::string DamageBase::_base_name

protectedinherited

Base name optionally used as prefix to material tensor names.

Definition at line 77 of file DamageBase.h.

Referenced by ScalarMaterialDamage::updateQpDamageIndex().

◆ _combination_type

const CombinationType CombinedScalarDamage::_combination_type

protected

Type of expansion.

Definition at line 42 of file CombinedScalarDamage.h.

Referenced by updateQpDamageIndex().

◆ _damage_index

MaterialProperty<Real>& ScalarDamageBase::_damage_index

protectedinherited

Material property that provides the damage index.

Definition at line 54 of file ScalarDamageBase.h.

Referenced by ScalarDamageBase::computeTimeStepLimit(), ScalarDamageBase::getQpDamageIndex(), ScalarDamageBase::initQpStatefulProperties(), ScalarDamageBase::updateJacobianMultForDamage(), ScalarMaterialDamage::updateQpDamageIndex(), updateQpDamageIndex(), and ScalarDamageBase::updateStressForDamage().

◆ _damage_index_name

const std::string ScalarDamageBase::_damage_index_name

protectedinherited

Name of the material property where the damage index is stored.

Definition at line 48 of file ScalarDamageBase.h.

Referenced by ScalarDamageBase::getDamageIndexName().

◆ _damage_index_old

const MaterialProperty<Real>& ScalarDamageBase::_damage_index_old

protectedinherited

Definition at line 55 of file ScalarDamageBase.h.

Referenced by ScalarDamageBase::computeTimeStepLimit(), ScalarDamageBase::computeUndamagedOldStress(), ScalarDamageBase::updateJacobianMultForDamage(), updateQpDamageIndex(), and ScalarDamageBase::updateStressForDamage().

◆ _damage_index_older

const MaterialProperty<Real>& ScalarDamageBase::_damage_index_older

protectedinherited

Definition at line 56 of file ScalarDamageBase.h.

Referenced by ScalarDamageBase::computeUndamagedOldStress().

◆ _damage_models

std::vector<ScalarDamageBase *> CombinedScalarDamage::_damage_models

protected

Definition at line 46 of file CombinedScalarDamage.h.

Referenced by initialSetup(), and updateQpDamageIndex().

◆ _damage_models_names

std::vector<MaterialName> CombinedScalarDamage::_damage_models_names

protected

Definition at line 44 of file CombinedScalarDamage.h.

Referenced by initialSetup().

◆ _maximum_damage_increment

const Real& ScalarDamageBase::_maximum_damage_increment

protectedinherited

Maximum damage increment allowed for the time step.

Definition at line 66 of file ScalarDamageBase.h.

Referenced by ScalarDamageBase::computeTimeStepLimit().

◆ _residual_stiffness_fraction

const Real& ScalarDamageBase::_residual_stiffness_fraction

protectedinherited

Residual fraction of stiffness used for material that is fully damaged.

Definition at line 63 of file ScalarDamageBase.h.

Referenced by ScalarDamageBase::updateJacobianMultForDamage().

◆ _use_old_damage

const bool ScalarDamageBase::_use_old_damage

protectedinherited

If true, use the damage index from the old state (rather than the current state)

Definition at line 60 of file ScalarDamageBase.h.

Referenced by ScalarDamageBase::computeUndamagedOldStress(), ScalarDamageBase::updateJacobianMultForDamage(), and ScalarDamageBase::updateStressForDamage().

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

tensor_mechanics/include/materials/CombinedScalarDamage.h
tensor_mechanics/src/materials/CombinedScalarDamage.C

Public Member Functions

Static Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

Member Enumeration Documentation

◆ CombinationType

Constructor & Destructor Documentation

◆ CombinedScalarDamage()

Member Function Documentation

◆ computeTimeStepLimit()

◆ computeUndamagedOldStress()

◆ finiteStrainRotation()

◆ getDamageIndexName()

◆ getQpDamageIndex()

◆ initialSetup()

◆ initQpStatefulProperties()

◆ resetProperties()

◆ resetQpProperties()

◆ setQp()

◆ updateDamage()

◆ updateJacobianMultForDamage()

◆ updateQpDamageIndex()

◆ updateStressForDamage()

◆ validParams()

Member Data Documentation

◆ _base_name

◆ _combination_type

◆ _damage_index

◆ _damage_index_name

◆ _damage_index_old

◆ _damage_index_older

◆ _damage_models

◆ _damage_models_names

◆ _maximum_damage_increment

◆ _residual_stiffness_fraction

◆ _use_old_damage