ComputeDynamicWeightedGapLMMechanicalContact.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 "ADMortarConstraint.h"

#include <unordered_map>

class ComputeDynamicWeightedGapLMMechanicalContact : public ADMortarConstraint
{
public:
  static InputParameters validParams();

  ComputeDynamicWeightedGapLMMechanicalContact(const InputParameters & parameters);

protected:
  virtual void computeQpProperties();

  virtual void computeQpIProperties();

  ADReal computeQpResidual(Moose::MortarType mortar_type) final;

  using ADMortarConstraint::computeResidual;
  void computeResidual(Moose::MortarType mortar_type) override;
  using ADMortarConstraint::computeJacobian;
  void computeJacobian(Moose::MortarType mortar_type) override;
  void residualSetup() override;
  void jacobianSetup() override final;

  void timestepSetup() override;

  virtual void post() override;

  virtual void
  incorrectEdgeDroppingPost(const std::unordered_set<const Node *> & inactive_lm_nodes) override;

  void communicateWear();

  virtual void enforceConstraintOnDof(const DofObject * const dof);

  const ADVariableValue & _secondary_disp_x;
  const ADVariableValue & _primary_disp_x;
  const ADVariableValue & _secondary_disp_y;
  const ADVariableValue & _primary_disp_y;

  const bool _has_disp_z;
  const ADVariableValue * const _secondary_disp_z;
  const ADVariableValue * const _primary_disp_z;

  const Real _c;

  ADReal _qp_gap;

  Real _qp_factor;

  bool _normalize_c;

  const bool _nodal;

  const MooseVariable * const _disp_x_var;
  const MooseVariable * const _disp_y_var;
  const MooseVariable * const _disp_z_var;

  ADRealVectorValue _qp_gap_nodal;

  std::unordered_map<const DofObject *, std::pair<ADReal, Real>> _dof_to_weighted_gap;

  const ADReal * _weighted_gap_ptr = nullptr;
  const Real * _normalization_ptr = nullptr;

  const Real _capture_tolerance;

  const ADVariableValue & _secondary_x_dot;
  const ADVariableValue & _primary_x_dot;
  const ADVariableValue & _secondary_y_dot;
  const ADVariableValue & _primary_y_dot;

  const ADVariableValue * _secondary_z_dot;
  const ADVariableValue * _primary_z_dot;

  const bool _has_wear;

  const VariableValue & _wear_depth;

  std::unordered_map<const DofObject *, ADReal> _dof_to_old_weighted_gap;

  ADRealVectorValue _qp_gap_nodal_dynamics;

  ADRealVectorValue _qp_velocity;

  std::unordered_map<const DofObject *, ADReal> _dof_to_weighted_gap_dynamics;

  std::unordered_map<const DofObject *, ADReal> _dof_to_velocity;
  std::unordered_map<const DofObject *, ADReal> _dof_to_old_velocity;

  // Newmark-beta beta parameter
  const Real _newmark_beta;

  // Newmark-beta gamma parameter
  const Real _newmark_gamma;

  std::unordered_map<const DofObject *, ADReal> _dof_to_nodal_wear_depth;

  std::unordered_map<const DofObject *, ADReal> _dof_to_nodal_old_wear_depth;

  ADRealVectorValue _relative_velocity;
};