www.mooseframework.org
Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
ADComputeAxisymmetricRZIncrementalStrain< compute_stage > Class Template Reference

ADComputeAxisymmetricRZIncrementalStrain defines a strain increment only for incremental strains in an Axisymmetric simulation. More...

#include <ADComputeAxisymmetricRZIncrementalStrain.h>

Inheritance diagram for ADComputeAxisymmetricRZIncrementalStrain< compute_stage >:
[legend]

Public Member Functions

 ADComputeAxisymmetricRZIncrementalStrain (const InputParameters &parameters)
 
void initialSetup () override
 
virtual void computeProperties () override
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

ADReal computeOutOfPlaneGradDisp () override
 Computes the current out-of-plane component of the displacement gradient; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes. More...
 
Real computeOutOfPlaneGradDispOld () override
 Computes the old out-of-plane component of the displacement gradient; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes. More...
 
virtual void displacementIntegrityCheck () override
 
virtual void computeTotalStrainIncrement (ADRankTwoTensor &total_strain_increment) override
 Computes the current and old deformation gradients with the assumptions for 2D geometries, including plane strain, generalized plane strain, and axisymmetric, and returns the total strain increment tensor. More...
 
virtual void initQpStatefulProperties () override
 
void subtractEigenstrainIncrementFromStrain (ADRankTwoTensor &strain)
 
 ADMaterialProperty (RankTwoTensor) &_strain_rate
 
 ADMaterialProperty (RankTwoTensor) &_strain_increment
 
 ADMaterialProperty (RankTwoTensor) &_rotation_increment
 
const ADMaterialProperty (RankTwoTensor) *_global_strain
 

Protected Attributes

const VariableValue & _disp_old_0
 the old value of the first component of the displacements vector More...
 
 usingCompute2DIncrementalStrainMembers
 
const unsigned int _out_of_plane_direction
 
 usingComputeIncrementalSmallStrainMembers
 
 usingComputeIncrementalStrainBaseMembers
 
std::vector< const VariableGradient * > _grad_disp_old
 
const MaterialProperty< RankTwoTensor > & _mechanical_strain_old
 
const MaterialProperty< RankTwoTensor > & _total_strain_old
 
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains_old
 
 usingComputeStrainBaseMembers
 
const unsigned int _ndisp
 Coupled displacement variables. More...
 
std::vector< const ADVariableValue * > _disp
 
std::vector< const ADVariableGradient * > _grad_disp
 
const std::string _base_name
 
std::vector< MaterialPropertyName > _eigenstrain_names
 
std::vector< const ADMaterialProperty(RankTwoTensor) * > _eigenstrains
 
const bool _volumetric_locking_correction
 
const Real & _current_elem_volume
 
 usingMaterialMembers
 

Detailed Description

template<ComputeStage compute_stage>
class ADComputeAxisymmetricRZIncrementalStrain< compute_stage >

ADComputeAxisymmetricRZIncrementalStrain defines a strain increment only for incremental strains in an Axisymmetric simulation.

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

Definition at line 15 of file ADComputeAxisymmetricRZIncrementalStrain.h.

Constructor & Destructor Documentation

◆ ADComputeAxisymmetricRZIncrementalStrain()

template<ComputeStage compute_stage>
ADComputeAxisymmetricRZIncrementalStrain< compute_stage >::ADComputeAxisymmetricRZIncrementalStrain ( const InputParameters &  parameters)

Definition at line 30 of file ADComputeAxisymmetricRZIncrementalStrain.C.

33  _disp_old_0(coupledValueOld("displacements", 0))
34 {
35 }

Member Function Documentation

◆ ADMaterialProperty() [1/4]

template<ComputeStage compute_stage>
ADComputeIncrementalStrainBase< compute_stage >::ADMaterialProperty ( RankTwoTensor  ) &
protectedinherited

◆ ADMaterialProperty() [2/4]

template<ComputeStage compute_stage>
ADComputeIncrementalStrainBase< compute_stage >::ADMaterialProperty ( RankTwoTensor  ) &
protectedinherited

