Kernel class for Form I of the horizon-stabilized peridynamic correspondence model for small strain. More...

#include <HorizonStabilizedFormISmallStrainMechanicsNOSPD.h>

Inheritance diagram for HorizonStabilizedFormISmallStrainMechanicsNOSPD:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	HorizonStabilizedFormISmallStrainMechanicsNOSPD (const InputParameters &parameters)

virtual void	computeOffDiagJacobian (unsigned int jvar) override

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	computeLocalJacobian () override

virtual void	computeNonlocalJacobian () override

virtual void	computeLocalOffDiagJacobian (unsigned int jvar_num, unsigned int coupled_component) override
	Function to compute local contribution to the off-diagonal Jacobian at the current nodes.

virtual void	computePDNonlocalOffDiagJacobian (unsigned int jvar_num, unsigned int coupled_component) override
	Function to compute nonlocal contribution to the off-diagonal Jacobian at the current nodes.

virtual RankTwoTensor	computeDSDU (unsigned int component, unsigned int nd)
	Function to compute derivative of stress with respect to displacements for small strain problems.

Protected Attributes
const unsigned int	_component
	The index of displacement component.

std::vector< MooseVariable * >	_disp_var
	displacement variables

unsigned int	_ndisp
	number of displacement components

const std::vector< RealGradient > *	_orientation
	Vector of bond in current configuration.

std::vector< dof_id_type >	_ivardofs
	Current variable dof numbers for nodes i and j.

std::vector< Real >	_weights
	weights used for the current element to obtain the nodal stress

RealGradient	_current_vec
	Vector of bond in current configuration.

RealGradient	_current_unit_vec
	Unit vector of bond in current configuration.


const MaterialProperty< Real > &	_multi
	Material point based material properties.

const MaterialProperty< RankTwoTensor > &	_stress

const MaterialProperty< RankTwoTensor > &	_shape2

const MaterialProperty< RankTwoTensor > &	_dgrad

const MaterialProperty< RankTwoTensor > &	_ddgraddu

const MaterialProperty< RankTwoTensor > &	_ddgraddv

const MaterialProperty< RankTwoTensor > &	_ddgraddw

const MaterialProperty< RankFourTensor > &	_Jacobian_mult

const std::vector< MaterialPropertyName >	_eigenstrain_names

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


const bool	_temp_coupled
	Temperature variable.

MooseVariable *	_temp_var


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

MooseVariable *	_out_of_plane_strain_var

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

Kernel class for Form I of the horizon-stabilized peridynamic correspondence model for small strain.

Definition at line 20 of file HorizonStabilizedFormISmallStrainMechanicsNOSPD.h.

Constructor & Destructor Documentation

◆ HorizonStabilizedFormISmallStrainMechanicsNOSPD()

HorizonStabilizedFormISmallStrainMechanicsNOSPD::HorizonStabilizedFormISmallStrainMechanicsNOSPD ( const InputParameters & parameters )

Definition at line 30 of file HorizonStabilizedFormISmallStrainMechanicsNOSPD.C.

  : MechanicsBaseNOSPD(parameters), _component(getParam<unsigned int>("component"))
{
}

Member Function Documentation

◆ computeDSDU()

RankTwoTensor MechanicsBaseNOSPD::computeDSDU	(	unsigned int	component,
		unsigned int	nd
	)

protectedvirtualinherited

Function to compute derivative of stress with respect to displacements for small strain problems.

Parameters

component	The index of displacement component
nd	The local index of element node (either 1 or 2 for Edge2 element)

Returns: The calculated derivative

Reimplemented in MechanicsFiniteStrainBaseNOSPD.

Definition at line 47 of file MechanicsBaseNOSPD.C.

{
  // compute the derivative of stress w.r.t the solution components for small strain
  RankTwoTensor dSdU;
  if (component == 0)
    dSdU = _Jacobian_mult[nd] * 0.5 * (_ddgraddu[nd].transpose() + _ddgraddu[nd]);
  else if (component == 1)
    dSdU = _Jacobian_mult[nd] * 0.5 * (_ddgraddv[nd].transpose() + _ddgraddv[nd]);
  else if (component == 2)
    dSdU = _Jacobian_mult[nd] * 0.5 * (_ddgraddw[nd].transpose() + _ddgraddw[nd]);
 
  return dSdU;
}

