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