opt_triangle.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       opt_triangle.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:57 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include "mg_maillage.h"
#include "ot_decalage_parametre.h"
#include "opt_triangle.h"
#include "fem_triangle3.h"
#include "opt_noeud.h"
 
 
OPT_TRIANGLE::OPT_TRIANGLE(FEM_TRIANGLE3* fem_tri3,OPT_NOEUD* noeud1,OPT_NOEUD* noeud2,OPT_NOEUD* noeud3):fem_triangle3(fem_tri3),opt_noeud1(noeud1),opt_noeud2(noeud2),opt_noeud3(noeud3)
{
    double uv[6];
    
    OT_VECTEUR_3D dNi_de(fem_tri3->get_fonction_derive_interpolation(1,1,uv),fem_tri3->get_fonction_derive_interpolation(2,1,uv),fem_tri3->get_fonction_derive_interpolation(3,1,uv));
    OT_VECTEUR_3D dNi_dn(fem_tri3->get_fonction_derive_interpolation(1,2,uv),fem_tri3->get_fonction_derive_interpolation(2,2,uv),fem_tri3->get_fonction_derive_interpolation(3,2,uv));
    OT_VECTEUR_3D dNi_dZ(0.,0.,0.);
    
    OT_MATRICE_3D dNi(dNi_de,dNi_dn,dNi_dZ);
    dNi_transpose=dNi.transpose();
}
 
OPT_TRIANGLE::OPT_TRIANGLE(OPT_TRIANGLE& mdd):fem_triangle3(mdd.fem_triangle3),opt_noeud1(mdd.opt_noeud1),opt_noeud2(mdd.opt_noeud2),opt_noeud3(mdd.opt_noeud3)
{
    
}
 
OPT_TRIANGLE::~OPT_TRIANGLE()
{
    
}
 
unsigned long OPT_TRIANGLE::get_id(void)
{
    return fem_triangle3->get_id();
}
 
void OPT_TRIANGLE::change_angle1_initial(double angle)
{
    a1_init=angle;
}
 
void OPT_TRIANGLE::change_angle2_initial(double angle)
{
    a2_init=angle;
}
 
void  OPT_TRIANGLE::change_angle3_initial(double angle)
{
    a3_init=angle;
}
 
void OPT_TRIANGLE::change_angle1(double angle)
{
    a1=angle;
}
 
void OPT_TRIANGLE::change_angle2(double angle)
{
    a2=angle;
}
 
void OPT_TRIANGLE::change_angle3(double angle)
{
    a3=angle;
}
 
void OPT_TRIANGLE::change_normale_unitaire(OT_VECTEUR_3D vec)
{
    normale_unitaire=vec;
}
 
OT_VECTEUR_3D OPT_TRIANGLE::get_normale_unitaire(void)
{
    return normale_unitaire;
}
 
OT_MATRICE_3D OPT_TRIANGLE::change_jacobien_inverse(void)
{
    double j[9];
    double uv[6];
    double unite;
    
    fem_triangle3->get_inverse_jacob(j,uv,unite);
    
    OT_VECTEUR_3D col1_ji(j[0],j[3],j[6]);
    OT_VECTEUR_3D col2_ji(j[1],j[4],j[7]);
    OT_VECTEUR_3D col3_ji(j[2],j[5],j[8]);
    
    OT_MATRICE_3D ji(col1_ji,col2_ji,col3_ji);
    jacobien_inverse=ji;
    return jacobien_inverse;
}
 
OT_MATRICE_3D OPT_TRIANGLE::get_jacobien_inverse(void)
{
    return jacobien_inverse;
}
 
OT_MATRICE_3D OPT_TRIANGLE::change_matrice_abc(void)
{
    matrice_abc=jacobien_inverse*dNi_transpose;
    return matrice_abc;
}
 
OT_MATRICE_3D OPT_TRIANGLE::get_matrice_abc(void)
{
    return matrice_abc;
}
 
void OPT_TRIANGLE::change_vecteur_deplacement_reel(void)
{
    OPT_NOEUD* opt_n1=get_noeud1();
    OPT_NOEUD* opt_n2=get_noeud2();
    OPT_NOEUD* opt_n3=get_noeud3();
    
    double d1=opt_n1->get_norme_orientee_deplacement();
    double d2=opt_n2->get_norme_orientee_deplacement();
    double d3=opt_n3->get_norme_orientee_deplacement();
    
    OT_VECTEUR_3D D(d1,d2,d3);
    vec_deplacement_reel=D;
}
 
OT_VECTEUR_3D OPT_TRIANGLE::get_vecteur_deplacement_reel(void)
{
    return vec_deplacement_reel;
}
 
