www.mooseframework.org
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
ComputeEigenstrainFromInitialStress Class Reference

ComputeEigenstrain computes an Eigenstrain that results from an initial stress The initial stress is defined in terms of Functions, which may be multiplied by optional AuxVariables. More...

#include <ComputeEigenstrainFromInitialStress.h>

Inheritance diagram for ComputeEigenstrainFromInitialStress:
[legend]

Public Member Functions

 ComputeEigenstrainFromInitialStress (const InputParameters &parameters)
 

Protected Member Functions

virtual void computeQpEigenstrain () override
 Compute the eigenstrain and store in _eigenstrain. More...
 
virtual void initQpStatefulProperties ()
 
virtual void computeQpProperties ()
 
Real computeVolumetricStrainComponent (const Real volumetric_strain) const
 Helper function for models that compute the eigenstrain based on a volumetric strain. More...
 

Protected Attributes

const std::string _base_name
 base_name for elasticity tensor to use to convert stress to strain More...
 
const MaterialProperty< RankFourTensor > & _elasticity_tensor
 elasticity tensor used to convert stress to strain More...
 
const MaterialProperty< RankTwoTensor > & _eigenstrain_old
 Stores the total eigenstrain in the previous step. More...
 
const bool _ini_aux_provided
 Whether the user has supplied AuxVariables representing the initial stress. More...
 
std::vector< Function * > _initial_stress_fcn
 initial stress components More...
 
std::vector< const VariableValue * > _ini_aux
 AuxVariables defining the initial stress. More...
 
std::string _eigenstrain_name
 Material property name for the eigenstrain tensor. More...
 
MaterialProperty< RankTwoTensor > & _eigenstrain
 Stores the current total eigenstrain. More...
 
bool & _step_zero
 Restartable data to check for the zeroth and first time steps for thermal calculations. More...
 

Detailed Description

ComputeEigenstrain computes an Eigenstrain that results from an initial stress The initial stress is defined in terms of Functions, which may be multiplied by optional AuxVariables.

Definition at line 26 of file ComputeEigenstrainFromInitialStress.h.

Constructor & Destructor Documentation

◆ ComputeEigenstrainFromInitialStress()

ComputeEigenstrainFromInitialStress::ComputeEigenstrainFromInitialStress ( const InputParameters &  parameters)

Definition at line 41 of file ComputeEigenstrainFromInitialStress.C.

43  : ComputeEigenstrainBase(parameters),
44  _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
45  _elasticity_tensor(getMaterialPropertyByName<RankFourTensor>(_base_name + "elasticity_tensor")),
46  _eigenstrain_old(getMaterialPropertyOld<RankTwoTensor>(_eigenstrain_name)),
47  _ini_aux_provided(isParamValid("initial_stress_aux"))
48 {
49  const std::vector<FunctionName> & fcn_names(
50  getParam<std::vector<FunctionName>>("initial_stress"));
51  const std::size_t num = fcn_names.size();
52 
53  if (num != LIBMESH_DIM * LIBMESH_DIM)
54  paramError(
55  "initial_stress",
56  "ComputeEigenstrainFromInitialStress: " + Moose::stringify(LIBMESH_DIM * LIBMESH_DIM) +
57  " initial stress functions must be provided. You supplied " + Moose::stringify(num) +
58  "\n");
59 
60  _initial_stress_fcn.resize(num);
61  for (unsigned i = 0; i < num; ++i)
62  _initial_stress_fcn[i] = &getFunctionByName(fcn_names[i]);
63 
65  {
66  const std::size_t aux_size = coupledComponents("initial_stress_aux");
67  if (aux_size != LIBMESH_DIM * LIBMESH_DIM)
68  paramError("initial_stress_aux",
69  "ComputeEigenstrainFromInitialStress: If you supply initial_stress_aux, " +
70  Moose::stringify(LIBMESH_DIM * LIBMESH_DIM) +
71  " values must be given. You supplied " + Moose::stringify(aux_size) + "\n");
72  _ini_aux.resize(0);
73  for (unsigned i = 0; i < aux_size; ++i)
74  _ini_aux.push_back(&coupledValue("initial_stress_aux", i));
75  }
76 }
std::vector< Function * > _initial_stress_fcn
initial stress components
const bool _ini_aux_provided
Whether the user has supplied AuxVariables representing the initial stress.
std::string _eigenstrain_name
Material property name for the eigenstrain tensor.
std::vector< const VariableValue * > _ini_aux
AuxVariables defining the initial stress.
const std::string _base_name
base_name for elasticity tensor to use to convert stress to strain
ComputeEigenstrainBase(const InputParameters &parameters)
const MaterialProperty< RankTwoTensor > & _eigenstrain_old
Stores the total eigenstrain in the previous step.
const MaterialProperty< RankFourTensor > & _elasticity_tensor
elasticity tensor used to convert stress to strain

Member Function Documentation

◆ computeQpEigenstrain()

void ComputeEigenstrainFromInitialStress::computeQpEigenstrain ( )
overrideprotectedvirtual

Compute the eigenstrain and store in _eigenstrain.

Implements ComputeEigenstrainBase.

Definition at line 79 of file ComputeEigenstrainFromInitialStress.C.

