ComputeDynamicFrictionalForceLMMechanicalContact.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 "ComputeDynamicWeightedGapLMMechanicalContact.h"
 
class ComputeDynamicFrictionalForceLMMechanicalContact
  : public ComputeDynamicWeightedGapLMMechanicalContact
{
public:
  static InputParameters validParams();
 
  ComputeDynamicFrictionalForceLMMechanicalContact(const InputParameters & parameters);
  void residualSetup() override;
  void post() override;
 
  void
  incorrectEdgeDroppingPost(const std::unordered_set<const Node *> & inactive_lm_nodes) override;
 
protected:
  virtual void computeQpProperties() override;
 
  virtual void computeQpIProperties() override;
 
  virtual void enforceConstraintOnDof(const DofObject * const dof) override;
 
  virtual void enforceConstraintOnDof3d(const DofObject * const dof);
 
  void timestepSetup() override;
 
  ADReal computeFrictionValue(const ADReal & contact_pressure,
                              const Real & tangential_vel,
                              const Real & tangential_vel_dir);
 
  std::unordered_map<const DofObject *, std::array<ADReal, 2>> _dof_to_weighted_tangential_velocity;
 
  std::unordered_map<const DofObject *, std::array<Real, 2>> _dof_to_real_tangential_velocity;
 
  std::unordered_map<const DofObject *, std::array<Real, 2>> _dof_to_old_real_tangential_velocity;
 
  std::array<const ADReal *, 2> _tangential_vel_ptr = {{nullptr, nullptr}};
 
  ADRealVectorValue _qp_tangential_velocity_nodal;
 
  ADRealVectorValue _qp_real_tangential_velocity_nodal;
 
  const Real _c_t;
 
  std::vector<MooseVariable *> _friction_vars;
 
  const ADVariableValue & _secondary_x_dot;
 
  const ADVariableValue & _primary_x_dot;
 
  const ADVariableValue & _secondary_y_dot;
 
  const ADVariableValue & _primary_y_dot;
 
  const ADVariableValue * const _secondary_z_dot;
 
  const ADVariableValue * const _primary_z_dot;
 
  const Real _epsilon;
 
  const Real _mu;
 
  const Function * const _function_friction;
 
  const bool _has_friction_function;
 
  bool _3d;
};