mgopt_mvt_normal.h

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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       mgopt_mvt_normal.h
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#ifndef mgoptmvtnormalH
#define mgoptmvtnormalH
 
 
 
#include "tpl_map_entite.h"
#include "ot_mathematique.h"
#include "ot_parametres.h"
#include "mg_gestionnaire.h"
#include "opt_noeud.h"
#include "fem_triangle3.h"
 
class OPT_NOEUD;
class OPT_TRIANGLE;
class FEM_TRIANGLE3;
 
class  MGOPT_MVT_NORMAL
{
public :
    
    MGOPT_MVT_NORMAL();
    MGOPT_MVT_NORMAL(class FEM_MAILLAGE* fem, class FEM_MAILLAGE* fem_visu);
    MGOPT_MVT_NORMAL(MGOPT_MVT_NORMAL& mdd);
    ~MGOPT_MVT_NORMAL();
 
    virtual void active_affichage(void (*fonc)(char*));
    virtual void init_parametre(void);
    virtual void ecrire_fichier_params(char *fichierparam);
    virtual void get_liste_noeud_deplacer_initiale(TPL_MAP_ENTITE< OPT_NOEUD* > listenoeud, TPL_MAP_ENTITE< OPT_NOEUD* >* lstnoeud_deplacer_initiale);
    virtual void lissage(TPL_MAP_ENTITE<OPT_NOEUD*> lst_noeud);
    virtual void lissage_controle(TPL_MAP_ENTITE<OPT_NOEUD*> lst_noeud, double ecart_max_pondere_norme_pos);
    virtual void optimisation(char* fichierin, char* fichierout, char* fichierparam, char* fichierparamaster);
    //virtual void lissage(char* fichier);
    //virtual void deplacement_normal(double dpl_max, double angle_min, int nbpas);
    //virtual void visualisation(int facteur);
    
    // Correspondance entre opt_noeud et fem_noeud
    TPL_MAP_ENTITE<class OPT_NOEUD*> listenoeud;
    class OPT_NOEUD* get_noeud(class FEM_NOEUD* fem_noeud);
    std::map<unsigned long, std::pair<class OPT_NOEUD*,class FEM_NOEUD*> > correspondance_noeud;
    
    void printf(const char* arg1);
    std::ostream fem_tri3();
    
protected:
    
    // Fonctions essentielles méthode du mouvement normal
    virtual void lire_params(char *fichier);
    //virtual void get_liste_noeud_deplacer_initiale(TPL_MAP_ENTITE< OPT_NOEUD* > listenoeud, TPL_MAP_ENTITE< OPT_NOEUD* >* lstnoeud_deplacer_initiale);
    //virtual void change_liste_noeud_voisin(TPL_MAP_ENTITE<OPT_NOEUD*> lstnoeud_deplacer_initiale);
    virtual void change_liste_noeud_voisin(OPT_NOEUD* opt_noeud);
    virtual double get_aire_fem_tri3(FEM_TRIANGLE3* fem_tri);
    virtual void calcul_normale_opt_noeud(class OPT_NOEUD* noeud);
    virtual double get_sigma_vm_max(MG_GESTIONNAIRE* gest, FEM_NOEUD* fem_noeud, int num_sol_inf, int num_sol_moy, int num_sol_sup);
    virtual double get_ecart_type_sigma_vm(MG_GESTIONNAIRE* gest, TPL_MAP_ENTITE<OPT_NOEUD*> lstnoeud, double moyenne, int numsolinf, int numsolmoy, int numsolsup);
    
    // Fonctions pour la gestion des déplacement virtuels et réels
    virtual void change_deplacement_virtuel_opt_noeud(OPT_NOEUD* opt_noeud, double d);
    virtual void deplace_opt_noeud(OPT_NOEUD* opt_noeud, double d);
    
    // Fonctions de visualisation des déplacements
    virtual void ecriture_fem_solution_deplacement(FEM_SOLUTION* sol, FEM_NOEUD* fem_noeud);
    
