GluedContactConstraint Class Reference

A GluedContactConstraint forces the value of a variable to be the same on both sides of an interface. More...

#include <GluedContactConstraint.h>

Inheritance diagram for GluedContactConstraint:

Public Member Functions
	GluedContactConstraint (const InputParameters &parameters)
virtual	~GluedContactConstraint ()
virtual void	timestepSetup ()
virtual void	jacobianSetup ()
virtual void	updateContactSet (bool beginning_of_step=false)
virtual Real	computeQpSlaveValue ()
virtual Real	computeQpResidual (Moose::ConstraintType type)
virtual Real	computeQpJacobian (Moose::ConstraintJacobianType type)
virtual Real	computeQpOffDiagJacobian (Moose::ConstraintJacobianType type, unsigned int jvar)
	Compute off-diagonal Jacobian entries. More...
bool	shouldApply ()
virtual void	getConnectedDofIndices ()
	Gets the indices for all dofs conected to the constraint Get indices for all the slave Jacobian columns, not only those based on the mesh connectivity. More...

Protected Attributes
const unsigned int	_component
const ContactModel	_model
const ContactFormulation	_formulation
const Real	_penalty
const Real	_friction_coefficient
const Real	_tension_release
bool	_updateContactSet
NumericVector< Number > &	_residual_copy
unsigned int	_x_var
unsigned int	_y_var
unsigned int	_z_var
std::vector< unsigned int >	_vars
MooseVariable *	_nodal_area_var
SystemBase &	_aux_system
const NumericVector< Number > *	_aux_solution

Detailed Description

A GluedContactConstraint forces the value of a variable to be the same on both sides of an interface.

Definition at line 27 of file GluedContactConstraint.h.

Constructor & Destructor Documentation

GluedContactConstraint()

GluedContactConstraint::GluedContactConstraint

(

const InputParameters &

parameters

)

Definition at line 62 of file GluedContactConstraint.C.

  : SparsityBasedContactConstraint(parameters),
    _component(getParam<unsigned int>("component")),
    _model(getParam<MooseEnum>("model").getEnum<ContactModel>()),
    _formulation(getParam<MooseEnum>("formulation").getEnum<ContactFormulation>()),
    _penalty(getParam<Real>("penalty")),
    _friction_coefficient(getParam<Real>("friction_coefficient")),
    _tension_release(getParam<Real>("tension_release")),
    _updateContactSet(true),
    _residual_copy(_sys.residualGhosted()),
    _vars(3, libMesh::invalid_uint),
    _nodal_area_var(getVar("nodal_area", 0)),
    _aux_system(_nodal_area_var->sys()),
    _aux_solution(_aux_system.currentSolution())
{
  if (parameters.isParamValid("tangential_tolerance"))
    _penetration_locator.setTangentialTolerance(getParam<Real>("tangential_tolerance"));
 
  if (parameters.isParamValid("normal_smoothing_distance"))
    _penetration_locator.setNormalSmoothingDistance(getParam<Real>("normal_smoothing_distance"));
 
  if (parameters.isParamValid("normal_smoothing_method"))
    _penetration_locator.setNormalSmoothingMethod(
        parameters.get<MooseEnum>("normal_smoothing_method"));
 
  if (isParamValid("displacements"))
  {
    // modern parameter scheme for displacements
    for (unsigned int i = 0; i < coupledComponents("displacements"); ++i)
      _vars[i] = coupled("displacements", i);
  }
  else
  {
    // Legacy parameter scheme for displacements
    if (isParamValid("disp_x"))
      _vars[0] = coupled("disp_x");
    if (isParamValid("disp_y"))
      _vars[1] = coupled("disp_y");
    if (isParamValid("disp_z"))
      _vars[2] = coupled("disp_z");
 
    mooseDeprecated("use the `displacements` parameter rather than the `disp_*` parameters (those "
                    "will go away with the deprecation of the Solid Mechanics module).");
  }
 
  _penetration_locator.setUpdate(false);
}

~GluedContactConstraint()

virtual GluedContactConstraint::~GluedContactConstraint ( )

inlinevirtual

Definition at line 31 of file GluedContactConstraint.h.

{}

Member Function Documentation

computeQpJacobian()

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

virtual

Reimplemented from SparsityBasedContactConstraint.

Definition at line 210 of file GluedContactConstraint.C.

