mg_calcule_proprietes.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_calcule_proprietes.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "mg_calcule_proprietes.h"
#include "mg_file.h"
#include "ot_parametres.h"
#include <fct_generateur_calibrage.h>
#include <math.h>
 
MG_CALCULE_PROPRIETES::MG_CALCULE_PROPRIETES():affichageactif(0)
{
}
 
MG_CALCULE_PROPRIETES::MG_CALCULE_PROPRIETES(MG_CALCULE_PROPRIETES& mdd):affichageactif(mdd.affichageactif)
{
}
 
MG_CALCULE_PROPRIETES::~MG_CALCULE_PROPRIETES()
{
}
 
void MG_CALCULE_PROPRIETES::active_affichage(void (*fonc)(char*))
{
affiche=fonc;
affichageactif=1;
}
 
char MG_CALCULE_PROPRIETES::get_seprateur_decimale(void)
{
OT_PARAMETRES param;
param.ajouter("Separateur_decimale",",",OT_PARAMETRES::STRING,"Seprateur décimale pour l'exportation des fichiers excel");
char fichier[500];
sprintf(fichier,"%s/.magic",getenv("HOME"));
param.lire(fichier);
param.enregistrer(fichier);
char c=(param.get_nom("Separateur_decimale").c_str())[0];
return c;
}
 
