#include <NormalNodalMechanicalContact.h>

Inheritance diagram for NormalNodalMechanicalContact:

Public Member Functions
	NormalNodalMechanicalContact (const InputParameters &parameters)

void	computeJacobian () override

void	computeOffDiagJacobian (unsigned int jvar) override

Protected Member Functions
virtual Real	computeQpSlaveValue () override

virtual Real	computeQpResidual (Moose::ConstraintType type) override

virtual Real	computeQpJacobian (Moose::ConstraintJacobianType type) override

virtual Real	computeQpOffDiagJacobian (Moose::ConstraintJacobianType type, unsigned jvar) override

Protected Attributes
const Real &	_lambda

const unsigned	_lambda_id

const Real	_epsilon

const MooseEnum	_component

Detailed Description

Definition at line 21 of file NormalNodalMechanicalContact.h.

Constructor & Destructor Documentation

◆ NormalNodalMechanicalContact()

NormalNodalMechanicalContact::NormalNodalMechanicalContact ( const InputParameters & parameters )

Definition at line 32 of file NormalNodalMechanicalContact.C.

   : NodeFaceConstraint(parameters),
     _lambda(getVar("lambda", 0)->nodalValue()),
     _lambda_id(coupled("lambda")),
     _epsilon(std::numeric_limits<Real>::epsilon()),
     _component(getParam<MooseEnum>("component"))
 {
   _overwrite_slave_residual = false;
 }

Member Function Documentation

◆ computeJacobian()

void NormalNodalMechanicalContact::computeJacobian ( )

override

Definition at line 49 of file NormalNodalMechanicalContact.C.

50 {

51 }

◆ computeOffDiagJacobian()

void NormalNodalMechanicalContact::computeOffDiagJacobian ( unsigned int jvar )

override

Definition at line 54 of file NormalNodalMechanicalContact.C.

 {
   // Our residual only strongly depends on the lagrange multiplier (the normal vector does indeed
   // depend on displacements but it's complicated and shouldn't be too strong of a dependence)
   if (jvar != _lambda_id)
     return;
  
   MooseVariableFEBase & var = _sys.getVariable(0, jvar);
   _connected_dof_indices.clear();
   _connected_dof_indices.push_back(var.nodalDofIndex());
  
   _qp = 0;
  
   _Kee.resize(1, 1);
   _Kee(0, 0) = computeQpOffDiagJacobian(Moose::SlaveSlave, jvar);
  
   _Kne.resize(_test_master.size(), 1);
  
   for (_i = 0; _i < _test_master.size(); ++_i)
     _Kne(_i, 0) = computeQpOffDiagJacobian(Moose::MasterSlave, jvar);
 }

◆ computeQpJacobian()

Real NormalNodalMechanicalContact::computeQpJacobian ( Moose::ConstraintJacobianType type )

overrideprotectedvirtual

Definition at line 104 of file NormalNodalMechanicalContact.C.

 {
   return 0;
 }

◆ computeQpOffDiagJacobian()

Real NormalNodalMechanicalContact::computeQpOffDiagJacobian	(	Moose::ConstraintJacobianType	type,
		unsigned	jvar
	)

overrideprotectedvirtual

Definition at line 110 of file NormalNodalMechanicalContact.C.

 {
   std::map<dof_id_type, PenetrationInfo *>::iterator found =
       _penetration_locator._penetration_info.find(_current_node->id());
   if (found != _penetration_locator._penetration_info.end())
   {
     PenetrationInfo * pinfo = found->second;
  
     // Latter check here is actually redundant because we don't call into this function unless jvar
     // == _lambda_id
     if (pinfo && jvar == _lambda_id)
     {
       switch (type)
       {
         case Moose::SlaveSlave:
           return -pinfo->_normal(_component);
         case Moose::MasterSlave:
           return _test_master[_i][_qp] * pinfo->_normal(_component);
         default:
           return 0;
       }
     }
   }
  
   return 0;
 }

Referenced by computeOffDiagJacobian().

◆ computeQpResidual()

Real NormalNodalMechanicalContact::computeQpResidual ( Moose::ConstraintType type )

overrideprotectedvirtual

Definition at line 77 of file NormalNodalMechanicalContact.C.

 {
   std::map<dof_id_type, PenetrationInfo *>::iterator found =
       _penetration_locator._penetration_info.find(_current_node->id());
   if (found != _penetration_locator._penetration_info.end())
   {
     PenetrationInfo * pinfo = found->second;
     if (pinfo)
     {
       switch (type)
       {
         case Moose::ConstraintType::Slave:
           // This normal appears to point in the opposite direction of the slave surface so we need
           // a negative sign
           return _lambda * -pinfo->_normal(_component);
  
         case Moose::ConstraintType::Master:
           return _test_master[_i][_qp] * _lambda * pinfo->_normal(_component);
  
         default:
           return 0;
       }
     }
   }
   return 0;
 }