ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/geometrie/src/xfem_tetra4.cpp
Revision: 378
Committed: Tue Nov 27 20:54:58 2012 UTC (12 years, 5 months ago) by francois
File size: 6620 byte(s)
Log Message: 
ajout de methode de construction pour les fem elements compatible avec le decoupage XFEM. + bug de format de fichier dans les fem elements non visible en dehors du xfem

File Contents

# User Rev Content
1 francois 340 //------------------------------------------------------------
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_element34.cpp
16     //
17     //------------------------------------------------------------
18     //------------------------------------------------------------
19     // COPYRIGHT 2000
20     // Version du 02/03/2006 � 11H22
21     //------------------------------------------------------------
22     //------------------------------------------------------------
23    
24    
25     #include "gestionversion.h"
26     #include "xfem_tetra4.h"
27     #include "fem_tetra4.h"
28     #include "fem_noeud.h"
29     #include "mg_element_maillage.h"
30     #include "fem_maillage.h"
31    
32    
33     XFEM_TETRA4::XFEM_TETRA4(unsigned long num,class FEM_TETRA4* tet,class FEM_NOEUD** tabnoeud):XFEM_ELEMENT3(num,tet->get_mg_element_maillage()),FEM_TEMPLATE_ELEMENT<4>(tabnoeud),tetra(tet)
34     {
35     tetra->ajouter_xfem(this);
36     }
37    
38     XFEM_TETRA4::XFEM_TETRA4(class FEM_TETRA4* tet,FEM_NOEUD** tabnoeud):XFEM_ELEMENT3(tet->get_mg_element_maillage()),FEM_TEMPLATE_ELEMENT<4>(tabnoeud),tetra(tet)
39     {
40     tetra->ajouter_xfem(this);
41     }
42    
43     XFEM_TETRA4::XFEM_TETRA4(XFEM_TETRA4& mdd):XFEM_ELEMENT3(mdd),FEM_TEMPLATE_ELEMENT<4>(mdd),tetra(mdd.tetra),n1(mdd.n1),n2(mdd.n2),n3(mdd.n3),n4(mdd.n4)
44     {
45     }
46    
47     XFEM_TETRA4::~XFEM_TETRA4()
48     {
49     tetra->supprimer_xfem(this);
50     }
51    
52    
53    
54     int XFEM_TETRA4::get_type_entite(void)
55     {
56     return IDXFEM_TETRA4;
57     }
58    
59     int XFEM_TETRA4::get_dimension(void)
60     {
61     return 3;
62     }
63    
64    
65     int XFEM_TETRA4::get_nb_fem_noeud(void)
66     {
67     return FEM_TEMPLATE_ELEMENT<4>::get_nb_fem_noeud();
68     }
69    
70     FEM_NOEUD* XFEM_TETRA4::get_fem_noeud(int num)
71     {
72     return FEM_TEMPLATE_ELEMENT<4>::get_fem_noeud(num);
73     }
74    
75     void XFEM_TETRA4::change_noeud(int num,FEM_NOEUD* noeud)
76     {
77     FEM_TEMPLATE_ELEMENT<4>::change_noeud(num,noeud);
78     }
79    
80    
81     BOITE_3D& XFEM_TETRA4::get_boite_3D(void)
82     {
83     return FEM_TEMPLATE_ELEMENT<4>::get_boite_3D();
84     }
85    
86     void XFEM_TETRA4::enregistrer(std::ostream& o)
87     {
88 francois 378 o << "%" << get_id() << "=XFEM_TETRA4($"<< tetra->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << "," << etat << ");" << std::endl;
89 francois 340
90     }
91    
92     FEM_ELEMENT_MAILLAGE* XFEM_TETRA4::dupliquer(FEM_MAILLAGE *femmai,long decalage)
93     {
94     tab[0]=femmai->get_fem_noeudid(n1);
95     tab[1]=femmai->get_fem_noeudid(n2);
96     tab[2]=femmai->get_fem_noeudid(n3);
97     tab[3]=femmai->get_fem_noeudid(n4);
98     }
99    
100    
101    
102     int XFEM_TETRA4::nb_fonction_interpolation(void)
103     {
104     return 4;
105     }
106     double XFEM_TETRA4::get_fonction_interpolation(int num,double *uv)
107     {
108     double val;
109     switch (num)
110     {
111     case 1:
112     val=1-uv[0]-uv[1]-uv[2];
113     break;
114     case 2:
115     val=uv[0];
116     break;
117     case 3:
118     val=uv[1];
119     break;
120     case 4:
121     val=uv[2];
122     break;
123     }
124     return val;
125     }
126     double XFEM_TETRA4::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
127     {
128     double val;
129    
130     switch (num)
131     {
132     case 1:
133     switch (num_variable)
134     {
135     case 1:
136     val=-1;
137     break;
138     case 2:
139     val=-1;
140     break;
141     case 3:
142     val=-1;
143     break;
144     } break;
145     case 2:
146     switch (num_variable)
147     {
148     case 1:
149     val=1;
150     break;
151     case 2:
152     val=0;
153     break;
154     case 3:
155     val=0;
156     break;
157     }break;
158     case 3:
159     switch (num_variable)
160     {
161     case 1:
162     val=0;
163     break;
164     case 2:
165     val=1;
166     break;
167     case 3:
168     val=0;
169     break;
170     }break;
171     case 4:
172     switch (num_variable)
173     {
174     case 1:
175     val=0;
176     break;
177     case 2:
178     val=0;
179     break;
180     case 3:
181     val=1;
182     break;
183     }break;
184    
185     }
186     return val;
187     }
188     double XFEM_TETRA4::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
189     {
190     li=3;
191     col=3;
192     double *xyz1=tab[0]->get_coord();
193     double *xyz2=tab[1]->get_coord();
194     double *xyz3=tab[2]->get_coord();
195     double *xyz4=tab[3]->get_coord();
196    
197     jac[0]=unite*(xyz2[0] - xyz1[0]) ;
198     jac[1]=unite*(xyz2[1] - xyz1[1]) ;
199     jac[2]=unite*(xyz2[2] - xyz1[2]) ;
200    
201     jac[3]=unite*(xyz3[0] - xyz1[0]) ;
202     jac[4]=unite*(xyz3[1] - xyz1[1]) ;
203     jac[5]=unite*(xyz3[2] - xyz1[2]) ;
204    
205     jac[6]=unite*(xyz4[0] - xyz1[0]) ;
206     jac[7]=unite*(xyz4[1] - xyz1[1]) ;
207     jac[8]=unite*(xyz4[2] - xyz1[2]) ;
208    
209     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]);
210    
211     return SIX_V;
212     }
213     void XFEM_TETRA4::get_inverse_jacob(double* j,double *uv,double unite)
214     {
215     double *xyz1=tab[0]->get_coord();
216     double *xyz2=tab[1]->get_coord();
217     double *xyz3=tab[2]->get_coord();
218     double *xyz4=tab[3]->get_coord();
219    
220     double jac[9];
221    
222     jac[0*3+0]= unite*(xyz2[0] - xyz1[0]) ;
223     jac[0*3+1]= unite*(xyz2[1] - xyz1[1]) ;
224     jac[0*3+2]= unite*(xyz2[2] - xyz1[2]) ;
225    
226     jac[1*3+0]= unite*(xyz3[0] - xyz1[0]) ;
227     jac[1*3+1]= unite*(xyz3[1] - xyz1[1]) ;
228     jac[1*3+2]= unite*(xyz3[2] - xyz1[2]) ;
229    
230     jac[2*3+0]= unite*(xyz4[0] - xyz1[0]) ;
231     jac[2*3+1]= unite*(xyz4[1] - xyz1[1]) ;
232     jac[2*3+2]= unite*(xyz4[2] - xyz1[2]) ;
233    
234     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]);
235    
236     j[0*3+0] = (jac[1*3+1]*jac[2*3+2]-jac[1*3+2]*jac[2*3+1])/(detj);
237     j[0*3+1] = -(jac[0*3+1]*jac[2*3+2]-jac[0*3+2]*jac[2*3+1])/(detj);
238     j[0*3+2] =-(-jac[0*3+1]*jac[1*3+2]+jac[0*3+2]*jac[1*3+1])/(detj);
239    
240     j[1*3+0] = -(jac[1*3+0]*jac[2*3+2]-jac[1*3+2]*jac[2*3+0])/(detj);
241     j[1*3+1] = (jac[0*3+0]*jac[2*3+2]-jac[0*3+2]*jac[2*3+0])/(detj);
242     j[1*3+2] = -(jac[0*3+0]*jac[1*3+2]-jac[0*3+2]*jac[1*3+0])/(detj);
243    
244     j[2*3+0] =-(-jac[1*3+0]*jac[2*3+1]+jac[1*3+1]*jac[2*3+0])/(detj);
245     j[2*3+1] = -(jac[0*3+0]*jac[2*3+1]-jac[0*3+1]*jac[2*3+0])/(detj);
246     j[2*3+2] = (jac[0*3+0]*jac[1*3+1]-jac[0*3+1]*jac[1*3+0])/(detj);
247    
248     }