GeneralizedPlaneStrainOffDiagOSPD Class Reference

Kernel class for coupled off diagonal Jacobian entries of ordinary state-based peridynamic generalized plane strain model. More...

#include <GeneralizedPlaneStrainOffDiagOSPD.h>

Inheritance diagram for GeneralizedPlaneStrainOffDiagOSPD:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	GeneralizedPlaneStrainOffDiagOSPD (const InputParameters &parameters)
virtual void	computeOffDiagJacobian (unsigned int jvar) override
virtual void	computeLocalOffDiagJacobian (unsigned int, unsigned int)
	Function to compute local contribution to the off-diagonal Jacobian at the current nodes. More...
virtual void	computePDNonlocalOffDiagJacobian (unsigned int, unsigned int)
	Function to compute nonlocal contribution to the off-diagonal Jacobian at the current nodes. More...
virtual void	initialSetup () override
virtual void	prepare () override
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 void	computeLocalResidual () override
virtual void	computeOffDiagJacobianScalar (unsigned int jvar_num) override
void	computeDispFullOffDiagJacobianScalar (unsigned int component, unsigned int jvar_num)
	Function to compute the full off diagonal Jacobian for coupling between displacements and scalar variable. More...
void	computeDispPartialOffDiagJacobianScalar (unsigned int component, unsigned int jvar_num)
	Function to compute partial off diagonal Jacobian for coupling between displacements and scalar variable. More...
void	computeTempOffDiagJacobianScalar (unsigned int jvar_num)
	Function to compute off disgonal Jacobian for coupling between temperature and scalar variable. More...

Protected Attributes
const MaterialProperty< RankFourTensor > &	_Cijkl
	Material point based material property. More...
const unsigned int	_scalar_out_of_plane_strain_var_num
	The variable number of the scalar out-of-plane strain variable. More...
DenseMatrix< Number >	_ken
DenseMatrix< Number >	_kne
std::vector< MooseVariable * >	_disp_var
	displacement variables More...
unsigned int	_ndisp
	number of displacement components More...
const std::vector< RealGradient > *	_orientation
	Vector of bond in current configuration. More...
std::vector< dof_id_type >	_ivardofs
	Current variable dof numbers for nodes i and j. More...
std::vector< Real >	_weights
	weights used for the current element to obtain the nodal stress More...
RealGradient	_current_vec
	Vector of bond in current configuration. More...
RealGradient	_current_unit_vec
	Unit vector of bond in current configuration. More...
const MaterialProperty< Real > &	_bond_local_dfdE
	Bond based material properties. More...
const MaterialProperty< Real > &	_bond_nonlocal_dfdE
const MaterialProperty< Real > &	_alpha
const bool	_temp_coupled
	Temperature variable. More...
MooseVariable *	_temp_var
const bool	_out_of_plane_strain_coupled
	Parameters for out-of-plane strain in weak plane stress formulation. More...
MooseVariable *	_out_of_plane_strain_var

Detailed Description

Kernel class for coupled off diagonal Jacobian entries of ordinary state-based peridynamic generalized plane strain model.

Definition at line 18 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Constructor & Destructor Documentation

GeneralizedPlaneStrainOffDiagOSPD()

GeneralizedPlaneStrainOffDiagOSPD::GeneralizedPlaneStrainOffDiagOSPD

(

const InputParameters &

parameters

)

Definition at line 31 of file GeneralizedPlaneStrainOffDiagOSPD.C.

  : MechanicsBasePD(parameters),
    _bond_local_dfdE(getMaterialProperty<Real>("bond_local_dfdE")),
    _bond_nonlocal_dfdE(getMaterialProperty<Real>("bond_nonlocal_dfdE")),
    _alpha(getMaterialProperty<Real>("thermal_expansion_coeff")),
    _Cijkl(getMaterialProperty<RankFourTensor>("elasticity_tensor")),
    _scalar_out_of_plane_strain_var_num(coupledScalar("scalar_out_of_plane_strain"))
{
  // Consistency check
  if (_disp_var.size() != 2)
    mooseError("GeneralizedPlaneStrain only works for two dimensional case!");
}

Member Function Documentation

computeDispFullOffDiagJacobianScalar()

void GeneralizedPlaneStrainOffDiagOSPD::computeDispFullOffDiagJacobianScalar ( unsigned int  component, unsigned int  jvar_num  )

