ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/geometrie/src/fem_triangle3.cpp
Revision: 754
Committed: Thu Oct 29 21:49:15 2015 UTC (9 years, 6 months ago) by francois
File size: 10523 byte(s)
Log Message: 
Centralisation et unification de la définition des points de gauss des elements finis et arrangement dans le meme ordre que Aster

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_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_triangle3.h"
27     #include "fem_maillage.h"
28     #include "fem_noeud.h"
29     #include "mg_element_maillage.h"
30 francois 754 #include "ot_quadrature_gauss.h"
31 francois 283
32    
33 francois 309 FEM_TRIANGLE3::FEM_TRIANGLE3(unsigned long num,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,mai),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
34 francois 283 {
35     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
36 francois 309 tab[0]->get_lien_element2()->ajouter(this);
37     tab[1]->get_lien_element2()->ajouter(this);
38     tab[2]->get_lien_element2()->ajouter(this);
39     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
40 francois 283 }
41    
42 francois 309 FEM_TRIANGLE3::FEM_TRIANGLE3(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(mai),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
43 francois 283 {
44     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
45 francois 309 tab[0]->get_lien_element2()->ajouter(this);
46     tab[1]->get_lien_element2()->ajouter(this);
47     tab[2]->get_lien_element2()->ajouter(this);
48     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
49 francois 283 }
50 francois 378 FEM_TRIANGLE3::FEM_TRIANGLE3(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,topo),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
51     {
52     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
53     tab[0]->get_lien_element2()->ajouter(this);
54     tab[1]->get_lien_element2()->ajouter(this);
55     tab[2]->get_lien_element2()->ajouter(this);
56     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
57     }
58 francois 283
59 francois 378 FEM_TRIANGLE3::FEM_TRIANGLE3(class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(topo),FEM_TEMPLATE_ELEMENT<3>(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     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
66     }
67     FEM_TRIANGLE3::FEM_TRIANGLE3(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<3>(tabnoeud)
68     {
69     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
70     tab[0]->get_lien_element2()->ajouter(this);
71     tab[1]->get_lien_element2()->ajouter(this);
72     tab[2]->get_lien_element2()->ajouter(this);
73     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
74     }
75    
76     FEM_TRIANGLE3::FEM_TRIANGLE3(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(topo,mai),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
77     {
78     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
79     tab[0]->get_lien_element2()->ajouter(this);
80     tab[1]->get_lien_element2()->ajouter(this);
81     tab[2]->get_lien_element2()->ajouter(this);
82     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
83     }
84 francois 309 FEM_TRIANGLE3::FEM_TRIANGLE3(FEM_TRIANGLE3& mdd):FEM_ELEMENT2(mdd),FEM_TEMPLATE_ELEMENT<3>(mdd)
85 francois 283 {
86     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
87 francois 309 tab[0]->get_lien_element2()->ajouter(this);
88     tab[1]->get_lien_element2()->ajouter(this);
89     tab[2]->get_lien_element2()->ajouter(this);
90     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
91 francois 283 }
92     FEM_TRIANGLE3::~FEM_TRIANGLE3()
93     {
94 francois 663 if (liaison_topologique!=NULL)
95     if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
96 francois 309 tab[0]->get_lien_element2()->supprimer(this);
97     tab[1]->get_lien_element2()->supprimer(this);
98     tab[2]->get_lien_element2()->supprimer(this);
99     get_fem_noeudpetitid()->get_lien_petit_element2()->supprimer(this);
100 francois 283 }
101    
102    
103     FEM_ELEMENT_MAILLAGE* FEM_TRIANGLE3::dupliquer(FEM_MAILLAGE *femmai,long decalage)
104     {
105     FEM_NOEUD* tabnoeud[3];
106     tabnoeud[0]=femmai->get_fem_noeudid(tab[0]->get_id()+decalage);
107     tabnoeud[1]=femmai->get_fem_noeudid(tab[1]->get_id()+decalage);
108     tabnoeud[2]=femmai->get_fem_noeudid(tab[2]->get_id()+decalage);
109     FEM_TRIANGLE3* tri=new FEM_TRIANGLE3(get_id()+decalage,maillage,tabnoeud);
110 francois 309 femmai->ajouter_fem_element2(tri);
111 francois 283 return tri;
112     }
113    
114    
115 francois 684 void FEM_TRIANGLE3::reinit_boite_3D(void)
116     {
117     FEM_TEMPLATE_ELEMENT<3>::reinit_boite_3D();
118     }
119 francois 283
120    
121     int FEM_TRIANGLE3::get_type_entite(void)
122     {
123     return IDFEM_TRIANGLE3;
124     }
125    
126     int FEM_TRIANGLE3::get_dimension(void)
127     {
128     return 2;
129     }
130    
131    
132    
133     int FEM_TRIANGLE3::get_nb_fem_noeud(void)
134     {
135     return FEM_TEMPLATE_ELEMENT<3>::get_nb_fem_noeud();
136     }
137    
138     FEM_NOEUD* FEM_TRIANGLE3::get_fem_noeud(int num)
139     {
140     return FEM_TEMPLATE_ELEMENT<3>::get_fem_noeud(num);
141     }
142    
143     void FEM_TRIANGLE3::change_noeud(int num,FEM_NOEUD* noeud)
144     {
145     FEM_TEMPLATE_ELEMENT<3>::change_noeud(num,noeud);
146     }
147    
148     BOITE_3D& FEM_TRIANGLE3::get_boite_3D(void)
149     {
150     return FEM_TEMPLATE_ELEMENT<3>::get_boite_3D();
151     }
152    
153    
154     void FEM_TRIANGLE3::enregistrer(std::ostream& o)
155     {
156 francois 378 if (maillage!=NULL)
157     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_TRIANGLE3($"<< get_lien_topologie()->get_id() << ",$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ");" << std::endl;
158     else o << "%" << get_id() << "=FEM_TRIANGLE3(NULL,$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ");" << std::endl;
159     else
160     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_TRIANGLE3($"<< get_lien_topologie()->get_id() << ",NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ");" << std::endl;
161     else o << "%" << get_id() << "=FEM_TRIANGLE3(NULL,NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ");" << std::endl;
162    
163 francois 637 }
164 francois 310
165 francois 637 int FEM_TRIANGLE3::get_nb_pt_gauss(int degre)
166     {
167 francois 754 return OT_POINTS_GAUSS::get_nb_point_tri(degre);
168 francois 637 }
169    
170     void FEM_TRIANGLE3::get_pt_gauss(int degre,int num,double &w,double *uv)
171     {
172 francois 754 return OT_POINTS_GAUSS::get_pt_gauss_tri(degre,num,w,uv);
173 francois 637 }
174    
175 francois 638 int FEM_TRIANGLE3::get_nb_fonction_interpolation(void)
176 francois 310 {
177     return 3;
178     }
179 francois 638 int FEM_TRIANGLE3::get_degremax_fonction_interpolation(void)
180     {
181     return 1;
182     }
183 francois 310 double FEM_TRIANGLE3::get_fonction_interpolation(int num,double *uv)
184     {
185     double val;
186     switch (num)
187     {
188     case 1:
189     val=1-uv[0]-uv[1];
190     break;
191     case 2:
192     val=uv[0];
193     break;
194     case 3:
195     val=uv[1];
196     break;
197    
198     }
199     return val;
200     }
201    
202     double FEM_TRIANGLE3::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
203     {
204     double val;
205    
206     switch (num)
207     {
208     case 1:
209     switch (num_variable)
210     {
211     case 1:
212     val=-1;
213     break;
214     case 2:
215     val=-1;
216     break;
217     } break;
218     case 2:
219     switch (num_variable)
220     {
221     case 1:
222     val=1;
223     break;
224     case 2:
225     val=0.;
226     break;
227     }break;
228     case 3:
229     switch (num_variable)
230     {
231     case 1:
232     val=0.;
233     break;
234     case 2:
235     val=1.;
236     break;
237     }break;
238    
239     }
240     return val;
241     }
242    
243     double FEM_TRIANGLE3::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
244 gervaislavoie 382 {
245 francois 647 int nb=get_nb_fem_noeud();
246     OT_VECTEUR_3D *point= new OT_VECTEUR_3D[nb];
247 gervaislavoie 382
248 francois 647 for (int i=0;i<nb;i++)
249     {
250     point[i].change_x(get_fem_noeud(i)->get_x());
251     point[i].change_y(get_fem_noeud(i)->get_y());
252     point[i].change_z(get_fem_noeud(i)->get_z());
253     }
254    
255     OT_MATRICE_3D J;
256     for (int i=0;i<2;i++)
257     {
258     J(i,0)=0.;J(i,1)=0.;J(i,2)=0.;
259     for (int k=0;k<nb;k++)
260     {
261     J(i,0)=J(i,0)+get_fonction_derive_interpolation(k+1,i+1,uv)*point[k].get_x();
262     J(i,1)=J(i,1)+get_fonction_derive_interpolation(k+1,i+1,uv)*point[k].get_y();
263     J(i,2)=J(i,2)+get_fonction_derive_interpolation(k+1,i+1,uv)*point[k].get_z();
264     }
265     }
266     OT_VECTEUR_3D J1(J(0,0),J(0,1),J(0,2));
267     OT_VECTEUR_3D J2(J(1,0),J(1,1),J(1,2));
268     OT_VECTEUR_3D J3=J1&J2;
269     J3.norme();
270     J(2,0)=J3(0);
271     J(2,1)=J3(1);
272     J(2,2)=J3(2);
273 gervaislavoie 389
274 francois 647 jac[0]=J(0,0);
275     jac[1]=J(0,1);
276     jac[2]=J(0,2);
277 gervaislavoie 389
278 francois 647 jac[3]=J(1,0);
279     jac[4]=J(1,1);
280     jac[5]=J(1,2);
281 gervaislavoie 389
282 francois 647 jac[6]=J(2,0);
283     jac[7]=J(2,1);
284     jac[8]=J(2,2);
285     li=3;col=3;
286 gervaislavoie 389
287 francois 647 return J.get_determinant();
288 francois 310 }
289    
290 francois 635
291     bool FEM_TRIANGLE3::valide_parametre_EF(double* uvw)
292     {
293     if (uvw[0]>=-1e-10)
294     if (uvw[1]>=-1e-10)
295     if (uvw[0]+uvw[1]<=1.+1e-10)
296     return true;
297     return false;
298     }
299    
300 francois 310 void FEM_TRIANGLE3::get_inverse_jacob(double* j,double *uv,double unite)
301     {
302 francois 647 double jac[9];
303     int li,col;
304     get_jacobien(jac,uv,li,col);
305     OT_MATRICE_3D J;
306     J(0,0)=jac[0];
307     J(0,1)=jac[1];
308     J(0,2)=jac[2];
309     J(1,0)=jac[3];
310     J(1,1)=jac[4];
311     J(1,2)=jac[5];
312     J(2,0)=jac[6];
313     J(2,1)=jac[7];
314     J(2,2)=jac[8];
315     OT_MATRICE_3D j_i=J.inverse();
316    
317 gervaislavoie 382
318     j[0]=j_i(0,0);
319 francois 647 j[1]=j_i(0,1);
320     j[2]=j_i(0,2);
321 gervaislavoie 382
322 francois 647 j[3]=j_i(1,0);
323 gervaislavoie 382 j[4]=j_i(1,1);
324 francois 647 j[5]=j_i(1,2);
325 gervaislavoie 382
326 francois 647 j[6]=j_i(2,0);
327     j[7]=j_i(2,1);
328 gervaislavoie 382 j[8]=j_i(2,2);
329 francois 310 }