solveur/src/sl_tetra11.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       sl_tetra11.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"

#pragma hdrstop
#include <math.h>
#include <assert.h>
#include "tpl_map_entite.h"
#include "fem_triangle.h"
#include "fem_segment.h"

#include "sl_tetra11.h"
#include "fem_tetra.h"
#include "fem_noeud.h"

#pragma package(smart_init)


SL_TETRA11::SL_TETRA11(FEM_TETRA* mgtet):SL_TETRA(mgtet)
{
}

SL_TETRA11::SL_TETRA11(SL_TETRA11& mdd):SL_TETRA(mdd)
{
}

SL_TETRA11::~SL_TETRA11()
{
}

int SL_TETRA11::nb_point_gauss(void)
{
return 1;
}


int SL_TETRA11::nb_fonction_interpolation(void)
{
return 4;
}

void SL_TETRA11::get_point_gauss(int num,double *xyz,double& poids)
{
xyz[0]=0.25;
xyz[1]=0.25;
xyz[2]=0.25;

poids=0.16666666666666666666666666666667;

}

double SL_TETRA11::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 SL_TETRA11::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 SL_TETRA11::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
{


li=3;
col=3;


FEM_NOEUD* noeud1=((FEM_TETRA*) mai)->get_fem_noeud(0);
FEM_NOEUD* noeud2=((FEM_TETRA*) mai)->get_fem_noeud(1);
FEM_NOEUD* noeud3=((FEM_TETRA*) mai)->get_fem_noeud(2);
FEM_NOEUD* noeud4=((FEM_TETRA*) mai)->get_fem_noeud(3);


   double *xyz1=noeud1->get_coord();
   double *xyz2=noeud2->get_coord();
   double *xyz3=noeud3->get_coord();
   double *xyz4=noeud4->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 SL_TETRA11::get_inverse_jacob(double* j,double *uv,double unite)
{

FEM_NOEUD* noeud1=((FEM_TETRA*) mai)->get_fem_noeud(0);
FEM_NOEUD* noeud2=((FEM_TETRA*) mai)->get_fem_noeud(1);
FEM_NOEUD* noeud3=((FEM_TETRA*) mai)->get_fem_noeud(2);
FEM_NOEUD* noeud4=((FEM_TETRA*) mai)->get_fem_noeud(3);

   double *xyz1=noeud1->get_coord();
   double *xyz2=noeud2->get_coord();
   double *xyz3=noeud3->get_coord();
   double *xyz4=noeud4->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);


}


int SL_TETRA11::nb_point_extrapolation()
{
return 4;
}


void SL_TETRA11::get_cord_ref_noauds(double* sommet)
{
sommet[0]=0;  sommet[1]=0;  sommet[2]=0;
sommet[3]=1;  sommet[4]=0;  sommet[5]=0;
sommet[6]=0;  sommet[7]=1;  sommet[8]=0;
sommet[9]=0;  sommet[10]=0; sommet[11]=1;

}

