ContactLaw.hpp

Go to the documentation of this file.

// This file is part of GranOO, a workbench for DEM simulation.
//
// Author(s)    : - Damien Andre       IRCER/UNILIM, Limoges France
//                  <damien.andre@unilim.fr>
//                - Jean-luc Charles   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jean-luc.charles@ensam.eu>
//                - Jeremie Girardot   Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <jeremie.girardot@ensam.eu>
//                - Cedric Hubert      LAMIH/UPHF, Valenciennes France
//                  <cedric.hubert@uphf.fr>
//                - Ivan Iordanoff     Arts et Metiers ParisTech, CNRS, I2M, Bordeaux France
//                  <ivan.iordanoff@ensam.eu>
//
// Copyright (C) 2008-2019 D. Andre, JL. Charles, J. Girardot, C. Hubert, I. Iordanoff
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program.  If not, see <http://www.gnu.org/licenses/>.
 
 
#ifndef _libDEM_ContactLaw_hpp_
#define _libDEM_ContactLaw_hpp_
 
#include <iostream>
 
#include "GranOO3/DEM/DiscreteElement.hpp"
#include "GranOO3/Collision/Manager.hpp"
#include "GranOO3/Physic/CriticalTimeStep.hpp"
 
#define APPLY_FORCE(F1, F2, energyLabel)                             \
  de1.force() += F1;                                                 \
  de2.force() += F2;                                                 \
  if (Physic::EnergyBalance::get().is_enable())                      \
  {                                                               \
    Physic::EnergyBalance& energy = Physic::EnergyBalance::get(); \
    energy.add_force_on(de1, energyLabel, F1);                            \
    energy.add_force_on(de2, energyLabel, F2);                            \
  }
 
 
 
 
#define APPLY_TORQUE(T1, T2, energyLabel)                            \
  de1.torque() += T1;                                                \
  de2.torque() += T2;                                                \
  if (Physic::EnergyBalance::get().is_enable())                      \
  {                                                               \
    Physic::EnergyBalance& energy = Physic::EnergyBalance::get(); \
    energy.add_torque_on(de1, energyLabel, T1);                   \
    energy.add_torque_on(de2, energyLabel, T2);                   \
  }
 
 
namespace GranOO3
{
  namespace DEM
  {
 
    template<class T>
    class ContactLaw : public Collision::CallBack<DEM::DiscreteElement, T>
    {
    public:
      static std::string class_ID() {return "DEM::ContactLaw_" + T::class_ID() + "_";}
 
    protected:
      //CONSTRUCTORS & DESTRUCTORS
      ContactLaw();
      virtual ~ContactLaw();
 
    public:
      //USEFULL
      virtual void parse_xml() = 0;
      virtual void init();
      virtual void info(std::ostream &) const;
      virtual void compute_reaction(DEM::DiscreteElement& de1, T& de2,
                                    const Geom::Vector& normal, const double& penetration) = 0;
 
      virtual void pre_contact_detection();
      virtual void post_contact_detection();
      void manage_collision(DEM::DiscreteElement&, T&, const Collision::Data&) {
      };                                                                          // TODO be coherent
 
    private:
      ContactLaw(const ContactLaw&) = delete;
      ContactLaw& operator=(const ContactLaw&) = delete;
 
    };
 
 
    template<class T>
    ContactLaw<T>::ContactLaw() {
    }
 
    template<class T>
    ContactLaw<T>::~ContactLaw() {
    }
 
    template<class T> void
    ContactLaw<T>::init() {
    }
 
    template<class T> void
    ContactLaw<T>::pre_contact_detection() {
    }
 
    template<class T> void
    ContactLaw<T>::post_contact_detection() {
    }
 
    template<class T> void
    ContactLaw<T>::info(std::ostream & out) const {
      out << "ContactLaw::Infos:  None" << granoo::endl;
    }
 
  }
}
 
 
#endif