https://mooseframework.inl.gov
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | Private Member Functions | Private Attributes | List of all members
ComputeCosseratIncrementalSmallStrain Class Reference

ComputeCosseratIncrementalSmallStrain defines various incremental versions of the Cossserat strain tensor, assuming small strains. More...

#include <ComputeCosseratIncrementalSmallStrain.h>

Inheritance diagram for ComputeCosseratIncrementalSmallStrain:
[legend]

Public Types

typedef DerivativeMaterialPropertyNameInterface::SymbolName SymbolName
 

Public Member Functions

 ComputeCosseratIncrementalSmallStrain (const InputParameters &parameters)
 
void initialSetup () override
 
const GenericMaterialProperty< U, is_ad > & getDefaultMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< U, is_ad > & getDefaultMaterialPropertyByName (const std::string &name)
 
void validateDerivativeMaterialPropertyBase (const std::string &base)
 
const MaterialPropertyName derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
GenericMaterialProperty< U, is_ad > & declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)
 
const GenericMaterialProperty< U, is_ad > & getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")
 
void validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
 
void validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)
 
void validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
void validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")
 
const MaterialPropertyName propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const
 
const MaterialPropertyName propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const
 
const MaterialPropertyName propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const
 
const MaterialPropertyName propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 
const MaterialPropertyName propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const
 

Static Public Member Functions

static InputParameters validParams ()
 

Protected Member Functions

virtual void computeQpProperties ()
 
virtual void initQpStatefulProperties ()
 
void subtractEigenstrainIncrementFromStrain (RankTwoTensor &strain)
 
virtual void displacementIntegrityCheck ()
 

Protected Attributes

MaterialProperty< RankTwoTensor > & _curvature
 the Cosserat curvature strain: curvature_ij = nabla_j CosseratRotation_i
 
const unsigned int _nrots
 the number of Cosserat rotation variables supplied by the user (must be 3 in current implementation)
 
const std::vector< const VariableValue * > _wc
 The Cosserat rotations.
 
const std::vector< const VariableValue * > _wc_old
 The Cosserat rotations.
 
const std::vector< const VariableGradient * > _grad_wc
 Grad(Cosserat rotation)
 
const std::vector< const VariableGradient * > _grad_wc_old
 Grad(Cosserat rotation)
 
const MaterialProperty< RankTwoTensor > & _curvature_old
 the Cosserat curvature strain: curvature_ij = nabla_j CosseratRotation_i
 
MaterialProperty< RankTwoTensor > & _curvature_increment
 _curvature_increment = (curvature - _curvature_old)
 
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.
 
std::vector< const VariableValue * > _disp
 Displacement variables.
 
std::vector< const VariableGradient * > _grad_disp
 Gradient of displacements.
 
const std::string _base_name
 Base name of the material system.
 
MaterialProperty< RankTwoTensor > & _mechanical_strain
 
MaterialProperty< RankTwoTensor > & _total_strain
 
std::vector< MaterialPropertyName > _eigenstrain_names
 
std::vector< const MaterialProperty< RankTwoTensor > * > _eigenstrains
 
const MaterialProperty< RankTwoTensor > *const _global_strain
 
const bool _volumetric_locking_correction
 
const Real & _current_elem_volume
 

Private Member Functions

bool haveMaterialProperty (const std::string &prop_name)
 
std::vector< VariableName > buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)
 
void validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)
 
bool isNotObjectVariable (const VariableName &name)
 

Private Attributes

FEProblemBase_dmi_fe_problem
 

Detailed Description

ComputeCosseratIncrementalSmallStrain defines various incremental versions of the Cossserat strain tensor, assuming small strains.

Definition at line 18 of file ComputeCosseratIncrementalSmallStrain.h.

Constructor & Destructor Documentation

◆ ComputeCosseratIncrementalSmallStrain()

ComputeCosseratIncrementalSmallStrain::ComputeCosseratIncrementalSmallStrain ( const InputParameters parameters)

Definition at line 28 of file ComputeCosseratIncrementalSmallStrain.C.

