occ_cg_forme_volume_ellipsoide.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       occ_cg_forme_volume_ellipsoide.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:55 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#ifdef ALL_OCC
#include "occ_cg_forme_volume_ellipsoide.h"
#include "mg_cg_modele.h"
#include "occ_fonction.h"
#include "occ_fonction_v2017.h"
#include <TopoDS_Shape.hxx>
#include <TopoDS_Edge.hxx>
#include <TopoDS_Wire.hxx>
#include <gp_Elips.hxx>
#include <gp_Ax2.hxx>
#include <gp_Dir.hxx>
#include <gp_Pnt.hxx>
#include <BRepBuilderAPI_MakeEdge.hxx>
#include <BRepBuilderAPI_MakeWire.hxx>
#include <BRepBuilderAPI_MakeFace.hxx>
#include <BRepPrimAPI_MakeRevol.hxx>
#include <math.h>
OCC_CG_FORME_VOLUME_ELLIPSOIDE::OCC_CG_FORME_VOLUME_ELLIPSOIDE(MG_CG_MODELE* mgcg_modele,
                                                               double centre_x, double centre_y, double centre_z,
                                                               double direction_x, double direction_y, double direction_z,
                                                               double rayon_majeur, double rayon_mineur): MG_CG_FORME_VOLUME_ELLIPSOIDE(mgcg_modele,centre_x, centre_y, centre_z, direction_x, direction_y, direction_z, rayon_majeur, rayon_mineur)
{
 
}
 
OCC_CG_FORME_VOLUME_ELLIPSOIDE::OCC_CG_FORME_VOLUME_ELLIPSOIDE(MG_CG_MODELE* mgcg_modele,
                                                               int etat,
                                                               long unsigned int num,
                                                               MG_VOLUME* mg_volume,
                                                               double centre_x, double centre_y, double centre_z,
                                                               double direction_x, double direction_y, double direction_z,
                                                               double rayon_majeur, double rayon_mineur): MG_CG_FORME_VOLUME_ELLIPSOIDE(mgcg_modele,etat, num, mg_volume, centre_x, centre_y, centre_z, direction_x, direction_y, direction_z, rayon_majeur, rayon_mineur)
{
 
}
 
OCC_CG_FORME_VOLUME_ELLIPSOIDE::OCC_CG_FORME_VOLUME_ELLIPSOIDE(OCC_CG_FORME_VOLUME_ELLIPSOIDE& mdd): MG_CG_FORME_VOLUME_ELLIPSOIDE(mdd)
{
 
}
 
OCC_CG_FORME_VOLUME_ELLIPSOIDE::~OCC_CG_FORME_VOLUME_ELLIPSOIDE()
{
 
}
 
int OCC_CG_FORME_VOLUME_ELLIPSOIDE::construire(void)
{
  if(m_mgcg_modele==NULL)
  {
    std::cout << "*** ERREUR : OCC_CG_FORME_VOLUME_ELLIPSOIDE::construire_forme -> MG_CG_MODELE NULL ***" << std::endl;
    return FAIL;
  }
  return construire(m_mgcg_modele->get_fusionner_entite_similaire(),
                    m_mgcg_modele->get_precision(),
                    m_mgcg_modele->get_importer_triangulation(),
                    m_mgcg_modele->get_epsilon_triangulation());
}
 
int OCC_CG_FORME_VOLUME_ELLIPSOIDE::construire(bool fusionner_entite_similaire, double precision, bool importer_triangulation, double epsilon_triangulation)
{
  if(m_mgcg_modele==NULL)
  {
    std::cout << "*** ERREUR : OCC_CG_FORME_VOLUME_ELLIPSOIDE::construire_forme -> MG_CG_MODELE NULL ***" << std::endl;
    return FAIL;
  }
  gp_Pnt centre_ellipsoide(m_centre[0],m_centre[1],m_centre[2]);
  gp_Dir direction_ellipsoide (m_direction[0],m_direction[1],m_direction[2]);  
  gp_Ax2 axe_ellipsoide;
  axe_ellipsoide.SetLocation(centre_ellipsoide);
  axe_ellipsoide.SetDirection(direction_ellipsoide);
  axe_ellipsoide.Rotate(gp_Ax1(centre_ellipsoide,axe_ellipsoide.YDirection()),M_PI/2.);
  gp_Elips ellipse;
  if(m_rayon_majeur>=m_rayon_mineur)
  {
    ellipse = gp_Elips(axe_ellipsoide,m_rayon_majeur,m_rayon_mineur);
  }
  else
  {
    ellipse = gp_Elips(axe_ellipsoide,m_rayon_mineur,m_rayon_majeur);
  } 
  TopoDS_Edge edge ;
  edge = BRepBuilderAPI_MakeEdge(ellipse,0.0,M_PI).Edge();
  TopoDS_Wire wire;
  wire = BRepBuilderAPI_MakeWire(edge).Wire();
  TopoDS_Face face;
  face = BRepBuilderAPI_MakeFace(wire,true);  
  gp_Dir direction_axe_revolution(m_direction[0],m_direction[1],m_direction[2]);
  gp_Ax1 axe_revolution;
  axe_revolution.SetLocation(centre_ellipsoide);
  axe_revolution.SetDirection(direction_axe_revolution);
  TopoDS_Shape  topods_shape = BRepPrimAPI_MakeRevol(face,axe_revolution,2.*M_PI).Shape();  
  BRepTools::Write(topods_shape,(char*)"test.brep");
  OCC_FONCTION* occ_fonction=&m_mgcg_modele->get_mg_geometrie()->get_occ_fonction();
  if(occ_fonction->get_version()=="OCCV2017")
  {
    OCC_FONCTION_V2017* occ_fonction_v2017 = (OCC_FONCTION_V2017*)occ_fonction;
    m_mg_volume =  occ_fonction_v2017->importer_TopoDS_Solid(TopoDS::Solid(topods_shape),
                                                             NULL,
                                                             NULL,
                                                             NULL,
                                                             m_mgcg_modele->get_mg_geometrie(),
                                                             fusionner_entite_similaire,
                                                             precision);
  }
  change_etat(MG_CG_FORME::ETAT_FORME::CONSTRUITE);
  m_mgcg_modele->mise_a_jout_lien_mg_element_mgcg_forme(this);
  return OK;
}
 
#endif