mailleur3d_peau.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mailleur3d_peau.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "mailleur3d_peau.h"
#include "fct_taille.h"
#include "mg_volume.h"
#include "m3d_triangle.h"
 
 
 
 
 
MAILLEUR3D_PEAU::MAILLEUR3D_PEAU(MG_MAILLAGE* mgmai,MG_GEOMETRIE *mggeo,MG_VOLUME* mgvol,FCT_TAILLE* fct_taille,bool save):MAILLEUR3D(mgmai,mggeo,fct_taille,save,mgvol,MAGIC::MAILLEURFRONTALETAT::TOUTDETRUIRE)
{
 
}
 
 
MAILLEUR3D_PEAU::~MAILLEUR3D_PEAU()
{
}
 
 
 
void MAILLEUR3D_PEAU::cree_front(MG_VOLUME* mgvol) 
{
bool bienori=orientation_peau();
char chaine[255];
strcpy(chaine,"       Bonne orientation de la peau");
if (!bienori) strcpy(chaine,"       Changement d'orientation de la peau");
affiche(chaine);
LISTE_MG_TRIANGLE::iterator it;
TPL_MAP_ENTITE<M3D_TRIANGLE*> lsttrim3d;
 
for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(it);tri!=NULL;tri=mg_maillage->get_suivant_triangle(it))
  {
  M3D_TRIANGLE* nvtri=new M3D_TRIANGLE(tri->get_id(),tri->get_lien_topologie(),tri->get_noeud1(),tri->get_noeud2(),tri->get_noeud3(),tri->get_segment1(),tri->get_segment2(),tri->get_segment3(),tri->get_origine());  
  nvtri->change_frontiere();
  //if (tri->get_lien_topologie()==NULL)
  if (tri->get_origine()==MAGIC::ORIGINE::MAILLEUR_AUTO)
    {
    /*nvtri->change_origine(MAGIC::ORIGINE::IMPOSE);
    nvtri->get_noeud1()->change_origine(MAGIC::ORIGINE::IMPOSE);;
    nvtri->get_noeud2()->change_origine(MAGIC::ORIGINE::IMPOSE);;
    nvtri->get_noeud3()->change_origine(MAGIC::ORIGINE::IMPOSE);;
    nvtri->get_segment1()->change_origine(MAGIC::ORIGINE::IMPOSE);;
    nvtri->get_segment2()->change_origine(MAGIC::ORIGINE::IMPOSE);;
    nvtri->get_segment3()->change_origine(MAGIC::ORIGINE::IMPOSE);;
    nvtri->get_segment1()->change_dimension_topo_null(2);
    nvtri->get_segment2()->change_dimension_topo_null(2);
    nvtri->get_segment3()->change_dimension_topo_null(2);*/
    nvtri->change_origine(MAGIC::ORIGINE::TRIANGULATION);
    nvtri->get_noeud1()->change_origine(MAGIC::ORIGINE::TRIANGULATION);;
    nvtri->get_noeud2()->change_origine(MAGIC::ORIGINE::TRIANGULATION);;
    nvtri->get_noeud3()->change_origine(MAGIC::ORIGINE::TRIANGULATION);;
    nvtri->get_segment1()->change_origine(MAGIC::ORIGINE::TRIANGULATION);;
    nvtri->get_segment2()->change_origine(MAGIC::ORIGINE::TRIANGULATION);;
    nvtri->get_segment3()->change_origine(MAGIC::ORIGINE::TRIANGULATION);;
    }
    if (tri->get_lien_topologie()==NULL)
    {
      nvtri->get_segment1()->change_dimension_topo_null(2);
      nvtri->get_segment2()->change_dimension_topo_null(2);
      nvtri->get_segment3()->change_dimension_topo_null(2);
    }
  lsttrim3d.ajouter(nvtri);
  }
 
TPL_MAP_ENTITE<M3D_TRIANGLE*>::ITERATEUR it2;
 
for (M3D_TRIANGLE* tri=lsttrim3d.get_premier(it2);tri!=NULL;tri=lsttrim3d.get_suivant(it2))
  {
  mg_maillage->supprimer_mg_triangleid(tri->get_id());
  mg_maillage->ajouter_mg_triangle(tri);  
  }
 
