RBMPresetOldValuePD Class Reference

Class to apply preset BC of old variable solution based on the singularity of nodal shape tensor. More...

#include <RBMPresetOldValuePD.h>

Inheritance diagram for RBMPresetOldValuePD:

Public Member Functions
	RBMPresetOldValuePD (const InputParameters &parameters)
virtual Real	computeQpValue () override
virtual bool	shouldApply () override

Protected Attributes
PeridynamicsMesh &	_pdmesh
	Peridynamic mesh. More...
const VariableValue &	_u_old
	Value of the unknown variable this BC is acting on at last time step. More...
MooseVariable *	_bond_status_var
	AuxVariable for number of intact bonds associated with each material point. More...

Detailed Description

Class to apply preset BC of old variable solution based on the singularity of nodal shape tensor.

Used to fix nodes with rigid body motion.

Definition at line 24 of file RBMPresetOldValuePD.h.

Constructor & Destructor Documentation

RBMPresetOldValuePD()

RBMPresetOldValuePD::RBMPresetOldValuePD

(

const InputParameters &

parameters

)

Definition at line 28 of file RBMPresetOldValuePD.C.

  : DirichletBCBase(parameters),
    _pdmesh(dynamic_cast<PeridynamicsMesh &>(_mesh)),
    _u_old(_var.dofValuesOld()),
    _bond_status_var(&_subproblem.getStandardVariable(_tid, "bond_status"))
{
}

Member Function Documentation

computeQpValue()

Real RBMPresetOldValuePD::computeQpValue ( )

overridevirtual

Definition at line 37 of file RBMPresetOldValuePD.C.

{
  return _u_old[_qp];
}

shouldApply()

bool RBMPresetOldValuePD::shouldApply ( )

overridevirtual

Definition at line 43 of file RBMPresetOldValuePD.C.

{
  bool should_apply = true;
 
  // check whether the node shape tensor is singular
  unsigned int dim = _pdmesh.dimension();
  Real horiz_size = _pdmesh.getHorizon(_current_node->id());
  std::vector<dof_id_type> neighbors = _pdmesh.getNeighbors(_current_node->id());
  std::vector<dof_id_type> bonds = _pdmesh.getBonds(_current_node->id());
 
  RankTwoTensor shape_tensor;
  if (dim == 2)
    shape_tensor(2, 2) = 1.0;
 
  RealGradient ori_vec(dim);
  for (unsigned int i = 0; i < neighbors.size(); ++i)
    if (_bond_status_var->getElementalValue(_pdmesh.elemPtr(bonds[i])) > 0.5)
    {
      Real vol_i = _pdmesh.getPDNodeVolume(neighbors[i]);
      ori_vec = *(_pdmesh.nodePtr(neighbors[i])) - *_current_node;
      for (unsigned int j = 0; j < dim; ++j)
        for (unsigned int k = 0; k < dim; ++k)
          shape_tensor(j, k) += horiz_size / ori_vec.norm() * ori_vec(j) * ori_vec(k) * vol_i;
    }
 
  if (!MooseUtils::absoluteFuzzyEqual(shape_tensor.det(), 0.0))
    should_apply = false;
 
  return should_apply;
}

Member Data Documentation

_bond_status_var

MooseVariable* RBMPresetOldValuePD::_bond_status_var

protected

AuxVariable for number of intact bonds associated with each material point.

Definition at line 40 of file RBMPresetOldValuePD.h.

Referenced by shouldApply().

_pdmesh

PeridynamicsMesh& RBMPresetOldValuePD::_pdmesh

protected

Peridynamic mesh.

Definition at line 34 of file RBMPresetOldValuePD.h.

Referenced by shouldApply().

_u_old

const VariableValue& RBMPresetOldValuePD::_u_old

protected

Value of the unknown variable this BC is acting on at last time step.

Definition at line 37 of file RBMPresetOldValuePD.h.

Referenced by computeQpValue().

---

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

• peridynamics/include/bcs/RBMPresetOldValuePD.h
• peridynamics/src/bcs/RBMPresetOldValuePD.C