LinearViscoelasticityBase Class Referenceabstract

This class is a base class for materials consisting of an assembly of linear springs and dashpots. More...

#include <LinearViscoelasticityBase.h>

Inheritance diagram for LinearViscoelasticityBase:
[legend]

Public Types

enum  IntegrationRule { IntegrationRule::BackwardEuler, IntegrationRule::MidPoint, IntegrationRule::Newmark, IntegrationRule::Zienkiewicz }
Determines how theta is calculated for the time-integration system. More...

Public Member Functions

LinearViscoelasticityBase (const InputParameters &parameters)

void recomputeQpApparentProperties (unsigned int qp)
Compute the apparent properties at a quadrature point. More...

bool hasGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

Protected Member Functions

virtual void initQpStatefulProperties () override

virtual void computeQpElasticityTensor () final
Inherited from ComputeElasticityTensorBase. More...

virtual void computeQpViscoelasticProperties ()=0
This method assigns the mechanical properties of each spring and dashpot in the system. More...

virtual void computeQpViscoelasticPropertiesInv ()
This method computes the inverse elasticity tensor of each spring in the system (if required). More...

virtual void computeQpApparentElasticityTensors ()=0
This method computes the apparent elasticity tensor used in the internal time-stepping scheme. More...

virtual void computeQpApparentCreepStrain ()=0
This method computes the apparent creep strain corresponding to the current viscous_strain of each dashpot. More...

virtual void updateQpViscousStrains ()=0
Update the internal viscous strains at a quadrature point. More...

void declareViscoelasticProperties ()
Declare all necessary MaterialProperties for the model. More...

Real computeTheta (Real dt, Real viscosity) const
Provides theta as a function of the time step and a viscosity. More...

virtual void computeQpProperties ()

void issueGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

void revokeGuarantee (const MaterialPropertyName &prop_name, Guarantee guarantee)

Protected Attributes

IntegrationRule _integration_rule
Determines how theta is computed. More...

Real _theta
User-defined value for theta. More...

MaterialProperty< RankFourTensor > & _apparent_elasticity_tensor
Apparent elasticity tensor. This is NOT the elasticity tensor of the material. More...

MaterialProperty< RankFourTensor > & _apparent_elasticity_tensor_inv
Inverse of the apparent elasticity tensor. More...

MaterialProperty< RankFourTensor > & _elasticity_tensor_inv
Instantaneous elasticity tensor. This IS the real elasticity tensor of the material. More...

bool _need_viscoelastic_properties_inverse
If active, indicates that we need to call computeQpViscoelasticPropertiesInv() More...

bool _has_longterm_dashpot
Indicates if the spring-dashpot assembly has a single dashpot not associated with a spring. More...

unsigned int _components
This is the number of internal variables required by the model. More...

const MaterialProperty< RankTwoTensor > & _elastic_strain_old
previous value of the elastic strain for update purposes More...

const MaterialProperty< RankTwoTensor > & _creep_strain_old
Previous value of the true creep strain for update purposes. More...

bool _has_driving_eigenstrain
Indicates if the model is only driven by the stress, or also by an additional eigenstrain. More...

std::string _driving_eigenstrain_name
Name of the eigenstrain that drives the additional creep strain. More...

bool _force_recompute_properties
If activated, the time-stepping scheme will be re-initialized at each step of the solver. More...

bool & _step_zero
checks whether we are at the first time step More...

std::string _base_name

std::string _elasticity_tensor_name

MaterialProperty< RankFourTensor > & _elasticity_tensor

Function *const _prefactor_function
prefactor function to multiply the elasticity tensor with More...

MaterialProperty< RankFourTensor > & _first_elasticity_tensor
Elasticity tensor of a stand-alone elastic spring in the chain. More...

MaterialProperty< RankFourTensor > * _first_elasticity_tensor_inv

std::vector< MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors
List of elasticity tensor of each subsequent spring in the chain. More...

std::vector< MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors_inv

std::vector< const MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors_inv_old

std::vector< MaterialProperty< Real > * > _dashpot_viscosities
List of viscosities of each subsequent dashpot in the chain. More...

std::vector< const MaterialProperty< Real > * > _dashpot_viscosities_old

std::vector< MaterialProperty< RankTwoTensor > * > _viscous_strains

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

MaterialProperty< RankTwoTensor > & _apparent_creep_strain
The apparent creep strain resulting from the internal viscous strains. More...

const MaterialProperty< RankTwoTensor > & _apparent_creep_strain_old