for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(it);tri!=NULL;tri=mg_maillage->get_suivant_triangle(it))
   if (bienori) ajouter_front_courant(MAGIC::MAILLEURFRONTALETAT::NONFORCE,tri->get_noeud1(),tri->get_noeud2(),tri->get_noeud3(),tri);  
   else ajouter_front_courant(MAGIC::MAILLEURFRONTALETAT::NONFORCE,tri->get_noeud1(),tri->get_noeud3(),tri->get_noeud2(),tri);
  
}
 
 
 
bool MAILLEUR3D_PEAU::orientation_peau(void)
{
LISTE_MG_TRIANGLE::iterator it,itp;
/*TPL_MAP_ENTITE<MG_NOEUD*> lstnoeud;
BOITE_3D boite=mg_maillage->get_premier_triangle(it)->get_noeud1()->get_boite_3D();
for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(it);tri!=NULL;tri=mg_maillage->get_suivant_triangle(it))
    {
    lstnoeud.ajouter(tri->get_noeud1());    
    lstnoeud.ajouter(tri->get_noeud2());    
    lstnoeud.ajouter(tri->get_noeud3());    
    boite=boite+tri->get_noeud1()->get_boite_3D();
    boite=boite+tri->get_noeud2()->get_boite_3D();
    boite=boite+tri->get_noeud3()->get_boite_3D();
    }
TPL_OCTREE<MG_TRIANGLE*,MG_NOEUD*> octree;
octree.initialiser(&lstnoeud,1,boite.get_xmin(),boite.get_ymin(),boite.get_zmin(),boite.get_xmax(),boite.get_ymax(),boite.get_zmax());
for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(it);tri!=NULL;tri=mg_maillage->get_suivant_triangle(it))
    octree.inserer(tri);
*/
 
for (MG_TRIANGLE* tri=mg_maillage->get_premier_triangle(itp);tri!=NULL;tri=mg_maillage->get_suivant_triangle(itp))
{
double *xyz1=tri->get_noeud1()->get_coord();
double *xyz2=tri->get_noeud2()->get_coord();
double *xyz3=tri->get_noeud3()->get_coord();
OT_VECTEUR_3D n1n2(xyz1,xyz2);
OT_VECTEUR_3D n1n3(xyz1,xyz3);
OT_VECTEUR_3D n=n1n3&n1n2;
n.norme();
double xyz[3];
xyz[0]=(xyz1[0]+xyz2[0]+xyz3[0])/3.;
xyz[1]=(xyz1[1]+xyz2[1]+xyz3[1])/3.;
xyz[2]=(xyz1[2]+xyz2[2]+xyz3[2])/3.;
double eps=1e-6;
bool pbdefinition=true;
int tpositif=0,tnegatif=0;
for (MG_TRIANGLE* tri2=mg_maillage->get_premier_triangle(it);tri2!=NULL;tri2=mg_maillage->get_suivant_triangle(it))
    {
    if (tri2==tri) continue;
    double *xyz1=tri2->get_noeud1()->get_coord();
    double *xyz2=tri2->get_noeud2()->get_coord();
    double *xyz3=tri2->get_noeud3()->get_coord();
    OT_VECTEUR_3D n1n2(xyz1,xyz2);
    OT_VECTEUR_3D n1n3(xyz1,xyz3);
    OT_VECTEUR_3D n2=n1n3&n1n2;
    n2.norme();
    if (fabs(n2*n)<1-eps)
        {
        OT_VECTEUR_3D vec3=(-n);
        OT_MATRICE_3D mat(n1n2,n1n3,vec3);
        OT_VECTEUR_3D b(xyz1,xyz);
        OT_MATRICE_3D mat1(b,n1n3,vec3);
        OT_MATRICE_3D mat2(n1n2,b,vec3);
        OT_MATRICE_3D mat3(n1n2,n1n3,b);
        double det=mat.get_determinant();
        double t1=mat1.get_determinant()/det;
        double t2=mat2.get_determinant()/det;
        double t=mat3.get_determinant()/det;
        if (t1>-eps)
            if (t2>-eps)
                if (t1+t2<1+eps)
                {
                if (t>0) tpositif++; else tnegatif++;    
                }
            
    
        }
    else 
        {
        OT_VECTEUR_3D n1g(xyz1,xyz);
        if (fabs(n1g*n2)<eps)
            {
            pbdefinition=false;
            break;
            }
        }
    }
    if (pbdefinition)
    {
        if (tnegatif%2==1) return true;
    return false;
    }
}
 
return true; // ne doit jamais passer ici
}