ADCompute2DIncrementalStrain defines a strain increment only for incremental strains in 2D geometries, handling the out of plane strains. More...

#include <ADCompute2DIncrementalStrain.h>

Inheritance diagram for ADCompute2DIncrementalStrain< compute_stage >:

[legend]

Public Member Functions
	ADCompute2DIncrementalStrain (const InputParameters &parameters)

void	initialSetup () override

virtual void	computeProperties () override

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
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 ADReal	computeOutOfPlaneGradDisp ()=0
	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...

virtual Real	computeOutOfPlaneGradDispOld ()=0
	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	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 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 ADCompute2DIncrementalStrain< compute_stage >

ADCompute2DIncrementalStrain defines a strain increment only for incremental strains in 2D geometries, handling the out of plane strains.

ADCompute2DIncrementalStrain contains a virtual method to define the out-of-plane strain as a general nonzero value in the inherited classes ComputePlaneIncrementalStrain and ComputeAxisymmetricRZIncrementalStrain.

Definition at line 21 of file ADCompute2DIncrementalStrain.h.

Constructor & Destructor Documentation

◆ ADCompute2DIncrementalStrain()

template<ComputeStage compute_stage>

ADCompute2DIncrementalStrain< compute_stage >::ADCompute2DIncrementalStrain ( const InputParameters & parameters )

Definition at line 30 of file ADCompute2DIncrementalStrain.C.

   : ADComputeIncrementalSmallStrain<compute_stage>(parameters),
     _out_of_plane_direction(getParam<MooseEnum>("out_of_plane_direction"))
 {
 }

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>

virtual ADReal ADCompute2DIncrementalStrain< compute_stage >::computeOutOfPlaneGradDisp ( )

protectedpure virtual

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.

Implemented in ADComputeAxisymmetricRZIncrementalStrain< compute_stage >.

◆ computeOutOfPlaneGradDispOld()

template<ComputeStage compute_stage>

virtual Real ADCompute2DIncrementalStrain< compute_stage >::computeOutOfPlaneGradDispOld ( )

protectedpure virtual

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.

Implemented in ADComputeAxisymmetricRZIncrementalStrain< compute_stage >.

◆ computeProperties()

template<ComputeStage compute_stage>

void ADComputeIncrementalSmallStrain< compute_stage >::computeProperties ( )

overridevirtualinherited

Definition at line 36 of file ADComputeIncrementalSmallStrain.C.

 {
   ADReal volumetric_strain = 0.0;
   for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
   {
     ADRankTwoTensor total_strain_increment;
     computeTotalStrainIncrement(total_strain_increment);
  
     _strain_increment[_qp] = total_strain_increment;
  
     if (_volumetric_locking_correction)
       volumetric_strain += total_strain_increment.trace() * _JxW[_qp] * _coord[_qp];
   }
  
   if (_volumetric_locking_correction)
     volumetric_strain /= _current_elem_volume;
  
   for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
   {
     if (_volumetric_locking_correction)
     {
       const auto correction = (volumetric_strain - _strain_increment[_qp].trace()) / 3.0;
       _strain_increment[_qp](0, 0) += correction;
       _strain_increment[_qp](1, 1) += correction;
       _strain_increment[_qp](2, 2) += correction;
     }
  
     _total_strain[_qp] = _strain_increment[_qp] + _total_strain_old[_qp];
  
     // Remove the Eigen strain increment
     subtractEigenstrainIncrementFromStrain(_strain_increment[_qp]);
  
     // strain rate
     if (_dt > 0)
       _strain_rate[_qp] = _strain_increment[_qp] / _dt;
     else
       _strain_rate[_qp].zero();
  
     // Update strain in intermediate configuration: rotations are not needed
     _mechanical_strain[_qp] = _strain_increment[_qp] + _mechanical_strain_old[_qp];
  
     // incremental small strain does not include rotation
     _rotation_increment[_qp].setToIdentity();
   }
  
   copyDualNumbersToValues();
 }

◆ computeTotalStrainIncrement()

template<ComputeStage compute_stage>

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

overrideprotectedvirtual

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.

 {
   // Deformation gradient calculation for 2D problems
   ADRankTwoTensor A(
       (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]); // Deformation gradient
   RankTwoTensor Fbar((*_grad_disp_old[0])[_qp],
                      (*_grad_disp_old[1])[_qp],
                      (*_grad_disp_old[2])[_qp]); // Old Deformation gradient
  
   // Compute the displacement gradient of the out of plane direction for plane strain,
   // generalized plane strain, or axisymmetric problems
   A(_out_of_plane_direction, _out_of_plane_direction) = computeOutOfPlaneGradDisp();
   Fbar(_out_of_plane_direction, _out_of_plane_direction) = computeOutOfPlaneGradDispOld();
  
   A -= Fbar; // very nearly A = gradU - gradUold
  
   total_strain_increment = 0.5 * (A + A.transpose());
 }

◆ displacementIntegrityCheck()

template<ComputeStage compute_stage>

void ADCompute2DIncrementalStrain< compute_stage >::displacementIntegrityCheck ( )

overrideprotectedvirtual

Reimplemented from ADComputeStrainBase< compute_stage >.

Definition at line 85 of file ADCompute2DIncrementalStrain.C.

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

◆ initialSetup()

template<ComputeStage compute_stage>

void ADCompute2DIncrementalStrain< compute_stage >::initialSetup ( )

override