protected

Function to compute the full off diagonal Jacobian for coupling between displacements and scalar variable.

Parameters

component	The index of displacement component
jvar_num	The coupled scalar variable number

Definition at line 68 of file GeneralizedPlaneStrainOffDiagOSPD.C.

Referenced by computeOffDiagJacobianScalar().

{
  prepareMatrixTag(_assembly, _var.number(), jvar_num, _ken);
  prepareMatrixTag(_assembly, jvar_num, _var.number(), _kne);
  MooseVariableScalar & jvar = _sys.getScalarVariable(_tid, jvar_num);

  // LOCAL jacobian contribution
  // fill in the column corresponding to the scalar variable from bond ij
  for (unsigned int i = 0; i < _nnodes; ++i)
    for (unsigned int j = 0; j < jvar.order(); ++j)
      _ken(i, j) +=
          (i == j ? -1 : 1) * _current_unit_vec(component) * _bond_local_dfdE[0] * _bond_status;

  // NONLOCAL jacobian contribution
  std::vector<RankTwoTensor> shape(_nnodes), dgrad(_nnodes);

  // for off-diagonal low components
  if (_bond_status > 0.5)
  {
    for (unsigned int nd = 0; nd < _nnodes; nd++)
    {
      if (_dim == 2)
        shape[nd](2, 2) = dgrad[nd](2, 2) = 1.0;
      // calculation of jacobian contribution to current node's neighbors
      std::vector<dof_id_type> ivardofs(_nnodes);
      ivardofs[nd] = _ivardofs[nd];
      std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_elem->node_id(nd));
      std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));

      Real vol_nb, weight_nb;
      RealGradient origin_vec_nb, current_vec_nb;

      for (unsigned int nb = 0; nb < neighbors.size(); nb++)
        if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[nb])) > 0.5)
        {
          ivardofs[1 - nd] =
              _pdmesh.nodePtr(neighbors[nb])->dof_number(_sys.number(), _var.number(), 0);
          vol_nb = _pdmesh.getNodeVolume(neighbors[nb]);

          // obtain bond nb's origin length and current orientation
          origin_vec_nb = _pdmesh.getNodeCoord(neighbors[nb]) -
                          _pdmesh.getNodeCoord(_current_elem->node_id(nd));

          for (unsigned int k = 0; k < _dim; k++)
            current_vec_nb(k) = origin_vec_nb(k) +
                                _disp_var[k]->getNodalValue(*_pdmesh.nodePtr(neighbors[nb])) -
                                _disp_var[k]->getNodalValue(*_current_elem->node_ptr(nd));

          weight_nb = _horizon_radius[nd] / origin_vec_nb.norm();
          // prepare shape tensor and deformation gradient tensor for current node
          for (unsigned int k = 0; k < _dim; k++)
            for (unsigned int l = 0; l < _dim; l++)
            {
              shape[nd](k, l) += weight_nb * origin_vec_nb(k) * origin_vec_nb(l) * vol_nb;
              dgrad[nd](k, l) += weight_nb * current_vec_nb(k) * origin_vec_nb(l) * vol_nb;
            }

          // cache the nonlocal jacobian contribution
          _local_ke.resize(_ken.m(), _ken.n());
          _local_ke.zero();
          _local_ke(0, 0) = (nd == 0 ? -1 : 1) * current_vec_nb(component) / current_vec_nb.norm() *
                            _bond_nonlocal_dfdE[nd] / origin_vec_nb.norm() * vol_nb * _bond_status;
          _local_ke(1, 0) = (nd == 0 ? 1 : -1) * current_vec_nb(component) / current_vec_nb.norm() *
                            _bond_nonlocal_dfdE[nd] / origin_vec_nb.norm() * vol_nb * _bond_status;

          addJacobian(_assembly, _local_ke, ivardofs, jvar.dofIndices(), _var.scalingFactor());
        }

      // finalize the shape tensor and deformation gradient tensor for node_i
      if (shape[nd].det() == 0.)
        mooseError("Singular shape tensor is detected in GeneralizedPlaneStrainOffDiagOSPD! Use "
                   "SingularShapeTensorEliminatorUserObjectPD to avoid singular shape tensor");

      shape[nd] = shape[nd].inverse();
      dgrad[nd] = dgrad[nd] * shape[nd];
    }
  }

  // off-diagonal jacobian entries on the row
  Real dEidUi =
      -_node_vol[1] * _horizon_radius[0] / _origin_vec.norm() *
      (_Cijkl[0](2, 2, 0, 0) * (_origin_vec(0) * shape[0](0, 0) + _origin_vec(1) * shape[0](1, 0)) *
           dgrad[0](component, 0) +
       _Cijkl[0](2, 2, 1, 1) * (_origin_vec(0) * shape[0](0, 1) + _origin_vec(1) * shape[0](1, 1)) *
           dgrad[0](component, 1));
  Real dEjdUj =
      _node_vol[0] * _horizon_radius[1] / _origin_vec.norm() *
      (_Cijkl[0](2, 2, 0, 0) * (_origin_vec(0) * shape[1](0, 0) + _origin_vec(1) * shape[1](1, 0)) *
           dgrad[1](component, 0) +
       _Cijkl[0](2, 2, 1, 1) * (_origin_vec(0) * shape[1](0, 1) + _origin_vec(1) * shape[1](1, 1)) *
           dgrad[1](component, 1));

  // fill in the row corresponding to the scalar variable
  _kne(0, 0) += (dEidUi * _node_vol[0] - dEjdUj * _node_vol[1]) * _bond_status; // node i
  _kne(0, 1) += (dEjdUj * _node_vol[1] - dEidUi * _node_vol[0]) * _bond_status; // node j

  accumulateTaggedLocalMatrix(_assembly, _var.number(), jvar_num, _ken);
  accumulateTaggedLocalMatrix(_assembly, jvar_num, _var.number(), _kne);
}

