MechanicalContactConstraint.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
// MOOSE includes
#include "NodeFaceConstraint.h"
#include "PenetrationLocator.h"
 
// Forward Declarations
class MechanicalContactConstraint;
class ContactLineSearchBase;
enum class ContactModel;
enum class ContactFormulation;
 
template <>
InputParameters validParams<MechanicalContactConstraint>();
 
class MechanicalContactConstraint : public NodeFaceConstraint
{
public:
  MechanicalContactConstraint(const InputParameters & parameters);
 
  virtual void timestepSetup() override;
  virtual void jacobianSetup() override;
  virtual void residualEnd() override;
 
  virtual bool AugmentedLagrangianContactConverged();
 
  virtual void updateAugmentedLagrangianMultiplier(bool beginning_of_step = false);
 
  virtual void updateContactStatefulData(bool beginning_of_step = false);
 
  virtual Real computeQpSlaveValue() override;
 
  virtual Real computeQpResidual(Moose::ConstraintType type) override;
 
  virtual void computeJacobian() override;
 
  virtual void computeOffDiagJacobian(unsigned int jvar) override;
 
  virtual Real computeQpJacobian(Moose::ConstraintJacobianType type) override;
 
  virtual Real computeQpOffDiagJacobian(Moose::ConstraintJacobianType type,
                                        unsigned int jvar) override;
 
  virtual void getConnectedDofIndices(unsigned int var_num) override;
 
  bool getCoupledVarComponent(unsigned int var_num, unsigned int & component);
 
  virtual bool addCouplingEntriesToJacobian() override { return _master_slave_jacobian; }
 
  bool shouldApply() override;
  void computeContactForce(PenetrationInfo * pinfo, bool update_contact_set);
 
protected:
  MooseSharedPointer<DisplacedProblem> _displaced_problem;
  Real nodalArea(PenetrationInfo & pinfo);
  Real getPenalty(PenetrationInfo & pinfo);
  Real getTangentialPenalty(PenetrationInfo & pinfo);
 
  const unsigned int _component;
  const ContactModel _model;
  const ContactFormulation _formulation;
  const bool _normalize_penalty;
 
  const Real _penalty;
  Real _penalty_tangential;
  const Real _friction_coefficient;
  const Real _tension_release;
  const Real _capture_tolerance;
  const unsigned int _stick_lock_iterations;
  const Real _stick_unlock_factor;
  bool _update_stateful_data;
 
  NumericVector<Number> & _residual_copy;
  //  std::map<Point, PenetrationInfo *> _point_to_info;
 
  const unsigned int _mesh_dimension;
 
  std::vector<unsigned int> _vars;
  std::vector<MooseVariable *> _var_objects;
 
  MooseVariable * _nodal_area_var;
  SystemBase & _aux_system;
  const NumericVector<Number> * _aux_solution;
 
  const bool _master_slave_jacobian;
  const bool _connected_slave_nodes_jacobian;
  const bool _non_displacement_vars_jacobian;
 
  Real _al_penetration_tolerance;
  Real _al_incremental_slip_tolerance;
  Real _al_frictional_force_tolerance;
 
  ContactLineSearchBase * _contact_linesearch;
  std::set<dof_id_type> _current_contact_state;
  std::set<dof_id_type> _old_contact_state;
 
  const bool _print_contact_nodes;
  static Threads::spin_mutex _contact_set_mutex;
};