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)
{
    o << "%" << get_id() << "=XFEM_SEGMENT2($"<< femmai->get_id() << ",$" << topo->get_id() << ",$"  << 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:	399
Committed:	Mon May 27 15:27:41 2013 UTC (11 years, 11 months ago) by francois
File size:	3505 byte(s)
Log Message:	Creation d'une structure XFEM et FEM complete avec des element0 et des xelement0
#	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	399	o << "%" << get_id() << "=XFEM_SEGMENT2($"<< femmai->get_id() << ",$" << topo->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << "," << etat << ");" << std::endl;
89	francois	380
90			}
91
92
93
94
95			int XFEM_SEGMENT2::nb_fonction_interpolation(void)
96			{
97			return 2;
98			}
99			double XFEM_SEGMENT2::get_fonction_interpolation(int num,double *uv)
100			{
101			double val;
102			switch (num)
103			{
104			case 1:
105			val=0.5*(1-uv[0]);
106			break;
107			case 2:
108			val=0.5*(1+uv[0]);
109			break;
110			}
111			return val;
112			}
113			double XFEM_SEGMENT2::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
114			{
115			double val;
116
117			switch (num)
118			{
119			case 1:
120			switch (num_variable)
121			{
122			case 1:
123			val=-0.5;
124			break;
125			} break;
126			case 2:
127			switch (num_variable)
128			{
129			case 1:
130			val=0.5;
131			break;
132			}break;
133
134			}
135			return val;
136			}
137			double XFEM_SEGMENT2::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
138			{
139			}
140			void XFEM_SEGMENT2::get_inverse_jacob(double* j,double *uv,double unite)
141			{
142			}