ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/geometrie/src/xfem_tetra4.cpp
Revision: 404
Committed: Thu May 30 13:05:01 2013 UTC (11 years, 11 months ago) by francois
File size: 6649 byte(s)
Log Message: 
Enregistrement des entites 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 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 404 if (get_lien_topologie()!=NULL) 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     else o << "%" << get_id() << "=XFEM_TETRA4($"<< femmai->get_id() << ",NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << "," << etat << ");" << std::endl;
90 francois 340
91     }
92    
93    
94    
95    
96     int XFEM_TETRA4::nb_fonction_interpolation(void)
97     {
98     return 4;
99     }
100     double XFEM_TETRA4::get_fonction_interpolation(int num,double *uv)
101     {
102     double val;
103     switch (num)
104     {
105     case 1:
106     val=1-uv[0]-uv[1]-uv[2];
107     break;
108     case 2:
109     val=uv[0];
110     break;
111     case 3:
112     val=uv[1];
113     break;
114     case 4:
115     val=uv[2];
116     break;
117     }
118     return val;
119     }
120     double XFEM_TETRA4::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
121     {
122     double val;
123    
124     switch (num)
125     {
126     case 1:
127     switch (num_variable)
128     {
129     case 1:
130     val=-1;
131     break;
132     case 2:
133     val=-1;
134     break;
135     case 3:
136     val=-1;
137     break;
138     } break;
139     case 2:
140     switch (num_variable)
141     {
142     case 1:
143     val=1;
144     break;
145     case 2:
146     val=0;
147     break;
148     case 3:
149     val=0;
150     break;
151     }break;
152     case 3:
153     switch (num_variable)
154     {
155     case 1:
156     val=0;
157     break;
158     case 2:
159     val=1;
160     break;
161     case 3:
162     val=0;
163     break;
164     }break;
165     case 4:
166     switch (num_variable)
167     {
168     case 1:
169     val=0;
170     break;
171     case 2:
172     val=0;
173     break;
174     case 3:
175     val=1;
176     break;
177     }break;
178    
179     }
180     return val;
181     }
182     double XFEM_TETRA4::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
183     {
184     li=3;
185     col=3;
186     double *xyz1=tab[0]->get_coord();
187     double *xyz2=tab[1]->get_coord();
188     double *xyz3=tab[2]->get_coord();
189     double *xyz4=tab[3]->get_coord();
190    
191     jac[0]=unite*(xyz2[0] - xyz1[0]) ;
192     jac[1]=unite*(xyz2[1] - xyz1[1]) ;
193     jac[2]=unite*(xyz2[2] - xyz1[2]) ;
194    
195     jac[3]=unite*(xyz3[0] - xyz1[0]) ;
196     jac[4]=unite*(xyz3[1] - xyz1[1]) ;
197     jac[5]=unite*(xyz3[2] - xyz1[2]) ;
198    
199     jac[6]=unite*(xyz4[0] - xyz1[0]) ;
200     jac[7]=unite*(xyz4[1] - xyz1[1]) ;
201     jac[8]=unite*(xyz4[2] - xyz1[2]) ;
202    
203     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]);
204    
205     return SIX_V;
206     }
207     void XFEM_TETRA4::get_inverse_jacob(double* j,double *uv,double unite)
208     {
209     double *xyz1=tab[0]->get_coord();
210     double *xyz2=tab[1]->get_coord();
211     double *xyz3=tab[2]->get_coord();
212     double *xyz4=tab[3]->get_coord();
213    
214     double jac[9];
215    
216     jac[0*3+0]= unite*(xyz2[0] - xyz1[0]) ;
217     jac[0*3+1]= unite*(xyz2[1] - xyz1[1]) ;
218     jac[0*3+2]= unite*(xyz2[2] - xyz1[2]) ;
219    
220     jac[1*3+0]= unite*(xyz3[0] - xyz1[0]) ;
221     jac[1*3+1]= unite*(xyz3[1] - xyz1[1]) ;
222     jac[1*3+2]= unite*(xyz3[2] - xyz1[2]) ;
223    
224     jac[2*3+0]= unite*(xyz4[0] - xyz1[0]) ;
225     jac[2*3+1]= unite*(xyz4[1] - xyz1[1]) ;
226     jac[2*3+2]= unite*(xyz4[2] - xyz1[2]) ;
227    
228     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]);
229    
230     j[0*3+0] = (jac[1*3+1]*jac[2*3+2]-jac[1*3+2]*jac[2*3+1])/(detj);
231     j[0*3+1] = -(jac[0*3+1]*jac[2*3+2]-jac[0*3+2]*jac[2*3+1])/(detj);
232     j[0*3+2] =-(-jac[0*3+1]*jac[1*3+2]+jac[0*3+2]*jac[1*3+1])/(detj);
233    
234     j[1*3+0] = -(jac[1*3+0]*jac[2*3+2]-jac[1*3+2]*jac[2*3+0])/(detj);
235     j[1*3+1] = (jac[0*3+0]*jac[2*3+2]-jac[0*3+2]*jac[2*3+0])/(detj);
236     j[1*3+2] = -(jac[0*3+0]*jac[1*3+2]-jac[0*3+2]*jac[1*3+0])/(detj);
237    
238     j[2*3+0] =-(-jac[1*3+0]*jac[2*3+1]+jac[1*3+1]*jac[2*3+0])/(detj);
239     j[2*3+1] = -(jac[0*3+0]*jac[2*3+1]-jac[0*3+1]*jac[2*3+0])/(detj);
240     j[2*3+2] = (jac[0*3+0]*jac[1*3+1]-jac[0*3+1]*jac[1*3+0])/(detj);
241    
242     }