ComputePlaneSmallStrain defines small strains under generalized plane strain and plane stress assumptions, where the out of plane strain can be uniformly or non-uniformly zero or nonzero. More...

#include <ComputePlaneSmallStrain.h>

Inheritance diagram for ComputePlaneSmallStrain:

[legend]

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ComputePlaneSmallStrain (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 Real	computeOutOfPlaneStrain ()

unsigned int	getCurrentSubblockIndex () const
	gets its subblock index for current element

virtual void	computeProperties () override

virtual void	displacementIntegrityCheck () override

virtual void	initQpStatefulProperties () override

Protected Attributes
const SubblockIndexProvider *const	_subblock_id_provider
	A Userobject that carries the subblock ID for all elements.

const unsigned int	_out_of_plane_direction

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
const bool	_scalar_out_of_plane_strain_coupled
	Whether out-of-plane strain scalar variables are coupled.

const bool	_out_of_plane_strain_coupled
	Whether an out-of-plane strain variable is coupled.

const VariableValue &	_out_of_plane_strain
	The out-of-plane strain variable.

std::vector< const VariableValue * >	_scalar_out_of_plane_strain
	The out-of-plane strain scalar variables.

FEProblemBase &	_dmi_fe_problem

Detailed Description

ComputePlaneSmallStrain defines small strains under generalized plane strain and plane stress assumptions, where the out of plane strain can be uniformly or non-uniformly zero or nonzero.

Definition at line 20 of file ComputePlaneSmallStrain.h.

Constructor & Destructor Documentation

◆ ComputePlaneSmallStrain()

ComputePlaneSmallStrain::ComputePlaneSmallStrain ( const InputParameters & parameters )

Definition at line 30 of file ComputePlaneSmallStrain.C.

  : Compute2DSmallStrain(parameters),
    _subblock_id_provider(isParamValid("subblock_index_provider")
                              ? &getUserObject<SubblockIndexProvider>("subblock_index_provider")
                              : nullptr),
    _scalar_out_of_plane_strain_coupled(isCoupledScalar("scalar_out_of_plane_strain")),
    _out_of_plane_strain_coupled(isCoupled("out_of_plane_strain")),
    _out_of_plane_strain(_out_of_plane_strain_coupled ? coupledValue("out_of_plane_strain") : _zero)
{
  if (_out_of_plane_strain_coupled && _scalar_out_of_plane_strain_coupled)
    paramError("_out_of_plane_strain_coupled",
               "Must define only one of out_of_plane_strain or scalar_out_of_plane_strain");
 
  if (_scalar_out_of_plane_strain_coupled)
  {
    const auto nscalar_strains = coupledScalarComponents("scalar_out_of_plane_strain");
    _scalar_out_of_plane_strain.resize(nscalar_strains);
    for (unsigned int i = 0; i < nscalar_strains; ++i)
      _scalar_out_of_plane_strain[i] = &coupledScalarValue("scalar_out_of_plane_strain", i);
  }
}

Member Function Documentation

◆ computeOutOfPlaneStrain()

Real ComputePlaneSmallStrain::computeOutOfPlaneStrain ( )

protectedvirtual

Implements Compute2DSmallStrain.

Definition at line 53 of file ComputePlaneSmallStrain.C.

{
  if (_scalar_out_of_plane_strain_coupled)
    return (*_scalar_out_of_plane_strain[getCurrentSubblockIndex()])[0];
  else
    return _out_of_plane_strain[_qp];
}

◆ computeProperties()

void Compute2DSmallStrain::computeProperties ( )

overrideprotectedvirtualinherited

Reimplemented from ComputeSmallStrain.

Definition at line 52 of file Compute2DSmallStrain.C.

{
  const auto o0 = _out_of_plane_direction;
  const auto o1 = (_out_of_plane_direction + 1) % 3;
  const auto o2 = (_out_of_plane_direction + 2) % 3;
 
  Real volumetric_strain = 0.0;
  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
  {
    _total_strain[_qp](o0, o0) = computeOutOfPlaneStrain();
    _total_strain[_qp](o1, o1) = (*_grad_disp[o1])[_qp](o1);
    _total_strain[_qp](o2, o2) = (*_grad_disp[o2])[_qp](o2);
    _total_strain[_qp](o1, o2) = ((*_grad_disp[o1])[_qp](o2) + (*_grad_disp[o2])[_qp](o1)) / 2.0;
    _total_strain[_qp](o2, o1) = _total_strain[_qp](o1, o2); // force the symmetrical strain tensor
 
    if (_volumetric_locking_correction)
      volumetric_strain += _total_strain[_qp].trace() * _JxW[_qp] * _coord[_qp];
  }
 
  if (_volumetric_locking_correction)
    volumetric_strain /= _current_elem_volume;
 
  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
  {
    if (_volumetric_locking_correction)
    {
      const Real correction = (volumetric_strain - _total_strain[_qp].trace()) / 3.0;
      _total_strain[_qp](0, 0) += correction;
      _total_strain[_qp](1, 1) += correction;
      _total_strain[_qp](2, 2) += correction;
    }
 
    _mechanical_strain[_qp] = _total_strain[_qp];
 
    // Remove the eigenstrains
    for (const auto es : _eigenstrains)
      _mechanical_strain[_qp] -= (*es)[_qp];
  }
}

◆ displacementIntegrityCheck()

void Compute2DSmallStrain::displacementIntegrityCheck ( )

overrideprotectedvirtualinherited

Reimplemented from ComputeStrainBase.

Definition at line 93 of file Compute2DSmallStrain.C.

{
  if (_out_of_plane_direction != 2 && _ndisp != 3)
    mooseError("For 2D simulations where the out-of-plane direction is x or y the number of "
               "supplied displacements must be three.");
  else if (_out_of_plane_direction == 2 && _ndisp != 2)
    mooseError("For 2D simulations where the out-of-plane direction is z the number of supplied "
               "displacements must be two.");
}

Referenced by Compute2DSmallStrain::initialSetup().

◆ getCurrentSubblockIndex()

unsigned int ComputePlaneSmallStrain::getCurrentSubblockIndex ( ) const

inlineprotected

gets its subblock index for current element

Definition at line 31 of file ComputePlaneSmallStrain.h.

  {
    return _subblock_id_provider ? _subblock_id_provider->getSubblockIndex(*_current_elem) : 0;
  };

Referenced by computeOutOfPlaneStrain().

◆ initialSetup()

void Compute2DSmallStrain::initialSetup ( )

overrideinherited

Definition at line 33 of file Compute2DSmallStrain.C.

{
  displacementIntegrityCheck();
  for (unsigned int i = 0; i < 3; ++i)
  {
    if (_out_of_plane_direction == i)
    {
      _disp[i] = &_zero;
      _grad_disp[i] = &_grad_zero;
    }
    else
    {
      _disp[i] = &coupledValue("displacements", i);
      _grad_disp[i] = &coupledGradient("displacements", i);
    }
  }
}

Referenced by ComputeAxisymmetricRZSmallStrain::initialSetup().

◆ initQpStatefulProperties()

void ComputeStrainBase::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Reimplemented in ComputeCosseratIncrementalSmallStrain, and ComputeIncrementalStrainBase.

Definition at line 93 of file ComputeStrainBase.C.

{
  _mechanical_strain[_qp].zero();
  _total_strain[_qp].zero();
}

◆ validParams()

InputParameters ComputePlaneSmallStrain::validParams ( )

static

Definition at line 16 of file ComputePlaneSmallStrain.C.

{
  InputParameters params = Compute2DSmallStrain::validParams();
  params.addClassDescription("Compute a small strain under generalized plane strain assumptions "
                             "where the out of plane strain is generally nonzero.");
  params.addParam<UserObjectName>("subblock_index_provider",
                                  "SubblockIndexProvider user object name");
  params.addCoupledVar("scalar_out_of_plane_strain",
                       "Scalar variable for generalized plane strain");
  params.addCoupledVar("out_of_plane_strain", "Nonlinear variable for plane stress condition");
 
  return params;
}

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

Definition at line 56 of file ComputeStrainBase.h.

Referenced by Compute2DFiniteStrain::computeProperties(), Compute2DSmallStrain::computeProperties(), ComputeFiniteStrain::computeProperties(), ComputeIncrementalStrain::computeProperties(), and ComputeSmallStrain::computeProperties().

◆ _disp

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

protectedinherited

Displacement variables.

Definition at line 38 of file ComputeStrainBase.h.

Referenced by ComputeAxisymmetric1DFiniteStrain::computeGradDispZZ(), ComputeAxisymmetric1DIncrementalStrain::computeGradDispZZ(), ComputeAxisymmetricRZFiniteStrain::computeOutOfPlaneGradDisp(), ComputeAxisymmetricRZIncrementalStrain::computeOutOfPlaneGradDisp(), ComputeAxisymmetricRZSmallStrain::computeOutOfPlaneStrain(), ComputeRSphericalFiniteStrain::computeProperties(), ComputeRSphericalSmallStrain::computeProperties(), ComputeStrainBase::ComputeStrainBase(), ComputeAxisymmetric1DSmallStrain::computeStrainZZ(), ComputeRSphericalIncrementalStrain::computeTotalStrainIncrement(), Compute2DFiniteStrain::initialSetup(), Compute2DIncrementalStrain::initialSetup(), and Compute2DSmallStrain::initialSetup().

◆ _eigenstrain_names

std::vector<MaterialPropertyName> ComputeStrainBase::_eigenstrain_names

protectedinherited

Definition at line 50 of file ComputeStrainBase.h.

Referenced by ComputeIncrementalStrainBase::ComputeIncrementalStrainBase(), and ComputeStrainBase::ComputeStrainBase().

◆ _eigenstrains

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

protectedinherited

Definition at line 51 of file ComputeStrainBase.h.

Referenced by Compute1DSmallStrain::computeProperties(), Compute2DSmallStrain::computeProperties(), ComputeRSphericalSmallStrain::computeProperties(), ComputeSmallStrain::computeProperties(), ComputeCosseratSmallStrain::computeQpProperties(), ComputeCrackTipEnrichmentSmallStrain::computeQpProperties(), ComputeStrainBase::ComputeStrainBase(), and ComputeIncrementalStrainBase::subtractEigenstrainIncrementFromStrain().

◆ _global_strain

const MaterialProperty<RankTwoTensor>* const ComputeStrainBase::_global_strain

protectedinherited

Definition at line 53 of file ComputeStrainBase.h.

Referenced by ComputeSmallStrain::computeProperties(), and ComputeFiniteStrain::computeQpStrain().

◆ _grad_disp

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

protectedinherited

Gradient of displacements.

Definition at line 41 of file ComputeStrainBase.h.

◆ _mechanical_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBase::_mechanical_strain

protectedinherited

Definition at line 46 of file ComputeStrainBase.h.

Referenced by Compute1DSmallStrain::computeProperties(), Compute2DSmallStrain::computeProperties(), ComputeIncrementalStrain::computeProperties(), ComputeRSphericalSmallStrain::computeProperties(), ComputeSmallStrain::computeProperties(), ComputeCosseratIncrementalSmallStrain::computeQpProperties(), ComputeCosseratSmallStrain::computeQpProperties(), ComputeCrackTipEnrichmentSmallStrain::computeQpProperties(), ComputeFiniteStrain::computeQpStrain(), ComputeIncrementalStrainBase::initQpStatefulProperties(), and ComputeStrainBase::initQpStatefulProperties().

◆ _ndisp

unsigned int ComputeStrainBase::_ndisp

protectedinherited

Coupled displacement variables.

Definition at line 35 of file ComputeStrainBase.h.

Referenced by ComputeCrackTipEnrichmentSmallStrain::ComputeCrackTipEnrichmentSmallStrain(), ComputeCrackTipEnrichmentSmallStrain::computeQpProperties(), ComputeStrainBase::ComputeStrainBase(), ComputeStrainBase::displacementIntegrityCheck(), Compute2DFiniteStrain::displacementIntegrityCheck(), Compute2DIncrementalStrain::displacementIntegrityCheck(), Compute2DSmallStrain::displacementIntegrityCheck(), and ComputeIncrementalStrainBase::initialSetup().

◆ _out_of_plane_direction

const unsigned int Compute2DSmallStrain::_out_of_plane_direction

protectedinherited

Definition at line 34 of file Compute2DSmallStrain.h.

Referenced by Compute2DSmallStrain::computeProperties(), Compute2DSmallStrain::displacementIntegrityCheck(), Compute2DSmallStrain::initialSetup(), and ComputeAxisymmetricRZSmallStrain::initialSetup().

◆ _out_of_plane_strain

const VariableValue& ComputePlaneSmallStrain::_out_of_plane_strain

private

The out-of-plane strain variable.

Definition at line 47 of file ComputePlaneSmallStrain.h.

Referenced by computeOutOfPlaneStrain().

◆ _out_of_plane_strain_coupled

const bool ComputePlaneSmallStrain::_out_of_plane_strain_coupled

private

Whether an out-of-plane strain variable is coupled.

Definition at line 44 of file ComputePlaneSmallStrain.h.

Referenced by ComputePlaneSmallStrain().

◆ _scalar_out_of_plane_strain

std::vector<const VariableValue *> ComputePlaneSmallStrain::_scalar_out_of_plane_strain

private

The out-of-plane strain scalar variables.

Definition at line 50 of file ComputePlaneSmallStrain.h.

Referenced by computeOutOfPlaneStrain(), and ComputePlaneSmallStrain().

◆ _scalar_out_of_plane_strain_coupled

const bool ComputePlaneSmallStrain::_scalar_out_of_plane_strain_coupled

private

Whether out-of-plane strain scalar variables are coupled.

Definition at line 41 of file ComputePlaneSmallStrain.h.

Referenced by computeOutOfPlaneStrain(), and ComputePlaneSmallStrain().

◆ _subblock_id_provider

const SubblockIndexProvider* const ComputePlaneSmallStrain::_subblock_id_provider

protected

A Userobject that carries the subblock ID for all elements.

Definition at line 37 of file ComputePlaneSmallStrain.h.

Referenced by getCurrentSubblockIndex().

◆ _total_strain

MaterialProperty<RankTwoTensor>& ComputeStrainBase::_total_strain

protectedinherited

Definition at line 48 of file ComputeStrainBase.h.

Referenced by Compute1DSmallStrain::computeProperties(), Compute2DSmallStrain::computeProperties(), ComputeIncrementalStrain::computeProperties(), ComputeRSphericalSmallStrain::computeProperties(), ComputeSmallStrain::computeProperties(), ComputeCosseratIncrementalSmallStrain::computeQpProperties(), ComputeCosseratSmallStrain::computeQpProperties(), ComputeCrackTipEnrichmentSmallStrain::computeQpProperties(), ComputeFiniteStrain::computeQpStrain(), ComputeIncrementalStrainBase::initQpStatefulProperties(), and ComputeStrainBase::initQpStatefulProperties().

◆ _volumetric_locking_correction

const bool ComputeStrainBase::_volumetric_locking_correction

protectedinherited

Definition at line 55 of file ComputeStrainBase.h.

Referenced by Compute2DFiniteStrain::computeProperties(), Compute2DSmallStrain::computeProperties(), ComputeFiniteStrain::computeProperties(), ComputeIncrementalStrain::computeProperties(), ComputeSmallStrain::computeProperties(), and ComputeStrainBase::ComputeStrainBase().

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

solid_mechanics/include/materials/ComputePlaneSmallStrain.h
solid_mechanics/src/materials/ComputePlaneSmallStrain.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ComputePlaneSmallStrain()

Member Function Documentation

◆ computeOutOfPlaneStrain()

◆ computeProperties()

◆ displacementIntegrityCheck()

◆ getCurrentSubblockIndex()

◆ initialSetup()

◆ initQpStatefulProperties()

◆ validParams()

Member Data Documentation

◆ _base_name

◆ _current_elem_volume

◆ _disp

◆ _eigenstrain_names

◆ _eigenstrains

◆ _global_strain

◆ _grad_disp

◆ _mechanical_strain

◆ _ndisp

◆ _out_of_plane_direction

◆ _out_of_plane_strain

◆ _out_of_plane_strain_coupled

◆ _scalar_out_of_plane_strain

◆ _scalar_out_of_plane_strain_coupled

◆ _subblock_id_provider

◆ _total_strain

◆ _volumetric_locking_correction