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

ComputeStrainBase is the base class for strain tensors. More...

#include <ComputeStrainBase.h>

Inheritance diagram for ComputeStrainBase:
[legend]

Public Member Functions

 ComputeStrainBase (const InputParameters &parameters)
 
void initialSetup () override
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

virtual void initQpStatefulProperties () override
 
virtual void displacementIntegrityCheck ()
 

Protected Attributes

unsigned int _ndisp
 Coupled displacement variables. More...
 
std::vector< const VariableValue * > _disp
 
std::vector< const VariableGradient * > _grad_disp
 
const std::string _base_name
 
MaterialProperty< RankTwoTensor > & _mechanical_strain
 
MaterialProperty< RankTwoTensor > & _total_strain
 
std::vector< MaterialPropertyName > _eigenstrain_names
 
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains
 
const MaterialProperty< RankTwoTensor > * _global_strain
 
const bool _volumetric_locking_correction
 
const Real & _current_elem_volume
 

Detailed Description

ComputeStrainBase is the base class for strain tensors.

Definition at line 26 of file ComputeStrainBase.h.

Constructor & Destructor Documentation

◆ ComputeStrainBase()

ComputeStrainBase::ComputeStrainBase ( const InputParameters &  parameters)

Definition at line 38 of file ComputeStrainBase.C.

39  : DerivativeMaterialInterface<Material>(parameters),
40  _ndisp(coupledComponents("displacements")),
41  _disp(3),
42  _grad_disp(3),
43  _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
44  _mechanical_strain(declareProperty<RankTwoTensor>(_base_name + "mechanical_strain")),
45  _total_strain(declareProperty<RankTwoTensor>(_base_name + "total_strain")),
46  _eigenstrain_names(getParam<std::vector<MaterialPropertyName>>("eigenstrain_names")),
48  _global_strain(isParamValid("global_strain")
49  ? &getMaterialProperty<RankTwoTensor>(_base_name + "global_strain")
50  : nullptr),
51  _volumetric_locking_correction(getParam<bool>("volumetric_locking_correction") &&
52  !isBoundaryMaterial()),
53  _current_elem_volume(_assembly.elemVolume())
54 {
55  for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
56  {
58  _eigenstrains[i] = &getMaterialProperty<RankTwoTensor>(_eigenstrain_names[i]);
59  }
60 
62  paramError("volumetric_locking_correction", "has to be set to false for 1-D problems.");
63 
64  if (getParam<bool>("use_displaced_mesh"))
65  paramError("use_displaced_mesh", "The strain calculator needs to run on the undisplaced mesh.");
66 
67  // Generate warning when volumetric locking correction is used with second order elements
68  if (_mesh.hasSecondOrderElements() && _volumetric_locking_correction)
69  mooseWarning("Volumetric locking correction is not required for second order elements. Using "
70  "volumetric locking with second order elements could cause zigzag patterns in "
71  "stresses and strains.");
72 }

Member Function Documentation

◆ displacementIntegrityCheck()

void ComputeStrainBase::displacementIntegrityCheck ( )
protectedvirtual

Reimplemented in Compute2DFiniteStrain, Compute2DIncrementalStrain, and Compute2DSmallStrain.

Definition at line 94 of file ComputeStrainBase.C.

95 {
96  // Checking for consistency between mesh size and length of the provided displacements vector
97  if (_ndisp != _mesh.dimension())
98  paramError(
99  "displacements",
100  "The number of variables supplied in 'displacements' must match the mesh dimension.");
101 }

Referenced by initialSetup().

◆ initialSetup()

void ComputeStrainBase::initialSetup ( )
override

Definition at line 75 of file ComputeStrainBase.C.

76 {
78  // fetch coupled variables and gradients (as stateful properties if necessary)
79  for (unsigned int i = 0; i < _ndisp; ++i)
80  {
81  _disp[i] = &coupledValue("displacements", i);
82  _grad_disp[i] = &coupledGradient("displacements", i);
83  }
84 
85  // set unused dimensions to zero
86  for (unsigned i = _ndisp; i < 3; ++i)
87  {
88  _disp[i] = &_zero;
89  _grad_disp[i] = &_grad_zero;
90  }
91 }

Referenced by ComputeIncrementalStrainBase::initialSetup(), and ComputeAxisymmetric1DSmallStrain::initialSetup().

◆ initQpStatefulProperties()

