mg_file_recuperation.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_file_recuperation.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
 
#pragma hdrstop
 
#include <string.h>
#include "mg_file_recuperation.h"
#include "mg_sommet_noeud.h"
#include "pars_argument.h"
#include "parse.h"
#ifdef BREP_STEP
#include "step_surface.h"
#include "step_courbe.h"
#include "step_point.h"
#endif
#ifdef BREP_SAT
#include "acis_surface.h"
#include "acis_courbe.h"
#include "acis_point.h"
#endif
#ifdef BREP_SLD
#include "sld_point.h"
#include "sld_courbe.h"
#include "sld_surface.h"
#include "sld_fonction.h"
#endif
#include "lc_point.h"
#include "m3d_noeud.h"
#include "m3d_triangle.h"
#include "m3d_tetra.h"
#include "fem_segment2.h"
#include "fem_segment3.h"
#include "fem_triangle3.h"
#include "fem_triangle6.h"
#include "fem_quadrangle4.h"
#include "fem_quadrangle8.h"
#include "fem_tetra4.h"
#include "xfem_segment2.h"
#include "xfem_triangle3.h"
#include "xfem_tetra4.h"
#include "fem_tetra10.h"
#include "fem_hexa8.h"
#include "fem_hexa20.h"
#ifdef BREP_OCC
#include "occ_surface.h"
#include "occ_point.h"
#include "occ_courbe.h"
#include "occ_fonction.h"
#include "TopoDS_Shape.hxx"
#include "TopoDS.hxx"
#endif
#ifdef WINDOWS_VERSION
#include "CAD4FE_MCEdge.h"
#include "CAD4FE_MCFace.h"
#include "CAD4FE_MCVertex.h"
#include "CAD4FE_PolySurface.h"
#include "CAD4FE_PolyCurve.h"
#include "CAD4FE_MCNode.h"
#include "CAD4FE_MCSegment.h"
#include "CAD4FE_MCTriangle.h"
#include "CAD4FE_m3d_MCTriangle.h"
#endif
 
 
#pragma package(smart_init)
 
 
MG_FILE_RECUPERATION::MG_FILE_RECUPERATION(char* chemin):MG_FILE()
{
 code_de_lecture=lire(chemin);
  
}
 
MG_FILE_RECUPERATION::~MG_FILE_RECUPERATION()
{
}
 
 
void MG_FILE_RECUPERATION::lire_NOEUD(long id,class PARSE& parse,class PARS_ARGUMENT (&param)[100],std::string& data,MG_MAILLAGE *mgmai,MG_GEOMETRIE *mggeo)
{
parse.decode(data.c_str(),"@,@,@,@,@",param+2);
long idtopo=cid(param[2].argument[0]);
double x=atof(param[3].argument[0].c_str());
double y=atof(param[4].argument[0].c_str());
double z=atof(param[5].argument[0].c_str());
int ori=atoi(param[6].argument[0].c_str());
MG_ELEMENT_TOPOLOGIQUE* topo;
if (idtopo>-1)
      {
      topo=mggeo->get_mg_sommetid(idtopo);
      if (topo==NULL) topo=mggeo->get_mg_areteid(idtopo);
      if (topo==NULL) topo=mggeo->get_mg_faceid(idtopo);
      if (topo==NULL) topo=mggeo->get_mg_volumeid(idtopo);
      }
else topo=NULL;
MG_NOEUD* noeud=new M3D_NOEUD(id,topo,x,y,z,ori);
mgmai->ajouter_mg_noeud(noeud);
}
 
            
 
