www.mooseframework.org
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
ComputeIncrementalSmallStrain Class Reference

ComputeIncrementalSmallStrain defines a strain increment and rotation increment (=1), for small strains. More...

#include <ComputeIncrementalSmallStrain.h>

Inheritance diagram for ComputeIncrementalSmallStrain:
[legend]

Public Member Functions

 ComputeIncrementalSmallStrain (const InputParameters &parameters)
 

Protected Member Functions

virtual void computeProperties () override
 
virtual void computeTotalStrainIncrement (RankTwoTensor &total_strain_increment)
 Computes the current and old deformation gradients and passes back the total strain increment tensor. More...
 
virtual void initQpStatefulProperties () override
 
void subtractEigenstrainIncrementFromStrain (RankTwoTensor &strain)
 

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
 
bool _volumetric_locking_correction
 
const Real & _current_elem_volume
 

Detailed Description

ComputeIncrementalSmallStrain defines a strain increment and rotation increment (=1), for small strains.

Definition at line 16 of file ComputeIncrementalSmallStrain.h.

Constructor & Destructor Documentation

ComputeIncrementalSmallStrain::ComputeIncrementalSmallStrain ( const InputParameters &  parameters)

Definition at line 21 of file ComputeIncrementalSmallStrain.C.

22  : ComputeIncrementalStrainBase(parameters)
23 {
24 }
ComputeIncrementalStrainBase(const InputParameters &parameters)

Member Function Documentation

void ComputeIncrementalSmallStrain::computeProperties ( )
overrideprotectedvirtual

Definition at line 27 of file ComputeIncrementalSmallStrain.C.

28 {
29  Real volumetric_strain = 0.0;
30  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
31  {
32  RankTwoTensor total_strain_increment;
33  computeTotalStrainIncrement(total_strain_increment);
34 
35  _strain_increment[_qp] = total_strain_increment;
36 
38  volumetric_strain += total_strain_increment.trace() * _JxW[_qp] * _coord[_qp];
39  }
41  volumetric_strain /= _current_elem_volume;
42 
43  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
44  {
45  Real trace = _strain_increment[_qp].trace();
47  {
48  _strain_increment[_qp](0, 0) += (volumetric_strain - trace) / 3.0;
49  _strain_increment[_qp](1, 1) += (volumetric_strain - trace) / 3.0;
50  _strain_increment[_qp](2, 2) += (volumetric_strain - trace) / 3.0;
51  }
52 
54 
55  // Remove the Eigen strain increment
56  subtractEigenstrainIncrementFromStrain(_strain_increment[_qp]);
57 
58  // strain rate
59  if (_dt > 0)
60  _strain_rate[_qp] = _strain_increment[_qp] / _dt;
61  else
62  _strain_rate[_qp].zero();
63 
64  // Update strain in intermediate configuration: rotations are not needed
65  _mechanical_strain[_qp] = _mechanical_strain_old[_qp] + _strain_increment[_qp];
66  }
67 }
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
void ComputeIncrementalSmallStrain::computeTotalStrainIncrement ( RankTwoTensor &  total_strain_increment)
protectedvirtual

Computes the current and old deformation gradients and passes back the total strain increment tensor.

Reimplemented in Compute2DIncrementalStrain, ComputeRSphericalIncrementalStrain, and Compute1DIncrementalStrain.

Definition at line 70 of file ComputeIncrementalSmallStrain.C.

Referenced by computeProperties().

71 {
72  // Deformation gradient
73  RankTwoTensor A(
74  (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]); // Deformation gradient
75  RankTwoTensor Fbar((*_grad_disp_old[0])[_qp],
76  (*_grad_disp_old[1])[_qp],
77  (*_grad_disp_old[2])[_qp]); // Old Deformation gradient
78 
79  _deformation_gradient[_qp] = A;
80  _deformation_gradient[_qp].addIa(1.0);
81 
82  A -= Fbar; // A = grad_disp - grad_disp_old
83 
84  total_strain_increment = 0.5 * (A + A.transpose());
85 }
MaterialProperty< RankTwoTensor > & _deformation_gradient
std::vector< const VariableGradient * > _grad_disp_old
std::vector< const VariableGradient * > _grad_disp
void ComputeIncrementalStrainBase::initQpStatefulProperties ( )
overrideprotectedvirtualinherited

Reimplemented from ComputeStrainBase.

Reimplemented in ComputeCosseratIncrementalSmallStrain.

Definition at line 44 of file ComputeIncrementalStrainBase.C.

Referenced by ComputeCosseratIncrementalSmallStrain::initQpStatefulProperties(), and ComputeIncrementalStrainBase::~ComputeIncrementalStrainBase().

45 {
46  _mechanical_strain[_qp].zero();
47  _total_strain[_qp].zero();
48  _deformation_gradient[_qp].zero();
49  _deformation_gradient[_qp].addIa(1.0);
50 
51  // Note that for some models (small strain), the rotation increment is
52  // never updated. Because we always have stateful properties, this method
53  // always gets called, so we can rely on this getting set here without
54  // setting it again when properties get computed.
55  _rotation_increment[_qp].zero();
56  _rotation_increment[_qp].addIa(1.0);
57 }
MaterialProperty< RankTwoTensor > & _deformation_gradient
MaterialProperty< RankTwoTensor > & _mechanical_strain
MaterialProperty< RankTwoTensor > & _rotation_increment
MaterialProperty< RankTwoTensor > & _total_strain
void ComputeIncrementalStrainBase::subtractEigenstrainIncrementFromStrain ( RankTwoTensor &  strain)
protectedinherited

Definition at line 60 of file ComputeIncrementalStrainBase.C.

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

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

Member Data Documentation

std::string ComputeStrainBase::_base_name
protectedinherited

Definition at line 33 of file ComputeStrainBase.h.

Referenced by ComputeStrainBase::ComputeStrainBase().

const Real& ComputeStrainBase::_current_elem_volume
protectedinherited
MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_deformation_gradient
protectedinherited
std::vector<const VariableValue *> ComputeStrainBase::_disp
protectedinherited
std::vector<MaterialPropertyName> ComputeStrainBase::_eigenstrain_names
protectedinherited
std::vector<const MaterialProperty<RankTwoTensor> *> ComputeStrainBase::_eigenstrains
protectedinherited
std::vector<const MaterialProperty<RankTwoTensor> *> ComputeIncrementalStrainBase::_eigenstrains_old
protectedinherited
std::vector<const VariableGradient *> ComputeStrainBase::_grad_disp
protectedinherited
std::vector<const VariableGradient *> ComputeIncrementalStrainBase::_grad_disp_old
protectedinherited
MaterialProperty<RankTwoTensor>& ComputeStrainBase::_mechanical_strain
protectedinherited
const MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_mechanical_strain_old
protectedinherited
unsigned int ComputeStrainBase::_ndisp
protectedinherited
MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_rotation_increment
protectedinherited
MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_strain_increment
protectedinherited
MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_strain_rate
protectedinherited
MaterialProperty<RankTwoTensor>& ComputeStrainBase::_total_strain
protectedinherited
const MaterialProperty<RankTwoTensor>& ComputeIncrementalStrainBase::_total_strain_old
protectedinherited
bool ComputeStrainBase::_volumetric_locking_correction
protectedinherited

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