mg_hexa.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mg_hexa.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:54 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include "mg_hexa.h"
#include "mg_maillage.h"
 
MG_HEXA::MG_HEXA(unsigned long num,MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD *mgnoeud1,class MG_NOEUD *mgnoeud2,class MG_NOEUD *mgnoeud3,class MG_NOEUD *mgnoeud4,class MG_NOEUD *mgnoeud5,class MG_NOEUD *mgnoeud6,class MG_NOEUD *mgnoeud7,class MG_NOEUD *mgnoeud8,class MG_QUADRANGLE* mgquadrangle1,class MG_QUADRANGLE* mgquadrangle2,class MG_QUADRANGLE* mgquadrangle3,class MG_QUADRANGLE* mgquadrangle4,class MG_QUADRANGLE* mgquadrangle5,class MG_QUADRANGLE* mgquadrangle6,int origine):MG_ELEMENT_MAILLAGE(num,topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),noeud6(mgnoeud6),noeud7(mgnoeud7),noeud8(mgnoeud8),quadrangle1(mgquadrangle1),quadrangle2(mgquadrangle2),quadrangle3(mgquadrangle3),quadrangle4(mgquadrangle4),quadrangle5(mgquadrangle5),quadrangle6(mgquadrangle6)
{
    init_hexa();
}
 
MG_HEXA::MG_HEXA(MG_ELEMENT_TOPOLOGIQUE* topo,class MG_NOEUD *mgnoeud1,class MG_NOEUD *mgnoeud2,class MG_NOEUD *mgnoeud3,class MG_NOEUD *mgnoeud4,class MG_NOEUD *mgnoeud5,class MG_NOEUD *mgnoeud6,class MG_NOEUD *mgnoeud7,class MG_NOEUD *mgnoeud8,class MG_QUADRANGLE* mgquadrangle1,class MG_QUADRANGLE* mgquadrangle2,class MG_QUADRANGLE* mgquadrangle3,class MG_QUADRANGLE* mgquadrangle4,class MG_QUADRANGLE* mgquadrangle5,class MG_QUADRANGLE* mgquadrangle6,int origine):MG_ELEMENT_MAILLAGE(topo,origine),noeud1(mgnoeud1),noeud2(mgnoeud2),noeud3(mgnoeud3),noeud4(mgnoeud4),noeud5(mgnoeud5),noeud6(mgnoeud6),noeud7(mgnoeud7),noeud8(mgnoeud8),quadrangle1(mgquadrangle1),quadrangle2(mgquadrangle2),quadrangle3(mgquadrangle3),quadrangle4(mgquadrangle4),quadrangle5(mgquadrangle5),quadrangle6(mgquadrangle6)
{
    init_hexa();
}
 
 
MG_HEXA::MG_HEXA(MG_HEXA& mdd):MG_ELEMENT_MAILLAGE(mdd),noeud1(mdd.noeud1),noeud2(mdd.noeud2),noeud3(mdd.noeud3),noeud4(mdd.noeud4),noeud5(mdd.noeud5),noeud6(mdd.noeud6),noeud7(mdd.noeud7),noeud8(mdd.noeud8),quadrangle1(mdd.quadrangle1),quadrangle2(mdd.quadrangle2),quadrangle3(mdd.quadrangle3),quadrangle4(mdd.quadrangle4),quadrangle5(mdd.quadrangle5),quadrangle6(mdd.quadrangle6)
{
    init_hexa();
}
 
 
MG_HEXA* MG_HEXA::dupliquer(MG_MAILLAGE *mgmai,long decalage)
{
    MG_NOEUD* nd1=mgmai->get_mg_noeudid(noeud1->get_id()+decalage);
    MG_NOEUD* nd2=mgmai->get_mg_noeudid(noeud2->get_id()+decalage);
    MG_NOEUD* nd3=mgmai->get_mg_noeudid(noeud3->get_id()+decalage);
    MG_NOEUD* nd4=mgmai->get_mg_noeudid(noeud4->get_id()+decalage);
    MG_NOEUD* nd5=mgmai->get_mg_noeudid(noeud5->get_id()+decalage);
    MG_NOEUD* nd6=mgmai->get_mg_noeudid(noeud6->get_id()+decalage);
    MG_NOEUD* nd7=mgmai->get_mg_noeudid(noeud7->get_id()+decalage);
    MG_NOEUD* nd8=mgmai->get_mg_noeudid(noeud8->get_id()+decalage);
    MG_QUADRANGLE* quad1=mgmai->get_mg_quadrangleid(quadrangle1->get_id()+decalage);
    MG_QUADRANGLE* quad2=mgmai->get_mg_quadrangleid(quadrangle2->get_id()+decalage);
    MG_QUADRANGLE* quad3=mgmai->get_mg_quadrangleid(quadrangle3->get_id()+decalage);
    MG_QUADRANGLE* quad4=mgmai->get_mg_quadrangleid(quadrangle4->get_id()+decalage);
    MG_QUADRANGLE* quad5=mgmai->get_mg_quadrangleid(quadrangle5->get_id()+decalage);
    MG_QUADRANGLE* quad6=mgmai->get_mg_quadrangleid(quadrangle6->get_id()+decalage);
    MG_HEXA* hex=new MG_HEXA(this->get_id()+decalage,liaison_topologique,nd1,nd2,nd3,nd4,nd5,nd6,nd7,nd8,quad1,quad2,quad3,quad4,quad5,quad6,MAGIC::ORIGINE::DUPLIQUER);
    mgmai->ajouter_mg_hexa(hex);
    return hex;
}
 
