ComputeAxisymmetric1DIncrementalStrain defines a strain increment only for incremental small strains in an Axisymmetric 1D problem. More...

#include <ComputeAxisymmetric1DIncrementalStrain.h>

Inheritance diagram for ComputeAxisymmetric1DIncrementalStrain:

[legend]

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ComputeAxisymmetric1DIncrementalStrain (const InputParameters &parameters)

void	initialSetup () override

virtual void	computeProperties () 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
Real	computeGradDispYY () override
	Computes the current dUy/dy for axisymmetric problems.

Real	computeGradDispYYOld () override
	Computes the old dUy/dy for axisymmetric problems.

Real	computeGradDispZZ () override
	Computes the current dUz/dz for axisymmetric problems, where \( \epsilon_{\theta} = \frac{u_r}{r} \).

Real	computeGradDispZZOld () override
	Computes the old dUz/dz for axisymmetric problems, where \( \epsilon_{\theta-old} = \frac{u_{r-old}}{r_{old}} \).

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

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.

virtual void	initQpStatefulProperties () override

void	subtractEigenstrainIncrementFromStrain (RankTwoTensor &strain)

virtual void	displacementIntegrityCheck ()

Protected Attributes
const VariableValue &	_disp_old_0
	the old value of the first component of the displacements vector

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

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

bool	_has_scalar_out_of_plane_strain
	Whether an out-of-plane strain scalar variable is coupled.

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


const VariableValue &	_out_of_plane_strain
	Current and old values of the out-of-plane strain variable.

const VariableValue &	_out_of_plane_strain_old


std::vector< const VariableValue * >	_scalar_out_of_plane_strain
	Current and old values of the out-of-plane strain scalar variable.

std::vector< const VariableValue * >	_scalar_out_of_plane_strain_old

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

ComputeAxisymmetric1DIncrementalStrain defines a strain increment only for incremental small strains in an Axisymmetric 1D problem.

The COORD_TYPE in the Problem block must be set to RZ.

Definition at line 20 of file ComputeAxisymmetric1DIncrementalStrain.h.

Constructor & Destructor Documentation

◆ ComputeAxisymmetric1DIncrementalStrain()

ComputeAxisymmetric1DIncrementalStrain::ComputeAxisymmetric1DIncrementalStrain ( const InputParameters & parameters )

Definition at line 29 of file ComputeAxisymmetric1DIncrementalStrain.C.

  : Compute1DIncrementalStrain(parameters),
    _disp_old_0(coupledValueOld("displacements", 0)),
    _subblock_id_provider(isParamValid("subblock_index_provider")
                              ? &getUserObject<SubblockIndexProvider>("subblock_index_provider")
                              : nullptr),
    _has_out_of_plane_strain(isCoupled("out_of_plane_strain")),
    _out_of_plane_strain(_has_out_of_plane_strain ? coupledValue("out_of_plane_strain") : _zero),
    _out_of_plane_strain_old(_has_out_of_plane_strain ? coupledValueOld("out_of_plane_strain")
                                                      : _zero),
    _has_scalar_out_of_plane_strain(isCoupledScalar("scalar_out_of_plane_strain"))
{
  if (_has_out_of_plane_strain && _has_scalar_out_of_plane_strain)
    mooseError("Must define only one of out_of_plane_strain or scalar_out_of_plane_strain");
 
  if (_has_scalar_out_of_plane_strain)
  {
    const auto nscalar_strains = coupledScalarComponents("scalar_out_of_plane_strain");
    _scalar_out_of_plane_strain.resize(nscalar_strains);
    _scalar_out_of_plane_strain_old.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);
      _scalar_out_of_plane_strain_old[i] = &coupledScalarValueOld("scalar_out_of_plane_strain", i);
    }
  }
}

Member Function Documentation

◆ computeGradDispYY()

Real ComputeAxisymmetric1DIncrementalStrain::computeGradDispYY ( )

overrideprotectedvirtual

