ElementT4.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 _LibFEM_ElementT4
#define _LibFEM_ElementT4
 
#include <string>
#include <iostream>
#include <fstream>
 
#include "GranOO3/FEM/ElementFactory.hpp"
#include "GranOO3/FEM/Element.hpp"
 
namespace GranOO3 {
 
  namespace FEM {
 
    class ElementT4 : public Element, public Core::Register<ElementT4> {
 
    public:
      GRANOO_CLASS(FEM::ElementT4, FEM::Element);
 
      ElementT4();
      ElementT4(std::vector<unsigned int> connectivity);
      virtual ~ElementT4();
 
      Eigen::VectorXd shape_function_at_coordinate(double xi, double eta, double zeta);
      Eigen::MatrixXd shape_function_derivative_at_coordinate(double xi, double eta, double zeta);
 
      static const unsigned int node_number = 4;
      inline unsigned int node_count() { return node_number; };
      inline unsigned int dof_per_node_count() { return 3; };
      inline unsigned int gauss_point_count() { return 1; };
      inline static std::string element_name() { return "T4"; };
      inline std::string name() { return "T4"; };
      inline std::string vtk_name() { return "10"; };
 
      // declare the element to the factory
      DECLARE_ELEMENT(ElementT4)
    };
  }
 
  GRANOO_CLASS_DECLARE_TPL(FEM::ElementT4);
}
 
#endif