◆ ADMaterialProperty() [3/4]

template<ComputeStage compute_stage>
ADComputeIncrementalStrainBase< compute_stage >::ADMaterialProperty ( RankTwoTensor  ) &
protectedinherited

◆ ADMaterialProperty() [4/4]

template<ComputeStage compute_stage>
const ADComputeStrainBase< compute_stage >::ADMaterialProperty ( RankTwoTensor  )
protectedinherited

◆ computeOutOfPlaneGradDisp()

template<ComputeStage compute_stage>
ADReal ADComputeAxisymmetricRZIncrementalStrain< compute_stage >::computeOutOfPlaneGradDisp ( )
overrideprotectedvirtual

Computes the current out-of-plane component of the displacement gradient; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes.

Implements ADCompute2DIncrementalStrain< compute_stage >.

Definition at line 53 of file ADComputeAxisymmetricRZIncrementalStrain.C.

54 {
55  if (!MooseUtils::absoluteFuzzyEqual(_q_point[_qp](0), 0.0))
56  return (*_disp[0])[_qp] / _q_point[_qp](0);
57  else
58  return 0.0;
59 }

◆ computeOutOfPlaneGradDispOld()

template<ComputeStage compute_stage>
Real ADComputeAxisymmetricRZIncrementalStrain< compute_stage >::computeOutOfPlaneGradDispOld ( )
overrideprotectedvirtual

Computes the old out-of-plane component of the displacement gradient; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes.

Implements ADCompute2DIncrementalStrain< compute_stage >.

Definition at line 63 of file ADComputeAxisymmetricRZIncrementalStrain.C.

64 {
65  if (!MooseUtils::absoluteFuzzyEqual(_q_point[_qp](0), 0.0))
66  return _disp_old_0[_qp] / _q_point[_qp](0);
67  else
68  return 0.0;
69 }

◆ computeProperties()

template<ComputeStage compute_stage>
void ADComputeIncrementalSmallStrain< compute_stage >::computeProperties ( )
overridevirtualinherited

Definition at line 36 of file ADComputeIncrementalSmallStrain.C.

37 {
38  ADReal volumetric_strain = 0.0;
39  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
40  {
41  ADRankTwoTensor total_strain_increment;
42  computeTotalStrainIncrement(total_strain_increment);
43 
44  _strain_increment[_qp] = total_strain_increment;
45 
47  volumetric_strain += total_strain_increment.trace() * _JxW[_qp] * _coord[_qp];
48  }
49 
51  volumetric_strain /= _current_elem_volume;
52 
53  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
54  {
56  {
57  const auto correction = (volumetric_strain - _strain_increment[_qp].trace()) / 3.0;
58  _strain_increment[_qp](0, 0) += correction;
59  _strain_increment[_qp](1, 1) += correction;
60  _strain_increment[_qp](2, 2) += correction;
61  }
62 
63  _total_strain[_qp] = _strain_increment[_qp] + _total_strain_old[_qp];
64 
65  // Remove the Eigen strain increment
66  subtractEigenstrainIncrementFromStrain(_strain_increment[_qp]);
67 
68  // strain rate
69  if (_dt > 0)
70  _strain_rate[_qp] = _strain_increment[_qp] / _dt;
71  else
72  _strain_rate[_qp].zero();
73 
74  // Update strain in intermediate configuration: rotations are not needed
75  _mechanical_strain[_qp] = _strain_increment[_qp] + _mechanical_strain_old[_qp];
76 
77  // incremental small strain does not include rotation
78  _rotation_increment[_qp].setToIdentity();
79  }
80 
81  copyDualNumbersToValues();
82 }

◆ computeTotalStrainIncrement()

template<ComputeStage compute_stage>
void ADCompute2DIncrementalStrain< compute_stage >::computeTotalStrainIncrement ( ADRankTwoTensor &  total_strain_increment)
overrideprotectedvirtualinherited

Computes the current and old deformation gradients with the assumptions for 2D geometries, including plane strain, generalized plane strain, and axisymmetric, and returns the total strain increment tensor.

