ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/mtu/src/fem_quadrangle4.cpp
Revision: 1158
Committed: Thu Jun 13 22:18:49 2024 UTC (11 months, 1 week ago) by francois
File size: 10544 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_quadrangle4.cpp
15     //####//
16     //####//------------------------------------------------------------
17     //####//------------------------------------------------------------
18     //####// COPYRIGHT 2000-2024
19     //####// jeu 13 jun 2024 11:58:54 EDT
20     //####//------------------------------------------------------------
21     //####//------------------------------------------------------------
22 francois 310
23    
24     #include "gestionversion.h"
25     #include "fem_quadrangle4.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 310
31    
32     FEM_QUADRANGLE4::FEM_QUADRANGLE4(unsigned long num,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,mai),FEM_TEMPLATE_ELEMENT<4>(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     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
40     }
41    
42     FEM_QUADRANGLE4::FEM_QUADRANGLE4(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(mai),FEM_TEMPLATE_ELEMENT<4>(tabnoeud)
43     {
44     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
45     tab[0]->get_lien_element2()->ajouter(this);
46     tab[1]->get_lien_element2()->ajouter(this);
47     tab[2]->get_lien_element2()->ajouter(this);
48     tab[3]->get_lien_element2()->ajouter(this);
49     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
50     }
51    
52 francois 378 FEM_QUADRANGLE4::FEM_QUADRANGLE4(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,topo),FEM_TEMPLATE_ELEMENT<4>(tabnoeud)
53     {
54     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
55     tab[0]->get_lien_element2()->ajouter(this);
56     tab[1]->get_lien_element2()->ajouter(this);
57     tab[2]->get_lien_element2()->ajouter(this);
58     tab[3]->get_lien_element2()->ajouter(this);
59     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
60     }
61    
62     FEM_QUADRANGLE4::FEM_QUADRANGLE4(class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(topo),FEM_TEMPLATE_ELEMENT<4>(tabnoeud)
63     {
64     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
65     tab[0]->get_lien_element2()->ajouter(this);
66     tab[1]->get_lien_element2()->ajouter(this);
67     tab[2]->get_lien_element2()->ajouter(this);
68     tab[3]->get_lien_element2()->ajouter(this);
69     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
70     }
71     FEM_QUADRANGLE4::FEM_QUADRANGLE4(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<4>(tabnoeud)
72     {
73     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
74     tab[0]->get_lien_element2()->ajouter(this);
75     tab[1]->get_lien_element2()->ajouter(this);
76     tab[2]->get_lien_element2()->ajouter(this);
77     tab[3]->get_lien_element2()->ajouter(this);
78     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
79     }
80    
81     FEM_QUADRANGLE4::FEM_QUADRANGLE4(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(topo,mai),FEM_TEMPLATE_ELEMENT<4>(tabnoeud)
82     {
83     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
84     tab[0]->get_lien_element2()->ajouter(this);
85     tab[1]->get_lien_element2()->ajouter(this);
86     tab[2]->get_lien_element2()->ajouter(this);
87     tab[3]->get_lien_element2()->ajouter(this);
88     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
89     }
90    
91    
92    
93    
94    
95    
96    
97    
98    
99    
100    
101 francois 310 FEM_QUADRANGLE4::FEM_QUADRANGLE4(FEM_QUADRANGLE4& mdd):FEM_ELEMENT2(mdd),FEM_TEMPLATE_ELEMENT<4>(mdd)
102     {
103     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
104     tab[0]->get_lien_element2()->ajouter(this);
105     tab[1]->get_lien_element2()->ajouter(this);
106     tab[2]->get_lien_element2()->ajouter(this);
107     tab[3]->get_lien_element2()->ajouter(this);
108     get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
109     }
110     FEM_QUADRANGLE4::~FEM_QUADRANGLE4()
111     {
112 francois 663 if (liaison_topologique!=NULL)
113     if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
114 francois 310 tab[0]->get_lien_element2()->supprimer(this);
115     tab[1]->get_lien_element2()->supprimer(this);
116     tab[2]->get_lien_element2()->supprimer(this);
117     tab[3]->get_lien_element2()->supprimer(this);
118     get_fem_noeudpetitid()->get_lien_petit_element2()->supprimer(this);
119     }
120    
121    
122     FEM_ELEMENT_MAILLAGE* FEM_QUADRANGLE4::dupliquer(FEM_MAILLAGE *femmai,long decalage)
123     {
124     FEM_NOEUD* tabnoeud[4];
125     tabnoeud[0]=femmai->get_fem_noeudid(tab[0]->get_id()+decalage);
126     tabnoeud[1]=femmai->get_fem_noeudid(tab[1]->get_id()+decalage);
127     tabnoeud[2]=femmai->get_fem_noeudid(tab[2]->get_id()+decalage);
128     tabnoeud[3]=femmai->get_fem_noeudid(tab[3]->get_id()+decalage);
129     FEM_QUADRANGLE4* quad=new FEM_QUADRANGLE4(get_id()+decalage,maillage,tabnoeud);
130     femmai->ajouter_fem_element2(quad);
131     return quad;
132     }
133    
134 francois 786 void FEM_QUADRANGLE4::get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin)
135     {
136     voisin.vide();
137     if (no==tab[0])
138     {
139     voisin.ajouter(tab[1]);
140     voisin.ajouter(tab[3]);
141     }
142     if (no==tab[1])
143     {
144     voisin.ajouter(tab[0]);
145     voisin.ajouter(tab[2]);
146     }
147     if (no==tab[2])
148     {
149     voisin.ajouter(tab[1]);
150     voisin.ajouter(tab[3]);
151     }
152     if (no==tab[3])
153     {
154     voisin.ajouter(tab[0]);
155     voisin.ajouter(tab[2]);
156     }
157     }
158 francois 310
159    
160    
161    
162     int FEM_QUADRANGLE4::get_type_entite(void)
163     {
164 francois 1150 return MAGIC::TYPE_ENTITE::IDFEM_QUADRANGLE4;
165 francois 310 }
166    
167     int FEM_QUADRANGLE4::get_dimension(void)
168     {
169     return 2;
170     }
171    
172    
173    
174     int FEM_QUADRANGLE4::get_nb_fem_noeud(void)
175     {
176     return FEM_TEMPLATE_ELEMENT<4>::get_nb_fem_noeud();
177     }
178    
179     FEM_NOEUD* FEM_QUADRANGLE4::get_fem_noeud(int num)
180     {
181     return FEM_TEMPLATE_ELEMENT<4>::get_fem_noeud(num);
182     }
183    
184     void FEM_QUADRANGLE4::change_noeud(int num,FEM_NOEUD* noeud)
185     {
186     FEM_TEMPLATE_ELEMENT<4>::change_noeud(num,noeud);
187     }
188    
189     BOITE_3D& FEM_QUADRANGLE4::get_boite_3D(void)
190     {
191     return FEM_TEMPLATE_ELEMENT<4>::get_boite_3D();
192     }
193    
194    
195 francois 763 void FEM_QUADRANGLE4::enregistrer(std::ostream& o,double version)
196 francois 310 {
197 francois 378 if (maillage!=NULL)
198     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_QUADRANGLE4($"<< get_lien_topologie()->get_id() << ",$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ",$" << tab[3]->get_id()<< ");" << std::endl;
199     else o << "%" << get_id() << "=FEM_QUADRANGLE4(NULL,$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ",$" << tab[3]->get_id()<< ");" << std::endl;
200     else
201     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_QUADRANGLE4($"<< get_lien_topologie()->get_id() << ",NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ",$" << tab[3]->get_id()<< ");" << std::endl;
202     else o << "%" << get_id() << "=FEM_QUADRANGLE4(NULL,NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id()<< ",$" << tab[3]->get_id()<< ");" << std::endl;
203    
204 francois 310 }
205 francois 637
206    
207     int FEM_QUADRANGLE4::get_nb_pt_gauss(int degre)
208     {
209 francois 755 return OT_POINTS_GAUSS::get_nb_point_quad_prod(degre);
210 francois 637 }
211    
212     void FEM_QUADRANGLE4::get_pt_gauss(int degre,int num,double &w,double *uv)
213     {
214 francois 755 return OT_POINTS_GAUSS::get_pt_gauss_qua_prod(degre,num,w,uv);
215 francois 637 }
216 francois 757
217    
218     int FEM_QUADRANGLE4::get_degre_gauss(int num)
219     {
220     return OT_POINTS_GAUSS::get_degre_gauss_quad_prod(num);
221     }
222    
223    
224 francois 638 int FEM_QUADRANGLE4::get_nb_fonction_interpolation(void)
225 francois 310 {
226     return 4;
227     }
228 francois 638 int FEM_QUADRANGLE4::get_degremax_fonction_interpolation(void)
229     {
230     return 2;
231     }
232 francois 310
233     double FEM_QUADRANGLE4::get_fonction_interpolation(int num,double *uv)
234     {
235     double val;
236     switch (num)
237     {
238     case 1:
239     val=0.25*(1-uv[0])*(1-uv[1]);
240     break;
241     case 2:
242     val=0.25*(1+uv[0])*(1-uv[1]);
243     break;
244     case 3:
245     val=0.25*(1+uv[0])*(1+uv[1]);
246     break;
247     case 4:
248     val=0.25*(1-uv[0])*(1+uv[1]);
249     break;
250     }
251     return val;
252     }
253    
254     double FEM_QUADRANGLE4::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
255     {
256     double val;
257    
258     switch (num)
259     {
260     case 1:
261     switch (num_variable)
262     {
263     case 1:
264     val=0.25*(-1+uv[1]);
265     break;
266     case 2:
267     val=0.25*(-1+uv[0]);
268     break;
269     } break;
270     case 2:
271     switch (num_variable)
272     {
273     case 1:
274     val=0.25*(1-uv[1]);
275     break;
276     case 2:
277     val=0.25*(-1-uv[0]);
278     break;
279     }break;
280     case 3:
281     switch (num_variable)
282     {
283     case 1:
284     val=0.25*(1+uv[1]);
285     break;
286     case 2:
287     val=0.25*(1+uv[0]);
288     break;
289     }break;
290     case 4:
291     switch (num_variable)
292     {
293     case 1:
294     val=0.25*(-1-uv[1]);
295     break;
296     case 2:
297     val=0.25*(1-uv[0]);
298     break;
299     }break;
300    
301     }
302     return val;
303     }
304    
305 francois 635 bool FEM_QUADRANGLE4::valide_parametre_EF(double* uvw)
306     {
307 francois 674 if (uvw[0]>=-1.-1e-10)
308     if (uvw[1]>=-1.-1e-10)
309 francois 635 if (uvw[0]<=1.+1e-10)
310     if (uvw[1]<=1.+1e-10)
311     return true;
312     return false;
313     }
314    
315 francois 684
316     void FEM_QUADRANGLE4::reinit_boite_3D(void)
317     {
318     FEM_TEMPLATE_ELEMENT<4>::reinit_boite_3D();
319     }