fem_maillage.h

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fem_maillage.h
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:54 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#ifndef __FEMMAILLAGE_
#define __FEMMAILLAGE_
 
 
 
 
 
#include <map>
#include <set>
 
 
 
#include "mg_identificateur.h"
#include "fem_noeud.h"
#include "fem_element0.h"
#include "fem_element1.h"
#include "fem_element2.h"
#include "fem_element3.h"
#include "xfem_element0.h"
#include "xfem_element1.h"
#include "xfem_element2.h"
#include "xfem_element3.h"
#include "ot_boite_3d.h"
#include "tpl_map_entite.h"
 
#include "mg_definition.h"
 
 
 
 
typedef std::set<MG_IDENTIFICATEUR*,inf<MG_IDENTIFICATEUR*> > LISTE_ENTITE;
typedef std::map<unsigned long,FEM_NOEUD*,std::less<unsigned long> > LISTE_FEM_NOEUD;
typedef std::map<unsigned long,FEM_ELEMENT0*,std::less<unsigned long> > LISTE_FEM_ELEMENT0;
typedef std::map<unsigned long,FEM_ELEMENT1*,std::less<unsigned long> > LISTE_FEM_ELEMENT1;
typedef std::map<unsigned long,FEM_ELEMENT2*,std::less<unsigned long> > LISTE_FEM_ELEMENT2;
typedef std::map<unsigned long,FEM_ELEMENT3*,std::less<unsigned long> > LISTE_FEM_ELEMENT3;
typedef std::map<unsigned long,XFEM_ELEMENT0*,std::less<unsigned long> > LISTE_XFEM_ELEMENT0;
typedef std::map<unsigned long,XFEM_ELEMENT1*,std::less<unsigned long> > LISTE_XFEM_ELEMENT1;
typedef std::map<unsigned long,XFEM_ELEMENT2*,std::less<unsigned long> > LISTE_XFEM_ELEMENT2;
typedef std::map<unsigned long,XFEM_ELEMENT3*,std::less<unsigned long> > LISTE_XFEM_ELEMENT3;
 
 std::ostream& operator << (std::ostream& o,class FEM_MAILLAGE& maillage);
 
 
class  FEM_MAILLAGE : public MG_IDENTIFICATEUR
{
public:
    FEM_MAILLAGE(class MG_GEOMETRIE* geo,class MG_MAILLAGE* mai,int deg);
    FEM_MAILLAGE(unsigned long num,MG_GEOMETRIE* geo,MG_MAILLAGE* mai,int deg);
    virtual ~FEM_MAILLAGE();
    MG_GEOMETRIE* get_mg_geometrie(void);
    MG_MAILLAGE* get_mg_maillage(void);
    void change_gestionnaire(class MG_GESTIONNAIRE* mggest);
    int get_degre(void);
    void adapte_degre_importation_maillage(int deg);
 
 
    FEM_NOEUD* ajouter_fem_noeud(MG_NOEUD* mgnoeud,unsigned long num=0);
    int ajouter_fem_noeud(FEM_NOEUD *femnoeud);
    FEM_NOEUD* get_fem_noeudid(unsigned long num);
    FEM_NOEUD* get_fem_noeud(unsigned int num);
    FEM_NOEUD* get_premier_noeud(LISTE_FEM_NOEUD::iterator &  it);
    FEM_NOEUD* get_suivant_noeud(LISTE_FEM_NOEUD::iterator &  it);
    unsigned int get_nb_fem_noeud(void);
    int supprimer_fem_noeudid(unsigned long num);
    int supprimer_fem_noeud(unsigned int num);
    void supprimer_tout_fem_noeud(void);
 
 
 
    int ajouter_fem_element0(FEM_ELEMENT0 *femnoeud);
    FEM_ELEMENT0* get_fem_element0id(unsigned long num);
    FEM_ELEMENT0* get_fem_element0(unsigned int num);
    FEM_ELEMENT0* get_premier_element0(LISTE_FEM_ELEMENT0::iterator &  it);
    FEM_ELEMENT0* get_suivant_element0(LISTE_FEM_ELEMENT0::iterator &  it);
    unsigned int get_nb_fem_element0(void);
    int supprimer_fem_element0id(unsigned long num);
    int supprimer_fem_element0(unsigned int num);
    void supprimer_tout_fem_element0(void);
    