void OPT_TRIANGLE::change_vecteur_sol_deplacement_virtuel(void )
{
    OPT_NOEUD* opt_n1=get_noeud1();
    OPT_NOEUD* opt_n2=get_noeud2();
    OPT_NOEUD* opt_n3=get_noeud3();
    
    double d1=opt_n1->get_sol_deplacement_virtuel();
    double d2=opt_n2->get_sol_deplacement_virtuel();
    double d3=opt_n3->get_sol_deplacement_virtuel();
    
    OT_VECTEUR_3D D(d1,d2,d3);
    vec_sol_deplacement_virtuel=D;
}
 
OT_VECTEUR_3D OPT_TRIANGLE::get_vecteur_sol_deplacement_virtuel(void )
{
    return vec_sol_deplacement_virtuel;
}
 
void OPT_TRIANGLE::change_norme_gradient_deplacement_reel(void)
{
    OT_VECTEUR_3D gradient_deplacement_reel=matrice_abc*vec_deplacement_reel;
    norme_gradient_deplacement_reel=gradient_deplacement_reel.get_longueur();
}
 
double OPT_TRIANGLE::get_norme_gradient_deplacement_reel(void)
{
    return norme_gradient_deplacement_reel;
}
 
void OPT_TRIANGLE::change_norme_gradient_deplacement_virtuel(void)
{
    /*OPT_NOEUD* opt_n1=get_noeud1();
    OPT_NOEUD* opt_n2=get_noeud2();
    OPT_NOEUD* opt_n3=get_noeud3();
 
    OT_MATRICE_3D abc=get_matrice_abc();
 
    double d1=opt_n1->get_sol_deplacement_virtuel();
    double d2=opt_n2->get_sol_deplacement_virtuel();
    double d3=opt_n3->get_sol_deplacement_virtuel();
 
    OT_VECTEUR_3D D(d1,d2,d3);
      
    OT_VECTEUR_3D gradient_deplacement=abc*D;*/
    OT_VECTEUR_3D gradient_deplacement_virtuel=matrice_abc*vec_sol_deplacement_virtuel;
    norme_gradient_deplacement_virtuel=gradient_deplacement_virtuel.get_longueur();
}
 
double OPT_TRIANGLE::get_norme_gradient_deplacement_virtuel(void)
{
    return norme_gradient_deplacement_virtuel;
}
 
void OPT_TRIANGLE::change_noeuds_fixes_3(int val)
{
    noeuds_fixes_3=val;
}
 
int OPT_TRIANGLE::get_noeuds_fixes_3(void)
{
    return noeuds_fixes_3;
}
 
void OPT_TRIANGLE::change_nb_noeuds_fixes_iter(int nb)
{
    nb_noeuds_fixes_iter=nb;
}
 
int OPT_TRIANGLE::get_nb_noeuds_fixes_iter(void)
{
    return nb_noeuds_fixes_iter;
}
 
double OPT_TRIANGLE::get_angle1_initial(void)
{
    return a1_init;
}
 
double OPT_TRIANGLE::get_angle2_initial(void)
{
    return a2_init;
}
 
double OPT_TRIANGLE::get_angle3_initial(void)
{
    return a3_init;
}
 
double OPT_TRIANGLE::get_angle1(void)
{
    return a1;
}
 
double OPT_TRIANGLE::get_angle2(void)
{
    return a2;
}
 
double OPT_TRIANGLE::get_angle3(void)
{
    return a3;
}
 
OPT_NOEUD* OPT_TRIANGLE::get_noeud1(void)
{
    return opt_noeud1;
}
 
OPT_NOEUD* OPT_TRIANGLE::get_noeud2(void)
{
    return opt_noeud2;
}
 
OPT_NOEUD* OPT_TRIANGLE::get_noeud3(void)
{
     return opt_noeud3;
}
 
int OPT_TRIANGLE::get_num(void)
{
    return num;
}
 
void OPT_TRIANGLE::change_num(int val)
{
    num=val;
}
 
OPT_TRIANGLE* OPT_TRIANGLE::get_voisin1(void)
{
    return voisin1;
}
 
OPT_TRIANGLE* OPT_TRIANGLE::get_voisin2(void)
{
    return voisin2;
}
 
OPT_TRIANGLE* OPT_TRIANGLE::get_voisin3(void)
{
    return voisin3;
}
 
void OPT_TRIANGLE::change_voisin1(OPT_TRIANGLE* tri)
{
    voisin1=tri;
}
 
void OPT_TRIANGLE::change_voisin2(OPT_TRIANGLE* tri)
{
    voisin2=tri;
}
 
void OPT_TRIANGLE::change_voisin3(OPT_TRIANGLE* tri)
{
    voisin3=tri;
}