ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/geometrie/src/fem_quadrangle8.cpp
Revision: 674
Committed: Tue Jun 9 18:30:54 2015 UTC (9 years, 11 months ago) by francois
File size: 14606 byte(s)
Log Message: 
Transformation d'un support de carte de taille non structure en structure (formule à preciser et a revoir)
Calcul du temps CPU pour toutes les applications en incluant le temps CPU de code aster
Changement dans la projection et la comparaison de solution pour la generaliser à tous les maillages
Divers bug petit sur maillage hex/quadra

File Contents

# User Rev Content
1 francois 310 //------------------------------------------------------------
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_triangle3.cpp
16     //
17     //------------------------------------------------------------
18     //------------------------------------------------------------
19     // COPYRIGHT 2000
20     // Version du 02/03/2006 � 11H22
21     //------------------------------------------------------------
22     //------------------------------------------------------------
23    
24    
25     #include "gestionversion.h"
26     #include "fem_quadrangle8.h"
27     #include "fem_maillage.h"
28     #include "fem_noeud.h"
29     #include "mg_element_maillage.h"
30    
31    
32     FEM_QUADRANGLE8::FEM_QUADRANGLE8(unsigned long num,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,mai),FEM_TEMPLATE_ELEMENT<8>(tabnoeud)
33     {
34     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
35     tab[0]->get_lien_element2()->ajouter(this);
36     tab[1]->get_lien_element2()->ajouter(this);
37     tab[2]->get_lien_element2()->ajouter(this);
38     tab[3]->get_lien_element2()->ajouter(this);
39     tab[4]->get_lien_element2()->ajouter(this);
40     tab[5]->get_lien_element2()->ajouter(this);
41     tab[6]->get_lien_element2()->ajouter(this);
42     tab[7]->get_lien_element2()->ajouter(this);
43     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
44     }
45    
46     FEM_QUADRANGLE8::FEM_QUADRANGLE8(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(mai),FEM_TEMPLATE_ELEMENT<8>(tabnoeud)
47     {
48     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
49     tab[0]->get_lien_element2()->ajouter(this);
50     tab[1]->get_lien_element2()->ajouter(this);
51     tab[2]->get_lien_element2()->ajouter(this);
52     tab[3]->get_lien_element2()->ajouter(this);
53     tab[4]->get_lien_element2()->ajouter(this);
54     tab[5]->get_lien_element2()->ajouter(this);
55     tab[6]->get_lien_element2()->ajouter(this);
56     tab[7]->get_lien_element2()->ajouter(this);
57     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
58     }
59 francois 378 FEM_QUADRANGLE8::FEM_QUADRANGLE8(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,topo),FEM_TEMPLATE_ELEMENT<8>(tabnoeud)
60     {
61     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
62     tab[0]->get_lien_element2()->ajouter(this);
63     tab[1]->get_lien_element2()->ajouter(this);
64     tab[2]->get_lien_element2()->ajouter(this);
65     tab[3]->get_lien_element2()->ajouter(this);
66     tab[4]->get_lien_element2()->ajouter(this);
67     tab[5]->get_lien_element2()->ajouter(this);
68     tab[6]->get_lien_element2()->ajouter(this);
69     tab[7]->get_lien_element2()->ajouter(this);
70     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
71     }
72 francois 310
73 francois 378 FEM_QUADRANGLE8::FEM_QUADRANGLE8(class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(topo),FEM_TEMPLATE_ELEMENT<8>(tabnoeud)
74     {
75     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
76     tab[0]->get_lien_element2()->ajouter(this);
77     tab[1]->get_lien_element2()->ajouter(this);
78     tab[2]->get_lien_element2()->ajouter(this);
79     tab[3]->get_lien_element2()->ajouter(this);
80     tab[4]->get_lien_element2()->ajouter(this);
81     tab[5]->get_lien_element2()->ajouter(this);
82     tab[6]->get_lien_element2()->ajouter(this);
83     tab[7]->get_lien_element2()->ajouter(this);
84     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
85     }
86     FEM_QUADRANGLE8::FEM_QUADRANGLE8(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,topo,mai),FEM_TEMPLATE_ELEMENT<8>(tabnoeud)
87     {
88     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
89     tab[0]->get_lien_element2()->ajouter(this);
90     tab[1]->get_lien_element2()->ajouter(this);
91     tab[2]->get_lien_element2()->ajouter(this);
92     tab[3]->get_lien_element2()->ajouter(this);
93     tab[4]->get_lien_element2()->ajouter(this);
94     tab[5]->get_lien_element2()->ajouter(this);
95     tab[6]->get_lien_element2()->ajouter(this);
96     tab[7]->get_lien_element2()->ajouter(this);
97     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
98     }
99    
100     FEM_QUADRANGLE8::FEM_QUADRANGLE8(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(topo,mai),FEM_TEMPLATE_ELEMENT<8>(tabnoeud)
101     {
102     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
103     tab[0]->get_lien_element2()->ajouter(this);
104     tab[1]->get_lien_element2()->ajouter(this);
105     tab[2]->get_lien_element2()->ajouter(this);
106     tab[3]->get_lien_element2()->ajouter(this);
107     tab[4]->get_lien_element2()->ajouter(this);
108     tab[5]->get_lien_element2()->ajouter(this);
109     tab[6]->get_lien_element2()->ajouter(this);
110     tab[7]->get_lien_element2()->ajouter(this);
111     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
112     }
113    
114    
115    
116    
117    
118    
119    
120    
121    
122 francois 310 FEM_QUADRANGLE8::FEM_QUADRANGLE8(FEM_QUADRANGLE8& mdd):FEM_ELEMENT2(mdd),FEM_TEMPLATE_ELEMENT<8>(mdd)
123     {
124     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
125     tab[0]->get_lien_element2()->ajouter(this);
126     tab[1]->get_lien_element2()->ajouter(this);
127     tab[2]->get_lien_element2()->ajouter(this);
128     tab[3]->get_lien_element2()->ajouter(this);
129     tab[4]->get_lien_element2()->ajouter(this);
130     tab[5]->get_lien_element2()->ajouter(this);
131     tab[6]->get_lien_element2()->ajouter(this);
132     tab[7]->get_lien_element2()->ajouter(this);
133     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
134     }
135     FEM_QUADRANGLE8::~FEM_QUADRANGLE8()
136     {
137 francois 663 if (liaison_topologique!=NULL)
138     if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
139 francois 310 tab[0]->get_lien_element2()->supprimer(this);
140     tab[1]->get_lien_element2()->supprimer(this);
141     tab[2]->get_lien_element2()->supprimer(this);
142     tab[3]->get_lien_element2()->supprimer(this);
143     tab[4]->get_lien_element2()->supprimer(this);
144     tab[5]->get_lien_element2()->supprimer(this);
145     tab[6]->get_lien_element2()->supprimer(this);
146     tab[7]->get_lien_element2()->supprimer(this);
147     get_fem_noeudpetitid()->get_lien_petit_element2()->supprimer(this);
148     }
149    
150    
151     FEM_ELEMENT_MAILLAGE* FEM_QUADRANGLE8::dupliquer(FEM_MAILLAGE *femmai,long decalage)
152     {
153     FEM_NOEUD* tabnoeud[8];
154     tabnoeud[0]=femmai->get_fem_noeudid(tab[0]->get_id()+decalage);
155     tabnoeud[1]=femmai->get_fem_noeudid(tab[1]->get_id()+decalage);
156     tabnoeud[2]=femmai->get_fem_noeudid(tab[2]->get_id()+decalage);
157     tabnoeud[3]=femmai->get_fem_noeudid(tab[3]->get_id()+decalage);
158     tabnoeud[4]=femmai->get_fem_noeudid(tab[4]->get_id()+decalage);
159     tabnoeud[5]=femmai->get_fem_noeudid(tab[5]->get_id()+decalage);
160     tabnoeud[6]=femmai->get_fem_noeudid(tab[6]->get_id()+decalage);
161     tabnoeud[7]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
162     FEM_QUADRANGLE8* quad=new FEM_QUADRANGLE8(get_id()+decalage,maillage,tabnoeud);
163     femmai->ajouter_fem_element2(quad);
164     return quad;
165     }
166    
167 francois 635 bool FEM_QUADRANGLE8::valide_parametre_EF(double* uvw)
168     {
169 francois 674 if (uvw[0]>=-1.-1e-10)
170     if (uvw[1]>=-1.-1e-10)
171 francois 635 if (uvw[0]<=1.+1e-10)
172     if (uvw[1]<=1.+1e-10)
173     return true;
174     return false;
175     }
176 francois 310
177    
178    
179     int FEM_QUADRANGLE8::get_type_entite(void)
180     {
181     return IDFEM_QUADRANGLE8;
182     }
183    
184     int FEM_QUADRANGLE8::get_dimension(void)
185     {
186     return 2;
187     }
188    
189    
190    
191     int FEM_QUADRANGLE8::get_nb_fem_noeud(void)
192     {
193     return FEM_TEMPLATE_ELEMENT<8>::get_nb_fem_noeud();
194     }
195    
196     FEM_NOEUD* FEM_QUADRANGLE8::get_fem_noeud(int num)
197     {
198     return FEM_TEMPLATE_ELEMENT<8>::get_fem_noeud(num);
199     }
200    
201     void FEM_QUADRANGLE8::change_noeud(int num,FEM_NOEUD* noeud)
202     {
203     FEM_TEMPLATE_ELEMENT<8>::change_noeud(num,noeud);
204     }
205    
206     BOITE_3D& FEM_QUADRANGLE8::get_boite_3D(void)
207     {
208     return FEM_TEMPLATE_ELEMENT<8>::get_boite_3D();
209     }
210    
211    
212     void FEM_QUADRANGLE8::enregistrer(std::ostream& o)
213     {
214 francois 378 if (maillage!=NULL)
215     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_QUADRANGLE8($"<< get_lien_topologie()->get_id() << ",$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ",$" << tab[3]->get_id()<< ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id()<< ",$" << tab[7]->get_id()<< ");" << std::endl;
216     else o << "%" << get_id() << "=FEM_QUADRANGLE8(NULL,$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ",$" << tab[3]->get_id()<< ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id()<< ",$" << tab[7]->get_id()<<");" << std::endl;
217     else
218     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_QUADRANGLE8($"<< get_lien_topologie()->get_id() << ",NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ",$" << tab[3]->get_id()<< ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id()<< ",$" << tab[7]->get_id()<< ");" << std::endl;
219     else o << "%" << get_id() << "=FEM_QUADRANGLE8(NULL,NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ",$" << tab[3]->get_id()<< ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id()<< ",$" << tab[7]->get_id()<<");" << std::endl;
220    
221 francois 310 }
222 francois 637
223     int FEM_QUADRANGLE8::get_nb_pt_gauss(int degre)
224     {
225     if (degre<3) return 3;
226     if (degre<4) return 4;
227     if (degre<6) return 7;
228     return 0;
229     }
230    
231     void FEM_QUADRANGLE8::get_pt_gauss(int degre,int num,double &w,double *uv)
232     {
233     if (degre<3)
234     {
235     if (num==0) {w=1.333333333333333;uv[0]=0.816496580927726;uv[1]=0.;}
236     if (num==1) {w=1.333333333333333;uv[0]=-0.408248290463863;uv[1]=-0.707106781186548;}
237     if (num==2) {w=1.333333333333333;uv[0]=-0.408248290463863;uv[1]=0.707106781186548;}
238     return;
239     }
240     if (degre<4)
241     {
242     if (num==0) {w=1.;uv[0]=0.816496580927726;uv[1]=0.;}
243     if (num==1) {w=1.;uv[0]=-0.816496580927726;uv[1]=0.;}
244     if (num==2) {w=1.;uv[0]=0.;uv[1]=0.816496580927726;}
245     if (num==3) {w=1.;uv[0]=0.;uv[1]=-0.816496580927726;}
246     return;
247     }
248     if (degre<6)
249     {
250     if (num==0) {w=1.142857142857140;uv[0]=0.;uv[1]=0.;}
251     if (num==1) {w=0.317460317460317;uv[0]=0.;uv[1]=0.966091783079296;}
252     if (num==2) {w=0.317460317460317;uv[0]=0.;uv[1]=-0.966091783079296;}
253     if (num==3) {w=0.555555555555556;uv[0]=0.774596669241483;uv[1]=-0.577350269189626;}
254     if (num==4) {w=0.555555555555556;uv[0]=-0.774596669241483;uv[1]=-0.577350269189626;}
255     if (num==5) {w=0.555555555555556;uv[0]=-0.774596669241483;uv[1]=0.577350269189626;}
256     if (num==6) {w=0.555555555555556;uv[0]=0.774596669241483;uv[1]=0.577350269189626;}
257     return;
258     }
259     }
260 francois 638 int FEM_QUADRANGLE8::get_nb_fonction_interpolation(void)
261 francois 310 {
262     return 8;
263     }
264    
265 francois 638 int FEM_QUADRANGLE8::get_degremax_fonction_interpolation(void)
266     {
267     return 3;
268     }
269 francois 310 double FEM_QUADRANGLE8::get_fonction_interpolation(int num,double *uv)
270     {
271     double val;
272     switch (num)
273     {
274     case 1:
275     val=-0.25*(1-uv[0])*(1-uv[1])*(1.+uv[0]+uv[1]);
276     break;
277     case 2:
278     val=0.5*(1-uv[0]*uv[0])*(1-uv[1]);
279     break;
280     case 3:
281     val=-0.25*(1+uv[0])*(1-uv[1])*(1.-uv[0]+uv[1]);
282     break;
283     case 4:
284     val=0.5*(1+uv[0])*(1-uv[1]*uv[1]);
285     break;
286     case 5:
287     val=-0.25*(1+uv[0])*(1+uv[1])*(1.-uv[0]-uv[1]);
288     break;
289     case 6:
290     val=0.5*(1-uv[0]*uv[0])*(1+uv[1]);
291     break;
292     case 7:
293     val=-0.25*(1-uv[0])*(1+uv[1])*(1.+uv[0]-uv[1]);
294     break;
295     case 8:
296     val=0.5*(1-uv[0])*(1-uv[1]*uv[1]);
297     break;
298     }
299     return val;
300     }
301    
302     double FEM_QUADRANGLE8::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
303     {
304     double val;
305    
306     switch (num)
307     {
308     case 1:
309     switch (num_variable)
310     {
311     case 1:
312     val=0.25*(1-uv[1])*(2*uv[0]+uv[1]);
313     break;
314     case 2:
315     val=0.25*(1-uv[0])*(uv[0]+2.*uv[1]);
316     break;
317     } break;
318     case 2:
319     switch (num_variable)
320     {
321     case 1:
322     val=-uv[0]*(1.-uv[1]);
323     break;
324     case 2:
325     val=-0.5*(1-uv[0]*uv[0]);
326     break;
327     }break;
328     case 3:
329     switch (num_variable)
330     {
331     case 1:
332     val=0.25*(1.-uv[1])*(2*uv[0]-uv[1]);
333     break;
334     case 2:
335     val=-0.25*(1.+uv[0])*(uv[0]-2.*uv[1]);
336     break;
337     }break;
338     case 4:
339     switch (num_variable)
340     {
341     case 1:
342     val=0.5*(1.-uv[1]*uv[1]);
343     break;
344     case 2:
345     val=-uv[1]*(1+uv[0]);
346     break;
347     }break;
348     case 5:
349     switch (num_variable)
350     {
351     case 1:
352     val=0.25*(1+uv[1])*(2*uv[0]+uv[1]);
353     break;
354     case 2:
355     val=0.25*(1+uv[0])*(uv[0]+2.*uv[1]);
356     break;
357     } break;
358     case 6:
359     switch (num_variable)
360     {
361     case 1:
362     val=-uv[0]*(1.+uv[1]);
363     break;
364     case 2:
365     val=-0.5*(1-uv[0]*uv[0]);
366     break;
367     }break;
368     case 7:
369     switch (num_variable)
370     {
371     case 1:
372     val=0.25*(1.+uv[1])*(2*uv[0]-uv[1]);
373     break;
374     case 2:
375     val=-0.25*(1.-uv[0])*(uv[0]-2.*uv[1]);
376     break;
377     }break;
378     case 8:
379     switch (num_variable)
380     {
381     case 1:
382     val=-0.5*(1.-uv[1]*uv[1]);
383     break;
384     case 2:
385     val=-uv[1]*(1+uv[0]);
386     break;
387     }break;
388    
389     }
390     return val;
391     }
392    
393     double FEM_QUADRANGLE8::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
394     {
395    
396     }
397    
398     void FEM_QUADRANGLE8::get_inverse_jacob(double* j,double *uv,double unite)
399     {
400    
401     }