Referenced by GeneralizedPlaneStrainOffDiagNOSPD::computeDispFullOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagNOSPD::computeDispPartialOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalJacobian(), computeLocalJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), and HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian().

◆ computeLocalJacobian()

void HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalJacobian ( )

overrideprotectedvirtual

Definition at line 53 of file HorizonStabilizedFormISmallStrainMechanicsNOSPD.C.

{
  // excludes dTi/dUj and dTj/dUi contributions, which were considered as nonlocal contribution
  for (unsigned int i = 0; i < _nnodes; ++i)
    for (unsigned int j = 0; j < _nnodes; ++j)
      _local_ke(i, j) += (i == 0 ? -1 : 1) * _multi[j] *
                         (computeDSDU(_component, j) * _shape2[j].inverse()).row(_component) *
                         _origin_vec * _bond_status;
}

◆ computeLocalOffDiagJacobian()

void HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian	(	unsigned int	,
		unsigned int
	)

overrideprotectedvirtual

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

Parameters

coupled_component The coupled variable number

Reimplemented from MechanicsBasePD.

Definition at line 138 of file HorizonStabilizedFormISmallStrainMechanicsNOSPD.C.

{
  _local_ke.zero();
  if (_temp_coupled && jvar_num == _temp_var->number()) // temperature is coupled
  {
    std::vector<RankTwoTensor> dSdT(_nnodes);
    for (unsigned int nd = 0; nd < _nnodes; ++nd)
      for (unsigned int es = 0; es < _deigenstrain_dT.size(); ++es)
        dSdT[nd] = -_Jacobian_mult[nd] * (*_deigenstrain_dT[es])[nd];
 
    for (unsigned int i = 0; i < _nnodes; ++i)
      for (unsigned int j = 0; j < _nnodes; ++j)
        _local_ke(i, j) += (i == 0 ? -1 : 1) * _multi[j] *
                           (dSdT[j] * _shape2[j].inverse()).row(_component) * _origin_vec *
                           _bond_status;
  }
  else if (_out_of_plane_strain_coupled &&
           jvar_num == _out_of_plane_strain_var
                           ->number()) // weak plane stress case, out_of_plane_strain is coupled
  {
    std::vector<RankTwoTensor> dSdE33(_nnodes);
    for (unsigned int nd = 0; nd < _nnodes; ++nd)
      for (unsigned int i = 0; i < 3; ++i)
        for (unsigned int j = 0; j < 3; ++j)
          dSdE33[nd](i, j) = _Jacobian_mult[nd](i, j, 2, 2);
 
    for (unsigned int i = 0; i < _nnodes; ++i)
      for (unsigned int j = 0; j < _nnodes; ++j)
        _local_ke(i, j) += (i == 0 ? -1 : 1) * _multi[j] *
                           (dSdE33[j] * _shape2[j].inverse()).row(_component) * _origin_vec *
                           _bond_status;
  }
  else // off-diagonal Jacobian with respect to other displacement variables
  {
    // ONLY consider the contributions to node i and j
    // contributions to their neighbors are considered as nonlocal off-diagonal
    for (unsigned int i = 0; i < _nnodes; ++i)
      for (unsigned int j = 0; j < _nnodes; ++j)
        _local_ke(i, j) +=
            (i == 0 ? -1 : 1) * _multi[j] *
            (computeDSDU(coupled_component, j) * _shape2[j].inverse()).row(_component) *
            _origin_vec * _bond_status;
  }
}

◆ computeLocalResidual()

void HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeLocalResidual ( )

overrideprotectedvirtual

Definition at line 37 of file HorizonStabilizedFormISmallStrainMechanicsNOSPD.C.

{
  // For small strain assumptions, stress measures, i.e., Cauchy stress (Sigma), the first
  // Piola-Kirchhoff stress (P), and the second Piola-Kirchhoff stress (S) are approximately the
  // same. Thus, the nodal force state tensors are calculated using the Cauchy stresses,
  // i.e., T = Sigma * inv(Shape) * xi * multi.
  // Cauchy stress is calculated as Sigma = C * E in the ComputeSmallStrainNOSPD material class.
 
  for (unsigned int nd = 0; nd < _nnodes; ++nd)
    for (unsigned int i = 0; i < _nnodes; ++i)
      _local_re(i) += (i == 0 ? -1 : 1) * _multi[nd] *
                      (_stress[nd] * _shape2[nd].inverse()).row(_component) * _origin_vec *
                      _bond_status;
}