MG_HEXA::~MG_HEXA()
{
    if (noeud1!=NULL) noeud1->get_lien_hexa()->supprimer(this);
    if (noeud2!=NULL) noeud2->get_lien_hexa()->supprimer(this);
    if (noeud3!=NULL) noeud3->get_lien_hexa()->supprimer(this);
    if (noeud4!=NULL) noeud4->get_lien_hexa()->supprimer(this);
    if (noeud5!=NULL) noeud5->get_lien_hexa()->supprimer(this);
    if (noeud6!=NULL) noeud6->get_lien_hexa()->supprimer(this);
    if (noeud7!=NULL) noeud7->get_lien_hexa()->supprimer(this);
    if (noeud8!=NULL) noeud8->get_lien_hexa()->supprimer(this);
    unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
    noeudminid=std::min(noeudminid,noeud3->get_id());
    noeudminid=std::min(noeudminid,noeud4->get_id());
    noeudminid=std::min(noeudminid,noeud5->get_id());
    noeudminid=std::min(noeudminid,noeud6->get_id());
    noeudminid=std::min(noeudminid,noeud7->get_id());
    noeudminid=std::min(noeudminid,noeud8->get_id());
    if (noeud1!=NULL)
        if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_hexa()->supprimer(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_hexa()->supprimer(this);
    if (noeud3!=NULL)
        if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_hexa()->supprimer(this);
    if (noeud4!=NULL)
        if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_hexa()->supprimer(this);
    if (noeud5!=NULL)
        if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_hexa()->supprimer(this);
    if (noeud6!=NULL)
        if (noeud6->get_id()==noeudminid) noeud6->get_lien_petit_hexa()->supprimer(this);
    if (noeud7!=NULL)
        if (noeud7->get_id()==noeudminid) noeud7->get_lien_petit_hexa()->supprimer(this);
    if (noeud8!=NULL)
        if (noeud8->get_id()==noeudminid) noeud8->get_lien_petit_hexa()->supprimer(this);
    if (quadrangle1!=NULL) {
        quadrangle1->decrementer();
        quadrangle1->get_lien_hexa()->supprimer(this);
    }
    if (quadrangle2!=NULL) {
        quadrangle2->decrementer();
        quadrangle2->get_lien_hexa()->supprimer(this);
    }
    if (quadrangle3!=NULL) {
        quadrangle3->decrementer();
        quadrangle3->get_lien_hexa()->supprimer(this);
    }
    if (quadrangle4!=NULL) {
        quadrangle4->decrementer();
        quadrangle4->get_lien_hexa()->supprimer(this);
    }
    if (quadrangle5!=NULL) {
        quadrangle5->decrementer();
        quadrangle5->get_lien_hexa()->supprimer(this);
    }
    if (quadrangle6!=NULL) {
        quadrangle6->decrementer();
        quadrangle6->get_lien_hexa()->supprimer(this);
    }
    if (liaison_topologique!=NULL)
      if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->supprimer(this);
}
 
void MG_HEXA::init_hexa(void)
{
    if (noeud1!=NULL) noeud1->get_lien_hexa()->ajouter(this);
    if (noeud2!=NULL) noeud2->get_lien_hexa()->ajouter(this);
    if (noeud3!=NULL) noeud3->get_lien_hexa()->ajouter(this);
    if (noeud4!=NULL) noeud4->get_lien_hexa()->ajouter(this);
    if (noeud5!=NULL) noeud5->get_lien_hexa()->ajouter(this);
    if (noeud6!=NULL) noeud6->get_lien_hexa()->ajouter(this);
    if (noeud7!=NULL) noeud7->get_lien_hexa()->ajouter(this);
    if (noeud8!=NULL) noeud8->get_lien_hexa()->ajouter(this);
    unsigned long noeudminid=std::min(noeud1->get_id(),noeud2->get_id());
    noeudminid=std::min(noeudminid,noeud3->get_id());
    noeudminid=std::min(noeudminid,noeud4->get_id());
    noeudminid=std::min(noeudminid,noeud5->get_id());
    noeudminid=std::min(noeudminid,noeud6->get_id());
    noeudminid=std::min(noeudminid,noeud7->get_id());
    noeudminid=std::min(noeudminid,noeud8->get_id());
    if (noeud1!=NULL)
        if (noeud1->get_id()==noeudminid) noeud1->get_lien_petit_hexa()->ajouter(this);
    if (noeud2!=NULL)
        if (noeud2->get_id()==noeudminid) noeud2->get_lien_petit_hexa()->ajouter(this);
    if (noeud3!=NULL)
        if (noeud3->get_id()==noeudminid) noeud3->get_lien_petit_hexa()->ajouter(this);
    if (noeud4!=NULL)
        if (noeud4->get_id()==noeudminid) noeud4->get_lien_petit_hexa()->ajouter(this);
    if (noeud5!=NULL)
        if (noeud5->get_id()==noeudminid) noeud5->get_lien_petit_hexa()->ajouter(this);
    if (noeud6!=NULL)
        if (noeud6->get_id()==noeudminid) noeud6->get_lien_petit_hexa()->ajouter(this);
    if (noeud7!=NULL)
        if (noeud7->get_id()==noeudminid) noeud7->get_lien_petit_hexa()->ajouter(this);
    if (noeud8!=NULL)
        if (noeud8->get_id()==noeudminid) noeud8->get_lien_petit_hexa()->ajouter(this);
    if (quadrangle1!=NULL) {
        quadrangle1->incrementer();
        quadrangle1->get_lien_hexa()->ajouter(this);
    }
    if (quadrangle2!=NULL) {
        quadrangle2->incrementer();
        quadrangle2->get_lien_hexa()->ajouter(this);
    }
    if (quadrangle3!=NULL) {
        quadrangle3->incrementer();
        quadrangle3->get_lien_hexa()->ajouter(this);
    }
    if (quadrangle4!=NULL) {
        quadrangle4->incrementer();
        quadrangle4->get_lien_hexa()->ajouter(this);
    }
    if (quadrangle5!=NULL) {
        quadrangle5->incrementer();
        quadrangle5->get_lien_hexa()->ajouter(this);
    }
    if (quadrangle6!=NULL) {
        quadrangle6->incrementer();
        quadrangle6->get_lien_hexa()->ajouter(this);
    }
    
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==3) liaison_topologique->get_lien_maillage()->ajouter(this);
}
 
 
int MG_HEXA::get_type_entite(void)
{
    return MAGIC::TYPE_ENTITE::IDMG_HEXA;
}
 
 
 
 
 
