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