Kernel class for ordinary state based peridynamic solid mechanics models for small strain. More...

#include <MechanicsOSPD.h>

Inheritance diagram for MechanicsOSPD:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	MechanicsOSPD (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	computeNonlocalResidual () override

virtual void	computeLocalJacobian () override

virtual void	computeNonlocalJacobian () override

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

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

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

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_force
	Bond based material properties. More...

const MaterialProperty< Real > &	_bond_nonlocal_force

const MaterialProperty< Real > &	_bond_local_dfdU

const MaterialProperty< Real > &	_bond_nonlocal_dfdU

const MaterialProperty< Real > &	_bond_local_dfdT

const MaterialProperty< Real > &	_bond_nonlocal_dfdT


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 ordinary state based peridynamic solid mechanics models for small strain.

Definition at line 17 of file MechanicsOSPD.h.

Constructor & Destructor Documentation

◆ MechanicsOSPD()

MechanicsOSPD::MechanicsOSPD ( const InputParameters & parameters )

Definition at line 30 of file MechanicsOSPD.C.

   : MechanicsBasePD(parameters),
     _bond_local_force(getMaterialProperty<Real>("bond_local_force")),
     _bond_nonlocal_force(getMaterialProperty<Real>("bond_nonlocal_force")),
     _bond_local_dfdU(getMaterialProperty<Real>("bond_dfdU")),
     _bond_nonlocal_dfdU(getMaterialProperty<Real>("bond_nonlocal_dfdU")),
     _bond_local_dfdT(getMaterialProperty<Real>("bond_dfdT")),
     _bond_nonlocal_dfdT(getMaterialProperty<Real>("bond_nonlocal_dfdT")),
     _component(getParam<unsigned int>("component"))
 {
 }

Member Function Documentation

◆ computeLocalJacobian()

void MechanicsOSPD::computeLocalJacobian ( )

overrideprotectedvirtual

Definition at line 111 of file MechanicsOSPD.C.

 {
   const Real val =
       _current_unit_vec(_component) * _current_unit_vec(_component) * _bond_local_dfdU[0] +
       _bond_local_force[0] * (1.0 - _current_unit_vec(_component) * _current_unit_vec(_component)) /
           _current_vec.norm();
 
   for (unsigned int i = 0; i < _nnodes; ++i)
     for (unsigned int j = 0; j < _nnodes; ++j)
       _local_ke(i, j) += (i == j ? 1 : -1) * val * _bond_status;
 }

◆ computeLocalOffDiagJacobian()

void MechanicsOSPD::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 124 of file MechanicsOSPD.C.

 {
   if (_temp_coupled && jvar_num == _temp_var->number())
   {
     for (unsigned int i = 0; i < _nnodes; ++i)
       for (unsigned int j = 0; j < _nnodes; ++j)
         _local_ke(i, j) +=
             (i == 1 ? 1 : -1) * _current_unit_vec(_component) * _bond_local_dfdT[0] * _bond_status;
   }
   else
   {
     const Real val =
         _current_unit_vec(_component) * _current_unit_vec(coupled_component) * _bond_local_dfdU[0] -
         _bond_local_force[0] * _current_unit_vec(_component) *
             _current_unit_vec(coupled_component) / _current_vec.norm();
 
     for (unsigned int i = 0; i < _nnodes; ++i)
       for (unsigned int j = 0; j < _nnodes; ++j)
         _local_ke(i, j) += (i == j ? 1 : -1) * val * _bond_status;
   }
 }

◆ computeLocalResidual()

void MechanicsOSPD::computeLocalResidual ( )

overrideprotectedvirtual

Definition at line 43 of file MechanicsOSPD.C.

 {
   // H term
   _local_re(0) = -_bond_local_force[0] * _current_unit_vec(_component) * _bond_status;
   _local_re(1) = -_local_re(0);
 }

◆ computeNonlocalJacobian()

void MechanicsOSPD::computeNonlocalJacobian ( )

overrideprotectedvirtual

Definition at line 147 of file MechanicsOSPD.C.

 {
   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[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;
     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]);
 
         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));
 
         current_vec_nb /= current_vec_nb.norm();
 
         const Real val = (nd == 0 ? 1 : -1) * current_vec_nb(_component) *
                          _current_unit_vec(_component) * _bond_nonlocal_dfdU[nd];
 
         _local_ke.zero();
         for (unsigned int i = 0; i < _nnodes; ++i)
           for (unsigned int j = 0; j < _nnodes; ++j)
             _local_ke(i, j) +=
                 (i == j ? 1 : -1) * val / origin_vec_nb.norm() * vol_nb * _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);
 
           diag(1 - nd) = 0;
 
           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[nd] = _current_elem->node_ptr(nd)->dof_number(
                 _diag_save_in[i]->sys().number(), _diag_save_in[i]->number(), 0);
             diag_save_in_dofs[1 - nd] =
                 _pdmesh.nodePtr(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);
           }
         }
 
         const Real val2 = _bond_nonlocal_force[nd] *
                           (1.0 - current_vec_nb(_component) * current_vec_nb(_component)) /
                           current_vec_nb.norm();
 
         _local_ke.zero();
         for (unsigned int i = 0; i < _nnodes; ++i)
           for (unsigned int j = 0; j < _nnodes; ++j)
             _local_ke(i, j) +=
                 (i == j ? 1 : -1) * val2 / origin_vec_nb.norm() * vol_nb * _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[nd] = _current_elem->node_ptr(nd)->dof_number(
                 _diag_save_in[i]->sys().number(), _diag_save_in[i]->number(), 0);
             diag_save_in_dofs[1 - nd] =
                 _pdmesh.nodePtr(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);
           }
         }
       }
   }
 }

