rec_face_plane.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       rec_face_plane.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:56 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
 
#pragma hdrstop
 
#include "occ_fonction.h"
#include "rec_face_plane.h"
#include "tpl_map_entite.h"
#include "tpl_liste_entite.h"
#include "rec_contour.h"
#include "fem_noeud.h"
#include "fem_maillage.h"
#include "ot_mathematique.h"
 
 
#pragma package(smart_init)
 
 
 
 
REC_FACE_PLANE::REC_FACE_PLANE(FEM_MAILLAGE* fem_mai, int Idealise, OCC_FONCTION_RECONSTRUCTION& fonc, int test, int* faces): REC_FACE(fem_mai,Idealise, fonc, test, faces)
{
 
}
 
REC_FACE_PLANE::~REC_FACE_PLANE()
{
 
}
 
void REC_FACE_PLANE::Reconstruire(MG_FACE* face, double coef)
{
        //OCC_FONCTION fonction;
 
 
        MG_SURFACE* surface = face->get_surface();
        TPL_LISTE_ENTITE <double> Parametres;
        int type_surface = surface->get_type_geometrique(Parametres);
 
        double Pnt_Insersion[3];
        Pnt_Insersion[0]= Parametres.get(0);
        Pnt_Insersion[1]= Parametres.get(1);
        Pnt_Insersion[2]= Parametres.get(2);
        double Vect_Normal[3];
        Vect_Normal[0]= Parametres.get(3);
        Vect_Normal[1]= Parametres.get(4);
        Vect_Normal[2]= Parametres.get(5);
 
        
        int nb_triangle = face->get_lien_fem_maillage()->get_nb();
        TPL_MAP_ENTITE<FEM_NOEUD*> liste_noeuds0;
 
        for (int j =0; j<nb_triangle;j++)
                {
                FEM_ELEMENT2* mgtri=(FEM_ELEMENT2*)face->get_lien_fem_maillage()->get(j);
                if (!fem_maillage->contient(mgtri)) continue;
                if (mgtri->get_nb_fem_noeud()==3)
                        {
                        FEM_NOEUD* noeud0=mgtri->get_fem_noeud(0);
                        FEM_NOEUD* noeud1=mgtri->get_fem_noeud(1);
                        FEM_NOEUD* noeud2=mgtri->get_fem_noeud(2);
                        liste_noeuds0.ajouter(noeud0);
                        liste_noeuds0.ajouter(noeud1);
                        liste_noeuds0.ajouter(noeud2);
                        }
                if (mgtri->get_nb_fem_noeud()==6)
                        {
                        FEM_NOEUD* noeud0=mgtri->get_fem_noeud(0);
                        FEM_NOEUD* noeud1=mgtri->get_fem_noeud(1);
                        FEM_NOEUD* noeud2=mgtri->get_fem_noeud(2);
                        FEM_NOEUD* noeud3=mgtri->get_fem_noeud(3);
                        FEM_NOEUD* noeud4=mgtri->get_fem_noeud(4);
                        FEM_NOEUD* noeud5=mgtri->get_fem_noeud(5);
                        liste_noeuds0.ajouter(noeud0);
                        liste_noeuds0.ajouter(noeud1);
                        liste_noeuds0.ajouter(noeud2);
                        liste_noeuds0.ajouter(noeud3);
                        liste_noeuds0.ajouter(noeud4);
                        liste_noeuds0.ajouter(noeud5);
                        }
                }
 
        double xmin = 1e308, ymin = 1e308, zmin = 1e308;
        double xmax = -1e308, ymax = -1e308, zmax = -1e308;
 
        for (int k =0; k < liste_noeuds0.get_nb();k++)
                {
                FEM_NOEUD* noeud=liste_noeuds0.get(k);
                double* xyz=noeud->get_coord();
                if (xyz[0]>xmax) xmax=xyz[0];
                if (xyz[0]<xmin) xmin=xyz[0];
                if (xyz[1]>ymax) ymax=xyz[1];
                if (xyz[1]<ymin) ymin=xyz[1];
                if (xyz[2]>zmax) zmax=xyz[2];
                if (xyz[2]<zmin) zmin=xyz[2];
                }
 
        double Epsilon = OT_VECTEUR_3D(xmax-xmin,ymax-ymin,zmax-zmin).get_longueur()/10.0;
 
        for (int k =0; k < liste_noeuds0.get_nb();k++)
                {
                FEM_NOEUD* noeud=liste_noeuds0.get(k);
                fonction.Ajouter_Point_Array_Ideal(noeud->get_x(coef),noeud->get_y(coef),noeud->get_z(coef),k+1);
                }
 
        fonction.Determiner_Plan(liste_noeuds0.get_nb(),Epsilon);
 
        MG_GEOMETRIE* mggeom = fem_maillage->get_mg_geometrie();
        int num_faceid = 0;
        int nb_face=mggeom->get_nb_mg_face();
        for (int l =0; l < nb_face;l++)
                {
                MG_FACE* F = mggeom->get_mg_face(l);
                if(faceso[l]==1) num_faceid = num_faceid+1;
                if(F == face )
                        {
                        //faceso[l] == 1;
                        fonction.Ajouter_Surface_liste_Ideal(num_faceid);
                        goto Face_Trouvee;
                        }
                }
        Face_Trouvee:
 
        //fonction.Creer_Face_Surface();
 
        //fonction.Creer_Surf_Plane(Pnt_Insersion,Vect_Normal);
 
        int nb_boucle =  face->get_nb_mg_boucle();
 
for (int v =0; v<nb_boucle;v++)
        {
        MG_BOUCLE* Boucle = face->get_mg_boucle(v);
        REC_CONTOUR Rec(fem_maillage, Idealiser, fonction, testo, faceso);
        Rec.Reconstruire(Boucle, coef);
 
        fonction.Ajouter_Contour_Face(v);
        } 
}
 
void REC_FACE_PLANE::Enregistrer(char *path)
{
        OCC_FONCTION_RECONSTRUCTION fonction;
        fonction.Enregistrer(path);
}