GlobalStrainUserObject Class Reference

#include <GlobalStrainUserObject.h>

Inheritance diagram for GlobalStrainUserObject:

Public Member Functions
	GlobalStrainUserObject (const InputParameters &parameters)
void	initialize () override
void	execute () override
void	threadJoin (const UserObject &uo) override
void	finalize () override
virtual const RankTwoTensor &	getResidual () const override
virtual const RankFourTensor &	getJacobian () const override
virtual const VectorValue< bool > &	getPeriodicDirections () const override
virtual void	computeAdditionalStress ()
	Calculate additional applied stresses. More...

Static Public Member Functions
static InputParameters	validParams ()

Protected Attributes
const std::string	_base_name
const MaterialProperty< RankFourTensor > &	_dstress_dstrain
const MaterialProperty< RankTwoTensor > &	_stress
RankTwoTensor	_applied_stress_tensor
RankTwoTensor	_residual
RankFourTensor	_jacobian
const unsigned int	_dim
const unsigned int	_ndisp
std::vector< unsigned int >	_disp_var
VectorValue< bool >	_periodic_dir

Detailed Description

Definition at line 23 of file GlobalStrainUserObject.h.

Constructor & Destructor Documentation

GlobalStrainUserObject()

GlobalStrainUserObject::GlobalStrainUserObject

(

const InputParameters &

parameters

)

Definition at line 34 of file GlobalStrainUserObject.C.

  : ElementUserObject(parameters),
    GlobalStrainUserObjectInterface(),
    _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
    _dstress_dstrain(getMaterialProperty<RankFourTensor>(_base_name + "Jacobian_mult")),
    _stress(getMaterialProperty<RankTwoTensor>(_base_name + "stress")),
    _dim(_mesh.dimension()),
    _ndisp(coupledComponents("displacements")),
    _disp_var(_ndisp),
    _periodic_dir()
{
  for (unsigned int i = 0; i < _ndisp; ++i)
    _disp_var[i] = coupled("displacements", i);
 
  for (unsigned int dir = 0; dir < _dim; ++dir)
  {
    _periodic_dir(dir) = _mesh.isTranslatedPeriodic(_disp_var[0], dir);
 
    for (unsigned int i = 1; i < _ndisp; ++i)
      if (_mesh.isTranslatedPeriodic(_disp_var[i], dir) != _periodic_dir(dir))
        mooseError("All the displacement components in a particular direction should have same "
                   "periodicity.");
  }
 
  if (isParamValid("applied_stress_tensor"))
    _applied_stress_tensor.fillFromInputVector(
        getParam<std::vector<Real>>("applied_stress_tensor"));
  else
    _applied_stress_tensor.zero();
}

Member Function Documentation

computeAdditionalStress()

virtual void GlobalStrainUserObject::computeAdditionalStress ( )

inlinevirtual

Calculate additional applied stresses.

Definition at line 41 of file GlobalStrainUserObject.h.

{};

Referenced by execute().

execute()

void GlobalStrainUserObject::execute ( )

override

Definition at line 73 of file GlobalStrainUserObject.C.

{
  computeAdditionalStress();
 
  for (unsigned int _qp = 0; _qp < _qrule->n_points(); _qp++)
  {
    // residual, integral of stress components
    _residual += _JxW[_qp] * _coord[_qp] * (_stress[_qp] - _applied_stress_tensor);
 
    // diagonal jacobian, integral of elasticity tensor components
    _jacobian += _JxW[_qp] * _coord[_qp] * _dstress_dstrain[_qp];
  }
}

finalize()

void GlobalStrainUserObject::finalize ( )

override

Definition at line 96 of file GlobalStrainUserObject.C.

{
  std::vector<Real> residual(9);
  std::vector<Real> jacobian(81);
 
  std::copy(&_residual(0, 0), &_residual(0, 0) + 9, residual.begin());
  std::copy(&_jacobian(0, 0, 0, 0), &_jacobian(0, 0, 0, 0) + 81, jacobian.begin());
 
  gatherSum(residual);
  gatherSum(jacobian);
 
  std::copy(residual.begin(), residual.end(), &_residual(0, 0));
  std::copy(jacobian.begin(), jacobian.end(), &_jacobian(0, 0, 0, 0));
}

getJacobian()

