#include <WeakPlaneStress.h>

Inheritance diagram for WeakPlaneStress:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	WeakPlaneStress (const InputParameters &parameters)

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
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Protected Attributes
const std::string	_base_name
	Base name of the material system that this kernel applies to.

const MaterialProperty< RankTwoTensor > &	_stress
	The stress tensor that provides the out-of-plane stress.

const MaterialProperty< RankFourTensor > &	_Jacobian_mult

const unsigned int	_direction
	The direction of the out-of-plane strain variable.

const bool	_disp_coupled
	Coupled displacement variables.

unsigned int	_ndisp
	Number of displacement variables.

std::vector< unsigned int >	_disp_var
	Variable numbers of the displacement variables.

const bool	_temp_coupled

const unsigned int	_temp_var

std::vector< const MaterialProperty< RankTwoTensor > * >	_deigenstrain_dT
	d(strain)/d(temperature), if computed by ComputeThermalExpansionEigenstrain

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

Definition at line 17 of file WeakPlaneStress.h.

Constructor & Destructor Documentation

◆ WeakPlaneStress()

WeakPlaneStress::WeakPlaneStress ( const InputParameters & parameters )

Definition at line 48 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(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)
{
  if (_disp_coupled)
    for (unsigned int i = 0; i < _ndisp; ++i)
      _disp_var[i] = coupled("displacements", i);
 
  if (_temp_coupled)
  {
    for (auto eigenstrain_name : getParam<std::vector<MaterialPropertyName>>("eigenstrain_names"))
      _deigenstrain_dT.push_back(&getMaterialPropertyDerivative<RankTwoTensor>(
          eigenstrain_name, coupledName("temperature", 0)));
  }
 
  // 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 83 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 90 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)
  {
    RankTwoTensor total_deigenstrain_dT;
    for (const auto deigenstrain_dT : _deigenstrain_dT)
      total_deigenstrain_dT += (*deigenstrain_dT)[_qp];
 
    Real sum = 0.0;
    for (unsigned int i = 0; i < 3; ++i)
      sum += _Jacobian_mult[_qp](_direction, _direction, i, i) * total_deigenstrain_dT(i, i);
 
    val = -sum * _test[_i][_qp] * _phi[_j][_qp];
  }
 
  return val;
}

◆ computeQpResidual()

Real WeakPlaneStress::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 77 of file WeakPlaneStress.C.

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

◆ validParams()

InputParameters WeakPlaneStress::validParams ( )

static

Definition at line 21 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::vector<MaterialPropertyName>>(
      "eigenstrain_names",
      {},
      "List of eigenstrains used in the strain calculation. Used for computing their derivaties "
      "for off-diagonal Jacobian terms.");
  params.addParam<std::string>("base_name", "Material property base name");
 
  MooseEnum direction("x y z", "z");
  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

Base name of the material system that this kernel applies to.

Definition at line 30 of file WeakPlaneStress.h.

◆ _deigenstrain_dT

std::vector<const MaterialProperty<RankTwoTensor> *> WeakPlaneStress::_deigenstrain_dT

protected

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

Definition at line 52 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian(), and WeakPlaneStress().

◆ _direction

const unsigned int WeakPlaneStress::_direction

protected

The direction of the out-of-plane strain variable.

Definition at line 37 of file WeakPlaneStress.h.

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

◆ _disp_coupled

const bool WeakPlaneStress::_disp_coupled

protected

Coupled displacement variables.

Definition at line 40 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian(), and WeakPlaneStress().

◆ _disp_var

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

protected

Variable numbers of the displacement variables.

Definition at line 46 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian(), and WeakPlaneStress().

◆ _Jacobian_mult

const MaterialProperty<RankFourTensor>& WeakPlaneStress::_Jacobian_mult

protected

Definition at line 34 of file WeakPlaneStress.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ _ndisp

unsigned int WeakPlaneStress::_ndisp

protected

Number of displacement variables.

Definition at line 43 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian(), and WeakPlaneStress().

◆ _stress

const MaterialProperty<RankTwoTensor>& WeakPlaneStress::_stress

protected

The stress tensor that provides the out-of-plane stress.

Definition at line 33 of file WeakPlaneStress.h.

Referenced by computeQpResidual().

◆ _temp_coupled

const bool WeakPlaneStress::_temp_coupled

protected

Definition at line 48 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian(), and WeakPlaneStress().

◆ _temp_var

const unsigned int WeakPlaneStress::_temp_var

protected

Definition at line 49 of file WeakPlaneStress.h.

Referenced by computeQpOffDiagJacobian().

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

solid_mechanics/include/kernels/WeakPlaneStress.h
solid_mechanics/src/kernels/WeakPlaneStress.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ WeakPlaneStress()

Member Function Documentation

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ validParams()

Member Data Documentation

◆ _base_name

◆ _deigenstrain_dT

◆ _direction

◆ _disp_coupled

◆ _disp_var

◆ _Jacobian_mult

◆ _ndisp

◆ _stress

◆ _temp_coupled

◆ _temp_var