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, 4 months 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.
#	Content
1	//------------------------------------------------------------
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	#include "ot_boite_3d.h"
39
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	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	FEM_NOEUD(FEM_NOEUD& mdd);
55	virtual ~FEM_NOEUD();
56
57	virtual void change_topologie_null(MG_ELEMENT_TOPOLOGIQUE* ele);
58
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	virtual double get_rx(void);
69	virtual double get_ry(void);
70	virtual double get_rz(void);
71	virtual double* get_coord(void);
72	virtual double* get_coord_ori(void);
73
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	virtual void change_rx(double xx);
84	virtual void change_ry(double yy);
85	virtual void change_rz(double zz);
86	virtual void change_coord(double *coo);
87	virtual void reinit_boite_3D(void);
88	virtual void reinit_boite_element_3D(void);
89
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	virtual void enregistrer(std::ostream& o,double version);
98	virtual void change_solution(double val,int num=0);
99	virtual double get_solution(int num=0);
100	virtual bool valide_parametre_EF(double *uvw);
101	// virtual int get_etat(void);
102	// virtual void change_etat(int num);
103	virtual void get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin);
104
105	TPL_LISTE_ENTITE<class FEM_ELEMENT1> get_lien_element1(void);
106	TPL_LISTE_ENTITE<class FEM_ELEMENT1> get_lien_petit_element1(void);
107
108	TPL_LISTE_ENTITE<class FEM_ELEMENT2> get_lien_element2(void);
109	TPL_LISTE_ENTITE<class FEM_ELEMENT2> get_lien_petit_element2(void);
110
111	TPL_LISTE_ENTITE<class FEM_ELEMENT3> get_lien_element3(void);
112	TPL_LISTE_ENTITE<class FEM_ELEMENT3> get_lien_petit_element3(void);
113
114	protected:
115	double xyz[3];
116	double dxyz[3];
117	double rxyz[3];
118	double xyzorigine[3];
119	int numero;
120	int numero_opt;
121	double solution[MAX_TYPE_SOLUTION];
122	int activedeplacement1;
123	int activedeplacement2;
124	int activedeplacement3;
125	int activereaction1;
126	int activereaction2;
127	int activereaction3;
128
129	// int etat;
130
131
132	TPL_LISTE_ENTITE<class FEM_ELEMENT1*> element1;
133	TPL_LISTE_ENTITE<class FEM_ELEMENT1*> petit_element1;
134
135	TPL_LISTE_ENTITE<class FEM_ELEMENT2*> element2;
136	TPL_LISTE_ENTITE<class FEM_ELEMENT2*> petit_element2;
137
138	TPL_LISTE_ENTITE<class FEM_ELEMENT3*> element3;
139	TPL_LISTE_ENTITE<class FEM_ELEMENT3*> petit_element3;
140
141	};
142
143
144
145
146
147
148	#endif
149
150