int MG_CALCULE_PROPRIETES::calcul_masse_humidite(char *fichiersol,char *Resultat,MG_GESTIONNAIRE &gest,double &resu,double tf,int pas,double m0,int Pasvariable,int division_pas,MG_GEOMETRIE *geo)
{
FEM_MAILLAGE* fem;
FEM_SOLUTION* sol;
char c=get_seprateur_decimale();
recupere_maillage(gest,&fem);
char nom1[500];
char nom2[500];
sprintf(nom1,"%s.csv",Resultat);
sprintf(nom2,"%s.txt",Resultat);
FILE* file_resu1=fopen(nom1,"wt");
FILE* file_resu2=fopen(nom2,"wt");
fprintf (file_resu1, "%s ; %s ; %s ; %s ; %s ; %s ; %s\n","Pas","Nombre de jours","Racine carrée du jour","Prise de masse en Kg","Prise de masse en g","masse echantillon","Prise de masse en pourcentage");
fprintf (file_resu2, "%s  %s  %s  %s  %s  %s  %s\n","Pas","Nombre_de_jours","Racine_carrée_ du_jour","Prise_de_masse_en_Kg","Prise_de_masse_en_g","masse_echantillon","Prise_de_masse_en_pourcentage");
int l=1;
for(int k=0;k<(pas+1);k++)
{
recupere_concentration(fichiersol,gest,fem,&sol,k);
sol->active_solution(0);
 
double masse=0.;
double jac[9];
double unite=fem->get_mg_geometrie()->get_valeur_unite(); 
LISTE_FEM_ELEMENT3::iterator it;
for (FEM_ELEMENT3* ele=fem->get_premier_element3(it);ele!=NULL;ele=fem->get_suivant_element3(it))
        {        
       double w;
       double uvw[3];
       double degmax=ele->get_degremax_fonction_interpolation();
       int nbfonc=ele->get_nb_fonction_interpolation();
       int nbgauss=ele->get_nb_pt_gauss(degmax); 
       double M;
       M=0;
       for(int i=0;i<nbgauss;i++)
        {
        ele->get_pt_gauss(degmax,i,w,uvw);
        double uvw1[3]={uvw[0],uvw[1],uvw[2]};
        //Determinant du Jacobien//
        double detJ=ele->get_jacobien(jac,uvw1,1.);
        //Les fonction d'interpolation//
        
        // Concentrations aux noeuds//
        if (nbfonc==4)
        {
        double N1=ele->get_fonction_interpolation(1,uvw1);
        double N2=ele->get_fonction_interpolation(2,uvw1);  
        double N3=ele->get_fonction_interpolation(3,uvw1);
        double N4=ele->get_fonction_interpolation(4,uvw1);
        double C1= (ele->get_fem_noeud(0)->get_solution())*(pow(unite,3));
        double C2= (ele->get_fem_noeud(1)->get_solution())*(pow(unite,3));
        double C3= (ele->get_fem_noeud(2)->get_solution())*(pow(unite,3));
        double C4= (ele->get_fem_noeud(3)->get_solution())*(pow(unite,3));
        double F=(N1*C1+N2*C2+N3*C3+N4*C4)*fabs(detJ);
        M=M+(w*F); 
        }
        if (nbfonc==8)
        {
        double N1=ele->get_fonction_interpolation(1,uvw1);
        double N2=ele->get_fonction_interpolation(2,uvw1);  
        double N3=ele->get_fonction_interpolation(3,uvw1);
        double N4=ele->get_fonction_interpolation(4,uvw1);
        double N5=ele->get_fonction_interpolation(5,uvw1);
        double N6=ele->get_fonction_interpolation(6,uvw1);
        double N7=ele->get_fonction_interpolation(7,uvw1);
        double N8=ele->get_fonction_interpolation(8,uvw1);
        double C1= (ele->get_fem_noeud(0)->get_solution())*(pow(unite,3));
        double C2= (ele->get_fem_noeud(1)->get_solution())*(pow(unite,3));
        double C3= (ele->get_fem_noeud(2)->get_solution())*(pow(unite,3));
        double C4= (ele->get_fem_noeud(3)->get_solution())*(pow(unite,3));
        double C5=(ele->get_fem_noeud(4)->get_solution())*(pow(unite,3));
        double C6= (ele->get_fem_noeud(5)->get_solution())*(pow(unite,3));
        double C7= (ele->get_fem_noeud(6)->get_solution())*(pow(unite,3));
        double C8= (ele->get_fem_noeud(7)->get_solution())*(pow(unite,3));
        double F=(N1*C1+N2*C2+N3*C3+N4*C4+N5*C5+N6*C6+N7*C7+N8*C8)*fabs(detJ);
        M=M+(w*F); 
        }
        
        
         } 
     masse=masse+M;
    
        }
resu=8*masse;
double jours;
if (Pasvariable==1)
{
 MG_GEOM_FONCTION *fonc;  
      int val_Cond1;
      int val_Cond2 = geo->get_nb_mg_arete();
      int val_Cond3 = geo->get_nb_mg_sommet();
      int val_Cond4 = geo->get_nb_mg_face();
      int val_Cond5 = geo->get_nb_mg_geom_fonction();
      if (val_Cond5==1) val_Cond1=3+val_Cond2+val_Cond3+val_Cond4;   //fonctionne avec la géométrie virtuelle
      if (val_Cond5==2) val_Cond1=4+val_Cond2+val_Cond3+val_Cond4;   //fonctionne avec la géométrie virtuelle
      
          fonc = geo->get_mg_geom_fonctionid((int)val_Cond1);
          double col[2];
          double col1[2];
          double diff;
          int s;
          fonc->get_point(0,col);
          jours=col[1]/(86400);
          if(k%division_pas!=0 && k!=0)
          {  
          s=(division_pas*l)-k;  
          fonc->get_point(l,col);
          fonc->get_point(l-1,col1);
          diff=(col[1]-col1[1])/(86400*division_pas);
          jours=(col[1]/86400)-s*diff;
          }
          if(k%division_pas==0 && k!=0)
          {
          fonc->get_point(l,col);
          jours=col[1]/(86400);
          l=l+1;
          }       
}
else
{
jours=((tf/pas)*k)/86400;
}
double prise_pourcentage=((resu*(1/unite))/m0)*100;
char chaine[5000];
if (k>=division_pas)
{
sprintf (chaine, "%d ; %lf ; %lf ; %lf ; %lf ; %lf ; %lf",k,jours,sqrt(jours),resu,resu*(1/unite),m0+resu*(1/unite),prise_pourcentage);  
for (int i=0;i<strlen(chaine);i++)
  if (chaine[i]=='.') chaine[i]=c;
fprintf(file_resu1,"%s\n",chaine);  
sprintf (chaine, "%d  %lf  %lf  %lf  %lf  %lf  %lf",k,jours,sqrt(jours),resu,resu*(1/unite),m0+resu*(1/unite),prise_pourcentage);   
fprintf(file_resu2,"%s\n",chaine); 
}
}
fclose(file_resu2);
fclose(file_resu1);
return 0;
}
int MG_CALCULE_PROPRIETES::recupere_maillage(MG_GESTIONNAIRE &gest,FEM_MAILLAGE **femmesh)
{
FEM_MAILLAGE *fem=gest.get_fem_maillage(0);
*femmesh=fem;
return 0;
}
 
int MG_CALCULE_PROPRIETES::recupere_concentration(char *fichiersol,MG_GESTIONNAIRE &gest,FEM_MAILLAGE *fem,FEM_SOLUTION **solmesh,int k)
{
FEM_SOLUTION *sol=gest.get_fem_solution(k);
*solmesh=sol;
return 0;
}