◆ computeNonlocalJacobian()

void HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeNonlocalJacobian ( )

overrideprotectedvirtual

Definition at line 64 of file HorizonStabilizedFormISmallStrainMechanicsNOSPD.C.

{
  // includes dTi/dUj and dTj/dUi contributions
  // excludes contributions to node i and j, which were considered as local contributions
  for (unsigned int nd = 0; nd < _nnodes; ++nd)
  {
    // calculation of jacobian contribution to current_node's neighbors
    std::vector<dof_id_type> ivardofs(_nnodes);
    ivardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), _var.number(), 0);
    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));
 
    dof_id_type nb_index =
        std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - nd)) -
        neighbors.begin();
    std::vector<dof_id_type> dg_neighbors =
        _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(nd), nb_index);
 
    Real vol_nb;
    RealGradient origin_vec_nb;
    RankTwoTensor dFdUk, dPxdUkx;
 
    for (unsigned int nb = 0; nb < dg_neighbors.size(); ++nb)
      if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb]])) > 0.5)
      {
        ivardofs[1] = _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])
                          ->dof_number(_sys.number(), _var.number(), 0);
        vol_nb = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb]]);
 
        origin_vec_nb = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb]]) -
                        _pdmesh.getNodeCoord(_current_elem->node_id(nd));
 
        dFdUk.zero();
        for (unsigned int i = 0; i < _dim; ++i)
          dFdUk(_component, i) =
              _horizon_radius[nd] / origin_vec_nb.norm() * origin_vec_nb(i) * vol_nb;
 
        dFdUk *= _shape2[nd].inverse();
 
        dPxdUkx = _Jacobian_mult[nd] * 0.5 * (dFdUk.transpose() + dFdUk);
 
        for (unsigned int i = 0; i < _nnodes; ++i)
          for (unsigned int j = 0; j < _nnodes; ++j)
            _local_ke(i, j) = (i == 0 ? -1 : 1) * (j == 0 ? 0 : 1) * _multi[nd] *
                              (dPxdUkx * _shape2[nd].inverse()).row(_component) * _origin_vec *
                              _bond_status;
 
        addJacobian(_assembly, _local_ke, _ivardofs, ivardofs, _var.scalingFactor());
 
        if (_has_diag_save_in)
        {
          unsigned int rows = _nnodes;
          DenseVector<Real> diag(rows);
          for (unsigned int i = 0; i < rows; ++i)
            diag(i) = _local_ke(i, i);
 
          Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
          for (unsigned int i = 0; i < _diag_save_in.size(); ++i)
          {
            std::vector<dof_id_type> diag_save_in_dofs(2);
            diag_save_in_dofs[0] = _current_elem->node_ptr(nd)->dof_number(
                _diag_save_in[i]->sys().number(), _diag_save_in[i]->number(), 0);
            diag_save_in_dofs[1] =
                _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])
                    ->dof_number(_diag_save_in[i]->sys().number(), _diag_save_in[i]->number(), 0);
 
            _diag_save_in[i]->sys().solution().add_vector(diag, diag_save_in_dofs);
          }
        }
      }
  }
}

◆ 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);
    }
  }
}

◆ computePDNonlocalOffDiagJacobian()

void HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian	(	unsigned int	,
		unsigned int
	)

overrideprotectedvirtual

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 from MechanicsBasePD.

Definition at line 185 of file HorizonStabilizedFormISmallStrainMechanicsNOSPD.C.

