geometrie/src/fem_tetra10.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�re et Vincent FRANCOIS
//  D�partement de G�nie M�canique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�partement
//  de g�nie m�canique de l'Universit� du Qu�bec �
//  Trois Rivi�res
//  Les librairies ne peuvent �tre utilis�es sans l'accord
//  des auteurs (contact : francois@uqtr.ca)
//------------------------------------------------------------
//------------------------------------------------------------
//
//       fem_element310.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 � 11H22
//------------------------------------------------------------
//------------------------------------------------------------


#include "gestionversion.h"
#include "fem_tetra10.h"
#include "fem_noeud.h"
#include "mg_element_maillage.h"
#include "fem_maillage.h"
#include "ot_quadrature_gauss.h"


FEM_TETRA10::FEM_TETRA10(unsigned long num,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT3(num,mai),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}

FEM_TETRA10::FEM_TETRA10(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT3(mai),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
FEM_TETRA10::FEM_TETRA10(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class FEM_NOEUD** tabnoeud):FEM_ELEMENT3(num,topo),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}

FEM_TETRA10::FEM_TETRA10(class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT3(topo),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
FEM_TETRA10::FEM_TETRA10(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT3(num,topo,mai),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}

FEM_TETRA10::FEM_TETRA10(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT3(topo,mai),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}


FEM_TETRA10::FEM_TETRA10(FEM_TETRA10& mdd):FEM_ELEMENT3(mdd),FEM_TEMPLATE_ELEMENT<10>(mdd)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element3()->ajouter(this);
    tab[1]->get_lien_element3()->ajouter(this);
    tab[2]->get_lien_element3()->ajouter(this);
    tab[3]->get_lien_element3()->ajouter(this);
    tab[4]->get_lien_element3()->ajouter(this);
    tab[5]->get_lien_element3()->ajouter(this);
    tab[6]->get_lien_element3()->ajouter(this);
    tab[7]->get_lien_element3()->ajouter(this);
    tab[8]->get_lien_element3()->ajouter(this);
    tab[9]->get_lien_element3()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
}
FEM_TETRA10::~FEM_TETRA10()
{
    if (liaison_topologique!=NULL) 
      if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
    tab[0]->get_lien_element3()->supprimer(this);
    tab[1]->get_lien_element3()->supprimer(this);
    tab[2]->get_lien_element3()->supprimer(this);
    tab[3]->get_lien_element3()->supprimer(this);
    tab[4]->get_lien_element3()->supprimer(this);
    tab[5]->get_lien_element3()->supprimer(this);
    tab[6]->get_lien_element3()->supprimer(this);
    tab[7]->get_lien_element3()->supprimer(this);
    tab[8]->get_lien_element3()->supprimer(this);
    tab[9]->get_lien_element3()->supprimer(this);
    get_fem_noeudpetitid()->get_lien_petit_element3()->supprimer(this);
}

FEM_ELEMENT_MAILLAGE* FEM_TETRA10::dupliquer(FEM_MAILLAGE *femmai,long decalage)
{
    FEM_NOEUD* tabnoeud[10];
    tabnoeud[0]=femmai->get_fem_noeudid(tab[0]->get_id()+decalage);
    tabnoeud[1]=femmai->get_fem_noeudid(tab[1]->get_id()+decalage);
    tabnoeud[2]=femmai->get_fem_noeudid(tab[2]->get_id()+decalage);
    tabnoeud[3]=femmai->get_fem_noeudid(tab[3]->get_id()+decalage);
    tabnoeud[4]=femmai->get_fem_noeudid(tab[4]->get_id()+decalage);
    tabnoeud[5]=femmai->get_fem_noeudid(tab[5]->get_id()+decalage);
    tabnoeud[6]=femmai->get_fem_noeudid(tab[6]->get_id()+decalage);
    tabnoeud[7]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
    tabnoeud[8]=femmai->get_fem_noeudid(tab[8]->get_id()+decalage);
    tabnoeud[9]=femmai->get_fem_noeudid(tab[9]->get_id()+decalage);
    FEM_TETRA10* tet=new FEM_TETRA10(get_id()+decalage,maillage,tabnoeud);
    femmai->ajouter_fem_element3(tet);
    return tet;
}
void FEM_TETRA10::get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin)
{
voisin.vide();
if (no==tab[0])
      {
      voisin.ajouter(tab[1]);
      voisin.ajouter(tab[5]);
      voisin.ajouter(tab[6]);
      }
if (no==tab[1])
      {
      voisin.ajouter(tab[0]);
      voisin.ajouter(tab[2]);
      }
if (no==tab[2])
      {
      voisin.ajouter(tab[1]);
      voisin.ajouter(tab[3]);
      voisin.ajouter(tab[7]);
      }
if (no==tab[3])
      {
      voisin.ajouter(tab[2]);
      voisin.ajouter(tab[4]);
      }
if (no==tab[4])
      {
      voisin.ajouter(tab[3]);
      voisin.ajouter(tab[5]);
      voisin.ajouter(tab[8]);
      }
if (no==tab[5])
      {
      voisin.ajouter(tab[0]);
      voisin.ajouter(tab[4]);
      }
if (no==tab[6])
      {
      voisin.ajouter(tab[0]);
      voisin.ajouter(tab[9]);
      }
if (no==tab[7])
      {
      voisin.ajouter(tab[2]);
      voisin.ajouter(tab[9]);
      }
if (no==tab[8])
      {
      voisin.ajouter(tab[4]);
      voisin.ajouter(tab[9]);
      }
if (no==tab[9])
      {
      voisin.ajouter(tab[6]);
      voisin.ajouter(tab[7]);
      voisin.ajouter(tab[8]);
      }
}


