FlatBond.hpp

Go to the documentation of this file.

 
// This file is part of GranOO, a workbench for DEM simulation.
//
// Author(s)    : - Jean-luc CHARLES I2M-DuMAS/ENSAM Talence France
//                  <jean-luc.charles@ensam.eu>
//                - Damien ANDRE     SPCTS/ENS Ceramique industrielle, Limoges France
//                  <damien.andre@unilim.fr>
//                - Jeremie GIRARDOT I2M-DuMAS/ENSAM Talence France
//                  <jeremie.girardot@ensam.eu>
//                - Cedric Hubert    LAMIH/UVHC Valenciennes France
//                  <cedric.hubert@univ-valenciennes.fr>
//                - Ivan IORDANOFF   I2M-DuMAS-MPI/ENSAM Talence France
//                  <ivan.iordanoff@ensam.eu>
//
// Copyright (C) 2008-2016 JL. CHARLES, D. ANDRE, I. IORDANOFF, J. GIRARDOT
//
//
//
//
//
// 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_FlatBond_hpp_
#define _libDEM_FlatBond_hpp_
 
#include <string>
#include <vector>
 
 
 
#include "GranOO3/DEM/Bond.hpp"
#include "GranOO3/Core/SetOf.hpp"
#include "GranOO3/DEM/Element.hpp"
#include "GranOO3/DEM/Prop/VoronoiBond.hpp"
 
#include "GranOO3/Geom/SymTensor.hpp"
#include "GranOO3/Geom/Point.hpp" 
 
namespace GranOO3
{
  namespace DEM
  {
 
    
    class FlatBond : public Bond,
             public Core::Register<FlatBond >,
             public Physic::CriticalTimeStep
    {
 
      GRANOO_CLASS(DEM::FlatBond, Bond);
 
    public:
      enum DrawMode{NO_INTERFACE, INTERFACE};
      static DrawMode draw_mode;
      static bool update_current_stress;
      static bool update_current_strain;
 
    public:
      //CONSTRUCTORS & DESTRUCTORS
      FlatBond(Element& el1, Element& el2, bool build_voronoi_bond = true);
      FlatBond(Element& el1, Element& el2, double young_modulus, double poisson_ratio, bool build_voronoi_bond = true);
      virtual ~FlatBond();
 
      //USEFULL
      void update_strain_tensor_and_compute_reaction_force_and_torque();
      virtual void compute_load();
      double get_linear_stiffness() const;
      double get_angular_stiffness() const;
      virtual double compute_critical_time_step() const;
      virtual std::string info() const;
      virtual void draw();
      
 
      //ACCESSORS
      GRANOO_ACCESS(manage_crack_closure      , bool           , _manage_crack_closure       );
      GRANOO_ACCESS(young_modulus             , double         , _young_modulus              );
      GRANOO_ACCESS(poisson_ratio             , double         , _poisson_ratio              );
      GRANOO_ACCESS(friction_coeff            , double         , _friction_coeff             );
      GRANOO_ACCESS(maximum_tensile_stress    , double         , _maximum_tensile_stress     );
      GRANOO_ACCESS(maximum_compressive_stress, double         , _maximum_compressive_stress );
      GRANOO_ACCESS(maximum_shear_stress      , double         , _maximum_shear_stress       );
      GRANOO_ACCESS(current_stress_vector     , Geom::Vector   , _current_stress_vector      );
      GRANOO_ACCESS(current_normal_stress     , double         , _current_normal_stress      );
      GRANOO_ACCESS(current_shear_stress      , double         , _current_shear_stress       );
      GRANOO_ACCESS(current_stress_tensor     , Geom::SymTensor, _current_stress_tensor      );
      GRANOO_ACCESS(current_strain_tensor     , Geom::SymTensor, _current_strain_tensor      );
 
      GRANOO_ACCESS_GET_COPY(current_normal_stress_abs, double  ,fabs( _current_normal_stress)    );
 
    private:
      void fail(const FailureMode mode);
      
    private:
      FlatBond() = delete;      
      FlatBond(const FlatBond&) = delete;
      FlatBond & operator=(const FlatBond &) = delete;
 
    private:
      //BOOST SERIALIZATION
      friend class boost::serialization::access;
      template<class Archive> void serialize(Archive&, const unsigned int);
      
    private:
      VoronoiBond*    _voronoi_bond;
      bool            _manage_crack_closure;
      double          _young_modulus;
      double          _poisson_ratio;
      double          _friction_coeff;
      double          _maximum_tensile_stress;
      double          _maximum_compressive_stress;
      double          _maximum_shear_stress;
      Geom::Point     _FC1;
      Geom::Point     _FC2;
      Geom::Vector    _current_stress_vector;
      double          _current_normal_stress;
      double          _current_shear_stress;
      Geom::SymTensor _current_stress_tensor;
      Geom::SymTensor _current_strain_tensor;
    };    
 
    //BOOST SERIALIZATION 
    template<class Archive> void 
    FlatBond::serialize(Archive & ar, const unsigned int version) {
      ar & boost::serialization::base_object<Bond>(*this);
      ar & _young_modulus;
      ar & _poisson_ratio;
      ar & _maximum_tensile_stress;
      ar & _maximum_compressive_stress;
      ar & _current_stress_vector;
      ar & _current_normal_stress;
      ar & _FC1;
      ar & _FC2;
      if (version >= 1)
    ar & _current_stress_tensor;
      if (version >= 2)
    ar & _current_strain_tensor;
      if (version >= 3)
    ar & _maximum_shear_stress;
      if (version >= 4)
    ar & _manage_crack_closure;
      if (version >= 5)
    ar & _friction_coeff;
      if (version >= 6)
    ar & _current_shear_stress;
      
    }
    
  } 
}
 
 
#include <boost/serialization/version.hpp>
BOOST_CLASS_VERSION(GranOO3::DEM::FlatBond, 6)
 
 
// 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::FlatBond * 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>
    inline void load_construct_data(Archive & ar, 
                    GranOO3::DEM::FlatBond * t, 
                    const unsigned int) {
      GranOO3::DEM::Element* de1 = nullptr;
      GranOO3::DEM::Element* de2 = nullptr;
      ar >> de1;
      ar >> de2;
      ::new(t)GranOO3::DEM::FlatBond(*de1, *de2, false);
    }
  }
} 
 
 
namespace GranOO3{
  GRANOO_CLASS_DECLARE_TPL(DEM::FlatBond);
}
 
 
#endif