Reimplemented from ADComputeIncrementalSmallStrain< compute_stage >.

Definition at line 63 of file ADCompute2DIncrementalStrain.C.

65 {
66  // Deformation gradient calculation for 2D problems
67  ADRankTwoTensor A(
68  (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]); // Deformation gradient
69  RankTwoTensor Fbar((*_grad_disp_old[0])[_qp],
70  (*_grad_disp_old[1])[_qp],
71  (*_grad_disp_old[2])[_qp]); // Old Deformation gradient
72 
73  // Compute the displacement gradient of the out of plane direction for plane strain,
74  // generalized plane strain, or axisymmetric problems
77 
78  A -= Fbar; // very nearly A = gradU - gradUold
79 
80  total_strain_increment = 0.5 * (A + A.transpose());
81 }

◆ displacementIntegrityCheck()

template<ComputeStage compute_stage>
void ADCompute2DIncrementalStrain< compute_stage >::displacementIntegrityCheck ( )
overrideprotectedvirtualinherited

Reimplemented from ADComputeStrainBase< compute_stage >.

Definition at line 85 of file ADCompute2DIncrementalStrain.C.

86 {
87  if (_out_of_plane_direction != 2 && _ndisp != 3)
88  mooseError("For 2D simulations where the out-of-plane direction is x or y the number of "
89  "supplied displacements must be three.");
90  else if (_out_of_plane_direction == 2 && _ndisp != 2)
91  mooseError("For 2D simulations where the out-of-plane direction is z the number of supplied "
92  "displacements must be two.");
93 }

◆ initialSetup()

template<ComputeStage compute_stage>
void ADComputeAxisymmetricRZIncrementalStrain< compute_stage >::initialSetup ( )
override

Definition at line 39 of file ADComputeAxisymmetricRZIncrementalStrain.C.

40 {
42 
43  if (getBlockCoordSystem() != Moose::COORD_RZ)
44  mooseError("The coordinate system must be set to RZ for Axisymmetric geometries.");
45 
46  if (_out_of_plane_direction != 2)
47  paramError("out_of_plane_direction",
48  "The out-of-plane direction for axisymmetric systems is currently restricted to z");
49 }

◆ initQpStatefulProperties()

template<ComputeStage compute_stage>
void ADComputeIncrementalStrainBase< compute_stage >::initQpStatefulProperties ( )
overrideprotectedvirtualinherited

Reimplemented from ADComputeStrainBase< compute_stage >.

Definition at line 55 of file ADComputeIncrementalStrainBase.C.

56 {
57  _mechanical_strain[_qp].zero();
58  _total_strain[_qp].zero();
59 }

◆ subtractEigenstrainIncrementFromStrain()

template<ComputeStage compute_stage>
void ADComputeIncrementalStrainBase< compute_stage >::subtractEigenstrainIncrementFromStrain ( ADRankTwoTensor &  strain)
protectedinherited

Definition at line 63 of file ADComputeIncrementalStrainBase.C.

65 {
66  for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
67  {
68  strain -= (*_eigenstrains[i])[_qp];
69  strain += (*_eigenstrains_old[i])[_qp];
70  }
71 }

◆ validParams()

template<ComputeStage compute_stage>
InputParameters ADComputeAxisymmetricRZIncrementalStrain< compute_stage >::validParams ( )
static

Definition at line 21 of file ADComputeAxisymmetricRZIncrementalStrain.C.

22 {
24  params.addClassDescription("Compute a strain increment and rotation increment for finite strains "
25  "under axisymmetric assumptions.");
26  return params;
27 }

Member Data Documentation

◆ _base_name

template<ComputeStage compute_stage>
const std::string ADComputeStrainBase< compute_stage >::_base_name
protectedinherited

◆ _current_elem_volume

template<ComputeStage compute_stage>
const Real& ADComputeStrainBase< compute_stage >::_current_elem_volume
protectedinherited

Definition at line 71 of file ADComputeStrainBase.h.

