ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/geometrie/src/xfem_tetra4.cpp
Revision: 399
Committed: Mon May 27 15:27:41 2013 UTC (11 years, 11 months ago) by francois
File size: 6392 byte(s)
Log Message: 
Creation d'une structure XFEM et FEM complete avec des element0 et des xelement0

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 francois 399 XFEM_TETRA4::XFEM_TETRA4(unsigned long num,class FEM_ELEMENT_MAILLAGE *mai,class MG_ELEMENT_TOPOLOGIQUE *mgtopo,class FEM_NOEUD** tabnoeud):XFEM_ELEMENT3(num,mai,mgtopo),FEM_TEMPLATE_ELEMENT<4>(tabnoeud)
34 francois 340 {
35 francois 399 femmai->ajouter_xfem(this);
36 francois 340 }
37    
38 francois 399 XFEM_TETRA4::XFEM_TETRA4(class FEM_ELEMENT_MAILLAGE *mai,class MG_ELEMENT_TOPOLOGIQUE *mgtopo,FEM_NOEUD** tabnoeud):XFEM_ELEMENT3(mai,mgtopo),FEM_TEMPLATE_ELEMENT<4>(tabnoeud)
39 francois 340 {
40 francois 399 femmai->ajouter_xfem(this);
41 francois 340 }
42    
43 francois 399 XFEM_TETRA4::XFEM_TETRA4(XFEM_TETRA4& mdd):XFEM_ELEMENT3(mdd),FEM_TEMPLATE_ELEMENT<4>(mdd)
44 francois 340 {
45     }
46    
47     XFEM_TETRA4::~XFEM_TETRA4()
48     {
49 francois 399 femmai->supprimer_xfem(this);
50 francois 340 }
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 399 o << "%" << get_id() << "=XFEM_TETRA4($"<< femmai->get_id() << ",$" << topo->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    
93    
94    
95     int XFEM_TETRA4::nb_fonction_interpolation(void)
96     {
97     return 4;
98     }
99     double XFEM_TETRA4::get_fonction_interpolation(int num,double *uv)
100     {
101     double val;
102     switch (num)
103     {
104     case 1:
105     val=1-uv[0]-uv[1]-uv[2];
106     break;
107     case 2:
108     val=uv[0];
109     break;
110     case 3:
111     val=uv[1];
112     break;
113     case 4:
114     val=uv[2];
115     break;
116     }
117     return val;
118     }
119     double XFEM_TETRA4::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
120     {
121     double val;
122    
123     switch (num)
124     {
125     case 1:
126     switch (num_variable)
127     {
128     case 1:
129     val=-1;
130     break;
131     case 2:
132     val=-1;
133     break;
134     case 3:
135     val=-1;
136     break;
137     } break;
138     case 2:
139     switch (num_variable)
140     {
141     case 1:
142     val=1;
143     break;
144     case 2:
145     val=0;
146     break;
147     case 3:
148     val=0;
149     break;
150     }break;
151     case 3:
152     switch (num_variable)
153     {
154     case 1:
155     val=0;
156     break;
157     case 2:
158     val=1;
159     break;
160     case 3:
161     val=0;
162     break;
163     }break;
164     case 4:
165     switch (num_variable)
166     {
167     case 1:
168     val=0;
169     break;
170     case 2:
171     val=0;
172     break;
173     case 3:
174     val=1;
175     break;
176     }break;
177    
178     }
179     return val;
180     }
181     double XFEM_TETRA4::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
182     {
183     li=3;
184     col=3;
185     double *xyz1=tab[0]->get_coord();
186     double *xyz2=tab[1]->get_coord();
187     double *xyz3=tab[2]->get_coord();
188     double *xyz4=tab[3]->get_coord();
189    
190     jac[0]=unite*(xyz2[0] - xyz1[0]) ;
191     jac[1]=unite*(xyz2[1] - xyz1[1]) ;
192     jac[2]=unite*(xyz2[2] - xyz1[2]) ;
193    
194     jac[3]=unite*(xyz3[0] - xyz1[0]) ;
195     jac[4]=unite*(xyz3[1] - xyz1[1]) ;
196     jac[5]=unite*(xyz3[2] - xyz1[2]) ;
197    
198     jac[6]=unite*(xyz4[0] - xyz1[0]) ;
199     jac[7]=unite*(xyz4[1] - xyz1[1]) ;
200     jac[8]=unite*(xyz4[2] - xyz1[2]) ;
201    
202     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]);
203    
204     return SIX_V;
205     }
206     void XFEM_TETRA4::get_inverse_jacob(double* j,double *uv,double unite)
207     {
208     double *xyz1=tab[0]->get_coord();
209     double *xyz2=tab[1]->get_coord();
210     double *xyz3=tab[2]->get_coord();
211     double *xyz4=tab[3]->get_coord();
212    
213     double jac[9];
214    
215     jac[0*3+0]= unite*(xyz2[0] - xyz1[0]) ;
216     jac[0*3+1]= unite*(xyz2[1] - xyz1[1]) ;
217     jac[0*3+2]= unite*(xyz2[2] - xyz1[2]) ;
218    
219     jac[1*3+0]= unite*(xyz3[0] - xyz1[0]) ;
220     jac[1*3+1]= unite*(xyz3[1] - xyz1[1]) ;
221     jac[1*3+2]= unite*(xyz3[2] - xyz1[2]) ;
222    
223     jac[2*3+0]= unite*(xyz4[0] - xyz1[0]) ;
224     jac[2*3+1]= unite*(xyz4[1] - xyz1[1]) ;
225     jac[2*3+2]= unite*(xyz4[2] - xyz1[2]) ;
226    
227     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]);
228    
229     j[0*3+0] = (jac[1*3+1]*jac[2*3+2]-jac[1*3+2]*jac[2*3+1])/(detj);
230     j[0*3+1] = -(jac[0*3+1]*jac[2*3+2]-jac[0*3+2]*jac[2*3+1])/(detj);
231     j[0*3+2] =-(-jac[0*3+1]*jac[1*3+2]+jac[0*3+2]*jac[1*3+1])/(detj);
232    
233     j[1*3+0] = -(jac[1*3+0]*jac[2*3+2]-jac[1*3+2]*jac[2*3+0])/(detj);
234     j[1*3+1] = (jac[0*3+0]*jac[2*3+2]-jac[0*3+2]*jac[2*3+0])/(detj);
235     j[1*3+2] = -(jac[0*3+0]*jac[1*3+2]-jac[0*3+2]*jac[1*3+0])/(detj);
236    
237     j[2*3+0] =-(-jac[1*3+0]*jac[2*3+1]+jac[1*3+1]*jac[2*3+0])/(detj);
238     j[2*3+1] = -(jac[0*3+0]*jac[2*3+1]-jac[0*3+1]*jac[2*3+0])/(detj);
239     j[2*3+2] = (jac[0*3+0]*jac[1*3+1]-jac[0*3+1]*jac[1*3+0])/(detj);
240    
241     }