opt_triangle.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       opt_triangle.h
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#ifndef _OPTTRIANGLE_
#define _OPTTRIANGLE_
 
 
 
 
#include "mg_element_maillage.h"
#include "mg_segment.h"
#include "mg_noeud.h"
#include "opt_noeud.h"
#include "ot_reference.h"
 
 
class  OPT_TRIANGLE
{
public:
    
    OPT_TRIANGLE(class FEM_TRIANGLE3* fem_tri3,class OPT_NOEUD* noeud1,class OPT_NOEUD* noeud2,class OPT_NOEUD* noeud3);
    OPT_TRIANGLE(OPT_TRIANGLE& mdd);
    virtual ~OPT_TRIANGLE();
 
    virtual class OPT_NOEUD* get_noeud1(void);
    virtual class OPT_NOEUD* get_noeud2(void);
    virtual class OPT_NOEUD* get_noeud3(void);
    
    virtual unsigned long get_id(void);
    virtual void change_normale_unitaire(OT_VECTEUR_3D vec);
        
    virtual void change_angle1_initial(double angle);
    virtual void change_angle2_initial(double angle);
    virtual void change_angle3_initial(double angle);
    virtual void change_angle1(double angle);
    virtual void change_angle2(double angle);
    virtual void change_angle3(double angle);
    
    virtual OT_VECTEUR_3D get_normale_unitaire(void);   //Enlever?
    
    virtual OT_MATRICE_3D get_jacobien_inverse(void);
    virtual OT_MATRICE_3D get_matrice_abc(void);
    
    virtual OT_MATRICE_3D change_jacobien_inverse(void);
    virtual OT_MATRICE_3D change_matrice_abc(void);
    
    //virtual OT_MATRICE_3D get_jacobien_inverse_initial(void);
    //virtual OT_MATRICE_3D get_matrice_abc_initiale(void);
    
    virtual OT_VECTEUR_3D get_vecteur_deplacement_reel(void);
    virtual void change_vecteur_deplacement_reel(void);
    
    virtual OT_VECTEUR_3D get_vecteur_sol_deplacement_virtuel(void);
    virtual void change_vecteur_sol_deplacement_virtuel(void);
    
    virtual double get_norme_gradient_deplacement_reel(void);
    virtual void change_norme_gradient_deplacement_reel(void);
    
    virtual double get_norme_gradient_deplacement_virtuel(void);
    virtual void change_norme_gradient_deplacement_virtuel(void);
    
    virtual int get_noeuds_fixes_3(void);
    virtual void change_noeuds_fixes_3(int val);
    
    virtual int get_nb_noeuds_fixes_iter(void);
    virtual void change_nb_noeuds_fixes_iter(int nb);
    
    virtual double get_angle1_initial(void);
    virtual double get_angle2_initial(void);
    virtual double get_angle3_initial(void);
    virtual double get_angle1(void);
    virtual double get_angle2(void);
    virtual double get_angle3(void);
    
    
    virtual int get_num(void);
    virtual void change_num(int val);
    
    virtual OPT_TRIANGLE* get_voisin1(void);
    virtual OPT_TRIANGLE* get_voisin2(void);
    virtual OPT_TRIANGLE* get_voisin3(void);
    virtual void change_voisin1(OPT_TRIANGLE* tri);
    virtual void change_voisin2(OPT_TRIANGLE* tri);
    virtual void change_voisin3(OPT_TRIANGLE* tri);
    
    
protected:
    
    class FEM_TRIANGLE3* fem_triangle3;
    OPT_NOEUD* opt_noeud1;
    OPT_NOEUD* opt_noeud2;
    OPT_NOEUD* opt_noeud3;
    
    OT_VECTEUR_3D normale_unitaire;
    OT_MATRICE_3D dNi_transpose;
    OT_MATRICE_3D jacobien_inverse;
    OT_MATRICE_3D matrice_abc;                  // Matrice résultant du produit de la matrice jacobienne inverse par la matrice des dérivées des fonctions d'interpolation
    OT_VECTEUR_3D vec_deplacement_reel; // Vecteur contenant les déplacements des noeuds par rapport à leur position initiale
    OT_VECTEUR_3D vec_sol_deplacement_virtuel;  // Vecteur contenant les solutions de déplacement virtuel des noeuds
    double norme_gradient_deplacement_reel;
    double norme_gradient_deplacement_virtuel;
    int noeuds_fixes_3;         // Si noeuds_fixes_3==1, le triangle a 3 noeuds fixes
    int nb_noeuds_fixes_iter;   // Nb de noeuds fixés au cours d'une même itération de la méthode du mvt normal
    
    int num;
    double a1_init;
    double a2_init;
    double a3_init;
    double a1;
    double a2;
    double a3;
    OPT_TRIANGLE* voisin1;
    OPT_TRIANGLE* voisin2;
    OPT_TRIANGLE* voisin3;
};
 
#endif