Stress.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_Prop_Stress_HPP
#define _LibDEM_Prop_Stress_HPP
 
 
 
#include "GranOO3/Core/Prop.hpp"
#include "GranOO3/DEM/Element.hpp"
#include "GranOO3/Geom/SymTensor.hpp"
 
 
#define GRANOO_STESS_GET_CRITERION(I)      \
  double get_criterion ## I() const {if (_criterion_list.size() > I) return _criterion_list[I]; return 0.;}
 
 
namespace GranOO3
{
  namespace DEM
  {
 
    class Stress :  public Core::Prop<Element>,
            public Core::Register<Stress>
    {
 
      GRANOO_CLASS (DEM::Stress, Core::Prop<Element>);
 
    public:
      //CONSTRUCTORS & DESTRUCTORS
      Stress();
      virtual ~Stress();
 
      //USEFULL
      virtual std::string info() const;
 
      //ACCESSORS
      GRANOO_ACCESS(neighbour_level , unsigned int   , _neighbour_level  );
      GRANOO_ACCESS(volume_fraction , double         , _volume_fraction  );
      GRANOO_ACCESS(include_velocity, bool           , _include_velocity );
      GRANOO_ACCESS(stress_tensor   , Geom::SymTensor, _stress_tensor    );
      GRANOO_ACCESS(criterion       , double         , _criterion        );
      GRANOO_ACCESS(periodic        , bool           , _periodic         );
      GRANOO_ACCESS_GET(principal_stress_1, double      , _principal_stress_1 );
      GRANOO_ACCESS_GET(principal_stress_2, double      , _principal_stress_2 );
      GRANOO_ACCESS_GET(principal_stress_3, double      , _principal_stress_3 );
      GRANOO_ACCESS_GET(principal_frame   , Geom::Frame , _principal_frame    );
      GRANOO_ACCESS_GET(principal_stress  , Geom::Vector, _principal_stress   );
      
      double get_hydrostatic_stress() const;
      double get_vonmises_stress() const;
      double get_tresca_stress() const;
 
      void add_criterion(double value);
      double& criterion(size_t i);
 
      GRANOO_STESS_GET_CRITERION(0);
      GRANOO_STESS_GET_CRITERION(1);
      GRANOO_STESS_GET_CRITERION(2);
      GRANOO_STESS_GET_CRITERION(3);
      GRANOO_STESS_GET_CRITERION(4);
      GRANOO_STESS_GET_CRITERION(5);
 
      //METHOD
      void compute_virial_stress(); 
      void compute_virial_stress_for(Element* de, unsigned int level, std::set<Element*>& set);
      void compute_lsq_stress();
      void update_principal_stress();
 
    private:
      //Serialization
      friend class boost::serialization::access;
      template<class Archive> void serialize(Archive& ar, const unsigned int );
 
    private:
      Geom::SymTensor  _stress_tensor;
      Geom::Point      _principal_point;
      Geom::Quaternion _principal_quat;
      Geom::Frame      _principal_frame;
      Geom::Vector     _principal_stress;
      double           _principal_stress_1;
      double           _principal_stress_2;
      double           _principal_stress_3;
      // failure criterion
      double              _criterion;
      std::vector<double> _criterion_list;
      // for virial stress computation
      unsigned int _neighbour_level;
      double       _volume_fraction;
      bool         _include_velocity;
      bool         _periodic;
    };
 
    template<class Archive> void
    Stress::serialize(Archive& ar, const unsigned int version) {
      ar & boost::serialization::base_object<Core::Prop<Element> >(*this);
      ar & _neighbour_level;
      ar & _volume_fraction;
      ar & _stress_tensor;
      ar & _principal_stress;
      ar & _principal_stress_1;
      ar & _principal_stress_2;
      ar & _principal_stress_3;
      ar & _criterion;
      ar & _criterion_list;
    }
 
  }
}
 
 
#include <boost/serialization/version.hpp>
BOOST_CLASS_VERSION(GranOO3::DEM::Stress,0)
 
 
namespace GranOO3
{
  GRANOO_CLASS_DECLARE_TPL(DEM::Stress);
}
 
 
#endif