geometrie/src/fem_noeud.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_noeud.cpp
//
//------------------------------------------------------------
//------------------------------------------------------------
//  COPYRIGHT 2000
//  Version du 02/03/2006 � 11H22
//------------------------------------------------------------
//------------------------------------------------------------


#include "gestionversion.h"
#include "fem_noeud.h"
#include "fem_maillage.h"
#include "mg_noeud.h"


FEM_NOEUD::FEM_NOEUD(MG_NOEUD* noeud):FEM_ELEMENT_MAILLAGE(noeud),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
{
    xyz[0]=noeud->get_x();
    xyz[1]=noeud->get_y();
    xyz[2]=noeud->get_z();
    xyzorigine[0]=noeud->get_x();
    xyzorigine[1]=noeud->get_y();
    xyzorigine[2]=noeud->get_z();
    dxyz[0]=0.;
    dxyz[1]=0.;
    dxyz[2]=0.;
    rxyz[0]=0.;
    rxyz[1]=0.;
    rxyz[2]=0.;
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
}

FEM_NOEUD::FEM_NOEUD(unsigned long num,MG_NOEUD* noeud):FEM_ELEMENT_MAILLAGE(num,noeud),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
{
    xyz[0]=noeud->get_x();
    xyz[1]=noeud->get_y();
    xyz[2]=noeud->get_z();
    xyzorigine[0]=noeud->get_x();
    xyzorigine[1]=noeud->get_y();
    xyzorigine[2]=noeud->get_z();
    dxyz[0]=0.;
    dxyz[1]=0.;
    dxyz[2]=0.;
    rxyz[0]=0.;
    rxyz[1]=0.;
    rxyz[2]=0.;
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
}

FEM_NOEUD::FEM_NOEUD(MG_ELEMENT_MAILLAGE* mai,double x,double y,double z,double xori,double yori,double zori):FEM_ELEMENT_MAILLAGE(mai),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
{
    xyz[0]=x;
    xyz[1]=y;
    xyz[2]=z;
    xyzorigine[0]=xori;
    xyzorigine[1]=yori;
    xyzorigine[2]=zori;
    dxyz[0]=0.;
    dxyz[1]=0.;
    dxyz[2]=0.;
    rxyz[0]=0.;
    rxyz[1]=0.;
    rxyz[2]=0.;
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
}

FEM_NOEUD* FEM_NOEUD::dupliquer(FEM_MAILLAGE *femmai,long decalage)
{
    FEM_NOEUD* noeud=new FEM_NOEUD(this->get_id()+decalage,liaison_topologique,xyz[0],xyz[1],xyz[2]);
    femmai->ajouter_fem_noeud(noeud);
    return noeud;
}
void FEM_NOEUD::reinit_boite_3D(void)
{
}


void FEM_NOEUD::reinit_boite_element_3D(void)
{
for (int i=0;i<get_lien_element3()->get_nb();i++)
      get_lien_element3()->get(i)->reinit_boite_3D();
    for (int i=0;i<get_lien_element2()->get_nb();i++)
      get_lien_element2()->get(i)->reinit_boite_3D();
    for (int i=0;i<get_lien_element1()->get_nb();i++)
      get_lien_element1()->get(i)->reinit_boite_3D();
}

FEM_NOEUD::FEM_NOEUD(unsigned long num,MG_ELEMENT_TOPOLOGIQUE* topo,double x,double y,double z,double xori,double yori,double zori):FEM_ELEMENT_MAILLAGE(num,topo),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
{
    xyz[0]=x;
    xyz[1]=y;
    xyz[2]=z;
    xyzorigine[0]=xori;
    xyzorigine[1]=yori;
    xyzorigine[2]=zori;
    dxyz[0]=0.;
    dxyz[1]=0.;
    dxyz[2]=0.;
    rxyz[0]=0.;
    rxyz[1]=0.;
    rxyz[2]=0.;
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
}

FEM_NOEUD::FEM_NOEUD(MG_ELEMENT_TOPOLOGIQUE* topo,double x,double y,double z,double xori,double yori,double zori):FEM_ELEMENT_MAILLAGE(topo),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
{
    xyz[0]=x;
    xyz[1]=y;
    xyz[2]=z;
    xyzorigine[0]=xori;
    xyzorigine[1]=yori;
    xyzorigine[2]=zori;
    dxyz[0]=0.;
    dxyz[1]=0.;
    dxyz[2]=0.;
    rxyz[0]=0.;
    rxyz[1]=0.;
    rxyz[2]=0.;
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
}

