mstruct_analyse_erosion.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mstruct_analyse_erosion.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "mstruct_analyse_erosion.h"
#include "mstruct_analyse_champ.h"
#include "mstruct_analyse_orientation.h"
#include "mstruct_analyse_cao.h"
#include "mstruct_analyse_mg_maillage.h"
#include "mstruct_analyse_fem_maillage.h"
#include "mstruct_analyse_modules_elasticite.h"
#include "mstruct_analyse_conductivite_thermique.h"
#include "mstruct_analyse_energie_hill.h"
#include "mstruct_definition.h"
#include "ot_boite_3d.h"
MSTRUCT_ANALYSE_EROSION::MSTRUCT_ANALYSE_EROSION(void): MSTRUCT_ANALYSE()
{
  m_epaisseur_couche=0;
}
 
MSTRUCT_ANALYSE_EROSION::MSTRUCT_ANALYSE_EROSION(std::string identifiant,
                                                 double epaisseur_couche,
                                                 std::string nom_groupe_forme,
                                                 BOITE_3D* boite_3d):MSTRUCT_ANALYSE(identifiant,nom_groupe_forme,boite_3d)
{
  m_epaisseur_couche=epaisseur_couche;
}
 
 
MSTRUCT_ANALYSE_EROSION::MSTRUCT_ANALYSE_EROSION(std::string identifiant,
                                                 MSTRUCT_ANALYSE* analyse_initiale,
                                                 long nb_couche,
                                                 double epaisseur_couche): MSTRUCT_ANALYSE(identifiant,
                                                                                           analyse_initiale->get_nom_groupe_forme(),
                                                                                           analyse_initiale->get_boite_analyse()),
                                                                                           m_epaisseur_couche(epaisseur_couche)
{
  std::string iden_ini=identifiant+"_0";
  analyse_initiale->change_identifiant(iden_ini);
  ajouter_analyse(analyse_initiale);
  int type_analyse = analyse_initiale->get_type();
  for(long i=1;i<nb_couche;i++)
  {
    MSTRUCT_ANALYSE* nouvelle_analyse;
    if(type_analyse==MSTRUCT::TYPE_ANALYSE::CHAMP)
    {
      MSTRUCT_ANALYSE_CHAMP* analyse_ini_champ = (MSTRUCT_ANALYSE_CHAMP*)analyse_initiale;
      MSTRUCT_ANALYSE_CHAMP* analyse_champ = new MSTRUCT_ANALYSE_CHAMP(*analyse_ini_champ);
      nouvelle_analyse=analyse_champ;      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::ORIENTATION)
    {
      MSTRUCT_ANALYSE_ORIENTATION* analyse_ini_ori = (MSTRUCT_ANALYSE_ORIENTATION*)analyse_initiale;
      MSTRUCT_ANALYSE_ORIENTATION* analyse_ori = new MSTRUCT_ANALYSE_ORIENTATION(*analyse_ini_ori);
      nouvelle_analyse=analyse_ori;  
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::CAO)
    {
      MSTRUCT_ANALYSE_CAO* analyse_ini_cao = (MSTRUCT_ANALYSE_CAO*)analyse_initiale;
      MSTRUCT_ANALYSE_CAO* analyse_cao = new MSTRUCT_ANALYSE_CAO(*analyse_ini_cao);
      nouvelle_analyse=analyse_cao;  
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::MAILLAGE_MG)
    {
      MSTRUCT_ANALYSE_MG_MAILLAGE* analyse_ini_mg_maill = (MSTRUCT_ANALYSE_MG_MAILLAGE*)analyse_initiale;
      MSTRUCT_ANALYSE_MG_MAILLAGE* analyse_mg_maill = new MSTRUCT_ANALYSE_MG_MAILLAGE(*analyse_ini_mg_maill);
      nouvelle_analyse=analyse_mg_maill;  
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::MAILLAGE_FEM)
    {
      MSTRUCT_ANALYSE_FEM_MAILLAGE* analyse_ini_fem_maill = (MSTRUCT_ANALYSE_FEM_MAILLAGE*)analyse_initiale;
      MSTRUCT_ANALYSE_FEM_MAILLAGE* analyse_fem_maill = new MSTRUCT_ANALYSE_FEM_MAILLAGE(*analyse_ini_fem_maill);
      nouvelle_analyse=analyse_fem_maill;  
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::CONDUCTIVITE_THERMIQUE)
    {
      MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE* analyse_ini_cond_ther = (MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE*)analyse_initiale;
      MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE* analyse_cond_ther = new MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE(*analyse_ini_cond_ther);
      nouvelle_analyse=analyse_cond_ther;  
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::ENERGIE_HILL)
    {
      MSTRUCT_ANALYSE_ENERGIE_HILL* analyse_ini_cond_ther = (MSTRUCT_ANALYSE_ENERGIE_HILL*)analyse_initiale;
      MSTRUCT_ANALYSE_ENERGIE_HILL* analyse_cond_ther = new MSTRUCT_ANALYSE_ENERGIE_HILL(*analyse_ini_cond_ther);
      nouvelle_analyse=analyse_cond_ther;  
    }
    char iden[500];
    sprintf(iden,"%s_%li",identifiant.c_str(),i);
    nouvelle_analyse->change_identifiant(iden);
    if(nouvelle_analyse->get_boite_analyse()==NULL)
    {
      nouvelle_analyse->change_boite_analyse(BOITE_3D(0.0,0.0,0.0,1.0,1.0,1.0));
    }
    double xmin,ymin,zmin,xmax,ymax,zmax;
    xmin=nouvelle_analyse->get_boite_analyse()->get_xmin();
    ymin=nouvelle_analyse->get_boite_analyse()->get_ymin();
    zmin=nouvelle_analyse->get_boite_analyse()->get_zmin();
    xmax=nouvelle_analyse->get_boite_analyse()->get_xmax();
    ymax=nouvelle_analyse->get_boite_analyse()->get_ymax();
    zmax=nouvelle_analyse->get_boite_analyse()->get_zmax();
    xmin=xmin+i*m_epaisseur_couche;
    ymin=ymin+i*m_epaisseur_couche;
    zmin=zmin+i*m_epaisseur_couche;
    xmax=xmax-i*m_epaisseur_couche;
    ymax=ymax-i*m_epaisseur_couche;
    zmax=zmax-i*m_epaisseur_couche;
    nouvelle_analyse->get_boite_analyse()->reinit(xmin,ymin,zmin,xmax,ymax,zmax);
    ajouter_analyse(nouvelle_analyse);
  }
}
 
 
MSTRUCT_ANALYSE_EROSION::~MSTRUCT_ANALYSE_EROSION(void)
{
  std::vector<MSTRUCT_ANALYSE*>::iterator it_analyse;
  for(MSTRUCT_ANALYSE* analyse=get_premiere_analyse(it_analyse);analyse!=NULL;analyse=get_suivante_analyse(it_analyse)) delete analyse;
}
 
