mstruct_analyse_modules_elasticite.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_modules_elasticite.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "mstruct_analyse_modules_elasticite.h"
#include "mstruct_analyse_champ.h"
#include "ot_boite_3d.h"
#include "mstruct_definition.h"
 
MSTRUCT_ANALYSE_MODULES_ELASTICITE::MSTRUCT_ANALYSE_MODULES_ELASTICITE(void):MSTRUCT_ANALYSE()
{
}
 
 
MSTRUCT_ANALYSE_MODULES_ELASTICITE::MSTRUCT_ANALYSE_MODULES_ELASTICITE(std::string identifiant,
                                                                       MSTRUCT_ANALYSE_CHAMP* epsilon_sph,
                                                                       MSTRUCT_ANALYSE_CHAMP* sigma_sph,
                                                                       MSTRUCT_ANALYSE_CHAMP* epsilon_dev,
                                                                       MSTRUCT_ANALYSE_CHAMP* sigma_dev): MSTRUCT_ANALYSE(identifiant)
{
  if(epsilon_sph->get_boite_analyse()!=NULL) change_boite_analyse(*epsilon_sph->get_boite_analyse());
  change_nom_groupe_forme(epsilon_sph->get_nom_groupe_forme());  
  m_epsilon_sph=epsilon_sph;
  m_sigma_sph=sigma_sph;
  m_epsilon_dev=epsilon_dev;
  m_sigma_dev=sigma_dev;
  m_Kapp=0;
  m_Gapp=0;
  m_Eapp=0;
  m_Nuapp=0;
}
 
MSTRUCT_ANALYSE_MODULES_ELASTICITE::~MSTRUCT_ANALYSE_MODULES_ELASTICITE(void)
{
}
 
double MSTRUCT_ANALYSE_MODULES_ELASTICITE::get_Kapp(void)
{
  return m_Kapp;
}
 
double MSTRUCT_ANALYSE_MODULES_ELASTICITE::get_Gapp(void)
{
  return m_Gapp;
}
 
double MSTRUCT_ANALYSE_MODULES_ELASTICITE::get_Nuapp(void)
{
  return m_Nuapp;
}
 
double MSTRUCT_ANALYSE_MODULES_ELASTICITE::get_Eapp(void)
{
  return m_Eapp;
}
 
long MSTRUCT_ANALYSE_MODULES_ELASTICITE::get_type()
{
  return MSTRUCT::TYPE_ANALYSE::MODULES_ELASTICITE;
}
 
void MSTRUCT_ANALYSE_MODULES_ELASTICITE::executer()
{
  m_Kapp = (m_sigma_sph->get_moyenne()[0]+m_sigma_sph->get_moyenne()[1]+m_sigma_sph->get_moyenne()[2])/(3.0*(m_epsilon_sph->get_moyenne()[0]+m_epsilon_sph->get_moyenne()[1]+m_epsilon_sph->get_moyenne()[2]));
  m_Gapp = (1./3.)*((m_sigma_dev->get_moyenne()[3]/(2.0*m_epsilon_dev->get_moyenne()[3]))+
                    (m_sigma_dev->get_moyenne()[4]/(2.0*m_epsilon_dev->get_moyenne()[4]))+
                    (m_sigma_dev->get_moyenne()[5]/(2.0*m_epsilon_dev->get_moyenne()[5])));
  m_Eapp = (9.0*m_Kapp*m_Gapp)/(3.0*m_Kapp+m_Gapp);
  m_Nuapp = (3.0*m_Kapp-2.0*m_Gapp)/(2.0*(3.0*m_Kapp+m_Gapp));
}
 
 
void MSTRUCT_ANALYSE_MODULES_ELASTICITE::exporter(std::ofstream& ofstrm, long i, bool avec_entete, bool avec_graph, char* prefix_graph)
{
  if(avec_entete) ofstrm << "#(1) Eapp(2) Nuapp(3) Gapp(4) Kapp(5)" << std::endl;
  ofstrm << i << " " 
         << m_Eapp << " "
         << m_Nuapp << " "
         << m_Gapp << " "
         << m_Kapp << " " << std::endl;    
}
 
void MSTRUCT_ANALYSE_MODULES_ELASTICITE::enregistrer(std::ofstream& ofstrm)
{
  long type_analyse=get_type();
  ofstrm.write((char*)&type_analyse,sizeof(long));
  MSTRUCT_ANALYSE::enregistrer(ofstrm);
  ofstrm.write((char*)&m_Eapp,sizeof(double));
  ofstrm.write((char*)&m_Nuapp,sizeof(double));
  ofstrm.write((char*)&m_Gapp,sizeof(double));
  ofstrm.write((char*)&m_Kapp,sizeof(double));
}
 
void MSTRUCT_ANALYSE_MODULES_ELASTICITE::ouvrir(std::ifstream& ifstrm)
{
  MSTRUCT_ANALYSE::ouvrir(ifstrm);
  ifstrm.read((char*)&m_Eapp,sizeof(double));
  ifstrm.read((char*)&m_Nuapp,sizeof(double));
  ifstrm.read((char*)&m_Gapp,sizeof(double));
  ifstrm.read((char*)&m_Kapp,sizeof(double));
}