mtu/src/xfem_tetra4.cpp

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       xfem_tetra4.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------


#include "gestionversion.h"
#include "xfem_tetra4.h"
#include "fem_tetra4.h"
#include "fem_noeud.h"
#include "mg_element_maillage.h"
#include "fem_maillage.h"


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)
{
femmai->ajouter_xfem(this);
}

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)
{
femmai->ajouter_xfem(this);
}

XFEM_TETRA4::XFEM_TETRA4(XFEM_TETRA4& mdd):XFEM_ELEMENT3(mdd),FEM_TEMPLATE_ELEMENT<4>(mdd)
{
}

XFEM_TETRA4::~XFEM_TETRA4()
{
femmai->supprimer_xfem(this);
}


int XFEM_TETRA4::get_type_entite(void)
{
    return   IDXFEM_TETRA4;
}

int XFEM_TETRA4::get_dimension(void)
{
    return 3;
}


int XFEM_TETRA4::get_nb_fem_noeud(void)
{
    return FEM_TEMPLATE_ELEMENT<4>::get_nb_fem_noeud();
}

FEM_NOEUD* XFEM_TETRA4::get_fem_noeud(int num)
{
    return FEM_TEMPLATE_ELEMENT<4>::get_fem_noeud(num);
}

void XFEM_TETRA4::change_noeud(int num,FEM_NOEUD* noeud)
{
    FEM_TEMPLATE_ELEMENT<4>::change_noeud(num,noeud);
}


BOITE_3D& XFEM_TETRA4::get_boite_3D(void)
{
    return FEM_TEMPLATE_ELEMENT<4>::get_boite_3D();
}

void XFEM_TETRA4::enregistrer(std::ostream& o,double version)
{
   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;
    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;

}