    // Fonctions pour le contrôle par gradient de déplacement
    virtual void get_liste_tri_2nd_fixe(FEM_MAILLAGE* fem,TPL_MAP_ENTITE<OPT_NOEUD*> lstnoeud_deplacer,TPL_MAP_ENTITE<class OPT_TRIANGLE*>* lst_tri_2nd_fixe);
    virtual void get_liste_noeud_a_partir_liste_tri(TPL_MAP_ENTITE< OPT_TRIANGLE* > lst_tri, TPL_MAP_ENTITE< OPT_NOEUD* >* lst_noeud);
    virtual void get_liste_noeud_mobile_tri_2nd_fixe(TPL_MAP_ENTITE<OPT_TRIANGLE*> lst_tri_2nd_fixe, TPL_MAP_ENTITE< OPT_NOEUD* > lstnoeud_deplacer, TPL_MAP_ENTITE<OPT_NOEUD*> *lst_noeud_mobile_tri_2nd_fixe);
    //virtual void get_liste_opt_tri_3nd_fixe(FEM_MAILLAGE* fem, TPL_MAP_ENTITE<OPT_NOEUD*> lstnoeud_deplacer, TPL_MAP_ENTITE<OPT_TRIANGLE*>* lst_opt_tri_3nd_fixe);    
    //virtual void get_noeud_mobile_opt_tri(OPT_TRIANGLE* opt_tri, TPL_MAP_ENTITE<OPT_NOEUD*> lstnoeud_deplacer, OPT_NOEUD** opt_noeud_mobile);
    virtual double evaluer_fonction_f(double norme_max_grad, TPL_MAP_ENTITE< OPT_TRIANGLE* > lst_tri, int affichage);//, TPL_MAP_ENTITE<OPT_NOEUD*> lst_noeud_mobile_tri_2nd_fixe);
    virtual void evaluer_gradient_f(int num_repl, double norme_max_grad, TPL_MAP_ENTITE<OPT_TRIANGLE*> lst_tri);
    virtual double get_norme_gradient_f(int num_repl, TPL_MAP_ENTITE<OPT_NOEUD*> lst_noeud);
    virtual int algorithme_gradient(int num_repl, double norme_max_grad, double pas, double epsilon, int iter_max_grad, TPL_MAP_ENTITE<OPT_TRIANGLE*> lst_tri, TPL_MAP_ENTITE<OPT_NOEUD*> lst_noeud, int affichage);
    
    // Fonction non utilisée pour calculer l'angle entre deux triangles
    //virtual double calcul_angle(class OPT_TRIANGLE* tri1,class OPT_TRIANGLE* tri2);
    
    // Correspondance entre opt_triangle et fem_triangle
    TPL_MAP_ENTITE<class OPT_TRIANGLE*> listetriangle;
    class OPT_TRIANGLE* get_triangle(class FEM_TRIANGLE3* fem_tri3);
    std::map<unsigned long, std::pair<class OPT_TRIANGLE*,class FEM_TRIANGLE3*> > correspondance_triangle;
    
    class MG_MAILLAGE *mai;
    class FEM_MAILLAGE *fem;            // Maillage servant pour l'optimisation
    class FEM_MAILLAGE* fem_visu;       // Maillage servant pour la visualisation
    class MG_GEOMETRIE *geo;
    int orientation;
    
    void (*affiche)(char *mess);
    int affichageactif;
    OT_PARAMETRES params;
    char etude[500];
    char etudesortie[500];
    
    int nb_tri_grad_depasse;    // Nb de triangles dont la norme du gradient de déplacement dépasse la valeur seuil lors de l'écriture de la fonction f
    //int nb_tri_3noeud_fixe;   // Nb de triangles qui ont 3 noeuds fixes et qui avaient au moins 1 noeud mobile au départ
};
 
#endif