NormalNodalLMMechanicalContact Class Reference

#include <NormalNodalLMMechanicalContact.h>

Inheritance diagram for NormalNodalLMMechanicalContact:

Public Member Functions
	NormalNodalLMMechanicalContact (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpSlaveValue () override
virtual void	computeResidual () override
virtual void	computeJacobian () override
virtual void	computeOffDiagJacobian (unsigned jvar) 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 unsigned	_disp_y_id
const unsigned	_disp_z_id
const Real	_c
const Real	_epsilon
const MooseEnum	_ncp_type

Detailed Description

Definition at line 21 of file NormalNodalLMMechanicalContact.h.

Constructor & Destructor Documentation

NormalNodalLMMechanicalContact()

NormalNodalLMMechanicalContact::NormalNodalLMMechanicalContact

(

const InputParameters &

parameters

)

Definition at line 44 of file NormalNodalLMMechanicalContact.C.

  : NodeFaceConstraint(parameters),
    _disp_y_id(coupled("disp_y")),
    _disp_z_id(coupled("disp_z")),
    _c(getParam<Real>("c")),
    _epsilon(std::numeric_limits<Real>::epsilon()),
    _ncp_type(getParam<MooseEnum>("ncp_function_type"))
 
{
  _overwrite_slave_residual = false;
}

Member Function Documentation

computeJacobian()

void NormalNodalLMMechanicalContact::computeJacobian ( )

overrideprotectedvirtual

Definition at line 73 of file NormalNodalLMMechanicalContact.C.

{
  _Kee.resize(1, 1);
  _connected_dof_indices.clear();
  _connected_dof_indices.push_back(_var.nodalDofIndex());
 
  _qp = 0;
 
  _Kee(0, 0) += computeQpJacobian(Moose::SlaveSlave);
}

Referenced by computeOffDiagJacobian().

computeOffDiagJacobian()

void NormalNodalLMMechanicalContact::computeOffDiagJacobian ( unsigned  jvar )

overrideprotectedvirtual

Definition at line 85 of file NormalNodalLMMechanicalContact.C.

{
  if (jvar == _var.number())
  {
    computeJacobian();
    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);
 
  DenseMatrix<Number> & Ken =
      _assembly.jacobianBlockNeighbor(Moose::ElementNeighbor, _var.number(), jvar);
 
  for (_j = 0; _j < _phi_master.size(); ++_j)
    Ken(0, _j) += computeQpOffDiagJacobian(Moose::SlaveMaster, jvar);
}

computeQpJacobian()

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

overrideprotectedvirtual

Definition at line 137 of file NormalNodalLMMechanicalContact.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 != NULL)
    {
      DualNumber<Real, Real> dual_u_slave(_u_slave[_qp]);
      dual_u_slave.derivatives() = 1.;
 
      auto a = -pinfo->_distance * _c;
      auto b = dual_u_slave;
 
      if (_ncp_type == "fb")
        return (a + b - std::sqrt(a * a + b * b + _epsilon)).derivatives();
      else
        return std::min(a, b).derivatives();
    }
  }
  return 0;
}

Referenced by computeJacobian().

computeQpOffDiagJacobian()

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

overrideprotectedvirtual

Definition at line 162 of file NormalNodalLMMechanicalContact.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 != NULL)
    {
      DualNumber<Real, Real> gap(-pinfo->_distance);
 
      unsigned comp;
      if (jvar == _master_var_num)
        comp = 0;
      else if (jvar == _disp_y_id)
        comp = 1;
      else if (jvar == _disp_z_id)
        comp = 2;
      else
        return 0;
 
      gap.derivatives() = pinfo->_normal(comp);
 
      switch (type)
      {
        case Moose::SlaveSlave:
          gap.derivatives() *= 1;
          break;
        case Moose::SlaveMaster:
          gap.derivatives() *= -_phi_master[_j][_qp];
          break;
        default:
          mooseError("LMs do not have a master contribution.");
      }
 
      auto a = gap * _c;
      auto b = _u_slave[_qp];
 
      if (_ncp_type == "fb")
        return (a + b - std::sqrt(a * a + b * b + _epsilon)).derivatives();
      else
        return std::min(a, b).derivatives();
    }
  }
  return 0.;
}

Referenced by computeOffDiagJacobian().

computeQpResidual()

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

overrideprotectedvirtual

Definition at line 109 of file NormalNodalLMMechanicalContact.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 != NULL)
    {
      Real a = -pinfo->_distance * _c;
      mooseAssert(
          _qp < _u_slave.size(),
          "qp is greater than the size of u_slave in NormalNodalLMMechanicalContact. Check and "
          "make sure that your Lagrange multiplier variable has the same order as the mesh");
      Real b = _u_slave[_qp];
 
      if (_ncp_type == "fb")
        return a + b - std::sqrt(a * a + b * b + _epsilon);
      else
        return std::min(a, b);
    }
  }
  return 0;
}

Referenced by computeResidual().

computeQpSlaveValue()

Real NormalNodalLMMechanicalContact::computeQpSlaveValue ( )

overrideprotectedvirtual

Definition at line 57 of file NormalNodalLMMechanicalContact.C.

{
  return _u_slave[_qp];
}

computeResidual()

void NormalNodalLMMechanicalContact::computeResidual ( )

overrideprotectedvirtual

Definition at line 63 of file NormalNodalLMMechanicalContact.C.

{
  DenseVector<Number> & re = _assembly.residualBlock(_var.number());
 
  _qp = 0;
 
  re(0) = computeQpResidual(Moose::Slave);
}

Member Data Documentation

_c

const Real NormalNodalLMMechanicalContact::_c

protected

Definition at line 39 of file NormalNodalLMMechanicalContact.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

_disp_y_id

const unsigned NormalNodalLMMechanicalContact::_disp_y_id

protected

Definition at line 37 of file NormalNodalLMMechanicalContact.h.

Referenced by computeQpOffDiagJacobian().

_disp_z_id

const unsigned NormalNodalLMMechanicalContact::_disp_z_id

protected

Definition at line 38 of file NormalNodalLMMechanicalContact.h.

Referenced by computeQpOffDiagJacobian().

_epsilon

const Real NormalNodalLMMechanicalContact::_epsilon

protected

Definition at line 40 of file NormalNodalLMMechanicalContact.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

_ncp_type

const MooseEnum NormalNodalLMMechanicalContact::_ncp_type

protected

Definition at line 42 of file NormalNodalLMMechanicalContact.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

---

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

• contact/include/constraints/NormalNodalLMMechanicalContact.h
• contact/src/constraints/NormalNodalLMMechanicalContact.C