void FEM_TETRA10::reinit_boite_3D(void)
{
FEM_TEMPLATE_ELEMENT<10>::reinit_boite_3D();  
}


int FEM_TETRA10::get_type_entite(void)
{
    return   IDFEM_TETRA10;
}

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


int FEM_TETRA10::get_nb_fem_noeud(void)
{
    return FEM_TEMPLATE_ELEMENT<10>::get_nb_fem_noeud();
}

FEM_NOEUD* FEM_TETRA10::get_fem_noeud(int num)
{
    return FEM_TEMPLATE_ELEMENT<10>::get_fem_noeud(num);
}

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

BOITE_3D& FEM_TETRA10::get_boite_3D(void)
{
    return FEM_TEMPLATE_ELEMENT<10>::get_boite_3D();
}


int FEM_TETRA10::get_nb_pt_gauss(int degre)
{
return OT_POINTS_GAUSS::get_nb_point_tetra(degre);
}
void FEM_TETRA10::get_pt_gauss(int degre,int num,double &w,double *uvw)
{
return OT_POINTS_GAUSS::get_pt_gauss_tet(degre,num,w,uvw);  
}

int FEM_TETRA10::get_degre_gauss(int num)
{
return OT_POINTS_GAUSS::get_degre_gauss_tetra(num);  
}

int FEM_TETRA10::get_nb_fonction_interpolation(void)
{
    return 10;
}

int FEM_TETRA10::get_degremax_fonction_interpolation(void)
{
     return 2;  
}
double FEM_TETRA10::get_fonction_interpolation(int num,double* uv)
{
    double val;
    double landa=1-uv[0]-uv[1]-uv[2];
    switch (num)
    {
    case 1:
        val=-landa*(1-2*landa);
        break;
    case 2:
        val=4*uv[0]*landa;
        break;
    case 3:
        val=-uv[0]*(1-2*uv[0]);
        break;
    case 4:
        val=4*uv[0]*uv[1];
        break;
    case 5:
        val=-uv[1]*(1-2*uv[1]);
        break;
    case 6:
        val=4*uv[1]*landa;
        break;
    case 7:
        val=4*uv[2]*landa;
        break;
    case 8:
        val=4*uv[0]*uv[2];
        break;
    case 9:
        val=4*uv[1]*uv[2];
        break;
    case 10:
        val=-uv[2]*(1-2*uv[2]);
        break;
    }
    return val;
}