    int ajouter_fem_element1(FEM_ELEMENT1 *femsegment);
    FEM_ELEMENT1* get_fem_element1id(unsigned long num);
    FEM_ELEMENT1* get_fem_element1(unsigned int num);
    FEM_ELEMENT1* get_premier_element1(LISTE_FEM_ELEMENT1::iterator &  it);
    FEM_ELEMENT1* get_suivant_element1(LISTE_FEM_ELEMENT1::iterator &  it);
    unsigned int get_nb_fem_element1(void);
    int supprimer_fem_element1id(unsigned long num);
    int supprimer_fem_element1(unsigned int num);
    void supprimer_tout_fem_element1(void);
 
 
 
    int ajouter_fem_element2(FEM_ELEMENT2 *femtriangle);
    FEM_ELEMENT2* get_fem_element2id(unsigned long num);
    FEM_ELEMENT2* get_fem_element2(unsigned int num);
    FEM_ELEMENT2* get_premier_element2(LISTE_FEM_ELEMENT2::iterator &  it) ;
    FEM_ELEMENT2* get_suivant_element2(LISTE_FEM_ELEMENT2::iterator &  it) ;
    unsigned int get_nb_fem_element2(void);
    int supprimer_fem_element2id(unsigned long num);
    int supprimer_fem_element2(unsigned int num);
    void supprimer_tout_fem_element2(void);
 
    int ajouter_fem_element3(FEM_ELEMENT3 *femelement3);
    FEM_ELEMENT3* get_fem_element3id(unsigned long num);
    FEM_ELEMENT3* get_fem_element3(unsigned int num);
    FEM_ELEMENT3* get_premier_element3(LISTE_FEM_ELEMENT3::iterator &  it);
    FEM_ELEMENT3* get_suivant_element3(LISTE_FEM_ELEMENT3::iterator &  it);
    unsigned int get_nb_fem_element3(void);
    int supprimer_fem_element3id(unsigned long num);
    int supprimer_fem_element3(unsigned int num);
    void supprimer_tout_fem_element3(void);
    
    int ajouter_xfem_element0(XFEM_ELEMENT0 *femelement0);
    XFEM_ELEMENT0* get_xfem_element0id(unsigned long num);
    XFEM_ELEMENT0* get_xfem_element0(unsigned int num);
    XFEM_ELEMENT0* get_premier_xelement0(LISTE_XFEM_ELEMENT0::iterator &  it);
    XFEM_ELEMENT0* get_suivant_xelement0(LISTE_XFEM_ELEMENT0::iterator &  it);
    unsigned int get_nb_xfem_element0(void);
    int supprimer_xfem_element0id(unsigned long num);
    int supprimer_xfem_element0(unsigned int num);
    void supprimer_tout_xfem_element0(void);
    
    int ajouter_xfem_element1(XFEM_ELEMENT1 *femelement1);
    XFEM_ELEMENT1* get_xfem_element1id(unsigned long num);
    XFEM_ELEMENT1* get_xfem_element1(unsigned int num);
    XFEM_ELEMENT1* get_premier_xelement1(LISTE_XFEM_ELEMENT1::iterator &  it);
    XFEM_ELEMENT1* get_suivant_xelement1(LISTE_XFEM_ELEMENT1::iterator &  it);
    unsigned int get_nb_xfem_element1(void);
    int supprimer_xfem_element1id(unsigned long num);
    int supprimer_xfem_element1(unsigned int num);
    void supprimer_tout_xfem_element1(void);
    