MG_NOEUD* MG_HEXA::get_noeud1(void)
{
    return noeud1;
}
 
MG_NOEUD* MG_HEXA::get_noeud2(void)
{
    return noeud2;
}
 
MG_NOEUD* MG_HEXA::get_noeud3(void)
{
    return noeud3;
}
 
MG_NOEUD* MG_HEXA::get_noeud4(void)
{
    return noeud4;
}
MG_NOEUD* MG_HEXA::get_noeud5(void)
{
    return noeud5;
}
MG_NOEUD* MG_HEXA::get_noeud6(void)
{
    return noeud6;
}
MG_NOEUD* MG_HEXA::get_noeud7(void)
{
    return noeud7;
}
MG_NOEUD* MG_HEXA::get_noeud8(void)
{
    return noeud8;
}
 
MG_QUADRANGLE* MG_HEXA::get_quadrangle1(void)
{
    return quadrangle1;
}
 
MG_QUADRANGLE* MG_HEXA::get_quadrangle2(void)
{
    return quadrangle2;
}
 
MG_QUADRANGLE* MG_HEXA::get_quadrangle3(void)
{
    return quadrangle3;
}
 
MG_QUADRANGLE* MG_HEXA::get_quadrangle4(void)
{
    return quadrangle4;
}
 