const MaterialProperty< RankTwoTensor > * _driving_eigenstrain
Pointer to the value of the driving eigenstrain. More...

const MaterialProperty< RankTwoTensor > * _driving_eigenstrain_old

Detailed Description

This class is a base class for materials consisting of an assembly of linear springs and dashpots.

It represents the arrangement of the system (typically, in parallel or in series), the mechanical * properties of each individual spring and dashpot in the model, and the internal strains associated with each dashpot.

To use a linear viscoelastic model, the user must provide a LinearViscoelasticityManager, that will initialize and update the internal time-stepping scheme at the beginning and the end of each time step.

To compute the stress associated with a linear viscoelastic model, the user must either use a ComputeLinearViscoelasticStress material (for total small strain formulations), or a ComputeMultipleInelasticStress material associated with a LinearViscoelasticStressUpdate material if there are multiple sources of inelastic strains (this requires to use incremental strains).

The rate-dependent problem is solved using an internal 1-step finite difference scheme. The scheme itself depends only on the arrangement of the system, and not the individual properties. The scheme acts upon a series of internal strain variables, stored in the _viscous_strains material properties. At the start of each time step, the apparent elastic properties of the assembly and the apparent creep strain are calculated. The rate-dependent problem becomes then equivalent to a purely elastic system. At the end of the time step, the internal strain variables are updated according to the current values of the stress and the strain.

The real system refers to stress = elasticity_tensor : (mechanical_strain - creep_strain)

The internal time-stepping scheme uses the following conversion: stress = apparent_elasticity_tensor : (mechanical_strain - apparent_creep_strain)

The creep strain is generally driven only by the mechanical stress, but it can be driven by an eigenstrain in addition to the stress. This means that the creep strain is calculated as if the eigenstrain were to contribute to the global mechanical stress of the system.

The determination of the apparent properties, as well as the update of the internal strain variables depend only on the arrangement of the system (in parralel or series), and not their actual properties. These calculations are generally only performed at the beginning of the time step, but the user (or derived classes) can enforce these to be performed at each step of the nonlinear solver if necessary.

The time-integration uses a 1-step Newmark finite difference scheme. This scheme is controlled by a parameter theta (between 0 and 1, default-value 1). Theta can be automatically calibrated depending on the value of the dashpot viscosities in order to reproduce the exact integral of exponential series (integration_rule = "zienkiewicz"). Theta cannot be set to 0 (purely explicit system), and setting theta < 0.5 may cause the scheme to diverge (the "backward_euler" "mid_point" and "zienkiewicz" schemes are unconditionally stable).

This class is virtual. Deriving this class must be done in two steps.

1. Classes inheriting directly from LinearViscoelasticityBase must represent a distinct spring-dashpot architecture (examples in GeneralizedKelvinVoigtBase and GeneralizedMaxwellBase), and must inherit only methods related to the time-stepping scheme itself.
2. Classes inheriting inderectly from LinearViscoelasticityBase must provide the actual properties of each spring and dashpot in the assembly (examples in GeneralizedKelvinVoigtModel and GeneralizedMaxwellModel). These classes must not override methods related to the time-stepping scheme.

This class does not store the true creep strain of the material, but an equivalent creep strain adapted for numerical efficiency. The true creep strain itself is declared by ComputeLinearViscoelasticStress or LinearViscoelasticStressUpdate.

Definition at line 82 of file LinearViscoelasticityBase.h.

◆ IntegrationRule

 strong

Determines how theta is calculated for the time-integration system.

Enumerator
BackwardEuler

theta = 1

MidPoint

theta = 0.5

Newmark

theta defined by the user

Zienkiewicz

theta automatically adjusted as a function of the time step and the viscosity

Definition at line 88 of file LinearViscoelasticityBase.h.

89  {
91  BackwardEuler,
93  MidPoint,
95  Newmark,
97  Zienkiewicz,
98  };

◆ LinearViscoelasticityBase()

 LinearViscoelasticityBase::LinearViscoelasticityBase ( const InputParameters & parameters )

Definition at line 49 of file LinearViscoelasticityBase.C.