80 {
81  if (_t_step == 1)
82  {
83  RankTwoTensor initial_stress;
84  for (unsigned i = 0; i < LIBMESH_DIM; ++i)
85  for (unsigned j = 0; j < LIBMESH_DIM; ++j)
86  {
87  initial_stress(i, j) = _initial_stress_fcn[i * LIBMESH_DIM + j]->value(_t, _q_point[_qp]);
89  initial_stress(i, j) *= (*_ini_aux[i * LIBMESH_DIM + j])[_qp];
90  }
91 
92  _eigenstrain[_qp] = -_elasticity_tensor[_qp].invSymm() * initial_stress;
93  }
94  else
95  _eigenstrain[_qp] = _eigenstrain_old[_qp];
96 }
std::vector< Function * > _initial_stress_fcn
initial stress components
const bool _ini_aux_provided
Whether the user has supplied AuxVariables representing the initial stress.
std::vector< const VariableValue * > _ini_aux
AuxVariables defining the initial stress.
MaterialProperty< RankTwoTensor > & _eigenstrain
Stores the current total eigenstrain.
const MaterialProperty< RankTwoTensor > & _eigenstrain_old
Stores the total eigenstrain in the previous step.
const MaterialProperty< RankFourTensor > & _elasticity_tensor
elasticity tensor used to convert stress to strain

◆ computeQpProperties()

void ComputeEigenstrainBase::computeQpProperties ( )
protectedvirtualinherited

Definition at line 53 of file ComputeEigenstrainBase.C.

54 {
55  if (_t_step >= 1)
56  _step_zero = false;
57 
58  // Skip the eigenstrain calculation in step zero because no solution is computed during
59  // the zeroth step, hence computing the eigenstrain in the zeroth step would result in
60  // an incorrect calculation of mechanical_strain, which is stateful.
61  if (!_step_zero)
63 }
virtual void computeQpEigenstrain()=0
Compute the eigenstrain and store in _eigenstrain.
bool & _step_zero
Restartable data to check for the zeroth and first time steps for thermal calculations.

◆ computeVolumetricStrainComponent()

Real ComputeEigenstrainBase::computeVolumetricStrainComponent ( const Real  volumetric_strain) const
protectedinherited

Helper function for models that compute the eigenstrain based on a volumetric strain.

This function computes the diagonal components of the eigenstrain tensor as logarithmic strains.

Parameters
volumetric_strainThe current volumetric strain to be applied
Returns
Current strain in one direction due to volumetric strain, expressed as a logarithmic strain

Definition at line 66 of file ComputeEigenstrainBase.C.

67 {
68  // The engineering strain in a given direction is:
69  // epsilon_eng = cbrt(volumetric_strain + 1.0) - 1.0
70  //
71  // We need to provide this as a logarithmic strain to be consistent with the strain measure
72  // used for finite strain:
73  // epsilon_log = log(1.0 + epsilon_eng)
74  //
75  // This can be simplified down to a more direct form:
76  // epsilon_log = log(cbrt(volumetric_strain + 1.0))
77  // or:
78  // epsilon_log = (1/3) log(volumetric_strain + 1.0)
79 
80  return std::log(volumetric_strain + 1.0) / 3.0;
81 }

◆ initQpStatefulProperties()

void ComputeEigenstrainBase::initQpStatefulProperties ( )
protectedvirtualinherited

Reimplemented in ComputeReducedOrderEigenstrain.

Definition at line 45 of file ComputeEigenstrainBase.C.

46 {
47  // This property can be promoted to be stateful by other models that use it,
48  // so it needs to be initalized.
49  _eigenstrain[_qp].zero();
50 }
MaterialProperty< RankTwoTensor > & _eigenstrain
Stores the current total eigenstrain.

Member Data Documentation

◆ _base_name

const std::string ComputeEigenstrainFromInitialStress::_base_name
protected

base_name for elasticity tensor to use to convert stress to strain

Definition at line 35 of file ComputeEigenstrainFromInitialStress.h.

◆ _eigenstrain

MaterialProperty<RankTwoTensor>& ComputeEigenstrainBase::_eigenstrain
protectedinherited

◆ _eigenstrain_name

std::string ComputeEigenstrainBase::_eigenstrain_name
protectedinherited

Material property name for the eigenstrain tensor.

Definition at line 43 of file ComputeEigenstrainBase.h.

◆ _eigenstrain_old

const MaterialProperty<RankTwoTensor>& ComputeEigenstrainFromInitialStress::_eigenstrain_old
protected

Stores the total eigenstrain in the previous step.

Definition at line 41 of file ComputeEigenstrainFromInitialStress.h.

Referenced by computeQpEigenstrain().

◆ _elasticity_tensor

const MaterialProperty<RankFourTensor>& ComputeEigenstrainFromInitialStress::_elasticity_tensor
protected

elasticity tensor used to convert stress to strain

Definition at line 38 of file ComputeEigenstrainFromInitialStress.h.

Referenced by computeQpEigenstrain().

◆ _ini_aux

std::vector<const VariableValue *> ComputeEigenstrainFromInitialStress::_ini_aux
protected

AuxVariables defining the initial stress.

Definition at line 50 of file ComputeEigenstrainFromInitialStress.h.

Referenced by ComputeEigenstrainFromInitialStress(), and computeQpEigenstrain().

◆ _ini_aux_provided

const bool ComputeEigenstrainFromInitialStress::_ini_aux_provided
protected

Whether the user has supplied AuxVariables representing the initial stress.

Definition at line 44 of file ComputeEigenstrainFromInitialStress.h.

Referenced by ComputeEigenstrainFromInitialStress(), and computeQpEigenstrain().

◆ _initial_stress_fcn

std::vector<Function *> ComputeEigenstrainFromInitialStress::_initial_stress_fcn
protected

initial stress components

Definition at line 47 of file ComputeEigenstrainFromInitialStress.h.

Referenced by ComputeEigenstrainFromInitialStress(), and computeQpEigenstrain().

◆ _step_zero

bool& ComputeEigenstrainBase::_step_zero
protectedinherited

Restartable data to check for the zeroth and first time steps for thermal calculations.

Definition at line 59 of file ComputeEigenstrainBase.h.

Referenced by ComputeEigenstrainBase::computeQpProperties().


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