ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/geometrie/src/fem_maillage.h
Revision: 399
Committed: Mon May 27 15:27:41 2013 UTC (11 years, 11 months ago) by francois
Content type: text/plain
File size: 9427 byte(s)
Log Message:
Creation d'une structure XFEM et FEM complete avec des element0 et des xelement0

File Contents

# User Rev Content
1 francois 283 //------------------------------------------------------------
2     //------------------------------------------------------------
3     // MAGiC
4     // Jean Christophe Cuilli�re et Vincent FRANCOIS
5     // D�partement de G�nie M�canique - UQTR
6     //------------------------------------------------------------
7     // Le projet MAGIC est un projet de recherche du d�partement
8     // de g�nie m�canique de l'Universit� du Qu�bec �
9     // Trois Rivi�res
10     // Les librairies ne peuvent �tre utilis�es sans l'accord
11     // des auteurs (contact : francois@uqtr.ca)
12     //------------------------------------------------------------
13     //------------------------------------------------------------
14     //
15     // fem_maillage.h
16     //
17     //------------------------------------------------------------
18     //------------------------------------------------------------
19     // COPYRIGHT 2000
20     // Version du 02/03/2006 � 11H22
21     //------------------------------------------------------------
22     //------------------------------------------------------------
23     #ifndef __FEMMAILLAGE_
24     #define __FEMMAILLAGE_
25    
26    
27     #ifdef WINDOWS_VERSION
28     #ifdef BUILT_DLL_GEOMETRIE
29     #define DLLPORTGEOMETRIE __declspec(dllexport)
30     #else
31     #define DLLPORTGEOMETRIE __declspec(dllimport)
32     #endif
33     #else
34     #define DLLPORTGEOMETRIE
35     #endif
36    
37     #ifdef BORLANDCPP
38     #include <map.h>
39     #include <set.h>
40     #else
41     #include <map>
42     #include <set>
43     #endif
44    
45    
46     #include "mg_identificateur.h"
47     #include "fem_noeud.h"
48 francois 399 #include "fem_element0.h"
49 francois 310 #include "fem_element1.h"
50     #include "fem_element2.h"
51     #include "fem_element3.h"
52 francois 399 #include "xfem_element0.h"
53 francois 380 #include "xfem_element1.h"
54     #include "xfem_element2.h"
55 francois 339 #include "xfem_element3.h"
56 francois 310
57 francois 375 #include "mg_definition.h"
58 francois 283
59    
60    
61    
62     typedef std::set<MG_IDENTIFICATEUR*,inf<MG_IDENTIFICATEUR*> > LISTE_ENTITE;
63     typedef std::map<unsigned long,FEM_NOEUD*,std::less<unsigned long> > LISTE_FEM_NOEUD;
64 francois 399 typedef std::map<unsigned long,FEM_ELEMENT0*,std::less<unsigned long> > LISTE_FEM_ELEMENT0;
65 francois 309 typedef std::map<unsigned long,FEM_ELEMENT1*,std::less<unsigned long> > LISTE_FEM_ELEMENT1;
66     typedef std::map<unsigned long,FEM_ELEMENT2*,std::less<unsigned long> > LISTE_FEM_ELEMENT2;
67     typedef std::map<unsigned long,FEM_ELEMENT3*,std::less<unsigned long> > LISTE_FEM_ELEMENT3;
68 francois 399 typedef std::map<unsigned long,XFEM_ELEMENT0*,std::less<unsigned long> > LISTE_XFEM_ELEMENT0;
69 francois 380 typedef std::map<unsigned long,XFEM_ELEMENT1*,std::less<unsigned long> > LISTE_XFEM_ELEMENT1;
70     typedef std::map<unsigned long,XFEM_ELEMENT2*,std::less<unsigned long> > LISTE_XFEM_ELEMENT2;
71 francois 339 typedef std::map<unsigned long,XFEM_ELEMENT3*,std::less<unsigned long> > LISTE_XFEM_ELEMENT3;
72 francois 283
73     DLLPORTGEOMETRIE std::ostream& operator << (std::ostream& o,class FEM_MAILLAGE& maillage);
74    
75    
76     class DLLPORTGEOMETRIE FEM_MAILLAGE : public MG_IDENTIFICATEUR
77     {
78     public:
79     FEM_MAILLAGE(class MG_GEOMETRIE* geo,class MG_MAILLAGE* mai,int deg);
80     FEM_MAILLAGE(unsigned long num,MG_GEOMETRIE* geo,MG_MAILLAGE* mai,int deg);
81     virtual ~FEM_MAILLAGE();
82     MG_GEOMETRIE* get_mg_geometrie(void);
83     MG_MAILLAGE* get_mg_maillage(void);
84     void change_gestionnaire(class MG_GESTIONNAIRE* mggest);
85     int get_degre(void);
86    
87    
88     FEM_NOEUD* ajouter_fem_noeud(MG_NOEUD* mgnoeud,unsigned long num=0);
89     int ajouter_fem_noeud(FEM_NOEUD *femnoeud);
90     FEM_NOEUD* get_fem_noeudid(unsigned long num);
91     FEM_NOEUD* get_fem_noeud(unsigned int num);
92     FEM_NOEUD* get_premier_noeud(LISTE_FEM_NOEUD::iterator & it);
93     FEM_NOEUD* get_suivant_noeud(LISTE_FEM_NOEUD::iterator & it);
94     unsigned int get_nb_fem_noeud(void);
95     int supprimer_fem_noeudid(unsigned long num);
96     int supprimer_fem_noeud(unsigned int num);
97     void supprimer_tout_fem_noeud(void);
98    
99    
100    
101 francois 399 int ajouter_fem_element0(FEM_ELEMENT0 *femnoeud);
102     FEM_ELEMENT0* get_fem_element0id(unsigned long num);
103     FEM_ELEMENT0* get_fem_element0(unsigned int num);
104     FEM_ELEMENT0* get_premier_element0(LISTE_FEM_ELEMENT0::iterator & it);
105     FEM_ELEMENT0* get_suivant_element0(LISTE_FEM_ELEMENT0::iterator & it);
106     unsigned int get_nb_fem_element0(void);
107     int supprimer_fem_element0id(unsigned long num);
108     int supprimer_fem_element0(unsigned int num);
109     void supprimer_tout_fem_element0(void);
110    
111 francois 309 int ajouter_fem_element1(FEM_ELEMENT1 *femsegment);
112     FEM_ELEMENT1* get_fem_element1id(unsigned long num);
113     FEM_ELEMENT1* get_fem_element1(unsigned int num);
114     FEM_ELEMENT1* get_premier_element1(LISTE_FEM_ELEMENT1::iterator & it);
115     FEM_ELEMENT1* get_suivant_element1(LISTE_FEM_ELEMENT1::iterator & it);
116     unsigned int get_nb_fem_element1(void);
117     int supprimer_fem_element1id(unsigned long num);
118     int supprimer_fem_element1(unsigned int num);
119     void supprimer_tout_fem_element1(void);
120 francois 283
121    
122    
123 francois 309 int ajouter_fem_element2(FEM_ELEMENT2 *femtriangle);
124     FEM_ELEMENT2* get_fem_element2id(unsigned long num);
125     FEM_ELEMENT2* get_fem_element2(unsigned int num);
126     FEM_ELEMENT2* get_premier_element2(LISTE_FEM_ELEMENT2::iterator & it) ;
127     FEM_ELEMENT2* get_suivant_element2(LISTE_FEM_ELEMENT2::iterator & it) ;
128     unsigned int get_nb_fem_element2(void);
129     int supprimer_fem_element2id(unsigned long num);
130     int supprimer_fem_element2(unsigned int num);
131     void supprimer_tout_fem_element2(void);
132 francois 283
133 francois 309 int ajouter_fem_element3(FEM_ELEMENT3 *femelement3);
134     FEM_ELEMENT3* get_fem_element3id(unsigned long num);
135     FEM_ELEMENT3* get_fem_element3(unsigned int num);
136     FEM_ELEMENT3* get_premier_element3(LISTE_FEM_ELEMENT3::iterator & it);
137     FEM_ELEMENT3* get_suivant_element3(LISTE_FEM_ELEMENT3::iterator & it);
138     unsigned int get_nb_fem_element3(void);
139     int supprimer_fem_element3id(unsigned long num);
140     int supprimer_fem_element3(unsigned int num);
141     void supprimer_tout_fem_element3(void);
142 francois 399
143     int ajouter_xfem_element0(XFEM_ELEMENT0 *femelement0);
144     XFEM_ELEMENT0* get_xfem_element0id(unsigned long num);
145     XFEM_ELEMENT0* get_xfem_element0(unsigned int num);
146     XFEM_ELEMENT0* get_premier_xelement0(LISTE_XFEM_ELEMENT0::iterator & it);
147     XFEM_ELEMENT0* get_suivant_xelement0(LISTE_XFEM_ELEMENT0::iterator & it);
148     unsigned int get_nb_xfem_element0(void);
149     int supprimer_xfem_element0id(unsigned long num);
150     int supprimer_xfem_element0(unsigned int num);
151     void supprimer_tout_xfem_element0(void);
152    
153 francois 380 int ajouter_xfem_element1(XFEM_ELEMENT1 *femelement1);
154     XFEM_ELEMENT1* get_xfem_element1id(unsigned long num);
155     XFEM_ELEMENT1* get_xfem_element1(unsigned int num);
156     XFEM_ELEMENT1* get_premier_xelement1(LISTE_XFEM_ELEMENT1::iterator & it);
157     XFEM_ELEMENT1* get_suivant_xelement1(LISTE_XFEM_ELEMENT1::iterator & it);
158     unsigned int get_nb_xfem_element1(void);
159     int supprimer_xfem_element1id(unsigned long num);
160     int supprimer_xfem_element1(unsigned int num);
161     void supprimer_tout_xfem_element1(void);
162    
163     int ajouter_xfem_element2(XFEM_ELEMENT2 *femelement2);
164     XFEM_ELEMENT2* get_xfem_element2id(unsigned long num);
165     XFEM_ELEMENT2* get_xfem_element2(unsigned int num);
166     XFEM_ELEMENT2* get_premier_xelement2(LISTE_XFEM_ELEMENT2::iterator & it);
167     XFEM_ELEMENT2* get_suivant_xelement2(LISTE_XFEM_ELEMENT2::iterator & it);
168     unsigned int get_nb_xfem_element2(void);
169     int supprimer_xfem_element2id(unsigned long num);
170     int supprimer_xfem_element2(unsigned int num);
171     void supprimer_tout_xfem_element2(void);
172    
173 francois 339 int ajouter_xfem_element3(XFEM_ELEMENT3 *femelement3);
174     XFEM_ELEMENT3* get_xfem_element3id(unsigned long num);
175     XFEM_ELEMENT3* get_xfem_element3(unsigned int num);
176     XFEM_ELEMENT3* get_premier_xelement3(LISTE_XFEM_ELEMENT3::iterator & it);
177     XFEM_ELEMENT3* get_suivant_xelement3(LISTE_XFEM_ELEMENT3::iterator & it);
178     unsigned int get_nb_xfem_element3(void);
179     int supprimer_xfem_element3id(unsigned long num);
180     int supprimer_xfem_element3(unsigned int num);
181     void supprimer_tout_xfem_element3(void);
182 francois 380
183    
184 francois 283 void construire(int num=0);
185     void exporter_cosmos(std::string& nomfic);
186     void optimise_numerotation(void);
187     void echange_numerotation(void);
188     bool contient ( MG_IDENTIFICATEUR * id );
189    
190 francois 339
191    
192    
193 francois 283 void calcul_deforme(class FEM_SOLUTION* solx,int numx,class FEM_SOLUTION* soly,int numy,class FEM_SOLUTION* solz,int numz);
194     int existe_deforme(void);
195     void active_deforme(void);
196 gervaislavoie 304 void desactive_deforme(void);
197 francois 339 int existe_reaction(void);
198     void active_reaction(void);
199     void desactive_reaction(void);
200 francois 283 double calcul_coef(void);
201 francois 325 void calcul_somme_reaction(double &rx,double &ry,double &rz);
202 francois 283 virtual void enregistrer(std::ostream& o);
203     FEM_MAILLAGE* dupliquer(MG_GESTIONNAIRE* gest);
204 francois 310 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);
205 francois 283 friend std::ostream& operator << (std::ostream& o,FEM_MAILLAGE& maillage);
206 francois 339
207 francois 283 private:
208 francois 310 void construire_lineaire(void);
209 francois 283 void construire_quadratique(void);
210     LISTE_ENTITE lst_entite;
211     LISTE_FEM_NOEUD lst_fem_noeud;
212 francois 399 LISTE_FEM_ELEMENT0 lst_fem_element0;
213 francois 309 LISTE_FEM_ELEMENT1 lst_fem_element1;
214     LISTE_FEM_ELEMENT2 lst_fem_element2;
215     LISTE_FEM_ELEMENT3 lst_fem_element3;
216 francois 399 LISTE_XFEM_ELEMENT0 lst_xfem_element0;
217 francois 380 LISTE_XFEM_ELEMENT1 lst_xfem_element1;
218     LISTE_XFEM_ELEMENT2 lst_xfem_element2;
219 francois 339 LISTE_XFEM_ELEMENT3 lst_xfem_element3;
220 francois 283 MG_GEOMETRIE* mg_geometrie;
221     MG_GESTIONNAIRE *gest;
222     MG_MAILLAGE* maillage;
223     int degre;
224     int deforme;
225 francois 339 int reaction;
226    
227 francois 283 };
228    
229    
230    
231     #endif
232    
233