fct_taille_fem_solution_generateur_isotrope.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fct_taille_fem_solution_generateur_isotrope.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:52 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "mg_gestionnaire.h"
#include <math.h>
#include "ot_mathematique.h"
#include "mailleur_fem.h"
#include <fem_solution.h>
#include "mailleur3d_structure.h"
#include "fct_taille_fem_solution_generateur_isotrope.h"
 
 
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(MG_GESTIONNAIRE *gesttmp,MG_GEOMETRIE* geotmp,double dg,int fechantillon,int nxtmp,int nytmp,int nztmp,char *nomsolution):FCT_TAILLE_FEM_SOLUTION(gesttmp,nxtmp,nytmp,nztmp),geo(geotmp),ecart_nodal(dg),fechan(fechantillon),facteur_grossissement(1.125),typedata(10),carteini(NULL)
{
strcpy(nom,nomsolution);
typedata=10;
}
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(MG_GESTIONNAIRE *gesttmp,MG_MAILLAGE* maitmp,double dg,int fechantillon,int nxtmp,int nytmp,int nztmp,char *nomsolution):FCT_TAILLE_FEM_SOLUTION(gesttmp,nxtmp,nytmp,nztmp),mai(maitmp),ecart_nodal(dg),fechan(fechantillon),facteur_grossissement(1.125),typedata(11),carteini(NULL)
{
strcpy(nom,nomsolution);
typedata=11;
}
 
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(MG_GESTIONNAIRE *gesttmp,MG_GEOMETRIE* geotmp,FCT_TAILLE* carte,int fechantillon,int nxtmp,int nytmp,int nztmp,char *nomsolution):FCT_TAILLE_FEM_SOLUTION(gesttmp,nxtmp,nytmp,nztmp),geo(geotmp),ecart_nodal(0.),fechan(fechantillon),facteur_grossissement(1.125),typedata(10),carteini(carte)
{
strcpy(nom,nomsolution);
typedata=10;
}
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(MG_GESTIONNAIRE *gesttmp,MG_MAILLAGE* maitmp,FCT_TAILLE* carte,int fechantillon,int nxtmp,int nytmp,int nztmp,char *nomsolution):FCT_TAILLE_FEM_SOLUTION(gesttmp,nxtmp,nytmp,nztmp),mai(maitmp),ecart_nodal(0.),fechan(fechantillon),facteur_grossissement(1.125),typedata(11),carteini(carte)
{
strcpy(nom,nomsolution);
typedata=11;
}
 
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE& mdd):FCT_TAILLE_FEM_SOLUTION(mdd),geo(mdd.geo),mai(mdd.mai),ecart_nodal(mdd.ecart_nodal),fechan(mdd.fechan),facteur_grossissement(mdd.facteur_grossissement),typedata(mdd.typedata),carteini(mdd.carteini)
{
strcpy(nom,mdd.nom);
}
 
 
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE::~FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE()
{
  
}
 
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE::construit(void)
{
if (gest==NULL)
  {
  gest=new MG_GESTIONNAIRE;
  gest_alloue=true;
  }
double xmin=1e308;
double ymin=1e308;
double zmin=1e308;
double xmax=-1e308;
double ymax=-1e308;
double zmax=-1e308;
 
if(typedata==10)
{
BOITE_3D boite;
bool affect=false;
int nbvol=geo->get_nb_mg_volume();
MG_ELEMENT_TOPOLOGIQUE::change_pas_echantillon(fechan);
for (int j=0;j<nbvol;j++)
    {
    MG_VOLUME* vol=geo->get_mg_volume(j);
    BOITE_3D tmp=vol->get_boite_3D();
    if (affect==false) boite=tmp;
    else boite=boite+tmp;
    affect=true;
    }
int nbcoque=geo->get_nb_mg_coque();
for (int j=0;j<nbcoque;j++)
   {
    MG_COQUE* coq=geo->get_mg_coque(j);
    BOITE_3D tmp=coq->get_boite_3D();
    if (affect==false) boite=tmp;
    else boite=boite+tmp;
    affect=true;
    }
MG_ELEMENT_TOPOLOGIQUE::restaure_pas_echantillon();
boite_englobante=boite;
}
if(typedata==11) // when there is no geometry, the bounding box is defined by the current mesh
{
  LISTE_MG_NOEUD::iterator it;
  for (MG_NOEUD* no=mai->get_premier_noeud(it);no!=NULL;no=mai->get_suivant_noeud(it))
  {
    if (no->get_x()<xmin) xmin=no->get_x();  
    if (no->get_y()<ymin) ymin=no->get_y();
    if (no->get_z()<zmin) zmin=no->get_z();
    if (no->get_x()>xmax) xmax=no->get_x();
    if (no->get_y()>ymax) ymax=no->get_y();
    if (no->get_z()>zmax) zmax=no->get_z();
  }
boite_englobante.reinit(xmin,ymin,zmin,xmax,ymax,zmax);
}
boite_englobante.change_grosseur(facteur_grossissement);
MAILLEUR3D_STRUCTURE mail(boite_englobante.get_xmin(),boite_englobante.get_xmax(),boite_englobante.get_ymin(),boite_englobante.get_ymax(),boite_englobante.get_zmin(),boite_englobante.get_zmax(),nx,ny,nz,gest);
mail.maille();
MG_MAILLAGE* mai2=gest->get_mg_maillage(gest->get_nb_mg_maillage()-1);
fem = new FEM_MAILLAGE(NULL, mai2,1);
gest->ajouter_fem_maillage(fem);
MAILLEUR_FEM m;
m.maille(fem);
sol=new FEM_SOLUTION(fem,1,nom,1,(char*)"Carte de taille",MAGIC::ENTITE_SOLUTION::ENTITE_NOEUD,MAGIC::TYPE_SOLUTION::TENSEUR);
gest->ajouter_fem_solution(sol);
sol->change_legende(0,"metrique");
pasx=(boite_englobante.get_xmax()-boite_englobante.get_xmin())/nx;
pasy=(boite_englobante.get_ymax()-boite_englobante.get_ymin())/ny;
pasz=(boite_englobante.get_zmax()-boite_englobante.get_zmin())/nz;
construit_particulier();
isotrope=true;
numchamp=0;
sol->active_solution(numchamp);
LISTE_FEM_ELEMENT3::iterator ite3;
for (FEM_ELEMENT3* ele=fem->get_premier_element3(ite3);ele!=NULL;ele=fem->get_suivant_element3(ite3)) 
  element_structure.push_back(ele);
}