geometrie/src/fem_noeud.h

//------------------------------------------------------------
//------------------------------------------------------------
//  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_noeud.h
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 � 11H22
//------------------------------------------------------------
//------------------------------------------------------------
#ifndef _FEMNOEUD_
#define _FEMNOEUD_

#ifdef WINDOWS_VERSION
#ifdef BUILT_DLL_GEOMETRIE
#define DLLPORTGEOMETRIE __declspec(dllexport)
#else
#define DLLPORTGEOMETRIE __declspec(dllimport)
#endif
#else
#define DLLPORTGEOMETRIE
#endif

#include "fem_element_maillage.h"
#include "tpl_liste_entite.h"
#include "ot_boite_3d.h"


const unsigned int IDFEM_NOEUD=1000;


class DLLPORTGEOMETRIE FEM_NOEUD : public FEM_ELEMENT_MAILLAGE
{
public:
    FEM_NOEUD(unsigned long num,class MG_NOEUD* noeud);
    FEM_NOEUD(MG_NOEUD* noeud);
    FEM_NOEUD(class MG_ELEMENT_MAILLAGE* mai,double x,double y,double z,double xori=0.,double yori=0.,double zori=0.);
    FEM_NOEUD(unsigned long num,MG_ELEMENT_MAILLAGE* mai,double x,double y,double z,double xori=0.,double yori=0.,double zori=0.);
    FEM_NOEUD(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,double x,double y,double z,double xori=0.,double yori=0.,double zori=0.);
    FEM_NOEUD(class MG_ELEMENT_TOPOLOGIQUE* topo,double x,double y,double z,double xori=0.,double yori=0.,double zori=0.);
    FEM_NOEUD(FEM_NOEUD& mdd);
    virtual ~FEM_NOEUD();
    
    virtual void change_topologie_null(MG_ELEMENT_TOPOLOGIQUE* ele);

    virtual int get_type_entite();
    virtual BOITE_3D get_boite_3D(void);

    virtual double get_x(double coef=0.);
    virtual double get_y(double coef=0.);
    virtual double get_z(double coef=0.);
    virtual double get_dx(double coef=1.);
    virtual double get_dy(double coef=1.);
    virtual double get_dz(double coef=1.);
    virtual double get_rx(void);
    virtual double get_ry(void);
    virtual double get_rz(void);
    virtual double* get_coord(void);
    virtual double* get_coord_ori(void);

    virtual FEM_NOEUD* dupliquer(FEM_MAILLAGE *femmai,long decalage);


    virtual void change_x(double xx);
    virtual void change_y(double yy);
    virtual void change_z(double zz);
    virtual void change_dx(double xx);
    virtual void change_dy(double yy);
    virtual void change_dz(double zz);
    virtual void change_rx(double xx);
    virtual void change_ry(double yy);
    virtual void change_rz(double zz);
    virtual void change_coord(double *coo);
    virtual void reinit_boite_3D(void);
    virtual void reinit_boite_element_3D(void);

    int get_numero(void);
    void change_numero(int num);
    int get_numero_opt(void);
    void change_numero_opt(int num);

    virtual int get_dimension(void);

    virtual void enregistrer(std::ostream& o,double version);
    virtual void change_solution(double val,int num=0);
    virtual double get_solution(int num=0);
    virtual bool valide_parametre_EF(double *uvw);
   // virtual int get_etat(void);
   // virtual void change_etat(int num);
    virtual void get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin);

    TPL_LISTE_ENTITE<class FEM_ELEMENT1*> *get_lien_element1(void);
    TPL_LISTE_ENTITE<class FEM_ELEMENT1*> *get_lien_petit_element1(void);

    TPL_LISTE_ENTITE<class FEM_ELEMENT2*> *get_lien_element2(void);
    TPL_LISTE_ENTITE<class FEM_ELEMENT2*> *get_lien_petit_element2(void);

    TPL_LISTE_ENTITE<class FEM_ELEMENT3*> *get_lien_element3(void);
    TPL_LISTE_ENTITE<class FEM_ELEMENT3*> *get_lien_petit_element3(void);

protected:
    double xyz[3];
    double dxyz[3];
    double rxyz[3];
    double xyzorigine[3];
    int numero;
    int numero_opt;
    double solution[MAX_TYPE_SOLUTION];
    int activedeplacement1;
    int activedeplacement2;
    int activedeplacement3;
    int activereaction1;
    int activereaction2;
    int activereaction3;

   // int etat;
   

    TPL_LISTE_ENTITE<class FEM_ELEMENT1*> element1;
    TPL_LISTE_ENTITE<class FEM_ELEMENT1*> petit_element1;

    TPL_LISTE_ENTITE<class FEM_ELEMENT2*> element2;
    TPL_LISTE_ENTITE<class FEM_ELEMENT2*> petit_element2;

    TPL_LISTE_ENTITE<class FEM_ELEMENT3*> element3;
    TPL_LISTE_ENTITE<class FEM_ELEMENT3*> petit_element3;

};


#endif