FEM_NOEUD::FEM_NOEUD(unsigned long num,MG_ELEMENT_MAILLAGE* mai,double x,double y,double z,double xori,double yori,double zori):FEM_ELEMENT_MAILLAGE(num,mai),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
{
    xyz[0]=x;
    xyz[1]=y;
    xyz[2]=z;
    xyzorigine[0]=xori;
    xyzorigine[1]=yori;
    xyzorigine[2]=zori;
    dxyz[0]=0.;
    dxyz[1]=0.;
    dxyz[2]=0.;
    rxyz[0]=0.;
    rxyz[1]=0.;
    rxyz[2]=0.;
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
}


FEM_NOEUD::FEM_NOEUD(FEM_NOEUD& mdd):FEM_ELEMENT_MAILLAGE(mdd),numero(-1),numero_opt(-1),activedeplacement1(mdd.activedeplacement1),activedeplacement2(mdd.activedeplacement2),activedeplacement3(mdd.activedeplacement3),activereaction1(mdd.activereaction1),activereaction2(mdd.activereaction2),activereaction3(mdd.activereaction3)
{
    xyz[0]=mdd.xyz[0];
    xyz[1]=mdd.xyz[1];
    xyz[2]=mdd.xyz[2];
    xyzorigine[0]=mdd.xyzorigine[0];
    xyzorigine[1]=mdd.xyzorigine[1];
    xyzorigine[2]=mdd.xyzorigine[2];
    dxyz[0]=mdd.dxyz[0];
    dxyz[1]=mdd.dxyz[1];
    dxyz[2]=mdd.dxyz[2];
    rxyz[0]=mdd.rxyz[0];
    rxyz[1]=mdd.rxyz[1];
    rxyz[2]=mdd.rxyz[2];
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
}

FEM_NOEUD::~FEM_NOEUD()
{
    if (liaison_topologique==NULL) return;
    if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
}
void FEM_NOEUD::get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin)
{
voisin.vide();
}

void FEM_NOEUD::change_topologie_null(MG_ELEMENT_TOPOLOGIQUE* ele)
{
if (liaison_topologique!=NULL) return;
liaison_topologique=ele;
if (liaison_topologique==NULL) return;
if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
}

int FEM_NOEUD::get_type_entite(void)
{
    return IDFEM_NOEUD;
}

double FEM_NOEUD::get_x(double coef)
{
    return xyz[0]+coef*dxyz[0];
}

double FEM_NOEUD::get_y(double coef)
{
    return xyz[1]+coef*dxyz[1];
}

double FEM_NOEUD::get_z(double coef)
{
    return xyz[2]+coef*dxyz[2];
}

double FEM_NOEUD::get_dx(double coef)
{
    return dxyz[0]*coef;
}

double FEM_NOEUD::get_dy(double coef)
{
    return dxyz[1]*coef;
}

double FEM_NOEUD::get_dz(double coef)
{
    return dxyz[2]*coef;
}


double FEM_NOEUD::get_rx(void)
{
return rxyz[0];
}
double FEM_NOEUD::get_ry(void)
{
return rxyz[1];
}

double FEM_NOEUD::get_rz(void)
{
return rxyz[2];
}


double* FEM_NOEUD::get_coord(void)
{
    return xyz;
}

double* FEM_NOEUD::get_coord_ori(void)
{
    return xyzorigine;
}


void FEM_NOEUD::change_x(double xx)
{
    xyz[0]=xx;
    reinit_boite_element_3D();
    
}

void FEM_NOEUD::change_y(double yy)
{
    xyz[1]=yy;    
    reinit_boite_element_3D();
}

void FEM_NOEUD::change_z(double zz)
{
    xyz[2]=zz;    
    reinit_boite_element_3D();
}
void FEM_NOEUD::change_rx(double xx)
{
    rxyz[0]=xx;
    activereaction1=1;
}

void FEM_NOEUD::change_ry(double yy)
{
    rxyz[1]=yy;
    activereaction2=1;
}

