Kernel class for bond based peridynamic solid mechanics models. More...

#include <MechanicsBPD.h>

Inheritance diagram for MechanicsBPD:

Public Member Functions
	MechanicsBPD (const InputParameters &parameters)

virtual void	computeOffDiagJacobian (MooseVariableFEBase &jvar) override

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

Protected Member Functions
virtual void	computeLocalResidual () override

virtual void	computeLocalJacobian () override

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

Protected Attributes
const unsigned int	_component
	The index of displcement 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_ij
	Current variable dof numbers for nodes i and j. More...

RealGradient	_cur_ori_ij
	Vector of bond in current configuration. More...

Real	_cur_len_ij
	Current bond length. More...


const MaterialProperty< Real > &	_bond_force_ij
	Bond based material properties. More...

const MaterialProperty< Real > &	_bond_dfdU_ij

const MaterialProperty< Real > &	_bond_dfdT_ij


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 bond based peridynamic solid mechanics models.

Definition at line 22 of file MechanicsBPD.h.

Constructor & Destructor Documentation

◆ MechanicsBPD()

MechanicsBPD::MechanicsBPD ( const InputParameters & parameters )

Definition at line 29 of file MechanicsBPD.C.

   : MechanicsBasePD(parameters),
     _bond_force_ij(getMaterialProperty<Real>("bond_force_ij")),
     _bond_dfdU_ij(getMaterialProperty<Real>("bond_dfdU_ij")),
     _bond_dfdT_ij(getMaterialProperty<Real>("bond_dfdT_ij")),
     _component(getParam<unsigned int>("component"))
 {
 }

Member Function Documentation

◆ computeLocalJacobian()

void MechanicsBPD::computeLocalJacobian ( )

overrideprotectedvirtual

Definition at line 46 of file MechanicsBPD.C.

 {
   Real diag =
       _cur_ori_ij(_component) * _cur_ori_ij(_component) * _bond_dfdU_ij[0] +
       _bond_force_ij[0] * (1.0 - _cur_ori_ij(_component) * _cur_ori_ij(_component)) / _cur_len_ij;
  
   for (_i = 0; _i < _test.size(); _i++)
     for (_j = 0; _j < _phi.size(); _j++)
       _local_ke(_i, _j) += (_i == _j ? 1 : -1) * diag * _bond_status_ij;
 }

◆ computeLocalOffDiagJacobian()