double FEM_TETRA10::get_fonction_derive_interpolation(int num,int num_variable,double* uv)
{
    double val;
    double landa=1-uv[0]-uv[1]-uv[2];
    switch (num)
    {
    case 1:
        switch (num_variable)
        {
        case 1:
            val=1-4*landa;
            break;
        case 2:
            val=1-4*landa;
            break;
        case 3:
            val=1-4*landa;
            break;
        } break;
    case 2:
        switch (num_variable)
        {
        case 1:
            val=4*(landa-uv[0]);
            break;
        case 2:
            val=-4*uv[0];
            break;
        case 3:
            val=-4*uv[0];
            break;
        }break;
    case 3:
        switch (num_variable)
        {
        case 1:
            val=-1+4*uv[0];
            break;
        case 2:
            val=0.;
            break;
        case 3:
            val=0.;
            break;
        }break;
    case 4:
        switch (num_variable)
        {
        case 1:
            val=4*uv[1];
            break;
        case 2:
            val=4*uv[0];
            break;
        case 3:
            val=0.;
            break;
        }break;
    case 5:
        switch (num_variable)
        {
        case 1:
            val=0.;
            break;
        case 2:
            val=-1+4*uv[1];
            break;
        case 3:
            val=0.;
            break;
        }break;
    case 6:
        switch (num_variable)
        {
        case 1:
            val=-4*uv[1];
            break;
        case 2:
            val=4*(landa-uv[1]);
            break;
        case 3:
            val=-4*uv[1];
            break;
        }break;
    case 7:
        switch (num_variable)
        {
        case 1:
            val=-4*uv[2];
            break;
        case 2:
            val=-4*uv[2];
            break;
        case 3:
            val=4*(landa-uv[2]);
            break;
        }break;
    case 8:
        switch (num_variable)
        {
        case 1:
            val=4*uv[2];
            break;
        case 2:
            val=0;
            break;
        case 3:
            val=4*uv[0];
            break;
        }break;
    case 9:
        switch (num_variable)
        {
        case 1:
            val=0.;
            break;
        case 2:
            val=4*uv[2];
            break;
        case 3:
            val=4*uv[1];
            break;
        }break;
    case 10:
        switch (num_variable)
        {
        case 1:
            val=0.;
            break;
        case 2:
            val=0.;
            break;
        case 3:
            val=-1+4*uv[2];
            break;
        }break;
    }
    return val;

}

