fem_triangle6.cpp

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fem_triangle6.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include "fem_triangle6.h"
#include "fem_noeud.h"
#include "mg_element_maillage.h"
#include "fem_maillage.h"
#include "math.h"
#include "ot_quadrature_gauss.h"
 
FEM_TRIANGLE6::FEM_TRIANGLE6(unsigned long num,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,mai),FEM_TEMPLATE_ELEMENT<6>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element2()->ajouter(this);
    tab[1]->get_lien_element2()->ajouter(this);
    tab[2]->get_lien_element2()->ajouter(this);
    tab[3]->get_lien_element2()->ajouter(this);
    tab[4]->get_lien_element2()->ajouter(this);
    tab[5]->get_lien_element2()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
}
 
FEM_TRIANGLE6::FEM_TRIANGLE6(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(mai),FEM_TEMPLATE_ELEMENT<6>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element2()->ajouter(this);
    tab[1]->get_lien_element2()->ajouter(this);
    tab[2]->get_lien_element2()->ajouter(this);
    tab[3]->get_lien_element2()->ajouter(this);
    tab[4]->get_lien_element2()->ajouter(this);
    tab[5]->get_lien_element2()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
}
 
FEM_TRIANGLE6::FEM_TRIANGLE6(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,topo),FEM_TEMPLATE_ELEMENT<6>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element2()->ajouter(this);
    tab[1]->get_lien_element2()->ajouter(this);
    tab[2]->get_lien_element2()->ajouter(this);
    tab[3]->get_lien_element2()->ajouter(this);
    tab[4]->get_lien_element2()->ajouter(this);
    tab[5]->get_lien_element2()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
}
 
FEM_TRIANGLE6::FEM_TRIANGLE6(class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(topo),FEM_TEMPLATE_ELEMENT<6>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element2()->ajouter(this);
    tab[1]->get_lien_element2()->ajouter(this);
    tab[2]->get_lien_element2()->ajouter(this);
    tab[3]->get_lien_element2()->ajouter(this);
    tab[4]->get_lien_element2()->ajouter(this);
    tab[5]->get_lien_element2()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
}
FEM_TRIANGLE6::FEM_TRIANGLE6(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,class FEM_NOEUD** tabnoeud):FEM_ELEMENT2(num,topo,mai),FEM_TEMPLATE_ELEMENT<6>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element2()->ajouter(this);
    tab[1]->get_lien_element2()->ajouter(this);
    tab[2]->get_lien_element2()->ajouter(this);
    tab[3]->get_lien_element2()->ajouter(this);
    tab[4]->get_lien_element2()->ajouter(this);
    tab[5]->get_lien_element2()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
}
 
FEM_TRIANGLE6::FEM_TRIANGLE6(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT2(topo,mai),FEM_TEMPLATE_ELEMENT<6>(tabnoeud)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element2()->ajouter(this);
    tab[1]->get_lien_element2()->ajouter(this);
    tab[2]->get_lien_element2()->ajouter(this);
    tab[3]->get_lien_element2()->ajouter(this);
    tab[4]->get_lien_element2()->ajouter(this);
    tab[5]->get_lien_element2()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
}
 
 
 
void FEM_TRIANGLE6::reinit_boite_3D(void)
{
FEM_TEMPLATE_ELEMENT<6>::reinit_boite_3D();  
}
 
 
 
 
FEM_TRIANGLE6::FEM_TRIANGLE6(FEM_TRIANGLE6& mdd):FEM_ELEMENT2(mdd),FEM_TEMPLATE_ELEMENT<6>(mdd)
{
    if (liaison_topologique!=NULL) if (liaison_topologique->get_dimension()==2) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
    tab[0]->get_lien_element2()->ajouter(this);
    tab[1]->get_lien_element2()->ajouter(this);
    tab[2]->get_lien_element2()->ajouter(this);
    tab[3]->get_lien_element2()->ajouter(this);
    tab[4]->get_lien_element2()->ajouter(this);
    tab[5]->get_lien_element2()->ajouter(this);
    get_fem_noeudpetitid()->get_lien_petit_element2()->ajouter(this);
}
FEM_TRIANGLE6::~FEM_TRIANGLE6()
{
    if (liaison_topologique!=NULL) 
        if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
    tab[0]->get_lien_element2()->supprimer(this);
    tab[1]->get_lien_element2()->supprimer(this);
    tab[2]->get_lien_element2()->supprimer(this);
    tab[3]->get_lien_element2()->supprimer(this);
    tab[4]->get_lien_element2()->supprimer(this);
    tab[5]->get_lien_element2()->supprimer(this);
    get_fem_noeudpetitid()->get_lien_petit_element2()->supprimer(this);
}
 
FEM_ELEMENT_MAILLAGE* FEM_TRIANGLE6::dupliquer(FEM_MAILLAGE *femmai,long decalage)
{
    FEM_NOEUD* tabnoeud[6];
    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);
    FEM_TRIANGLE6* tri=new FEM_TRIANGLE6(get_id()+decalage,maillage,tabnoeud);
    femmai->ajouter_fem_element2(tri);
    return tri;
}
 
 
void FEM_TRIANGLE6::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]);
      }