void MechanicsBPD::computeLocalOffDiagJacobian ( 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 58 of file MechanicsBPD.C.

 {
   if (coupled_component == 3)
   {
     for (_i = 0; _i < _test.size(); _i++)
       for (_j = 0; _j < _phi.size(); _j++)
         _local_ke(_i, _j) +=
             (_i == 1 ? 1 : -1) * _cur_ori_ij(_component) * _bond_dfdT_ij[_j] * _bond_status_ij;
   }
   else
   {
     Real val =
         _cur_ori_ij(_component) * _cur_ori_ij(coupled_component) * _bond_dfdU_ij[0] -
         _bond_force_ij[0] * _cur_ori_ij(_component) * _cur_ori_ij(coupled_component) / _cur_len_ij;
     for (_i = 0; _i < _test.size(); _i++)
       for (_j = 0; _j < _phi.size(); _j++)
         _local_ke(_i, _j) += (_i == _j ? 1 : -1) * val * _bond_status_ij;
   }
 }

◆ computeLocalResidual()

void MechanicsBPD::computeLocalResidual ( )

overrideprotectedvirtual

Definition at line 39 of file MechanicsBPD.C.

 {
   _local_re(0) = -_bond_force_ij[0] * _cur_ori_ij(_component) * _bond_status_ij;
   _local_re(1) = -_local_re(0);
 }

◆ computeOffDiagJacobian()

void MechanicsBasePD::computeOffDiagJacobian ( MooseVariableFEBase & jvar )

overridevirtualinherited

Definition at line 71 of file MechanicsBasePD.C.

 {
   prepare();
  
   if (jvar.number() == _var.number())
     computeJacobian();
   else
   {
     unsigned int coupled_component = 0;
     bool active = false;
  
     for (unsigned int i = 0; i < _dim; ++i)
       if (jvar.number() == _disp_var[i]->number())
       {
         coupled_component = i;
         active = true;
       }
  
     if (_temp_coupled && jvar.number() == _temp_var->number())
     {
       coupled_component = 3;
       active = true;
     }
  
     if (_out_of_plane_strain_coupled && jvar.number() == _out_of_plane_strain_var->number())
     {
       coupled_component = 4;
       active = true;
     }
  
     if (active)
     {
       DenseMatrix<Number> & ke = _assembly.jacobianBlock(_var.number(), jvar.number());
       _local_ke.resize(ke.m(), ke.n());
       _local_ke.zero();
  
       computeLocalOffDiagJacobian(coupled_component);
  
       ke += _local_ke;
  
       if (_use_full_jacobian)
         computePDNonlocalOffDiagJacobian(jvar.number(), coupled_component);
     }
   }
 }

◆ 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 MechanicsOSPD, FiniteStrainMechanicsNOSPD, ForceStabilizedSmallStrainMechanicsNOSPD, SmallStrainMechanicsNOSPD, and WeakPlaneStressNOSPD.

Definition at line 42 of file MechanicsBasePD.h.

43 {};

Referenced by MechanicsBasePD::computeOffDiagJacobian().

◆ initialSetup()

void MechanicsBasePD::initialSetup ( )

overridevirtualinherited

Definition at line 47 of file MechanicsBasePD.C.

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

◆ prepare()

void MechanicsBasePD::prepare ( )

overridevirtualinherited

Definition at line 53 of file MechanicsBasePD.C.

 {
   PeridynamicsKernelBase::prepare();
  
   _ivardofs_ij.resize(_nnodes);
  
   for (unsigned int i = 0; i < _nnodes; ++i)
     _ivardofs_ij[i] = _current_elem->node_ptr(i)->dof_number(_sys.number(), _var.number(), 0);
  
   for (unsigned int i = 0; i < _dim; ++i)
     _cur_ori_ij(i) = _origin_vec_ij(i) + _disp_var[i]->getNodalValue(*_current_elem->node_ptr(1)) -
                      _disp_var[i]->getNodalValue(*_current_elem->node_ptr(0));
  
   _cur_len_ij = _cur_ori_ij.norm();
   _cur_ori_ij /= _cur_len_ij;
 }

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

Member Data Documentation

◆ _bond_dfdT_ij

const MaterialProperty<Real>& MechanicsBPD::_bond_dfdT_ij

protected

Definition at line 35 of file MechanicsBPD.h.

Referenced by computeLocalOffDiagJacobian().

◆ _bond_dfdU_ij

const MaterialProperty<Real>& MechanicsBPD::_bond_dfdU_ij

protected

Definition at line 34 of file MechanicsBPD.h.

Referenced by computeLocalJacobian(), and computeLocalOffDiagJacobian().

◆ _bond_force_ij

const MaterialProperty<Real>& MechanicsBPD::_bond_force_ij

protected

Bond based material properties.

Definition at line 33 of file MechanicsBPD.h.

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

◆ _component

const unsigned int MechanicsBPD::_component

protected

The index of displcement component.

Definition at line 39 of file MechanicsBPD.h.

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

◆ _cur_len_ij

Real MechanicsBasePD::_cur_len_ij

protectedinherited

Current bond length.

Definition at line 75 of file MechanicsBasePD.h.

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

◆ _cur_ori_ij

RealGradient MechanicsBasePD::_cur_ori_ij

protectedinherited

Vector of bond in current configuration.

Definition at line 72 of file MechanicsBasePD.h.

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

◆ _disp_var

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

protectedinherited

displacement variables

Definition at line 50 of file MechanicsBasePD.h.

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

◆ _ivardofs_ij

std::vector<dof_id_type> MechanicsBasePD::_ivardofs_ij

protectedinherited

Current variable dof numbers for nodes i and j.

Definition at line 69 of file MechanicsBasePD.h.

Referenced by GeneralizedPlaneStrainOffDiagOSPD::computeDispFullOffDiagJacobianScalar(), FiniteStrainMechanicsNOSPD::computeNonlocalJacobian(), SmallStrainMechanicsNOSPD::computeNonlocalJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computeNonlocalJacobian(), MechanicsOSPD::computeNonlocalJacobian(), MechanicsOSPD::computeNonlocalResidual(), WeakPlaneStressNOSPD::computePDNonlocalOffDiagJacobian(), SmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), FiniteStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), ForceStabilizedSmallStrainMechanicsNOSPD::computePDNonlocalOffDiagJacobian(), MechanicsOSPD::computePDNonlocalOffDiagJacobian(), and MechanicsBasePD::prepare().

◆ _ndisp

unsigned int MechanicsBasePD::_ndisp

protectedinherited

number of displacement components

Definition at line 58 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 66 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 61 of file MechanicsBasePD.h.

Referenced by MechanicsBasePD::computeOffDiagJacobian().

◆ _out_of_plane_strain_var

MooseVariable* MechanicsBasePD::_out_of_plane_strain_var

protectedinherited

Definition at line 62 of file MechanicsBasePD.h.

Referenced by MechanicsBasePD::computeOffDiagJacobian().

◆ _temp_coupled

const bool MechanicsBasePD::_temp_coupled

protectedinherited

Temperature variable.

Definition at line 53 of file MechanicsBasePD.h.

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

◆ _temp_var

MooseVariable* MechanicsBasePD::_temp_var

protectedinherited

Definition at line 54 of file MechanicsBasePD.h.

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

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

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

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MechanicsBPD()

Member Function Documentation

◆ computeLocalJacobian()

◆ computeLocalOffDiagJacobian()

◆ computeLocalResidual()

◆ computeOffDiagJacobian()

◆ computePDNonlocalOffDiagJacobian()

◆ initialSetup()

◆ prepare()

Member Data Documentation

◆ _bond_dfdT_ij

◆ _bond_dfdU_ij

◆ _bond_force_ij

◆ _component

◆ _cur_len_ij

◆ _cur_ori_ij

◆ _disp_var

◆ _ivardofs_ij

◆ _ndisp

◆ _orientation

◆ _out_of_plane_strain_coupled

◆ _out_of_plane_strain_var

◆ _temp_coupled

◆ _temp_var