geometrie/src/fem_element3.cpp

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


#include "gestionversion.h"
#include <math.h>
#include "fem_element3.h"
#include "xfem_element3.h"
#include "fem_noeud.h"
#include "fem_maillage.h"
#include "mg_element_maillage.h"
#include "ot_boite_3D.h"

FEM_ELEMENT3::FEM_ELEMENT3(unsigned long num,class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(num,mai),etat(0)
{
}
FEM_ELEMENT3::FEM_ELEMENT3(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo):FEM_ELEMENT_MAILLAGE(num,topo),etat(0)
{
}
FEM_ELEMENT3::FEM_ELEMENT3(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(num,topo,mai),etat(0)
{
}

FEM_ELEMENT3::FEM_ELEMENT3(class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(mai),etat(0)
{
}

FEM_ELEMENT3::FEM_ELEMENT3(class MG_ELEMENT_TOPOLOGIQUE* topo):FEM_ELEMENT_MAILLAGE(topo),etat(0)
{
}
FEM_ELEMENT3::FEM_ELEMENT3(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(topo,mai),etat(0)
{
}


FEM_ELEMENT3::FEM_ELEMENT3(FEM_ELEMENT3& mdd):FEM_ELEMENT_MAILLAGE(mdd),etat(mdd.etat)
{
}


FEM_ELEMENT3::~FEM_ELEMENT3()
{
}


void FEM_ELEMENT3::extrapoler_solution_noeud(void)
{
    int nb=get_nb_fem_noeud();
    for (int i=0;i<nb;i++)
        for (int j=0;j<MAX_TYPE_SOLUTION;j++)  
          get_fem_noeud(i)->change_solution(solution[j],j);
}


int FEM_ELEMENT3::get_etat(void)
{
    return etat;
}

void FEM_ELEMENT3::change_etat(int num)
{
    etat=num;
}

void FEM_ELEMENT3::ajouter_xfem(XFEM_ELEMENT3* xele)
{
xfem.push_back(xele);
}

void FEM_ELEMENT3::supprimer_xfem(XFEM_ELEMENT3* xele)
{
std::vector<XFEM_ELEMENT3*>::iterator it=find(xfem.begin(),xfem.end(),xele);
XFEM_ELEMENT3* ele=(*it);
xfem.erase(it);
}

int FEM_ELEMENT3::get_nb_xfem(void)
{
return xfem.size();
}

XFEM_ELEMENT3* FEM_ELEMENT3::get_xfem(int num)
{
return xfem[num];
}


Revision:	378
Committed:	Tue Nov 27 20:54:58 2012 UTC (12 years, 5 months ago) by francois
File size:	2805 byte(s)
Log Message:	ajout de methode de construction pour les fem elements compatible avec le decoupage XFEM. + bug de format de fichier dans les fem elements non visible en dehors du xfem
#	User	Rev	Content
1	francois	283	//------------------------------------------------------------
2			//------------------------------------------------------------
3			// MAGiC
4			// Jean Christophe Cuilli�e et Vincent FRANCOIS
5			// D�artement de G�ie M�anique - UQTR
6			//------------------------------------------------------------
7			// Le projet MAGIC est un projet de recherche du d�artement
8			// de g�ie m�anique de l'Universit�du Qu�ec �
9			// Trois Rivi�es
10			// Les librairies ne peuvent �re utilis�s sans l'accord
11			// des auteurs (contact : francois@uqtr.ca)
12			//------------------------------------------------------------
13			//------------------------------------------------------------
14			//
15			// fem_tetra.cpp
16			//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20			// Version du 02/03/2006 �11H22
21			//------------------------------------------------------------
22			//------------------------------------------------------------
23
24
25			#include "gestionversion.h"
26			#include <math.h>
27	francois	309	#include "fem_element3.h"
28	francois	339	#include "xfem_element3.h"
29	francois	283	#include "fem_noeud.h"
30	francois	339	#include "fem_maillage.h"
31	francois	283	#include "mg_element_maillage.h"
32			#include "ot_boite_3D.h"
33
34	francois	309	FEM_ELEMENT3::FEM_ELEMENT3(unsigned long num,class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(num,mai),etat(0)
35	francois	283	{
36			}
37	francois	378	FEM_ELEMENT3::FEM_ELEMENT3(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo):FEM_ELEMENT_MAILLAGE(num,topo),etat(0)
38			{
39			}
40			FEM_ELEMENT3::FEM_ELEMENT3(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(num,topo,mai),etat(0)
41			{
42			}
43	francois	283
44	francois	309	FEM_ELEMENT3::FEM_ELEMENT3(class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(mai),etat(0)
45	francois	283	{
46			}
47
48	francois	378	FEM_ELEMENT3::FEM_ELEMENT3(class MG_ELEMENT_TOPOLOGIQUE* topo):FEM_ELEMENT_MAILLAGE(topo),etat(0)
49			{
50			}
51			FEM_ELEMENT3::FEM_ELEMENT3(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai):FEM_ELEMENT_MAILLAGE(topo,mai),etat(0)
52			{
53			}
54	francois	283
55	francois	378
56	francois	309	FEM_ELEMENT3::FEM_ELEMENT3(FEM_ELEMENT3& mdd):FEM_ELEMENT_MAILLAGE(mdd),etat(mdd.etat)
57	francois	283	{
58			}
59
60
61
62	francois	309	FEM_ELEMENT3::~FEM_ELEMENT3()
63	francois	283	{
64			}
65
66
67	francois	309	void FEM_ELEMENT3::extrapoler_solution_noeud(void)
68	francois	283	{
69			int nb=get_nb_fem_noeud();
70			for (int i=0;i<nb;i++)
71	francois	377	for (int j=0;j<MAX_TYPE_SOLUTION;j++)
72	francois	375	get_fem_noeud(i)->change_solution(solution[j],j);
73	francois	283	}
74
75
76
77	francois	309	int FEM_ELEMENT3::get_etat(void)
78	francois	283	{
79			return etat;
80			}
81
82	francois	309	void FEM_ELEMENT3::change_etat(int num)
83	francois	283	{
84			etat=num;
85			}
86
87	francois	339	void FEM_ELEMENT3::ajouter_xfem(XFEM_ELEMENT3* xele)
88			{
89			xfem.push_back(xele);
90			}
91	francois	283
92	francois	339	void FEM_ELEMENT3::supprimer_xfem(XFEM_ELEMENT3* xele)
93			{
94			std::vector<XFEM_ELEMENT3*>::iterator it=find(xfem.begin(),xfem.end(),xele);
95			XFEM_ELEMENT3* ele=(*it);
96			xfem.erase(it);
97			}
98
99			int FEM_ELEMENT3::get_nb_xfem(void)
100			{
101			return xfem.size();
102			}
103
104			XFEM_ELEMENT3* FEM_ELEMENT3::get_xfem(int num)
105			{
106			return xfem[num];
107			}
108
109