Computes the current dUy/dy for axisymmetric problems.

Implements Compute1DIncrementalStrain.

Definition at line 68 of file ComputeAxisymmetric1DIncrementalStrain.C.

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

◆ computeGradDispYYOld()

Real ComputeAxisymmetric1DIncrementalStrain::computeGradDispYYOld ( )

overrideprotectedvirtual

Computes the old dUy/dy for axisymmetric problems.

Implements Compute1DIncrementalStrain.

Definition at line 77 of file ComputeAxisymmetric1DIncrementalStrain.C.

{
  if (_has_scalar_out_of_plane_strain)
    return (*_scalar_out_of_plane_strain_old[getCurrentSubblockIndex()])[0];
  else
    return _out_of_plane_strain_old[_qp];
}

◆ computeGradDispZZ()

Real ComputeAxisymmetric1DIncrementalStrain::computeGradDispZZ ( )

overrideprotectedvirtual

Computes the current dUz/dz for axisymmetric problems, where \( \epsilon_{\theta} = \frac{u_r}{r} \).

Implements Compute1DIncrementalStrain.

Definition at line 86 of file ComputeAxisymmetric1DIncrementalStrain.C.

{
  if (!MooseUtils::absoluteFuzzyEqual(_q_point[_qp](0), 0.0))
    return (*_disp[0])[_qp] / _q_point[_qp](0);
  else
    return 0.0;
}

◆ computeGradDispZZOld()

Real ComputeAxisymmetric1DIncrementalStrain::computeGradDispZZOld ( )

overrideprotectedvirtual

Computes the old dUz/dz for axisymmetric problems, where \( \epsilon_{\theta-old} = \frac{u_{r-old}}{r_{old}} \).

Implements Compute1DIncrementalStrain.

Definition at line 95 of file ComputeAxisymmetric1DIncrementalStrain.C.

{
  if (!MooseUtils::absoluteFuzzyEqual(_q_point[_qp](0), 0.0))
    return _disp_old_0[_qp] / _q_point[_qp](0);
  else
    return 0.0;
}

◆ computeProperties()

void ComputeIncrementalStrain::computeProperties ( )

overridevirtualinherited

Definition at line 35 of file ComputeIncrementalStrain.C.

{
  Real volumetric_strain = 0.0;
  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
  {
    RankTwoTensor total_strain_increment;
    computeTotalStrainIncrement(total_strain_increment);
 
    _strain_increment[_qp] = total_strain_increment;
 
    if (_volumetric_locking_correction)
      volumetric_strain += total_strain_increment.trace() * _JxW[_qp] * _coord[_qp];
  }
  if (_volumetric_locking_correction)
    volumetric_strain /= _current_elem_volume;
 
  for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
  {
    Real trace = _strain_increment[_qp].trace();
    if (_volumetric_locking_correction)
    {
      _strain_increment[_qp](0, 0) += (volumetric_strain - trace) / 3.0;
      _strain_increment[_qp](1, 1) += (volumetric_strain - trace) / 3.0;
      _strain_increment[_qp](2, 2) += (volumetric_strain - trace) / 3.0;
    }
 
    _total_strain[_qp] = _total_strain_old[_qp] + _strain_increment[_qp];
 
    // Remove the Eigen strain increment
    subtractEigenstrainIncrementFromStrain(_strain_increment[_qp]);
 
    // strain rate
    if (_dt > 0)
      _strain_rate[_qp] = _strain_increment[_qp] / _dt;
    else
      _strain_rate[_qp].zero();
 
    // Update strain in intermediate configuration: rotations are not needed
    _mechanical_strain[_qp] = _mechanical_strain_old[_qp] + _strain_increment[_qp];
 
    // incremental small strain does not include rotation
    _rotation_increment[_qp].setToIdentity();
  }
}

◆ computeTotalStrainIncrement()

void Compute1DIncrementalStrain::computeTotalStrainIncrement ( RankTwoTensor & total_strain_increment )

