TangentialNodalLMMechanicalContact Class Reference

#include <TangentialNodalLMMechanicalContact.h>

Inheritance diagram for TangentialNodalLMMechanicalContact:

Public Member Functions
	TangentialNodalLMMechanicalContact (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 Real &	_contact_pressure
const unsigned	_contact_pressure_id
const Real &	_disp_x_dot
const Real &	_disp_y_dot
const unsigned	_disp_y_id
const VariableValue &	_du_dot_du
const Real	_mu
const Real	_epsilon
const MooseEnum	_ncp_type

Detailed Description

Definition at line 21 of file TangentialNodalLMMechanicalContact.h.

Constructor & Destructor Documentation

TangentialNodalLMMechanicalContact()

TangentialNodalLMMechanicalContact::TangentialNodalLMMechanicalContact

(

const InputParameters &

parameters

)

Definition at line 52 of file TangentialNodalLMMechanicalContact.C.

  : NodeFaceConstraint(parameters),
    _contact_pressure(getVar("contact_pressure", 0)->nodalValue()),
    _contact_pressure_id(coupled("contact_pressure")),
    _disp_x_dot(_master_var.nodalValueDot()),
    _disp_y_dot(getVar("disp_y", 0)->nodalValueDot()),
    _disp_y_id(coupled("disp_y")),
    _du_dot_du(_master_var.dofValuesDuDotDu()),
 
    _mu(getParam<Real>("mu")),
    _epsilon(std::numeric_limits<Real>::epsilon()),
    _ncp_type(getParam<MooseEnum>("ncp_function_type"))
{
  _overwrite_slave_residual = false;
}

Member Function Documentation

computeJacobian()

void TangentialNodalLMMechanicalContact::computeJacobian ( )

overrideprotectedvirtual

Definition at line 86 of file TangentialNodalLMMechanicalContact.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 TangentialNodalLMMechanicalContact::computeOffDiagJacobian ( unsigned  jvar )

overrideprotectedvirtual

Definition at line 98 of file TangentialNodalLMMechanicalContact.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);
}

computeQpJacobian()

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

overrideprotectedvirtual

Definition at line 154 of file TangentialNodalLMMechanicalContact.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)
    {
      if (_contact_pressure < _epsilon)
        return 1.;
      else
      {
        RealVectorValue tangent_vec(pinfo->_normal(1), -pinfo->_normal(0), 0);
        Real v_dot_tan = RealVectorValue(_disp_x_dot, _disp_y_dot, 0) * tangent_vec;
 
        // NCP part 1: requirement that either there is no slip **or** slip velocity and
        // frictional force exerted **by** the slave side are in the same direction
        Real a;
        if (v_dot_tan * _u_slave[_qp] < 0)
          a = -std::abs(v_dot_tan);
        else
          a = std::abs(v_dot_tan);
 
        // NCP part 2: require that the frictional force can never exceed the frictional
        // coefficient times the normal force.
        DualNumber<Real, Real> dual_u_slave(_u_slave[_qp]);
        dual_u_slave.derivatives() = 1;
 
        auto b = _mu * _contact_pressure - std::abs(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 TangentialNodalLMMechanicalContact::computeQpOffDiagJacobian ( Moose::ConstraintJacobianType  type, unsigned  jvar  )

overrideprotectedvirtual

Definition at line 196 of file TangentialNodalLMMechanicalContact.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)
    {
      if (_contact_pressure < _epsilon)
        return 0.;
 
      // Our local dual number is going to depend on only three degrees of freedom: the slave nodal
      // dofs for disp_x (index 0), disp_y (index 1), and the contact pressure (index 2). The latter
      // of course exists only on the slave side
      typedef DualNumber<Real, DNDerivativeSize<3>> LocalDN;
 
      RealVectorValue tangent_vec(pinfo->_normal(1), -pinfo->_normal(0), 0);
 
      LocalDN dual_disp_x_dot(_disp_x_dot);
      // We index the zeroth entry of the _du_dot_du member variable because there is only one
      // degree of freedom on the node
      Moose::derivInsert(dual_disp_x_dot.derivatives(), 0, _du_dot_du[0]);
 
      LocalDN dual_disp_y_dot(_disp_y_dot);
      Moose::derivInsert(dual_disp_y_dot.derivatives(), 1, _du_dot_du[0]);
 
      LocalDN dual_contact_pressure(_contact_pressure);
      Moose::derivInsert(dual_contact_pressure.derivatives(), 2, 1);
 
      auto v_dot_tan = VectorValue<LocalDN>(dual_disp_x_dot, dual_disp_y_dot, 0) * tangent_vec;
 
      // NCP part 1: requirement that either there is no slip **or** slip velocity and
      // frictional force exerted **by** the slave side are in the same direction
      LocalDN a;
      if (v_dot_tan * _u_slave[_qp] < 0)
        a = -std::abs(v_dot_tan);
      else
        a = std::abs(v_dot_tan);
 
      // NCP part 2: require that the frictional force can never exceed the frictional
      // coefficient times the normal force.
      auto b = _mu * dual_contact_pressure - std::abs(_u_slave[_qp]);
 
      LocalDN ncp_value;
      if (_ncp_type == "fb")
        ncp_value = a + b - std::sqrt(a * a + b * b + _epsilon);
      else
        ncp_value = std::min(a, b);
 
      if (jvar == _contact_pressure_id)
        return ncp_value.derivatives()[2];
      else if (jvar == _master_var.number())
        return ncp_value.derivatives()[0];
      else if (jvar == _disp_y_id)
        return ncp_value.derivatives()[1];
    }
  }
  return 0.;
}