Revision:	883
Committed:	Thu Apr 20 13:38:18 2017 UTC (8 years ago) by francois
Content type:	text/plain
File size:	4929 byte(s)
Log Message:	Creation d'un mailleur FEM pour pouvoir avoir des stratégies paramétrées de maillage. L'ancienne méthode construit disparait et est remplacée par l'utilisation d'un MAILLEUR_FEM. Stratégie de quadratisation mise en place : déplacer les noeuds pour s'assurer que les tetras quadratiques présentent une distortion au dessu d'une borne inférieure. Mais ces noeuds quittent la géométrie. Les paramètres dans ~/.magic.
#	User	Rev	Content
1	francois	283	//------------------------------------------------------------
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_noeud.h
16			//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20			// Version du 02/03/2006 � 11H22
21			//------------------------------------------------------------
22			//------------------------------------------------------------
23			#ifndef _FEMNOEUD_
24			#define _FEMNOEUD_
25
26			#ifdef WINDOWS_VERSION
27			#ifdef BUILT_DLL_GEOMETRIE
28			#define DLLPORTGEOMETRIE __declspec(dllexport)
29			#else
30			#define DLLPORTGEOMETRIE __declspec(dllimport)
31			#endif
32			#else
33			#define DLLPORTGEOMETRIE
34			#endif
35
36			#include "fem_element_maillage.h"
37			#include "tpl_liste_entite.h"
38	francois	481	#include "ot_boite_3d.h"
39	francois	283
40
41			const unsigned int IDFEM_NOEUD=1000;
42
43
44
45			class DLLPORTGEOMETRIE FEM_NOEUD : public FEM_ELEMENT_MAILLAGE
46			{
47			public:
48			FEM_NOEUD(unsigned long num,class MG_NOEUD* noeud);
49			FEM_NOEUD(MG_NOEUD* noeud);
50	francois	883	FEM_NOEUD(class MG_ELEMENT_MAILLAGE* mai,double x,double y,double z,double xori=0.,double yori=0.,double zori=0.);
51			FEM_NOEUD(unsigned long num,MG_ELEMENT_MAILLAGE* mai,double x,double y,double z,double xori=0.,double yori=0.,double zori=0.);
52			FEM_NOEUD(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,double x,double y,double z,double xori=0.,double yori=0.,double zori=0.);
53			FEM_NOEUD(class MG_ELEMENT_TOPOLOGIQUE* topo,double x,double y,double z,double xori=0.,double yori=0.,double zori=0.);
54	francois	283	FEM_NOEUD(FEM_NOEUD& mdd);
55			virtual ~FEM_NOEUD();
56	francois	414
57			virtual void change_topologie_null(MG_ELEMENT_TOPOLOGIQUE* ele);
58	francois	283
59			virtual int get_type_entite();
60			virtual BOITE_3D get_boite_3D(void);
61
62			virtual double get_x(double coef=0.);
63			virtual double get_y(double coef=0.);
64			virtual double get_z(double coef=0.);
65			virtual double get_dx(double coef=1.);
66			virtual double get_dy(double coef=1.);
67			virtual double get_dz(double coef=1.);
68	francois	325	virtual double get_rx(void);
69			virtual double get_ry(void);
70			virtual double get_rz(void);
71	francois	283	virtual double* get_coord(void);
72	francois	883	virtual double* get_coord_ori(void);
73	francois	283
74			virtual FEM_NOEUD* dupliquer(FEM_MAILLAGE *femmai,long decalage);
75
76
77			virtual void change_x(double xx);
78			virtual void change_y(double yy);
79			virtual void change_z(double zz);
80			virtual void change_dx(double xx);
81			virtual void change_dy(double yy);
82			virtual void change_dz(double zz);
83	francois	325	virtual void change_rx(double xx);
84			virtual void change_ry(double yy);
85			virtual void change_rz(double zz);
86	francois	283	virtual void change_coord(double *coo);
87	francois	684	virtual void reinit_boite_3D(void);
88			virtual void reinit_boite_element_3D(void);
89	francois	283
90			int get_numero(void);
91			void change_numero(int num);
92			int get_numero_opt(void);
93			void change_numero_opt(int num);
94
95			virtual int get_dimension(void);
96
97	francois	763	virtual void enregistrer(std::ostream& o,double version);
98	francois	375	virtual void change_solution(double val,int num=0);
99			virtual double get_solution(int num=0);
100	francois	635	virtual bool valide_parametre_EF(double *uvw);
101	francois	410	// virtual int get_etat(void);
102			// virtual void change_etat(int num);
103	francois	786	virtual void get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin);
104	francois	283
105	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT1> get_lien_element1(void);
106			TPL_LISTE_ENTITE<class FEM_ELEMENT1> get_lien_petit_element1(void);
107	francois	283
108	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT2> get_lien_element2(void);
109			TPL_LISTE_ENTITE<class FEM_ELEMENT2> get_lien_petit_element2(void);
110	francois	283
111	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT3> get_lien_element3(void);
112			TPL_LISTE_ENTITE<class FEM_ELEMENT3> get_lien_petit_element3(void);
113	francois	283
114			protected:
115			double xyz[3];
116			double dxyz[3];
117	francois	325	double rxyz[3];
118	francois	883	double xyzorigine[3];
119	francois	283	int numero;
120			int numero_opt;
121	francois	375	double solution[MAX_TYPE_SOLUTION];
122	francois	283	int activedeplacement1;
123			int activedeplacement2;
124			int activedeplacement3;
125	francois	339	int activereaction1;
126			int activereaction2;
127			int activereaction3;
128	francois	283
129	francois	410	// int etat;
130	francois	339
131	francois	283
132	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT1*> element1;
133			TPL_LISTE_ENTITE<class FEM_ELEMENT1*> petit_element1;
134	francois	283
135	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT2*> element2;
136			TPL_LISTE_ENTITE<class FEM_ELEMENT2*> petit_element2;
137	francois	283
138	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT3*> element3;
139			TPL_LISTE_ENTITE<class FEM_ELEMENT3*> petit_element3;
140	francois	283
141			};
142
143
144
145
146
147
148			#endif
149
150