contact_interaction.hpp

Go to the documentation of this file.

// Copyright (c) Lawrence Livermore National Security, LLC and
// other Smith Project Developers. See the top-level LICENSE file for
// details.
//
// SPDX-License-Identifier: (BSD-3-Clause)
 
#pragma once
 
#include "smith/smith_config.hpp"
 
#ifdef SMITH_USE_TRIBOL
 
#include <memory>
#include <set>
 
#include "mfem.hpp"
#include "tribol/common/Parameters.hpp"
 
#include "smith/physics/contact/contact_config.hpp"
#include "smith/physics/state/finite_element_dual.hpp"
#include "smith/physics/state/finite_element_state.hpp"
 
namespace smith {
 
class ContactInteraction {
 public:
  ContactInteraction(int interaction_id, const mfem::ParMesh& mesh, const std::set<int>& bdry_attr_surf1,
                     const std::set<int>& bdry_attr_surf2, const mfem::ParGridFunction& current_coords,
                     ContactOptions contact_opts);
 
  int getInteractionId() const { return interaction_id_; }
 
  const ContactOptions& getContactOptions() const { return contact_opts_; }
 
  FiniteElementDual forces() const;
 
  FiniteElementState pressure() const;
 
  FiniteElementDual gaps() const;
 
  std::unique_ptr<mfem::BlockOperator> jacobian() const;
 
  std::unique_ptr<mfem::BlockOperator> jacobianContribution() const;
 
  mfem::ParFiniteElementSpace& pressureSpace() const;
 
  void setPressure(const FiniteElementState& pressure) const;
 
  int numPressureDofs() const;
 
  const mfem::Array<int>& inactiveDofs() const;
 
  const mfem::Array<int>& inactiveDofs(const FiniteElementState& pressure) const;
 
  void evalJacobian(bool eval) const;
 
 private:
  tribol::ContactMethod getMethod() const;
 
  int interaction_id_;
 
  ContactOptions contact_opts_;
 
  const mfem::ParGridFunction& current_coords_;
 
  mutable mfem::Array<int> inactive_tdofs_;
};
 
}  // namespace smith
 
#endif