const RankFourTensor & GlobalStrainUserObject::getJacobian ( ) const

overridevirtual

Implements GlobalStrainUserObjectInterface.

Definition at line 118 of file GlobalStrainUserObject.C.

{
  return _jacobian;
}

getPeriodicDirections()

const VectorValue< bool > & GlobalStrainUserObject::getPeriodicDirections ( ) const

overridevirtual

Implements GlobalStrainUserObjectInterface.

Definition at line 124 of file GlobalStrainUserObject.C.

{
  return _periodic_dir;
}

getResidual()

const RankTwoTensor & GlobalStrainUserObject::getResidual ( ) const

overridevirtual

Implements GlobalStrainUserObjectInterface.

Definition at line 112 of file GlobalStrainUserObject.C.

{
  return _residual;
}

initialize()

void GlobalStrainUserObject::initialize ( )

override

Definition at line 66 of file GlobalStrainUserObject.C.

{
  _residual.zero();
  _jacobian.zero();
}

threadJoin()

void GlobalStrainUserObject::threadJoin ( const UserObject &  uo )

override

Definition at line 88 of file GlobalStrainUserObject.C.

{
  const GlobalStrainUserObject & pstuo = static_cast<const GlobalStrainUserObject &>(uo);
  _residual += pstuo._residual;
  _jacobian += pstuo._jacobian;
}

validParams()

InputParameters GlobalStrainUserObject::validParams ( )

static

Definition at line 19 of file GlobalStrainUserObject.C.

{
  InputParameters params = ElementUserObject::validParams();
  params.addClassDescription(
      "Global Strain UserObject to provide Residual and diagonal Jacobian entry");
  params.addParam<std::vector<Real>>("applied_stress_tensor",
                                     "Vector of values defining the constant applied stress "
                                     "to add, in order 11, 22, 33, 23, 13, 12");
  params.addParam<std::string>("base_name", "Material properties base name");
  params.addCoupledVar("displacements", "The name of the displacement variables");
  params.set<ExecFlagEnum>("execute_on") = EXEC_LINEAR;
 
  return params;
}

Member Data Documentation

_applied_stress_tensor

RankTwoTensor GlobalStrainUserObject::_applied_stress_tensor

protected

Definition at line 49 of file GlobalStrainUserObject.h.

Referenced by execute(), and GlobalStrainUserObject().

_base_name

const std::string GlobalStrainUserObject::_base_name

protected

Definition at line 41 of file GlobalStrainUserObject.h.

_dim

const unsigned int GlobalStrainUserObject::_dim

protected

Definition at line 53 of file GlobalStrainUserObject.h.

Referenced by GlobalStrainUserObject().

_disp_var

std::vector<unsigned int> GlobalStrainUserObject::_disp_var

protected

Definition at line 55 of file GlobalStrainUserObject.h.

Referenced by GlobalStrainUserObject().

_dstress_dstrain

const MaterialProperty<RankFourTensor>& GlobalStrainUserObject::_dstress_dstrain

protected

Definition at line 46 of file GlobalStrainUserObject.h.

Referenced by execute().

_jacobian

RankFourTensor GlobalStrainUserObject::_jacobian

protected

Definition at line 51 of file GlobalStrainUserObject.h.

Referenced by execute(), finalize(), getJacobian(), initialize(), and threadJoin().

_ndisp

const unsigned int GlobalStrainUserObject::_ndisp

protected

Definition at line 54 of file GlobalStrainUserObject.h.

Referenced by GlobalStrainUserObject().

_periodic_dir

VectorValue<bool> GlobalStrainUserObject::_periodic_dir

protected

Definition at line 56 of file GlobalStrainUserObject.h.

Referenced by getPeriodicDirections(), and GlobalStrainUserObject().

_residual

RankTwoTensor GlobalStrainUserObject::_residual

protected

Definition at line 50 of file GlobalStrainUserObject.h.

Referenced by execute(), finalize(), getResidual(), initialize(), and threadJoin().

_stress

const MaterialProperty<RankTwoTensor>& GlobalStrainUserObject::_stress

protected

Definition at line 47 of file GlobalStrainUserObject.h.

Referenced by execute().

---

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

• tensor_mechanics/include/userobjects/GlobalStrainUserObject.h
• tensor_mechanics/src/userobjects/GlobalStrainUserObject.C