MG_QUADRANGLE* MG_HEXA::get_quadrangle5(void)
{
    return quadrangle5;
}
 
MG_QUADRANGLE* MG_HEXA::get_quadrangle6(void)
{
    return quadrangle6;
}
 
void MG_HEXA::inverse_sens(void)
{
MG_NOEUD* noeud=noeud2;
noeud2=noeud4;
noeud4=noeud;
noeud=noeud6;
noeud6=noeud8;
noeud8=noeud;
MG_QUADRANGLE* quadrangle=quadrangle1;
quadrangle1=quadrangle4;
quadrangle4=quadrangle;
quadrangle=quadrangle2;
quadrangle2=quadrangle3;
quadrangle3=quadrangle;
}
 
 
int MG_HEXA::get_dimension(void)
{
    return 3;
}
 
 
BOITE_3D MG_HEXA::get_boite_3D(void)
{
    double coo1[3];
    double coo2[3];
    double coo3[3];
    double coo4[3];
    double coo5[3];
    double coo6[3];
    double coo7[3];
    double coo8[3];
    coo1[0]=noeud1->get_x();
    coo1[1]=noeud1->get_y();
    coo1[2]=noeud1->get_z();
    coo2[0]=noeud2->get_x();
    coo2[1]=noeud2->get_y();
    coo2[2]=noeud2->get_z();
    coo3[0]=noeud3->get_x();
    coo3[1]=noeud3->get_y();
    coo3[2]=noeud3->get_z();
    coo4[0]=noeud4->get_x();
    coo4[1]=noeud4->get_y();
    coo4[2]=noeud4->get_z();
    coo5[0]=noeud5->get_x();
    coo5[1]=noeud5->get_y();
    coo5[2]=noeud5->get_z();
    coo6[0]=noeud6->get_x();
    coo6[1]=noeud6->get_y();
    coo6[2]=noeud6->get_z();
    coo7[0]=noeud7->get_x();
    coo7[1]=noeud7->get_y();
    coo7[2]=noeud7->get_z();
    coo8[0]=noeud8->get_x();
    coo8[1]=noeud8->get_y();
    coo8[2]=noeud8->get_z();
    double xmin,ymin,zmin,xmax,ymax,zmax;
    xmin=std::min(coo1[0],coo2[0]);
    xmin=std::min(xmin,coo3[0]);
    xmin=std::min(xmin,coo4[0]);
    xmin=std::min(xmin,coo5[0]);
    xmin=std::min(xmin,coo6[0]);
    xmin=std::min(xmin,coo7[0]);
    xmin=std::min(xmin,coo8[0]);
    xmax=std::max(coo1[0],coo2[0]);
    xmax=std::max(xmax,coo3[0]);
    xmax=std::max(xmax,coo4[0]);
    xmax=std::max(xmax,coo5[0]);
    xmax=std::max(xmax,coo6[0]);
    xmax=std::max(xmax,coo7[0]);
    xmax=std::max(xmax,coo8[0]);
    ymin=std::min(coo1[1],coo2[1]);
    ymin=std::min(ymin,coo3[1]);
    ymin=std::min(ymin,coo4[1]);
    ymin=std::min(ymin,coo5[1]);
    ymin=std::min(ymin,coo6[1]);
    ymin=std::min(ymin,coo7[1]);
    ymin=std::min(ymin,coo8[1]);
    ymax=std::max(coo1[1],coo2[1]);
    ymax=std::max(ymax,coo3[1]);
    ymax=std::max(ymax,coo4[1]);
    ymax=std::max(ymax,coo5[1]);
    ymax=std::max(ymax,coo6[1]);
    ymax=std::max(ymax,coo7[1]);
    ymax=std::max(ymax,coo8[1]);
    zmin=std::min(coo1[2],coo2[2]);
    zmin=std::min(zmin,coo3[2]);
    zmin=std::min(zmin,coo4[2]);
    zmin=std::min(zmin,coo5[2]);
    zmin=std::min(zmin,coo6[2]);
    zmin=std::min(zmin,coo7[2]);
    zmin=std::min(zmin,coo8[2]);
    zmax=std::max(coo1[2],coo2[2]);
    zmax=std::max(zmax,coo3[2]);
    zmax=std::max(zmax,coo4[2]);
    zmax=std::max(zmax,coo5[2]);
    zmax=std::max(zmax,coo6[2]);
    zmax=std::max(zmax,coo7[2]);
    zmax=std::max(zmax,coo8[2]);
    BOITE_3D boite(xmin,ymin,zmin,xmax,ymax,zmax);
    return boite;
}
void MG_HEXA::extrapoler_solution_noeud(void)
{
    for (int j=0;j<MAX_TYPE_SOLUTION;j++)  
      {
        noeud1->change_solution(solution[j],j);
        noeud2->change_solution(solution[j],j);
        noeud3->change_solution(solution[j],j);
        noeud4->change_solution(solution[j],j);
        noeud5->change_solution(solution[j],j);
        noeud6->change_solution(solution[j],j);
        noeud7->change_solution(solution[j],j);
        noeud8->change_solution(solution[j],j);
      }
}
 
void MG_HEXA::enregistrer(std::ostream& o,double version)
{
    if ((liaison_topologique==NULL) && (MG_MAILLAGE::DIMENSIONMAILLAGESANSTOPO!=3)) return;
    if (liaison_topologique==NULL) o << "%" << get_id() << "=HEXAEDRE(NULL,$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id()  << ",$" << noeud4->get_id()  << ",$" << noeud5->get_id() << ",$" << noeud6->get_id() << ",$" << noeud7->get_id()  << ",$" << noeud8->get_id()  << "," << origine_creation << ");" << std::endl;
    else if (liaison_topologique->get_dimension()==3)  o << "%" << get_id() << "=HEXAEDRE($"<< get_lien_topologie()->get_id() << ",$" << noeud1->get_id() << ",$" << noeud2->get_id() << ",$" << noeud3->get_id()  << ",$" << noeud4->get_id()  << ",$" <<noeud5->get_id() << ",$" << noeud6->get_id() << ",$" << noeud7->get_id()  << ",$" << noeud8->get_id()  << "," << origine_creation << ");" << std::endl;
}