ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/mtu/src/xfem_triangle3.cpp
Revision: 1158
Committed: Thu Jun 13 22:18:49 2024 UTC (11 months, 1 week ago) by francois
File size: 4192 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     //####// xfem_triangle3.cpp
15     //####//
16     //####//------------------------------------------------------------
17     //####//------------------------------------------------------------
18     //####// COPYRIGHT 2000-2024
19     //####// jeu 13 jun 2024 11:58:53 EDT
20     //####//------------------------------------------------------------
21     //####//------------------------------------------------------------
22 francois 410
23    
24     #include "gestionversion.h"
25     #include "xfem_triangle3.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     XFEM_TRIANGLE3::XFEM_TRIANGLE3(unsigned long num,class FEM_ELEMENT_MAILLAGE *mai,class MG_ELEMENT_TOPOLOGIQUE *mgtopo,class FEM_NOEUD** tabnoeud):XFEM_ELEMENT2(num,mai,mgtopo),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
33     {
34     femmai->ajouter_xfem(this);
35     }
36    
37     XFEM_TRIANGLE3::XFEM_TRIANGLE3(class FEM_ELEMENT_MAILLAGE *mai,class MG_ELEMENT_TOPOLOGIQUE *mgtopo,FEM_NOEUD** tabnoeud):XFEM_ELEMENT2(mai,mgtopo),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
38     {
39     femmai->ajouter_xfem(this);
40     }
41    
42     XFEM_TRIANGLE3::XFEM_TRIANGLE3(XFEM_TRIANGLE3& mdd):XFEM_ELEMENT2(mdd),FEM_TEMPLATE_ELEMENT<3>(mdd)
43     {
44     }
45    
46     XFEM_TRIANGLE3::~XFEM_TRIANGLE3()
47     {
48     femmai->supprimer_xfem(this);
49     }
50    
51    
52    
53     int XFEM_TRIANGLE3::get_type_entite(void)
54     {
55     return IDXFEM_TRIANGLE3;
56     }
57    
58     int XFEM_TRIANGLE3::get_dimension(void)
59     {
60     return 2;
61     }
62    
63    
64     int XFEM_TRIANGLE3::get_nb_fem_noeud(void)
65     {
66     return FEM_TEMPLATE_ELEMENT<3>::get_nb_fem_noeud();
67     }
68    
69     FEM_NOEUD* XFEM_TRIANGLE3::get_fem_noeud(int num)
70     {
71     return FEM_TEMPLATE_ELEMENT<3>::get_fem_noeud(num);
72     }
73    
74     void XFEM_TRIANGLE3::change_noeud(int num,FEM_NOEUD* noeud)
75     {
76     FEM_TEMPLATE_ELEMENT<3>::change_noeud(num,noeud);
77     }
78    
79    
80     BOITE_3D& XFEM_TRIANGLE3::get_boite_3D(void)
81     {
82     return FEM_TEMPLATE_ELEMENT<3>::get_boite_3D();
83     }
84    
85 francois 763 void XFEM_TRIANGLE3::enregistrer(std::ostream& o,double version)
86 francois 410 {
87     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=XFEM_TRIANGLE3($"<< femmai->get_id() << ",$" << topo->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << "," << etat << ");" << std::endl;
88     else o << "%" << get_id() << "=XFEM_TRIANGLE3($"<< femmai->get_id() << ",NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << "," << etat << ");" << std::endl;
89    
90     }
91    
92    
93    
94    
95     int XFEM_TRIANGLE3::nb_fonction_interpolation(void)
96     {
97     return 3;
98     }
99     double XFEM_TRIANGLE3::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];
106     break;
107     case 2:
108     val=uv[0];
109     break;
110     case 3:
111     val=uv[1];
112     break;
113     }
114     return val;
115     }
116     double XFEM_TRIANGLE3::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
117     {
118     double val;
119    
120     switch (num)
121     {
122     case 1:
123     switch (num_variable)
124     {
125     case 1:
126     val=-1;
127     break;
128     case 2:
129     val=-1;
130     break;
131     } break;
132     case 2:
133     switch (num_variable)
134     {
135     case 1:
136     val=1;
137     break;
138     case 2:
139     val=0;
140     break;
141     }break;
142     case 3:
143     switch (num_variable)
144     {
145     case 1:
146     val=0;
147     break;
148     case 2:
149     val=1;
150     break;
151     }break;
152    
153     }
154     return val;
155     }
156     double XFEM_TRIANGLE3::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
157     {
158 francois 1075 return 0.;
159 francois 410 }
160 francois 1075
161 francois 410 void XFEM_TRIANGLE3::get_inverse_jacob(double* j,double *uv,double unite)
162     {
163 francois 1075 }