50  : ComputeElasticityTensorBase(parameters),
51  _integration_rule(getParam<MooseEnum>("integration_rule").getEnum<IntegrationRule>()),
52  _theta(getParam<Real>("theta")),
54  declareProperty<RankFourTensor>(_base_name + "apparent_elasticity_tensor")),
56  declareProperty<RankFourTensor>(_base_name + "apparent_elasticity_tensor_inv")),
57  _elasticity_tensor_inv(declareProperty<RankFourTensor>(_elasticity_tensor_name + "_inv")),
58  _need_viscoelastic_properties_inverse(getParam<bool>("need_viscoelastic_properties_inverse")),
59  _has_longterm_dashpot(false),
60  _components(0),
62  declareProperty<RankFourTensor>(_base_name + "spring_elasticity_tensor_0")),
65  ? &declareProperty<RankFourTensor>(_base_name + "spring_elasticity_tensor_0_inv")
66  : nullptr),
67  _apparent_creep_strain(declareProperty<RankTwoTensor>(_base_name + "apparent_creep_strain")),
69  getMaterialPropertyOld<RankTwoTensor>(_base_name + "apparent_creep_strain")),
71  getMaterialPropertyOld<RankTwoTensor>(getParam<std::string>("elastic_strain_name"))),
73  getMaterialPropertyOld<RankTwoTensor>(getParam<std::string>("creep_strain_name"))),
74  _has_driving_eigenstrain(isParamValid("driving_eigenstrain")),
76  _has_driving_eigenstrain ? getParam<std::string>("driving_eigenstrain") : ""),
78  ? &getMaterialPropertyByName<RankTwoTensor>(_driving_eigenstrain_name)
79  : nullptr),
81  ? &getMaterialPropertyOld<RankTwoTensor>(_driving_eigenstrain_name)
82  : nullptr),
83  _force_recompute_properties(getParam<bool>("force_recompute_properties")),
84  _step_zero(declareRestartableData<bool>("step_zero", true))
85 {
86  if (_theta < 0.5)
87  mooseWarning("theta parameter for LinearViscoelasticityBase is below 0.5; time integration may "
88  "not converge!");
89
90  // force material properties to be considered stateful
91  getMaterialPropertyOld<RankFourTensor>(_base_name + "apparent_elasticity_tensor");
92  getMaterialPropertyOld<RankFourTensor>(_base_name + "apparent_elasticity_tensor_inv");
93  getMaterialPropertyOld<RankFourTensor>(_elasticity_tensor_name);
94  getMaterialPropertyOld<RankFourTensor>(_elasticity_tensor_name + "_inv");
95  getMaterialPropertyOld<RankFourTensor>(_base_name + "spring_elasticity_tensor_0");
97  getMaterialPropertyOld<RankFourTensor>(_base_name + "spring_elasticity_tensor_0_inv");
98 }
bool & _step_zero
checks whether we are at the first time step
const MaterialProperty< RankTwoTensor > & _apparent_creep_strain_old
bool _has_longterm_dashpot
Indicates if the spring-dashpot assembly has a single dashpot not associated with a spring...
bool _force_recompute_properties
If activated, the time-stepping scheme will be re-initialized at each step of the solver...
bool _need_viscoelastic_properties_inverse
If active, indicates that we need to call computeQpViscoelasticPropertiesInv()
MaterialProperty< RankFourTensor > & _apparent_elasticity_tensor
Apparent elasticity tensor. This is NOT the elasticity tensor of the material.
const MaterialProperty< RankTwoTensor > & _elastic_strain_old
previous value of the elastic strain for update purposes
MaterialProperty< RankFourTensor > & _first_elasticity_tensor
Elasticity tensor of a stand-alone elastic spring in the chain.
MaterialProperty< RankFourTensor > * _first_elasticity_tensor_inv
const MaterialProperty< RankTwoTensor > * _driving_eigenstrain_old
Real _theta
User-defined value for theta.
std::string _driving_eigenstrain_name
Name of the eigenstrain that drives the additional creep strain.
unsigned int _components
This is the number of internal variables required by the model.
MaterialProperty< RankTwoTensor > & _apparent_creep_strain
The apparent creep strain resulting from the internal viscous strains.
const MaterialProperty< RankTwoTensor > * _driving_eigenstrain
Pointer to the value of the driving eigenstrain.
ComputeElasticityTensorBase(const InputParameters &parameters)
bool _has_driving_eigenstrain
Indicates if the model is only driven by the stress, or also by an additional eigenstrain.
MaterialProperty< RankFourTensor > & _elasticity_tensor_inv
Instantaneous elasticity tensor. This IS the real elasticity tensor of the material.
IntegrationRule _integration_rule
Determines how theta is computed.
const MaterialProperty< RankTwoTensor > & _creep_strain_old
Previous value of the true creep strain for update purposes.
MaterialProperty< RankFourTensor > & _apparent_elasticity_tensor_inv
Inverse of the apparent elasticity tensor.

