mstruct_analyse_conductivite_thermique.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_conductivite_thermique.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:56 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "mstruct_analyse_conductivite_thermique.h"
#include "mstruct_definition.h"
#include "ot_boite_3d.h"
#include "fem_solution.h"
#include "mstruct_ves.h"
#include "mg_gestionnaire.h"
#include "mg_element_topologique.h"
#include <math.h>
#include <vector>
MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE():MSTRUCT_ANALYSE()
{
  m_id_fem_solution_temp=-1;
  m_id_fem_solution_flux=-1;
  m_direction=-1;
  m_lambda_x=NULL;
  m_lambda_y=NULL;
  m_lambda_z=NULL;
}
 
MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE(MSTRUCT_VES* ves,
                                                                               std::string identifiant,
                                                                               long id_fem_solution_temp,
                                                                               long id_fem_solution_flux,
                                                                               int direction,
                                                                               BOITE_3D* boite_3d):MSTRUCT_ANALYSE(ves,identifiant)
{
  if(boite_3d==NULL) change_boite_analyse(ves->get_boite3d_ves());
  m_id_fem_solution_temp=id_fem_solution_temp;
  m_id_fem_solution_flux=id_fem_solution_flux;
  m_direction=direction;
  m_lambda_x=NULL;
  m_lambda_y=NULL;
  m_lambda_z=NULL;
}
 
 
MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE(std::string identifiant,
                                                                               MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE* lambda_x,
                                                                               MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE* lambda_y,
                                                                               MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE* lambda_z):MSTRUCT_ANALYSE(identifiant)
{
  m_id_fem_solution_temp=-1;
  m_id_fem_solution_flux=-1;
  m_direction=-1;
  m_lambda_x=lambda_x;
  m_lambda_y=lambda_y;
  m_lambda_z=lambda_z;
}
 
MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE(MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE& mdd, bool cpy_data):MSTRUCT_ANALYSE(mdd,cpy_data)
{
  m_id_fem_solution_temp=mdd.m_id_fem_solution_temp;
  m_id_fem_solution_flux=mdd.m_id_fem_solution_flux;
  m_direction=mdd.m_direction;
  m_lambda_x=mdd.m_lambda_x;
  m_lambda_y=mdd.m_lambda_y;
  m_lambda_z=mdd.m_lambda_z;
}
 
 
MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::~MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE()
{
}
 
long MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::get_type()
{
  return MSTRUCT::TYPE_ANALYSE::CONDUCTIVITE_THERMIQUE;
}
 
double MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::get_Lambda()
{
  return m_Lambda;
}
 
 
void MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::executer()
{
  if(m_lambda_x==NULL && m_lambda_y==NULL && m_lambda_z==NULL)
  {
    FEM_SOLUTION* fem_sol_flux = m_ves->get_mg_gestionnaire()->get_fem_solutionid(m_id_fem_solution_flux);
    LISTE_FEM_ELEMENT3::iterator it;
    long k=0;
    FEM_MAILLAGE* fem_mail = fem_sol_flux->get_maillage();
    
    double volume=0;
    double grad_temp=0;
    double flux=0;
    for(FEM_ELEMENT3*ele=fem_mail->get_premier_element3(it);ele!=NULL;ele=fem_mail->get_suivant_element3(it))
    {
      BOITE_3D boite_ele=ele->get_boite_3D();
      if(m_boite_analyse!=NULL)
        if(!m_boite_analyse->contient(boite_ele.get_xcentre(),boite_ele.get_ycentre(),boite_ele.get_zcentre()))
        {
          k++;
          continue;
        }
      MG_ELEMENT_TOPOLOGIQUE* topo=ele->get_lien_topologie();
      double lambda_mat;
      topo->get_valeur_ccf((char*)"Kx",lambda_mat);
      double unite=fem_mail->get_mg_geometrie()->get_valeur_unite();
      double val=0.;
      int nbgauss=fem_sol_flux->get_nb_gauss(k);
      for (int r=0;r<nbgauss;r++)
      {
        int degre=ele->get_degre_gauss(nbgauss);
        double w;
        double uvw[3];
        ele->get_pt_gauss(degre,r,w,uvw);
        double jac[9];
        double det=ele->get_jacobien(jac,uvw,unite);
        double valeur=fem_sol_flux->lire(k,m_direction,0,r)*fabs(det);
        val=val+w*valeur;          
        volume=volume+w*fabs(det);          
      }
      grad_temp+=-val/lambda_mat;
      flux+=val;
      k++;
    }
    grad_temp=grad_temp/volume;
    flux=flux/volume;
    m_Lambda=-flux/grad_temp;
 
  }
  else if(m_lambda_x!=NULL && m_lambda_y!=NULL && m_lambda_z!=NULL)
  {
    m_Lambda=(m_lambda_x->get_Lambda()+m_lambda_y->get_Lambda()+m_lambda_z->get_Lambda())/3.0;
  }
}
 
 
 
void MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::exporter(std::ofstream& ofstrm, long i, bool avec_entete, bool avec_graph, char* prefix_graph)
{
  if(avec_entete) ofstrm << "#(1) Lambda(2)" << std::endl;
  ofstrm << i << " " 
         << m_Lambda << " " << std::endl;   
}
 
void MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::enregistrer(std::ofstream& ofstrm)
{
  long type_analyse=get_type();
  ofstrm.write((char*)&type_analyse,sizeof(long));
  MSTRUCT_ANALYSE::enregistrer(ofstrm);
  ofstrm.write((char*)&m_Lambda,sizeof(double));
}
 
void MSTRUCT_ANALYSE_CONDUCTIVITE_THERMIQUE::ouvrir(std::ifstream& ifstrm)
{
  MSTRUCT_ANALYSE::ouvrir(ifstrm);
  ifstrm.read((char*)&m_Lambda,sizeof(double));
}