GlobalStrain Class Reference

#include <GlobalStrain.h>

Inheritance diagram for GlobalStrain:

Public Member Functions
	GlobalStrain (const InputParameters &parameters)
virtual void	reinit ()
virtual void	computeResidual ()
virtual void	computeJacobian ()
virtual void	computeOffDiagJacobian (unsigned int)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual void	assignComponentIndices (Order var_order)

Protected Attributes
const GlobalStrainUserObjectInterface &	_pst
const RankTwoTensor &	_pst_residual
const RankFourTensor &	_pst_jacobian
const VectorValue< bool > &	_periodic_dir
std::vector< std::pair< unsigned int, unsigned int > >	_components
const unsigned int	_dim

Detailed Description

Definition at line 27 of file GlobalStrain.h.

Constructor & Destructor Documentation

GlobalStrain()

GlobalStrain::GlobalStrain

(

const InputParameters &

parameters

)

Definition at line 35 of file GlobalStrain.C.

  : ScalarKernel(parameters),
    _pst(getUserObject<GlobalStrainUserObjectInterface>("global_strain_uo")),
    _pst_residual(_pst.getResidual()),
    _pst_jacobian(_pst.getJacobian()),
    _periodic_dir(_pst.getPeriodicDirections()),
    _components(_var.order()),
    _dim(_mesh.dimension())
{
  if ((_dim == 1 && _var.order() != FIRST) || (_dim == 2 && _var.order() != THIRD) ||
      (_dim == 3 && _var.order() != SIXTH))
    mooseError("PerdiodicStrain ScalarKernel is only compatible with scalar variables of order "
               "FIRST in 1D, THIRD in 2D, and SIXTH in 3D. Please change the order of the scalar"
               "variable according to the mesh dimension.");
 
  assignComponentIndices(_var.order());
}

Member Function Documentation

assignComponentIndices()

void GlobalStrain::assignComponentIndices ( Order  var_order )

protectedvirtual

Definition at line 83 of file GlobalStrain.C.

{
  switch (order)
  {
    case 1:
      _components[0].first = 0;
      _components[0].second = 0;
      break;
 
    case 3:
      _components[0].first = 0;
      _components[0].second = 0;
      _components[1].first = 1;
      _components[1].second = 1;
      _components[2].first = 0;
      _components[2].second = 1;
      break;
 
    case 6:
      _components[0].first = 0;
      _components[0].second = 0;
      _components[1].first = 1;
      _components[1].second = 1;
      _components[2].first = 2;
      _components[2].second = 2;
      _components[3].first = 1;
      _components[3].second = 2;
      _components[4].first = 0;
      _components[4].second = 2;
      _components[5].first = 0;
      _components[5].second = 1;
      break;
 
    default:
      mooseError("PerdiodicStrain ScalarKernel is only compatible with FIRST, THIRD, and SIXTH "
                 "order scalar variables.");
  }
}

Referenced by GlobalStrain().

computeJacobian()

void GlobalStrain::computeJacobian ( )

virtual

Definition at line 65 of file GlobalStrain.C.

{
  DenseMatrix<Number> & ke = _assembly.jacobianBlock(_var.number(), _var.number());
  for (_i = 0; _i < ke.m(); ++_i)
    for (_j = 0; _j < ke.m(); ++_j)
      // periodic direction check is not done for jacobian calculations to avoid zero pivot error
      ke(_i, _j) += _pst_jacobian(_components[_i].first,
                                  _components[_i].second,
                                  _components[_j].first,
                                  _components[_j].second);
}

computeOffDiagJacobian()

void GlobalStrain::computeOffDiagJacobian ( unsigned int  )

virtual

Definition at line 78 of file GlobalStrain.C.

{
}

computeResidual()

void GlobalStrain::computeResidual ( )

virtual

Definition at line 54 of file GlobalStrain.C.

{
  DenseVector<Number> & re = _assembly.residualBlock(_var.number());
  for (_i = 0; _i < re.size(); ++_i)
  {
    if (_periodic_dir(_components[_i].first) || _periodic_dir(_components[_i].second))
      re(_i) += _pst_residual(_components[_i].first, _components[_i].second);
  }
}

reinit()

virtual void GlobalStrain::reinit ( )

inlinevirtual

Definition at line 34 of file GlobalStrain.h.

{};

validParams()

InputParameters GlobalStrain::validParams ( )

static

Definition at line 25 of file GlobalStrain.C.

{
  InputParameters params = ScalarKernel::validParams();
  params.addClassDescription("Scalar Kernel to solve for the global strain");
  params.addRequiredParam<UserObjectName>("global_strain_uo",
                                          "The name of the GlobalStrainUserObject");
 
  return params;
}

Member Data Documentation

_components

std::vector<std::pair<unsigned int, unsigned int> > GlobalStrain::_components

protected

Definition at line 47 of file GlobalStrain.h.

Referenced by assignComponentIndices(), computeJacobian(), and computeResidual().

_dim

const unsigned int GlobalStrain::_dim

protected

Definition at line 48 of file GlobalStrain.h.

Referenced by GlobalStrain().

_periodic_dir

const VectorValue<bool>& GlobalStrain::_periodic_dir

protected

Definition at line 45 of file GlobalStrain.h.

Referenced by computeResidual().

_pst

const GlobalStrainUserObjectInterface& GlobalStrain::_pst

protected

Definition at line 42 of file GlobalStrain.h.

_pst_jacobian

const RankFourTensor& GlobalStrain::_pst_jacobian

protected

Definition at line 44 of file GlobalStrain.h.

Referenced by computeJacobian().

_pst_residual

const RankTwoTensor& GlobalStrain::_pst_residual

protected

Definition at line 43 of file GlobalStrain.h.

Referenced by computeResidual().

---

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

• tensor_mechanics/include/scalarkernels/GlobalStrain.h
• tensor_mechanics/src/scalarkernels/GlobalStrain.C