computeDispPartialOffDiagJacobianScalar()

void GeneralizedPlaneStrainOffDiagOSPD::computeDispPartialOffDiagJacobianScalar ( unsigned int  component, unsigned int  jvar_num  )

protected

Function to compute partial off diagonal Jacobian for coupling between displacements and scalar variable.

Parameters

component	The index of displacement component
jvar_num	The coupled scalar variable number

Definition at line 170 of file GeneralizedPlaneStrainOffDiagOSPD.C.

Referenced by computeOffDiagJacobianScalar().

{
  prepareMatrixTag(_assembly, _var.number(), jvar_num, _ken);
  prepareMatrixTag(_assembly, jvar_num, _var.number(), _kne);
  MooseVariableScalar & jvar = _sys.getScalarVariable(_tid, jvar_num);

  // fill in the column corresponding to the scalar variable from bond ij
  for (unsigned int i = 0; i < _nnodes; ++i)
    for (unsigned int j = 0; j < jvar.order(); ++j)
    {
      _ken(i, j) +=
          (i == j ? -1 : 1) * _current_unit_vec(component) * _bond_local_dfdE[0] * _bond_status;
      _kne(j, i) +=
          (i == j ? -1 : 1) * _current_unit_vec(component) * _bond_local_dfdE[0] * _bond_status;
    }

  accumulateTaggedLocalMatrix(_assembly, _var.number(), jvar_num, _ken);
  accumulateTaggedLocalMatrix(_assembly, jvar_num, _var.number(), _kne);
}

computeLocalOffDiagJacobian()

virtual void MechanicsBasePD::computeLocalOffDiagJacobian ( unsigned int  , unsigned int    )

inlinevirtualinherited

Function to compute local contribution to the off-diagonal Jacobian at the current nodes.

Parameters

coupled_component

The coupled variable number

Reimplemented in HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD, HorizonStabilizedFormIISmallStrainMechanicsNOSPD, HorizonStabilizedFormIFiniteStrainMechanicsNOSPD, HorizonStabilizedFormISmallStrainMechanicsNOSPD, MechanicsOSPD, ForceStabilizedSmallStrainMechanicsNOSPD, MechanicsBPD, and WeakPlaneStressNOSPD.

Definition at line 31 of file MechanicsBasePD.h.

Referenced by MechanicsBasePD::computeOffDiagJacobian().

                                                          {};

computeLocalResidual()

virtual void GeneralizedPlaneStrainOffDiagOSPD::computeLocalResidual ( )

inlineoverrideprotectedvirtual

Definition at line 26 of file GeneralizedPlaneStrainOffDiagOSPD.h.