void FEM_NOEUD::change_rz(double zz)
{
    rxyz[2]=zz;
    activereaction3=1;

}

void FEM_NOEUD::change_dx(double xx)
{
    dxyz[0]=xx;
    activedeplacement1=1;
}

void FEM_NOEUD::change_dy(double yy)
{
    dxyz[1]=yy;
    activedeplacement2=1;
}

void FEM_NOEUD::change_dz(double zz)
{
    dxyz[2]=zz;
    activedeplacement3=1;
}

void FEM_NOEUD::change_coord(double *coo)
{
    xyz[0]=coo[0];
    xyz[1]=coo[1];
    xyz[2]=coo[2];    
    reinit_boite_element_3D();
}

int FEM_NOEUD::get_numero(void)
{
    return numero;
}

void FEM_NOEUD::change_numero(int num)
{
    numero=num;
}

int FEM_NOEUD::get_numero_opt(void)
{
    return numero_opt;
}

void FEM_NOEUD::change_numero_opt(int num)
{
    numero_opt=num;
}

bool FEM_NOEUD::valide_parametre_EF(double *uvw)
{
return true;  
}

int FEM_NOEUD::get_dimension(void)
{
    return 0;
}

void FEM_NOEUD::change_solution(double val,int num)
{
    solution[num]=val;
}


double FEM_NOEUD::get_solution(int num)
{
    return solution[num];
}

BOITE_3D FEM_NOEUD::get_boite_3D(void)
{
    BOITE_3D boite(xyz[0],xyz[1],xyz[2],xyz[0],xyz[1],xyz[2]);
    return boite;

}

void FEM_NOEUD::enregistrer(std::ostream& o,double version)
{
    if ((get_mg_element_maillage()!=NULL) && (get_lien_topologie()!=NULL)) o << "%" << get_id() << "=FEM_NOEUD($"<< get_lien_topologie()->get_id() << ",$" << maillage->get_id() << "," << xyz[0] << "," << xyz[1] << "," << xyz[2] << "," << numero  << ","  << numero_opt << ");" << std::endl;
    else if ((get_mg_element_maillage()!=NULL) && (get_lien_topologie()==NULL)) o << "%" << get_id() << "=FEM_NOEUD(NULL,$" << maillage->get_id() << "," << xyz[0] << "," << xyz[1] << "," << xyz[2] <<   "," << numero  << ","  << numero_opt << ");" << std::endl;
    else if ((get_mg_element_maillage()==NULL) && (get_lien_topologie()!=NULL)) o << "%" << get_id() << "=FEM_NOEUD($" << get_lien_topologie()->get_id()  << ",NULL," << xyz[0] << "," << xyz[1] << "," << xyz[2] << "," <<    numero  << ","  << numero_opt << ");" << std::endl;
    else  o << "%" << get_id() << "=FEM_NOEUD(NULL,NULL," << xyz[0] << "," << xyz[1] << "," << xyz[2]   << "," << numero  << ","  << numero_opt <<  ");" << std::endl;
    if (activereaction1==1)
        if (activereaction2==1)
            if (activereaction3==1)
          o << "%" << get_id() <<  "=FEM_NOEUD_REAC(" << rxyz[0] << "," << rxyz[1] << "," <<  rxyz[2] << ");" << std::endl;
    if (activedeplacement1==1)
        if (activedeplacement2==1)
            if (activedeplacement3==1)
                o << "%" << get_id() <<  "=FEM_NOEUD_DEF(" << dxyz[0] << "," << dxyz[1] << "," <<  dxyz[2] << ");" << std::endl;
}


TPL_LISTE_ENTITE<class FEM_ELEMENT1*>* FEM_NOEUD::get_lien_element1(void)
{
    return &element1;
}


TPL_LISTE_ENTITE<class FEM_ELEMENT1*>*  FEM_NOEUD::get_lien_petit_element1(void)
{
    return &petit_element1;
}


TPL_LISTE_ENTITE<class FEM_ELEMENT2*>* FEM_NOEUD::get_lien_element2(void)
{
    return &element2;
}


TPL_LISTE_ENTITE<class FEM_ELEMENT2*>* FEM_NOEUD::get_lien_petit_element2(void)
{
    return &petit_element2;
}