void ComputeStrainBase::initQpStatefulProperties ( )
overrideprotectedvirtual

Reimplemented in ComputeIncrementalStrainBase, and ComputeCosseratIncrementalSmallStrain.

Definition at line 104 of file ComputeStrainBase.C.

105 {
106  _mechanical_strain[_qp].zero();
107  _total_strain[_qp].zero();
108 }

◆ validParams()

InputParameters ComputeStrainBase::validParams ( )
static

Definition at line 17 of file ComputeStrainBase.C.

18 {
19  InputParameters params = Material::validParams();
20  params.addRequiredCoupledVar(
21  "displacements",
22  "The displacements appropriate for the simulation geometry and coordinate system");
23  params.addParam<std::string>("base_name",
24  "Optional parameter that allows the user to define "
25  "multiple mechanics material systems on the same "
26  "block, i.e. for multiple phases");
27  params.addParam<bool>(
28  "volumetric_locking_correction", false, "Flag to correct volumetric locking");
29  params.addParam<std::vector<MaterialPropertyName>>(
30  "eigenstrain_names", "List of eigenstrains to be applied in this strain calculation");
31  params.addParam<MaterialPropertyName>("global_strain",
32  "Optional material property holding a global strain "
33  "tensor applied to the mesh as a whole");
34  params.suppressParameter<bool>("use_displaced_mesh");
35  return params;
36 }

Referenced by ComputeCosseratSmallStrain::validParams(), ComputeIncrementalStrainBase::validParams(), and ComputeSmallStrain::validParams().

Member Data Documentation

◆ _base_name

const std::string ComputeStrainBase::_base_name
protected

Definition at line 44 of file ComputeStrainBase.h.

Referenced by ComputeStrainBase().

◆ _current_elem_volume

const Real& ComputeStrainBase::_current_elem_volume
protected

◆ _disp

std::vector<const VariableValue *> ComputeStrainBase::_disp
protected

◆ _eigenstrain_names

std::vector<MaterialPropertyName> ComputeStrainBase::_eigenstrain_names
protected

◆ _eigenstrains

std::vector<const MaterialProperty<RankTwoTensor> *> ComputeStrainBase::_eigenstrains
protected

◆ _global_strain

const MaterialProperty<RankTwoTensor>* ComputeStrainBase::_global_strain
protected

◆ _grad_disp

std::vector<const VariableGradient *> ComputeStrainBase::_grad_disp
protected

◆ _mechanical_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBase::_mechanical_strain
protected

◆ _ndisp

unsigned int ComputeStrainBase::_ndisp
protected

◆ _total_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBase::_total_strain
protected

◆ _volumetric_locking_correction

const bool ComputeStrainBase::_volumetric_locking_correction
protected

The documentation for this class was generated from the following files:
ComputeStrainBase::displacementIntegrityCheck
virtual void displacementIntegrityCheck()
Definition: ComputeStrainBase.C:94
ComputeStrainBase::_current_elem_volume
const Real & _current_elem_volume
Definition: ComputeStrainBase.h:56
ComputeStrainBase::_disp
std::vector< const VariableValue * > _disp
Definition: ComputeStrainBase.h:41
ComputeStrainBase::_base_name
const std::string _base_name
Definition: ComputeStrainBase.h:44
ComputeStrainBase::_ndisp
unsigned int _ndisp
Coupled displacement variables.
Definition: ComputeStrainBase.h:40
ComputeStrainBase::_mechanical_strain
MaterialProperty< RankTwoTensor > & _mechanical_strain
Definition: ComputeStrainBase.h:46
validParams
InputParameters validParams()
ComputeStrainBase::_global_strain
const MaterialProperty< RankTwoTensor > * _global_strain
Definition: ComputeStrainBase.h:53
ComputeStrainBase::_eigenstrain_names
std::vector< MaterialPropertyName > _eigenstrain_names
Definition: ComputeStrainBase.h:50
ComputeStrainBase::_eigenstrains
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains
Definition: ComputeStrainBase.h:51
ComputeStrainBase::_total_strain
MaterialProperty< RankTwoTensor > & _total_strain
Definition: ComputeStrainBase.h:48
ComputeStrainBase::_volumetric_locking_correction
const bool _volumetric_locking_correction
Definition: ComputeStrainBase.h:55
ComputeStrainBase::_grad_disp
std::vector< const VariableGradient * > _grad_disp
Definition: ComputeStrainBase.h:42