{};

computeOffDiagJacobian()

void MechanicsBasePD::computeOffDiagJacobian ( unsigned int  jvar )

overridevirtualinherited

Definition at line 73 of file MechanicsBasePD.C.

{
  prepare();

  if (jvar_num == _var.number())
    computeJacobian();
  else
  {
    unsigned int coupled_component = 0;
    bool active = false;

    for (unsigned int i = 0; i < _dim; ++i)
      if (jvar_num == _disp_var[i]->number())
      {
        coupled_component = i;
        active = true;
      }

    if (_temp_coupled && jvar_num == _temp_var->number())
      active = true;

    if (_out_of_plane_strain_coupled && jvar_num == _out_of_plane_strain_var->number())
      active = true;

    if (active)
    {
      prepareMatrixTag(_assembly, _var.number(), jvar_num);
      computeLocalOffDiagJacobian(jvar_num, coupled_component);
      accumulateTaggedLocalMatrix();

      if (_use_full_jacobian)
        computePDNonlocalOffDiagJacobian(jvar_num, coupled_component);
    }
  }
}

computeOffDiagJacobianScalar()

void GeneralizedPlaneStrainOffDiagOSPD::computeOffDiagJacobianScalar ( unsigned int  jvar_num )

overrideprotectedvirtual

Definition at line 46 of file GeneralizedPlaneStrainOffDiagOSPD.C.

{
  if (jvar_num == _scalar_out_of_plane_strain_var_num)
  {
    prepare();

    if (_var.number() == _disp_var[0]->number())
      if (_use_full_jacobian)
        computeDispFullOffDiagJacobianScalar(0, jvar_num);
      else
        computeDispPartialOffDiagJacobianScalar(0, jvar_num);
    else if (_var.number() == _disp_var[1]->number())
      if (_use_full_jacobian)
        computeDispFullOffDiagJacobianScalar(1, jvar_num);
      else
        computeDispPartialOffDiagJacobianScalar(1, jvar_num);
    else if (_temp_coupled ? _var.number() == _temp_var->number() : 0)
      computeTempOffDiagJacobianScalar(jvar_num);
  }
}

computePDNonlocalOffDiagJacobian()

virtual void MechanicsBasePD::computePDNonlocalOffDiagJacobian ( unsigned int  , unsigned int    )

inlinevirtualinherited

Function to compute nonlocal contribution to the off-diagonal Jacobian at the current nodes.

Parameters

jvar_num	The number of the first coupled variable
coupled_component	The component number of the second coupled variable

Reimplemented in HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD, HorizonStabilizedFormIISmallStrainMechanicsNOSPD, HorizonStabilizedFormIFiniteStrainMechanicsNOSPD, HorizonStabilizedFormISmallStrainMechanicsNOSPD, MechanicsOSPD, ForceStabilizedSmallStrainMechanicsNOSPD, and WeakPlaneStressNOSPD.

Definition at line 39 of file MechanicsBasePD.h.

Referenced by MechanicsBasePD::computeOffDiagJacobian().

                                                               {};

computeTempOffDiagJacobianScalar()

void GeneralizedPlaneStrainOffDiagOSPD::computeTempOffDiagJacobianScalar ( unsigned int  jvar_num )

protected

Function to compute off disgonal Jacobian for coupling between temperature and scalar variable.

Parameters

jvar_num

The coupled scalar variable number

Definition at line 192 of file GeneralizedPlaneStrainOffDiagOSPD.C.

Referenced by computeOffDiagJacobianScalar().

{
  // off-diagonal jacobian entries on the row
  prepareMatrixTag(_assembly, jvar_num, _var.number(), _kne);

  // calculate number of active neighbors for node i and j
  std::vector<unsigned int> active_neighbors(_nnodes, 0);
  for (unsigned int nd = 0; nd < _nnodes; nd++)
  {
    std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_elem->node_id(nd));
    for (unsigned int nb = 0; nb < bonds.size(); ++nb)
      if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[nb])) > 0.5)
        active_neighbors[nd]++;

    if (active_neighbors[nd] == 0) // avoid dividing by zero
      active_neighbors[nd] = 1;
  }

  //  one-way coupling between the strain_zz and temperature. fill in the row corresponding to the
  //  scalar variable strain_zz
  _kne(0, 0) += -_alpha[0] *
                (_Cijkl[0](2, 2, 0, 0) + _Cijkl[0](2, 2, 1, 1) + _Cijkl[0](2, 2, 2, 2)) *
                _node_vol[0] / active_neighbors[0] * _bond_status; // node i
  _kne(0, 1) += -_alpha[0] *
                (_Cijkl[0](2, 2, 0, 0) + _Cijkl[0](2, 2, 1, 1) + _Cijkl[0](2, 2, 2, 2)) *
                _node_vol[1] / active_neighbors[1] * _bond_status; // node j

  accumulateTaggedLocalMatrix(_assembly, jvar_num, _var.number(), _kne);
}

