mstruct_analyse_energie_hill.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_energie_hill.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "mstruct_analyse_energie_hill.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_ENERGIE_HILL::MSTRUCT_ANALYSE_ENERGIE_HILL():MSTRUCT_ANALYSE()
{
    m_id_fem_epsilon=-1;
    m_id_fem_sigma=-1;
}
 
MSTRUCT_ANALYSE_ENERGIE_HILL::MSTRUCT_ANALYSE_ENERGIE_HILL(MSTRUCT_VES* ves, std::string identifiant, long id_fem_epsilon, long id_fem_sigma,BOITE_3D* boite_3d):MSTRUCT_ANALYSE(ves,identifiant)
{
  if(boite_3d==NULL) change_boite_analyse(ves->get_boite3d_ves());
  m_id_fem_epsilon=id_fem_epsilon;
  m_id_fem_sigma=id_fem_sigma;
}
 
 
MSTRUCT_ANALYSE_ENERGIE_HILL::MSTRUCT_ANALYSE_ENERGIE_HILL(MSTRUCT_ANALYSE_ENERGIE_HILL& mdd, bool cpy_data):MSTRUCT_ANALYSE(mdd,cpy_data)
{
  m_id_fem_epsilon=mdd.m_id_fem_epsilon;
  m_id_fem_sigma=mdd.m_id_fem_sigma;
}
 
MSTRUCT_ANALYSE_ENERGIE_HILL::~MSTRUCT_ANALYSE_ENERGIE_HILL()
{
}
 
double MSTRUCT_ANALYSE_ENERGIE_HILL::get_energie()
{
    return m_energie;
}
 
void MSTRUCT_ANALYSE_ENERGIE_HILL::exporter(std::ofstream& ofstrm, long i, bool avec_entete, bool avec_graph, char* prefix_graph)
{
    if(avec_entete) ofstrm << "#(1) Energie(2)" << std::endl;
  ofstrm << i << " " 
         << m_energie << " " << std::endl;   
}
 
long MSTRUCT_ANALYSE_ENERGIE_HILL::get_type()
{
    return MSTRUCT::TYPE_ANALYSE::ENERGIE_HILL;
}
 
void MSTRUCT_ANALYSE_ENERGIE_HILL::executer()
{
 
    double volume_total=0.0;
    m_energie=0.0;
    FEM_SOLUTION* fem_sol_epsilon = m_ves->get_mg_gestionnaire()->get_fem_solutionid(m_id_fem_epsilon);
    FEM_SOLUTION* fem_sol_sigma = m_ves->get_mg_gestionnaire()->get_fem_solutionid(m_id_fem_sigma);
    
    LISTE_FEM_ELEMENT3::iterator it;
    long k=0;
    FEM_MAILLAGE* fem_mail = fem_sol_epsilon->get_maillage();    
    int nb_champ = fem_sol_epsilon->get_nb_champ();
    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;
        }
        double epsilon[6];
        double sigma[6];
        double volume=0;
        for(int i=0;i<nb_champ;i++)
        {
            double unite=fem_mail->get_mg_geometrie()->get_valeur_unite();
            double val_eps=0.;
            double val_sigma=0.;
            volume=0;
            int nbgauss=fem_sol_epsilon->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_eps=fem_sol_epsilon->lire(k,i,0,r)*fabs(det);
                double valeur_sigma=fem_sol_sigma->lire(k,i,0,r)*fabs(det);
                val_eps=val_eps+w*valeur_eps;       
                val_sigma=val_sigma+w*valeur_sigma;  
                volume=volume+w*fabs(det);          
            }
            epsilon[i]=val_eps;
            sigma[i]=val_sigma;
        }
        k++;
 
      
      double energie=0.0;
      energie+=epsilon[0]*sigma[0];
      energie+=epsilon[1]*sigma[1];
      energie+=epsilon[2]*sigma[2];
      energie+=2*epsilon[3]*sigma[3];
      energie+=2*epsilon[4]*sigma[4];
      energie+=2*epsilon[5]*sigma[5];
      energie/=volume;
      
      m_energie+=energie;
      volume_total+=volume;
    }
    m_energie=m_energie/volume_total;
}
 
void MSTRUCT_ANALYSE_ENERGIE_HILL::enregistrer(std::ofstream& ofstrm)
{
    long type_analyse=get_type();
  ofstrm.write((char*)&type_analyse,sizeof(long));
  MSTRUCT_ANALYSE::enregistrer(ofstrm);
  ofstrm.write((char*)&m_energie,sizeof(double));
}
 
void MSTRUCT_ANALYSE_ENERGIE_HILL::ouvrir(std::ifstream& ifstrm)
{
    MSTRUCT_ANALYSE::ouvrir(ifstrm);
    ifstrm.read((char*)&m_energie,sizeof(double));
}