PenaltyFrictionUserObject.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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
 
#include "PenaltyWeightedGapUserObject.h"
#include "WeightedVelocitiesUserObject.h"
#include "AugmentedLagrangeInterface.h"
#include "TwoVector.h"
 
class PenaltyFrictionUserObject : virtual public PenaltyWeightedGapUserObject,
                                  virtual public WeightedVelocitiesUserObject
{
public:
  static InputParameters validParams();
 
  PenaltyFrictionUserObject(const InputParameters & parameters);
 
  virtual const ADVariableValue & contactTangentialPressureDirOne() const override;
  virtual const ADVariableValue & contactTangentialPressureDirTwo() const override;
 
  virtual void initialize() override;
  virtual void finalize() override;
  virtual void reinit() override;
  virtual ADReal nodalTangentialPressure(const Node & node, unsigned int direction) const override;
  virtual void timestepSetup() override;
 
  virtual Real getFrictionalContactPressure(const Node * const node,
                                            const unsigned int component) const override;
  virtual Real getAccumulatedSlip(const Node * const node,
                                  const unsigned int component) const override;
  virtual Real getTangentialVelocity(const Node * const node,
                                     const unsigned int component) const override;
  virtual Real getDeltaTangentialLagrangeMultiplier(const Node * const node,
                                                    const unsigned int component) const override;
 
  virtual bool isAugmentedLagrangianConverged() override;
  virtual void updateAugmentedLagrangianMultipliers() override;
 
protected:
  virtual const VariableTestValue & test() const override;
  virtual bool constrainedByOwner() const override { return false; }
 
  const Real _penalty;
 
  const Real _penalty_friction;
 
  const Real _slip_tolerance;
 
  const Real _friction_coefficient;
 
  // The maps below are keyed by node id so that they can be stored as restartable data. Recovery
  // must restore their keys as well as their values: the augmented Lagrange updates loop over every
  // node a map holds, including nodes that were in contact only in earlier timesteps.
 
  std::unordered_map<dof_id_type, std::pair<TwoVector, TwoVector>> & _dof_to_step_slip;
 
  std::unordered_map<dof_id_type, std::pair<TwoVector, TwoVector>> & _dof_to_accumulated_slip;
 
  std::unordered_map<dof_id_type, std::pair<ADTwoVector, TwoVector>> & _dof_to_tangential_traction;
 
  ADVariableValue _frictional_contact_traction_one;
 
  ADVariableValue _frictional_contact_traction_two;
 
  std::unordered_map<dof_id_type, TwoVector> & _dof_to_frictional_lagrange_multipliers;
 
  std::unordered_map<dof_id_type, Real> & _dof_to_local_penalty_friction;
 
  const Real _penalty_multiplier_friction;
 
  const enum class AdaptivityFrictionalPenalty { SIMPLE, FRICTION_LIMIT } _adaptivity_friction;
 
  const Real _epsilon_tolerance;
 
  int & _t_step_old_friction;
 
  std::unordered_map<dof_id_type, std::pair<TwoVector, TwoVector>> _dof_to_step_slip_start;
  std::unordered_map<dof_id_type, std::pair<TwoVector, TwoVector>> _dof_to_accumulated_slip_start;
  std::unordered_map<dof_id_type, std::pair<ADTwoVector, TwoVector>>
      _dof_to_tangential_traction_start;
  std::unordered_map<dof_id_type, TwoVector> _dof_to_frictional_lagrange_multipliers_start;
  std::unordered_map<dof_id_type, Real> _dof_to_local_penalty_friction_start;
};