fct_taille_fem_solution_generateur_gradient.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_gradient.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "gestionversion.h"
#include "fct_taille_fem_solution_generateur_gradient.h"
#include "ot_geometrie.h"
#include "tpl_grille.h"
 
 FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT(MG_GESTIONNAIRE *gesttmp,                                                   MG_GEOMETRIE* geotmp,double dg,int fechantillon,int nnx,int nny,int nnz,char *nomsolution,std::vector<double> lst,double gradient):FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(gesttmp,geotmp,dg,fechantillon,nnx,nny,nnz,nomsolution),gradientmax(gradient),lstraffine(lst)
 {
     
 }
                                                     
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT(MG_GESTIONNAIRE *gesttmp,                                                    MG_GEOMETRIE* geotmp,FCT_TAILLE* carte,int fechantillon,int nnx,int nny,int nnz,char *nomsolution,std::vector<double> lst,double gradient):FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(gesttmp,geotmp,carte,fechantillon,nnx,nny,nnz,nomsolution),gradientmax(gradient),lstraffine(lst)
 {
     
 }
 
 FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT(MG_GESTIONNAIRE *gesttmp,                                                   MG_MAILLAGE* maitmp,double dg,int fechantillon,int nnx,int nny,int nnz,char *nomsolution,std::vector<double> lst,double gradient):FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(gesttmp,maitmp,dg,fechantillon,nnx,nny,nnz,nomsolution),gradientmax(gradient),lstraffine(lst)
 {
     
 }
                                                     
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT(MG_GESTIONNAIRE *gesttmp,                                                    MG_MAILLAGE* maitmp,FCT_TAILLE* carte,int fechantillon,int nnx,int nny,int nnz,char *nomsolution,std::vector<double> lst,double gradient):FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(gesttmp,maitmp,carte,fechantillon,nnx,nny,nnz,nomsolution),gradientmax(gradient),lstraffine(lst)
 {
     
 }
                                                     
                                                     
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT(FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT& mdd):FCT_TAILLE_FEM_SOLUTION_GENERATEUR_ISOTROPE(mdd),gradientmax(mdd.gradientmax),lstraffine(mdd.lstraffine)
{
}
                                                     
                                                     
                                                     
                                                     
 
 
FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::~FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT()
{
}
 
 
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::construit_particulier()
{
cree_grille();
affiche((char*)"    Points a raffiner");
echantillonnage();
affiche((char*)"    Lissage");
lisser_solution();
affiche((char*)"    Extrapolation aux noeuds");
extrapoler_solution_noeud();
}
 
int FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::echantillonnage()
{
    LISTE_FEM_ELEMENT3::iterator it_ele3;
    for(FEM_ELEMENT3* ele3=fem->get_premier_element3(it_ele3);ele3!=NULL;ele3=fem->get_suivant_element3(it_ele3))
        if (carteini==NULL) ele3->change_solution(ecart_nodal);
        else
            {
             double xyz[3];
             ele3->get_boite_3D().get_centre(xyz);
             double tab[9];
             carteini->evaluer(xyz,tab);
             double ecart=1./sqrt(tab[0]);
             ele3->change_solution(ecart);
            }
    int nb=lstraffine.size();
    for (int i=0;i<nb/4;i=i+1)
            {
            double x=lstraffine[4*i];    
            double y=lstraffine[4*i+1];    
            double z=lstraffine[4*i+2];    
            double val=lstraffine[4*i+3];    
            TPL_CELLULE_GRILLE<FEM_ELEMENT3*> *cell=m_grille_ele3.get_cellule(x,y,z);
            FEM_ELEMENT3* ele3=cell->get_entite(0);
            double valcell=ele3->get_solution();
            if (val<valcell)
                {
                ele3->change_solution(val);
                std::pair<FEM_ELEMENT3*,FEM_ELEMENT3*> tmp(ele3,ele3);
                m_map_fem_ele3.insert(tmp);
                }
            }
 
 
    std::map<FEM_ELEMENT3*,FEM_ELEMENT3*>::iterator it_map_ele3;
    for(it_map_ele3=m_map_fem_ele3.begin();it_map_ele3!=m_map_fem_ele3.end();it_map_ele3++)
    {
      FEM_ELEMENT3* ele3=it_map_ele3->second;
      double val=ele3->get_solution();
      std::pair<double,FEM_ELEMENT3*> tmp(val,ele3);
      m_map_fem_classe.insert(tmp);
    }  
return 0;
}
 
 
 
 
 
 
 
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::lisser_solution(void)
{
 
  while(!m_map_fem_classe.empty())
  {
    FEM_ELEMENT3* ele3=m_map_fem_classe.begin()->second;
    double val=ele3->get_solution();
    double centre[3];
    ele3->get_boite_3D().get_centre(centre);
    double en_centre=ele3->get_solution(0);
    TPL_CELLULE_GRILLE<FEM_ELEMENT3*> *tab_cell[26];
    for(int i=0;i<26;i++)tab_cell[i]=NULL;
    
    tab_cell[0] = m_grille_ele3.get_cellule(centre[0]-pasx,centre[1]-pasy,centre[2]-pasz);
    tab_cell[1] = m_grille_ele3.get_cellule(centre[0],centre[1]-pasy,centre[2]-pasz);
    tab_cell[2] = m_grille_ele3.get_cellule(centre[0]+pasx,centre[1]-pasy,centre[2]-pasz);
    tab_cell[3] = m_grille_ele3.get_cellule(centre[0]-pasx,centre[1],centre[2]-pasz);
    tab_cell[4] = m_grille_ele3.get_cellule(centre[0],centre[1],centre[2]-pasz);
    tab_cell[5] = m_grille_ele3.get_cellule(centre[0]+pasx,centre[1],centre[2]-pasz);
    tab_cell[6] = m_grille_ele3.get_cellule(centre[0]-pasx,centre[1]+pasy,centre[2]-pasz);
    tab_cell[7] = m_grille_ele3.get_cellule(centre[0],centre[1]+pasy,centre[2]-pasz);
    tab_cell[8] = m_grille_ele3.get_cellule(centre[0]+pasx,centre[1]+pasy,centre[2]-pasz);
    
    tab_cell[9] = m_grille_ele3.get_cellule(centre[0]-pasx,centre[1]-pasy,centre[2]);
    tab_cell[10] = m_grille_ele3.get_cellule(centre[0],centre[1]-pasy,centre[2]);
    tab_cell[11] = m_grille_ele3.get_cellule(centre[0]+pasx,centre[1]-pasy,centre[2]);
    tab_cell[12] = m_grille_ele3.get_cellule(centre[0]-pasx,centre[1],centre[2]);
    tab_cell[13] = m_grille_ele3.get_cellule(centre[0]+pasx,centre[1],centre[2]);
    tab_cell[14] = m_grille_ele3.get_cellule(centre[0]-pasx,centre[1]+pasy,centre[2]);
    tab_cell[15] = m_grille_ele3.get_cellule(centre[0],centre[1]+pasy,centre[2]);
    tab_cell[16] = m_grille_ele3.get_cellule(centre[0]+pasx,centre[1]+pasy,centre[2]);
    
    tab_cell[17] = m_grille_ele3.get_cellule(centre[0]-pasx,centre[1]-pasy,centre[2]+pasz);
    tab_cell[18] = m_grille_ele3.get_cellule(centre[0],centre[1]-pasy,centre[2]+pasz);
    tab_cell[19] = m_grille_ele3.get_cellule(centre[0]+pasx,centre[1]-pasy,centre[2]+pasz);
    tab_cell[20] = m_grille_ele3.get_cellule(centre[0]-pasx,centre[1],centre[2]+pasz);
    tab_cell[21] = m_grille_ele3.get_cellule(centre[0],centre[1],centre[2]+pasz);
    tab_cell[22] = m_grille_ele3.get_cellule(centre[0]+pasx,centre[1],centre[2]+pasz);
    tab_cell[23] = m_grille_ele3.get_cellule(centre[0]-pasx,centre[1]+pasy,centre[2]+pasz);
    tab_cell[24] = m_grille_ele3.get_cellule(centre[0],centre[1]+pasy,centre[2]+pasz);
    tab_cell[25] = m_grille_ele3.get_cellule(centre[0]+pasx,centre[1]+pasy,centre[2]+pasz);
    for(int i=0;i<26;i++)
    {
      if(tab_cell[i]!=NULL)
        if(tab_cell[i]->get_nb_entite()==1)
            {
            FEM_ELEMENT3* ele3_voisin = tab_cell[i]->get_entite(0);
            double valvoisin=ele3_voisin->get_solution();
            double centrevoisin[3];
            ele3_voisin->get_boite_3D().get_centre(centrevoisin);
            OT_VECTEUR_3D vec(centre,centrevoisin);
            double valvoisinmax=val+gradientmax*vec.get_longueur();
            if(valvoisinmax<valvoisin)
                    {              
                    ele3_voisin->change_solution(valvoisinmax); 
                    std::map<double,FEM_ELEMENT3*>::iterator j=m_map_fem_classe.lower_bound(valvoisin);
                    int ok=0;
                    if (j!=m_map_fem_classe.end())
                      while (ok==0)
                        {
                        if ((*j).first>valvoisin) ok=1;
                        if (ele3_voisin==(*j).second) {m_map_fem_classe.erase(j);ok=1;}
                        if (ok==0) j++;    
                        }
                    std::pair<double,FEM_ELEMENT3*> tmp(valvoisinmax,ele3_voisin);
                    m_map_fem_classe.insert(tmp);    
                    }
            }
    }
   m_map_fem_classe.erase(m_map_fem_classe.begin());
  }  
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::extrapoler_solution_noeud(void)
{
  int i=0;
  LISTE_FEM_NOEUD::iterator it_fem_nd;
  for(FEM_NOEUD* femnd=fem->get_premier_noeud(it_fem_nd);femnd!=NULL;femnd=fem->get_suivant_noeud(it_fem_nd))
  {
    double ecart_nodal_noeud=0.0;
    TPL_LISTE_ENTITE<FEM_ELEMENT3*> *liste_fem_element3;
    liste_fem_element3 = femnd->get_lien_element3();    
    for(int i=0;i<liste_fem_element3->get_nb();i++)
    {
      ecart_nodal_noeud+=liste_fem_element3->get(i)->get_solution(0);
    }
    ecart_nodal_noeud=ecart_nodal_noeud/liste_fem_element3->get_nb();
  /*  femnd->change_solution(1.0/ecart_nodal_noeud/ecart_nodal_noeud,0);
    femnd->change_solution(1.0/ecart_nodal_noeud/ecart_nodal_noeud,4);
    femnd->change_solution(1.0/ecart_nodal_noeud/ecart_nodal_noeud,8);
    femnd->change_solution(0.,1);
    femnd->change_solution(0.,2);
    femnd->change_solution(0.,3);
    femnd->change_solution(0.,5);
    femnd->change_solution(0.,6);
    femnd->change_solution(0.,7);*/
    sol->ecrire(1.0/ecart_nodal_noeud/ecart_nodal_noeud,i,0,0);
    sol->ecrire(0.0,i,0,1);
    sol->ecrire(0.0,i,0,2);
    sol->ecrire(0.0,i,0,3);
    sol->ecrire(1./ecart_nodal_noeud/ecart_nodal_noeud,i,0,4);
    sol->ecrire(0.0,i,0,5);
    sol->ecrire(0.0,i,0,6);
    sol->ecrire(0.0,i,0,7);
    sol->ecrire(1./ecart_nodal_noeud/ecart_nodal_noeud,i,0,8);
    i++;
    
  }
}
 
void FCT_TAILLE_FEM_SOLUTION_GENERATEUR_GRADIENT::cree_grille(void)
{
    double xmin=boite_englobante.get_xmin();
    double ymin=boite_englobante.get_ymin();
    double zmin=boite_englobante.get_zmin();
    double xmax=boite_englobante.get_xmax();
    double ymax=boite_englobante.get_ymax();
    double zmax=boite_englobante.get_zmax();
       
    m_grille_ele3.initialiser(xmin,
                          ymin,
                          zmin,
                          xmax,
                          ymax,
                          zmax,
                          nx,ny,nz);
  LISTE_FEM_ELEMENT3::iterator it_ele3;
  for(FEM_ELEMENT3* ele3=fem->get_premier_element3(it_ele3);ele3!=NULL;ele3=fem->get_suivant_element3(it_ele3))
  {
    double centre[3];
    ele3->get_boite_3D().get_centre(centre);
    TPL_CELLULE_GRILLE<FEM_ELEMENT3*>* cellule = m_grille_ele3.get_cellule(centre[0],centre[1],centre[2]);
    cellule->lst_entite_A.ajouter(ele3);
  }  
}