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_SEGMENT2* seg,class FEM_NOEUD** tabnoeud):XFEM_ELEMENT1(num,seg->get_mg_element_maillage()),FEM_TEMPLATE_ELEMENT<2>(tabnoeud),segment(seg)
{
segment->ajouter_xfem(this);
}

XFEM_SEGMENT2::XFEM_SEGMENT2(class FEM_SEGMENT2* seg,FEM_NOEUD** tabnoeud):XFEM_ELEMENT1(seg->get_mg_element_maillage()),FEM_TEMPLATE_ELEMENT<2>(tabnoeud),segment(seg)
{
segment->ajouter_xfem(this);
}

XFEM_SEGMENT2::XFEM_SEGMENT2(XFEM_SEGMENT2& mdd):XFEM_ELEMENT1(mdd),FEM_TEMPLATE_ELEMENT<2>(mdd),segment(mdd.segment),n1(mdd.n1),n2(mdd.n2)
{
}

XFEM_SEGMENT2::~XFEM_SEGMENT2()
{
segment->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($"<< segment->get_id() << ",$"  << tab[0]->get_id() << ",$" << tab[1]->get_id() << "," << etat << ");" << std::endl;

}

FEM_ELEMENT_MAILLAGE* XFEM_SEGMENT2::dupliquer(FEM_MAILLAGE *femmai,long decalage)
{
tab[0]=femmai->get_fem_noeudid(n1);
tab[1]=femmai->get_fem_noeudid(n2);
}


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:	380
Committed:	Thu Nov 29 15:40:22 2012 UTC (12 years, 5 months ago) by francois
File size:	3661 byte(s)
Log Message:	gestion des entités xfem 1d et 2d
#	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			XFEM_SEGMENT2::XFEM_SEGMENT2(unsigned long num,class FEM_SEGMENT2* seg,class FEM_NOEUD** tabnoeud):XFEM_ELEMENT1(num,seg->get_mg_element_maillage()),FEM_TEMPLATE_ELEMENT<2>(tabnoeud),segment(seg)
34			{
35			segment->ajouter_xfem(this);
36			}
37
38			XFEM_SEGMENT2::XFEM_SEGMENT2(class FEM_SEGMENT2* seg,FEM_NOEUD** tabnoeud):XFEM_ELEMENT1(seg->get_mg_element_maillage()),FEM_TEMPLATE_ELEMENT<2>(tabnoeud),segment(seg)
39			{
40			segment->ajouter_xfem(this);
41			}
42
43			XFEM_SEGMENT2::XFEM_SEGMENT2(XFEM_SEGMENT2& mdd):XFEM_ELEMENT1(mdd),FEM_TEMPLATE_ELEMENT<2>(mdd),segment(mdd.segment),n1(mdd.n1),n2(mdd.n2)
44			{
45			}
46
47			XFEM_SEGMENT2::~XFEM_SEGMENT2()
48			{
49			segment->supprimer_xfem(this);
50			}
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			o << "%" << get_id() << "=XFEM_SEGMENT2($"<< segment->get_id() << ",$" << tab[0]->get_id() << ",$" << tab[1]->get_id() << "," << etat << ");" << std::endl;
89
90			}
91
92			FEM_ELEMENT_MAILLAGE* XFEM_SEGMENT2::dupliquer(FEM_MAILLAGE *femmai,long decalage)
93			{
94			tab[0]=femmai->get_fem_noeudid(n1);
95			tab[1]=femmai->get_fem_noeudid(n2);
96			}
97
98
99
100			int XFEM_SEGMENT2::nb_fonction_interpolation(void)
101			{
102			return 2;
103			}
104			double XFEM_SEGMENT2::get_fonction_interpolation(int num,double *uv)
105			{
106			double val;
107			switch (num)
108			{
109			case 1:
110			val=0.5*(1-uv[0]);
111			break;
112			case 2:
113			val=0.5*(1+uv[0]);
114			break;
115			}
116			return val;
117			}
118			double XFEM_SEGMENT2::get_fonction_derive_interpolation(int num,int num_variable,double *uv)
119			{
120			double val;
121
122			switch (num)
123			{
124			case 1:
125			switch (num_variable)
126			{
127			case 1:
128			val=-0.5;
129			break;
130			} break;
131			case 2:
132			switch (num_variable)
133			{
134			case 1:
135			val=0.5;
136			break;
137			}break;
138
139			}
140			return val;
141			}
142			double XFEM_SEGMENT2::get_jacobien(double* jac,double *uv,int& li,int& col,double unite)
143			{
144			}
145			void XFEM_SEGMENT2::get_inverse_jacob(double* j,double *uv,double unite)
146			{
147			}