PeridynamicsKernelBase Class Reference

Base kernel class for peridynamic models. More...

#include <PeridynamicsKernelBase.h>

Inheritance diagram for PeridynamicsKernelBase:

Public Member Functions
	PeridynamicsKernelBase (const InputParameters &parameters)
virtual void	computeResidual () override
virtual void	computeJacobian () override
virtual Real	computeQpResidual () override

Protected Member Functions
virtual void	computeLocalResidual ()=0
	Function to compute local contribution to the residual at the current nodes. More...
virtual void	computeNonlocalResidual ()
	Function to compute nonlocal contribution to the residual at the current nodes. More...
virtual void	computeLocalJacobian ()
	Function to compute local contribution to the diagonal Jacobian at the current nodes. More...
virtual void	prepare ()
	Function to precalculate data which will be used in the derived classes. More...

Protected Attributes
MooseVariable *	_bond_status_var
	Bond_status variable. More...
const bool	_use_full_jacobian
	Option to use full jacobian including nonlocal constribution or not. More...
RealGradient	_origin_vec_ij
	Vector for current bond under undefored configuration. More...
Real	_bond_status_ij
	Bond status of current bond/edge2. More...
PeridynamicsMesh &	_pdmesh
	Parameters for peridynamic mesh information. More...
const unsigned int	_dim
const unsigned int	_nnodes
std::vector< Real >	_vols_ij
std::vector< Real >	_dg_vol_frac_ij
std::vector< Real >	_horiz_rad

Detailed Description

Base kernel class for peridynamic models.

Definition at line 25 of file PeridynamicsKernelBase.h.

Constructor & Destructor Documentation

PeridynamicsKernelBase()

PeridynamicsKernelBase::PeridynamicsKernelBase

(

const InputParameters &

parameters

)

Definition at line 29 of file PeridynamicsKernelBase.C.

  : Kernel(parameters),
    _bond_status_var(&_subproblem.getStandardVariable(_tid, "bond_status")),
    _use_full_jacobian(getParam<bool>("full_jacobian")),
    _pdmesh(dynamic_cast<PeridynamicsMesh &>(_mesh)),
    _dim(_pdmesh.dimension()),
    _nnodes(2),
    _vols_ij(_nnodes),
    _dg_vol_frac_ij(_nnodes),
    _horiz_rad(_nnodes)
{
}

Member Function Documentation

computeJacobian()

void PeridynamicsKernelBase::computeJacobian ( )

overridevirtual

Definition at line 87 of file PeridynamicsKernelBase.C.

{
  prepare();
 
  DenseMatrix<Number> & ke = _assembly.jacobianBlock(_var.number(), _var.number());
  _local_ke.resize(ke.m(), ke.n());
  _local_ke.zero();
 
  computeLocalJacobian();
 
  ke += _local_ke;
 
  if (_has_diag_save_in)
  {
    unsigned int rows = ke.m();
    DenseVector<Number> 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)
      _diag_save_in[i]->sys().solution().add_vector(diag, _diag_save_in[i]->dofIndices());
  }
 
  _local_ke.zero();
 
  if (_use_full_jacobian)
    computeNonlocalJacobian();
}

computeLocalJacobian()

virtual void PeridynamicsKernelBase::computeLocalJacobian ( )

inlineprotectedvirtual

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

Reimplemented in HeatConductionBPD.

Definition at line 48 of file PeridynamicsKernelBase.h.

{};

Referenced by computeJacobian().

computeLocalResidual()

virtual void PeridynamicsKernelBase::computeLocalResidual ( )

protectedpure virtual

Function to compute local contribution to the residual at the current nodes.

Implemented in HeatConductionBPD, and HeatSourceBPD.

Referenced by computeResidual().

computeNonlocalResidual()

virtual void PeridynamicsKernelBase::computeNonlocalResidual ( )

inlineprotectedvirtual

Function to compute nonlocal contribution to the residual at the current nodes.

Definition at line 43 of file PeridynamicsKernelBase.h.

{};

Referenced by computeResidual().

computeQpResidual()

virtual Real PeridynamicsKernelBase::computeQpResidual ( )

inlineoverridevirtual

Definition at line 32 of file PeridynamicsKernelBase.h.

{ return 0.0; }

computeResidual()

void PeridynamicsKernelBase::computeResidual ( )

