PenaltySimpleCohesiveZoneModel.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 "TwoVector.h"
#include "PenaltyWeightedGapUserObject.h"
#include "WeightedVelocitiesUserObject.h"
#include "AugmentedLagrangeInterface.h"

class PenaltySimpleCohesiveZoneModel : virtual public PenaltyWeightedGapUserObject,
                                       virtual public WeightedVelocitiesUserObject
{
public:
  static InputParameters validParams();

  PenaltySimpleCohesiveZoneModel(const InputParameters & parameters);

  virtual const ADVariableValue & contactTangentialPressureDirOne() const override;
  virtual const ADVariableValue & contactTangentialPressureDirTwo() const override;

  virtual void initialize() override;
  virtual void finalize() override;
  virtual void reinit() override;
  virtual void timestepSetup() override;

protected:
  virtual const VariableTestValue & test() const override;
  virtual bool constrainedByOwner() const override { return false; }

  // Compute CZM kinematics.
  virtual void prepareJumpKinematicQuantities() {}
  virtual void computeFandR(const Node * const /*node*/) {}

  virtual void computeBilinearMixedModeTraction(const Node * const /*node*/) {}

  virtual void computeGlobalTraction(const Node * const /*node*/) {}

  const Real _penalty;

  const Real _penalty_friction;

  const Real _friction_coefficient;

  std::unordered_map<const DofObject *, std::pair<TwoVector, TwoVector>> _dof_to_step_slip;

  std::unordered_map<const DofObject *, std::pair<TwoVector, TwoVector>> _dof_to_accumulated_slip;

  std::unordered_map<const DofObject *, std::pair<ADTwoVector, TwoVector>>
      _dof_to_tangential_traction;

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

  ADVariableValue _frictional_contact_traction_one;

  ADVariableValue _frictional_contact_traction_two;

  std::unordered_map<const DofObject *, TwoVector> _dof_to_frictional_lagrange_multipliers;

  std::unordered_map<const DofObject *, Real> _dof_to_local_penalty_friction;

  const Real _epsilon_tolerance;
};