30 : ComputeIncrementalStrainBase(parameters),
31 _curvature(declareProperty<RankTwoTensor>("curvature")),
32 _nrots(coupledComponents("Cosserat_rotations")),
33 _wc(coupledValues("Cosserat_rotations")),
34 _wc_old(coupledValuesOld("Cosserat_rotations")),
35 _grad_wc(coupledGradients("Cosserat_rotations")),
36 _grad_wc_old(coupledGradientsOld("Cosserat_rotations")),
37 _curvature_old(getMaterialPropertyOld<RankTwoTensor>("curvature")),
38 _curvature_increment(declareProperty<RankTwoTensor>("curvature_increment"))
39{
40 if (_nrots != 3)
41 mooseError("ComputeCosseratSmallStrain: This Material is only defined for 3-dimensional "
42 "simulations so 3 Cosserat rotation variables are needed");
43}
void mooseError(Args &&... args)
MaterialProperty< RankTwoTensor > & _curvature
the Cosserat curvature strain: curvature_ij = nabla_j CosseratRotation_i
const unsigned int _nrots
the number of Cosserat rotation variables supplied by the user (must be 3 in current implementation)
const std::vector< const VariableGradient * > _grad_wc
Grad(Cosserat rotation)
const std::vector< const VariableValue * > _wc
The Cosserat rotations.
MaterialProperty< RankTwoTensor > & _curvature_increment
_curvature_increment = (curvature - _curvature_old)
const MaterialProperty< RankTwoTensor > & _curvature_old
the Cosserat curvature strain: curvature_ij = nabla_j CosseratRotation_i
const std::vector< const VariableGradient * > _grad_wc_old
Grad(Cosserat rotation)
const std::vector< const VariableValue * > _wc_old
The Cosserat rotations.
ComputeIncrementalStrainBase is the base class for strain tensors using incremental formulations.

Member Function Documentation

◆ computeQpProperties()

void ComputeCosseratIncrementalSmallStrain::computeQpProperties ( )
protectedvirtual

Definition at line 54 of file ComputeCosseratIncrementalSmallStrain.C.

55{
57 (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]);
58 auto strain_old = RankTwoTensor::initializeFromRows(
59 (*_grad_disp_old[0])[_qp], (*_grad_disp_old[1])[_qp], (*_grad_disp_old[2])[_qp]);
60 RealVectorValue wc_vector((*_wc[0])[_qp], (*_wc[1])[_qp], (*_wc[2])[_qp]);
61 RealVectorValue wc_vector_old((*_wc_old[0])[_qp], (*_wc_old[1])[_qp], (*_wc_old[2])[_qp]);
62
63 for (unsigned i = 0; i < LIBMESH_DIM; ++i)
64 for (unsigned j = 0; j < LIBMESH_DIM; ++j)
65 for (unsigned k = 0; k < LIBMESH_DIM; ++k)
66 {
67 strain(i, j) += PermutationTensor::eps(i, j, k) * wc_vector(k);
68 strain_old(i, j) += PermutationTensor::eps(i, j, k) * wc_vector_old(k);
69 }
70
71 // Gauss point deformation gradient
72 _deformation_gradient[_qp] = strain;
73 _deformation_gradient[_qp].addIa(1.0);
74
75 const RankTwoTensor total_strain_increment = strain - strain_old;
76
77 _strain_increment[_qp] = total_strain_increment;
78
79 // Remove the eigenstrain increment
81
82 _strain_rate[_qp] = _strain_increment[_qp] / _dt;
83
84 _total_strain[_qp] = _total_strain_old[_qp] + total_strain_increment;
86
88 (*_grad_wc[0])[_qp], (*_grad_wc[1])[_qp], (*_grad_wc[2])[_qp]);
90 (*_grad_wc_old[0])[_qp], (*_grad_wc_old[1])[_qp], (*_grad_wc_old[2])[_qp]);
91 _curvature_increment[_qp] = curv - curv_old;
93
94 // incremental small strain does not include rotation
95 _rotation_increment[_qp].setToIdentity();
96}
MaterialProperty< RankTwoTensor > & _strain_rate
std::vector< const VariableGradient * > _grad_disp_old
const MaterialProperty< RankTwoTensor > & _mechanical_strain_old
void subtractEigenstrainIncrementFromStrain(RankTwoTensor &strain)
MaterialProperty< RankTwoTensor > & _rotation_increment
MaterialProperty< RankTwoTensor > & _strain_increment
MaterialProperty< RankTwoTensor > & _deformation_gradient
const MaterialProperty< RankTwoTensor > & _total_strain_old
std::vector< const VariableGradient * > _grad_disp
Gradient of displacements.
MaterialProperty< RankTwoTensor > & _mechanical_strain
MaterialProperty< RankTwoTensor > & _total_strain
static RankTwoTensorTempl initializeFromRows(const libMesh::TypeVector< T > &row0, const libMesh::TypeVector< T > &row1, const libMesh::TypeVector< T > &row2)
int eps(unsigned int i, unsigned int j)
VectorValue< Real > RealVectorValue

◆ displacementIntegrityCheck()

void ComputeStrainBase::displacementIntegrityCheck ( )
protectedvirtualinherited

Reimplemented in Compute2DFiniteStrain, Compute2DIncrementalStrain, and Compute2DSmallStrain.

