ExplicitDynamicsContactConstraint.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

// MOOSE includes
#include "NodeFaceConstraint.h"
#include "PenetrationLocator.h"
#include "TwoMaterialPropertyInterface.h"
#include "Coupleable.h"

// Forward Declarations
enum class ExplicitDynamicsContactModel;
class DisplacedProblem;

class ExplicitDynamicsContactConstraint : public NodeFaceConstraint,
                                          public TwoMaterialPropertyInterface
{
public:
  static InputParameters validParams();

  ExplicitDynamicsContactConstraint(const InputParameters & parameters);

  virtual void timestepSetup() override;
  virtual void jacobianSetup() override {}
  virtual void residualEnd() override {}

  virtual void updateContactStatefulData(bool beginning_of_step);
  virtual Real computeQpSecondaryValue() override;
  virtual Real computeQpResidual(Moose::ConstraintType type) override;

  virtual void computeJacobian() override {}

  virtual void computeOffDiagJacobian(unsigned int /*jvar*/) override {}

  virtual Real computeQpJacobian(Moose::ConstraintJacobianType /*type*/) override { return 0.0; }

  virtual Real computeQpOffDiagJacobian(Moose::ConstraintJacobianType /*type*/,
                                        unsigned int /*jvar*/) override
  {
    return 0.0;
  }

  bool shouldApply() override;
  void computeContactForce(const Node & node, PenetrationInfo * pinfo, bool update_contact_set);
  virtual bool isExplicitConstraint() const override { return true; }
  virtual bool addCouplingEntriesToJacobian() override { return false; }

  virtual const std::unordered_set<unsigned int> & getMatPropDependencies() const override;

  virtual void overwriteBoundaryVariables(NumericVector<Number> &, const Node &) const override
  {
    return;
  };

protected:
  void solveImpactEquations(const Node & node,
                            PenetrationInfo * pinfo,
                            const RealVectorValue & distance_gap);

  MooseSharedPointer<DisplacedProblem> _displaced_problem;
  Real gapOffset(const Node & node);
  Real getPenalty(const Node & node);

  const unsigned int _component;
  const ExplicitDynamicsContactModel _model;

  bool _update_stateful_data;

  const unsigned int _mesh_dimension;

  std::vector<unsigned int> _vars;
  std::vector<MooseVariable *> _var_objects;

  const bool _has_secondary_gap_offset;
  const MooseVariable * const _secondary_gap_offset_var;
  const bool _has_mapped_primary_gap_offset;
  const MooseVariable * const _mapped_primary_gap_offset_var;

  const Real _penalty;

  const bool _print_contact_nodes;

  const static unsigned int _no_iterations;

  NumericVector<Number> & _residual_copy;

  MooseWritableVariable * _gap_rate;

  const VariableValue & _neighbor_vel_x;
  const VariableValue & _neighbor_vel_y;
  const VariableValue & _neighbor_vel_z;

private:
  std::unordered_map<dof_id_type, Real> _dof_to_position;

  Real computeFaceMass(const NumericVector<Real> & lumped_mass);
};

inline const std::unordered_set<unsigned int> &
ExplicitDynamicsContactConstraint::getMatPropDependencies() const
{
  return TwoMaterialPropertyInterface::getMatPropDependencies();
}