◆ computeNonlocalResidual()

void MechanicsOSPD::computeNonlocalResidual ( )

overrideprotectedvirtual

Definition at line 51 of file MechanicsOSPD.C.

 {
   // P and Q terms
   for (unsigned int nd = 0; nd < _nnodes; ++nd)
   {
     // calculation of residual 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;
     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]);
 
         origin_vec_nb =
             _pdmesh.getNodeCoord(neighbors[nb]) - _pdmesh.getNodeCoord(_current_elem->node_id(nd));
 
         for (unsigned int i = 0; i < _dim; ++i)
           current_vec_nb(i) = origin_vec_nb(i) +
                               _disp_var[i]->getNodalValue(*_pdmesh.nodePtr(neighbors[nb])) -
                               _disp_var[i]->getNodalValue(*_current_elem->node_ptr(nd));
 
         current_vec_nb /= current_vec_nb.norm();
 
         _local_re(0) = (nd == 0 ? -1 : 1) * _bond_nonlocal_force[nd] * vol_nb /
                        origin_vec_nb.norm() * current_vec_nb(_component) * _bond_status;
         _local_re(1) = -_local_re(0);
 
         // cache the residual contribution to node_i and its neighbor nb using their global dof
         // indices
         addResiduals(_assembly, _local_re, ivardofs, _var.scalingFactor());
 
         // save in the displacement residuals
         if (_has_save_in)
         {
           Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
           for (unsigned int i = 0; i < _save_in.size(); ++i)
           {
             std::vector<dof_id_type> save_in_dofs(2);
             save_in_dofs[nd] = _current_elem->node_ptr(nd)->dof_number(
                 _save_in[i]->sys().number(), _save_in[i]->number(), 0);
             save_in_dofs[1 - nd] =
                 _pdmesh.nodePtr(neighbors[nb])
                     ->dof_number(_save_in[i]->sys().number(), _save_in[i]->number(), 0);
 
             _save_in[i]->sys().solution().add_vector(_local_re, 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 MechanicsOSPD::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 248 of file MechanicsOSPD.C.

 {
   std::vector<dof_id_type> jvardofs_ij(_nnodes);
   for (unsigned int nd = 0; nd < _nnodes; ++nd)
     jvardofs_ij[nd] = _current_elem->node_ptr(nd)->dof_number(_sys.number(), jvar_num, 0);
 
   for (unsigned int nd = 0; nd < _nnodes; ++nd)
   {
     // calculation of jacobian contribution to current_node's neighbors
     std::vector<dof_id_type> ivardofs(_nnodes), jvardofs(_nnodes);
     ivardofs[nd] = _ivardofs[nd];
     jvardofs[nd] = jvardofs_ij[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;
     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);
         jvardofs[1 - nd] = _pdmesh.nodePtr(neighbors[nb])->dof_number(_sys.number(), jvar_num, 0);
         vol_nb = _pdmesh.getNodeVolume(neighbors[nb]);
 
         origin_vec_nb =
             _pdmesh.getNodeCoord(neighbors[nb]) - _pdmesh.getNodeCoord(_current_elem->node_id(nd));
 
         for (unsigned int i = 0; i < _dim; ++i)
           current_vec_nb(i) = origin_vec_nb(i) +
                               _disp_var[i]->getNodalValue(*_pdmesh.nodePtr(neighbors[nb])) -
                               _disp_var[i]->getNodalValue(*_current_elem->node_ptr(nd));
 
         current_vec_nb /= current_vec_nb.norm();
 
         _local_ke.zero();
         if (_temp_coupled && jvar_num == _temp_var->number())
         {
           const Real val =
               current_vec_nb(_component) * _bond_nonlocal_dfdT[nd] / origin_vec_nb.norm() * vol_nb;
 
           _local_ke(0, nd) += (nd == 0 ? -1 : 1) * val * _bond_status;
           _local_ke(1, nd) += (nd == 0 ? 1 : -1) * val * _bond_status;
         }
         else
         {
           const Real val = (nd == 0 ? 1 : -1) * current_vec_nb(_component) *
                            _current_unit_vec(coupled_component) * _bond_nonlocal_dfdU[nd] /
                            origin_vec_nb.norm() * vol_nb;
 
           for (unsigned int i = 0; i < _nnodes; ++i)
             for (unsigned int j = 0; j < _nnodes; ++j)
               _local_ke(i, j) += (i == j ? 1 : -1) * val * _bond_status;
         }
 
         addJacobian(_assembly, _local_ke, ivardofs, jvardofs_ij, _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.

Referenced by MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagOSPD::computeOffDiagJacobianScalar(), and GeneralizedPlaneStrainOffDiagNOSPD::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 MechanicsOSPD::validParams ( )

static

Definition at line 16 of file MechanicsOSPD.C.

 {
   InputParameters params = MechanicsBasePD::validParams();
   params.addClassDescription(
       "Class for calculating the residual and Jacobian for the ordinary state-based "
       "peridynamic mechanics formulation");
 
   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

◆ _bond_local_dfdT

const MaterialProperty<Real>& MechanicsOSPD::_bond_local_dfdT

protected

Definition at line 41 of file MechanicsOSPD.h.

Referenced by computeLocalOffDiagJacobian().

◆ _bond_local_dfdU

const MaterialProperty<Real>& MechanicsOSPD::_bond_local_dfdU

protected

Definition at line 39 of file MechanicsOSPD.h.

Referenced by computeLocalJacobian(), and computeLocalOffDiagJacobian().

◆ _bond_local_force

const MaterialProperty<Real>& MechanicsOSPD::_bond_local_force

protected

Bond based material properties.

Definition at line 37 of file MechanicsOSPD.h.

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

◆ _bond_nonlocal_dfdT

const MaterialProperty<Real>& MechanicsOSPD::_bond_nonlocal_dfdT

protected

Definition at line 42 of file MechanicsOSPD.h.

Referenced by computePDNonlocalOffDiagJacobian().

◆ _bond_nonlocal_dfdU

const MaterialProperty<Real>& MechanicsOSPD::_bond_nonlocal_dfdU

protected

Definition at line 40 of file MechanicsOSPD.h.

Referenced by computeNonlocalJacobian(), and computePDNonlocalOffDiagJacobian().

◆ _bond_nonlocal_force

const MaterialProperty<Real>& MechanicsOSPD::_bond_nonlocal_force

protected

Definition at line 38 of file MechanicsOSPD.h.

Referenced by computeNonlocalJacobian(), and computeNonlocalResidual().

◆ _component

const unsigned int MechanicsOSPD::_component

protected

The index of displacement component.

Definition at line 46 of file MechanicsOSPD.h.

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

◆ _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 GeneralizedPlaneStrainOffDiagOSPD::computeDispFullOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagOSPD::computeDispPartialOffDiagJacobianScalar(), MechanicsBPD::computeLocalJacobian(), computeLocalJacobian(), MechanicsBPD::computeLocalOffDiagJacobian(), computeLocalOffDiagJacobian(), computeLocalResidual(), MechanicsBPD::computeLocalResidual(), computeNonlocalJacobian(), 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(), computeLocalJacobian(), MechanicsBPD::computeLocalOffDiagJacobian(), 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 GeneralizedPlaneStrainOffDiagOSPD::computeDispFullOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeLocalResidual(), computeNonlocalJacobian(), computeNonlocalResidual(), MechanicsBasePD::computeOffDiagJacobian(), GeneralizedPlaneStrainOffDiagOSPD::computeOffDiagJacobianScalar(), GeneralizedPlaneStrainOffDiagNOSPD::computeOffDiagJacobianScalar(), computePDNonlocalOffDiagJacobian(), GeneralizedPlaneStrainOffDiagNOSPD::GeneralizedPlaneStrainOffDiagNOSPD(), GeneralizedPlaneStrainOffDiagOSPD::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 GeneralizedPlaneStrainOffDiagOSPD::computeDispFullOffDiagJacobianScalar(), ForceStabilizedSmallStrainMechanicsNOSPD::computeNonlocalJacobian(), computeNonlocalJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computeNonlocalJacobian(), computeNonlocalResidual(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeNonlocalResidual(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeNonlocalResidual(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBasePD::prepare().

◆ _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(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::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(), HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computeLocalOffDiagJacobian(), MechanicsBasePD::computeOffDiagJacobian(), HorizonStabilizedFormIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), HorizonStabilizedFormIISmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), and HorizonStabilizedFormIIFiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian().

◆ _temp_coupled

const bool MechanicsBasePD::_temp_coupled

protectedinherited

Temperature variable.

Definition at line 50 of file MechanicsBasePD.h.

◆ _temp_var

MooseVariable* MechanicsBasePD::_temp_var

protectedinherited

Definition at line 51 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 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/MechanicsOSPD.h
peridynamics/src/kernels/MechanicsOSPD.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MechanicsOSPD()

Member Function Documentation

◆ computeLocalJacobian()

◆ computeLocalOffDiagJacobian()

◆ computeLocalResidual()

◆ computeNonlocalJacobian()

◆ computeNonlocalResidual()

◆ computeOffDiagJacobian()

◆ computePDNonlocalOffDiagJacobian()

◆ initialSetup()

◆ prepare()

◆ validParams()

Member Data Documentation

◆ _bond_local_dfdT

◆ _bond_local_dfdU

◆ _bond_local_force

◆ _bond_nonlocal_dfdT

◆ _bond_nonlocal_dfdU

◆ _bond_nonlocal_force

◆ _component

◆ _current_unit_vec

◆ _current_vec

◆ _disp_var

◆ _ivardofs

◆ _ndisp

◆ _orientation

◆ _out_of_plane_strain_coupled

◆ _out_of_plane_strain_var

◆ _temp_coupled

◆ _temp_var

◆ _weights