ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/geometrie/src/fem_element0.cpp
Revision: 786
Committed: Wed Mar 9 21:10:49 2016 UTC (9 years, 2 months ago) by francois
File size: 4220 byte(s)
Log Message: 
Possibilité d'obtenir la matrice de rigidite de code aster sous form de matrice bande.
Pour cela ajouter -opt dans l'operation -fem à la création du maillage FEM et dans -calculaster utiliser les nouveaux types d'étude rigidite.
Attention pour cela il faut une version de code aster modifiée
De manière automatique elle peut etre optenu en faisant ./instal_aster 12.5mod
Uniquement implantée avec la version 12.5 de code aster.

File Contents

# User Rev Content
1 francois 399 //------------------------------------------------------------
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_segment.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     #include "fem_element0.h"
28     #include "fem_noeud.h"
29     #include "fem_maillage.h"
30     #include "mg_element_maillage.h"
31    
32     FEM_ELEMENT0::FEM_ELEMENT0(unsigned long num,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT_MAILLAGE(num,mai),FEM_TEMPLATE_ELEMENT<1>(tabnoeud)
33     {
34     }
35    
36     FEM_ELEMENT0::FEM_ELEMENT0(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT_MAILLAGE(num,topo),FEM_TEMPLATE_ELEMENT<1>(tabnoeud)
37     {
38     }
39    
40    
41     FEM_ELEMENT0::FEM_ELEMENT0(unsigned long num,class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT_MAILLAGE(num,topo,mai),FEM_TEMPLATE_ELEMENT<1>(tabnoeud)
42     {
43     }
44    
45     FEM_ELEMENT0::FEM_ELEMENT0(class MG_ELEMENT_TOPOLOGIQUE* topo,FEM_NOEUD** tabnoeud):FEM_ELEMENT_MAILLAGE(topo),FEM_TEMPLATE_ELEMENT<1>(tabnoeud)
46     {
47     }
48    
49    
50     FEM_ELEMENT0::FEM_ELEMENT0(class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT_MAILLAGE(mai),FEM_TEMPLATE_ELEMENT<1>(tabnoeud)
51     {
52     }
53    
54     FEM_ELEMENT0::FEM_ELEMENT0(class MG_ELEMENT_TOPOLOGIQUE* topo,class MG_ELEMENT_MAILLAGE* mai,FEM_NOEUD** tabnoeud):FEM_ELEMENT_MAILLAGE(topo,mai),FEM_TEMPLATE_ELEMENT<1>(tabnoeud)
55     {
56     }
57    
58     FEM_ELEMENT0::FEM_ELEMENT0(FEM_ELEMENT0& mdd):FEM_ELEMENT_MAILLAGE(mdd),FEM_TEMPLATE_ELEMENT<1>(mdd)
59     {
60     }
61    
62 francois 786 void FEM_ELEMENT0::get_voisin_noeud(class FEM_NOEUD* no,TPL_LISTE_ENTITE<FEM_NOEUD*> &voisin)
63     {
64     voisin.vide();
65     }
66 francois 399
67    
68     FEM_ELEMENT0::~FEM_ELEMENT0()
69     {
70     }
71    
72     FEM_ELEMENT_MAILLAGE* FEM_ELEMENT0::dupliquer(FEM_MAILLAGE *femmai,long decalage)
73     {
74     FEM_NOEUD* tabnoeud[1];
75     tabnoeud[0]=femmai->get_fem_noeudid(tab[0]->get_id()+decalage);
76     FEM_ELEMENT0* fem=new FEM_ELEMENT0(get_id()+decalage,maillage,tabnoeud);
77     femmai->ajouter_fem_element0(fem);
78     return fem;
79     }
80    
81     int FEM_ELEMENT0::get_type_entite(void)
82     {
83     return IDFEM_ELEMENT0;
84     }
85    
86     int FEM_ELEMENT0::get_dimension(void)
87     {
88     return 0;
89     }
90    
91 francois 684 void FEM_ELEMENT0::reinit_boite_3D(void)
92     {
93     FEM_TEMPLATE_ELEMENT<1>::reinit_boite_3D();
94     }
95 francois 399
96     int FEM_ELEMENT0::get_nb_fem_noeud(void)
97     {
98     return FEM_TEMPLATE_ELEMENT<1>::get_nb_fem_noeud();
99     }
100    
101     FEM_NOEUD* FEM_ELEMENT0::get_fem_noeud(int num)
102     {
103     return FEM_TEMPLATE_ELEMENT<1>::get_fem_noeud(num);
104     }
105    
106     void FEM_ELEMENT0::change_noeud(int num,FEM_NOEUD* noeud)
107     {
108     FEM_TEMPLATE_ELEMENT<1>::change_noeud(num,noeud);
109     }
110    
111     BOITE_3D& FEM_ELEMENT0::get_boite_3D(void)
112     {
113     return FEM_TEMPLATE_ELEMENT<1>::get_boite_3D();
114     }
115    
116 francois 635 bool FEM_ELEMENT0::valide_parametre_EF(double *uvw)
117     {
118     return true;
119     }
120 francois 399
121 francois 763 void FEM_ELEMENT0::enregistrer(std::ostream& o,double version)
122 francois 399 {
123     if (maillage!=NULL)
124     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_ELEMENT_NOEUD($"<< get_lien_topologie()->get_id() << ",$"<< maillage->get_id() << ",$" << tab[0]->get_id() << ");" << std::endl;
125     else o << "%" << get_id() << "=FEM_ELEMENT_NOEUD(NULL,$" << maillage->get_id() <<",$" << tab[0]->get_id() << ");" << std::endl;
126     if (maillage==NULL)
127     if (get_lien_topologie()!=NULL) o << "%" << get_id() << "=FEM_SEGMENT2($"<< get_lien_topologie()->get_id() << ",NULL,$" << tab[0]->get_id() << ");" << std::endl;
128     else o << "%" << get_id() << "=FEM_ELEMENT_NOEUD(NULL,NULL,$" << tab[0]->get_id() << ");" << std::endl;
129     }
130    
131    
132    
133