fem_element1.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       fem_element1.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include <math.h>
#include "fem_element1.h"
#include "fem_noeud.h"
#include "mg_element_maillage.h"
#include "ot_mathematique.h"
 
FEM_ELEMENT1::FEM_ELEMENT1(unsigned long num,class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(num,mai)
{
}
 
FEM_ELEMENT1::FEM_ELEMENT1(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo):FEM_ELEMENT_MAILLAGE(num,topo)
{
}
 
 
FEM_ELEMENT1::FEM_ELEMENT1(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(num,topo,mai)
{
}
 
FEM_ELEMENT1::FEM_ELEMENT1(class MG_ELEMENT_TOPOLOGIQUE* topo):FEM_ELEMENT_MAILLAGE(topo)
{
}
 
 
FEM_ELEMENT1::FEM_ELEMENT1(class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(mai)
{
}
 
FEM_ELEMENT1::FEM_ELEMENT1(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(topo,mai)
{
}
 
FEM_ELEMENT1::FEM_ELEMENT1(FEM_ELEMENT1& mdd):FEM_ELEMENT_MAILLAGE(mdd)
{
}
 
 
 
FEM_ELEMENT1::~FEM_ELEMENT1()
{
}
 
 
void FEM_ELEMENT1::extrapoler_solution_noeud(void)
{
    int nb=get_nb_fem_noeud();
    for (int i=0;i<nb;i++)
      for (int j=0;j<MAX_TYPE_SOLUTION;j++)  
        get_fem_noeud(i)->change_solution(solution[j],j);
}
 
bool FEM_ELEMENT1::est_un_mini_element1(void)
{
return false;  
}
 
double FEM_ELEMENT1::get_jacobien(double* jac,double *uv,double unite)
{   
    int nb=get_nb_fem_noeud();
      
    OT_MATRICE_3D jacobien;
    jacobien(0,0)=0.;jacobien(0,1)=0.;jacobien(0,2)=0.;
      for (int k=0;k<nb;k++)
          {
      double valderiv=get_fonction_derive_interpolation(k+1,1,uv);
      double *xyz=get_fem_noeud(k)->get_coord();
          jacobien(0,0)=jacobien(0,0)+valderiv*xyz[0]*unite;
          jacobien(0,1)=jacobien(0,1)+valderiv*xyz[1]*unite;
          jacobien(0,2)=jacobien(0,2)+valderiv*xyz[2]*unite;
          }
        
    OT_VECTEUR_3D vec1,vec2,vec3;
    vec1.change_x(jacobien(0,0));
    vec1.change_y(jacobien(0,1));
    vec1.change_z(jacobien(0,2));
    
    if ((OPERATEUR::egal(vec1.get_x(),0.,1e-12)) && (OPERATEUR::egal(vec1.get_y(),0.,1e-12)))
        {
        vec2.change_x(0.);    
        vec2.change_y(vec1.get_z());    
        vec2.change_z(0.);    
        }
    else
        {
        vec2.change_x(vec1.get_y());    
        vec2.change_y(-vec1.get_x());    
        vec2.change_z(0.);    
        }
    vec2.norme();
    vec3=vec1&vec2;
    
    jacobien(1,0)=vec2.get_x();
    jacobien(1,1)=vec2.get_y();
    jacobien(1,2)=vec2.get_z();
    jacobien(2,0)=vec3.get_x();
    jacobien(2,1)=vec3.get_y();
    jacobien(2,2)=vec3.get_z();
    
    
    jac[0]=jacobien(0,0);
    jac[1]=jacobien(0,1);
    jac[2]=jacobien(0,2);
    
    jac[3]=jacobien(1,0);
    jac[4]=jacobien(1,1);
    jac[5]=jacobien(1,2);
    
    jac[6]=jacobien(2,0);
    jac[7]=jacobien(2,1);
    jac[8]=jacobien(2,2);
    
    double det=jacobien.get_determinant();
    return det;
}
 
void FEM_ELEMENT1::get_inverse_jacob(double* j,double *uv,double unite)
{
    double jac[9];
    get_jacobien(jac,uv,unite);
    OT_MATRICE_3D J;
    J(0,0)=jac[0];
    J(0,1)=jac[1];
    J(0,2)=jac[2];
    J(1,0)=jac[3];
    J(1,1)=jac[4];
    J(1,2)=jac[5];
    J(2,0)=jac[6];
    J(2,1)=jac[7];
    J(2,2)=jac[8];
    OT_MATRICE_3D j_i=J.inverse();
  
    
    j[0]=j_i(0,0);
    j[1]=j_i(0,1);
    j[2]=j_i(0,2);
    
    j[3]=j_i(1,0);
    j[4]=j_i(1,1);
    j[5]=j_i(1,2);
    
    j[6]=j_i(2,0);
    j[7]=j_i(2,1);
    j[8]=j_i(2,2);
}