double MSTRUCT_ANALYSE_EROSION::get_epaisseur_couche(void)
{
  return m_epaisseur_couche;
}
 
long int MSTRUCT_ANALYSE_EROSION::get_nb_analyse(void)
{
  return m_vector_analyse.size();
}
 
int MSTRUCT_ANALYSE_EROSION::ajouter_analyse(MSTRUCT_ANALYSE* analyse)
{
  char iden[500];
  sprintf(iden,"%s_%li",m_identifiant.c_str(),get_nb_analyse());
  analyse->change_identifiant(iden);
  m_vector_analyse.push_back(analyse);
  return 0;
}
 
MSTRUCT_ANALYSE* MSTRUCT_ANALYSE_EROSION::get_premiere_analyse(std::vector<MSTRUCT_ANALYSE*>::iterator& it)
{
  it=m_vector_analyse.begin();
  if(it!=m_vector_analyse.end()) return *it;
  else return NULL;
}
 
MSTRUCT_ANALYSE* MSTRUCT_ANALYSE_EROSION::get_suivante_analyse(std::vector<MSTRUCT_ANALYSE*>::iterator& it)
{
  it++;
  if(it!=m_vector_analyse.end()) return *it;
  else return NULL;
}
 
MSTRUCT_ANALYSE* MSTRUCT_ANALYSE_EROSION::get_analyse(long num)
{
  return m_vector_analyse.at(num);
}
 
long int MSTRUCT_ANALYSE_EROSION::get_type(void)
{
  return MSTRUCT::TYPE_ANALYSE::EROSION;
}
 