◆ computeQpApparentCreepStrain()

 virtual void LinearViscoelasticityBase::computeQpApparentCreepStrain ( )
protectedpure virtual

This method computes the apparent creep strain corresponding to the current viscous_strain of each dashpot.

It must be called after the apparent elasticity tensors have been calculated.

This method is purely virtual. Inherited classes must override it.

This method is related to the internal time-stepping scheme. It should only be overwritten by classes that inherit directly from LinearViscoelasticityBase, and that represent a different spring-dashpot assembly. See GeneralizedKelvinVoigtBase for example.

Implemented in GeneralizedKelvinVoigtBase, and GeneralizedMaxwellBase.

Referenced by recomputeQpApparentProperties().

◆ computeQpApparentElasticityTensors()

 virtual void LinearViscoelasticityBase::computeQpApparentElasticityTensors ( )
protectedpure virtual

This method computes the apparent elasticity tensor used in the internal time-stepping scheme.

It is called after the mechanical properties have been set, and before the apparent creep strains are calculated.

This method is also responsible for calculating the instantaneous elasticity tensor, and the inverse of both the apparent and instantaneous elasticity tensors.

This method is purely virtual. Inherited classes must override it.

This method is related to the internal time-stepping scheme. It should only be overwritten by classes that inherit directly from LinearViscoelasticityBase, and that represent a different spring-dashpot assembly. See GeneralizedKelvinVoigtBase for example.

Implemented in GeneralizedKelvinVoigtBase, and GeneralizedMaxwellBase.

Referenced by recomputeQpApparentProperties().

◆ computeQpElasticityTensor()

 void LinearViscoelasticityBase::computeQpElasticityTensor ( )
finalprotectedvirtual

Inherited from ComputeElasticityTensorBase.

Implements ComputeElasticityTensorBase.

Definition at line 191 of file LinearViscoelasticityBase.C.

192 {
195 }
void recomputeQpApparentProperties(unsigned int qp)
Compute the apparent properties at a quadrature point.
bool _force_recompute_properties
If activated, the time-stepping scheme will be re-initialized at each step of the solver...

◆ computeQpProperties()

 void ComputeElasticityTensorBase::computeQpProperties ( )
protectedvirtualinherited

Definition at line 41 of file ComputeElasticityTensorBase.C.

42 {
44
45  // Multiply by prefactor
47  _elasticity_tensor[_qp] *= _prefactor_function->value(_t, _q_point[_qp]);
48 }
Function *const _prefactor_function
prefactor function to multiply the elasticity tensor with
virtual void computeQpElasticityTensor()=0
MaterialProperty< RankFourTensor > & _elasticity_tensor

◆ computeQpViscoelasticProperties()

 virtual void LinearViscoelasticityBase::computeQpViscoelasticProperties ( )
protectedpure virtual

This method assigns the mechanical properties of each spring and dashpot in the system.

This method is purely virtual. Inherited classes must override it.

This method is related to the storage of the mechanical properties of each spring and dashpot in the system, and not the internal time-stepping procedure. Only end-user classes should override it. See GeneralizedKelvinVoigModel for example.

Implemented in GeneralizedMaxwellModel, and GeneralizedKelvinVoigtModel.

Referenced by recomputeQpApparentProperties().

◆ computeQpViscoelasticPropertiesInv()

 void LinearViscoelasticityBase::computeQpViscoelasticPropertiesInv ( )
protectedvirtual

This method computes the inverse elasticity tensor of each spring in the system (if required).

This method is virtual. Its default behavior computes the inverse of each tensor. It must be inherited only if there is a faster way to compute this inverse (for example, if they are known).

This method is related to the storage of the mechanical properties of each spring and dashpot in the system, and not the internal time-stepping procedure. Only end-user classes should override it. See GeneralizedKelvinVoigtModel for example.

Reimplemented in GeneralizedMaxwellModel, and GeneralizedKelvinVoigtModel.

Definition at line 198 of file LinearViscoelasticityBase.C.

Referenced by recomputeQpApparentProperties().

