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 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;

  std::unordered_map<const DofObject *, std::pair<TwoVector, TwoVector>> _dof_to_step_slip;

  std::unordered_map<const DofObject *, std::pair<TwoVector, TwoVector>> _dof_to_accumulated_slip;

  std::unordered_map<const DofObject *, std::pair<ADTwoVector, TwoVector>>
      _dof_to_tangential_traction;

  ADVariableValue _frictional_contact_traction_one;

  ADVariableValue _frictional_contact_traction_two;

  std::unordered_map<const DofObject *, TwoVector> _dof_to_frictional_lagrange_multipliers;

  std::unordered_map<const DofObject *, Real> _dof_to_local_penalty_friction;

  const Real _penalty_multiplier_friction;

  const enum class AdaptivityFrictionalPenalty { SIMPLE, FRICTION_LIMIT } _adaptivity_friction;

  const Real _epsilon_tolerance;
};