Base kernel class for peridynamic solid mechanics models. More...

#include <MechanicsBasePD.h>

Inheritance diagram for MechanicsBasePD:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	MechanicsBasePD (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 Attributes
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 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

Base kernel class for peridynamic solid mechanics models.

Definition at line 18 of file MechanicsBasePD.h.

Constructor & Destructor Documentation

◆ MechanicsBasePD()

MechanicsBasePD::MechanicsBasePD ( const InputParameters & parameters )

Definition at line 27 of file MechanicsBasePD.C.

   : DerivativeMaterialInterface<PeridynamicsKernelBase>(parameters),
     _temp_coupled(isCoupled("temperature")),
     _temp_var(_temp_coupled ? getVar("temperature", 0) : nullptr),
     _ndisp(coupledComponents("displacements")),
     _out_of_plane_strain_coupled(isCoupled("out_of_plane_strain")),
     _out_of_plane_strain_var(_out_of_plane_strain_coupled ? getVar("out_of_plane_strain", 0)
                                                           : nullptr),
     _orientation(nullptr)
 {
   if (_ndisp != _dim)
     mooseError("Number of displacements should be consistent with mesh dimension!");
 
   for (unsigned int i = 0; i < _ndisp; ++i)
     _disp_var.push_back(getVar("displacements", i));
 }

Member Function Documentation

◆ computeLocalOffDiagJacobian()

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

inlinevirtual

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 computeOffDiagJacobian().

32 {};

◆ computeOffDiagJacobian()

void MechanicsBasePD::computeOffDiagJacobian ( unsigned int jvar )

overridevirtual

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

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

inlinevirtual

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 computeOffDiagJacobian().

40 {};

◆ initialSetup()

void MechanicsBasePD::initialSetup ( )

overridevirtual

Definition at line 45 of file MechanicsBasePD.C.

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

◆ prepare()

void MechanicsBasePD::prepare ( )

overridevirtual

Definition at line 51 of file MechanicsBasePD.C.

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

static

Definition at line 13 of file MechanicsBasePD.C.

Referenced by MechanicsBaseNOSPD::validParams(), MechanicsOSPD::validParams(), MechanicsBPD::validParams(), and GeneralizedPlaneStrainOffDiagOSPD::validParams().

 {
   InputParameters params = PeridynamicsKernelBase::validParams();
   params.addClassDescription(
       "Base class for calculating the residual and Jacobian for the peridynamic mechanic kernels");
 
   params.addCoupledVar("temperature", "Nonlinear variable name for the temperature");
   params.addCoupledVar("out_of_plane_strain",
                        "Nonlinear variable name for the out_of_plane strain for "
                        "plane stress analysis using NOSPD formulation");
 
   return params;
 }

Member Data Documentation

◆ _current_unit_vec

RealGradient MechanicsBasePD::_current_unit_vec

protected

Unit vector of bond in current configuration.

Definition at line 75 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 prepare().

◆ _current_vec

RealGradient MechanicsBasePD::_current_vec

protected

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 prepare().

◆ _disp_var

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

protected

displacement variables

Definition at line 47 of file MechanicsBasePD.h.

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

◆ _ivardofs

std::vector<dof_id_type> MechanicsBasePD::_ivardofs

protected

Current variable dof numbers for nodes i and j.

Definition at line 66 of file MechanicsBasePD.h.

Referenced by GeneralizedPlaneStrainOffDiagOSPD::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 prepare().

◆ _ndisp

unsigned int MechanicsBasePD::_ndisp

protected

number of displacement components

Definition at line 55 of file MechanicsBasePD.h.

Referenced by MechanicsBasePD().

◆ _orientation

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

protected

Vector of bond in current configuration.

Definition at line 63 of file MechanicsBasePD.h.

Referenced by initialSetup().

◆ _out_of_plane_strain_coupled

const bool MechanicsBasePD::_out_of_plane_strain_coupled

protected

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

Definition at line 58 of file MechanicsBasePD.h.

◆ _out_of_plane_strain_var

MooseVariable* MechanicsBasePD::_out_of_plane_strain_var

protected

Definition at line 59 of file MechanicsBasePD.h.

◆ _temp_coupled

const bool MechanicsBasePD::_temp_coupled

protected

Temperature variable.

Definition at line 50 of file MechanicsBasePD.h.

◆ _temp_var

MooseVariable* MechanicsBasePD::_temp_var

protected

Definition at line 51 of file MechanicsBasePD.h.

◆ _weights

std::vector<Real> MechanicsBasePD::_weights

protected

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 prepare().

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

peridynamics/include/kernels/MechanicsBasePD.h
peridynamics/src/kernels/MechanicsBasePD.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MechanicsBasePD()

Member Function Documentation

◆ computeLocalOffDiagJacobian()

◆ computeOffDiagJacobian()

◆ computePDNonlocalOffDiagJacobian()

◆ initialSetup()

◆ prepare()

◆ validParams()

Member Data Documentation

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