    int ajouter_xfem_element2(XFEM_ELEMENT2 *femelement2);
    XFEM_ELEMENT2* get_xfem_element2id(unsigned long num);
    XFEM_ELEMENT2* get_xfem_element2(unsigned int num);
    XFEM_ELEMENT2* get_premier_xelement2(LISTE_XFEM_ELEMENT2::iterator &  it);
    XFEM_ELEMENT2* get_suivant_xelement2(LISTE_XFEM_ELEMENT2::iterator &  it);
    unsigned int get_nb_xfem_element2(void);
    int supprimer_xfem_element2id(unsigned long num);
    int supprimer_xfem_element2(unsigned int num);
    void supprimer_tout_xfem_element2(void);
    
    int ajouter_xfem_element3(XFEM_ELEMENT3 *femelement3);
    XFEM_ELEMENT3* get_xfem_element3id(unsigned long num);
    XFEM_ELEMENT3* get_xfem_element3(unsigned int num);
    XFEM_ELEMENT3* get_premier_xelement3(LISTE_XFEM_ELEMENT3::iterator &  it);
    XFEM_ELEMENT3* get_suivant_xelement3(LISTE_XFEM_ELEMENT3::iterator &  it);
    unsigned int get_nb_xfem_element3(void);
    int supprimer_xfem_element3id(unsigned long num);
    int supprimer_xfem_element3(unsigned int num);
    void supprimer_tout_xfem_element3(void);
    
    
    //void construire(int num=0);
    void exporter_cosmos(std::string& nomfic);
    //void optimise_numerotation(void);
    //void copie_numerotation_opt(void);
    bool contient ( MG_IDENTIFICATEUR * id );
 
    void recherche_connexion_multidimension(TPL_MAP_ENTITE<class MG_SEGMENT*> &mini);
    
    
    void calcul_deforme(class FEM_SOLUTION* solx,int numx,class FEM_SOLUTION* soly,int numy,class FEM_SOLUTION* solz,int numz);
    int existe_deforme(void);
    void active_deforme(void);
    void desactive_deforme(void);
    int existe_reaction(void);
    void active_reaction(void);
    void desactive_reaction(void);
    double calcul_coef(void);
    void calcul_somme_reaction(double &rx,double &ry,double &rz);
    virtual void enregistrer(std::ostream& o,double version);
    virtual void get_fichier_dependant(std::vector<std::string> &liste_fichier);
    FEM_MAILLAGE* dupliquer(MG_GESTIONNAIRE* gest);
    void decompte(int &nb_noeud,int &nb_seg2,int &nb_seg3,int &nb_tri3,int &nb_tri6,int &nb_quad4,int &nb_quad8,int &nb_tet4,int &nb_tet10,int &nb_hex8,int &nb_hex20);
    friend std::ostream& operator << (std::ostream& o,FEM_MAILLAGE& maillage);
    void change_param_structure(BOITE_3D btmp,int nxtmp,int nytmp,int nztmp);
    bool get_param_structure(BOITE_3D &btmp,int &nxtmp,int &nytmp,int &nztmp);
    bool get_est_structure(void);
 
private:
    //void construire_lineaire(TPL_MAP_ENTITE<MG_SEGMENT*> &mini);
    //void construire_quadratique(void);
    LISTE_ENTITE lst_entite;
    LISTE_FEM_NOEUD lst_fem_noeud;
    LISTE_FEM_ELEMENT0 lst_fem_element0;
    LISTE_FEM_ELEMENT1 lst_fem_element1;
    LISTE_FEM_ELEMENT2 lst_fem_element2;
    LISTE_FEM_ELEMENT3 lst_fem_element3;
    LISTE_XFEM_ELEMENT0 lst_xfem_element0;
    LISTE_XFEM_ELEMENT1 lst_xfem_element1;
    LISTE_XFEM_ELEMENT2 lst_xfem_element2;
    LISTE_XFEM_ELEMENT3 lst_xfem_element3;
    MG_GEOMETRIE* mg_geometrie;
    MG_GESTIONNAIRE *gest;
    MG_MAILLAGE* maillage;
    int degre;
    int deforme;
    int reaction;
    
    
    BOITE_3D boite_englobante;
    int nx,ny,nz;
    bool est_structure;
  
};
 
 
 
#endif