HertzSpring.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_HertzSpring_hpp_
#define _libDEM_HertzSpring_hpp_
 
#include <string>
#include <vector>
 
#include "GranOO3/DEM/Bond.hpp"
#include "GranOO3/Core/SetOf.hpp"
#include "GranOO3/DEM/Element.hpp"
#include "GranOO3/Geom/Point.hpp"
 
namespace GranOO3
{
  namespace DEM
  {
    class HertzSpring : public Bond,
            public Core::Register<HertzSpring>,
            public Physic::CriticalTimeStep {
 
      GRANOO_CLASS(DEM::HertzSpring, Bond);
 
    public:
      //CONSTRUCTORS & DESTRUCTORS
      HertzSpring(Element& el1, Element& el2);
      HertzSpring(Element& el1, Element& el2, double young_modulus, double poisson_ratio);
      virtual ~HertzSpring();
 
      //USEFULL
      void init();
      void process();
      virtual void update();
      virtual void compute_load();
      virtual double get_linear_stiffness() const;
      virtual double get_angular_stiffness() const;
      virtual double compute_critical_time_step() const;
      virtual std::string info() const;
 
      //ACCESSORS
      GRANOO_ACCESS    (stiffness        , double       , _stiffness         );
      GRANOO_ACCESS    (young_modulus    , double       , _young_modulus     );
      GRANOO_ACCESS_GET(poisson_ratio    , double       , _poisson_ratio     );
      GRANOO_ACCESS_GET(direction        , Geom::Vector , _direction         );
      GRANOO_ACCESS_GET(force            , Geom::Vector , _force             );
      GRANOO_ACCESS_GET(equivalent_radius, double       , _equivalent_radius );
      GRANOO_ACCESS_GET(contact_radius   , double       , _contact_radius    );
      GRANOO_ACCESS_GET(contact_surface  , double       , _contact_surface   );
 
      
    private:
      HertzSpring() = delete;
      HertzSpring(const HertzSpring&) = delete;
      HertzSpring & operator=(const HertzSpring &) = delete;
 
    private:
      // BOOST SERIALIZATION
      friend class boost::serialization::access;
      template<class Archive> void serialize(Archive &, const unsigned int);
      
    private:
      //Initial Parameter
      double _young_modulus;
      double _poisson_ratio;
      //Parameter
      Geom::Vector _direction;
      Geom::Vector _force;
      double _equivalent_radius;
      double _stiffness;
      double _contact_radius;
      double _contact_surface;
 
   
    };
 
    // BOOST SERIALIZATION
    template<class Archive> void
    HertzSpring::serialize(Archive &ar, const unsigned int version) {
      ar & boost::serialization::base_object<Bond>(*this);
      ar & _young_modulus;
      ar & _poisson_ratio;
    }
  }
}
 
 
#include <boost/serialization/version.hpp>
BOOST_CLASS_VERSION(GranOO3::DEM::HertzSpring, 0)
 
 
// need placement new snippet because no default constructor exist
namespace boost
{
  namespace serialization
  {
    template<class Archive>
    inline void save_construct_data(Archive &ar, const GranOO3::DEM::HertzSpring *t, const unsigned int) {
 
      const GranOO3::DEM::Element *el1 = &t->get_element1();
      const GranOO3::DEM::Element *el2 = &t->get_element2();
      ar << el1;
      ar << el2;
    }
 
    template<class Archive>
    inline void load_construct_data(Archive &ar, GranOO3::DEM::HertzSpring *t, const unsigned int) {
 
      GranOO3::DEM::Element* el1 = nullptr;
      GranOO3::DEM::Element* el2 = nullptr;
      ar >> el1;
      ar >> el2;
      ::new(t)GranOO3::DEM::HertzSpring(*el1,*el2);
    }
  }
} // namespace ...
 
 
namespace GranOO3 {
  GRANOO_CLASS_DECLARE_TPL(DEM::HertzSpring);
}
 
 
#endif