d8/d42/fem__maillage_8h_source.html

//####//------------------------------------------------------------

//####//------------------------------------------------------------

//####//  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