ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/mecanique/src/mc_gestionnaire.cpp
Revision: 706
Committed: Mon Aug 31 13:35:48 2015 UTC (9 years, 8 months ago) by francois
File size: 9006 byte(s)
Log Message:
réérciture de la procédure pour ecrire des fichiers code aster. Changement -mecanique devient -elastique et un fichier de parametre obligatoire est maintenant necessaire pour faire un calcul

File Contents

# 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     // mt_gestionnaire.cpp
16     //
17     //------------------------------------------------------------
18     //------------------------------------------------------------
19     // COPYRIGHT 2000
20     // Version du 02/03/2006 � 11H23
21     //------------------------------------------------------------
22     //------------------------------------------------------------
23    
24    
25     #include "gestionversion.h"
26    
27     #pragma hdrstop
28    
29 francois 296 #include "mc_gestionnaire.h"
30     #include "mc_materiau.h"
31 francois 283 #include "parse.h"
32     #include "pars_argument.h"
33     #include <string.h>
34     //---------------------------------------------------------------------------
35 francois 296 MC_GESTIONNAIRE::MC_GESTIONNAIRE(char* path)
36 francois 283 {
37     chemin=path;
38     lire();
39     }
40    
41 francois 296
42     MC_GESTIONNAIRE::MC_GESTIONNAIRE()
43 francois 283 {
44 francois 296 std::pair<std::string,std::string> tmp("Dx","Déplacement en X");
45     mecaliste.insert(tmp);
46     std::pair<std::string,std::string> tmp2("Dy","Déplacement en Y");
47     mecaliste.insert(tmp2);
48     std::pair<std::string,std::string> tmp3("Dz","Déplacement en z");
49     mecaliste.insert(tmp3);
50 gervaislavoie 302 std::pair<std::string,std::string> tmp4("Da","Tous les déplacements et toutes les rotations");
51 francois 296 mecaliste.insert(tmp4);
52     std::pair<std::string,std::string> tmp5("Fx","Force concentrée en X");
53     mecaliste.insert(tmp5);
54     std::pair<std::string,std::string> tmp6("Fy","Force concentrée en Y");
55     mecaliste.insert(tmp6);
56     std::pair<std::string,std::string> tmp7("Fz","Force concentrée en Z");
57     mecaliste.insert(tmp7);
58     std::pair<std::string,std::string> tmp8("Px","Pression en x");
59     mecaliste.insert(tmp8);
60     std::pair<std::string,std::string> tmp9("Py","Pression en y");
61     mecaliste.insert(tmp9);
62     std::pair<std::string,std::string> tmp10("Pz","Pression en z");
63     mecaliste.insert(tmp10);
64     std::pair<std::string,std::string> tmp11("Pn","Pression normale");
65     mecaliste.insert(tmp11);
66 francois 706 /*std::pair<std::string,std::string> tmp12("Tp","Temperature de reference");
67     mecaliste.insert(tmp12);*/
68 francois 296 std::pair<std::string,std::string> tmp13("Cc","Couleur pour algorithme de comparaison");
69     mecaliste.insert(tmp13);
70 gervaislavoie 304 std::pair<std::string,std::string> tmp14("Em","Module d'Young isotrope");
71 francois 296 mecaliste.insert(tmp14);
72 gervaislavoie 304 std::pair<std::string,std::string> tmp15("nu","Coefficient de Poisson isotrope");
73 francois 296 mecaliste.insert(tmp15);
74     std::pair<std::string,std::string> tmp16("Ro","Densité du matériau");
75     mecaliste.insert(tmp16);
76 francois 300 std::pair<std::string,std::string> tmp17("ep","Epaisseur d'une coque");
77     mecaliste.insert(tmp17);
78 gervaislavoie 301 std::pair<std::string,std::string> tmp18("Rx","Rotation en x");
79     mecaliste.insert(tmp18);
80     std::pair<std::string,std::string> tmp19("Ry","Rotation en y");
81     mecaliste.insert(tmp19);
82     std::pair<std::string,std::string> tmp20("Rz","Rotation en z");
83     mecaliste.insert(tmp20);
84 gervaislavoie 302 std::pair<std::string,std::string> tmp21("Rt","Toutes les rotations");
85     mecaliste.insert(tmp21);
86     std::pair<std::string,std::string> tmp22("Dt","Tous les déplacements");
87     mecaliste.insert(tmp22);
88 francois 316 std::pair<std::string,std::string> tmp23("Cx","Abcisse d'orientation du point d'orientation d'une coque");
89     mecaliste.insert(tmp23);
90     std::pair<std::string,std::string> tmp24("Cy","Ordonnee d'orientation du point d'orientation d'une coque");
91     mecaliste.insert(tmp24);
92     std::pair<std::string,std::string> tmp25("Cz","Cote d'orientation du point d'orientation d'une coque");
93     mecaliste.insert(tmp25);
94 gervaislavoie 332 std::pair<std::string,std::string> tmp26("Nd","Topologie à ne pas optimiser");
95     mecaliste.insert(tmp26);
96 francois 576 std::pair<std::string,std::string> tmp27("Tn","Temparature à un noeud pour le problème thermique");
97     mecaliste.insert(tmp27);
98 francois 706 /*std::pair<std::string,std::string> tmp28("fx","Flux en x pour le problème thermique");
99 francois 576 mecaliste.insert(tmp28);
100 francois 592 std::pair<std::string,std::string> tmp29("fy","Flux en y pour le problème thermique");
101 francois 576 mecaliste.insert(tmp29);
102 francois 592 std::pair<std::string,std::string> tmp30("fz","Flux en z pour le problème thermique");
103 francois 706 mecaliste.insert(tmp30);*/
104 francois 592 std::pair<std::string,std::string> tmp31("fn","Flux normal pour le problème thermique");
105 francois 576 mecaliste.insert(tmp31);
106 francois 706 std::pair<std::string,std::string> tmp32("Kx","Conductivité thermique selon x ou isotrope");
107 francois 586 mecaliste.insert(tmp32);
108 francois 706 /*std::pair<std::string,std::string> tmp33("Ky","Conductivité thermique selon y");
109 francois 586 mecaliste.insert(tmp33);
110     std::pair<std::string,std::string> tmp34("Kz","Conductivité thermique selon z");
111 francois 706 mecaliste.insert(tmp34);*/
112 francois 586 std::pair<std::string,std::string> tmp35("Cp","Chaleur specifique");
113     mecaliste.insert(tmp35);
114 francois 706 /*
115 mejrim 597 std::pair<std::string,std::string> tmp36("ti","Temps initial");
116     mecaliste.insert(tmp36);
117     std::pair<std::string,std::string> tmp37("tf","Temps final");
118     mecaliste.insert(tmp37);
119     std::pair<std::string,std::string> tmp38("pa","pas de temps");
120     mecaliste.insert(tmp38);
121 mejrim 606 std::pair<std::string,std::string> tmp39("PE","Paramètre de la méthode d'EULER");
122 francois 706 mecaliste.insert(tmp39);*/
123 francois 296 }
124    
125    
126     int MC_GESTIONNAIRE::existe(std::string code)
127     {
128     std::map<std::string,std::string,std::less<std::string> >::iterator it=mecaliste.find(code);
129     if (it==mecaliste.end()) return 0;
130     return 1;
131     }
132    
133    
134     std::string MC_GESTIONNAIRE::get_description(std::string code)
135     {
136     std::map<std::string,std::string,std::less<std::string> >::iterator it=mecaliste.find(code);
137     if (it==mecaliste.end()) return "";
138     return (*it).second;
139     }
140    
141     int MC_GESTIONNAIRE::get_description(std::string &code,std::string &description,int init)
142     {
143     static std::map<std::string,std::string,std::less<std::string> >::iterator it;
144     static int passe=0;
145     if (passe==0) {init=1;passe=1;}
146     if (it==mecaliste.end()) init=1;
147     if (init==1) it=mecaliste.begin();
148     if (init==0) it++;
149     if (it!=mecaliste.end())
150     {
151     code=(*it).first;
152     description=(*it).second;
153     return 1;
154     }
155     return 0;
156     }
157    
158    
159     MC_GESTIONNAIRE::~MC_GESTIONNAIRE()
160     {
161 francois 283 int nb=lst_mat.get_nb();
162     while (nb!=0)
163     {
164 francois 296 MC_MATERIAU* mat=lst_mat.get(0);
165 francois 283 lst_mat.supprimer(mat);
166     delete mat;
167     nb=lst_mat.get_nb();
168     }
169     }
170    
171 francois 296 class MC_MATERIAU* MC_GESTIONNAIRE::ajouter_materiau(std::string nom,std::string descrip)
172 francois 283 {
173 francois 296 MC_MATERIAU* mat=new MC_MATERIAU((char*)nom.c_str(),(char*)descrip.c_str());
174 francois 283 lst_mat.ajouter(mat);
175     return mat;
176     }
177    
178 francois 296 int MC_GESTIONNAIRE::get_nb_materiau()
179 francois 283 {
180     return lst_mat.get_nb();
181     }
182    
183 francois 296 class MC_MATERIAU* MC_GESTIONNAIRE::get_materiau(int num)
184 francois 283 {
185     return lst_mat.get(num);
186     }
187    
188 francois 296 char* MC_GESTIONNAIRE::get_chemin()
189 francois 283 {
190     return (char*)chemin.c_str();
191     }
192    
193    
194 francois 296 void MC_GESTIONNAIRE::lire(void)
195 francois 283 {
196     FILE *in=fopen(chemin.c_str(),"rt");
197     if (in==NULL) return ;
198     PARS_ARGUMENT param[20];
199     PARSE parse;
200     char motcle[255];
201     do
202     {
203     char mess[255];
204     char* res=fgets(mess,255,in);
205     std::string ligne=mess;
206     parse.decode(ligne.c_str(),(char*)"@\n",param);
207     sscanf(param[0].argument[0].c_str(),"%s",motcle);
208     if (strcmp(motcle,"END")!=0)
209     {
210     parse.decode(ligne.c_str(),(char*)"@:@\n",param);
211     sscanf(param[0].argument[0].c_str(),"%s",motcle);
212     if (strcmp(motcle,"MATL")==0)
213     {
214     parse.decode(ligne.c_str(),(char*)"@:@:@:@\n",param);
215 francois 296 MC_MATERIAU* mate=ajouter_materiau(param[1].argument[0],param[3].argument[0].c_str());
216 francois 283 do
217     {
218     char mess[255];
219     char *res=fgets(mess,255,in);
220     std::string ligne=mess;
221     parse.decode(ligne.c_str(),(char*)"@\n",param);
222     sscanf(param[0].argument[0].c_str(),"%s",motcle);
223     if (strcmp(motcle,"ENDMATL")!=0)
224     {
225     char nom[20],uni1[50],uni2[50],uni3[50];
226     double val1,val2,val3;
227     int nb=sscanf(ligne.c_str(),"%s %lf %s %lf %s %lf %s",nom,&val1,uni1,&val2,uni2,&val3,uni3);
228     if (nb==7)
229     mate->ajouter_propriete(nom,val1,uni1,val2,uni2,val3,uni3);
230     else
231     {
232     sscanf(ligne.c_str(),"%s %lf %lf %lf",nom,&val1,&val2,&val3);
233     mate->ajouter_propriete(nom,val1,(char*)"",val2,(char*)"",val3,(char*)"");
234     }
235    
236     }
237     sscanf(param[0].argument[0].c_str(),"%s",motcle);
238     }
239     while (strcmp(motcle,"ENDMATL")!=0);
240     }
241     }
242     sscanf(param[0].argument[0].c_str(),"%s",motcle);
243     }
244     while (strcmp(motcle,"END")!=0);
245     fclose(in);
246     }
247    
248