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