199 {
200  if (MooseUtils::absoluteFuzzyEqual(_first_elasticity_tensor[_qp].L2norm(), 0.0))
201  (*_first_elasticity_tensor_inv)[_qp].zero();
202  else
204
205  for (unsigned int i = 0; i < _springs_elasticity_tensors.size(); ++i)
206  {
207  if (MooseUtils::absoluteFuzzyEqual((*_springs_elasticity_tensors[i])[_qp].L2norm(), 0.0))
208  (*_springs_elasticity_tensors_inv[i])[_qp].zero();
209  else
210  (*_springs_elasticity_tensors_inv[i])[_qp] = (*_springs_elasticity_tensors[i])[_qp].invSymm();
211  }
212 }
std::vector< MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors_inv
std::vector< MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors
List of elasticity tensor of each subsequent spring in the chain.
MaterialProperty< RankFourTensor > & _first_elasticity_tensor
Elasticity tensor of a stand-alone elastic spring in the chain.
MaterialProperty< RankFourTensor > * _first_elasticity_tensor_inv
Real L2norm(const RankTwoTensor &r2tensor)

◆ computeTheta()

 Real LinearViscoelasticityBase::computeTheta ( Real dt, Real viscosity ) const
protected

Provides theta as a function of the time step and a viscosity.

Definition at line 215 of file LinearViscoelasticityBase.C.

216 {
217  if (MooseUtils::absoluteFuzzyEqual(dt, 0.0))
218  mooseError("linear viscoelasticity cannot be integrated over a dt of ", dt);
219
220  switch (_integration_rule)
221  {
223  return 1.;
225  return 0.5;
227  return _theta;
229  return 1. / (1. - std::exp(-dt / viscosity)) - viscosity / dt;
230  default:
231  return 1.;
232  }
233  return 1.;
234 }
theta automatically adjusted as a function of the time step and the viscosity
Real _theta
User-defined value for theta.
IntegrationRule _integration_rule
Determines how theta is computed.

◆ declareViscoelasticProperties()

 void LinearViscoelasticityBase::declareViscoelasticProperties ( )
protected

Declare all necessary MaterialProperties for the model.

This method must be called once at the end of the constructor of a final inherited class, after _components has been set. See GeneralizedKelvinVoigtModel or GeneralizedMaxwell model for example.

Definition at line 101 of file LinearViscoelasticityBase.C.

102 {
103  for (unsigned int i = 0; i < _components; ++i)
104  {
105  std::string ith = Moose::stringify(i + 1);
106
107  if (!_has_longterm_dashpot || (_components > 0 && i < _components - 1))
108  {
109  _springs_elasticity_tensors.push_back(
110  &declareProperty<RankFourTensor>(_base_name + "spring_elasticity_tensor_" + ith));
111  getMaterialPropertyOld<RankFourTensor>(_base_name + "spring_elasticity_tensor_" + ith);
112  }
113
114  _dashpot_viscosities.push_back(&declareProperty<Real>(_base_name + "dashpot_viscosity_" + ith));
115  _dashpot_viscosities_old.push_back(
116  &getMaterialPropertyOld<Real>(_base_name + "dashpot_viscosity_" + ith));
117
118  _viscous_strains.push_back(
119  &declareProperty<RankTwoTensor>(_base_name + "viscous_strain_" + ith));
120  _viscous_strains_old.push_back(
121  &getMaterialPropertyOld<RankTwoTensor>(_base_name + "viscous_strain_" + ith));
122
124  {
125  _springs_elasticity_tensors_inv.push_back(&declareProperty<RankFourTensor>(
126  _base_name + "spring_elasticity_tensor_" + ith + "_inv"));
127  _springs_elasticity_tensors_inv_old.push_back(&getMaterialPropertyOld<RankFourTensor>(
128  _base_name + "spring_elasticity_tensor_" + ith + "_inv"));
129  }
130  }
131 }
std::vector< MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors_inv
std::vector< MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors
List of elasticity tensor of each subsequent spring in the chain.
std::vector< const MaterialProperty< RankTwoTensor > * > _viscous_strains_old
std::vector< const MaterialProperty< Real > * > _dashpot_viscosities_old
bool _has_longterm_dashpot
Indicates if the spring-dashpot assembly has a single dashpot not associated with a spring...
bool _need_viscoelastic_properties_inverse
If active, indicates that we need to call computeQpViscoelasticPropertiesInv()
std::vector< MaterialProperty< Real > * > _dashpot_viscosities
List of viscosities of each subsequent dashpot in the chain.
std::vector< const MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors_inv_old
unsigned int _components
This is the number of internal variables required by the model.
std::vector< MaterialProperty< RankTwoTensor > * > _viscous_strains

◆ hasGuarantee()

 bool GuaranteeProvider::hasGuarantee ( const MaterialPropertyName & prop_name, Guarantee guarantee )
inherited

Definition at line 16 of file GuaranteeProvider.C.

