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

ComputeAxisymmetricRZSmallStrain defines small strains in an Axisymmetric system. More...

#include <ComputeAxisymmetricRZSmallStrain.h>

Inheritance diagram for ComputeAxisymmetricRZSmallStrain:
[legend]

Public Member Functions

 ComputeAxisymmetricRZSmallStrain (const InputParameters &parameters)
 
virtual void initialSetup () override
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

virtual Real computeOutOfPlaneStrain () override
 
virtual void computeProperties () override
 
virtual void displacementIntegrityCheck () override
 
virtual void initQpStatefulProperties () override
 

Protected Attributes

const unsigned int _out_of_plane_direction
 
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

ComputeAxisymmetricRZSmallStrain defines small strains in an Axisymmetric system.

The COORD_TYPE in the Problem block must be set to RZ.

Definition at line 23 of file ComputeAxisymmetricRZSmallStrain.h.

Constructor & Destructor Documentation

◆ ComputeAxisymmetricRZSmallStrain()

ComputeAxisymmetricRZSmallStrain::ComputeAxisymmetricRZSmallStrain ( const InputParameters &  parameters)

Definition at line 26 of file ComputeAxisymmetricRZSmallStrain.C.

28  : Compute2DSmallStrain(parameters)
29 {
30 }

Member Function Documentation

◆ computeOutOfPlaneStrain()

Real ComputeAxisymmetricRZSmallStrain::computeOutOfPlaneStrain ( )
overrideprotectedvirtual

Implements Compute2DSmallStrain.

Definition at line 46 of file ComputeAxisymmetricRZSmallStrain.C.

47 {
48  if (!MooseUtils::absoluteFuzzyEqual(_q_point[_qp](0), 0.0))
49  return (*_disp[0])[_qp] / _q_point[_qp](0);
50  else
51  return 0.0;
52 }

◆ computeProperties()

void Compute2DSmallStrain::computeProperties ( )
overrideprotectedvirtualinherited

Reimplemented from ComputeSmallStrain.

Definition at line 54 of file Compute2DSmallStrain.C.

55 {
56  const auto o0 = _out_of_plane_direction;
57  const auto o1 = (_out_of_plane_direction + 1) % 3;
58  const auto o2 = (_out_of_plane_direction + 2) % 3;
59 
60  Real volumetric_strain = 0.0;
61  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
62  {
63  _total_strain[_qp](o0, o0) = computeOutOfPlaneStrain();
64  _total_strain[_qp](o1, o1) = (*_grad_disp[o1])[_qp](o1);
65  _total_strain[_qp](o2, o2) = (*_grad_disp[o2])[_qp](o2);
66  _total_strain[_qp](o1, o2) = ((*_grad_disp[o1])[_qp](o2) + (*_grad_disp[o2])[_qp](o1)) / 2.0;
67  _total_strain[_qp](o2, o1) = _total_strain[_qp](o1, o2); // force the symmetrical strain tensor
68 
70  volumetric_strain += _total_strain[_qp].trace() * _JxW[_qp] * _coord[_qp];
71  }
72 
74  volumetric_strain /= _current_elem_volume;
75 
76  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
77  {
79  {
80  const Real correction = (volumetric_strain - _total_strain[_qp].trace()) / 3.0;
81  _total_strain[_qp](0, 0) += correction;
82  _total_strain[_qp](1, 1) += correction;
83  _total_strain[_qp](2, 2) += correction;
84  }
85 
87 
88  // Remove the eigenstrains
89  for (const auto es : _eigenstrains)
90  _mechanical_strain[_qp] -= (*es)[_qp];
91  }
92 }

◆ displacementIntegrityCheck()

void Compute2DSmallStrain::displacementIntegrityCheck ( )
overrideprotectedvirtualinherited

Reimplemented from ComputeStrainBase.

Definition at line 95 of file Compute2DSmallStrain.C.

96 {
97  if (_out_of_plane_direction != 2 && _ndisp != 3)
98  mooseError("For 2D simulations where the out-of-plane direction is x or y the number of "
99  "supplied displacements must be three.");
100  else if (_out_of_plane_direction == 2 && _ndisp != 2)
101  mooseError("For 2D simulations where the out-of-plane direction is z the number of supplied "
102  "displacements must be two.");
103 }

Referenced by Compute2DSmallStrain::initialSetup().

◆ initialSetup()

void ComputeAxisymmetricRZSmallStrain::initialSetup ( )
overridevirtual

Definition at line 33 of file ComputeAxisymmetricRZSmallStrain.C.

34 {
36 
37  if (getBlockCoordSystem() != Moose::COORD_RZ)
38  mooseError("The coordinate system must be set to RZ for Axisymmetric geometries.");
39 
40  if (_out_of_plane_direction != 2)
41  paramError("out_of_plane_direction",
42  "The out-of-plane direction for axisymmetric systems is currently restricted to z");
43 }

◆ initQpStatefulProperties()

void ComputeStrainBase::initQpStatefulProperties ( )
overrideprotectedvirtualinherited

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 ComputeAxisymmetricRZSmallStrain::validParams ( )
static

Definition at line 19 of file ComputeAxisymmetricRZSmallStrain.C.

20 {
21  InputParameters params = Compute2DSmallStrain::validParams();
22  params.addClassDescription("Compute a small strain in an Axisymmetric geometry");
23  return params;
24 }

Member Data Documentation

◆ _base_name

const std::string ComputeStrainBase::_base_name
protectedinherited

Definition at line 44 of file ComputeStrainBase.h.

Referenced by ComputeStrainBase::ComputeStrainBase().

◆ _current_elem_volume

const Real& ComputeStrainBase::_current_elem_volume
protectedinherited

◆ _disp

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

◆ _eigenstrain_names

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

◆ _eigenstrains

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

◆ _global_strain

const MaterialProperty<RankTwoTensor>* ComputeStrainBase::_global_strain
protectedinherited

◆ _grad_disp

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

◆ _mechanical_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBase::_mechanical_strain
protectedinherited

◆ _ndisp

unsigned int ComputeStrainBase::_ndisp
protectedinherited

◆ _out_of_plane_direction

const unsigned int Compute2DSmallStrain::_out_of_plane_direction
protectedinherited

◆ _total_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBase::_total_strain
protectedinherited

◆ _volumetric_locking_correction

const bool ComputeStrainBase::_volumetric_locking_correction
protectedinherited

The documentation for this class was generated from the following files:
Compute2DSmallStrain::Compute2DSmallStrain
Compute2DSmallStrain(const InputParameters &parameters)
Definition: Compute2DSmallStrain.C:28
ComputeStrainBase::_current_elem_volume
const Real & _current_elem_volume
Definition: ComputeStrainBase.h:56
ComputeStrainBase::_disp
std::vector< const VariableValue * > _disp
Definition: ComputeStrainBase.h:41
Compute2DSmallStrain::computeOutOfPlaneStrain
virtual Real computeOutOfPlaneStrain()=0
Compute2DSmallStrain::_out_of_plane_direction
const unsigned int _out_of_plane_direction
Definition: Compute2DSmallStrain.h:39
ComputeStrainBase::_ndisp
unsigned int _ndisp
Coupled displacement variables.
Definition: ComputeStrainBase.h:40
Compute2DSmallStrain::validParams
static InputParameters validParams()
Definition: Compute2DSmallStrain.C:17
ComputeStrainBase::_mechanical_strain
MaterialProperty< RankTwoTensor > & _mechanical_strain
Definition: ComputeStrainBase.h:46
Compute2DSmallStrain::initialSetup
void initialSetup() override
Definition: Compute2DSmallStrain.C:35
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