Postprocessor class to compute L2 norm of displacements difference between prediction and analytical solution for peridynamic model. More...

#include <NodalDisplacementDifferenceL2NormPD.h>

Inheritance diagram for NodalDisplacementDifferenceL2NormPD:

[legend]

Public Member Functions
	NodalDisplacementDifferenceL2NormPD (const InputParameters &parameters)

virtual Real	getValue () override

virtual void	initialize () override

virtual void	execute () override

virtual void	threadJoin (const UserObject &uo) override

Protected Member Functions
virtual Real	computeNodalValue () override
	Function to evaluate the given function at each material point. More...

Protected Attributes
unsigned int	_n_disps
	Displacement variables. More...

std::vector< MooseVariable * >	_disp_var

std::vector< const Function * >	_funcs
	Known analytic displacement functions. More...

Real	_integral_value
	nodal area integral result More...

PeridynamicsMesh &	_pdmesh
	Reference to peridynamic mesh object. More...

const unsigned int	_dim
	Mesh dimension. More...

Detailed Description

Postprocessor class to compute L2 norm of displacements difference between prediction and analytical solution for peridynamic model.

Definition at line 23 of file NodalDisplacementDifferenceL2NormPD.h.

Constructor & Destructor Documentation

◆ NodalDisplacementDifferenceL2NormPD()

NodalDisplacementDifferenceL2NormPD::NodalDisplacementDifferenceL2NormPD ( const InputParameters & parameters )

Definition at line 32 of file NodalDisplacementDifferenceL2NormPD.C.

   : NodalIntegralPostprocessorBasePD(parameters)
 {
   const std::vector<NonlinearVariableName> & nl_vnames(
       getParam<std::vector<NonlinearVariableName>>("displacements"));
  
   const std::vector<FunctionName> & func_names(
       getParam<std::vector<FunctionName>>("analytic_functions"));
  
   _n_disps = nl_vnames.size();
   if (_n_disps > _dim)
     mooseError("Number of displacements components should not greater than problem dimension!");
  
   if (_n_disps != func_names.size())
     mooseError("Number of analytic_functions components should be the same as the number of "
                "displacements components!");
  
   for (unsigned int i = 0; i < _n_disps; ++i)
   {
     _disp_var.push_back(&_subproblem.getStandardVariable(_tid, nl_vnames[i]));
     _funcs.push_back(&getFunctionByName(func_names[i]));
   }
 }

Member Function Documentation

◆ computeNodalValue()

Real NodalDisplacementDifferenceL2NormPD::computeNodalValue ( )

overrideprotectedvirtual

Function to evaluate the given function at each material point.

Returns: The computed nodal value

Implements NodalIntegralPostprocessorBasePD.

Definition at line 64 of file NodalDisplacementDifferenceL2NormPD.C.

 {
   Real diff = 0;
   for (unsigned int i = 0; i < _n_disps; ++i)
     diff += (_disp_var[i]->getNodalValue(*_current_node) - _funcs[i]->value(_t, *_current_node)) *
             (_disp_var[i]->getNodalValue(*_current_node) - _funcs[i]->value(_t, *_current_node));
  
   return diff;
 }

◆ execute()

void NodalIntegralPostprocessorBasePD::execute ( )

overridevirtualinherited

Definition at line 35 of file NodalIntegralPostprocessorBasePD.C.

 {
   _integral_value += computeNodalValue() * _pdmesh.getPDNodeVolume(_current_node->id());
 }

◆ getValue()

Real NodalDisplacementDifferenceL2NormPD::getValue ( )

overridevirtual

Reimplemented from NodalIntegralPostprocessorBasePD.

Definition at line 58 of file NodalDisplacementDifferenceL2NormPD.C.

 {
   return std::sqrt(NodalIntegralPostprocessorBasePD::getValue());
 }

◆ initialize()

void NodalIntegralPostprocessorBasePD::initialize ( )

overridevirtualinherited

Definition at line 29 of file NodalIntegralPostprocessorBasePD.C.

 {
   _integral_value = 0;
 }

◆ threadJoin()

void NodalIntegralPostprocessorBasePD::threadJoin ( const UserObject & uo )

overridevirtualinherited

Definition at line 49 of file NodalIntegralPostprocessorBasePD.C.

 {
   const NodalIntegralPostprocessorBasePD & pps =
       static_cast<const NodalIntegralPostprocessorBasePD &>(uo);
   _integral_value += pps._integral_value;
 }