ADComputeRSphericalIncrementalStrain defines a strain increment only for small strains in 1D spherical symmetry geometries. More...

#include <ADComputeRSphericalIncrementalStrain.h>

Inheritance diagram for ADComputeRSphericalIncrementalStrain< compute_stage >:

[legend]

Public Member Functions
	ADComputeRSphericalIncrementalStrain (const InputParameters &parameters)

virtual void	initialSetup () override

virtual void	computeProperties () override

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	computeTotalStrainIncrement (ADRankTwoTensor &total_strain_increment) override
	Computes the current and old deformation gradients with the assumptions for 1D spherical symmetry geometries: \( \epsilon_{\theta} = \epsilon_{\phi} = \frac{u_r}{r} \). 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

virtual void	displacementIntegrityCheck ()

Protected Attributes
const VariableValue &	_disp_old_0
	the old value of the first component of the displacements vector More...

	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 ADComputeRSphericalIncrementalStrain< compute_stage >

ADComputeRSphericalIncrementalStrain defines a strain increment only for small strains in 1D spherical symmetry geometries.

The strains in the polar and azimuthal directions are functions of the radial displacement.

Definition at line 15 of file ADComputeRSphericalIncrementalStrain.h.

Constructor & Destructor Documentation

◆ ADComputeRSphericalIncrementalStrain()

template<ComputeStage compute_stage>

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

Definition at line 32 of file ADComputeRSphericalIncrementalStrain.C.

   : ADComputeIncrementalSmallStrain<compute_stage>(parameters),
     _disp_old_0(coupledValueOld("displacements", 0))
 {
 }

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

◆ 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 ADComputeRSphericalIncrementalStrain< compute_stage >::computeTotalStrainIncrement ( ADRankTwoTensor & total_strain_increment )

overrideprotectedvirtual

Computes the current and old deformation gradients with the assumptions for 1D spherical symmetry geometries: \( \epsilon_{\theta} = \epsilon_{\phi} = \frac{u_r}{r} \).

Reimplemented from ADComputeIncrementalSmallStrain< compute_stage >.

Definition at line 53 of file ADComputeRSphericalIncrementalStrain.C.

 {
   // Deformation gradient calculation in cylindrical coordinates
   ADRankTwoTensor A;  // Deformation gradient
   RankTwoTensor Fbar; // Old Deformation gradient
  
   // Step through calculating the current and old deformation gradients
   // Only diagonal components are nonzero because this is a 1D material
   // Note: x_disp is the radial displacement
   A(0, 0) = (*_grad_disp[0])[_qp](0);
   Fbar(0, 0) = (*_grad_disp_old[0])[_qp](0);
  
   // The polar and azimuthal strains are functions of radial displacement
   if (!MooseUtils::relativeFuzzyEqual(_q_point[_qp](0), 0.0))
   {
     A(1, 1) = (*_disp[0])[_qp] / _q_point[_qp](0);
     Fbar(1, 1) = _disp_old_0[_qp] / _q_point[_qp](0);
   }
  
   // The polar and azimuthal strains are equivalent in this 1D problem
   A(2, 2) = A(1, 1);
   Fbar(2, 2) = Fbar(1, 1);
  
   // very nearly A = gradU - gradUold, adapted to cylindrical coords
   A -= Fbar;
  
   total_strain_increment = 0.5 * (A + A.transpose());
 }

◆ displacementIntegrityCheck()

template<ComputeStage compute_stage>

void ADComputeStrainBase< compute_stage >::displacementIntegrityCheck ( )

protectedvirtualinherited

Reimplemented in ADCompute2DFiniteStrain< compute_stage >, ADCompute2DIncrementalStrain< compute_stage >, and ADCompute2DSmallStrain< compute_stage >.

Definition at line 93 of file ADComputeStrainBase.C.

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

◆ initialSetup()

template<ComputeStage compute_stage>

void ADComputeRSphericalIncrementalStrain< compute_stage >::initialSetup ( )

overridevirtual

Definition at line 41 of file ADComputeRSphericalIncrementalStrain.C.

 {
   ADComputeIncrementalStrainBase<compute_stage>::initialSetup();
  
   const auto & subdomainIDs = _mesh.meshSubdomains();
   for (auto subdomainID : subdomainIDs)
     if (_fe_problem.getCoordSystem(subdomainID) != Moose::COORD_RSPHERICAL)
       mooseError("The coordinate system must be set to RSPHERICAL for 1D R spherical simulations.");
 }

◆ 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 ADComputeRSphericalIncrementalStrain< compute_stage >::validParams ( )

static

Definition at line 23 of file ADComputeRSphericalIncrementalStrain.C.

 {
   InputParameters params = ADComputeIncrementalSmallStrain<compute_stage>::validParams();
   params.addClassDescription(
       "Compute a strain increment for incremental strains in 1D spherical symmetry problems.");
   return params;
 }

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.

◆ _disp_old_0

template<ComputeStage compute_stage>

const VariableValue& ADComputeRSphericalIncrementalStrain< compute_stage >::_disp_old_0

protected

the old value of the first component of the displacements vector

Definition at line 41 of file ADComputeRSphericalIncrementalStrain.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().

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

ADComputeRSphericalIncrementalStrain< compute_stage >::usingComputeIncrementalSmallStrainMembers

protected

Definition at line 43 of file ADComputeRSphericalIncrementalStrain.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/ADComputeRSphericalIncrementalStrain.h
tensor_mechanics/src/materials/ADComputeRSphericalIncrementalStrain.C

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

template<ComputeStage compute_stage> class ADComputeRSphericalIncrementalStrain< compute_stage >

Constructor & Destructor Documentation

◆ ADComputeRSphericalIncrementalStrain()

Member Function Documentation

◆ ADMaterialProperty() [1/4]

◆ ADMaterialProperty() [2/4]

◆ ADMaterialProperty() [3/4]

◆ ADMaterialProperty() [4/4]

◆ computeProperties()

◆ computeTotalStrainIncrement()

◆ displacementIntegrityCheck()

◆ initialSetup()

◆ initQpStatefulProperties()

◆ subtractEigenstrainIncrementFromStrain()

◆ validParams()

Member Data Documentation

◆ _base_name

◆ _current_elem_volume

◆ _disp

◆ _disp_old_0

◆ _eigenstrain_names

◆ _eigenstrains

◆ _eigenstrains_old

◆ _grad_disp

◆ _grad_disp_old

◆ _mechanical_strain_old

◆ _ndisp

◆ _total_strain_old

◆ _volumetric_locking_correction

◆ usingComputeIncrementalSmallStrainMembers

◆ usingComputeIncrementalStrainBaseMembers

◆ usingComputeStrainBaseMembers

◆ usingMaterialMembers

template<ComputeStage compute_stage>
class ADComputeRSphericalIncrementalStrain< compute_stage >