UserObject class to compute the residual and diagonal Jacobian components for scalar out-of-plane strain variable of generalized plane strain model based on ordinary state-based peridynamic model. More...

#include <GeneralizedPlaneStrainUserObjectOSPD.h>

Inheritance diagram for GeneralizedPlaneStrainUserObjectOSPD:

[legend]

Public Member Functions
	GeneralizedPlaneStrainUserObjectOSPD (const InputParameters &parameters)

virtual void	execute () override

virtual void	initialize () override

virtual void	threadJoin (const UserObject &uo) override

virtual void	finalize () override

Real	returnResidual () const
	Function to return the computed residual. More...

Real	returnJacobian () const
	Function to return the computed diagonal Jacobian. More...

Protected Attributes
MooseVariable *	_out_of_plane_stress_var
	Variable for out-of-plane stress component. More...

const MooseEnum	_strain
	Option of strain formulation: SMALL or FINITE. More...

const MaterialProperty< RankFourTensor > &	_Cijkl
	Elasticity tensor. More...

Real	_residual
	Residual parameter. More...

Real	_jacobian
	Jacobian parameter. More...

MooseVariable *	_bond_status_var
	Bond status aux variable. More...

AuxiliarySystem &	_aux
	Reference to auxiliary system. More...

NumericVector< Number > &	_aux_sln
	Solution vector for aux variables. More...

PeridynamicsMesh &	_pdmesh
	Reference to Peridynamic mesh. More...


const Function &	_pressure
	Applied out-of-plane force parameters. More...

const Real	_factor

Detailed Description

UserObject class to compute the residual and diagonal Jacobian components for scalar out-of-plane strain variable of generalized plane strain model based on ordinary state-based peridynamic model.

Definition at line 24 of file GeneralizedPlaneStrainUserObjectOSPD.h.

Constructor & Destructor Documentation

◆ GeneralizedPlaneStrainUserObjectOSPD()

GeneralizedPlaneStrainUserObjectOSPD::GeneralizedPlaneStrainUserObjectOSPD ( const InputParameters & parameters )

Definition at line 30 of file GeneralizedPlaneStrainUserObjectOSPD.C.

   : GeneralizedPlaneStrainUserObjectBasePD(parameters),
     _out_of_plane_stress_var(getVar("out_of_plane_stress_variable", 0))
 {
 }

Member Function Documentation

◆ execute()

void GeneralizedPlaneStrainUserObjectOSPD::execute ( )

overridevirtual

Definition at line 38 of file GeneralizedPlaneStrainUserObjectOSPD.C.

 {
   // dof_id_type for node i and j
   dof_id_type node_i = _current_elem->node_id(0);
   dof_id_type node_j = _current_elem->node_id(1);
  
   // coordinates for node i and j
   Point coord_i = *_pdmesh.nodePtr(node_i);
   Point coord_j = *_pdmesh.nodePtr(node_j);
  
   // nodal area for node i and j
   Real nv_i = _pdmesh.getPDNodeVolume(node_i);
   Real nv_j = _pdmesh.getPDNodeVolume(node_j);
  
   // number of neighbors for node i and j, used to avoid repeated accounting nodal stress in
   // element-wise loop
   unsigned int nn_i = _pdmesh.getNeighbors(node_i).size();
   unsigned int nn_j = _pdmesh.getNeighbors(node_j).size();
  
   // residual
   _residual += (_out_of_plane_stress_var->getNodalValue(*_current_elem->node_ptr(0)) -
                 _pressure.value(_t, coord_i) * _factor) *
                nv_i / nn_i;
   _residual += (_out_of_plane_stress_var->getNodalValue(*_current_elem->node_ptr(1)) -
                 _pressure.value(_t, coord_j) * _factor) *
                nv_j / nn_j;
  
   // diagonal jacobian
   _jacobian += _Cijkl[0](2, 2, 2, 2) * nv_i / nn_i + _Cijkl[0](2, 2, 2, 2) * nv_j / nn_j;
 }

◆ finalize()

void GeneralizedPlaneStrainUserObjectBasePD::finalize ( )

overridevirtualinherited

Definition at line 60 of file GeneralizedPlaneStrainUserObjectBasePD.C.

 {
   gatherSum(_residual);
   gatherSum(_jacobian);
 }

◆ initialize()

void GeneralizedPlaneStrainUserObjectBasePD::initialize ( )

overridevirtualinherited

Definition at line 44 of file GeneralizedPlaneStrainUserObjectBasePD.C.

 {
   _residual = 0;
   _jacobian = 0;
 }

◆ returnJacobian()

Real GeneralizedPlaneStrainUserObjectBasePD::returnJacobian ( ) const

inherited

Function to return the computed diagonal Jacobian.

Returns: The computed Jacobian

Definition at line 73 of file GeneralizedPlaneStrainUserObjectBasePD.C.

 {
   return _jacobian;
 }

Referenced by GeneralizedPlaneStrainPD::computeJacobian().

◆ returnResidual()

Real GeneralizedPlaneStrainUserObjectBasePD::returnResidual ( ) const

inherited

Function to return the computed residual.

Returns: The computed residual

Definition at line 67 of file GeneralizedPlaneStrainUserObjectBasePD.C.

 {
   return _residual;
 }

Referenced by GeneralizedPlaneStrainPD::computeResidual().

◆ threadJoin()

void GeneralizedPlaneStrainUserObjectBasePD::threadJoin ( const UserObject & uo )

overridevirtualinherited

Definition at line 51 of file GeneralizedPlaneStrainUserObjectBasePD.C.

 {
   const GeneralizedPlaneStrainUserObjectBasePD & gpsuo =
       static_cast<const GeneralizedPlaneStrainUserObjectBasePD &>(uo);
   _residual += gpsuo._residual;
   _jacobian += gpsuo._jacobian;
 }