initialSetup()

void MechanicsBasePD::initialSetup ( )

overridevirtualinherited

Definition at line 45 of file MechanicsBasePD.C.

{
  _orientation = &_assembly.getFE(FEType(), 1)->get_dxyzdxi();
}

prepare()

void MechanicsBasePD::prepare ( )

overridevirtualinherited

Definition at line 51 of file MechanicsBasePD.C.

Referenced by MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), and computeOffDiagJacobianScalar().

{
  PeridynamicsKernelBase::prepare();

  _ivardofs.resize(_nnodes);
  _weights.resize(_nnodes);
  for (unsigned int nd = 0; nd < _nnodes; ++nd)
  {
    _ivardofs[nd] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
    _weights[nd] = _pdmesh.getNeighborWeight(
        _current_elem->node_id(nd),
        _pdmesh.getNeighborIndex(_current_elem->node_id(nd), _current_elem->node_id(1 - nd)));
  }

  for (unsigned int i = 0; i < _dim; ++i)
    _current_vec(i) = _origin_vec(i) + _disp_var[i]->getNodalValue(*_current_elem->node_ptr(1)) -
                      _disp_var[i]->getNodalValue(*_current_elem->node_ptr(0));

  _current_unit_vec = _current_vec / _current_vec.norm();
}

validParams()

InputParameters GeneralizedPlaneStrainOffDiagOSPD::validParams ( )

static

Definition at line 17 of file GeneralizedPlaneStrainOffDiagOSPD.C.

{
  InputParameters params = MechanicsBasePD::validParams();
  params.addClassDescription(
      "Class for calculating the off-diagonal Jacobian corresponding to "
      "coupling between displacements (or temperature) and the scalar out-of-plane strain for "
      "the generalized plane strain using the OSPD formulation");

  params.addCoupledVar("scalar_out_of_plane_strain",
                       "Scalar variable for strain in the out-of-plane direction");

  return params;
}

Member Data Documentation

_alpha

const MaterialProperty<Real>& GeneralizedPlaneStrainOffDiagOSPD::_alpha

protected

Definition at line 54 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeTempOffDiagJacobianScalar().

_bond_local_dfdE

const MaterialProperty<Real>& GeneralizedPlaneStrainOffDiagOSPD::_bond_local_dfdE

protected

Bond based material properties.

Definition at line 52 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), and computeDispPartialOffDiagJacobianScalar().

_bond_nonlocal_dfdE

const MaterialProperty<Real>& GeneralizedPlaneStrainOffDiagOSPD::_bond_nonlocal_dfdE

protected

Definition at line 53 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar().

_Cijkl

const MaterialProperty<RankFourTensor>& GeneralizedPlaneStrainOffDiagOSPD::_Cijkl

protected

Material point based material property.

Definition at line 58 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), and computeTempOffDiagJacobianScalar().

_current_unit_vec

RealGradient MechanicsBasePD::_current_unit_vec

protectedinherited

Unit vector of bond in current configuration.

Definition at line 75 of file MechanicsBasePD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), computeDispPartialOffDiagJacobianScalar(), MechanicsBPD::computeLocalJacobian(), MechanicsOSPD::computeLocalJacobian(), MechanicsBPD::computeLocalOffDiagJacobian(), MechanicsOSPD::computeLocalOffDiagJacobian(), MechanicsBPD::computeLocalResidual(), MechanicsOSPD::computeLocalResidual(), MechanicsOSPD::computeNonlocalJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBasePD::prepare().

_current_vec

RealGradient MechanicsBasePD::_current_vec

protectedinherited

Vector of bond in current configuration.