17 {
18  auto it = _guarantees.find(prop_name);
19  if (it == _guarantees.end())
20  return false;
21
22  auto it2 = it->second.find(guarantee);
23  return it2 != it->second.end();
24 }
std::map< MaterialPropertyName, std::set< Guarantee > > _guarantees

◆ initQpStatefulProperties()

 void LinearViscoelasticityBase::initQpStatefulProperties ( )
overrideprotectedvirtual

Definition at line 134 of file LinearViscoelasticityBase.C.

135 {
136  if (_components != _viscous_strains.size())
137  mooseError(
138  "inconsistent numbers of dashpots and viscous strains in LinearViscoelasticityBase;"
139  " Make sure declareViscoelasticProperties has been called in the viscoelastic model");
140
141  _apparent_creep_strain[_qp].zero();
142  _apparent_elasticity_tensor[_qp].zero();
144  _elasticity_tensor_inv[_qp].zero();
145  _first_elasticity_tensor[_qp].zero();
147  (*_first_elasticity_tensor_inv)[_qp].zero();
148
149  for (unsigned int i = 0; i < _components; ++i)
150  {
151  if (!_has_longterm_dashpot || (_components > 0 && i < _components - 1))
152  {
153  (*_springs_elasticity_tensors[i])[_qp].zero();
155  (*_springs_elasticity_tensors_inv[i])[_qp].zero();
156  }
157
158  (*_dashpot_viscosities[i])[_qp] = 0.0;
159  (*_viscous_strains[i])[_qp].zero();
160  }
161 }
std::vector< MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors_inv
std::vector< MaterialProperty< RankFourTensor > * > _springs_elasticity_tensors
List of elasticity tensor of each subsequent spring in the chain.
bool _has_longterm_dashpot
Indicates if the spring-dashpot assembly has a single dashpot not associated with a spring...
bool _need_viscoelastic_properties_inverse
If active, indicates that we need to call computeQpViscoelasticPropertiesInv()
MaterialProperty< RankFourTensor > & _apparent_elasticity_tensor
Apparent elasticity tensor. This is NOT the elasticity tensor of the material.
std::vector< MaterialProperty< Real > * > _dashpot_viscosities
List of viscosities of each subsequent dashpot in the chain.
MaterialProperty< RankFourTensor > & _first_elasticity_tensor
Elasticity tensor of a stand-alone elastic spring in the chain.
unsigned int _components
This is the number of internal variables required by the model.
MaterialProperty< RankTwoTensor > & _apparent_creep_strain
The apparent creep strain resulting from the internal viscous strains.
MaterialProperty< RankFourTensor > & _elasticity_tensor_inv
Instantaneous elasticity tensor. This IS the real elasticity tensor of the material.
std::vector< MaterialProperty< RankTwoTensor > * > _viscous_strains
MaterialProperty< RankFourTensor > & _apparent_elasticity_tensor_inv
Inverse of the apparent elasticity tensor.

◆ issueGuarantee()

 void GuaranteeProvider::issueGuarantee ( const MaterialPropertyName & prop_name, Guarantee guarantee )
protectedinherited

Definition at line 27 of file GuaranteeProvider.C.

28 {
29  // intentional insertion
30  _guarantees[prop_name].insert(guarantee);
31 }
std::map< MaterialPropertyName, std::set< Guarantee > > _guarantees

◆ recomputeQpApparentProperties()

 void LinearViscoelasticityBase::recomputeQpApparentProperties ( unsigned int qp )

Compute the apparent properties at a quadrature point.

This initializes the internal time-stepping scheme, and must be called at the beginning of the time step.

This method is called by LinearViscoelasticityManager.

Definition at line 164 of file LinearViscoelasticityBase.C.

Referenced by computeQpElasticityTensor().

165 {
166  unsigned int qp_prev = _qp;
167  _qp = qp;
168
169  if (_t_step >= 1)
170  _step_zero = false;
171
172  // 1. we get the viscoelastic properties and their inverse if needed
176
177  // 2. we update the internal viscous strains from the previous time step
179
180  // 3. we compute the apparent elasticity tensor
182
183  // 4. we transform the internal viscous strains in an apparent creep strain
184  if (!_step_zero)
186
187  _qp = qp_prev;
188 }
bool & _step_zero
checks whether we are at the first time step
bool _need_viscoelastic_properties_inverse
If active, indicates that we need to call computeQpViscoelasticPropertiesInv()
virtual void computeQpApparentElasticityTensors()=0
This method computes the apparent elasticity tensor used in the internal time-stepping scheme...
virtual void computeQpViscoelasticProperties()=0
This method assigns the mechanical properties of each spring and dashpot in the system.
virtual void computeQpApparentCreepStrain()=0
This method computes the apparent creep strain corresponding to the current viscous_strain of each da...
virtual void computeQpViscoelasticPropertiesInv()
This method computes the inverse elasticity tensor of each spring in the system (if required)...
virtual void updateQpViscousStrains()=0
Update the internal viscous strains at a quadrature point.