double FEM_TETRA10::get_jacobien(double* jac,double *uv,int& lig,int& col,double unite)
{
    lig=3;
    col=3;


    double nb_fct_int=this->get_nb_fonction_interpolation();

    double xyz[30];

    int nb_nd=get_nb_fem_noeud();

    for (int i=0;i<nb_nd;i++)
    {
        FEM_NOEUD* fmnd=get_fem_noeud(i);
        double *xyz1=fmnd->get_coord();

        xyz[3*i]=unite*xyz1[0];
        xyz[3*i+1]=unite*xyz1[1];
        xyz[3*i+2]=unite*xyz1[2];

    }

    for (int i=0;i<lig;i++)
    {
        for (int j=0;j<col;j++)
        {
            jac[i*col+j]=0.0;
            for (int k=0;k<nb_fct_int;k++)
            {
                double NDERI=this->get_fonction_derive_interpolation(k+1,i+1,uv);

                jac[i*col+j]+=NDERI*xyz[3*k+j];
            }
        }
    }

    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 FEM_TETRA10::get_inverse_jacob(double* j,double *uv,double unite)
{
    double jac[9];
    int lig,col;
    double detj=this->get_jacobien(jac,uv,lig,col,unite);


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


}

bool FEM_TETRA10::valide_parametre_EF(double* uvw)
{
if (uvw[0]>=-1e-10)
   if (uvw[1]>=-1e-10)
      if (uvw[2]>=-1e-10)
           if (uvw[0]+uvw[1]+uvw[2]<=1.+1e-10)  
          return true;
return false;
}


void FEM_TETRA10::enregistrer(std::ostream& o,double version)
{
    if (maillage!=NULL)
      if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_TETRA10($"<< get_lien_topologie()->get_id() << ",$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id() << ",$" << tab[7]->get_id() << ",$" << tab[8]->get_id() << ",$" << tab[9]->get_id()<<    ");" << std::endl;
      else o << "%" << get_id() << "=FEM_TETRA10(NULL,$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id() << ",$" << tab[7]->get_id() << ",$" << tab[8]->get_id() << ",$" << tab[9]->get_id()  << ");" << std::endl;
    else
      if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_TETRA10($"<< get_lien_topologie()->get_id() << ",NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id() << ",$" << tab[7]->get_id() << ",$" << tab[8]->get_id() << ",$" << tab[9]->get_id()<<    ");" << std::endl;
      else o << "%" << get_id() << "=FEM_TETRA10(NULL,NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id() << ",$" << tab[7]->get_id() << ",$" << tab[8]->get_id() << ",$" << tab[9]->get_id()  << ");" << std::endl;
      
}


int FEM_TETRA10::verifie_validite_decoupage_xfem(double *vol)
{
  
}
Revision:	786
Committed:	Wed Mar 9 21:10:49 2016 UTC (9 years, 2 months ago) by francois
File size:	18234 byte(s)
Log Message:	Possibilité d'obtenir la matrice de rigidite de code aster sous form de matrice bande. Pour cela ajouter -opt dans l'operation -fem à la création du maillage FEM et dans -calculaster utiliser les nouveaux types d'étude rigidite. Attention pour cela il faut une version de code aster modifiée De manière automatique elle peut etre optenu en faisant ./instal_aster 12.5mod Uniquement implantée avec la version 12.5 de code aster.
#	User	Rev	Content
1	francois	283	//------------------------------------------------------------
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	francois	309	// fem_element310.cpp
16	francois	283	//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20			// Version du 02/03/2006 � 11H22
21			//------------------------------------------------------------
22			//------------------------------------------------------------
23
24
25			#include "gestionversion.h"
26			#include "fem_tetra10.h"
27			#include "fem_noeud.h"
28			#include "mg_element_maillage.h"
29			#include "fem_maillage.h"
30	francois	754	#include "ot_quadrature_gauss.h"
31	francois	283
32
33	francois	309	FEM_TETRA10::FEM_TETRA10(unsigned long num,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT3(num,mai),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
34	francois	283	{
35			if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
36	francois	309	tab[0]->get_lien_element3()->ajouter(this);
37			tab[1]->get_lien_element3()->ajouter(this);
38			tab[2]->get_lien_element3()->ajouter(this);
39			tab[3]->get_lien_element3()->ajouter(this);
40			tab[4]->get_lien_element3()->ajouter(this);
41			tab[5]->get_lien_element3()->ajouter(this);
42			tab[6]->get_lien_element3()->ajouter(this);
43			tab[7]->get_lien_element3()->ajouter(this);
44			tab[8]->get_lien_element3()->ajouter(this);
45			tab[9]->get_lien_element3()->ajouter(this);
46			get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
47	francois	283	}
48
49	francois	309	FEM_TETRA10::FEM_TETRA10(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT3(mai),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
50	francois	283	{
51			if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
52	francois	309	tab[0]->get_lien_element3()->ajouter(this);
53			tab[1]->get_lien_element3()->ajouter(this);
54			tab[2]->get_lien_element3()->ajouter(this);
55			tab[3]->get_lien_element3()->ajouter(this);
56			tab[4]->get_lien_element3()->ajouter(this);
57			tab[5]->get_lien_element3()->ajouter(this);
58			tab[6]->get_lien_element3()->ajouter(this);
59			tab[7]->get_lien_element3()->ajouter(this);
60			tab[8]->get_lien_element3()->ajouter(this);
61			tab[9]->get_lien_element3()->ajouter(this);
62			get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
63	francois	283	}
64	francois	378	FEM_TETRA10::FEM_TETRA10(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class FEM_NOEUD** tabnoeud):FEM_ELEMENT3(num,topo),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
65			{
66			if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
67			tab[0]->get_lien_element3()->ajouter(this);
68			tab[1]->get_lien_element3()->ajouter(this);
69			tab[2]->get_lien_element3()->ajouter(this);
70			tab[3]->get_lien_element3()->ajouter(this);
71			tab[4]->get_lien_element3()->ajouter(this);
72			tab[5]->get_lien_element3()->ajouter(this);
73			tab[6]->get_lien_element3()->ajouter(this);
74			tab[7]->get_lien_element3()->ajouter(this);
75			tab[8]->get_lien_element3()->ajouter(this);
76			tab[9]->get_lien_element3()->ajouter(this);
77			get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
78			}
79	francois	283
80	francois	378	FEM_TETRA10::FEM_TETRA10(class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT3(topo),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
81			{
82			if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
83			tab[0]->get_lien_element3()->ajouter(this);
84			tab[1]->get_lien_element3()->ajouter(this);
85			tab[2]->get_lien_element3()->ajouter(this);
86			tab[3]->get_lien_element3()->ajouter(this);
87			tab[4]->get_lien_element3()->ajouter(this);
88			tab[5]->get_lien_element3()->ajouter(this);
89			tab[6]->get_lien_element3()->ajouter(this);
90			tab[7]->get_lien_element3()->ajouter(this);
91			tab[8]->get_lien_element3()->ajouter(this);
92			tab[9]->get_lien_element3()->ajouter(this);
93			get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
94			}
95			FEM_TETRA10::FEM_TETRA10(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT3(num,topo,mai),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
96			{
97			if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
98			tab[0]->get_lien_element3()->ajouter(this);
99			tab[1]->get_lien_element3()->ajouter(this);
100			tab[2]->get_lien_element3()->ajouter(this);
101			tab[3]->get_lien_element3()->ajouter(this);
102			tab[4]->get_lien_element3()->ajouter(this);
103			tab[5]->get_lien_element3()->ajouter(this);
104			tab[6]->get_lien_element3()->ajouter(this);
105			tab[7]->get_lien_element3()->ajouter(this);
106			tab[8]->get_lien_element3()->ajouter(this);
107			tab[9]->get_lien_element3()->ajouter(this);
108			get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
109			}
110
111			FEM_TETRA10::FEM_TETRA10(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT3(topo,mai),FEM_TEMPLATE_ELEMENT<10>(tabnoeud)
112			{
113			if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
114			tab[0]->get_lien_element3()->ajouter(this);
115			tab[1]->get_lien_element3()->ajouter(this);
116			tab[2]->get_lien_element3()->ajouter(this);
117			tab[3]->get_lien_element3()->ajouter(this);
118			tab[4]->get_lien_element3()->ajouter(this);
119			tab[5]->get_lien_element3()->ajouter(this);
120			tab[6]->get_lien_element3()->ajouter(this);
121			tab[7]->get_lien_element3()->ajouter(this);
122			tab[8]->get_lien_element3()->ajouter(this);
123			tab[9]->get_lien_element3()->ajouter(this);
124			get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
125			}
126
127
128
129
130
131
132
133	francois	309	FEM_TETRA10::FEM_TETRA10(FEM_TETRA10& mdd):FEM_ELEMENT3(mdd),FEM_TEMPLATE_ELEMENT<10>(mdd)
134	francois	283	{
135			if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
136	francois	309	tab[0]->get_lien_element3()->ajouter(this);
137			tab[1]->get_lien_element3()->ajouter(this);
138			tab[2]->get_lien_element3()->ajouter(this);
139			tab[3]->get_lien_element3()->ajouter(this);
140			tab[4]->get_lien_element3()->ajouter(this);
141			tab[5]->get_lien_element3()->ajouter(this);
142			tab[6]->get_lien_element3()->ajouter(this);
143			tab[7]->get_lien_element3()->ajouter(this);
144			tab[8]->get_lien_element3()->ajouter(this);
145			tab[9]->get_lien_element3()->ajouter(this);
146			get_fem_noeudpetitid()->get_lien_petit_element3()->ajouter(this);
147	francois	283	}
148			FEM_TETRA10::~FEM_TETRA10()
149			{
150	francois	663	if (liaison_topologique!=NULL)
151			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
152	francois	309	tab[0]->get_lien_element3()->supprimer(this);
153			tab[1]->get_lien_element3()->supprimer(this);
154			tab[2]->get_lien_element3()->supprimer(this);
155			tab[3]->get_lien_element3()->supprimer(this);
156			tab[4]->get_lien_element3()->supprimer(this);
157			tab[5]->get_lien_element3()->supprimer(this);
158			tab[6]->get_lien_element3()->supprimer(this);
159			tab[7]->get_lien_element3()->supprimer(this);
160			tab[8]->get_lien_element3()->supprimer(this);
161			tab[9]->get_lien_element3()->supprimer(this);
162			get_fem_noeudpetitid()->get_lien_petit_element3()->supprimer(this);
163	francois	283	}
164
165			FEM_ELEMENT_MAILLAGE* FEM_TETRA10::dupliquer(FEM_MAILLAGE *femmai,long decalage)
166			{
167			FEM_NOEUD* tabnoeud[10];
168			tabnoeud[0]=femmai->get_fem_noeudid(tab[0]->get_id()+decalage);
169			tabnoeud[1]=femmai->get_fem_noeudid(tab[1]->get_id()+decalage);
170			tabnoeud[2]=femmai->get_fem_noeudid(tab[2]->get_id()+decalage);
171			tabnoeud[3]=femmai->get_fem_noeudid(tab[3]->get_id()+decalage);
172			tabnoeud[4]=femmai->get_fem_noeudid(tab[4]->get_id()+decalage);
173			tabnoeud[5]=femmai->get_fem_noeudid(tab[5]->get_id()+decalage);
174			tabnoeud[6]=femmai->get_fem_noeudid(tab[6]->get_id()+decalage);
175			tabnoeud[7]=femmai->get_fem_noeudid(tab[7]->get_id()+decalage);
176			tabnoeud[8]=femmai->get_fem_noeudid(tab[8]->get_id()+decalage);
177			tabnoeud[9]=femmai->get_fem_noeudid(tab[9]->get_id()+decalage);
178			FEM_TETRA10* tet=new FEM_TETRA10(get_id()+decalage,maillage,tabnoeud);
179	francois	309	femmai->ajouter_fem_element3(tet);
180	francois	283	return tet;
181			}
182	francois	786	void FEM_TETRA10::get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin)
183			{
184			voisin.vide();
185			if (no==tab[0])
186			{
187			voisin.ajouter(tab[1]);
188			voisin.ajouter(tab[5]);
189			voisin.ajouter(tab[6]);
190			}
191			if (no==tab[1])
192			{
193			voisin.ajouter(tab[0]);
194			voisin.ajouter(tab[2]);
195			}
196			if (no==tab[2])
197			{
198			voisin.ajouter(tab[1]);
199			voisin.ajouter(tab[3]);
200			voisin.ajouter(tab[7]);
201			}
202			if (no==tab[3])
203			{
204			voisin.ajouter(tab[2]);
205			voisin.ajouter(tab[4]);
206			}
207			if (no==tab[4])
208			{
209			voisin.ajouter(tab[3]);
210			voisin.ajouter(tab[5]);
211			voisin.ajouter(tab[8]);
212			}
213			if (no==tab[5])
214			{
215			voisin.ajouter(tab[0]);
216			voisin.ajouter(tab[4]);
217			}
218			if (no==tab[6])
219			{
220			voisin.ajouter(tab[0]);
221			voisin.ajouter(tab[9]);
222			}
223			if (no==tab[7])
224			{
225			voisin.ajouter(tab[2]);
226			voisin.ajouter(tab[9]);
227			}
228			if (no==tab[8])
229			{
230			voisin.ajouter(tab[4]);
231			voisin.ajouter(tab[9]);
232			}
233			if (no==tab[9])
234			{
235			voisin.ajouter(tab[6]);
236			voisin.ajouter(tab[7]);
237			voisin.ajouter(tab[8]);
238			}
239			}
240	francois	283
241
242	francois	684	void FEM_TETRA10::reinit_boite_3D(void)
243			{
244			FEM_TEMPLATE_ELEMENT<10>::reinit_boite_3D();
245			}
246	francois	283
247
248			int FEM_TETRA10::get_type_entite(void)
249			{
250			return IDFEM_TETRA10;
251			}
252
253			int FEM_TETRA10::get_dimension(void)
254			{
255			return 3;
256			}
257
258
259			int FEM_TETRA10::get_nb_fem_noeud(void)
260			{
261			return FEM_TEMPLATE_ELEMENT<10>::get_nb_fem_noeud();
262			}
263
264			FEM_NOEUD* FEM_TETRA10::get_fem_noeud(int num)
265			{
266			return FEM_TEMPLATE_ELEMENT<10>::get_fem_noeud(num);
267			}
268
269			void FEM_TETRA10::change_noeud(int num,FEM_NOEUD* noeud)
270			{
271			FEM_TEMPLATE_ELEMENT<10>::change_noeud(num,noeud);
272			}
273
274			BOITE_3D& FEM_TETRA10::get_boite_3D(void)
275			{
276			return FEM_TEMPLATE_ELEMENT<10>::get_boite_3D();
277			}
278
279
280	francois	637	int FEM_TETRA10::get_nb_pt_gauss(int degre)
281			{
282	francois	754	return OT_POINTS_GAUSS::get_nb_point_tetra(degre);
283	francois	637	}
284			void FEM_TETRA10::get_pt_gauss(int degre,int num,double &w,double *uvw)
285			{
286	francois	754	return OT_POINTS_GAUSS::get_pt_gauss_tet(degre,num,w,uvw);
287	francois	637	}
288	francois	283
289	francois	757	int FEM_TETRA10::get_degre_gauss(int num)
290			{
291			return OT_POINTS_GAUSS::get_degre_gauss_tetra(num);
292			}
293	francois	637
294	francois	638	int FEM_TETRA10::get_nb_fonction_interpolation(void)
295	francois	283	{
296			return 10;
297			}
298
299	francois	638	int FEM_TETRA10::get_degremax_fonction_interpolation(void)
300			{
301			return 2;
302			}
303	francois	283	double FEM_TETRA10::get_fonction_interpolation(int num,double* uv)
304			{
305			double val;
306			double landa=1-uv[0]-uv[1]-uv[2];
307			switch (num)
308			{
309			case 1:
310			val=-landa(1-2landa);
311			break;
312			case 2:
313			val=4uv[0]landa;
314			break;
315			case 3:
316			val=-uv[0](1-2uv[0]);
317			break;
318			case 4:
319			val=4uv[0]uv[1];
320			break;
321			case 5:
322			val=-uv[1](1-2uv[1]);
323			break;
324			case 6:
325			val=4uv[1]landa;
326			break;
327			case 7:
328			val=4uv[2]landa;
329			break;
330			case 8:
331			val=4uv[0]uv[2];
332			break;
333			case 9:
334			val=4uv[1]uv[2];
335			break;
336			case 10:
337			val=-uv[2](1-2uv[2]);
338			break;
339			}
340			return val;
341			}
342
343			double FEM_TETRA10::get_fonction_derive_interpolation(int num,int num_variable,double* uv)
344			{
345			double val;
346			double landa=1-uv[0]-uv[1]-uv[2];
347			switch (num)
348			{
349			case 1:
350			switch (num_variable)
351			{
352			case 1:
353			val=1-4*landa;
354			break;
355			case 2:
356			val=1-4*landa;
357			break;
358			case 3:
359			val=1-4*landa;
360			break;
361			} break;
362			case 2:
363			switch (num_variable)
364			{
365			case 1:
366			val=4*(landa-uv[0]);
367			break;
368			case 2:
369			val=-4*uv[0];
370			break;
371			case 3:
372			val=-4*uv[0];
373			break;
374			}break;
375			case 3:
376			switch (num_variable)
377			{
378			case 1:
379			val=-1+4*uv[0];
380			break;
381			case 2:
382			val=0.;
383			break;
384			case 3:
385			val=0.;
386			break;
387			}break;
388			case 4:
389			switch (num_variable)
390			{
391			case 1:
392			val=4*uv[1];
393			break;
394			case 2:
395			val=4*uv[0];
396			break;
397			case 3:
398			val=0.;
399			break;
400			}break;
401			case 5:
402			switch (num_variable)
403			{
404			case 1:
405			val=0.;
406			break;
407			case 2:
408			val=-1+4*uv[1];
409			break;
410			case 3:
411			val=0.;
412			break;
413			}break;
414			case 6:
415			switch (num_variable)
416			{
417			case 1:
418			val=-4*uv[1];
419			break;
420			case 2:
421			val=4*(landa-uv[1]);
422			break;
423			case 3:
424			val=-4*uv[1];
425			break;
426			}break;
427			case 7:
428			switch (num_variable)
429			{
430			case 1:
431			val=-4*uv[2];
432			break;
433			case 2:
434			val=-4*uv[2];
435			break;
436			case 3:
437			val=4*(landa-uv[2]);
438			break;
439			}break;
440			case 8:
441			switch (num_variable)
442			{
443			case 1:
444			val=4*uv[2];
445			break;
446			case 2:
447			val=0;
448			break;
449			case 3:
450			val=4*uv[0];
451			break;
452			}break;
453			case 9:
454			switch (num_variable)
455			{
456			case 1:
457			val=0.;
458			break;
459			case 2:
460			val=4*uv[2];
461			break;
462			case 3:
463			val=4*uv[1];
464			break;
465			}break;
466			case 10:
467			switch (num_variable)
468			{
469			case 1:
470			val=0.;
471			break;
472			case 2:
473			val=0.;
474			break;
475			case 3:
476			val=-1+4*uv[2];
477			break;
478			}break;
479			}
480			return val;
481
482			}
483
484			double FEM_TETRA10::get_jacobien(double* jac,double *uv,int& lig,int& col,double unite)
485			{
486			lig=3;
487			col=3;
488
489
490	francois	638	double nb_fct_int=this->get_nb_fonction_interpolation();
491	francois	283
492			double xyz[30];
493
494			int nb_nd=get_nb_fem_noeud();
495
496			for (int i=0;i<nb_nd;i++)
497			{
498			FEM_NOEUD* fmnd=get_fem_noeud(i);
499			double *xyz1=fmnd->get_coord();
500
501			xyz[3i]=unitexyz1[0];
502			xyz[3i+1]=unitexyz1[1];
503			xyz[3i+2]=unitexyz1[2];
504
505			}
506
507			for (int i=0;i<lig;i++)
508			{
509			for (int j=0;j<col;j++)
510			{
511			jac[i*col+j]=0.0;
512			for (int k=0;k<nb_fct_int;k++)
513			{
514			double NDERI=this->get_fonction_derive_interpolation(k+1,i+1,uv);
515
516			jac[icol+j]+=NDERIxyz[3*k+j];
517			}
518			}
519			}
520
521			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]);
522
523			return SIX_V;
524
525			}
526
527
528			void FEM_TETRA10::get_inverse_jacob(double* j,double *uv,double unite)
529			{
530			double jac[9];
531			int lig,col;
532			double detj=this->get_jacobien(jac,uv,lig,col,unite);
533
534
535			j[03+0] = (jac[13+1]jac[23+2]-jac[13+2]jac[2*3+1])/(detj);
536			j[03+1] = -(jac[03+1]jac[23+2]-jac[03+2]jac[2*3+1])/(detj);
537			j[03+2] =-(-jac[03+1]jac[13+2]+jac[03+2]jac[1*3+1])/(detj);
538
539			j[13+0] = -(jac[13+0]jac[23+2]-jac[13+2]jac[2*3+0])/(detj);
540			j[13+1] = (jac[03+0]jac[23+2]-jac[03+2]jac[2*3+0])/(detj);
541			j[13+2] = -(jac[03+0]jac[13+2]-jac[03+2]jac[1*3+0])/(detj);
542
543			j[23+0] =-(-jac[13+0]jac[23+1]+jac[13+1]jac[2*3+0])/(detj);
544			j[23+1] = -(jac[03+0]jac[23+1]-jac[03+1]jac[2*3+0])/(detj);
545			j[23+2] = (jac[03+0]jac[13+1]-jac[03+1]jac[1*3+0])/(detj);
546
547
548			}
549
550	francois	635	bool FEM_TETRA10::valide_parametre_EF(double* uvw)
551			{
552			if (uvw[0]>=-1e-10)
553			if (uvw[1]>=-1e-10)
554			if (uvw[2]>=-1e-10)
555			if (uvw[0]+uvw[1]+uvw[2]<=1.+1e-10)
556			return true;
557			return false;
558			}
559	francois	283
560
561	francois	763	void FEM_TETRA10::enregistrer(std::ostream& o,double version)
562	francois	283	{
563	francois	378	if (maillage!=NULL)
564			if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_TETRA10($"<< get_lien_topologie()->get_id() << ",$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id() << ",$" << tab[7]->get_id() << ",$" << tab[8]->get_id() << ",$" << tab[9]->get_id()<< ");" << std::endl;
565			else o << "%" << get_id() << "=FEM_TETRA10(NULL,$" << maillage->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id() << ",$" << tab[7]->get_id() << ",$" << tab[8]->get_id() << ",$" << tab[9]->get_id() << ");" << std::endl;
566			else
567			if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_TETRA10($"<< get_lien_topologie()->get_id() << ",NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id() << ",$" << tab[7]->get_id() << ",$" << tab[8]->get_id() << ",$" << tab[9]->get_id()<< ");" << std::endl;
568			else o << "%" << get_id() << "=FEM_TETRA10(NULL,NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << ",$" << tab[3]->get_id() << ",$" << tab[4]->get_id() << ",$" << tab[5]->get_id() << ",$" << tab[6]->get_id() << ",$" << tab[7]->get_id() << ",$" << tab[8]->get_id() << ",$" << tab[9]->get_id() << ");" << std::endl;
569
570	francois	283	}
571	francois	405
572
573	francois	406	int FEM_TETRA10::verifie_validite_decoupage_xfem(double *vol)
574	francois	405	{
575
576			}