ComputeWeightedGapCartesianLMMechanicalContact.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 ComputeWeightedGapCartesianLMMechanicalContact : public ADMortarConstraint
{
public:
  static InputParameters validParams();
 
  ComputeWeightedGapCartesianLMMechanicalContact(const InputParameters & parameters);
  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 post() override;
 
  void
  incorrectEdgeDroppingPost(const std::unordered_set<const Node *> & inactive_lm_nodes) override;
 
protected:
  ADReal computeQpResidual(Moose::MortarType mortar_type) final;
 
  virtual void computeQpProperties();
 
  virtual void computeQpIProperties();
 
  void timestepSetup() override;
  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;
 
  std::vector<MooseVariable *> _lm_vars;
 
  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;
 
  std::unordered_map<const DofObject *, RealVectorValue> _dof_to_normal_vector;
 
  std::unordered_map<const DofObject *, RealVectorValue> _dof_to_old_normal_vector;
 
  std::unordered_map<const DofObject *, std::array<RealVectorValue, 2>> _dof_to_tangent_vectors;
 
  const ADReal * _weighted_gap_ptr = nullptr;
  const Real * _normalization_ptr = nullptr;
};