◆ _disp

template<ComputeStage compute_stage>
std::vector<const ADVariableValue *> ADComputeStrainBase< compute_stage >::_disp
protectedinherited

Definition at line 57 of file ADComputeStrainBase.h.

◆ _disp_old_0

template<ComputeStage compute_stage>
const VariableValue& ADComputeAxisymmetricRZIncrementalStrain< compute_stage >::_disp_old_0
protected

the old value of the first component of the displacements vector

Definition at line 40 of file ADComputeAxisymmetricRZIncrementalStrain.h.

◆ _eigenstrain_names

template<ComputeStage compute_stage>
std::vector<MaterialPropertyName> ADComputeStrainBase< compute_stage >::_eigenstrain_names
protectedinherited

◆ _eigenstrains

template<ComputeStage compute_stage>
std::vector<const ADMaterialProperty(RankTwoTensor) *> ADComputeStrainBase< compute_stage >::_eigenstrains
protectedinherited

◆ _eigenstrains_old

template<ComputeStage compute_stage>
std::vector<const MaterialProperty<RankTwoTensor> *> ADComputeIncrementalStrainBase< compute_stage >::_eigenstrains_old
protectedinherited

◆ _grad_disp

template<ComputeStage compute_stage>
std::vector<const ADVariableGradient *> ADComputeStrainBase< compute_stage >::_grad_disp
protectedinherited

Definition at line 58 of file ADComputeStrainBase.h.

◆ _grad_disp_old

template<ComputeStage compute_stage>
std::vector<const VariableGradient *> ADComputeIncrementalStrainBase< compute_stage >::_grad_disp_old
protectedinherited

Definition at line 49 of file ADComputeIncrementalStrainBase.h.

◆ _mechanical_strain_old

template<ComputeStage compute_stage>
const MaterialProperty<RankTwoTensor>& ADComputeIncrementalStrainBase< compute_stage >::_mechanical_strain_old
protectedinherited

Definition at line 55 of file ADComputeIncrementalStrainBase.h.

◆ _ndisp

template<ComputeStage compute_stage>
const unsigned int ADComputeStrainBase< compute_stage >::_ndisp
protectedinherited

Coupled displacement variables.

Definition at line 56 of file ADComputeStrainBase.h.

Referenced by ADComputeStrainBase< compute_stage >::ADComputeStrainBase().

◆ _out_of_plane_direction

template<ComputeStage compute_stage>
const unsigned int ADCompute2DIncrementalStrain< compute_stage >::_out_of_plane_direction
protectedinherited

Definition at line 64 of file ADCompute2DIncrementalStrain.h.

◆ _total_strain_old

template<ComputeStage compute_stage>
const MaterialProperty<RankTwoTensor>& ADComputeIncrementalStrainBase< compute_stage >::_total_strain_old
protectedinherited

Definition at line 56 of file ADComputeIncrementalStrainBase.h.

◆ _volumetric_locking_correction

template<ComputeStage compute_stage>
const bool ADComputeStrainBase< compute_stage >::_volumetric_locking_correction
protectedinherited

◆ usingCompute2DIncrementalStrainMembers

template<ComputeStage compute_stage>
ADComputeAxisymmetricRZIncrementalStrain< compute_stage >::usingCompute2DIncrementalStrainMembers
protected

Definition at line 42 of file ADComputeAxisymmetricRZIncrementalStrain.h.

◆ usingComputeIncrementalSmallStrainMembers

template<ComputeStage compute_stage>
ADCompute2DIncrementalStrain< compute_stage >::usingComputeIncrementalSmallStrainMembers
protectedinherited

Definition at line 66 of file ADCompute2DIncrementalStrain.h.

◆ usingComputeIncrementalStrainBaseMembers

template<ComputeStage compute_stage>
ADComputeIncrementalSmallStrain< compute_stage >::usingComputeIncrementalStrainBaseMembers
protectedinherited

Definition at line 44 of file ADComputeIncrementalSmallStrain.h.

◆ usingComputeStrainBaseMembers

