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