Definition at line 72 of file MechanicsBasePD.h.

Referenced by MechanicsBPD::computeLocalJacobian(), MechanicsOSPD::computeLocalJacobian(), MechanicsBPD::computeLocalOffDiagJacobian(), MechanicsOSPD::computeLocalOffDiagJacobian(), and MechanicsBasePD::prepare().

_disp_var

std::vector<MooseVariable *> MechanicsBasePD::_disp_var

protectedinherited

displacement variables

Definition at line 47 of file MechanicsBasePD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalResidual(), MechanicsOSPD::computeNonlocalJacobian(), MechanicsOSPD::computeNonlocalResidual(), MechanicsBasePD::computeOffDiagJacobian(), computeOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::GeneralizedPlaneStrainOffDiagNOSPD(), GeneralizedPlaneStrainOffDiagOSPD(), MechanicsBasePD::MechanicsBasePD(), and MechanicsBasePD::prepare().

_ivardofs

std::vector<dof_id_type> MechanicsBasePD::_ivardofs

protectedinherited

Current variable dof numbers for nodes i and j.

Definition at line 66 of file MechanicsBasePD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeNonlocalJacobian(), MechanicsOSPD::computeNonlocalJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeNonlocalJacobian(), MechanicsOSPD::computeNonlocalResidual(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalResidual(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeNonlocalResidual(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBasePD::prepare().

_ken

DenseMatrix<Number> GeneralizedPlaneStrainOffDiagOSPD::_ken

protected

Definition at line 63 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), and computeDispPartialOffDiagJacobianScalar().

_kne

DenseMatrix<Number> GeneralizedPlaneStrainOffDiagOSPD::_kne

protected

Definition at line 64 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeDispFullOffDiagJacobianScalar(), computeDispPartialOffDiagJacobianScalar(), and computeTempOffDiagJacobianScalar().

_ndisp

unsigned int MechanicsBasePD::_ndisp

protectedinherited

number of displacement components

Definition at line 55 of file MechanicsBasePD.h.

Referenced by MechanicsBasePD::MechanicsBasePD().

_orientation

const std::vector<RealGradient>* MechanicsBasePD::_orientation

protectedinherited

Vector of bond in current configuration.

Definition at line 63 of file MechanicsBasePD.h.

Referenced by MechanicsBasePD::initialSetup().

_out_of_plane_strain_coupled

const bool MechanicsBasePD::_out_of_plane_strain_coupled

protectedinherited

Parameters for out-of-plane strain in weak plane stress formulation.

Definition at line 58 of file MechanicsBasePD.h.

Referenced by HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian().

_out_of_plane_strain_var

MooseVariable* MechanicsBasePD::_out_of_plane_strain_var

protectedinherited

Definition at line 59 of file MechanicsBasePD.h.

Referenced by HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian().

_scalar_out_of_plane_strain_var_num

const unsigned int GeneralizedPlaneStrainOffDiagOSPD::_scalar_out_of_plane_strain_var_num

protected

The variable number of the scalar out-of-plane strain variable.

Definition at line 61 of file GeneralizedPlaneStrainOffDiagOSPD.h.

Referenced by computeOffDiagJacobianScalar().

_temp_coupled

const bool MechanicsBasePD::_temp_coupled

protectedinherited

Temperature variable.

Definition at line 50 of file MechanicsBasePD.h.

Referenced by MechanicsBPD::computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), computeOffDiagJacobianScalar(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBaseNOSPD::MechanicsBaseNOSPD().

_temp_var

MooseVariable* MechanicsBasePD::_temp_var

protectedinherited

Definition at line 51 of file MechanicsBasePD.h.

Referenced by MechanicsBPD::computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), computeOffDiagJacobianScalar(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBaseNOSPD::MechanicsBaseNOSPD().

_weights

std::vector<Real> MechanicsBasePD::_weights

protectedinherited

weights used for the current element to obtain the nodal stress

Definition at line 69 of file MechanicsBasePD.h.

Referenced by WeakPlaneStressNOSPD::computeLocalJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalResidual(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBasePD::prepare().

---

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

• peridynamics/include/kernels/GeneralizedPlaneStrainOffDiagOSPD.h
• peridynamics/src/kernels/GeneralizedPlaneStrainOffDiagOSPD.C