void MSTRUCT_ANALYSE_EROSION::executer(void)
{
  std::vector<MSTRUCT_ANALYSE*>::iterator it_analyse;
  for(MSTRUCT_ANALYSE* analyse=get_premiere_analyse(it_analyse);analyse!=NULL;analyse=get_suivante_analyse(it_analyse)) analyse->executer();
}
 
void MSTRUCT_ANALYSE_EROSION::exporter(std::ofstream& ofstrm, long i, bool avec_entete, bool avec_graph, char* prefix_graph)
{
  std::vector<MSTRUCT_ANALYSE*>::iterator it_analyse;
  for(MSTRUCT_ANALYSE* analyse=get_premiere_analyse(it_analyse);analyse!=NULL;analyse=get_suivante_analyse(it_analyse))
  {
    ofstrm << analyse->get_identifiant() << std::endl;
    analyse->exporter(ofstrm,i,avec_entete,avec_graph,prefix_graph);
  }
}
 
void MSTRUCT_ANALYSE_EROSION::enregistrer(std::ofstream& ofstrm)
{
  long type_analyse=get_type();
  ofstrm.write((char*)&type_analyse,sizeof(long));
  MSTRUCT_ANALYSE::enregistrer(ofstrm);
  ofstrm.write((char*)&m_epaisseur_couche,sizeof(double));
  long nb_analyse=get_nb_analyse();
  ofstrm.write((char*)&nb_analyse,sizeof(long));
  std::vector<MSTRUCT_ANALYSE*>::iterator it_analyse;
  for(MSTRUCT_ANALYSE* analyse=get_premiere_analyse(it_analyse);analyse!=NULL;analyse=get_suivante_analyse(it_analyse)) analyse->enregistrer(ofstrm);  
}
 
void MSTRUCT_ANALYSE_EROSION::ouvrir(std::ifstream& ifstrm)
{
  MSTRUCT_ANALYSE::ouvrir(ifstrm);
  ifstrm.read((char*)&m_epaisseur_couche,sizeof(double));
  long nb_analyse;
  ifstrm.read((char*)&nb_analyse,sizeof(long));
  for(long i=0;i<nb_analyse;i++)
  {
    long type_analyse;
    ifstrm.read((char*)&type_analyse,sizeof(long));
    MSTRUCT_ANALYSE *analyse;
    if(type_analyse==MSTRUCT::TYPE_ANALYSE::CHAMP)
    {
      analyse = new MSTRUCT_ANALYSE_CHAMP;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::ORIENTATION)
    {
      analyse = new MSTRUCT_ANALYSE_ORIENTATION;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::CAO)
    {
      analyse = new MSTRUCT_ANALYSE_CAO;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::MAILLAGE_MG)
    {
      analyse = new MSTRUCT_ANALYSE_MG_MAILLAGE;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::MAILLAGE_FEM)
    {
      analyse = new MSTRUCT_ANALYSE_FEM_MAILLAGE;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::MODULES_ELASTICITE)
    {
      analyse = new MSTRUCT_ANALYSE_MODULES_ELASTICITE;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::CONDUCTIVITE_THERMIQUE)
    {
      analyse = new MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE;
      analyse->ouvrir(ifstrm);      
    }
    else if(type_analyse==MSTRUCT::TYPE_ANALYSE::ENERGIE_HILL)
    {
      analyse = new MSTRUCT_ANALYSE_ENERGIE_HILL;
      analyse->ouvrir(ifstrm);      
    }
    ajouter_analyse(analyse);
  }
}
 
void MSTRUCT_ANALYSE_EROSION::affiche_contenu(fonction_affiche* fonc)
{
  MSTRUCT_ANALYSE::affiche_contenu(fonc);
  char ligne[5000];
  sprintf(ligne,"MSTRUCT_ANALYSE_EROSION");fonc(ligne);
  sprintf(ligne,"-> Nb analyse : %li",get_nb_analyse());fonc(ligne);
  sprintf(ligne,"-> Epaisseur couche  : %lf",m_epaisseur_couche);fonc(ligne);
  std::vector<MSTRUCT_ANALYSE*>::iterator it_analyse;
  for(MSTRUCT_ANALYSE* analyse=get_premiere_analyse(it_analyse);analyse!=NULL;analyse=get_suivante_analyse(it_analyse))
  {
    analyse->affiche_contenu(fonc);
  }
}