ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/mtu/src/fem_triangle3.cpp
Revision: 1158
Committed: Thu Jun 13 22:18:49 2024 UTC (11 months, 1 week ago) by francois
File size: 9665 byte(s)
Log Message: 
compatibilité Ubuntu 22.04
Suppression des refeences à Windows
Ajout d'une banière

File Contents

# User Rev Content
1 francois 1158 //####//------------------------------------------------------------
2     //####//------------------------------------------------------------
3     //####// MAGiC
4     //####// Jean Christophe Cuilliere et Vincent FRANCOIS
5     //####// Departement de Genie Mecanique - UQTR
6     //####//------------------------------------------------------------
7     //####// MAGIC est un projet de recherche de l equipe ERICCA
8     //####// du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
9     //####// http://www.uqtr.ca/ericca
10     //####// http://www.uqtr.ca/
11     //####//------------------------------------------------------------
12     //####//------------------------------------------------------------
13     //####//
14     //####// fem_triangle3.cpp
15     //####//
16     //####//------------------------------------------------------------
17     //####//------------------------------------------------------------
18     //####// COPYRIGHT 2000-2024
19     //####// jeu 13 jun 2024 11:58:55 EDT
20     //####//------------------------------------------------------------
21     //####//------------------------------------------------------------
22 francois 283
23    
24     #include "gestionversion.h"
25     #include "fem_triangle3.h"
26     #include "fem_maillage.h"
27     #include "fem_noeud.h"
28     #include "mg_element_maillage.h"
29 francois 754 #include "ot_quadrature_gauss.h"
30 francois 283
31    
32 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)
33 francois 283 {
34     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
35 francois 309 tab[0]->get_lien_element2()->ajouter(this);
36     tab[1]->get_lien_element2()->ajouter(this);
37     tab[2]->get_lien_element2()->ajouter(this);
38     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
39 francois 283 }
40    
41 francois 309 FEM_TRIANGLE3::FEM_TRIANGLE3(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(mai),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
42 francois 283 {
43     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
44 francois 309 tab[0]->get_lien_element2()->ajouter(this);
45     tab[1]->get_lien_element2()->ajouter(this);
46     tab[2]->get_lien_element2()->ajouter(this);
47     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
48 francois 283 }
49 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)
50     {
51     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
52     tab[0]->get_lien_element2()->ajouter(this);
53     tab[1]->get_lien_element2()->ajouter(this);
54     tab[2]->get_lien_element2()->ajouter(this);
55     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
56     }
57 francois 283
58 francois 378 FEM_TRIANGLE3::FEM_TRIANGLE3(class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(topo),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
59     {
60     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
61     tab[0]->get_lien_element2()->ajouter(this);
62     tab[1]->get_lien_element2()->ajouter(this);
63     tab[2]->get_lien_element2()->ajouter(this);
64     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
65     }
66     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)
67     {
68     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
69     tab[0]->get_lien_element2()->ajouter(this);
70     tab[1]->get_lien_element2()->ajouter(this);
71     tab[2]->get_lien_element2()->ajouter(this);
72     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
73     }
74    
75     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)
76     {
77     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
78     tab[0]->get_lien_element2()->ajouter(this);
79     tab[1]->get_lien_element2()->ajouter(this);
80     tab[2]->get_lien_element2()->ajouter(this);
81     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
82     }
83 francois 309 FEM_TRIANGLE3::FEM_TRIANGLE3(FEM_TRIANGLE3& mdd):FEM_ELEMENT2(mdd),FEM_TEMPLATE_ELEMENT<3>(mdd)
84 francois 283 {
85     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
86 francois 309 tab[0]->get_lien_element2()->ajouter(this);
87     tab[1]->get_lien_element2()->ajouter(this);
88     tab[2]->get_lien_element2()->ajouter(this);
89     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
90 francois 283 }
91     FEM_TRIANGLE3::~FEM_TRIANGLE3()
92     {
93 francois 663 if (liaison_topologique!=NULL)
94     if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
95 francois 309 tab[0]->get_lien_element2()->supprimer(this);
96     tab[1]->get_lien_element2()->supprimer(this);
97     tab[2]->get_lien_element2()->supprimer(this);
98     get_fem_noeudpetitid()->get_lien_petit_element2()->supprimer(this);
99 francois 283 }
100    
101    
102     FEM_ELEMENT_MAILLAGE* FEM_TRIANGLE3::dupliquer(FEM_MAILLAGE *femmai,long decalage)
103     {
104     FEM_NOEUD* tabnoeud[3];
105     tabnoeud[0]=femmai->get_fem_noeudid(tab[0]->get_id()+decalage);
106     tabnoeud[1]=femmai->get_fem_noeudid(tab[1]->get_id()+decalage);
107     tabnoeud[2]=femmai->get_fem_noeudid(tab[2]->get_id()+decalage);
108     FEM_TRIANGLE3* tri=new FEM_TRIANGLE3(get_id()+decalage,maillage,tabnoeud);
109 francois 309 femmai->ajouter_fem_element2(tri);
110 francois 283 return tri;
111     }
112    
113 francois 786 void FEM_TRIANGLE3::get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin)
114     {
115     voisin.vide();
116     if (no==tab[0])
117     {
118     voisin.ajouter(tab[1]);
119     voisin.ajouter(tab[2]);
120     }
121     if (no==tab[1])
122     {
123     voisin.ajouter(tab[0]);
124     voisin.ajouter(tab[2]);
125     }
126     if (no==tab[2])
127     {
128     voisin.ajouter(tab[0]);
129     voisin.ajouter(tab[1]);
130     }
131     }
132 francois 283
133 francois 684 void FEM_TRIANGLE3::reinit_boite_3D(void)
134     {
135     FEM_TEMPLATE_ELEMENT<3>::reinit_boite_3D();
136     }
137 francois 283
138    
139     int FEM_TRIANGLE3::get_type_entite(void)
140     {
141 francois 1150 return MAGIC::TYPE_ENTITE::IDFEM_TRIANGLE3;
142 francois 283 }
143    
144     int FEM_TRIANGLE3::get_dimension(void)
145     {
146     return 2;
147     }
148    
149    
150     int FEM_TRIANGLE3::get_nb_fem_noeud(void)
151     {
152     return FEM_TEMPLATE_ELEMENT<3>::get_nb_fem_noeud();
153     }
154    
155     FEM_NOEUD* FEM_TRIANGLE3::get_fem_noeud(int num)
156     {
157     return FEM_TEMPLATE_ELEMENT<3>::get_fem_noeud(num);
158     }
159    
160     void FEM_TRIANGLE3::change_noeud(int num,FEM_NOEUD* noeud)
161     {
162     FEM_TEMPLATE_ELEMENT<3>::change_noeud(num,noeud);
163     }
164    
165     BOITE_3D& FEM_TRIANGLE3::get_boite_3D(void)
166     {
167     return FEM_TEMPLATE_ELEMENT<3>::get_boite_3D();
168     }
169    
170    
171 francois 763 void FEM_TRIANGLE3::enregistrer(std::ostream& o,double version)
172 francois 283 {
173 francois 378 if (maillage!=NULL)
174     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;
175     else o << "%" << get_id() << "=FEM_TRIANGLE3(NULL,$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ");" << std::endl;
176     else
177     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;
178     else o << "%" << get_id() << "=FEM_TRIANGLE3(NULL,NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ");" << std::endl;
179    
180 francois 637 }
181 francois 310
182 francois 637 int FEM_TRIANGLE3::get_nb_pt_gauss(int degre)
183     {
184 francois 754 return OT_POINTS_GAUSS::get_nb_point_tri(degre);
185 francois 637 }
186    
187     void FEM_TRIANGLE3::get_pt_gauss(int degre,int num,double &w,double *uv)
188     {
189 francois 754 return OT_POINTS_GAUSS::get_pt_gauss_tri(degre,num,w,uv);
190 francois 637 }
191 francois 757
192     int FEM_TRIANGLE3::get_degre_gauss(int num)
193     {
194     return OT_POINTS_GAUSS::get_degre_gauss_tri(num);
195     }
196 francois 637
197 francois 638 int FEM_TRIANGLE3::get_nb_fonction_interpolation(void)
198 francois 310 {
199     return 3;
200     }
201 francois 638 int FEM_TRIANGLE3::get_degremax_fonction_interpolation(void)
202     {
203     return 1;
204     }
205 francois 310 double FEM_TRIANGLE3::get_fonction_interpolation(int num,double *uv)
206     {
207     double val;
208     switch (num)
209     {
210     case 1:
211     val=1-uv[0]-uv[1];
212     break;
213     case 2:
214     val=uv[0];
215     break;
216     case 3:
217     val=uv[1];
218     break;
219    
220     }
221     return val;
222     }
223    
224     double FEM_TRIANGLE3::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
225     {
226     double val;
227    
228     switch (num)
229     {
230     case 1:
231     switch (num_variable)
232     {
233     case 1:
234     val=-1;
235     break;
236     case 2:
237     val=-1;
238     break;
239     } break;
240     case 2:
241     switch (num_variable)
242     {
243     case 1:
244     val=1;
245     break;
246     case 2:
247     val=0.;
248     break;
249     }break;
250     case 3:
251     switch (num_variable)
252     {
253     case 1:
254     val=0.;
255     break;
256     case 2:
257     val=1.;
258     break;
259     }break;
260    
261     }
262     return val;
263     }
264    
265 francois 1118 bool FEM_TRIANGLE3::get_param_element_fini_2D(double *xyz,double *uv)
266     {
267     double *xyz1=get_fem_noeud(0)->get_coord();
268     double *xyz2=get_fem_noeud(1)->get_coord();
269     double *xyz3=get_fem_noeud(2)->get_coord();
270     OPERATEUR::calculeparamEFdansletriangle(xyz1,xyz2,xyz3,xyz[0],xyz[1],xyz[2],uv[0],uv[1]);
271     return true;
272     }
273 francois 310
274 francois 635 bool FEM_TRIANGLE3::valide_parametre_EF(double* uvw)
275     {
276     if (uvw[0]>=-1e-10)
277     if (uvw[1]>=-1e-10)
278     if (uvw[0]+uvw[1]<=1.+1e-10)
279     return true;
280     return false;
281     }