{
  switch (type)
  {
    case Moose::SlaveSlave:
      return _penalty * _phi_slave[_j][_qp] * _test_slave[_i][_qp];
 
    case Moose::SlaveMaster:
      return _penalty * -_phi_master[_j][_qp] * _test_slave[_i][_qp];
 
    case Moose::MasterSlave:
    {
      const double slave_jac = (*_jacobian)(_current_node->dof_number(0, _vars[_component], 0),
                                            _connected_dof_indices[_j]);
      return slave_jac * _test_master[_i][_qp];
    }
 
    case Moose::MasterMaster:
      return 0.0;
 
    default:
      mooseError("Unhandled ConstraintJacobianType");
  }
 
  return 0.0;
}

computeQpOffDiagJacobian()

Real GluedContactConstraint::computeQpOffDiagJacobian ( Moose::ConstraintJacobianType  type, unsigned int  jvar  )

virtual

Compute off-diagonal Jacobian entries.

Parameters

type	The type of coupling
jvar	The index of the coupled variable

Definition at line 238 of file GluedContactConstraint.C.

{
  Real retVal = 0.0;
 
  switch (type)
  {
    case Moose::SlaveSlave:
    case Moose::SlaveMaster:
    case Moose::MasterMaster:
      break;
 
    case Moose::MasterSlave:
    {
      const double slave_jac = (*_jacobian)(_current_node->dof_number(0, _vars[_component], 0),
                                            _connected_dof_indices[_j]);
      retVal = slave_jac * _test_master[_i][_qp];
      break;
    }
 
    default:
      mooseError("Unhandled ConstraintJacobianType");
  }
 
  return retVal;
}

computeQpResidual()

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

virtual

Reimplemented from SparsityBasedContactConstraint.

Definition at line 175 of file GluedContactConstraint.C.

{
  switch (type)
  {
    case Moose::Slave:
    {
      PenetrationInfo * pinfo = _penetration_locator._penetration_info[_current_node->id()];
      Real distance = (*_current_node)(_component)-pinfo->_closest_point(_component);
      Real pen_force = _penalty * distance;
 
      Real resid = pen_force;
      pinfo->_contact_force(_component) = resid;
      pinfo->_mech_status = PenetrationInfo::MS_STICKING;
 
      return _test_slave[_i][_qp] * resid;
    }
 
    case Moose::Master:
    {
      PenetrationInfo * pinfo = _penetration_locator._penetration_info[_current_node->id()];
 
      dof_id_type dof_number = _current_node->dof_number(0, _vars[_component], 0);
      Real resid = _residual_copy(dof_number);
 
      pinfo->_contact_force(_component) = -resid;
      pinfo->_mech_status = PenetrationInfo::MS_STICKING;
 
      return _test_master[_i][_qp] * resid;
    }
  }
 
  return 0.0;
}

computeQpSlaveValue()

Real GluedContactConstraint::computeQpSlaveValue ( )

virtual

Reimplemented from SparsityBasedContactConstraint.

Definition at line 169 of file GluedContactConstraint.C.

{
  return _u_slave[_qp];
}

getConnectedDofIndices()

void SparsityBasedContactConstraint::getConnectedDofIndices ( )

virtualinherited

Gets the indices for all dofs conected to the constraint Get indices for all the slave Jacobian columns, not only those based on the mesh connectivity.

Definition at line 29 of file SparsityBasedContactConstraint.C.

{
#if defined(LIBMESH_HAVE_PETSC) && !PETSC_VERSION_LESS_THAN(3, 3, 0)
  _connected_dof_indices.clear();
 
  // An ugly hack: have to extract the row and look at it's sparsity structure, since otherwise I
  // won't get the dofs connected to this row by virtue of intervariable coupling
  // This is easier than sifting through all coupled variables, selecting those active on the
  // current subdomain, dealing with the scalar variables, etc.
  // Also, importantly, this will miss coupling to variables that might have introduced by prior
  // constraints similar to this one!
  PetscMatrix<Number> * petsc_jacobian = dynamic_cast<PetscMatrix<Number> *>(_jacobian);
  mooseAssert(petsc_jacobian, "Expected a PETSc matrix");
  Mat jac = petsc_jacobian->mat();
  PetscErrorCode ierr;
  PetscInt ncols;
  const PetscInt * cols;
  ierr = MatGetRow(jac, static_cast<PetscInt>(_var.nodalDofIndex()), &ncols, &cols, PETSC_NULL);
  CHKERRABORT(_communicator.get(), ierr);
  bool debug = false;
  if (debug)
  {
    libMesh::out << "_connected_dof_indices: adding " << ncols << " dofs from Jacobian row["
                 << _var.nodalDofIndex() << "] = [";
  }
  for (PetscInt i = 0; i < ncols; ++i)
  {
    if (debug)
    {
      libMesh::out << cols[i] << " ";
    }
    _connected_dof_indices.push_back(cols[i]);
  }
  if (debug)
  {
    libMesh::out << "]\n";
  }
  ierr = MatRestoreRow(jac, static_cast<PetscInt>(_var.nodalDofIndex()), &ncols, &cols, PETSC_NULL);
  CHKERRABORT(_communicator.get(), ierr);
#else
  NodeFaceConstraint::getConnectedDofIndices();
#endif
}