Definition at line 39 of file ADCompute2DIncrementalStrain.C.

 {
   for (unsigned int i = 0; i < 3; ++i)
   {
     if (_out_of_plane_direction == i)
     {
       _disp[i] = &adZeroValue();
       _grad_disp[i] = &adZeroGradient();
     }
     else
     {
       _disp[i] = &adCoupledValue("displacements", i);
       _grad_disp[i] = &adCoupledGradient("displacements", i);
     }
  
     if (_fe_problem.isTransient() && i != _out_of_plane_direction)
       _grad_disp_old[i] = &coupledGradientOld("displacements", i);
     else
       _grad_disp_old[i] = &_grad_zero;
   }
 }

◆ initQpStatefulProperties()

template<ComputeStage compute_stage>

void ADComputeIncrementalStrainBase< compute_stage >::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Reimplemented from ADComputeStrainBase< compute_stage >.

Definition at line 55 of file ADComputeIncrementalStrainBase.C.

 {
   _mechanical_strain[_qp].zero();
   _total_strain[_qp].zero();
 }

◆ subtractEigenstrainIncrementFromStrain()

template<ComputeStage compute_stage>

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

protectedinherited

Definition at line 63 of file ADComputeIncrementalStrainBase.C.

 {
   for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
   {
     strain -= (*_eigenstrains[i])[_qp];
     strain += (*_eigenstrains_old[i])[_qp];
   }
 }

◆ validParams()

template<ComputeStage compute_stage>

InputParameters ADCompute2DIncrementalStrain< compute_stage >::validParams ( )

static

Definition at line 18 of file ADCompute2DIncrementalStrain.C.

 {
   InputParameters params = ADComputeIncrementalSmallStrain<compute_stage>::validParams();
   params.addClassDescription("Compute strain increment for incremental strains in 2D geometries.");
  
   MooseEnum outOfPlaneDirection("x y z", "z");
   params.addParam<MooseEnum>(
       "out_of_plane_direction", outOfPlaneDirection, "The direction of the out-of-plane strain.");
   return params;
 }

Referenced by ADComputeAxisymmetricRZIncrementalStrain< compute_stage >::validParams().

Member Data Documentation

◆ _base_name

template<ComputeStage compute_stage>

const std::string ADComputeStrainBase< compute_stage >::_base_name

protectedinherited

Definition at line 60 of file ADComputeStrainBase.h.

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

◆ _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.

◆ _eigenstrain_names

template<ComputeStage compute_stage>

std::vector<MaterialPropertyName> ADComputeStrainBase< compute_stage >::_eigenstrain_names

protectedinherited

Definition at line 65 of file ADComputeStrainBase.h.

Referenced by ADComputeIncrementalStrainBase< compute_stage >::ADComputeIncrementalStrainBase(), and ADComputeStrainBase< compute_stage >::ADComputeStrainBase().

◆ _eigenstrains

template<ComputeStage compute_stage>

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

protectedinherited

Definition at line 66 of file ADComputeStrainBase.h.

Referenced by ADComputeGreenLagrangeStrain< compute_stage >::ADComputeGreenLagrangeStrain(), and ADComputeStrainBase< compute_stage >::ADComputeStrainBase().

◆ _eigenstrains_old

template<ComputeStage compute_stage>

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

protectedinherited

Definition at line 58 of file ADComputeIncrementalStrainBase.h.

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

◆ _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

protected

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

Definition at line 70 of file ADComputeStrainBase.h.

Referenced by ADComputeGreenLagrangeStrain< compute_stage >::ADComputeGreenLagrangeStrain(), and ADComputeStrainBase< compute_stage >::ADComputeStrainBase().

◆ usingComputeIncrementalSmallStrainMembers

template<ComputeStage compute_stage>

ADCompute2DIncrementalStrain< compute_stage >::usingComputeIncrementalSmallStrainMembers

protected

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:

tensor_mechanics/include/materials/ADCompute2DIncrementalStrain.h
tensor_mechanics/src/materials/ADCompute2DIncrementalStrain.C

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

template<ComputeStage compute_stage> class ADCompute2DIncrementalStrain< compute_stage >

Constructor & Destructor Documentation

◆ ADCompute2DIncrementalStrain()

Member Function Documentation

◆ ADMaterialProperty() [1/4]

◆ ADMaterialProperty() [2/4]

◆ ADMaterialProperty() [3/4]

◆ ADMaterialProperty() [4/4]

◆ computeOutOfPlaneGradDisp()

◆ computeOutOfPlaneGradDispOld()

◆ computeProperties()

◆ computeTotalStrainIncrement()

◆ displacementIntegrityCheck()

◆ initialSetup()

◆ initQpStatefulProperties()

◆ subtractEigenstrainIncrementFromStrain()

◆ validParams()

Member Data Documentation

◆ _base_name

◆ _current_elem_volume

◆ _disp

◆ _eigenstrain_names

◆ _eigenstrains

◆ _eigenstrains_old

◆ _grad_disp

◆ _grad_disp_old

◆ _mechanical_strain_old

◆ _ndisp

◆ _out_of_plane_direction

◆ _total_strain_old

◆ _volumetric_locking_correction

◆ usingComputeIncrementalSmallStrainMembers

◆ usingComputeIncrementalStrainBaseMembers

◆ usingComputeStrainBaseMembers

◆ usingMaterialMembers

template<ComputeStage compute_stage>
class ADCompute2DIncrementalStrain< compute_stage >