ComputeAxisymmetric1DSmallStrain defines small strains in an Axisymmetric 1D problem. More...

#include <ComputeAxisymmetric1DSmallStrain.h>

Inheritance diagram for ComputeAxisymmetric1DSmallStrain:

Public Member Functions
	ComputeAxisymmetric1DSmallStrain (const InputParameters &parameters)

void	initialSetup () override

void	computeProperties () override

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
Real	computeStrainYY () override
	Computes the strain_yy for axisymmetric problems. More...

Real	computeStrainZZ () override
	Computes the strain_zz for axisymmetric problems, where \( \epsilon_{\theta} = \frac{u_r}{r} \). More...

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

virtual void	initQpStatefulProperties () override

virtual void	displacementIntegrityCheck ()

Protected Attributes
const SubblockIndexProvider *	_subblock_id_provider

const bool	_has_out_of_plane_strain

const VariableValue &	_out_of_plane_strain

const bool	_has_scalar_out_of_plane_strain

unsigned int	_nscalar_strains

std::vector< const VariableValue * >	_scalar_out_of_plane_strain

unsigned int	_ndisp
	Coupled displacement variables. More...

std::vector< const VariableValue * >	_disp

std::vector< const VariableGradient * >	_grad_disp

const 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

const bool	_volumetric_locking_correction

const Real &	_current_elem_volume

Detailed Description

ComputeAxisymmetric1DSmallStrain defines small strains in an Axisymmetric 1D problem.

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

Definition at line 24 of file ComputeAxisymmetric1DSmallStrain.h.

Constructor & Destructor Documentation

◆ ComputeAxisymmetric1DSmallStrain()

ComputeAxisymmetric1DSmallStrain::ComputeAxisymmetric1DSmallStrain ( const InputParameters & parameters )

Definition at line 29 of file ComputeAxisymmetric1DSmallStrain.C.

   : Compute1DSmallStrain(parameters),
     _subblock_id_provider(isParamValid("subblock_index_provider")
                               ? &getUserObject<SubblockIndexProvider>("subblock_index_provider")
                               : nullptr),
     _has_out_of_plane_strain(isParamValid("out_of_plane_strain")),
     _out_of_plane_strain(_has_out_of_plane_strain ? coupledValue("out_of_plane_strain") : _zero),
     _has_scalar_out_of_plane_strain(isParamValid("scalar_out_of_plane_strain")),
     _nscalar_strains(coupledScalarComponents("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)
   {
     _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

◆ computeProperties()

void Compute1DSmallStrain::computeProperties ( )

overridevirtualinherited

Reimplemented from ComputeSmallStrain.

Definition at line 30 of file Compute1DSmallStrain.C.

 {
   for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
   {
     _total_strain[_qp](0, 0) = (*_grad_disp[0])[_qp](0);
     _total_strain[_qp](1, 1) = computeStrainYY();
     _total_strain[_qp](2, 2) = computeStrainZZ();
  
     _mechanical_strain[_qp] = _total_strain[_qp];
  
     // Remove the eigenstrain
     for (auto es : _eigenstrains)
       _mechanical_strain[_qp] -= (*es)[_qp];
   }
 }

◆ computeStrainYY()

Real ComputeAxisymmetric1DSmallStrain::computeStrainYY ( )

overrideprotectedvirtual

Computes the strain_yy for axisymmetric problems.

Implements Compute1DSmallStrain.

Definition at line 61 of file ComputeAxisymmetric1DSmallStrain.C.

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

◆ computeStrainZZ()

Real ComputeAxisymmetric1DSmallStrain::computeStrainZZ ( )

overrideprotectedvirtual

Computes the strain_zz for axisymmetric problems, where \( \epsilon_{\theta} = \frac{u_r}{r} \).

Implements Compute1DSmallStrain.

Definition at line 70 of file ComputeAxisymmetric1DSmallStrain.C.

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

◆ displacementIntegrityCheck()

void ComputeStrainBase::displacementIntegrityCheck ( )

protectedvirtualinherited

Reimplemented in Compute2DFiniteStrain, Compute2DIncrementalStrain, and Compute2DSmallStrain.

Definition at line 94 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 ComputeAxisymmetric1DSmallStrain::getCurrentSubblockIndex ( ) const

inlineprotected

gets its subblock index for current element

Definition at line 42 of file ComputeAxisymmetric1DSmallStrain.h.

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

Referenced by computeStrainYY().

◆ initialSetup()

void ComputeAxisymmetric1DSmallStrain::initialSetup ( )

override

Definition at line 52 of file ComputeAxisymmetric1DSmallStrain.C.

 {
   ComputeStrainBase::initialSetup();
  
   if (getBlockCoordSystem() != Moose::COORD_RZ)
     mooseError("The coordinate system must be set to RZ for Axisymmetric geometries.");
 }

◆ initQpStatefulProperties()

void ComputeStrainBase::initQpStatefulProperties ( )

overrideprotectedvirtualinherited

Reimplemented in ComputeIncrementalStrainBase, and ComputeCosseratIncrementalSmallStrain.

Definition at line 104 of file ComputeStrainBase.C.

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

◆ validParams()

InputParameters ComputeAxisymmetric1DSmallStrain::validParams ( )

static

Definition at line 17 of file ComputeAxisymmetric1DSmallStrain.C.

 {
   InputParameters params = Compute1DSmallStrain::validParams();
   params.addClassDescription("Compute a small strain 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;
 }