Definition at line 83 of file ComputeStrainBase.C.

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

Referenced by ComputeStrainBase::initialSetup().

◆ initialSetup()

void ComputeIncrementalStrainBase::initialSetup ( )
overrideinherited

Definition at line 36 of file ComputeIncrementalStrainBase.C.

37{
39 for (unsigned int i = 0; i < 3; ++i)
40 {
41 if (_fe_problem.isTransient() && i < _ndisp)
42 _grad_disp_old[i] = &coupledGradientOld("displacements", i);
43 else
44 _grad_disp_old[i] = &_grad_zero;
45 }
46}
void initialSetup() override

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

◆ initQpStatefulProperties()

void ComputeCosseratIncrementalSmallStrain::initQpStatefulProperties ( )
protectedvirtual

Reimplemented from ComputeIncrementalStrainBase.

Definition at line 46 of file ComputeCosseratIncrementalSmallStrain.C.

◆ subtractEigenstrainIncrementFromStrain()

void ComputeIncrementalStrainBase::subtractEigenstrainIncrementFromStrain ( RankTwoTensor strain)
protectedinherited

Definition at line 58 of file ComputeIncrementalStrainBase.C.

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

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

◆ validParams()

InputParameters ComputeCosseratIncrementalSmallStrain::validParams ( )
static

Definition at line 20 of file ComputeCosseratIncrementalSmallStrain.C.

21{
23 params.addClassDescription("Compute incremental small Cosserat strains");
24 params.addRequiredCoupledVar("Cosserat_rotations", "The 3 Cosserat rotation variables");
25 return params;
26}
void addRequiredCoupledVar(const std::string &name, const std::string &doc_string)
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _base_name

const std::string ComputeStrainBase::_base_name
protectedinherited

Base name of the material system.

Definition at line 44 of file ComputeStrainBase.h.

Referenced by ComputeStrainBase::ComputeStrainBase().

◆ _current_elem_volume

const Real& ComputeStrainBase::_current_elem_volume
protectedinherited

◆ _curvature

MaterialProperty<RankTwoTensor>& ComputeCosseratIncrementalSmallStrain::_curvature
protected

the Cosserat curvature strain: curvature_ij = nabla_j CosseratRotation_i

Definition at line 31 of file ComputeCosseratIncrementalSmallStrain.h.

Referenced by computeQpProperties(), and initQpStatefulProperties().

◆ _curvature_increment

MaterialProperty<RankTwoTensor>& ComputeCosseratIncrementalSmallStrain::_curvature_increment
protected

_curvature_increment = (curvature - _curvature_old)

Definition at line 52 of file ComputeCosseratIncrementalSmallStrain.h.

Referenced by computeQpProperties().

◆ _curvature_old

const MaterialProperty<RankTwoTensor>& ComputeCosseratIncrementalSmallStrain::_curvature_old
protected

the Cosserat curvature strain: curvature_ij = nabla_j CosseratRotation_i

Definition at line 49 of file ComputeCosseratIncrementalSmallStrain.h.

Referenced by computeQpProperties().

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

◆ _grad_wc

const std::vector<const VariableGradient *> ComputeCosseratIncrementalSmallStrain::_grad_wc
protected

Grad(Cosserat rotation)

Definition at line 43 of file ComputeCosseratIncrementalSmallStrain.h.

Referenced by computeQpProperties().

◆ _grad_wc_old

const std::vector<const VariableGradient *> ComputeCosseratIncrementalSmallStrain::_grad_wc_old
protected

Grad(Cosserat rotation)

Definition at line 46 of file ComputeCosseratIncrementalSmallStrain.h.

Referenced by computeQpProperties().

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

◆ _nrots

const unsigned int ComputeCosseratIncrementalSmallStrain::_nrots
protected

the number of Cosserat rotation variables supplied by the user (must be 3 in current implementation)

Definition at line 34 of file ComputeCosseratIncrementalSmallStrain.h.

Referenced by ComputeCosseratIncrementalSmallStrain().

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

const bool ComputeStrainBase::_volumetric_locking_correction
protectedinherited

◆ _wc

const std::vector<const VariableValue *> ComputeCosseratIncrementalSmallStrain::_wc
protected

The Cosserat rotations.

Definition at line 37 of file ComputeCosseratIncrementalSmallStrain.h.

Referenced by computeQpProperties().

◆ _wc_old

const std::vector<const VariableValue *> ComputeCosseratIncrementalSmallStrain::_wc_old
protected

The Cosserat rotations.

Definition at line 40 of file ComputeCosseratIncrementalSmallStrain.h.

Referenced by computeQpProperties().


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