Referenced by computeOffDiagJacobian().

computeQpResidual()

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

overrideprotectedvirtual

Definition at line 116 of file TangentialNodalLMMechanicalContact.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)
    {
      if (_contact_pressure < _epsilon)
        return _u_slave[_qp];
      else
      {
        RealVectorValue tangent_vec(pinfo->_normal(1), -pinfo->_normal(0), 0);
        Real v_dot_tan = RealVectorValue(_disp_x_dot, _disp_y_dot, 0) * tangent_vec;
 
        // NCP part 1: requirement that either there is no slip **or** slip velocity and
        // frictional force exerted **by** the slave side are in the same direction
        Real a;
        if (v_dot_tan * _u_slave[_qp] < 0)
          a = -std::abs(v_dot_tan);
        else
          a = std::abs(v_dot_tan);
 
        // NCP part 2: require that the frictional force can never exceed the frictional
        // coefficient times the normal force.
        auto b = _mu * _contact_pressure - std::abs(_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 TangentialNodalLMMechanicalContact::computeQpSlaveValue ( )

overrideprotectedvirtual

Definition at line 70 of file TangentialNodalLMMechanicalContact.C.

{
  return _u_slave[_qp];
}

computeResidual()

void TangentialNodalLMMechanicalContact::computeResidual ( )

overrideprotectedvirtual

Definition at line 76 of file TangentialNodalLMMechanicalContact.C.

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

Member Data Documentation

_contact_pressure

const Real& TangentialNodalLMMechanicalContact::_contact_pressure

protected

Definition at line 37 of file TangentialNodalLMMechanicalContact.h.

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

_contact_pressure_id

const unsigned TangentialNodalLMMechanicalContact::_contact_pressure_id

protected

Definition at line 38 of file TangentialNodalLMMechanicalContact.h.

Referenced by computeQpOffDiagJacobian().

_disp_x_dot

const Real& TangentialNodalLMMechanicalContact::_disp_x_dot

protected

Definition at line 40 of file TangentialNodalLMMechanicalContact.h.

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

_disp_y_dot

const Real& TangentialNodalLMMechanicalContact::_disp_y_dot

protected

Definition at line 41 of file TangentialNodalLMMechanicalContact.h.

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

_disp_y_id

const unsigned TangentialNodalLMMechanicalContact::_disp_y_id

protected

Definition at line 43 of file TangentialNodalLMMechanicalContact.h.

Referenced by computeQpOffDiagJacobian().

_du_dot_du

const VariableValue& TangentialNodalLMMechanicalContact::_du_dot_du

protected

Definition at line 45 of file TangentialNodalLMMechanicalContact.h.

Referenced by computeQpOffDiagJacobian().

_epsilon

const Real TangentialNodalLMMechanicalContact::_epsilon

protected

Definition at line 48 of file TangentialNodalLMMechanicalContact.h.

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

_mu

const Real TangentialNodalLMMechanicalContact::_mu

protected

Definition at line 47 of file TangentialNodalLMMechanicalContact.h.

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

_ncp_type

const MooseEnum TangentialNodalLMMechanicalContact::_ncp_type

protected

Definition at line 50 of file TangentialNodalLMMechanicalContact.h.

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

---

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

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