{
  if (_temp_coupled && jvar_num == _temp_var->number())
  {
    // no nonlocal contribution from temperature
  }
  else if (_out_of_plane_strain_coupled && jvar_num == _out_of_plane_strain_var->number())
  {
    // no nonlocal contribution from out of plane strain
  }
  else
  {
    for (unsigned int nd = 0; nd < _nnodes; ++nd)
    {
      // calculation of jacobian contribution to current_node's neighbors
      // NOT including the contribution to nodes i and j, which is considered as local off-diagonal
      std::vector<dof_id_type> jvardofs(_nnodes);
      jvardofs[0] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), jvar_num, 0);
      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));
 
      dof_id_type nb_index =
          std::find(neighbors.begin(), neighbors.end(), _current_elem->node_id(1 - nd)) -
          neighbors.begin();
      std::vector<dof_id_type> dg_neighbors =
          _pdmesh.getBondDeformationGradientNeighbors(_current_elem->node_id(nd), nb_index);
 
      Real vol_nb;
      RealGradient origin_vec_nb;
      RankTwoTensor dFdUk, dPxdUky;
 
      for (unsigned int nb = 0; nb < dg_neighbors.size(); ++nb)
        if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[dg_neighbors[nb]])) > 0.5)
        {
          jvardofs[1] =
              _pdmesh.nodePtr(neighbors[dg_neighbors[nb]])->dof_number(_sys.number(), jvar_num, 0);
          vol_nb = _pdmesh.getNodeVolume(neighbors[dg_neighbors[nb]]);
 
          origin_vec_nb = _pdmesh.getNodeCoord(neighbors[dg_neighbors[nb]]) -
                          _pdmesh.getNodeCoord(_current_elem->node_id(nd));
 
          dFdUk.zero();
          for (unsigned int i = 0; i < _dim; ++i)
            dFdUk(coupled_component, i) =
                _horizon_radius[nd] / origin_vec_nb.norm() * origin_vec_nb(i) * vol_nb;
 
          dFdUk *= _shape2[nd].inverse();
 
          dPxdUky = _Jacobian_mult[nd] * 0.5 * (dFdUk.transpose() + dFdUk);
 
          for (unsigned int i = 0; i < _nnodes; ++i)
            for (unsigned int j = 0; j < _nnodes; ++j)
              _local_ke(i, j) = (i == 0 ? -1 : 1) * (j == 0 ? 0 : 1) * _multi[nd] *
                                (dPxdUky * _shape2[nd].inverse()).row(_component) * _origin_vec *
                                _bond_status;
 
          addJacobian(_assembly, _local_ke, _ivardofs, jvardofs, _var.scalingFactor());
        }
    }
  }
}

◆ 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.

{
  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();
}

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

◆ validParams()

InputParameters HorizonStabilizedFormISmallStrainMechanicsNOSPD::validParams ( )

static

Definition at line 16 of file HorizonStabilizedFormISmallStrainMechanicsNOSPD.C.

{
  InputParameters params = MechanicsBaseNOSPD::validParams();
  params.addClassDescription(
      "Class for calculating the residual and the Jacobian for Form I of the horizon-stabilized"
      "peridynamic correspondence model under small strain assumptions");
 
  params.addRequiredParam<unsigned int>(
      "component",
      "An integer corresponding to the variable this kernel acts on (0 for x, 1 for y, 2 for z)");
 
  return params;
}

Member Data Documentation

◆ _component

const unsigned int HorizonStabilizedFormISmallStrainMechanicsNOSPD::_component

protected

The index of displacement component.

Definition at line 39 of file HorizonStabilizedFormISmallStrainMechanicsNOSPD.h.

Referenced by computeLocalJacobian(), computeLocalOffDiagJacobian(), computeLocalResidual(), computeNonlocalJacobian(), and computePDNonlocalOffDiagJacobian().

◆ _current_unit_vec

RealGradient MechanicsBasePD::_current_unit_vec

protectedinherited

Unit vector of bond in current configuration.

Definition at line 79 of file MechanicsBasePD.h.

Referenced by GeneralizedPlaneStrainOffDiagOSPD::computeDispFullOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagOSPD::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 76 of file MechanicsBasePD.h.

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

◆ _ddgraddu

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_ddgraddu

protectedinherited

Definition at line 39 of file MechanicsBaseNOSPD.h.

Referenced by MechanicsFiniteStrainBaseNOSPD::computeDinvFTDU(), MechanicsFiniteStrainBaseNOSPD::computeDJDU(), MechanicsBaseNOSPD::computeDSDU(), and MechanicsFiniteStrainBaseNOSPD::computeDSDU().

◆ _ddgraddv

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_ddgraddv

protectedinherited

Definition at line 40 of file MechanicsBaseNOSPD.h.

Referenced by MechanicsFiniteStrainBaseNOSPD::computeDinvFTDU(), MechanicsFiniteStrainBaseNOSPD::computeDJDU(), MechanicsBaseNOSPD::computeDSDU(), and MechanicsFiniteStrainBaseNOSPD::computeDSDU().

◆ _ddgraddw

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_ddgraddw

protectedinherited

