WeakPlaneStress Class Reference

#include <WeakPlaneStress.h>

Inheritance diagram for WeakPlaneStress:

Public Member Functions
	WeakPlaneStress (const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpResidual () override
virtual Real	computeQpJacobian () override
virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Protected Attributes
const std::string	_base_name
const MaterialProperty< RankTwoTensor > &	_stress
const MaterialProperty< RankFourTensor > &	_Jacobian_mult
const unsigned int	_direction
const bool	_disp_coupled
	Coupled displacement variables. More...
unsigned int	_ndisp
std::vector< unsigned int >	_disp_var
const bool	_temp_coupled
const unsigned int	_temp_var
const MaterialProperty< RankTwoTensor > *const	_deigenstrain_dT
	d(strain)/d(temperature), if computed by ComputeThermalExpansionEigenstrain More...

Detailed Description

Definition at line 26 of file WeakPlaneStress.h.

Constructor & Destructor Documentation

WeakPlaneStress()

WeakPlaneStress::WeakPlaneStress

(

const InputParameters &

parameters

)

Definition at line 54 of file WeakPlaneStress.C.

  : DerivativeMaterialInterface<Kernel>(parameters),
    _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
    _stress(getMaterialProperty<RankTwoTensor>(_base_name + "stress")),
    _Jacobian_mult(getMaterialProperty<RankFourTensor>(_base_name + "Jacobian_mult")),
    _direction(parameters.isParamSetByUser("direction")
                   ? getParam<MooseEnum>("direction")
                   : getParam<MooseEnum>("out_of_plane_strain_direction")),
    _disp_coupled(isCoupled("displacements")),
    _ndisp(_disp_coupled ? coupledComponents("displacements") : 0),
    _disp_var(_ndisp),
    _temp_coupled(isCoupled("temperature")),
    _temp_var(_temp_coupled ? coupled("temperature") : 0),
    _deigenstrain_dT(_temp_coupled ? &getMaterialPropertyDerivative<RankTwoTensor>(
                                         getParam<std::string>("thermal_eigenstrain_name"),
                                         getVar("temperature", 0)->name())
                                   : nullptr)
{
  if (_disp_coupled)
    for (unsigned int i = 0; i < _ndisp; ++i)
      _disp_var[i] = coupled("displacements", i);
 
  if (parameters.isParamSetByUser("direction") &&
      parameters.isParamSetByUser("out_of_plane_strain_direction"))
    mooseError("Cannot specify both 'direction' and 'out_of_plane_strain_direction'! Use "
               "'out_of_plane_strain_direction'.");
 
  // Checking for consistency between mesh size and length of the provided displacements vector
  if (_disp_coupled && _ndisp != _mesh.dimension())
    mooseError("The number of displacement variables supplied must match the mesh dimension.");
}

Member Function Documentation

computeQpJacobian()

Real WeakPlaneStress::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 93 of file WeakPlaneStress.C.

{
  return _Jacobian_mult[_qp](_direction, _direction, _direction, _direction) * _test[_i][_qp] *
         _phi[_j][_qp];
}

computeQpOffDiagJacobian()

Real WeakPlaneStress::computeQpOffDiagJacobian ( unsigned int  jvar )

overrideprotectedvirtual

Definition at line 100 of file WeakPlaneStress.C.

{
  Real val = 0.0;
 
  // off-diagonal Jacobian with respect to a coupled displacement component
  if (_disp_coupled)
  {
    for (unsigned int coupled_direction = 0; coupled_direction < _ndisp; ++coupled_direction)
    {
      if (jvar == _disp_var[coupled_direction])
      {
        unsigned int coupled_direction_index = 0;
        switch (_direction)
        {
          case 0: // x
          {
            if (coupled_direction == 0)
              coupled_direction_index = 1;
            else
              coupled_direction_index = 2;
            break;
          }
          case 1: // y
          {
            if (coupled_direction == 0)
              coupled_direction_index = 0;
            else
              coupled_direction_index = 2;
            break;
          }
          default: // z
          {
            coupled_direction_index = coupled_direction;
            break;
          }
        }
 
        val = _Jacobian_mult[_qp](
                  _direction, _direction, coupled_direction_index, coupled_direction_index) *
              _test[_i][_qp] * _grad_phi[_j][_qp](coupled_direction_index);
      }
    }
  }
 
  // off-diagonal Jacobian with respect to a coupled temperature variable
  if (_temp_coupled && jvar == _temp_var)
  {
    Real sum = 0.0;
    for (unsigned int i = 0; i < 3; ++i)
      sum += _Jacobian_mult[_qp](_direction, _direction, i, i) * (*_deigenstrain_dT)[_qp](i, i);
 
    val = -sum * _test[_i][_qp] * _phi[_j][_qp];
  }
 
  return val;
}

computeQpResidual()

Real WeakPlaneStress::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 87 of file WeakPlaneStress.C.

{
  return _stress[_qp](_direction, _direction) * _test[_i][_qp];
}

validParams()

InputParameters WeakPlaneStress::validParams ( )

static

Definition at line 23 of file WeakPlaneStress.C.

{
  InputParameters params = Kernel::validParams();
  params.addClassDescription("Plane stress kernel to provide out-of-plane strain contribution.");
  params.addCoupledVar("displacements",
                       "The string of displacements suitable for the problem statement");
  params.addCoupledVar("temperature",
                       "The name of the temperature variable used in the "
                       "ComputeThermalExpansionEigenstrain.  (Not required for "
                       "simulations without temperature coupling.)");
  params.addParam<std::string>(
      "thermal_eigenstrain_name",
      "thermal_eigenstrain",
      "The eigenstrain_name used in the ComputeThermalExpansionEigenstrain.");
  params.addParam<std::string>("base_name", "Material property base name");
 
  MooseEnum direction("x y z", "z");
  params.addDeprecatedParam<MooseEnum>(
      "direction",
      direction,
      "The direction of the out-of-plane strain variable",
      "Use the new parameter name 'out_of_plane_strain_direction'");
  params.addParam<MooseEnum>("out_of_plane_strain_direction",
                             direction,
                             "The direction of the out-of-plane strain variable");
 
  params.set<bool>("use_displaced_mesh") = false;
 
  return params;
}

Member Data Documentation

_base_name

const std::string WeakPlaneStress::_base_name

protected

Definition at line 38 of file WeakPlaneStress.h.

_deigenstrain_dT

const MaterialProperty<RankTwoTensor>* const WeakPlaneStress::_deigenstrain_dT

protected

d(strain)/d(temperature), if computed by ComputeThermalExpansionEigenstrain

Definition at line 54 of file WeakPlaneStress.h.

_direction

const unsigned int WeakPlaneStress::_direction

protected

Definition at line 43 of file WeakPlaneStress.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

_disp_coupled

const bool WeakPlaneStress::_disp_coupled

protected

Coupled displacement variables.

Definition at line 46 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian(), and WeakPlaneStress().

_disp_var

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

protected

Definition at line 48 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian(), and WeakPlaneStress().

_Jacobian_mult

const MaterialProperty<RankFourTensor>& WeakPlaneStress::_Jacobian_mult

protected

Definition at line 41 of file WeakPlaneStress.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

_ndisp

unsigned int WeakPlaneStress::_ndisp

protected

Definition at line 47 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian(), and WeakPlaneStress().

_stress

const MaterialProperty<RankTwoTensor>& WeakPlaneStress::_stress

protected

Definition at line 40 of file WeakPlaneStress.h.

Referenced by computeQpResidual().

_temp_coupled

const bool WeakPlaneStress::_temp_coupled

protected

Definition at line 50 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian().

_temp_var

const unsigned int WeakPlaneStress::_temp_var

protected

Definition at line 51 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian().

---

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

• tensor_mechanics/include/kernels/WeakPlaneStress.h
• tensor_mechanics/src/kernels/WeakPlaneStress.C