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

#include <GeneralizedPlaneStrainUserObjectNOSPD.h>

Inheritance diagram for GeneralizedPlaneStrainUserObjectNOSPD:

[legend]

Public Member Functions
	GeneralizedPlaneStrainUserObjectNOSPD (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
const MaterialProperty< RankTwoTensor > &	_stress
	Materials property stress. 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 self-stablized non-ordinary state-based peridynamic model.

Definition at line 24 of file GeneralizedPlaneStrainUserObjectNOSPD.h.

Constructor & Destructor Documentation

◆ GeneralizedPlaneStrainUserObjectNOSPD()

GeneralizedPlaneStrainUserObjectNOSPD::GeneralizedPlaneStrainUserObjectNOSPD ( const InputParameters & parameters )

Definition at line 28 of file GeneralizedPlaneStrainUserObjectNOSPD.C.

   : GeneralizedPlaneStrainUserObjectBasePD(parameters),
     _stress(getMaterialProperty<RankTwoTensor>("stress"))
 {
 }

Member Function Documentation

◆ execute()

void GeneralizedPlaneStrainUserObjectNOSPD::execute ( )

overridevirtual

Definition at line 36 of file GeneralizedPlaneStrainUserObjectNOSPD.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);
  
   // sum of volumes of material points used in bond-associated deformation gradient calculation
   dof_id_type id_j_in_i = _pdmesh.getNeighborIndex(node_i, node_j);
   dof_id_type id_i_in_j = _pdmesh.getNeighborIndex(node_j, node_i);
  
   Real dg_vol_frac_i = _pdmesh.getDefGradVolFraction(node_i, id_j_in_i);
   Real dg_vol_frac_j = _pdmesh.getDefGradVolFraction(node_j, id_i_in_j);
  
   // residual
   _residual += (_stress[0](2, 2) - _pressure.value(_t, coord_i) * _factor) * nv_i * dg_vol_frac_i;
   _residual += (_stress[1](2, 2) - _pressure.value(_t, coord_j) * _factor) * nv_j * dg_vol_frac_j;
  
   // diagonal jacobian
   _jacobian +=
       _Cijkl[0](2, 2, 2, 2) * nv_i * dg_vol_frac_i + _Cijkl[1](2, 2, 2, 2) * nv_j * dg_vol_frac_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;
 }