geometrie/src/xfem_segment2.cpp

//------------------------------------------------------------
//------------------------------------------------------------
//  MAGiC
//  Jean Christophe Cuilli�re et Vincent FRANCOIS
//  D�partement de G�nie M�canique - UQTR
//------------------------------------------------------------
//  Le projet MAGIC est un projet de recherche du d�partement
//  de g�nie m�canique de l'Universit� du Qu�bec �
//  Trois Rivi�res
//  Les librairies ne peuvent �tre utilis�es sans l'accord
//  des auteurs (contact : francois@uqtr.ca)
//------------------------------------------------------------
//------------------------------------------------------------
//
//       fem_element34.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 � 11H22
//------------------------------------------------------------
//------------------------------------------------------------


#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)
{
   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)
{
}
void XFEM_SEGMENT2::get_inverse_jacob(double* j,double *uv,double unite)
{
}
Revision:	404
Committed:	Thu May 30 13:05:01 2013 UTC (11 years, 11 months ago) by francois
File size:	3706 byte(s)
Log Message:	Enregistrement des entites XFEM
#	User	Rev	Content
1	francois	380	//------------------------------------------------------------
2			//------------------------------------------------------------
3			// MAGiC
4			// Jean Christophe Cuilli�re et Vincent FRANCOIS
5			// D�partement de G�nie M�canique - UQTR
6			//------------------------------------------------------------
7			// Le projet MAGIC est un projet de recherche du d�partement
8			// de g�nie m�canique de l'Universit� du Qu�bec �
9			// Trois Rivi�res
10			// Les librairies ne peuvent �tre utilis�es sans l'accord
11			// des auteurs (contact : francois@uqtr.ca)
12			//------------------------------------------------------------
13			//------------------------------------------------------------
14			//
15			// fem_element34.cpp
16			//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20			// Version du 02/03/2006 � 11H22
21			//------------------------------------------------------------
22			//------------------------------------------------------------
23
24
25			#include "gestionversion.h"
26			#include "xfem_segment2.h"
27			#include "fem_segment2.h"
28			#include "fem_noeud.h"
29			#include "mg_element_maillage.h"
30			#include "fem_maillage.h"
31
32
33	francois	399	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)
34	francois	380	{
35	francois	399	femmai->ajouter_xfem(this);
36	francois	380	}
37
38	francois	399	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)
39	francois	380	{
40	francois	399	femmai->ajouter_xfem(this);
41	francois	380	}
42
43	francois	399	XFEM_SEGMENT2::XFEM_SEGMENT2(XFEM_SEGMENT2& mdd):XFEM_ELEMENT1(mdd),FEM_TEMPLATE_ELEMENT<2>(mdd)
44	francois	380	{
45			}
46
47			XFEM_SEGMENT2::~XFEM_SEGMENT2()
48			{
49	francois	399	femmai->supprimer_xfem(this);
50	francois	380	}
51
52
53
54			int XFEM_SEGMENT2::get_type_entite(void)
55			{
56			return IDXFEM_SEGMENT2;
57			}
58
59			int XFEM_SEGMENT2::get_dimension(void)
60			{
61			return 1;
62			}
63
64
65			int XFEM_SEGMENT2::get_nb_fem_noeud(void)
66			{
67			return FEM_TEMPLATE_ELEMENT<2>::get_nb_fem_noeud();
68			}
69
70			FEM_NOEUD* XFEM_SEGMENT2::get_fem_noeud(int num)
71			{
72			return FEM_TEMPLATE_ELEMENT<2>::get_fem_noeud(num);
73			}
74
75			void XFEM_SEGMENT2::change_noeud(int num,FEM_NOEUD* noeud)
76			{
77			FEM_TEMPLATE_ELEMENT<2>::change_noeud(num,noeud);
78			}
79
80
81			BOITE_3D& XFEM_SEGMENT2::get_boite_3D(void)
82			{
83			return FEM_TEMPLATE_ELEMENT<2>::get_boite_3D();
84			}
85
86			void XFEM_SEGMENT2::enregistrer(std::ostream& o)
87			{
88	francois	404	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;
89			else o << "%" << get_id() << "=XFEM_SEGMENT2($"<< femmai->get_id() << ",NULL,$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << "," << etat << ");" << std::endl;
90	francois	380
91			}
92
93
94
95
96			int XFEM_SEGMENT2::nb_fonction_interpolation(void)
97			{
98			return 2;
99			}
100			double XFEM_SEGMENT2::get_fonction_interpolation(int num,double *uv)
101			{
102			double val;
103			switch (num)
104			{
105			case 1:
106			val=0.5*(1-uv[0]);
107			break;
108			case 2:
109			val=0.5*(1+uv[0]);
110			break;
111			}
112			return val;
113			}
114			double XFEM_SEGMENT2::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
115			{
116			double val;
117
118			switch (num)
119			{
120			case 1:
121			switch (num_variable)
122			{
123			case 1:
124			val=-0.5;
125			break;
126			} break;
127			case 2:
128			switch (num_variable)
129			{
130			case 1:
131			val=0.5;
132			break;
133			}break;
134
135			}
136			return val;
137			}
138			double XFEM_SEGMENT2::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
139			{
140			}
141			void XFEM_SEGMENT2::get_inverse_jacob(double* j,double *uv,double unite)
142			{
143			}