www.mooseframework.org
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
Compute1DIncrementalStrain Class Referenceabstract

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

#include <Compute1DIncrementalStrain.h>

Inheritance diagram for Compute1DIncrementalStrain:
[legend]

Public Member Functions

 Compute1DIncrementalStrain (const InputParameters &parameters)
 

Protected Member Functions

void computeTotalStrainIncrement (RankTwoTensor &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 Real 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 Real 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 computeProperties () override
 
void initialSetup () override
 
virtual void initQpStatefulProperties () override
 
void subtractEigenstrainIncrementFromStrain (RankTwoTensor &strain)
 
virtual void displacementIntegrityCheck ()
 

Protected Attributes

std::vector< const VariableGradient * > _grad_disp_old
 
MaterialProperty< RankTwoTensor > & _strain_rate
 
MaterialProperty< RankTwoTensor > & _strain_increment
 
MaterialProperty< RankTwoTensor > & _rotation_increment
 
MaterialProperty< RankTwoTensor > & _deformation_gradient
 
const MaterialProperty< RankTwoTensor > & _mechanical_strain_old
 
const MaterialProperty< RankTwoTensor > & _total_strain_old
 
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains_old
 
unsigned int _ndisp
 Coupled displacement variables. More...
 
std::vector< const VariableValue * > _disp
 
std::vector< const VariableGradient * > _grad_disp
 
std::string _base_name
 
MaterialProperty< RankTwoTensor > & _mechanical_strain
 
MaterialProperty< RankTwoTensor > & _total_strain
 
std::vector< MaterialPropertyName > _eigenstrain_names
 
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains
 
const MaterialProperty< RankTwoTensor > * _global_strain
 
bool _volumetric_locking_correction
 
const Real & _current_elem_volume
 

Detailed Description

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

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

Definition at line 26 of file Compute1DIncrementalStrain.h.

Constructor & Destructor Documentation

◆ Compute1DIncrementalStrain()

Compute1DIncrementalStrain::Compute1DIncrementalStrain ( const InputParameters &  parameters)

Definition at line 24 of file Compute1DIncrementalStrain.C.

25  : ComputeIncrementalSmallStrain(parameters)
26 {
27 }
ComputeIncrementalSmallStrain(const InputParameters &parameters)

Member Function Documentation

◆ computeGradDispYY()

virtual Real Compute1DIncrementalStrain::computeGradDispYY ( )
protectedpure virtual

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

Implemented in ComputeAxisymmetric1DIncrementalStrain.

Referenced by computeTotalStrainIncrement().

◆ computeGradDispYYOld()

virtual Real Compute1DIncrementalStrain::computeGradDispYYOld ( )
protectedpure virtual

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

Implemented in ComputeAxisymmetric1DIncrementalStrain.

Referenced by computeTotalStrainIncrement().

◆ computeGradDispZZ()

virtual Real Compute1DIncrementalStrain::computeGradDispZZ ( )
protectedpure virtual

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

Implemented in ComputeAxisymmetric1DIncrementalStrain.

Referenced by computeTotalStrainIncrement().

◆ computeGradDispZZOld()

virtual Real Compute1DIncrementalStrain::computeGradDispZZOld ( )
protectedpure virtual

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

Implemented in ComputeAxisymmetric1DIncrementalStrain.

Referenced by computeTotalStrainIncrement().

◆ computeProperties()

void ComputeIncrementalSmallStrain::computeProperties ( )
overrideprotectedvirtualinherited

Definition at line 32 of file ComputeIncrementalSmallStrain.C.

33 {
34  Real volumetric_strain = 0.0;
35  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
36  {
37  RankTwoTensor total_strain_increment;
38  computeTotalStrainIncrement(total_strain_increment);
39 
40  _strain_increment[_qp] = total_strain_increment;
41 
43  volumetric_strain += total_strain_increment.trace() * _JxW[_qp] * _coord[_qp];
44  }
46  volumetric_strain /= _current_elem_volume;
47 
48  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
49  {
50  Real trace = _strain_increment[_qp].trace();
52  {
53  _strain_increment[_qp](0, 0) += (volumetric_strain - trace) / 3.0;
54  _strain_increment[_qp](1, 1) += (volumetric_strain - trace) / 3.0;
55  _strain_increment[_qp](2, 2) += (volumetric_strain - trace) / 3.0;
56  }
57 
59 
60  // Remove the Eigen strain increment
62 
63  // strain rate
64  if (_dt > 0)
65  _strain_rate[_qp] = _strain_increment[_qp] / _dt;
66  else
67  _strain_rate[_qp].zero();
68 
69  // Update strain in intermediate configuration: rotations are not needed
71  }
72 }
const Real & _current_elem_volume
const MaterialProperty< RankTwoTensor > & _total_strain_old
MaterialProperty< RankTwoTensor > & _strain_increment
MaterialProperty< RankTwoTensor > & _mechanical_strain
virtual void computeTotalStrainIncrement(RankTwoTensor &total_strain_increment)
Computes the current and old deformation gradients and passes back the total strain increment tensor...
MaterialProperty< RankTwoTensor > & _strain_rate
void subtractEigenstrainIncrementFromStrain(RankTwoTensor &strain)
const MaterialProperty< RankTwoTensor > & _mechanical_strain_old
MaterialProperty< RankTwoTensor > & _total_strain

◆ computeTotalStrainIncrement()

void Compute1DIncrementalStrain::computeTotalStrainIncrement ( RankTwoTensor &  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 ComputeIncrementalSmallStrain.

Definition at line 30 of file Compute1DIncrementalStrain.C.

31 {
32  // Deformation gradient calculation for 1D problems
33  RankTwoTensor A(
34  (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]); // Deformation gradient
35  RankTwoTensor Fbar((*_grad_disp_old[0])[_qp],
36  (*_grad_disp_old[1])[_qp],
37  (*_grad_disp_old[2])[_qp]); // Old Deformation gradient
38 
39  // Compute the displacement gradient dUy/dy and dUz/dz value for 1D problems
40  A(1, 1) = computeGradDispYY();
41  A(2, 2) = computeGradDispZZ();
42 
43  Fbar(1, 1) = computeGradDispYYOld();
44  Fbar(2, 2) = computeGradDispZZOld();
45 
46  // Gauss point deformation gradient
47  _deformation_gradient[_qp] = A;
48  _deformation_gradient[_qp].addIa(1.0);
49 
50  A -= Fbar; // very nearly A = gradU - gradUold, adapted to cylindrical coords
51 
52  total_strain_increment = 0.5 * (A + A.transpose());
53 }
MaterialProperty< RankTwoTensor > & _deformation_gradient
virtual Real computeGradDispYYOld()=0
Computes the old dUy/dY; as a virtual function, this function is overwritten for the specific geometr...
virtual Real computeGradDispZZOld()=0
Computes the old dUz/dZ; as a virtual function, this function is overwritten for the specific geometr...
virtual Real computeGradDispZZ()=0
Computes the current dUz/dZ; as a virtual function, this function is overwritten for the specific geo...
std::vector< const VariableGradient * > _grad_disp_old
virtual Real computeGradDispYY()=0
Computes the current dUy/dY; as a virtual function, this function is overwritten for the specific geo...
std::vector< const VariableGradient * > _grad_disp

◆ displacementIntegrityCheck()

void ComputeStrainBase::displacementIntegrityCheck ( )
protectedvirtualinherited

Reimplemented in Compute2DFiniteStrain, Compute2DIncrementalStrain, and Compute2DSmallStrain.

Definition at line 86 of file ComputeStrainBase.C.

Referenced by ComputeStrainBase::initialSetup().

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

◆ initialSetup()

void ComputeIncrementalStrainBase::initialSetup ( )
overrideprotectedinherited

Definition at line 37 of file ComputeIncrementalStrainBase.C.

Referenced by ComputeAxisymmetric1DFiniteStrain::initialSetup(), ComputeAxisymmetricRZIncrementalStrain::initialSetup(), ComputeRSphericalIncrementalStrain::initialSetup(), ComputeRSphericalFiniteStrain::initialSetup(), and ComputeAxisymmetric1DIncrementalStrain::initialSetup().

38 {
40  for (unsigned int i = 0; i < 3; ++i)
41  {
42  if (_fe_problem.isTransient() && i < _ndisp)
43  _grad_disp_old[i] = &coupledGradientOld("displacements", i);
44  else
45  _grad_disp_old[i] = &_grad_zero;
46  }
47 }
unsigned int _ndisp
Coupled displacement variables.
void initialSetup() override
std::vector< const VariableGradient * > _grad_disp_old

◆ initQpStatefulProperties()

void ComputeIncrementalStrainBase::initQpStatefulProperties ( )
overrideprotectedvirtualinherited

Reimplemented from ComputeStrainBase.

Reimplemented in ComputeCosseratIncrementalSmallStrain.

Definition at line 50 of file ComputeIncrementalStrainBase.C.

Referenced by ComputeCosseratIncrementalSmallStrain::initQpStatefulProperties().

51 {
52  _mechanical_strain[_qp].zero();
53  _total_strain[_qp].zero();
54  _deformation_gradient[_qp].zero();
55  _deformation_gradient[_qp].addIa(1.0);
56 
57  // Note that for some models (small strain), the rotation increment is
58  // never updated. Because we always have stateful properties, this method
59  // always gets called, so we can rely on this getting set here without
60  // setting it again when properties get computed.
61  _rotation_increment[_qp].zero();
62  _rotation_increment[_qp].addIa(1.0);
63 }
MaterialProperty< RankTwoTensor > & _deformation_gradient
MaterialProperty< RankTwoTensor > & _mechanical_strain
MaterialProperty< RankTwoTensor > & _rotation_increment
MaterialProperty< RankTwoTensor > & _total_strain

◆ subtractEigenstrainIncrementFromStrain()

void ComputeIncrementalStrainBase::subtractEigenstrainIncrementFromStrain ( RankTwoTensor &  strain)
protectedinherited

Definition at line 66 of file ComputeIncrementalStrainBase.C.

Referenced by ComputeIncrementalSmallStrain::computeProperties(), ComputeCosseratIncrementalSmallStrain::computeQpProperties(), and ComputeFiniteStrain::computeQpStrain().

67 {
68  for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
69  {
70  strain -= (*_eigenstrains[i])[_qp];
71  strain += (*_eigenstrains_old[i])[_qp];
72  }
73 }
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains_old
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains

Member Data Documentation

◆ _base_name

std::string ComputeStrainBase::_base_name
protectedinherited

Definition at line 43 of file ComputeStrainBase.h.

Referenced by ComputeStrainBase::ComputeStrainBase().

◆ _current_elem_volume

const Real& ComputeStrainBase::_current_elem_volume
protectedinherited

◆ _deformation_gradient

MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_deformation_gradient
protectedinherited

◆ _disp

std::vector<const VariableValue *> ComputeStrainBase::_disp
protectedinherited

◆ _eigenstrain_names

std::vector<MaterialPropertyName> ComputeStrainBase::_eigenstrain_names
protectedinherited

◆ _eigenstrains

std::vector<const MaterialProperty<RankTwoTensor> *> ComputeStrainBase::_eigenstrains
protectedinherited

◆ _eigenstrains_old

std::vector<const MaterialProperty<RankTwoTensor> *> ComputeIncrementalStrainBase::_eigenstrains_old
protectedinherited

◆ _global_strain

const MaterialProperty<RankTwoTensor>* ComputeStrainBase::_global_strain
protectedinherited

◆ _grad_disp

std::vector<const VariableGradient *> ComputeStrainBase::_grad_disp
protectedinherited

◆ _grad_disp_old

std::vector<const VariableGradient *> ComputeIncrementalStrainBase::_grad_disp_old
protectedinherited

◆ _mechanical_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBase::_mechanical_strain
protectedinherited

◆ _mechanical_strain_old

const MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_mechanical_strain_old
protectedinherited

◆ _ndisp

unsigned int ComputeStrainBase::_ndisp
protectedinherited

◆ _rotation_increment

MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_rotation_increment
protectedinherited

◆ _strain_increment

MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_strain_increment
protectedinherited

◆ _strain_rate

MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_strain_rate
protectedinherited

◆ _total_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBase::_total_strain
protectedinherited

◆ _total_strain_old

const MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_total_strain_old
protectedinherited

◆ _volumetric_locking_correction

bool ComputeStrainBase::_volumetric_locking_correction
protectedinherited

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