void MG_FILE_RECUPERATION::lire_TRIANGLE(long id,class PARSE& parse,class PARS_ARGUMENT (&param)[100],std::string& data,MG_MAILLAGE *mgmai,MG_GEOMETRIE *mggeo)
{
parse.decode(data.c_str(),"@,@,@,@,@",param+2);
long idtopo=cid(param[2].argument[0]);
long idn1=cid(param[3].argument[0]);
long idn2=cid(param[4].argument[0]);
long idn3=cid(param[5].argument[0]);
int ori=atoi(param[6].argument[0].c_str());
MG_NOEUD* noeud1=mgmai->get_mg_noeudid(idn1);
MG_NOEUD* noeud2=mgmai->get_mg_noeudid(idn2);
MG_NOEUD* noeud3=mgmai->get_mg_noeudid(idn3);
MG_ELEMENT_TOPOLOGIQUE* topo=NULL;
if (idtopo>-1)
    {
    topo=mggeo->get_mg_faceid(idtopo);
    if (topo==NULL) topo=mggeo->get_mg_volumeid(idtopo);
    }
MG_SEGMENT* mgsegment1=mgmai->get_mg_segment(noeud1->get_id(),noeud2->get_id());
MG_SEGMENT* mgsegment2=mgmai->get_mg_segment(noeud2->get_id(),noeud3->get_id());
MG_SEGMENT* mgsegment3=mgmai->get_mg_segment(noeud3->get_id(),noeud1->get_id());
if (mgsegment1==NULL) mgsegment1=mgmai->ajouter_mg_segment(topo,noeud1,noeud2,ori);
if (mgsegment2==NULL) mgsegment2=mgmai->ajouter_mg_segment(topo,noeud2,noeud3,ori);
if (mgsegment3==NULL) mgsegment3=mgmai->ajouter_mg_segment(topo,noeud3,noeud1,ori);
M3D_TRIANGLE* mgtriangle=new M3D_TRIANGLE(id,topo,noeud1,noeud2,noeud3,mgsegment1,mgsegment2,mgsegment3,ori);
mgtriangle->change_frontiere();
mgtriangle->get_segment1()->change_dimension_topo_null(2);
mgtriangle->get_segment2()->change_dimension_topo_null(2);
mgtriangle->get_segment3()->change_dimension_topo_null(2);
mgmai->ajouter_mg_triangle(mgtriangle);
}
 
 
void MG_FILE_RECUPERATION::lire_TETRA(long id,class PARSE& parse,class PARS_ARGUMENT (&param)[100],std::string& data,MG_MAILLAGE *mgmai,MG_GEOMETRIE *mggeo)
{
parse.decode(data.c_str(),"@,@,@,@,@,@",param+2);
long idtopo=cid(param[2].argument[0]);
long idn1=cid(param[3].argument[0]);
long idn2=cid(param[4].argument[0]);
long idn3=cid(param[5].argument[0]);
long idn4=cid(param[6].argument[0]);
int ori=atoi(param[7].argument[0].c_str());
MG_NOEUD* noeud1=mgmai->get_mg_noeudid(idn1);
MG_NOEUD* noeud2=mgmai->get_mg_noeudid(idn2);
MG_NOEUD* noeud3=mgmai->get_mg_noeudid(idn3);
MG_NOEUD* noeud4=mgmai->get_mg_noeudid(idn4);
MG_ELEMENT_TOPOLOGIQUE* topo;    
if (idtopo>-1) topo=mggeo->get_mg_volumeid(idtopo);
else topo=NULL;
MG_TRIANGLE* triangle1=mgmai->get_mg_triangle(noeud1->get_id(),noeud3->get_id(),noeud2->get_id());
MG_TRIANGLE* triangle2=mgmai->get_mg_triangle(noeud1->get_id(),noeud2->get_id(),noeud4->get_id());
MG_TRIANGLE* triangle3=mgmai->get_mg_triangle(noeud2->get_id(),noeud3->get_id(),noeud4->get_id());
MG_TRIANGLE* triangle4=mgmai->get_mg_triangle(noeud1->get_id(),noeud4->get_id(),noeud3->get_id());
if (triangle1==NULL) triangle1=insere_triangle(mgmai,topo,noeud1,noeud3,noeud2,MAGIC::ORIGINE::MAILLEUR_AUTO);
if (triangle2==NULL) triangle2=insere_triangle(mgmai,topo,noeud1,noeud2,noeud4,MAGIC::ORIGINE::MAILLEUR_AUTO);
if (triangle3==NULL) triangle3=insere_triangle(mgmai,topo,noeud2,noeud3,noeud4,MAGIC::ORIGINE::MAILLEUR_AUTO);
if (triangle4==NULL) triangle4=insere_triangle(mgmai,topo,noeud1,noeud4,noeud3,MAGIC::ORIGINE::MAILLEUR_AUTO);
M3D_TETRA* tetra=new M3D_TETRA(topo,noeud1,noeud2,noeud3,noeud4,triangle1,triangle2,triangle3,triangle4,ori);
mgmai->ajouter_mg_tetra(tetra);
M3D_TRIANGLE* tri=(M3D_TRIANGLE*)triangle1;
tri->ajouter_voisin(tetra);
tri=(M3D_TRIANGLE*)triangle2;
tri->ajouter_voisin(tetra);
tri=(M3D_TRIANGLE*)triangle3;
tri->ajouter_voisin(tetra);
tri=(M3D_TRIANGLE*)triangle4;
tri->ajouter_voisin(tetra);
}
 
 
class MG_TRIANGLE* MG_FILE_RECUPERATION::insere_triangle(MG_MAILLAGE* mgmai, MG_ELEMENT_TOPOLOGIQUE* mgvol,MG_NOEUD* noeud1,MG_NOEUD* noeud2,MG_NOEUD* noeud3,int origine)
{
    MG_SEGMENT* segment1=mgmai->get_mg_segment(noeud1->get_id(),noeud2->get_id());
    MG_SEGMENT* segment2=mgmai->get_mg_segment(noeud1->get_id(),noeud3->get_id());
    MG_SEGMENT* segment3=mgmai->get_mg_segment(noeud2->get_id(),noeud3->get_id());
    if (segment1==NULL) segment1=mgmai->ajouter_mg_segment(mgvol,noeud1,noeud2,origine);
    if (segment2==NULL) segment2=mgmai->ajouter_mg_segment(mgvol,noeud1,noeud3,origine);
    if (segment3==NULL) segment3=mgmai->ajouter_mg_segment(mgvol,noeud2,noeud3,origine);
    M3D_TRIANGLE* tri=new M3D_TRIANGLE(mgvol,noeud1,noeud2,noeud3,segment1,segment2,segment3,origine);
    mgmai->ajouter_mg_triangle(tri);
    return tri;
}