xfem_triangle3.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/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       xfem_triangle3.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:58:53 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
 
 
#include "gestionversion.h"
#include "xfem_triangle3.h"
#include "fem_tetra4.h"
#include "fem_noeud.h"
#include "mg_element_maillage.h"
#include "fem_maillage.h"
 
 
XFEM_TRIANGLE3::XFEM_TRIANGLE3(unsigned long num,class FEM_ELEMENT_MAILLAGE *mai,class MG_ELEMENT_TOPOLOGIQUE *mgtopo,class FEM_NOEUD** tabnoeud):XFEM_ELEMENT2(num,mai,mgtopo),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
{
femmai->ajouter_xfem(this);
}
 
XFEM_TRIANGLE3::XFEM_TRIANGLE3(class FEM_ELEMENT_MAILLAGE *mai,class MG_ELEMENT_TOPOLOGIQUE *mgtopo,FEM_NOEUD** tabnoeud):XFEM_ELEMENT2(mai,mgtopo),FEM_TEMPLATE_ELEMENT<3>(tabnoeud)
{
femmai->ajouter_xfem(this);
}
 
XFEM_TRIANGLE3::XFEM_TRIANGLE3(XFEM_TRIANGLE3& mdd):XFEM_ELEMENT2(mdd),FEM_TEMPLATE_ELEMENT<3>(mdd)
{
}
 
XFEM_TRIANGLE3::~XFEM_TRIANGLE3()
{
femmai->supprimer_xfem(this);
}
 
 
 
int XFEM_TRIANGLE3::get_type_entite(void)
{
    return   IDXFEM_TRIANGLE3;
}
 
int XFEM_TRIANGLE3::get_dimension(void)
{
    return 2;
}
 
 
int XFEM_TRIANGLE3::get_nb_fem_noeud(void)
{
    return FEM_TEMPLATE_ELEMENT<3>::get_nb_fem_noeud();
}
 
FEM_NOEUD* XFEM_TRIANGLE3::get_fem_noeud(int num)
{
    return FEM_TEMPLATE_ELEMENT<3>::get_fem_noeud(num);
}
 
void XFEM_TRIANGLE3::change_noeud(int num,FEM_NOEUD* noeud)
{
    FEM_TEMPLATE_ELEMENT<3>::change_noeud(num,noeud);
}
 
 
BOITE_3D& XFEM_TRIANGLE3::get_boite_3D(void)
{
    return FEM_TEMPLATE_ELEMENT<3>::get_boite_3D();
}
 
void XFEM_TRIANGLE3::enregistrer(std::ostream& o,double version)
{
   if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=XFEM_TRIANGLE3($"<< femmai->get_id() << ",$" << topo->get_id() << ",$"  << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << "," << etat << ");" << std::endl;
    else o << "%" << get_id() << "=XFEM_TRIANGLE3($"<< femmai->get_id() << ",NULL,$"  << tab[0]->get_id() << ",$" << tab[1]->get_id() << ",$" << tab[2]->get_id() << "," << etat << ");" << std::endl;
 
}
 
 
 
 
int XFEM_TRIANGLE3::nb_fonction_interpolation(void)
{
return 3;
}
double XFEM_TRIANGLE3::get_fonction_interpolation(int num,double *uv)
{
double val;
    switch (num)
    {
    case 1:
        val=1-uv[0]-uv[1];
        break;
    case 2:
        val=uv[0];
        break;
    case 3:
        val=uv[1];
        break;
    }
    return val;
}
double XFEM_TRIANGLE3::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
{
double val;
 
    switch (num)
    {
    case 1:
        switch (num_variable)
        {
        case 1:
            val=-1;
            break;
        case 2:
            val=-1;
            break;
        } break;
    case 2:
        switch (num_variable)
        {
        case 1:
            val=1;
            break;
        case 2:
            val=0;
            break;
        }break;
    case 3:
        switch (num_variable)
        {
        case 1:
            val=0;
            break;
        case 2:
            val=1;
            break;
        }break;
    
    }
    return val;
}
double XFEM_TRIANGLE3::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
{
return 0.;
}
 
void XFEM_TRIANGLE3::get_inverse_jacob(double* j,double *uv,double unite)
{
}