Revision:	1158
Committed:	Thu Jun 13 22:18:49 2024 UTC (11 months ago) by francois
File size:	6245 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			//####// sl_tetra11.cpp
15			//####//
16			//####//------------------------------------------------------------
17			//####//------------------------------------------------------------
18			//####// COPYRIGHT 2000-2024
19			//####// jeu 13 jun 2024 11:58:57 EDT
20			//####//------------------------------------------------------------
21			//####//------------------------------------------------------------
22		5	#include "gestionversion.h"
23
24			#pragma hdrstop
25			#include <math.h>
26			#include <assert.h>
27			#include "tpl_map_entite.h"
28			#include "fem_triangle.h"
29			#include "fem_segment.h"
30
31			#include "sl_tetra11.h"
32			#include "fem_tetra.h"
33			#include "fem_noeud.h"
34
35			#pragma package(smart_init)
36
37
38			SL_TETRA11::SL_TETRA11(FEM_TETRA* mgtet):SL_TETRA(mgtet)
39			{
40			}
41
42			SL_TETRA11::SL_TETRA11(SL_TETRA11& mdd):SL_TETRA(mdd)
43			{
44			}
45
46			SL_TETRA11::~SL_TETRA11()
47			{
48			}
49
50			int SL_TETRA11::nb_point_gauss(void)
51			{
52			return 1;
53			}
54
55
56			int SL_TETRA11::nb_fonction_interpolation(void)
57			{
58			return 4;
59			}
60
61			void SL_TETRA11::get_point_gauss(int num,double *xyz,double& poids)
62			{
63			xyz[0]=0.25;
64			xyz[1]=0.25;
65			xyz[2]=0.25;
66
67			poids=0.16666666666666666666666666666667;
68
69			}
70
71			double SL_TETRA11::get_fonction_interpolation(int num,double *uv)
72			{
73			double val;
74
75			switch (num)
76			{
77			case 1:val=1-uv[0]-uv[1]-uv[2];break;
78			case 2:val=uv[0];break;
79			case 3:val=uv[1];break;
80			case 4:val=uv[2];break;
81
82			}
83			return val;
84			}
85
86			double SL_TETRA11::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
87			{
88			double val;
89
90			switch (num)
91			{
92			case 1:switch (num_variable)
93			{
94			case 1:val=-1;break;
95			case 2:val=-1;break;
96			case 3:val=-1;break;
97			} break;
98			case 2:switch (num_variable)
99			{
100			case 1:val=1;break;
101			case 2:val=0;break;
102			case 3:val=0;break;
103			}break;
104			case 3:switch (num_variable)
105			{
106			case 1:val=0;break;
107			case 2:val=1;break;
108			case 3:val=0;break;
109			}break;
110			case 4:switch (num_variable)
111			{
112			case 1:val=0;break;
113			case 2:val=0;break;
114			case 3:val=1;break;
115			}break;
116
117			}
118			return val;
119
120			}
121
122			double SL_TETRA11::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
123			{
124
125
126			li=3;
127			col=3;
128
129
130			FEM_NOEUD* noeud1=((FEM_TETRA*) mai)->get_fem_noeud(0);
131			FEM_NOEUD* noeud2=((FEM_TETRA*) mai)->get_fem_noeud(1);
132			FEM_NOEUD* noeud3=((FEM_TETRA*) mai)->get_fem_noeud(2);
133			FEM_NOEUD* noeud4=((FEM_TETRA*) mai)->get_fem_noeud(3);
134
135
136			double *xyz1=noeud1->get_coord();
137			double *xyz2=noeud2->get_coord();
138			double *xyz3=noeud3->get_coord();
139			double *xyz4=noeud4->get_coord();
140
141			jac[0]=unite*(xyz2[0] - xyz1[0]) ;
142			jac[1]=unite*(xyz2[1] - xyz1[1]) ;
143			jac[2]=unite*(xyz2[2] - xyz1[2]) ;
144
145			jac[3]=unite*(xyz3[0] - xyz1[0]) ;
146			jac[4]=unite*(xyz3[1] - xyz1[1]) ;
147			jac[5]=unite*(xyz3[2] - xyz1[2]) ;
148
149			jac[6]=unite*(xyz4[0] - xyz1[0]) ;
150			jac[7]=unite*(xyz4[1] - xyz1[1]) ;
151			jac[8]=unite*(xyz4[2] - xyz1[2]) ;
152
153			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]);
154
155			return SIX_V;
156			}
157
158			void SL_TETRA11::get_inverse_jacob(double* j,double *uv,double unite)
159			{
160
161			FEM_NOEUD* noeud1=((FEM_TETRA*) mai)->get_fem_noeud(0);
162			FEM_NOEUD* noeud2=((FEM_TETRA*) mai)->get_fem_noeud(1);
163			FEM_NOEUD* noeud3=((FEM_TETRA*) mai)->get_fem_noeud(2);
164			FEM_NOEUD* noeud4=((FEM_TETRA*) mai)->get_fem_noeud(3);
165
166			double *xyz1=noeud1->get_coord();
167			double *xyz2=noeud2->get_coord();
168			double *xyz3=noeud3->get_coord();
169			double *xyz4=noeud4->get_coord();
170
171			double jac[9];
172
173			jac[03+0]= unite(xyz2[0] - xyz1[0]) ;
174			jac[03+1]= unite(xyz2[1] - xyz1[1]) ;
175			jac[03+2]= unite(xyz2[2] - xyz1[2]) ;
176
177			jac[13+0]= unite(xyz3[0] - xyz1[0]) ;
178			jac[13+1]= unite(xyz3[1] - xyz1[1]) ;
179			jac[13+2]= unite(xyz3[2] - xyz1[2]) ;
180
181			jac[23+0]= unite(xyz4[0] - xyz1[0]) ;
182			jac[23+1]= unite(xyz4[1] - xyz1[1]) ;
183			jac[23+2]= unite(xyz4[2] - xyz1[2]) ;
184
185			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]);
186
187			j[03+0] = (jac[13+1]jac[23+2]-jac[13+2]jac[2*3+1])/(detj);
188			j[03+1] = -(jac[03+1]jac[23+2]-jac[03+2]jac[2*3+1])/(detj);
189			j[03+2] =-(-jac[03+1]jac[13+2]+jac[03+2]jac[1*3+1])/(detj);
190
191			j[13+0] = -(jac[13+0]jac[23+2]-jac[13+2]jac[2*3+0])/(detj);
192			j[13+1] = (jac[03+0]jac[23+2]-jac[03+2]jac[2*3+0])/(detj);
193			j[13+2] = -(jac[03+0]jac[13+2]-jac[03+2]jac[1*3+0])/(detj);
194
195			j[23+0] =-(-jac[13+0]jac[23+1]+jac[13+1]jac[2*3+0])/(detj);
196			j[23+1] = -(jac[03+0]jac[23+1]-jac[03+1]jac[2*3+0])/(detj);
197			j[23+2] = (jac[03+0]jac[13+1]-jac[03+1]jac[1*3+0])/(detj);
198
199
200
201			}
202
203
204
205			int SL_TETRA11::nb_point_extrapolation()
206			{
207			return 4;
208			}
209
210
211			void SL_TETRA11::get_cord_ref_noauds(double* sommet)
212			{
213			sommet[0]=0; sommet[1]=0; sommet[2]=0;
214			sommet[3]=1; sommet[4]=0; sommet[5]=0;
215			sommet[6]=0; sommet[7]=1; sommet[8]=0;
216			sommet[9]=0; sommet[10]=0; sommet[11]=1;
217
218			}
219