jacobianSetup()

void GluedContactConstraint::jacobianSetup ( )

virtual

Definition at line 121 of file GluedContactConstraint.C.

{
  if (_component == 0)
  {
    if (_updateContactSet)
      updateContactSet();
 
    _updateContactSet = true;
  }
}

shouldApply()

bool GluedContactConstraint::shouldApply

(

)

Definition at line 162 of file GluedContactConstraint.C.

{
  auto hpit = _penetration_locator._has_penetrated.find(_current_node->id());
  return hpit != _penetration_locator._has_penetrated.end();
}

timestepSetup()

void GluedContactConstraint::timestepSetup ( )

virtual

Definition at line 111 of file GluedContactConstraint.C.

{
  if (_component == 0)
  {
    updateContactSet(true);
    _updateContactSet = false;
  }
}

updateContactSet()

void GluedContactConstraint::updateContactSet ( bool  beginning_of_step = false )

virtual

Definition at line 133 of file GluedContactConstraint.C.

{
  auto & has_penetrated = _penetration_locator._has_penetrated;
 
  for (auto & pinfo_pair : _penetration_locator._penetration_info)
  {
    PenetrationInfo * pinfo = pinfo_pair.second;
 
    // Skip this pinfo if there are no DOFs on this node.
    if (!pinfo || pinfo->_node->n_comp(_sys.number(), _vars[_component]) < 1)
      continue;
 
    const dof_id_type slave_node_num = pinfo_pair.first;
    auto hpit = has_penetrated.find(slave_node_num);
 
    if (beginning_of_step)
    {
      pinfo->_starting_elem = pinfo->_elem;
      pinfo->_starting_side_num = pinfo->_side_num;
      pinfo->_starting_closest_point_ref = pinfo->_closest_point_ref;
    }
 
    if (pinfo->_distance >= 0)
      if (hpit == has_penetrated.end())
        has_penetrated.insert(slave_node_num);
  }
}

Referenced by jacobianSetup(), and timestepSetup().

Member Data Documentation

_aux_solution

const NumericVector<Number>* GluedContactConstraint::_aux_solution

protected

Definition at line 73 of file GluedContactConstraint.h.

_aux_system

SystemBase& GluedContactConstraint::_aux_system

protected

Definition at line 72 of file GluedContactConstraint.h.

_component

const unsigned int GluedContactConstraint::_component

protected

Definition at line 54 of file GluedContactConstraint.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), computeQpResidual(), jacobianSetup(), timestepSetup(), and updateContactSet().

_formulation

const ContactFormulation GluedContactConstraint::_formulation

protected

Definition at line 56 of file GluedContactConstraint.h.

_friction_coefficient

const Real GluedContactConstraint::_friction_coefficient

protected

Definition at line 59 of file GluedContactConstraint.h.

_model

const ContactModel GluedContactConstraint::_model

protected

Definition at line 55 of file GluedContactConstraint.h.

_nodal_area_var

MooseVariable* GluedContactConstraint::_nodal_area_var

protected

Definition at line 71 of file GluedContactConstraint.h.

_penalty

const Real GluedContactConstraint::_penalty

protected

Definition at line 58 of file GluedContactConstraint.h.

Referenced by computeQpJacobian(), and computeQpResidual().

_residual_copy

NumericVector<Number>& GluedContactConstraint::_residual_copy

protected

Definition at line 63 of file GluedContactConstraint.h.

Referenced by computeQpResidual().

_tension_release

const Real GluedContactConstraint::_tension_release

protected

Definition at line 60 of file GluedContactConstraint.h.

_updateContactSet

bool GluedContactConstraint::_updateContactSet

protected

Definition at line 61 of file GluedContactConstraint.h.

Referenced by jacobianSetup(), and timestepSetup().

_vars

std::vector<unsigned int> GluedContactConstraint::_vars

protected

Definition at line 69 of file GluedContactConstraint.h.

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

_x_var

unsigned int GluedContactConstraint::_x_var

protected

Definition at line 65 of file GluedContactConstraint.h.

_y_var

unsigned int GluedContactConstraint::_y_var

protected

Definition at line 66 of file GluedContactConstraint.h.

_z_var

unsigned int GluedContactConstraint::_z_var

protected

Definition at line 67 of file GluedContactConstraint.h.

---

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

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