SinterNeck_Martin.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_SinterNeck_Martin_hpp_
#define _libDEM_SinterNeck_Martin_hpp_
 
#include <string>
#include <vector>
 
 
 
 
#include "GranOO3/DEM/Bond/SinterNeck.hpp"
 
namespace GranOO3
{
  namespace DEM
  {
 
 
    class SinterNeck_Martin : public SinterNeck,
                 public Core::Register< SinterNeck_Martin >
    {
 
      GRANOO_CLASS(DEM::SinterNeck_Martin, SinterNeck);
 
 
    public:
      //CONSTRUCTORS & DESTRUCTORS
      SinterNeck_Martin(Element& el1, Element& el2);
      virtual ~SinterNeck_Martin();
 
      //USEFULL
      virtual void compute_load();
      virtual std::string info() const;
 
      // UPDATE PARAMETERS
      void update_delta_b();
      
      //ACCESSORS
      GRANOO_ACCESS(delta_b_x_D0b, double, _delta_b_x_D0b ); // initial diffusion flux
      GRANOO_ACCESS(Qb           , double, _Qb            ); // activation energy 
      GRANOO_ACCESS(beta         , double, _beta          ); // model parameter
      GRANOO_ACCESS(alpha        , double, _alpha         ); // model parameter
      GRANOO_ACCESS(vt           , double, _vt            ); // tangential relative velocity
      GRANOO_ACCESS(eta          , double, _eta           ); // dimensionless viscous parameter
      GRANOO_ACCESS(Ft           , double, _Ft            ); // tangential viscous force
      GRANOO_ACCESS(t      , Geom::Vector, _t             ); // tangential direction
      
    private:
      SinterNeck_Martin() = delete;
      SinterNeck_Martin(const SinterNeck_Martin& frame) = delete;
      SinterNeck_Martin & operator=(const SinterNeck_Martin &) = delete;
 
    private:
      //SERIALIZATION
      friend class boost::serialization::access;
      template<class Archive> void save(Archive&, const unsigned int ) const;
      template<class Archive> void load(Archive&, const unsigned int);
      BOOST_SERIALIZATION_SPLIT_MEMBER();
 
    protected:
      double _delta_b_x_D0b = 0.;  // initial diffusion flux (δb.D0b: Martin)
      double _Qb            = 0.;  // activation energy      (    Qb: Martin)
 
      double _beta  = 4.0; // model parameter (α: Martin)
      double _alpha = 4.5; // model parameter (β: Martin)
 
      double _vt = 0.; // tangential relative velocity (du/dt: Martin)
      
      double       _eta = 0.;                     // dimensionless viscous parameter (η: Martin)
      double       _Ft  = 0.;                     // tangential viscous force
 
      Geom::Vector _t = Geom::Vector(0.,0.,0.);   // tangential direction
      
    };
 
 
    template<class Archive> void
    SinterNeck_Martin::save(Archive& ar, const unsigned int v) const {
      ar << boost::serialization::base_object<SinterNeck>(*this);
      ar << _delta_b_x_D0b;
      ar << _Qb;
      ar << _beta;
      ar << _alpha;
      ar << _vt;
      ar << _eta;
      ar << _Ft;
      ar << _t;
    }
 
    template<class Archive> void
    SinterNeck_Martin::load(Archive& ar, const unsigned int v) {
      ar >> boost::serialization::base_object<SinterNeck>(*this);
      ar >> _delta_b_x_D0b;
      ar >> _Qb;
      ar >> _beta;
      ar >> _alpha;
      ar >> _vt;
      ar >> _eta;
      ar >> _Ft;
      ar >> _t;
    }
 
 
  }
}
 
 
#include <boost/serialization/version.hpp>
BOOST_CLASS_VERSION(GranOO3::DEM::SinterNeck_Martin,0)
 
// need placement new snippet because no default constructor exist
namespace boost
{
  namespace serialization
  {
 
    template<class Archive>
    void save_construct_data(Archive & ar,
                             const GranOO3::DEM::SinterNeck_Martin * t,
                             const unsigned int) {
      const GranOO3::DEM::Element* de1 = &t->get_element1();
      const GranOO3::DEM::Element* de2 = &t->get_element2();
      ar << de1;
      ar << de2;
    }
 
    template<class Archive>
    void load_construct_data(Archive & ar,
                             GranOO3::DEM::SinterNeck_Martin * t,
                             const unsigned int) {
      GranOO3::DEM::Element* de1 = nullptr;
      GranOO3::DEM::Element* de2 = nullptr;
      ar >> de1;
      ar >> de2;
      ::new(t)GranOO3::DEM::SinterNeck_Martin(*de1,*de2);
    }
 
  }
} // namespace ...
 
namespace GranOO3
{
  GRANOO_CLASS_DECLARE_TPL(DEM::SinterNeck_Martin);
}
 
 
#endif