ApplyFracture.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 _GranOO_PlugIn_ApplyFracture_hpp_
#define _GranOO_PlugIn_ApplyFracture_hpp_
 
#include <string>
#include <iostream>
#include <fstream>
 
#include "GranOO3/Core/NeedSetOf.hpp"
#include "GranOO3/DEM/Element.hpp"
#include "GranOO3/DEM/Bond.hpp"
#include "GranOO3/DEM/Prop/Stress.hpp"
#include "GranOO3/Core/PlugIn.hpp"
 
#include "GranOO3/Math/Stat.hpp"
#include "GranOO3/Math/Expression.hpp"
 
 
namespace GranOO3
{
  namespace PlugIn
  {
 
 
 
    class ApplyFracture : public Core::PlugInInterface<ApplyFracture>,
              public Core::NeedSetOf< DEM::Element >
    {
 
      DECLARE_STD_GRANOO_PLUGIN(APPLY-FRACTURE);
 
    public:
      enum ApplyCriterionOn {ALL, HALF_PLANE, CLOSEST};
      enum FailureMode {BRITTLE, PLASTIC};
 
    public:
      ApplyFracture();
      void init();
      void run();
 
      void parse_xml();
      void write(std::ofstream&, const std::vector<double>&, const std::string&);
      void apply_failure_criterion();
 
 
      GRANOO_ACCESS(von_mises_criterion    ,  double      , _von_mises_criterion            );
      GRANOO_ACCESS(max_principal_criterion,  double      , _max_principal_stress_criterion );
      GRANOO_ACCESS(mode                   ,  FailureMode , _failure_mode                   );
      GRANOO_ACCESS(broken_bond_number     ,  unsigned int, _broken_bond_number              );
 
    private:
      ~ApplyFracture() {};
 
      void apply_failure_criterion(DEM::Element& de);
      void apply_failure_criterion(Core::SetOfBase< DEM::Bond >&);
      Geom::Vector compute_failure_direction(DEM::Element& de);
 
    private:
      const std::string _crack_set_ID;
      ApplyCriterionOn _apply_on;
      FailureMode _failure_mode;
      bool _add_element_to_crack_set;
      bool _delete_broken_bond;
      double _half_plane_coeff;
      double _hydrostatic_criterion;
      double _von_mises_criterion;
      double _tresca_criterion;
      double _max_principal_stress_criterion;
      double _s1, _s2, _s3;                         // store principal stress value for math::function
      std::vector<Math::Variable*>   _criterion_var;
      std::vector<double*>           _criterion_var_value;
      std::vector<Math::Expression*> _criterion_var_expr;
      std::vector<Math::Expression*> _criterion_law_expr;
      unsigned int _broken_bond_number;
    };
  }
}
#endif