int XFEM_TETRA4::nb_fonction_interpolation(void)
{
return 4;
}
double XFEM_TETRA4::get_fonction_interpolation(int num,double *uv)
{
double val;
    switch (num)
    {
    case 1:
        val=1-uv[0]-uv[1]-uv[2];
        break;
    case 2:
        val=uv[0];
        break;
    case 3:
        val=uv[1];
        break;
    case 4:
        val=uv[2];
        break;
    }
    return val;
}
double XFEM_TETRA4::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
{
double val;

    switch (num)
    {
    case 1:
        switch (num_variable)
        {
        case 1:
            val=-1;
            break;
        case 2:
            val=-1;
            break;
        case 3:
            val=-1;
            break;
        } break;
    case 2:
        switch (num_variable)
        {
        case 1:
            val=1;
            break;
        case 2:
            val=0;
            break;
        case 3:
            val=0;
            break;
        }break;
    case 3:
        switch (num_variable)
        {
        case 1:
            val=0;
            break;
        case 2:
            val=1;
            break;
        case 3:
            val=0;
            break;
        }break;
    case 4:
        switch (num_variable)
        {
        case 1:
            val=0;
            break;
        case 2:
            val=0;
            break;
        case 3:
            val=1;
            break;
        }break;

    }
    return val;
}
double XFEM_TETRA4::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
{
li=3;
    col=3;
    double *xyz1=tab[0]->get_coord();
    double *xyz2=tab[1]->get_coord();
    double *xyz3=tab[2]->get_coord();
    double *xyz4=tab[3]->get_coord();

    jac[0]=unite*(xyz2[0] - xyz1[0])  ;
    jac[1]=unite*(xyz2[1] - xyz1[1])  ;
    jac[2]=unite*(xyz2[2] - xyz1[2])  ;

    jac[3]=unite*(xyz3[0] - xyz1[0])  ;
    jac[4]=unite*(xyz3[1] - xyz1[1])  ;
    jac[5]=unite*(xyz3[2] - xyz1[2])  ;

    jac[6]=unite*(xyz4[0] - xyz1[0])  ;
    jac[7]=unite*(xyz4[1] - xyz1[1])  ;
    jac[8]=unite*(xyz4[2] - xyz1[2])  ;

    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]);

    return SIX_V;
}
void XFEM_TETRA4::get_inverse_jacob(double* j,double *uv,double unite)
{
double *xyz1=tab[0]->get_coord();
    double *xyz2=tab[1]->get_coord();
    double *xyz3=tab[2]->get_coord();
    double *xyz4=tab[3]->get_coord();

    double jac[9];

    jac[0*3+0]=  unite*(xyz2[0] - xyz1[0])  ;
    jac[0*3+1]=  unite*(xyz2[1] - xyz1[1])  ;
    jac[0*3+2]=  unite*(xyz2[2] - xyz1[2])  ;

    jac[1*3+0]=  unite*(xyz3[0] - xyz1[0])  ;
    jac[1*3+1]=  unite*(xyz3[1] - xyz1[1])  ;
    jac[1*3+2]=  unite*(xyz3[2] - xyz1[2])  ;

    jac[2*3+0]=  unite*(xyz4[0] - xyz1[0])  ;
    jac[2*3+1]=  unite*(xyz4[1] - xyz1[1])  ;
    jac[2*3+2]=  unite*(xyz4[2] - xyz1[2])  ;

    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]);

    j[0*3+0] =  (jac[1*3+1]*jac[2*3+2]-jac[1*3+2]*jac[2*3+1])/(detj);
    j[0*3+1] = -(jac[0*3+1]*jac[2*3+2]-jac[0*3+2]*jac[2*3+1])/(detj);
    j[0*3+2] =-(-jac[0*3+1]*jac[1*3+2]+jac[0*3+2]*jac[1*3+1])/(detj);

    j[1*3+0] = -(jac[1*3+0]*jac[2*3+2]-jac[1*3+2]*jac[2*3+0])/(detj);
    j[1*3+1] =  (jac[0*3+0]*jac[2*3+2]-jac[0*3+2]*jac[2*3+0])/(detj);
    j[1*3+2] = -(jac[0*3+0]*jac[1*3+2]-jac[0*3+2]*jac[1*3+0])/(detj);

    j[2*3+0] =-(-jac[1*3+0]*jac[2*3+1]+jac[1*3+1]*jac[2*3+0])/(detj);
    j[2*3+1] = -(jac[0*3+0]*jac[2*3+1]-jac[0*3+1]*jac[2*3+0])/(detj);
    j[2*3+2] =  (jac[0*3+0]*jac[1*3+1]-jac[0*3+1]*jac[1*3+0])/(detj);

}
Revision:	1158
Committed:	Thu Jun 13 22:18:49 2024 UTC (11 months, 1 week ago) by francois
File size:	6722 byte(s)
Log Message:	compatibilité Ubuntu 22.04 Suppression des refeences à Windows Ajout d'une banière
#	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_tetra4.cpp
15			//####//
16			//####//------------------------------------------------------------
17			//####//------------------------------------------------------------
18			//####// COPYRIGHT 2000-2024
19			//####// jeu 13 jun 2024 11:58:53 EDT
20			//####//------------------------------------------------------------
21			//####//------------------------------------------------------------
22	francois	340
23
24			#include "gestionversion.h"
25			#include "xfem_tetra4.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	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)
33	francois	340	{
34	francois	399	femmai->ajouter_xfem(this);
35	francois	340	}
36
37	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)
38	francois	340	{
39	francois	399	femmai->ajouter_xfem(this);
40	francois	340	}
41
42	francois	399	XFEM_TETRA4::XFEM_TETRA4(XFEM_TETRA4& mdd):XFEM_ELEMENT3(mdd),FEM_TEMPLATE_ELEMENT<4>(mdd)
43	francois	340	{
44			}
45
46			XFEM_TETRA4::~XFEM_TETRA4()
47			{
48	francois	399	femmai->supprimer_xfem(this);
49	francois	340	}
50
51
52
53			int XFEM_TETRA4::get_type_entite(void)
54			{
55			return IDXFEM_TETRA4;
56			}
57
58			int XFEM_TETRA4::get_dimension(void)
59			{
60			return 3;
61			}
62
63
64			int XFEM_TETRA4::get_nb_fem_noeud(void)
65			{
66			return FEM_TEMPLATE_ELEMENT<4>::get_nb_fem_noeud();
67			}
68
69			FEM_NOEUD* XFEM_TETRA4::get_fem_noeud(int num)
70			{
71			return FEM_TEMPLATE_ELEMENT<4>::get_fem_noeud(num);
72			}
73
74			void XFEM_TETRA4::change_noeud(int num,FEM_NOEUD* noeud)
75			{
76			FEM_TEMPLATE_ELEMENT<4>::change_noeud(num,noeud);
77			}
78
79
80			BOITE_3D& XFEM_TETRA4::get_boite_3D(void)
81			{
82			return FEM_TEMPLATE_ELEMENT<4>::get_boite_3D();
83			}
84
85	francois	763	void XFEM_TETRA4::enregistrer(std::ostream& o,double version)
86	francois	340	{
87	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;
88			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;
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[03+0]= unite(xyz2[0] - xyz1[0]) ;
216			jac[03+1]= unite(xyz2[1] - xyz1[1]) ;
217			jac[03+2]= unite(xyz2[2] - xyz1[2]) ;
218
219			jac[13+0]= unite(xyz3[0] - xyz1[0]) ;
220			jac[13+1]= unite(xyz3[1] - xyz1[1]) ;
221			jac[13+2]= unite(xyz3[2] - xyz1[2]) ;
222
223			jac[23+0]= unite(xyz4[0] - xyz1[0]) ;
224			jac[23+1]= unite(xyz4[1] - xyz1[1]) ;
225			jac[23+2]= unite(xyz4[2] - xyz1[2]) ;
226
227			double detj=(jac[03+0]jac[13+1]jac[23+2]-jac[03+0]jac[13+2]jac[23+1]-jac[13+0]jac[03+1]jac[23+2]+jac[13+0]jac[03+2]jac[23+1]+jac[23+0]jac[03+1]jac[13+2]-jac[23+0]jac[03+2]jac[13+1]);
228
229			j[03+0] = (jac[13+1]jac[23+2]-jac[13+2]jac[2*3+1])/(detj);
230			j[03+1] = -(jac[03+1]jac[23+2]-jac[03+2]jac[2*3+1])/(detj);
231			j[03+2] =-(-jac[03+1]jac[13+2]+jac[03+2]jac[1*3+1])/(detj);
232
233			j[13+0] = -(jac[13+0]jac[23+2]-jac[13+2]jac[2*3+0])/(detj);
234			j[13+1] = (jac[03+0]jac[23+2]-jac[03+2]jac[2*3+0])/(detj);
235			j[13+2] = -(jac[03+0]jac[13+2]-jac[03+2]jac[1*3+0])/(detj);
236
237			j[23+0] =-(-jac[13+0]jac[23+1]+jac[13+1]jac[2*3+0])/(detj);
238			j[23+1] = -(jac[03+0]jac[23+1]-jac[03+1]jac[2*3+0])/(detj);
239			j[23+2] = (jac[03+0]jac[13+1]-jac[03+1]jac[1*3+0])/(detj);
240
241	francois	1158	}