if (no==tab[1])
      {
      voisin.ajouter(tab[0]);
      voisin.ajouter(tab[2]);
      }
if (no==tab[2])
      {
      voisin.ajouter(tab[1]);
      voisin.ajouter(tab[3]);
      }
if (no==tab[3])
      {
      voisin.ajouter(tab[2]);
      voisin.ajouter(tab[4]);
      }
if (no==tab[4])
      {
      voisin.ajouter(tab[3]);
      voisin.ajouter(tab[5]);
      }
if (no==tab[5])
      {
      voisin.ajouter(tab[4]);
      voisin.ajouter(tab[0]);
      }
}
bool FEM_TRIANGLE6::valide_parametre_EF(double* uvw)
{
if (uvw[0]>=-1e-10)
   if (uvw[1]>=-1e-10)
     if (uvw[0]+uvw[1]<=1.+1e-10)  
          return true;
return false;
}
 
 
int FEM_TRIANGLE6::get_type_entite(void)
{
    return   MAGIC::TYPE_ENTITE::IDFEM_TRIANGLE6;
}
 
int FEM_TRIANGLE6::get_dimension(void)
{
    return 2;
}
 
 
int FEM_TRIANGLE6::get_nb_fem_noeud(void)
{
    return FEM_TEMPLATE_ELEMENT<6>::get_nb_fem_noeud();
}
 
FEM_NOEUD* FEM_TRIANGLE6::get_fem_noeud(int num)
{
    return FEM_TEMPLATE_ELEMENT<6>::get_fem_noeud(num);
}
 
void FEM_TRIANGLE6::change_noeud(int num,FEM_NOEUD* noeud)
{
    FEM_TEMPLATE_ELEMENT<6>::change_noeud(num,noeud);
}
 
BOITE_3D& FEM_TRIANGLE6::get_boite_3D(void)
{
    return FEM_TEMPLATE_ELEMENT<6>::get_boite_3D();
}
 
 
void FEM_TRIANGLE6::enregistrer(std::ostream& o,double version)
{
    if (maillage!=NULL)
      if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_TRIANGLE6($"<< 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()  <<    ");" << std::endl;
      else o << "%" << get_id() << "=FEM_TRIANGLE6(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()  << ");" << std::endl;
    else
      if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_TRIANGLE6($"<< 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()  <<    ");" << std::endl;
      else o << "%" << get_id() << "=FEM_TRIANGLE6(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()  << ");" << std::endl;
 
}
 
int FEM_TRIANGLE6::get_nb_pt_gauss(int degre)
{
return OT_POINTS_GAUSS::get_nb_point_tri(degre);
}
 
void FEM_TRIANGLE6::get_pt_gauss(int degre,int num,double &w,double *uv)
{
return OT_POINTS_GAUSS::get_pt_gauss_tri(degre,num,w,uv);
}
 
int FEM_TRIANGLE6::get_degre_gauss(int num)
{
return OT_POINTS_GAUSS::get_degre_gauss_tri(num);  
}
 
int FEM_TRIANGLE6::get_nb_fonction_interpolation(void)
{
    return 6;
}
int FEM_TRIANGLE6::get_degremax_fonction_interpolation(void)
{
     return 2;  
}
double FEM_TRIANGLE6::get_fonction_interpolation(int num,double *uv)
{
    double val;
    switch (num)
    {
    case 1:
        val=(-1+uv[0]+uv[1])*(-1+2*uv[0]+2*uv[1]);
        break;
    case 2:
        val=4.*uv[0]*(1-uv[0]-uv[1]);
        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]*(1-uv[0]-uv[1]);
        break;
    
    }
    return val;
}
 
double FEM_TRIANGLE6::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
{
    double val;
 
    switch (num)
    {
    case 1:
        switch (num_variable)
        {
        case 1:
          val=-3.+4*uv[0]+4*uv[1];
            break;
        case 2:
            val=-3.+4*uv[0]+4*uv[1];
            break;
        } break;
    case 2:
        switch (num_variable)
        {
        case 1:
            val=4.-8*uv[0]-4*uv[1];
            break;
        case 2:
            val=-4.*uv[0];
            break;
        }break;
    case 3:
        switch (num_variable)
        {
        case 1:
            val=-1.+4*uv[0];
            break;
        case 2:
            val=0.;
            break;
        }break;
    case 4:
        switch (num_variable)
        {
        case 1:
            val=4*uv[1];
            break;
        case 2:
            val=4*uv[0];
            break;
        }break;
    case 5:
        switch (num_variable)
        {
        case 1:
            val=0.;
            break;
        case 2:
            val=-1+4.*uv[1];
            break;
        }break;
    case 6:
        switch (num_variable)
        {
        case 1:
            val=-4*uv[1];
            break;
        case 2:
            val=4-4*uv[0]-8*uv[1];
            break;
        }break;
    
    }
    return val;
}