template<ComputeStage compute_stage>
ADComputeIncrementalStrainBase< compute_stage >::usingComputeStrainBaseMembers
protectedinherited

Definition at line 60 of file ADComputeIncrementalStrainBase.h.

◆ usingMaterialMembers

template<ComputeStage compute_stage>
ADComputeStrainBase< compute_stage >::usingMaterialMembers
protectedinherited

Definition at line 73 of file ADComputeStrainBase.h.


The documentation for this class was generated from the following files:
ADCompute2DIncrementalStrain
ADCompute2DIncrementalStrain defines a strain increment only for incremental strains in 2D geometries...
Definition: ADCompute2DIncrementalStrain.h:21
ADComputeIncrementalStrainBase::_total_strain_old
const MaterialProperty< RankTwoTensor > & _total_strain_old
Definition: ADComputeIncrementalStrainBase.h:56
ADComputeIncrementalSmallStrain::computeTotalStrainIncrement
virtual void computeTotalStrainIncrement(ADRankTwoTensor &total_strain_increment)
Computes the current and old deformation gradients and passes back the total strain increment tensor.
Definition: ADComputeIncrementalSmallStrain.C:86
ADComputeIncrementalStrainBase::subtractEigenstrainIncrementFromStrain
void subtractEigenstrainIncrementFromStrain(ADRankTwoTensor &strain)
Definition: ADComputeIncrementalStrainBase.C:63
ADComputeIncrementalStrainBase::_mechanical_strain_old
const MaterialProperty< RankTwoTensor > & _mechanical_strain_old
Definition: ADComputeIncrementalStrainBase.h:55
ADCompute2DIncrementalStrain::computeOutOfPlaneGradDisp
virtual ADReal computeOutOfPlaneGradDisp()=0
Computes the current out-of-plane component of the displacement gradient; as a virtual function,...
ADComputeStrainBase::_current_elem_volume
const Real & _current_elem_volume
Definition: ADComputeStrainBase.h:71
ADComputeStrainBase::_eigenstrains
std::vector< const ADMaterialProperty(RankTwoTensor) * > _eigenstrains
Definition: ADComputeStrainBase.h:66
ADCompute2DIncrementalStrain::_out_of_plane_direction
const unsigned int _out_of_plane_direction
Definition: ADCompute2DIncrementalStrain.h:64
ADCompute2DIncrementalStrain::validParams
static InputParameters validParams()
Definition: ADCompute2DIncrementalStrain.C:18
ADComputeStrainBase::_disp
std::vector< const ADVariableValue * > _disp
Definition: ADComputeStrainBase.h:57
ADComputeIncrementalStrainBase::_eigenstrains_old
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains_old
Definition: ADComputeIncrementalStrainBase.h:58
ADCompute2DIncrementalStrain::computeOutOfPlaneGradDispOld
virtual Real computeOutOfPlaneGradDispOld()=0
Computes the old out-of-plane component of the displacement gradient; as a virtual function,...
ADComputeIncrementalStrainBase::initialSetup
void initialSetup() override
Definition: ADComputeIncrementalStrainBase.C:41
ADComputeStrainBase::_grad_disp
std::vector< const ADVariableGradient * > _grad_disp
Definition: ADComputeStrainBase.h:58
RankTwoTensorTempl< Real >
ADComputeIncrementalStrainBase::_grad_disp_old
std::vector< const VariableGradient * > _grad_disp_old
Definition: ADComputeIncrementalStrainBase.h:49
ADComputeStrainBase::_ndisp
const unsigned int _ndisp
Coupled displacement variables.
Definition: ADComputeStrainBase.h:56
ADComputeAxisymmetricRZIncrementalStrain::_disp_old_0
const VariableValue & _disp_old_0
the old value of the first component of the displacements vector
Definition: ADComputeAxisymmetricRZIncrementalStrain.h:40
ADComputeStrainBase::_volumetric_locking_correction
const bool _volumetric_locking_correction
Definition: ADComputeStrainBase.h:70