Element.hpp

Go to the documentation of this file.

// This file is part of GranOO, a workbench for FEM 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 _LibFEM_Element
#define _LibFEM_Element
 
#include <string>
#include <iostream>
#include <fstream>
 
#include "GranOO3/Physic/ExtendedSetOf/Node.hpp"
#include "GranOO3/Physic/Node.hpp"
#include "GranOO3/3rdParty/Eigen/Sparse"
#include "GranOO3/3rdParty/Eigen/Dense"
#include "GranOO3/Common.hpp"
 
namespace GranOO3 {
 
  namespace FEM {
 
    class Element : public Core::Base, public Core::Register<Element> {
 
    public:
      GRANOO_CLASS(FEM::Element, Core::Base);
 
      Element();
 
      virtual ~Element();
 
      virtual Eigen::VectorXd shape_function_at_coordinate(double xi, double eta, double zeta) = 0;
      virtual Eigen::MatrixXd shape_function_derivative_at_coordinate(double xi, double eta, double zeta) = 0;
 
      virtual unsigned int node_count() = 0;
      virtual unsigned int dof_per_node_count() = 0;
      virtual unsigned int gauss_point_count() = 0;
      virtual std::string name() = 0;
      virtual std::string vtk_name() = 0;
      virtual Eigen::MatrixXd d_ndx_at_gauss_point(unsigned int gaussPointIndex);
      virtual Eigen::MatrixXd shape_function_at_gauss_point(unsigned int gaussPointIndex);
      virtual Eigen::MatrixXd shape_function_derivative_at_gauss_point(unsigned int gaussPointIndex);
      virtual Eigen::VectorXd strain_at_gauss_point(unsigned int gaussPointIndex);
      virtual Eigen::VectorXd internal_force();
      virtual Eigen::VectorXd internal_force_at_gauss_point(unsigned int gaussPointIndex);
      virtual Eigen::MatrixXd mass_matrix();
      virtual Eigen::MatrixXd mass_matrix_at_gauss_point(unsigned int gaussPointIndex);
 
      void setup_element();
      void reset_variable();
      Eigen::MatrixXd b_matrix_for_d_ndX(Eigen::MatrixXd dNdx);
      void update_state_at_gauss_point(unsigned int gaussPointIndex, const Eigen::VectorXd &dStrain);
      void update_variable();
 
      inline std::vector<GranOO3::Physic::Node *> &node_list() { return _node_list; };
 
      inline const Eigen::MatrixXd &gauss_point() { return _gauss_point; };
      inline const Eigen::MatrixXd &node_coord() { return _node_coord; };
      inline const std::vector<unsigned int> &connectivity() { return _connectivity; };
      inline const Eigen::VectorXd & stress_at_gauss_point(unsigned int igaussPoint) { return _oldStress[igaussPoint]; };
 
    public:
      double _youngModulus;
      double _poissonRatio;
      double _density;
 
    protected:
      Eigen::MatrixXd _gauss_point;
      Eigen::VectorXd _weights;
      std::vector<GranOO3::Physic::Node *> _node_list;
 
      Eigen::MatrixXd _node_coord;
      std::vector<Eigen::MatrixXd> _strainDisplacementMatrices;
      std::vector<double> _jacobians;
      std::vector<Eigen::VectorXd> _oldStrain;
      std::vector<Eigen::VectorXd> _newStrain;
      std::vector<Eigen::VectorXd> _oldStress;
      std::vector<Eigen::VectorXd> _newStress;
      std::vector<Eigen::VectorXd> _oldStateVariables;
      std::vector<Eigen::VectorXd> _newStateVariables;
      std::vector<Eigen::MatrixXd> _ctMatrix;
      std::vector<unsigned int> _connectivity;
    };
  }
 
  GRANOO_CLASS_DECLARE_TPL(FEM::Element);
}
 
#endif