ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/geometrie/src/fem_tetra4.cpp
Revision: 283
Committed: Tue Sep 13 21:11:20 2011 UTC (13 years, 8 months ago) by francois
File size: 8526 byte(s)
Log Message: 
structure de l'écriture

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_tetra4.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_tetra4.h"
27     #include "fem_noeud.h"
28     #include "mg_element_maillage.h"
29     #include "fem_maillage.h"
30    
31    
32     FEM_TETRA4::FEM_TETRA4(unsigned long num,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_TETRA(num,mai),FEM_TEMPLATE_ELEMENT<4>(tabnoeud)
33     {
34     if (liaison_topologique!=NULL)
35     if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
36     tab[0]->get_lien_tetra()->ajouter(this);
37     tab[1]->get_lien_tetra()->ajouter(this);
38     tab[2]->get_lien_tetra()->ajouter(this);
39     tab[3]->get_lien_tetra()->ajouter(this);
40     get_fem_noeudpetitid()->get_lien_petit_tetra()->ajouter(this);
41     }
42    
43     FEM_TETRA4::FEM_TETRA4(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_TETRA(mai),FEM_TEMPLATE_ELEMENT<4>(tabnoeud)
44     {
45     if (liaison_topologique!=NULL)
46     if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
47     tab[0]->get_lien_tetra()->ajouter(this);
48     tab[1]->get_lien_tetra()->ajouter(this);
49     tab[2]->get_lien_tetra()->ajouter(this);
50     tab[3]->get_lien_tetra()->ajouter(this);
51     get_fem_noeudpetitid()->get_lien_petit_tetra()->ajouter(this);
52     }
53    
54     FEM_TETRA4::FEM_TETRA4(FEM_TETRA4& mdd):FEM_TETRA(mdd),FEM_TEMPLATE_ELEMENT<4>(mdd)
55     {
56     if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
57     tab[0]->get_lien_tetra()->ajouter(this);
58     tab[1]->get_lien_tetra()->ajouter(this);
59     tab[2]->get_lien_tetra()->ajouter(this);
60     tab[3]->get_lien_tetra()->ajouter(this);
61     get_fem_noeudpetitid()->get_lien_petit_tetra()->ajouter(this);
62     }
63    
64     FEM_TETRA4::~FEM_TETRA4()
65     {
66     if (liaison_topologique==NULL) return;
67     if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
68     tab[0]->get_lien_tetra()->supprimer(this);
69     tab[1]->get_lien_tetra()->supprimer(this);
70     tab[2]->get_lien_tetra()->supprimer(this);
71     tab[3]->get_lien_tetra()->supprimer(this);
72     get_fem_noeudpetitid()->get_lien_petit_tetra()->supprimer(this);
73     }
74    
75    
76     FEM_ELEMENT_MAILLAGE* FEM_TETRA4::dupliquer(FEM_MAILLAGE *femmai,long decalage)
77     {
78     FEM_NOEUD* tabnoeud[4];
79     tabnoeud[0]=femmai->get_fem_noeudid(tab[0]->get_id()+decalage);
80     tabnoeud[1]=femmai->get_fem_noeudid(tab[1]->get_id()+decalage);
81     tabnoeud[2]=femmai->get_fem_noeudid(tab[2]->get_id()+decalage);
82     tabnoeud[3]=femmai->get_fem_noeudid(tab[3]->get_id()+decalage);
83     FEM_TETRA4* tet=new FEM_TETRA4(get_id()+decalage,maillage,tabnoeud);
84     femmai->ajouter_fem_tetra(tet);
85     return tet;
86     }
87    
88    
89     int FEM_TETRA4::get_type_entite(void)
90     {
91     return IDFEM_TETRA4;
92     }
93    
94     int FEM_TETRA4::get_dimension(void)
95     {
96     return 3;
97     }
98    
99    
100     int FEM_TETRA4::get_nb_fem_noeud(void)
101     {
102     return FEM_TEMPLATE_ELEMENT<4>::get_nb_fem_noeud();
103     }
104    
105     FEM_NOEUD* FEM_TETRA4::get_fem_noeud(int num)
106     {
107     return FEM_TEMPLATE_ELEMENT<4>::get_fem_noeud(num);
108     }
109    
110     void FEM_TETRA4::change_noeud(int num,FEM_NOEUD* noeud)
111     {
112     FEM_TEMPLATE_ELEMENT<4>::change_noeud(num,noeud);
113     }
114    
115    
116     BOITE_3D& FEM_TETRA4::get_boite_3D(void)
117     {
118     return FEM_TEMPLATE_ELEMENT<4>::get_boite_3D();
119     }
120     int FEM_TETRA4::nb_fonction_interpolation(void)
121     {
122     return 4;
123     }
124    
125    
126     double FEM_TETRA4::get_fonction_interpolation(int num,double *uv)
127     {
128     double val;
129     switch (num)
130     {
131     case 1:
132     val=1-uv[0]-uv[1]-uv[2];
133     break;
134     case 2:
135     val=uv[0];
136     break;
137     case 3:
138     val=uv[1];
139     break;
140     case 4:
141     val=uv[2];
142     break;
143     }
144     return val;
145     }
146    
147     double FEM_TETRA4::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
148     {
149     double val;
150    
151     switch (num)
152     {
153     case 1:
154     switch (num_variable)
155     {
156     case 1:
157     val=-1;
158     break;
159     case 2:
160     val=-1;
161     break;
162     case 3:
163     val=-1;
164     break;
165     } break;
166     case 2:
167     switch (num_variable)
168     {
169     case 1:
170     val=1;
171     break;
172     case 2:
173     val=0;
174     break;
175     case 3:
176     val=0;
177     break;
178     }break;
179     case 3:
180     switch (num_variable)
181     {
182     case 1:
183     val=0;
184     break;
185     case 2:
186     val=1;
187     break;
188     case 3:
189     val=0;
190     break;
191     }break;
192     case 4:
193     switch (num_variable)
194     {
195     case 1:
196     val=0;
197     break;
198     case 2:
199     val=0;
200     break;
201     case 3:
202     val=1;
203     break;
204     }break;
205    
206     }
207     return val;
208     }
209    
210     double FEM_TETRA4::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
211     {
212     li=3;
213     col=3;
214     double *xyz1=tab[0]->get_coord();
215     double *xyz2=tab[1]->get_coord();
216     double *xyz3=tab[2]->get_coord();
217     double *xyz4=tab[3]->get_coord();
218    
219     jac[0]=unite*(xyz2[0] - xyz1[0]) ;
220     jac[1]=unite*(xyz2[1] - xyz1[1]) ;
221     jac[2]=unite*(xyz2[2] - xyz1[2]) ;
222    
223     jac[3]=unite*(xyz3[0] - xyz1[0]) ;
224     jac[4]=unite*(xyz3[1] - xyz1[1]) ;
225     jac[5]=unite*(xyz3[2] - xyz1[2]) ;
226    
227     jac[6]=unite*(xyz4[0] - xyz1[0]) ;
228     jac[7]=unite*(xyz4[1] - xyz1[1]) ;
229     jac[8]=unite*(xyz4[2] - xyz1[2]) ;
230    
231     double SIX_V= jac[0]*(jac[4]*jac[8]-jac[7]*jac[5])+jac[1]*(jac[5]*jac[6]-jac[8]*jac[3])+jac[2]*(jac[3]*jac[7]-jac[6]*jac[4]);
232    
233     return SIX_V;
234     }
235    
236     void FEM_TETRA4::get_inverse_jacob(double* j,double *uv,double unite)
237     {
238     double *xyz1=tab[0]->get_coord();
239     double *xyz2=tab[1]->get_coord();
240     double *xyz3=tab[2]->get_coord();
241     double *xyz4=tab[3]->get_coord();
242    
243     double jac[9];
244    
245     jac[0*3+0]= unite*(xyz2[0] - xyz1[0]) ;
246     jac[0*3+1]= unite*(xyz2[1] - xyz1[1]) ;
247     jac[0*3+2]= unite*(xyz2[2] - xyz1[2]) ;
248    
249     jac[1*3+0]= unite*(xyz3[0] - xyz1[0]) ;
250     jac[1*3+1]= unite*(xyz3[1] - xyz1[1]) ;
251     jac[1*3+2]= unite*(xyz3[2] - xyz1[2]) ;
252    
253     jac[2*3+0]= unite*(xyz4[0] - xyz1[0]) ;
254     jac[2*3+1]= unite*(xyz4[1] - xyz1[1]) ;
255     jac[2*3+2]= unite*(xyz4[2] - xyz1[2]) ;
256    
257     double detj=(jac[0*3+0]*jac[1*3+1]*jac[2*3+2]-jac[0*3+0]*jac[1*3+2]*jac[2*3+1]-jac[1*3+0]*jac[0*3+1]*jac[2*3+2]+jac[1*3+0]*jac[0*3+2]*jac[2*3+1]+jac[2*3+0]*jac[0*3+1]*jac[1*3+2]-jac[2*3+0]*jac[0*3+2]*jac[1*3+1]);
258    
259     j[0*3+0] = (jac[1*3+1]*jac[2*3+2]-jac[1*3+2]*jac[2*3+1])/(detj);
260     j[0*3+1] = -(jac[0*3+1]*jac[2*3+2]-jac[0*3+2]*jac[2*3+1])/(detj);
261     j[0*3+2] =-(-jac[0*3+1]*jac[1*3+2]+jac[0*3+2]*jac[1*3+1])/(detj);
262    
263     j[1*3+0] = -(jac[1*3+0]*jac[2*3+2]-jac[1*3+2]*jac[2*3+0])/(detj);
264     j[1*3+1] = (jac[0*3+0]*jac[2*3+2]-jac[0*3+2]*jac[2*3+0])/(detj);
265     j[1*3+2] = -(jac[0*3+0]*jac[1*3+2]-jac[0*3+2]*jac[1*3+0])/(detj);
266    
267     j[2*3+0] =-(-jac[1*3+0]*jac[2*3+1]+jac[1*3+1]*jac[2*3+0])/(detj);
268     j[2*3+1] = -(jac[0*3+0]*jac[2*3+1]-jac[0*3+1]*jac[2*3+0])/(detj);
269     j[2*3+2] = (jac[0*3+0]*jac[1*3+1]-jac[0*3+1]*jac[1*3+0])/(detj);
270    
271    
272    
273     }
274    
275    
276     void FEM_TETRA4::enregistrer(std::ostream& o)
277     {
278     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_TETRA4($"<< 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;
279     else o << "%" << get_id() << "=FEM_TETRA4(NULL,$"<< maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << ");" << std::endl;
280    
281     }