overridevirtual

Definition at line 60 of file PeridynamicsKernelBase.C.

{
  prepare();
 
  DenseVector<Number> & re = _assembly.residualBlock(_var.number());
  mooseAssert(re.size() == _nnodes,
              "PeridynamicsKernelBase is designed to only work with EDGE2 elements");
  _local_re.resize(re.size());
  _local_re.zero();
 
  computeLocalResidual();
 
  re += _local_re;
 
  if (_has_save_in)
  {
    Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
    for (unsigned int i = 0; i < _save_in.size(); ++i)
      _save_in[i]->sys().solution().add_vector(_local_re, _save_in[i]->dofIndices());
  }
 
  _local_re.zero();
 
  computeNonlocalResidual();
}

prepare()

void PeridynamicsKernelBase::prepare ( )

protectedvirtual

Function to precalculate data which will be used in the derived classes.

Definition at line 43 of file PeridynamicsKernelBase.C.

{
  for (unsigned int nd = 0; nd < _nnodes; ++nd)
  {
    _vols_ij[nd] = _pdmesh.getPDNodeVolume(_current_elem->node_id(nd));
    dof_id_type id_ji_in_ij =
        _pdmesh.getNeighborIndex(_current_elem->node_id(nd), _current_elem->node_id(1 - nd));
    _dg_vol_frac_ij[nd] = _pdmesh.getDefGradVolFraction(_current_elem->node_id(nd), id_ji_in_ij);
    _horiz_rad[nd] = _pdmesh.getHorizon(_current_elem->node_id(nd));
  }
 
  _origin_vec_ij = _pdmesh.getPDNodeCoord(_current_elem->node_id(1)) -
                   _pdmesh.getPDNodeCoord(_current_elem->node_id(0));
  _bond_status_ij = _bond_status_var->getElementalValue(_current_elem);
}

Referenced by computeJacobian(), computeResidual(), and MechanicsBasePD::prepare().

Member Data Documentation

_bond_status_ij

Real PeridynamicsKernelBase::_bond_status_ij

protected

Bond status of current bond/edge2.

Definition at line 74 of file PeridynamicsKernelBase.h.

Referenced by HeatConductionBPD::computeLocalJacobian(), HeatConductionBPD::computeLocalResidual(), and prepare().

_bond_status_var

MooseVariable* PeridynamicsKernelBase::_bond_status_var

protected

Bond_status variable.

Definition at line 56 of file PeridynamicsKernelBase.h.

Referenced by prepare().

_dg_vol_frac_ij

std::vector<Real> PeridynamicsKernelBase::_dg_vol_frac_ij

protected

Definition at line 66 of file PeridynamicsKernelBase.h.

Referenced by prepare().

_dim

const unsigned int PeridynamicsKernelBase::_dim

protected

Definition at line 63 of file PeridynamicsKernelBase.h.

_horiz_rad

std::vector<Real> PeridynamicsKernelBase::_horiz_rad

protected

Definition at line 67 of file PeridynamicsKernelBase.h.

Referenced by prepare().

_nnodes

const unsigned int PeridynamicsKernelBase::_nnodes

protected

Definition at line 64 of file PeridynamicsKernelBase.h.

Referenced by HeatSourceBPD::computeLocalResidual(), computeResidual(), and prepare().

_origin_vec_ij

RealGradient PeridynamicsKernelBase::_origin_vec_ij

protected

Vector for current bond under undefored configuration.

Definition at line 71 of file PeridynamicsKernelBase.h.

Referenced by prepare().

_pdmesh

PeridynamicsMesh& PeridynamicsKernelBase::_pdmesh

protected

Parameters for peridynamic mesh information.

Definition at line 62 of file PeridynamicsKernelBase.h.

Referenced by HeatSourceBPD::computeLocalResidual(), and prepare().

_use_full_jacobian

const bool PeridynamicsKernelBase::_use_full_jacobian

protected

Option to use full jacobian including nonlocal constribution or not.

Definition at line 59 of file PeridynamicsKernelBase.h.

Referenced by computeJacobian().

_vols_ij

std::vector<Real> PeridynamicsKernelBase::_vols_ij

protected

Definition at line 65 of file PeridynamicsKernelBase.h.

Referenced by HeatSourceBPD::computeLocalResidual(), and prepare().

---

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

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