◆ revokeGuarantee()

 void GuaranteeProvider::revokeGuarantee ( const MaterialPropertyName & prop_name, Guarantee guarantee )
protectedinherited

Definition at line 34 of file GuaranteeProvider.C.

Referenced by ComputeElasticityTensorCP::ComputeElasticityTensorCP().

35 {
36  auto it = _guarantees.find(prop_name);
37  if (it != _guarantees.end())
38  it->second.erase(guarantee);
39 }
std::map< MaterialPropertyName, std::set< Guarantee > > _guarantees

◆ updateQpViscousStrains()

 virtual void LinearViscoelasticityBase::updateQpViscousStrains ( )
protectedpure virtual

Update the internal viscous strains at a quadrature point.

Calling this method is required at the end of each time step to update the internal time-stepping scheme correctly.

This method is pure virtual. Inherited classes must override it.

This method is related to the internal time-stepping scheme. It should only be overwritten by classes that inherit directly from LinearViscoelasticityBase, and that represent a different spring-dashpot assembly. See GeneralizedKelvinVoigtBase or GeneralizedMaxwellBase for example.

Implemented in GeneralizedKelvinVoigtBase, and GeneralizedMaxwellBase.

Referenced by recomputeQpApparentProperties().

◆ _apparent_creep_strain

 MaterialProperty& LinearViscoelasticityBase::_apparent_creep_strain
protected

The apparent creep strain resulting from the internal viscous strains.

Definition at line 240 of file LinearViscoelasticityBase.h.

◆ _apparent_creep_strain_old

 const MaterialProperty& LinearViscoelasticityBase::_apparent_creep_strain_old
protected

Definition at line 241 of file LinearViscoelasticityBase.h.

◆ _apparent_elasticity_tensor

 MaterialProperty& LinearViscoelasticityBase::_apparent_elasticity_tensor
protected

Apparent elasticity tensor. This is NOT the elasticity tensor of the material.

Definition at line 196 of file LinearViscoelasticityBase.h.

◆ _apparent_elasticity_tensor_inv

 MaterialProperty& LinearViscoelasticityBase::_apparent_elasticity_tensor_inv
protected

Inverse of the apparent elasticity tensor.

Definition at line 198 of file LinearViscoelasticityBase.h.

◆ _base_name

 std::string ComputeElasticityTensorBase::_base_name
protectedinherited

◆ _components

 unsigned int LinearViscoelasticityBase::_components
protected

This is the number of internal variables required by the model.

This must be set in the constructor of an inherited class. See GeneralizedKelvinVoigtModel for example.

Definition at line 213 of file LinearViscoelasticityBase.h.

◆ _creep_strain_old

 const MaterialProperty& LinearViscoelasticityBase::_creep_strain_old
protected

Previous value of the true creep strain for update purposes.

This is calculated by a ComputeLinearViscoelasticStress or a LinearViscoelasticStressUpdate material.

Definition at line 251 of file LinearViscoelasticityBase.h.

Referenced by GeneralizedMaxwellBase::updateQpViscousStrains().

◆ _dashpot_viscosities

 std::vector *> LinearViscoelasticityBase::_dashpot_viscosities
protected

List of viscosities of each subsequent dashpot in the chain.

Definition at line 227 of file LinearViscoelasticityBase.h.

◆ _dashpot_viscosities_old

 std::vector *> LinearViscoelasticityBase::_dashpot_viscosities_old
protected

Definition at line 228 of file LinearViscoelasticityBase.h.

◆ _driving_eigenstrain

 const MaterialProperty* LinearViscoelasticityBase::_driving_eigenstrain
protected

Pointer to the value of the driving eigenstrain.

Definition at line 258 of file LinearViscoelasticityBase.h.

◆ _driving_eigenstrain_name

 std::string LinearViscoelasticityBase::_driving_eigenstrain_name
protected

Name of the eigenstrain that drives the additional creep strain.

Definition at line 256 of file LinearViscoelasticityBase.h.

