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