TPL_LISTE_ENTITE<class FEM_ELEMENT3*>* FEM_NOEUD::get_lien_element3(void)
{
    return &element3;
}


TPL_LISTE_ENTITE<class FEM_ELEMENT3*>* FEM_NOEUD::get_lien_petit_element3(void)
{
    return &petit_element3;
}
/*
int FEM_NOEUD::get_etat(void)
{
    return etat;
}
void FEM_NOEUD::change_etat(int num)
{
    etat=num;
}
*/
Revision:	883
Committed:	Thu Apr 20 13:38:18 2017 UTC (8 years ago) by francois
File size:	11750 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.cpp
16			//
17			//------------------------------------------------------------
18			//------------------------------------------------------------
19			// COPYRIGHT 2000
20			// Version du 02/03/2006 � 11H22
21			//------------------------------------------------------------
22			//------------------------------------------------------------
23
24
25			#include "gestionversion.h"
26			#include "fem_noeud.h"
27			#include "fem_maillage.h"
28			#include "mg_noeud.h"
29
30
31
32
33
34	francois	410	FEM_NOEUD::FEM_NOEUD(MG_NOEUD* noeud):FEM_ELEMENT_MAILLAGE(noeud),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
35	francois	283	{
36			xyz[0]=noeud->get_x();
37			xyz[1]=noeud->get_y();
38			xyz[2]=noeud->get_z();
39	francois	883	xyzorigine[0]=noeud->get_x();
40			xyzorigine[1]=noeud->get_y();
41			xyzorigine[2]=noeud->get_z();
42	francois	283	dxyz[0]=0.;
43			dxyz[1]=0.;
44			dxyz[2]=0.;
45	francois	325	rxyz[0]=0.;
46			rxyz[1]=0.;
47			rxyz[2]=0.;
48	francois	283	if (liaison_topologique==NULL) return;
49			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
50			}
51
52	francois	410	FEM_NOEUD::FEM_NOEUD(unsigned long num,MG_NOEUD* noeud):FEM_ELEMENT_MAILLAGE(num,noeud),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
53	francois	283	{
54			xyz[0]=noeud->get_x();
55			xyz[1]=noeud->get_y();
56			xyz[2]=noeud->get_z();
57	francois	883	xyzorigine[0]=noeud->get_x();
58			xyzorigine[1]=noeud->get_y();
59			xyzorigine[2]=noeud->get_z();
60	francois	283	dxyz[0]=0.;
61			dxyz[1]=0.;
62			dxyz[2]=0.;
63	francois	325	rxyz[0]=0.;
64			rxyz[1]=0.;
65			rxyz[2]=0.;
66	francois	283	if (liaison_topologique==NULL) return;
67			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
68			}
69
70	francois	883	FEM_NOEUD::FEM_NOEUD(MG_ELEMENT_MAILLAGE* mai,double x,double y,double z,double xori,double yori,double zori):FEM_ELEMENT_MAILLAGE(mai),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
71	francois	283	{
72			xyz[0]=x;
73			xyz[1]=y;
74			xyz[2]=z;
75	francois	883	xyzorigine[0]=xori;
76			xyzorigine[1]=yori;
77			xyzorigine[2]=zori;
78	francois	283	dxyz[0]=0.;
79			dxyz[1]=0.;
80			dxyz[2]=0.;
81	francois	325	rxyz[0]=0.;
82			rxyz[1]=0.;
83			rxyz[2]=0.;
84	francois	283	if (liaison_topologique==NULL) return;
85			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
86			}
87
88			FEM_NOEUD* FEM_NOEUD::dupliquer(FEM_MAILLAGE *femmai,long decalage)
89			{
90			FEM_NOEUD* noeud=new FEM_NOEUD(this->get_id()+decalage,liaison_topologique,xyz[0],xyz[1],xyz[2]);
91			femmai->ajouter_fem_noeud(noeud);
92			return noeud;
93			}
94	francois	684	void FEM_NOEUD::reinit_boite_3D(void)
95			{
96			}
97	francois	283
98
99	francois	684	void FEM_NOEUD::reinit_boite_element_3D(void)
100			{
101			for (int i=0;i<get_lien_element3()->get_nb();i++)
102			get_lien_element3()->get(i)->reinit_boite_3D();
103			for (int i=0;i<get_lien_element2()->get_nb();i++)
104			get_lien_element2()->get(i)->reinit_boite_3D();
105			for (int i=0;i<get_lien_element1()->get_nb();i++)
106			get_lien_element1()->get(i)->reinit_boite_3D();
107			}
108
109	francois	883	FEM_NOEUD::FEM_NOEUD(unsigned long num,MG_ELEMENT_TOPOLOGIQUE* topo,double x,double y,double z,double xori,double yori,double zori):FEM_ELEMENT_MAILLAGE(num,topo),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
110	francois	283	{
111			xyz[0]=x;
112			xyz[1]=y;
113			xyz[2]=z;
114	francois	883	xyzorigine[0]=xori;
115			xyzorigine[1]=yori;
116			xyzorigine[2]=zori;
117	francois	283	dxyz[0]=0.;
118			dxyz[1]=0.;
119			dxyz[2]=0.;
120	francois	325	rxyz[0]=0.;
121			rxyz[1]=0.;
122			rxyz[2]=0.;
123	francois	283	if (liaison_topologique==NULL) return;
124			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
125			}
126
127	francois	883	FEM_NOEUD::FEM_NOEUD(MG_ELEMENT_TOPOLOGIQUE* topo,double x,double y,double z,double xori,double yori,double zori):FEM_ELEMENT_MAILLAGE(topo),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
128	francois	378	{
129			xyz[0]=x;
130			xyz[1]=y;
131			xyz[2]=z;
132	francois	883	xyzorigine[0]=xori;
133			xyzorigine[1]=yori;
134			xyzorigine[2]=zori;
135	francois	378	dxyz[0]=0.;
136			dxyz[1]=0.;
137			dxyz[2]=0.;
138			rxyz[0]=0.;
139			rxyz[1]=0.;
140			rxyz[2]=0.;
141			if (liaison_topologique==NULL) return;
142			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
143			}
144
145	francois	883	FEM_NOEUD::FEM_NOEUD(unsigned long num,MG_ELEMENT_MAILLAGE* mai,double x,double y,double z,double xori,double yori,double zori):FEM_ELEMENT_MAILLAGE(num,mai),numero(-1),numero_opt(-1),activedeplacement1(0),activedeplacement2(0),activedeplacement3(0),activereaction1(0),activereaction2(0),activereaction3(0)
146	francois	283	{
147			xyz[0]=x;
148			xyz[1]=y;
149			xyz[2]=z;
150	francois	883	xyzorigine[0]=xori;
151			xyzorigine[1]=yori;
152			xyzorigine[2]=zori;
153	francois	283	dxyz[0]=0.;
154			dxyz[1]=0.;
155			dxyz[2]=0.;
156	francois	325	rxyz[0]=0.;
157			rxyz[1]=0.;
158			rxyz[2]=0.;
159	francois	283	if (liaison_topologique==NULL) return;
160			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
161			}
162
163
164
165	francois	410	FEM_NOEUD::FEM_NOEUD(FEM_NOEUD& mdd):FEM_ELEMENT_MAILLAGE(mdd),numero(-1),numero_opt(-1),activedeplacement1(mdd.activedeplacement1),activedeplacement2(mdd.activedeplacement2),activedeplacement3(mdd.activedeplacement3),activereaction1(mdd.activereaction1),activereaction2(mdd.activereaction2),activereaction3(mdd.activereaction3)
166	francois	283	{
167			xyz[0]=mdd.xyz[0];
168			xyz[1]=mdd.xyz[1];
169			xyz[2]=mdd.xyz[2];
170	francois	883	xyzorigine[0]=mdd.xyzorigine[0];
171			xyzorigine[1]=mdd.xyzorigine[1];
172			xyzorigine[2]=mdd.xyzorigine[2];
173	francois	283	dxyz[0]=mdd.dxyz[0];
174			dxyz[1]=mdd.dxyz[1];
175			dxyz[2]=mdd.dxyz[2];
176	francois	325	rxyz[0]=mdd.rxyz[0];
177			rxyz[1]=mdd.rxyz[1];
178			rxyz[2]=mdd.rxyz[2];
179	francois	283	if (liaison_topologique==NULL) return;
180			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
181			}
182
183			FEM_NOEUD::~FEM_NOEUD()
184			{
185			if (liaison_topologique==NULL) return;
186			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->supprimer(this);
187			}
188	francois	786	void FEM_NOEUD::get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin)
189			{
190			voisin.vide();
191			}
192	francois	283
193	francois	414	void FEM_NOEUD::change_topologie_null(MG_ELEMENT_TOPOLOGIQUE* ele)
194			{
195			if (liaison_topologique!=NULL) return;
196			liaison_topologique=ele;
197			if (liaison_topologique==NULL) return;
198			if (liaison_topologique->get_dimension()==0) liaison_topologique->get_lien_fem_maillage()->ajouter(this);
199			}
200	francois	283
201			int FEM_NOEUD::get_type_entite(void)
202			{
203			return IDFEM_NOEUD;
204			}
205
206			double FEM_NOEUD::get_x(double coef)
207			{
208			return xyz[0]+coef*dxyz[0];
209			}
210
211			double FEM_NOEUD::get_y(double coef)
212			{
213			return xyz[1]+coef*dxyz[1];
214			}
215
216			double FEM_NOEUD::get_z(double coef)
217			{
218			return xyz[2]+coef*dxyz[2];
219			}
220
221			double FEM_NOEUD::get_dx(double coef)
222			{
223			return dxyz[0]*coef;
224			}
225
226			double FEM_NOEUD::get_dy(double coef)
227			{
228			return dxyz[1]*coef;
229			}
230
231			double FEM_NOEUD::get_dz(double coef)
232			{
233			return dxyz[2]*coef;
234			}
235
236	francois	325
237
238			double FEM_NOEUD::get_rx(void)
239			{
240			return rxyz[0];
241			}
242			double FEM_NOEUD::get_ry(void)
243			{
244			return rxyz[1];
245			}
246
247			double FEM_NOEUD::get_rz(void)
248			{
249			return rxyz[2];
250			}
251
252
253
254
255	francois	283	double* FEM_NOEUD::get_coord(void)
256			{
257			return xyz;
258			}
259
260	francois	883	double* FEM_NOEUD::get_coord_ori(void)
261			{
262			return xyzorigine;
263			}
264	francois	283
265	francois	883
266	francois	283	void FEM_NOEUD::change_x(double xx)
267			{
268			xyz[0]=xx;
269	francois	684	reinit_boite_element_3D();
270
271	francois	283	}
272
273			void FEM_NOEUD::change_y(double yy)
274			{
275	francois	684	xyz[1]=yy;
276			reinit_boite_element_3D();
277	francois	283	}
278
279			void FEM_NOEUD::change_z(double zz)
280			{
281	francois	684	xyz[2]=zz;
282			reinit_boite_element_3D();
283	francois	283	}
284	francois	325	void FEM_NOEUD::change_rx(double xx)
285			{
286			rxyz[0]=xx;
287	francois	339	activereaction1=1;
288	francois	325	}
289	francois	283
290	francois	325	void FEM_NOEUD::change_ry(double yy)
291			{
292			rxyz[1]=yy;
293	francois	339	activereaction2=1;
294	francois	325	}
295	francois	283
296	francois	325	void FEM_NOEUD::change_rz(double zz)
297			{
298			rxyz[2]=zz;
299	francois	339	activereaction3=1;
300
301	francois	325	}
302
303	francois	283	void FEM_NOEUD::change_dx(double xx)
304			{
305			dxyz[0]=xx;
306			activedeplacement1=1;
307			}
308
309			void FEM_NOEUD::change_dy(double yy)
310			{
311			dxyz[1]=yy;
312			activedeplacement2=1;
313			}
314
315			void FEM_NOEUD::change_dz(double zz)
316			{
317			dxyz[2]=zz;
318			activedeplacement3=1;
319			}
320
321			void FEM_NOEUD::change_coord(double *coo)
322			{
323			xyz[0]=coo[0];
324			xyz[1]=coo[1];
325	francois	684	xyz[2]=coo[2];
326			reinit_boite_element_3D();
327	francois	283	}
328
329			int FEM_NOEUD::get_numero(void)
330			{
331			return numero;
332			}
333
334			void FEM_NOEUD::change_numero(int num)
335			{
336			numero=num;
337			}
338
339			int FEM_NOEUD::get_numero_opt(void)
340			{
341			return numero_opt;
342			}
343
344			void FEM_NOEUD::change_numero_opt(int num)
345			{
346			numero_opt=num;
347			}
348
349	francois	635	bool FEM_NOEUD::valide_parametre_EF(double *uvw)
350			{
351			return true;
352			}
353	francois	283
354			int FEM_NOEUD::get_dimension(void)
355			{
356			return 0;
357			}
358
359	francois	375	void FEM_NOEUD::change_solution(double val,int num)
360	francois	283	{
361	francois	375	solution[num]=val;
362	francois	283	}
363
364
365	francois	375	double FEM_NOEUD::get_solution(int num)
366	francois	283	{
367	francois	375	return solution[num];
368	francois	283	}
369
370			BOITE_3D FEM_NOEUD::get_boite_3D(void)
371			{
372			BOITE_3D boite(xyz[0],xyz[1],xyz[2],xyz[0],xyz[1],xyz[2]);
373			return boite;
374
375			}
376
377	francois	763	void FEM_NOEUD::enregistrer(std::ostream& o,double version)
378	francois	283	{
379			if ((get_mg_element_maillage()!=NULL) && (get_lien_topologie()!=NULL)) o << "%" << get_id() << "=FEM_NOEUD($"<< get_lien_topologie()->get_id() << ",$" << maillage->get_id() << "," << xyz[0] << "," << xyz[1] << "," << xyz[2] << "," << numero << "," << numero_opt << ");" << std::endl;
380			else if ((get_mg_element_maillage()!=NULL) && (get_lien_topologie()==NULL)) o << "%" << get_id() << "=FEM_NOEUD(NULL,$" << maillage->get_id() << "," << xyz[0] << "," << xyz[1] << "," << xyz[2] << "," << numero << "," << numero_opt << ");" << std::endl;
381			else if ((get_mg_element_maillage()==NULL) && (get_lien_topologie()!=NULL)) o << "%" << get_id() << "=FEM_NOEUD($" << get_lien_topologie()->get_id() << ",NULL," << xyz[0] << "," << xyz[1] << "," << xyz[2] << "," << numero << "," << numero_opt << ");" << std::endl;
382			else o << "%" << get_id() << "=FEM_NOEUD(NULL,NULL," << xyz[0] << "," << xyz[1] << "," << xyz[2] << "," << numero << "," << numero_opt << ");" << std::endl;
383	francois	339	if (activereaction1==1)
384			if (activereaction2==1)
385			if (activereaction3==1)
386			o << "%" << get_id() << "=FEM_NOEUD_REAC(" << rxyz[0] << "," << rxyz[1] << "," << rxyz[2] << ");" << std::endl;
387	francois	283	if (activedeplacement1==1)
388			if (activedeplacement2==1)
389			if (activedeplacement3==1)
390			o << "%" << get_id() << "=FEM_NOEUD_DEF(" << dxyz[0] << "," << dxyz[1] << "," << dxyz[2] << ");" << std::endl;
391			}
392
393
394
395	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT1> FEM_NOEUD::get_lien_element1(void)
396	francois	283	{
397	francois	309	return &element1;
398	francois	283	}
399
400
401	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT1> FEM_NOEUD::get_lien_petit_element1(void)
402	francois	283	{
403	francois	309	return &petit_element1;
404	francois	283	}
405
406
407	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT2> FEM_NOEUD::get_lien_element2(void)
408	francois	283	{
409	francois	309	return &element2;
410	francois	283	}
411
412
413	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT2> FEM_NOEUD::get_lien_petit_element2(void)
414	francois	283	{
415	francois	309	return &petit_element2;
416	francois	283	}
417
418	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT3> FEM_NOEUD::get_lien_element3(void)
419	francois	283	{
420	francois	309	return &element3;
421	francois	283	}
422
423
424	francois	309	TPL_LISTE_ENTITE<class FEM_ELEMENT3> FEM_NOEUD::get_lien_petit_element3(void)
425	francois	283	{
426	francois	309	return &petit_element3;
427	francois	283	}
428	francois	410	/*
429	francois	283	int FEM_NOEUD::get_etat(void)
430			{
431			return etat;
432			}
433			void FEM_NOEUD::change_etat(int num)
434			{
435			etat=num;
436			}
437	francois	410	*/