overrideprotectedvirtualinherited

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

Reimplemented from ComputeIncrementalStrain.

Definition at line 29 of file Compute1DIncrementalStrain.C.

{
  // Deformation gradient calculation for 1D problems
  auto A = RankTwoTensor::initializeFromRows(
      (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]);
 
  // Old Deformation gradient
  auto Fbar = RankTwoTensor::initializeFromRows(
      (*_grad_disp_old[0])[_qp], (*_grad_disp_old[1])[_qp], (*_grad_disp_old[2])[_qp]);
 
  // Compute the displacement gradient dUy/dy and dUz/dz value for 1D problems
  A(1, 1) = computeGradDispYY();
  A(2, 2) = computeGradDispZZ();
 
  Fbar(1, 1) = computeGradDispYYOld();
  Fbar(2, 2) = computeGradDispZZOld();
 
  // Gauss point deformation gradient
  _deformation_gradient[_qp] = A;
  _deformation_gradient[_qp].addIa(1.0);
 
  A -= Fbar; // very nearly A = gradU - gradUold, adapted to cylindrical coords
 
  total_strain_increment = 0.5 * (A + A.transpose());
}

◆ displacementIntegrityCheck()

void ComputeStrainBase::displacementIntegrityCheck ( )

protectedvirtualinherited

Reimplemented in Compute2DFiniteStrain, Compute2DIncrementalStrain, and Compute2DSmallStrain.

Definition at line 83 of file ComputeStrainBase.C.

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

Referenced by ComputeStrainBase::initialSetup().

◆ getCurrentSubblockIndex()

unsigned int ComputeAxisymmetric1DIncrementalStrain::getCurrentSubblockIndex ( ) const

inlineprotected

gets its subblock index for current element

Definition at line 45 of file ComputeAxisymmetric1DIncrementalStrain.h.

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

Referenced by computeGradDispYY(), and computeGradDispYYOld().

◆ initialSetup()

void ComputeAxisymmetric1DIncrementalStrain::initialSetup ( )

override

Definition at line 59 of file ComputeAxisymmetric1DIncrementalStrain.C.

{
  ComputeIncrementalStrainBase::initialSetup();
 
  if (getBlockCoordSystem() != Moose::COORD_RZ)
    mooseError("The coordinate system must be set to RZ for Axisymmetric 1D simulations");
}

◆ initQpStatefulProperties()

void ComputeIncrementalStrainBase::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Reimplemented from ComputeStrainBase.

Reimplemented in ComputeCosseratIncrementalSmallStrain.

Definition at line 49 of file ComputeIncrementalStrainBase.C.

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

Referenced by ComputeCosseratIncrementalSmallStrain::initQpStatefulProperties().

◆ subtractEigenstrainIncrementFromStrain()

void ComputeIncrementalStrainBase::subtractEigenstrainIncrementFromStrain ( RankTwoTensor & strain )

protectedinherited

Definition at line 58 of file ComputeIncrementalStrainBase.C.

{
  for (unsigned int i = 0; i < _eigenstrains.size(); ++i)
  {
    strain -= (*_eigenstrains[i])[_qp];
    strain += (*_eigenstrains_old[i])[_qp];
  }
}

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

◆ validParams()

InputParameters ComputeAxisymmetric1DIncrementalStrain::validParams ( )

static

Definition at line 16 of file ComputeAxisymmetric1DIncrementalStrain.C.

{
  InputParameters params = Compute1DIncrementalStrain::validParams();
  params.addClassDescription(
      "Compute strain increment for small strains in an axisymmetric 1D problem");
  params.addParam<UserObjectName>("subblock_index_provider",
                                  "SubblockIndexProvider user object name");
  params.addCoupledVar("scalar_out_of_plane_strain", "Scalar variable for axisymmetric 1D problem");
  params.addCoupledVar("out_of_plane_strain", "Nonlinear variable for axisymmetric 1D problem");
 
  return params;
}