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

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
 

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
 
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 16 of file ADComputeRSphericalIncrementalStrain.h.

Constructor & Destructor Documentation

◆ ADComputeRSphericalIncrementalStrain()

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

Definition at line 26 of file ADComputeRSphericalIncrementalStrain.C.

29  _disp_old_0(coupledValueOld("displacements", 0))
30 {
31 }
ADComputeIncrementalSmallStrain defines a strain increment and rotation increment (=1)...
const VariableValue & _disp_old_0
the old value of the first component of the displacements vector

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 29 of file ADComputeIncrementalSmallStrain.C.

30 {
31  ADReal volumetric_strain = 0.0;
32  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
33  {
34  ADRankTwoTensor total_strain_increment;
35  computeTotalStrainIncrement(total_strain_increment);
36 
37  _strain_increment[_qp] = total_strain_increment;
38 
40  volumetric_strain += total_strain_increment.trace() * _JxW[_qp] * _coord[_qp];
41  }
42 
44  volumetric_strain /= _current_elem_volume;
45 
46  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
47  {
49  {
50  const auto correction = (volumetric_strain - _strain_increment[_qp].trace()) / 3.0;
51  _strain_increment[_qp](0, 0) += correction;
52  _strain_increment[_qp](1, 1) += correction;
53  _strain_increment[_qp](2, 2) += correction;
54  }
55 
56  _total_strain[_qp] = _strain_increment[_qp] + _total_strain_old[_qp];
57 
58  // Remove the Eigen strain increment
59  subtractEigenstrainIncrementFromStrain(_strain_increment[_qp]);
60 
61  // strain rate
62  if (_dt > 0)
63  _strain_rate[_qp] = _strain_increment[_qp] / _dt;
64  else
65  _strain_rate[_qp].zero();
66 
67  // Update strain in intermediate configuration: rotations are not needed
68  _mechanical_strain[_qp] = _strain_increment[_qp] + _mechanical_strain_old[_qp];
69 
70  // incremental small strain does not include rotation
71  _rotation_increment[_qp].setToIdentity();
72  }
73 
74  copyDualNumbersToValues();
75 }
const bool _volumetric_locking_correction
virtual void computeTotalStrainIncrement(ADRankTwoTensor &total_strain_increment)
Computes the current and old deformation gradients and passes back the total strain increment tensor...
void subtractEigenstrainIncrementFromStrain(ADRankTwoTensor &strain)
const Real & _current_elem_volume
const MaterialProperty< RankTwoTensor > & _mechanical_strain_old
const MaterialProperty< RankTwoTensor > & _total_strain_old

◆ 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 47 of file ADComputeRSphericalIncrementalStrain.C.

49 {
50  // Deformation gradient calculation in cylindrical coordinates
51  ADRankTwoTensor A; // Deformation gradient
52  RankTwoTensor Fbar; // Old Deformation gradient
53 
54  // Step through calculating the current and old deformation gradients
55  // Only diagonal components are nonzero because this is a 1D material
56  // Note: x_disp is the radial displacement
57  A(0, 0) = (*_grad_disp[0])[_qp](0);
58  Fbar(0, 0) = (*_grad_disp_old[0])[_qp](0);
59 
60  // The polar and azimuthal strains are functions of radial displacement
61  if (!MooseUtils::relativeFuzzyEqual(_q_point[_qp](0), 0.0))
62  {
63  A(1, 1) = (*_disp[0])[_qp] / _q_point[_qp](0);
64  Fbar(1, 1) = _disp_old_0[_qp] / _q_point[_qp](0);
65  }
66 
67  // The polar and azimuthal strains are equivalent in this 1D problem
68  A(2, 2) = A(1, 1);
69  Fbar(2, 2) = Fbar(1, 1);
70 
71  // very nearly A = gradU - gradUold, adapted to cylindrical coords
72  A -= Fbar;
73 
74  total_strain_increment = 0.5 * (A + A.transpose());
75 }
const VariableValue & _disp_old_0
the old value of the first component of the displacements vector
std::vector< const ADVariableValue * > _disp
std::vector< const ADVariableGradient * > _grad_disp
std::vector< const VariableGradient * > _grad_disp_old

◆ 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 87 of file ADComputeStrainBase.C.

88 {
89  // Checking for consistency between mesh size and length of the provided displacements vector
90  if (_ndisp != _mesh.dimension())
91  paramError(
92  "displacements",
93  "The number of variables supplied in 'displacements' must match the mesh dimension.");
94 }
const unsigned int _ndisp
Coupled displacement variables.

◆ initialSetup()

template<ComputeStage compute_stage>
void ADComputeRSphericalIncrementalStrain< compute_stage >::initialSetup ( )
overridevirtual

Definition at line 35 of file ADComputeRSphericalIncrementalStrain.C.

36 {
38 
39  const auto & subdomainIDs = _mesh.meshSubdomains();
40  for (auto subdomainID : subdomainIDs)
41  if (_fe_problem.getCoordSystem(subdomainID) != Moose::COORD_RSPHERICAL)
42  mooseError("The coordinate system must be set to RSPHERICAL for 1D R spherical simulations.");
43 }

◆ initQpStatefulProperties()

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

Reimplemented from ADComputeStrainBase< compute_stage >.

Definition at line 48 of file ADComputeIncrementalStrainBase.C.

49 {
50  _mechanical_strain[_qp].zero();
51  _total_strain[_qp].zero();
52 }

◆ subtractEigenstrainIncrementFromStrain()

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

Definition at line 56 of file ADComputeIncrementalStrainBase.C.

58 {
59  for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
60  {
61  strain -= (*_eigenstrains[i])[_qp];
62  strain += (*_eigenstrains_old[i])[_qp];
63  }
64 }
std::vector< const ADMaterialProperty(RankTwoTensor) * > _eigenstrains
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains_old

Member Data Documentation

◆ _base_name

template<ComputeStage compute_stage>
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 70 of file ADComputeStrainBase.h.

◆ _disp

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

Definition at line 56 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 40 of file ADComputeRSphericalIncrementalStrain.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 57 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 48 of file ADComputeIncrementalStrainBase.h.

◆ _mechanical_strain_old

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

Definition at line 54 of file ADComputeIncrementalStrainBase.h.

◆ _ndisp

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

Coupled displacement variables.

Definition at line 55 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 55 of file ADComputeIncrementalStrainBase.h.

◆ _volumetric_locking_correction

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

◆ usingComputeIncrementalSmallStrainMembers

template<ComputeStage compute_stage>
ADComputeRSphericalIncrementalStrain< compute_stage >::usingComputeIncrementalSmallStrainMembers
protected

Definition at line 42 of file ADComputeRSphericalIncrementalStrain.h.

◆ usingComputeIncrementalStrainBaseMembers

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

Definition at line 43 of file ADComputeIncrementalSmallStrain.h.

◆ usingComputeStrainBaseMembers

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

Definition at line 59 of file ADComputeIncrementalStrainBase.h.

◆ usingMaterialMembers

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

Definition at line 72 of file ADComputeStrainBase.h.


The documentation for this class was generated from the following files: