DeformPeriodicBox.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_ProcessPlugIn_DeformPeriodicBox_hpp_
#define _libDEM_ProcessPlugIn_DeformPeriodicBox_hpp_
 
#include "GranOO3/Core/PlugIn.hpp"
 
#include "GranOO3/Math/ExpressionArray.hpp"
#include "GranOO3/Math/Pid.hpp"
#include "GranOO3/Math/Variable.hpp"
 
#include "GranOO3/Core/SetOf.hpp"
#include "GranOO3/Core/Domain.hpp"
#include "GranOO3/Core/NeedSetOf.hpp"
#include "GranOO3/Core/Problem.hpp"
#include "GranOO3/Physic/Time.hpp"
#include "GranOO3/Physic/Node.hpp"
#include "GranOO3/Geom/SymTensor.hpp"
#include "GranOO3/DEM/Bond.hpp"
#include "GranOO3/DEM/PeriodicBoundary/Full3D.hpp"
 
namespace GranOO3
{
  namespace PlugIn
  {
    class DeformPeriodicBox : public Core::PlugInInterface<DeformPeriodicBox> {
 
      DECLARE_STD_GRANOO_PLUGIN(DEFORM-PERIODIC-BOX);
 
    public:
      DeformPeriodicBox();
      ~DeformPeriodicBox();
 
      GRANOO_ACCESS(stress, Geom::SymTensor, _stress);
      GRANOO_ACCESS(strain, Geom::SymTensor, _strain);
 
      GRANOO_ACCESS_REF(driven_stress , Math::ExpressionArray<6>, _driven_stress );
      GRANOO_ACCESS_REF(applied_strain, Math::ExpressionArray<6>, _applied_strain);
      
      void init_periodic_bond_set();
      void update_stress_tensor();
      void update_strain_tensor();
      
    public:
      void parse_xml();
      void run();
      void init();
 
    private:
      static void init_math_variable();
      static Geom::SymTensor _stress;
      static Geom::SymTensor _strain;
      static Math::Variable*  _stress_xx_var;
      static Math::Variable*  _stress_yy_var;
      static Math::Variable*  _stress_zz_var;
      static bool             _math_variable_initilized;
      
    private:
      bool _init_sensor;
      Core::SetOf<DEM::Bond>* _xmin;
      Core::SetOf<DEM::Bond>* _xmax;
      Core::SetOf<DEM::Bond>* _ymin;
      Core::SetOf<DEM::Bond>* _ymax;
      Core::SetOf<DEM::Bond>* _zmin;
      Core::SetOf<DEM::Bond>* _zmax;
      DEM::Full3D* _periodic_boundary;
      double _initial_length_x;
      double _initial_length_y;
      double _initial_length_z;
      Math::ExpressionArray<6> _applied_strain;
      Math::ExpressionArray<6> _driven_stress;
      Math::Pid _driven_stress_pid_xx;
      Math::Pid _driven_stress_pid_yy;
      Math::Pid _driven_stress_pid_zz;
      Math::Pid _driven_stress_pid_xy;
      Math::Pid _driven_stress_pid_xz;
      Math::Pid _driven_stress_pid_yz;   
    };
  }
}
 
#endif