◆ _driving_eigenstrain_old

 const MaterialProperty* LinearViscoelasticityBase::_driving_eigenstrain_old
protected

Definition at line 259 of file LinearViscoelasticityBase.h.

◆ _elastic_strain_old

 const MaterialProperty& LinearViscoelasticityBase::_elastic_strain_old
protected

previous value of the elastic strain for update purposes

Definition at line 245 of file LinearViscoelasticityBase.h.

◆ _elasticity_tensor_inv

 MaterialProperty& LinearViscoelasticityBase::_elasticity_tensor_inv
protected

Instantaneous elasticity tensor. This IS the real elasticity tensor of the material.

Inverse of the instaneous elasticity tensor

Definition at line 203 of file LinearViscoelasticityBase.h.

◆ _elasticity_tensor_name

 std::string ComputeElasticityTensorBase::_elasticity_tensor_name
protectedinherited

◆ _first_elasticity_tensor

 MaterialProperty& LinearViscoelasticityBase::_first_elasticity_tensor
protected

Elasticity tensor of a stand-alone elastic spring in the chain.

Definition at line 216 of file LinearViscoelasticityBase.h.

◆ _first_elasticity_tensor_inv

 MaterialProperty* LinearViscoelasticityBase::_first_elasticity_tensor_inv
protected

Definition at line 217 of file LinearViscoelasticityBase.h.

Referenced by computeQpViscoelasticPropertiesInv().

◆ _force_recompute_properties

 bool LinearViscoelasticityBase::_force_recompute_properties
protected

If activated, the time-stepping scheme will be re-initialized at each step of the solver.

This may be required for models in which the mechanical properties vary following other variables. If the mechanical properties are constant through the time step, this can be set to false.

Definition at line 267 of file LinearViscoelasticityBase.h.

Referenced by computeQpElasticityTensor().

◆ _has_driving_eigenstrain

 bool LinearViscoelasticityBase::_has_driving_eigenstrain
protected

Indicates if the model is only driven by the stress, or also by an additional eigenstrain.

Definition at line 254 of file LinearViscoelasticityBase.h.

◆ _has_longterm_dashpot

 bool LinearViscoelasticityBase::_has_longterm_dashpot
protected

Indicates if the spring-dashpot assembly has a single dashpot not associated with a spring.

Definition at line 208 of file LinearViscoelasticityBase.h.

◆ _integration_rule

 IntegrationRule LinearViscoelasticityBase::_integration_rule
protected

Determines how theta is computed.

Definition at line 191 of file LinearViscoelasticityBase.h.

Referenced by computeTheta().

◆ _need_viscoelastic_properties_inverse

 bool LinearViscoelasticityBase::_need_viscoelastic_properties_inverse
protected

If active, indicates that we need to call computeQpViscoelasticPropertiesInv()

Definition at line 206 of file LinearViscoelasticityBase.h.

◆ _prefactor_function

 Function* const ComputeElasticityTensorBase::_prefactor_function
protectedinherited

prefactor function to multiply the elasticity tensor with

Definition at line 42 of file ComputeElasticityTensorBase.h.

◆ _springs_elasticity_tensors

 std::vector *> LinearViscoelasticityBase::_springs_elasticity_tensors
protected

◆ _springs_elasticity_tensors_inv

 std::vector *> LinearViscoelasticityBase::_springs_elasticity_tensors_inv
protected

◆ _springs_elasticity_tensors_inv_old

 std::vector *> LinearViscoelasticityBase::_springs_elasticity_tensors_inv_old
protected

Definition at line 223 of file LinearViscoelasticityBase.h.

◆ _step_zero

 bool& LinearViscoelasticityBase::_step_zero
protected

checks whether we are at the first time step

Definition at line 270 of file LinearViscoelasticityBase.h.

Referenced by recomputeQpApparentProperties().

◆ _theta

 Real LinearViscoelasticityBase::_theta
protected

User-defined value for theta.

Definition at line 193 of file LinearViscoelasticityBase.h.

Referenced by computeTheta(), and LinearViscoelasticityBase().

◆ _viscous_strains

 std::vector *> LinearViscoelasticityBase::_viscous_strains
protected

The internal strain variables required by the time-stepping procedure (must be on a one-on-one basis with the number of dashpot).

Definition at line 235 of file LinearViscoelasticityBase.h.

◆ _viscous_strains_old

 std::vector *> LinearViscoelasticityBase::_viscous_strains_old
protected

Definition at line 236 of file LinearViscoelasticityBase.h.

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