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

ADCompute1DIncrementalStrain defines a strain increment only for incremental small strains in 1D problems, handling strains in other two directions. More...

#include <ADCompute1DIncrementalStrain.h>

Inheritance diagram for ADCompute1DIncrementalStrain< compute_stage >:
[legend]

Public Member Functions

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

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

void computeTotalStrainIncrement (ADRankTwoTensor &total_strain_increment) override
 Computes the current and old deformation gradients with the assumptions for axisymmetric 1D problems, and returns the total strain increment tensor. More...
 
virtual ADReal computeGradDispYY ()=0
 Computes the current dUy/dY; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes. More...
 
virtual Real computeGradDispYYOld ()=0
 Computes the old dUy/dY; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes. More...
 
virtual ADReal computeGradDispZZ ()=0
 Computes the current dUz/dZ; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes. More...
 
virtual Real computeGradDispZZOld ()=0
 Computes the old dUz/dZ; 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
 
virtual void displacementIntegrityCheck ()
 

Protected Attributes

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

ADCompute1DIncrementalStrain defines a strain increment only for incremental small strains in 1D problems, handling strains in other two directions.

ADCompute1DIncrementalStrain contains virtual methods to define the displacement gradients as a general nonzero value.

Definition at line 22 of file ADCompute1DIncrementalStrain.h.

Constructor & Destructor Documentation

◆ ADCompute1DIncrementalStrain()

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

Definition at line 26 of file ADCompute1DIncrementalStrain.C.

29 {
30 }

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

◆ computeGradDispYY()

template<ComputeStage compute_stage>
virtual ADReal ADCompute1DIncrementalStrain< compute_stage >::computeGradDispYY ( )
protectedpure virtual

Computes the current dUy/dY; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes.

◆ computeGradDispYYOld()

template<ComputeStage compute_stage>
virtual Real ADCompute1DIncrementalStrain< compute_stage >::computeGradDispYYOld ( )
protectedpure virtual

Computes the old dUy/dY; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes.

◆ computeGradDispZZ()

template<ComputeStage compute_stage>
virtual ADReal ADCompute1DIncrementalStrain< compute_stage >::computeGradDispZZ ( )
protectedpure virtual

Computes the current dUz/dZ; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes.

◆ computeGradDispZZOld()

template<ComputeStage compute_stage>
virtual Real ADCompute1DIncrementalStrain< compute_stage >::computeGradDispZZOld ( )
protectedpure virtual

Computes the old dUz/dZ; as a virtual function, this function is overwritten for the specific geometries defined by inheriting classes.

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

Computes the current and old deformation gradients with the assumptions for axisymmetric 1D problems, and returns the total strain increment tensor.

Reimplemented from ADComputeIncrementalSmallStrain< compute_stage >.

Definition at line 34 of file ADCompute1DIncrementalStrain.C.

36 {
37  // Deformation gradient calculation for 1D problems
38  ADRankTwoTensor A(
39  (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]); // Deformation gradient
40  RankTwoTensor Fbar((*_grad_disp_old[0])[_qp],
41  (*_grad_disp_old[1])[_qp],
42  (*_grad_disp_old[2])[_qp]); // Old Deformation gradient
43 
44  // Compute the displacement gradient dUy/dy and dUz/dz value for 1D problems
45  A(1, 1) = computeGradDispYY();
46  A(2, 2) = computeGradDispZZ();
47 
48  Fbar(1, 1) = computeGradDispYYOld();
49  Fbar(2, 2) = computeGradDispZZOld();
50 
51  A -= Fbar; // very nearly A = gradU - gradUold, adapted to cylindrical coords
52 
53  total_strain_increment = 0.5 * (A + A.transpose());
54 }

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

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

◆ initialSetup()

template<ComputeStage compute_stage>
void ADComputeIncrementalStrainBase< compute_stage >::initialSetup ( )
overrideinherited

Definition at line 41 of file ADComputeIncrementalStrainBase.C.

42 {
44  for (unsigned int i = 0; i < 3; ++i)
45  {
46  if (_fe_problem.isTransient() && i < _ndisp)
47  _grad_disp_old[i] = &coupledGradientOld("displacements", i);
48  else
49  _grad_disp_old[i] = &_grad_zero;
50  }
51 }

Referenced by ADComputeAxisymmetricRZIncrementalStrain< compute_stage >::initialSetup(), ADComputeRSphericalFiniteStrain< compute_stage >::initialSetup(), and ADComputeRSphericalIncrementalStrain< compute_stage >::initialSetup().

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

Definition at line 18 of file ADCompute1DIncrementalStrain.C.

19 {
21  params.addClassDescription("Compute strain increment for small strains in 1D problems.");
22  return params;
23 }

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.

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

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

◆ usingComputeIncrementalSmallStrainMembers

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

Definition at line 71 of file ADCompute1DIncrementalStrain.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:
ADComputeStrainBase::initialSetup
void initialSetup() override
Definition: ADComputeStrainBase.C:73
ADComputeIncrementalStrainBase::_total_strain_old
const MaterialProperty< RankTwoTensor > & _total_strain_old
Definition: ADComputeIncrementalStrainBase.h:56
ADCompute1DIncrementalStrain::computeGradDispYY
virtual ADReal computeGradDispYY()=0
Computes the current dUy/dY; as a virtual function, this function is overwritten for the specific geo...
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
ADCompute1DIncrementalStrain::computeGradDispZZOld
virtual Real computeGradDispZZOld()=0
Computes the old dUz/dZ; as a virtual function, this function is overwritten for the specific geometr...
ADCompute1DIncrementalStrain::computeGradDispZZ
virtual ADReal computeGradDispZZ()=0
Computes the current dUz/dZ; as a virtual function, this function is overwritten for the specific geo...
ADComputeIncrementalStrainBase::_mechanical_strain_old
const MaterialProperty< RankTwoTensor > & _mechanical_strain_old
Definition: ADComputeIncrementalStrainBase.h:55
ADComputeIncrementalSmallStrain::validParams
static InputParameters validParams()
Definition: ADComputeIncrementalSmallStrain.C:19
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
ADComputeIncrementalStrainBase::_eigenstrains_old
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains_old
Definition: ADComputeIncrementalStrainBase.h:58
ADCompute1DIncrementalStrain::computeGradDispYYOld
virtual Real computeGradDispYYOld()=0
Computes the old dUy/dY; as a virtual function, this function is overwritten for the specific geometr...
ADComputeIncrementalSmallStrain
ADComputeIncrementalSmallStrain defines a strain increment and rotation increment (=1),...
Definition: ADComputeIncrementalSmallStrain.h:19
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
ADComputeStrainBase::_volumetric_locking_correction
const bool _volumetric_locking_correction
Definition: ADComputeStrainBase.h:70