Definition at line 41 of file MechanicsBaseNOSPD.h.

Referenced by MechanicsFiniteStrainBaseNOSPD::computeDinvFTDU(), MechanicsFiniteStrainBaseNOSPD::computeDJDU(), MechanicsBaseNOSPD::computeDSDU(), and MechanicsFiniteStrainBaseNOSPD::computeDSDU().

◆ _deigenstrain_dT

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

protectedinherited

Definition at line 44 of file MechanicsBaseNOSPD.h.

Referenced by ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), computeLocalOffDiagJacobian(), WeakPlaneStressNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::computeTempOffDiagJacobianScalar(), and MechanicsBaseNOSPD::MechanicsBaseNOSPD().

◆ _dgrad

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_dgrad

protectedinherited

Definition at line 38 of file MechanicsBaseNOSPD.h.

◆ _disp_var

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

protectedinherited

displacement variables

Definition at line 51 of file MechanicsBasePD.h.

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

◆ _eigenstrain_names

const std::vector<MaterialPropertyName> MechanicsBaseNOSPD::_eigenstrain_names

protectedinherited

Definition at line 43 of file MechanicsBaseNOSPD.h.

Referenced by MechanicsBaseNOSPD::MechanicsBaseNOSPD().

◆ _ivardofs

std::vector<dof_id_type> MechanicsBasePD::_ivardofs

protectedinherited

Current variable dof numbers for nodes i and j.

Definition at line 70 of file MechanicsBasePD.h.

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

◆ _Jacobian_mult

const MaterialProperty<RankFourTensor>& MechanicsBaseNOSPD::_Jacobian_mult

protectedinherited

Definition at line 42 of file MechanicsBaseNOSPD.h.

◆ _multi

const MaterialProperty<Real>& MechanicsBaseNOSPD::_multi

protectedinherited

Material point based material properties.

Definition at line 35 of file MechanicsBaseNOSPD.h.

◆ _ndisp

unsigned int MechanicsBasePD::_ndisp

protectedinherited

number of displacement components

Definition at line 59 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 67 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 62 of file MechanicsBasePD.h.

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

◆ _out_of_plane_strain_var

MooseVariable* MechanicsBasePD::_out_of_plane_strain_var

protectedinherited

Definition at line 63 of file MechanicsBasePD.h.

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

◆ _shape2

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_shape2

protectedinherited

Definition at line 37 of file MechanicsBaseNOSPD.h.

◆ _stress

const MaterialProperty<RankTwoTensor>& MechanicsBaseNOSPD::_stress

protectedinherited

Definition at line 36 of file MechanicsBaseNOSPD.h.

◆ _temp_coupled

const bool MechanicsBasePD::_temp_coupled

protectedinherited

Temperature variable.

Definition at line 54 of file MechanicsBasePD.h.

◆ _temp_var

MooseVariable* MechanicsBasePD::_temp_var

protectedinherited

Definition at line 55 of file MechanicsBasePD.h.

◆ _weights

std::vector<Real> MechanicsBasePD::_weights

protectedinherited

weights used for the current element to obtain the nodal stress

Definition at line 73 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/HorizonStabilizedFormISmallStrainMechanicsNOSPD.h
peridynamics/src/kernels/HorizonStabilizedFormISmallStrainMechanicsNOSPD.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

◆ HorizonStabilizedFormISmallStrainMechanicsNOSPD()

Member Function Documentation

◆ computeDSDU()

◆ computeLocalJacobian()

◆ computeLocalOffDiagJacobian()

◆ computeLocalResidual()

◆ computeNonlocalJacobian()

◆ computeOffDiagJacobian()

◆ computePDNonlocalOffDiagJacobian()

◆ initialSetup()

◆ prepare()

◆ validParams()

Member Data Documentation

◆ _component

◆ _current_unit_vec

◆ _current_vec

◆ _ddgraddu

◆ _ddgraddv

◆ _ddgraddw

◆ _deigenstrain_dT

◆ _dgrad

◆ _disp_var

◆ _eigenstrain_names

◆ _ivardofs

◆ _Jacobian_mult

◆ _multi

◆ _ndisp

◆ _orientation

◆ _out_of_plane_strain_coupled

◆ _out_of_plane_strain_var

◆ _shape2

◆ _stress

◆ _temp_coupled

◆ _temp_var

◆ _weights