ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/microstructure/src/mstruct_parametres.cpp
Revision: 930
Committed: Tue May 15 15:50:32 2018 UTC (7 years ago) by couturad
File size: 53560 byte(s)
Log Message: 
Mise a jour des fonctionnalites de MICROSTRUCTURE

File Contents

# User Rev Content
1 couturad 919 #include "mstruct_parametres.h"
2     #include "mstruct_generateur_rsa.h"
3 couturad 926 #include "mg_export.h"
4 couturad 919 #include <string.h>
5     #include <fstream>
6     #include <cstdlib>
7    
8     using namespace MICROSTRUCTURE;
9    
10     void PARAMETRES::ini_param_generateur_rsa(OT_PARAMETRES* params, int type_inclusion,char* nom_groupe_forme)
11     {
12     params->ajouter("Type_generateur","0",OT_PARAMETRES::DOUBLE,"Type de generateur : RSA [1], ");
13     params->ajouter("Nom_mgcg_modele","VES",OT_PARAMETRES::STRING,"Nom du MG_CG_MODELE");
14     params->ajouter("Nom_groupe_inclusion",nom_groupe_forme,OT_PARAMETRES::STRING,"Nom du groupe d'inclusion (Ex: Sphere)");
15     params->ajouter("Boite3D_distribution_Xmin","-0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D de distribution d'inclusion");
16     params->ajouter("Boite3D_distribution_Ymin","-0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D de distribution d'inclusion");
17     params->ajouter("Boite3D_distribution_Zmin","-0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D de distribution d'inclusion");
18     params->ajouter("Boite3D_distribution_Xmax","1.1",OT_PARAMETRES::DOUBLE,"Xmax boite3D de distribution d'inclusion");
19     params->ajouter("Boite3D_distribution_Ymax","1.1",OT_PARAMETRES::DOUBLE,"Ymax boite3D de distribution d'inclusion");
20     params->ajouter("Boite3D_distribution_Zmax","1.1",OT_PARAMETRES::DOUBLE,"Zmax boite3D de distribution d'inclusion");
21     params->ajouter("Nb_pas_X","20",OT_PARAMETRES::DOUBLE,"Nombre de pas selon X pout la creation d'un octree");
22     params->ajouter("Nb_pas_Y","20",OT_PARAMETRES::DOUBLE,"Nombre de pas selon Y pout la creation d'un octree");
23     params->ajouter("Nb_pas_Z","20",OT_PARAMETRES::DOUBLE,"Nombre de pas selon Z pout la creation d'un octree");
24 couturad 929 params->ajouter("Fraction_volumique_cible","0.05",OT_PARAMETRES::DOUBLE,"Fraction volumique cible d'inclusion [0.0, 1.0]");
25     params->ajouter("Eps_fraction_volumique","0.001",OT_PARAMETRES::DOUBLE,"Epsilon fraction volumique");
26 couturad 919 if(type_inclusion==GENERATEUR_RSA::TYPE_INCLUSION::SPHERE)
27     {
28     params->ajouter("Type_inclusion","0",OT_PARAMETRES::DOUBLE,"Type d'inclusion : Sphere [0], Cylindre [1]");
29     params->ajouter("Mu_rayon","0.2",OT_PARAMETRES::DOUBLE,"Moyenne du rayon de la sphere");
30     params->ajouter("Sigma_rayon","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du rayon de la sphere");
31     }
32     if(type_inclusion==GENERATEUR_RSA::TYPE_INCLUSION::CYLINDRE)
33     {
34     params->ajouter("Type_inclusion","1",OT_PARAMETRES::DOUBLE,"Type d'inclusion : Sphere [0], Cylindre [1]");
35 couturad 926 params->ajouter("Mu_rayon","0.05",OT_PARAMETRES::DOUBLE,"Moyenne du rayon du cylindre");
36 couturad 919 params->ajouter("Sigma_rayon","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du rayon du cylindre");
37 couturad 926 params->ajouter("Mu_longueur","0.4",OT_PARAMETRES::DOUBLE,"Moyenne de la longueur du cylindre");
38 couturad 919 params->ajouter("Sigma_longueur","0.0",OT_PARAMETRES::DOUBLE,"Ecart type de la longueur du cylindre");
39 couturad 929 params->ajouter("Mu_theta","0.0",OT_PARAMETRES::DOUBLE,"Moyenne de l'angle theta");
40     params->ajouter("Sigma_theta","0.0",OT_PARAMETRES::DOUBLE,"Ecart type de l'angle theta");
41     params->ajouter("Type_distribution_theta","0",OT_PARAMETRES::DOUBLE,"Type de distribution de l'angle theta : Uniforme[0], Normale[1]");
42     params->ajouter("Mu_phi","0.0",OT_PARAMETRES::DOUBLE,"Moyenne de l'angle phi");
43     params->ajouter("Sigma_phi","0.0",OT_PARAMETRES::DOUBLE,"Ecart type de l'angle phi");
44     params->ajouter("Type_distribution_phi","0",OT_PARAMETRES::DOUBLE,"Type de distribution de l'angle phi : Uniforme[0], Normale[1]");
45 couturad 919 }
46     params->ajouter("Porosite","0",OT_PARAMETRES::DOUBLE,"Considerer l'inclusion comme etant une porosite (vide) [1]");
47 couturad 926 params->ajouter("Distance_min_inter_volume","0.04",OT_PARAMETRES::DOUBLE,"Distance minimale qui separe deux inclusions");
48     params->ajouter("Volume_min","0.0008",OT_PARAMETRES::DOUBLE,"Volume minimal admissible de l'inclusion");
49     params->ajouter("Aire_min","0.001963",OT_PARAMETRES::DOUBLE,"Aire minimal admissible d'une face de l'inclusion");
50 couturad 919 params->ajouter("Longueur_min","0.02",OT_PARAMETRES::DOUBLE,"Longueur minimal admissible d'une arete ");
51     params->ajouter("Angle_min","0.7854",OT_PARAMETRES::DOUBLE,"Angle minimal admissible entre deux faces");
52     params->ajouter("Nb_volume_max","500",OT_PARAMETRES::DOUBLE,"Nombre maximal d'inclusion");
53     params->ajouter("Nb_iteration_max","2000",OT_PARAMETRES::DOUBLE,"Nombre maximal d'iteration");
54     }
55    
56     void PARAMETRES::ini_param_materiau(OT_PARAMETRES* params,char* nom_groupe_forme)
57     {
58     params->ajouter("Nom_groupe_forme",nom_groupe_forme,OT_PARAMETRES::STRING,"Nom du groupe de forme (Ex: Sphere, Matrice)");
59     params->ajouter("Mu_E","69.0e9",OT_PARAMETRES::DOUBLE,"Moyenne du module d'elasticite");
60     params->ajouter("Sigma_E","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du module d'elasticite");
61     params->ajouter("Mu_nu","0.3",OT_PARAMETRES::DOUBLE,"Moyenne du coefficient de Poisson");
62     params->ajouter("Sigma_nu","0.0",OT_PARAMETRES::DOUBLE,"Ecart type du coefficient de Poisson");
63     }
64    
65     void PARAMETRES::ini_param_carte_taille(OT_PARAMETRES* params)
66     {
67     params->ajouter("Nom_carte","carte",OT_PARAMETRES::STRING,"Nom de la carte de taille");
68 couturad 926 params->ajouter("Ecart_nodal","0.04",OT_PARAMETRES::DOUBLE,"Ecart nodal");
69 couturad 919 params->ajouter("Fechantillonnage","20",OT_PARAMETRES::DOUBLE,"Nombre d'échantillon par entité");
70 couturad 926 params->ajouter("Nb_cellule_direction","50",OT_PARAMETRES::DOUBLE,"Nombre de cellules de la FEM_SOLUTION dans les directions x,y,z");
71 couturad 919 params->ajouter("Nom_fem_solution","solcarte",OT_PARAMETRES::STRING,"Nom de la FEM_SOLUTION");
72     params->ajouter("Nb_couche_min","2",OT_PARAMETRES::DOUBLE,"Nombre de couches entre chaque solide");
73     params->ajouter("Nb_pas","32",OT_PARAMETRES::DOUBLE,"Nombre de pas de balayage");
74     params->ajouter("Facteur_augmentation","25",OT_PARAMETRES::DOUBLE,"Facteur d'augmentation de l'ecart nodal d'une cellule a l'autre");
75     }
76    
77     void PARAMETRES::ini_param_mailleur(OT_PARAMETRES* params)
78     {
79     params->ajouter("Niveau","3",OT_PARAMETRES::DOUBLE,"Niveau du maillage : 1D [1], 2D [2], 3D [3]");
80     params->ajouter("Niveau_opt_2d","5",OT_PARAMETRES::DOUBLE,"Limite supérieure des mailles à optimiser");
81     params->ajouter("Niveau_opt_3d","2",OT_PARAMETRES::DOUBLE,"Limite supérieure des mailles à optimiser");
82     params->ajouter("Priorite_metrique","0.65",OT_PARAMETRES::DOUBLE,"Valeur de la prise en compte de la métrique versus la qualité");
83     params->ajouter("Analyse","0",OT_PARAMETRES::DOUBLE,"Analyse de la qualité du maillage après le maillage : Analyse [1]");
84     }
85    
86     void PARAMETRES::ini_param_mailleur_fem(OT_PARAMETRES* params)
87     {
88     params->ajouter("Degre","2",OT_PARAMETRES::DOUBLE,"Degre du maillage a creer");
89     params->ajouter("Analyse","0",OT_PARAMETRES::DOUBLE,"Analyse de la qualité du maillage après le maillage : Analyse [1]");
90     params->ajouter("Optimisation_num_noeud","0",OT_PARAMETRES::DOUBLE,"Optimisation de la numerotation des noeuds : Optimisation [1]");
91     }
92    
93     void PARAMETRES::ini_param_etude(OT_PARAMETRES* params, int type_etude)
94     {
95     params->ajouter("Nom_etude","etude_1",OT_PARAMETRES::STRING,"Nom de l'etude");
96     params->ajouter("Fichier_param_aster","param_aster.txt",OT_PARAMETRES::STRING,"Nom du fichier de parametres aster");
97     params->ajouter("Type_etude",type_etude,OT_PARAMETRES::DOUBLE,"Type d'etude : Mecanique [0], Thermique [1]");
98     params->ajouter("Type_calcul","0",OT_PARAMETRES::DOUBLE,"Type de calcul : Elastique [0]");
99     params->ajouter("Code_resu","11111111",OT_PARAMETRES::STRING,"Code de resolution");
100 couturad 926 if(type_etude==MSTRUCT_VES::TYPE_ETUDE::MECANIQUE)
101 couturad 919 {
102     params->ajouter("Type_Chargement","0",OT_PARAMETRES::DOUBLE,"Type de chargement : Spherique [0], Deviatorique [1]");
103     params->ajouter("Type_CL","0",OT_PARAMETRES::DOUBLE,"Type de conditions aux limites : CLDH [0], CLCH[1]");
104     params->ajouter("Valeur_CL","0.0001",OT_PARAMETRES::DOUBLE,"Valeur de la condition aux limites");
105     }
106     }
107    
108     void PARAMETRES::ini_param_homogeneisation(OT_PARAMETRES* params)
109     {
110     params->ajouter("Fraction_volumique_CAD",0,OT_PARAMETRES::DOUBLE,"Calculer la fraction volumique CAD d'un groupe de forme [1]");
111     params->ajouter("Fraction_volumique_mg_maillage",0,OT_PARAMETRES::DOUBLE,"Calculer la fraction volumique du mg_maillage d'un groupe de forme [1]");
112     params->ajouter("Fraction_volumique_fem_maillage",0,OT_PARAMETRES::DOUBLE,"Calculer la fraction volumique du fem_maillage d'un groupe de forme [1]");
113     params->ajouter("Tenseur_orientation",0,OT_PARAMETRES::DOUBLE,"Calculer le tenseur d'orientation d'un groupe de forme [1]");
114     params->ajouter("Comparer_champ_theorie",0,OT_PARAMETRES::DOUBLE,"Comparer un champ par rapport a un champ theorique [1]");
115     params->ajouter("Critere_comparaison",0,OT_PARAMETRES::DOUBLE,"Erreur moyenne absolue [0], Erreur moyenne au carre [2]");
116     params->ajouter("Moyenne_volumique_champ",0,OT_PARAMETRES::DOUBLE,"Calculer la moyenne volumique d'un champ [1]");
117     params->ajouter("Num_solution",0,OT_PARAMETRES::DOUBLE,"Numero de la solution");
118     params->ajouter("Nom_champ","Sigma",OT_PARAMETRES::STRING,"Nom du champ");
119 couturad 926 params->ajouter("Nom_groupe_forme","ALL",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [ALL] prend le ves au complet");
120 couturad 919 params->ajouter("Boite_analyse",0,OT_PARAMETRES::DOUBLE,"Calculer à partir des elements donc le centre de la boite englobante figure dans la boite d'analyse");
121     params->ajouter("Boite3D_analyse_Xmin","0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D d'analyse");
122     params->ajouter("Boite3D_analyse_Ymin","0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D d'analyse");
123     params->ajouter("Boite3D_analyse_Zmin","0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D d'analyse");
124     params->ajouter("Boite3D_analyse_Xmax","0.9",OT_PARAMETRES::DOUBLE,"Xmax boite3D d'analyse");
125     params->ajouter("Boite3D_analyse_Ymax","0.9",OT_PARAMETRES::DOUBLE,"Ymax boite3D d'analyse");
126     params->ajouter("Boite3D_analyse_Zmax","0.9",OT_PARAMETRES::DOUBLE,"Zmax boite3D d'analyse");
127     params->ajouter("Exporter_liste_resultat","0",OT_PARAMETRES::DOUBLE,"Exporter une liste de resultats");
128     params->ajouter("Fichier_liste_resultat","fichier_liste_resultat.txt",OT_PARAMETRES::STRING,"Exporter une liste de resultats");
129     params->ajouter("Val_champ_theorique_XX","0.0",OT_PARAMETRES::DOUBLE,"Valeur d'un champ theorique XX");
130     params->ajouter("Val_champ_theorique_YY","0.0",OT_PARAMETRES::DOUBLE,"Valeur d'un champ theorique YY");
131     params->ajouter("Val_champ_theorique_ZZ","0.0",OT_PARAMETRES::DOUBLE,"Valeur d'un champ theorique ZZ");
132     params->ajouter("Val_champ_theorique_XY","0.0",OT_PARAMETRES::DOUBLE,"Valeur d'un champ theorique XY");
133     params->ajouter("Val_champ_theorique_XZ","0.0",OT_PARAMETRES::DOUBLE,"Valeur d'un champ theorique XZ");
134     params->ajouter("Val_champ_theorique_YZ","0.0",OT_PARAMETRES::DOUBLE,"Valeur d'un champ theorique YZ");
135     }
136    
137 couturad 930 void PARAMETRES::ini_param_analyse(OT_PARAMETRES* params,int type_analyse,bool erosion)
138 couturad 926 {
139     if(type_analyse==MSTRUCT_ANALYSE::TYPE_ANALYSE::CHAMP)
140     {
141     params->ajouter("Identifiant","Champ_SIGMA",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
142     params->ajouter("Type_analyse",0,OT_PARAMETRES::DOUBLE,"Type d'analyse : Champ [0], Orientation [1], CAO [2], MAILLAGE_MG [3], MAILLAGE_FEM [4]");
143     params->ajouter("Num_solution",2,OT_PARAMETRES::DOUBLE,"Numero de la FEM_SOLUTION solution (0,1,2,...)");
144     params->ajouter("Nom_groupe_forme","ALL",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [ALL] prend le ves au complet");
145     params->ajouter("Boite_analyse",0,OT_PARAMETRES::DOUBLE,"Calculer à partir des elements donc le centre de la boite englobante figure dans la boite d'analyse [1]");
146     params->ajouter("Boite3D_analyse_Xmin","0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D d'analyse");
147     params->ajouter("Boite3D_analyse_Ymin","0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D d'analyse");
148     params->ajouter("Boite3D_analyse_Zmin","0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D d'analyse");
149     params->ajouter("Boite3D_analyse_Xmax","0.9",OT_PARAMETRES::DOUBLE,"Xmax boite3D d'analyse");
150     params->ajouter("Boite3D_analyse_Ymax","0.9",OT_PARAMETRES::DOUBLE,"Ymax boite3D d'analyse");
151     params->ajouter("Boite3D_analyse_Zmax","0.9",OT_PARAMETRES::DOUBLE,"Zmax boite3D d'analyse");
152     params->ajouter("Largeur_colonne_distribution","0.1",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
153     }
154     else if(type_analyse==MSTRUCT_ANALYSE::TYPE_ANALYSE::ORIENTATION)
155     {
156     params->ajouter("Identifiant","Orientation_Cylindre",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
157     params->ajouter("Type_analyse",1,OT_PARAMETRES::DOUBLE,"Type d'analyse : Champ [0], Orientation [1], CAO [2], MAILLAGE_MG [3], MAILLAGE_FEM [4]");
158     params->ajouter("Nom_groupe_forme","Cylindre",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion.");
159     params->ajouter("Boite_analyse",0,OT_PARAMETRES::DOUBLE,"Calculer à partir des elements donc le centre de la boite englobante figure dans la boite d'analyse [1]");
160     params->ajouter("Boite3D_analyse_Xmin","0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D d'analyse");
161     params->ajouter("Boite3D_analyse_Ymin","0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D d'analyse");
162     params->ajouter("Boite3D_analyse_Zmin","0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D d'analyse");
163     params->ajouter("Boite3D_analyse_Xmax","0.9",OT_PARAMETRES::DOUBLE,"Xmax boite3D d'analyse");
164     params->ajouter("Boite3D_analyse_Ymax","0.9",OT_PARAMETRES::DOUBLE,"Ymax boite3D d'analyse");
165     params->ajouter("Boite3D_analyse_Zmax","0.9",OT_PARAMETRES::DOUBLE,"Zmax boite3D d'analyse");
166     }
167     else if(type_analyse==MSTRUCT_ANALYSE::TYPE_ANALYSE::CAO)
168     {
169     params->ajouter("Identifiant","Cao",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
170     params->ajouter("Type_analyse",2,OT_PARAMETRES::DOUBLE,"Type d'analyse : Champ [0], Orientation [1], CAO [2], MAILLAGE_MG [3], MAILLAGE_FEM [4]");
171     params->ajouter("Nom_groupe_forme","ALL",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [NULL] prend le ves au complet");
172     params->ajouter("Largeur_colonne_distribution_nb_volume","1.0",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
173     params->ajouter("Largeur_colonne_distribution_nb_forme","1.0",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
174     params->ajouter("Largeur_colonne_distribution_volume","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
175     params->ajouter("Largeur_colonne_distribution_fraction_volumique","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
176     params->ajouter("Largeur_colonne_distribution_volume_forme","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
177     }
178     else if(type_analyse==MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_MG)
179     {
180     params->ajouter("Identifiant","Maillage_MG",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
181     params->ajouter("Type_analyse",3,OT_PARAMETRES::DOUBLE,"Type d'analyse : Champ [0], Orientation [1], CAO [2], MAILLAGE_MG [3], MAILLAGE_FEM [4]");
182     params->ajouter("Nom_groupe_forme","ALL",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [ALL] prend le ves au complet");
183     params->ajouter("Boite_analyse",0,OT_PARAMETRES::DOUBLE,"Calculer à partir des elements donc le centre de la boite englobante figure dans la boite d'analyse [1]");
184     params->ajouter("Boite3D_analyse_Xmin","0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D d'analyse");
185     params->ajouter("Boite3D_analyse_Ymin","0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D d'analyse");
186     params->ajouter("Boite3D_analyse_Zmin","0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D d'analyse");
187     params->ajouter("Boite3D_analyse_Xmax","0.9",OT_PARAMETRES::DOUBLE,"Xmax boite3D d'analyse");
188     params->ajouter("Boite3D_analyse_Ymax","0.9",OT_PARAMETRES::DOUBLE,"Ymax boite3D d'analyse");
189     params->ajouter("Boite3D_analyse_Zmax","0.9",OT_PARAMETRES::DOUBLE,"Zmax boite3D d'analyse");
190     params->ajouter("Largeur_colonne_distribution_nb_element_2D","100",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
191     params->ajouter("Largeur_colonne_distribution_nb_element_3D","100",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
192     params->ajouter("Largeur_colonne_distribution_qualite_2D","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
193     params->ajouter("Largeur_colonne_distribution_qualite_3D","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
194     params->ajouter("Largeur_colonne_distribution_taille_2D","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
195     params->ajouter("Largeur_colonne_distribution_taille_3D","0.0001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
196     }
197     else if(type_analyse==MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_FEM)
198     {
199     params->ajouter("Identifiant","Maillage_FEM",OT_PARAMETRES::STRING,"Identifiant de l'analyse");
200     params->ajouter("Type_analyse",4,OT_PARAMETRES::DOUBLE,"Type d'analyse : Champ [0], Orientation [1], CAO [2], MAILLAGE_MG [3], MAILLAGE_FEM [4]");
201     params->ajouter("Nom_groupe_forme","ALL",OT_PARAMETRES::STRING,"Nom du groupe d'inclusion. Si [ALL] prend le ves au complet");
202     params->ajouter("Boite_analyse",0,OT_PARAMETRES::DOUBLE,"Calculer à partir des elements donc le centre de la boite englobante figure dans la boite d'analyse [1]");
203     params->ajouter("Boite3D_analyse_Xmin","0.1",OT_PARAMETRES::DOUBLE,"Xmin boite3D d'analyse");
204     params->ajouter("Boite3D_analyse_Ymin","0.1",OT_PARAMETRES::DOUBLE,"Ymin boite3D d'analyse");
205     params->ajouter("Boite3D_analyse_Zmin","0.1",OT_PARAMETRES::DOUBLE,"Zmin boite3D d'analyse");
206     params->ajouter("Boite3D_analyse_Xmax","0.9",OT_PARAMETRES::DOUBLE,"Xmax boite3D d'analyse");
207     params->ajouter("Boite3D_analyse_Ymax","0.9",OT_PARAMETRES::DOUBLE,"Ymax boite3D d'analyse");
208     params->ajouter("Boite3D_analyse_Zmax","0.9",OT_PARAMETRES::DOUBLE,"Zmax boite3D d'analyse");
209     params->ajouter("Largeur_colonne_distribution_nb_element_2D","100",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
210     params->ajouter("Largeur_colonne_distribution_nb_element_3D","100",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
211     params->ajouter("Largeur_colonne_distribution_jacobien_2D_min","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
212     params->ajouter("Largeur_colonne_distribution_jacobien_2D_max","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
213     params->ajouter("Largeur_colonne_distribution_jacobien_3D_min","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
214     params->ajouter("Largeur_colonne_distribution_jacobien_3D_max","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
215     params->ajouter("Largeur_colonne_distribution_distortion_2D","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
216     params->ajouter("Largeur_colonne_distribution_distortion_3D","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
217 couturad 930 }
218     if(erosion)
219     {
220     params->ajouter("Analyse_erosion","1",OT_PARAMETRES::DOUBLE,"Analyse d'erosion [1], sans erosion [0]");
221     params->ajouter("Nb_couche","40",OT_PARAMETRES::DOUBLE,"Nombre de couche d'erosion");
222     params->ajouter("Epaisseur_couche","0.01",OT_PARAMETRES::DOUBLE,"Epaisseur d'une couche d'erosion");
223 couturad 926 }
224 couturad 930 if(!erosion)
225     {
226     params->ajouter("Analyse_erosion","0",OT_PARAMETRES::DOUBLE,"Analyse d'erosion [1], sans erosion [0]");
227     }
228 couturad 926 }
229    
230     void PARAMETRES::ini_param_ver(OT_PARAMETRES* params)
231     {
232     params->ajouter("Identifiant_epsilon","Champ_EPSILON",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ epsilon");
233     params->ajouter("Identifiant_sigma","Champ_SIGMA",OT_PARAMETRES::STRING,"Identifiant de l'analyse du champ sigma");
234 couturad 927 params->ajouter("Largeur_colonne_distribution_module_Kapp","100000",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
235     params->ajouter("Largeur_colonne_distribution_module_Gapp","100000",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
236     params->ajouter("Largeur_colonne_distribution_module_Eapp","100000",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
237     params->ajouter("Largeur_colonne_distribution_module_Nuapp","0.001",OT_PARAMETRES::DOUBLE,"Largeur d'une colonne pour graph de distribution");
238 couturad 930 params->ajouter("Analyse_erosion","0",OT_PARAMETRES::DOUBLE,"Analyse d'erosion [1], sans erosion [0]");
239 couturad 926 }
240    
241    
242    
243    
244 couturad 919 int PARAMETRES::generer_fichier_parametres(char* nom_dossier_param)
245     {
246     int mkdir_err=0;
247     std::string sys_command = "mkdir "; sys_command.append(nom_dossier_param);
248     mkdir_err = system(sys_command.c_str());
249     if(mkdir_err !=0)
250     {
251     std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametres! ***" << std::endl;
252     return FAIL;
253     }
254     sys_command.clear(); sys_command.append("mkdir -p ");sys_command.append(nom_dossier_param); sys_command.append("/geometrie");
255     mkdir_err = system(sys_command.c_str());
256     if(mkdir_err !=0)
257     {
258     std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametres/geometrie ! ***" << std::endl;
259     return FAIL;
260     }
261     sys_command.clear(); sys_command.append("mkdir -p ");sys_command.append(nom_dossier_param); sys_command.append("/materiau");
262     mkdir_err = system(sys_command.c_str());
263     if(mkdir_err !=0)
264     {
265     std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/materiau ! ***" << std::endl;
266     return FAIL;
267     }
268     sys_command.clear(); sys_command.append("mkdir -p ");sys_command.append(nom_dossier_param); sys_command.append("/carte");
269     mkdir_err = system(sys_command.c_str());
270     if(mkdir_err !=0)
271     {
272     std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametres/carte ! ***" << std::endl;
273     return FAIL;
274     }
275     sys_command.clear(); sys_command.append("mkdir -p ");sys_command.append(nom_dossier_param); sys_command.append("/maillage");
276     mkdir_err = system(sys_command.c_str());
277     if(mkdir_err !=0)
278     {
279     std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/maillage ! ***" << std::endl;
280     return FAIL;
281     }
282     sys_command.clear(); sys_command.append("mkdir -p ");sys_command.append(nom_dossier_param); sys_command.append("/fem_maillage");
283     mkdir_err = system(sys_command.c_str());
284     if(mkdir_err !=0)
285     {
286     std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/fem_maillage ! ***" << std::endl;
287     return FAIL;
288     }
289     sys_command.clear(); sys_command.append("mkdir -p ");sys_command.append(nom_dossier_param); sys_command.append("/etude");
290     mkdir_err = system(sys_command.c_str());
291     if(mkdir_err !=0)
292     {
293     std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/etude ! ***" << std::endl;
294     return FAIL;
295     }
296 couturad 926 sys_command.clear(); sys_command.append("mkdir -p ");sys_command.append(nom_dossier_param); sys_command.append("/analyse");
297 couturad 919 mkdir_err = system(sys_command.c_str());
298     if(mkdir_err !=0)
299     {
300 couturad 926 std::cerr << "*** PARAMETRES::generer_fichier_parametres -> Erreur lors de la creation des dossiers parametre/analyse ! ***" << std::endl;
301 couturad 919 return FAIL;
302     }
303     std::ofstream f;
304     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/geometrie/lstparam_geo.txt");
305     f.open(sys_command.c_str(),ios::out);
306     f << "// Fichier de parametres GEOMETRIE" << std::endl;
307     f << "// [Nombre de fichier de parametres]" << std::endl;
308 couturad 926 f << "NB_FICHIER_PARAM 1" << std::endl;
309 couturad 919 f << "param_geo_cylindre.txt" << std::endl;
310     f.close();
311     OT_PARAMETRES params_geo_1;
312 couturad 926 ini_param_generateur_rsa(&params_geo_1,GENERATEUR_RSA::TYPE_INCLUSION::CYLINDRE,(char*)"Cylindre");
313     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/geometrie/param_geo_cylindre.txt");
314 couturad 919 params_geo_1.enregistrer((char*)sys_command.c_str());
315    
316     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/materiau/lstparam_mat.txt");
317     f.open(sys_command.c_str(),ios::out);
318     f << "// Fichier de parametres MATERIAU" << std::endl;
319     f << "// [Nombre de fichier de parametres]" << std::endl;
320 couturad 926 f << "NB_FICHIER_PARAM 2" << std::endl;
321 couturad 919 f << "param_mat_cylindre.txt" << std::endl;
322     f << "param_mat_matrice.txt" << std::endl;
323     f.close();
324     OT_PARAMETRES params_mat_cylindre;
325     ini_param_materiau(&params_mat_cylindre,(char*)"Cylindre");
326 couturad 926 params_mat_cylindre.change_valeur("Mu_E",75.0e9);
327     params_mat_cylindre.change_valeur("Mu_nu",0.22);
328 couturad 919 sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/materiau/param_mat_cylindre.txt");
329     params_mat_cylindre.enregistrer((char*)sys_command.c_str());
330     OT_PARAMETRES params_mat_matrice;
331     ini_param_materiau(&params_mat_matrice,(char*)"Matrice");
332 couturad 926 params_mat_matrice.change_valeur("Mu_E",3.5e9);
333     params_mat_matrice.change_valeur("Mu_nu",0.33);
334 couturad 919 sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/materiau/param_mat_matrice.txt");
335     params_mat_matrice.enregistrer((char*)sys_command.c_str());
336    
337     OT_PARAMETRES params_carte;
338     ini_param_carte_taille(&params_carte);
339     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/carte/param_carte.txt");
340     params_carte.enregistrer((char*)sys_command.c_str());
341    
342     OT_PARAMETRES params_maill;
343     ini_param_mailleur(&params_maill);
344     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/maillage/param_maill.txt");
345     params_maill.enregistrer((char*)sys_command.c_str());
346    
347     OT_PARAMETRES params_fem_maillage;
348     ini_param_mailleur_fem(&params_fem_maillage);
349     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/fem_maillage/param_fem_maill.txt");
350     params_fem_maillage.enregistrer((char*)sys_command.c_str());
351    
352 couturad 926 OT_PARAMETRES params_etude_CLDH_sph;
353     ini_param_etude(&params_etude_CLDH_sph,MSTRUCT_VES::TYPE_ETUDE::MECANIQUE);
354     params_etude_CLDH_sph.change_valeur("Nom_etude","etude_CLDH_sph");
355     params_etude_CLDH_sph.change_valeur("Type_Chargement",0);
356     params_etude_CLDH_sph.change_valeur("Type_CL",0);
357     params_etude_CLDH_sph.change_valeur("Valeur_CL",0.0001);
358     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/etude/param_etude_CLDH_sph.txt");
359     params_etude_CLDH_sph.enregistrer((char*)sys_command.c_str());
360 couturad 919
361 couturad 926 OT_PARAMETRES params_etude_CLCH_sph;
362     ini_param_etude(&params_etude_CLCH_sph,MSTRUCT_VES::TYPE_ETUDE::MECANIQUE);
363     params_etude_CLCH_sph.change_valeur("Nom_etude","etude_CLCH_sph");
364     params_etude_CLCH_sph.change_valeur("Type_Chargement",0);
365     params_etude_CLCH_sph.change_valeur("Type_CL",1);
366     params_etude_CLCH_sph.change_valeur("Valeur_CL",1000);
367     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/etude/param_etude_CLCH_sph.txt");
368     params_etude_CLCH_sph.enregistrer((char*)sys_command.c_str());
369    
370     OT_PARAMETRES params_etude_CLDH_dev;
371     ini_param_etude(&params_etude_CLDH_dev,MSTRUCT_VES::TYPE_ETUDE::MECANIQUE);
372     params_etude_CLDH_dev.change_valeur("Nom_etude","etude_CLDH_dev");
373     params_etude_CLDH_dev.change_valeur("Type_Chargement",1);
374     params_etude_CLDH_dev.change_valeur("Type_CL",0);
375     params_etude_CLDH_dev.change_valeur("Valeur_CL",0.0001);
376     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/etude/param_etude_CLDH_dev.txt");
377     params_etude_CLDH_dev.enregistrer((char*)sys_command.c_str());
378    
379     OT_PARAMETRES params_etude_CLCH_dev;
380     ini_param_etude(&params_etude_CLCH_dev,MSTRUCT_VES::TYPE_ETUDE::MECANIQUE);
381     params_etude_CLCH_dev.change_valeur("Nom_etude","etude_CLCH_dev");
382     params_etude_CLCH_dev.change_valeur("Type_Chargement",1);
383     params_etude_CLCH_dev.change_valeur("Type_CL",1);
384     params_etude_CLCH_dev.change_valeur("Valeur_CL",1000);
385     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/etude/param_etude_CLCH_dev.txt");
386     params_etude_CLCH_dev.enregistrer((char*)sys_command.c_str());
387    
388     MG_EXPORT exp;
389     OT_PARAMETRES *params_aster=exp.get_param_aster();
390     params_aster->change_valeur("Memoire",3072);
391     params_aster->change_valeur("Temps_max",250000);
392     params_aster->change_valeur("Nb_CPU",8);
393     params_aster->change_valeur("Noeud_ele",2);
394 couturad 929 params_aster->change_valeur("Base_num_mail",36);
395 couturad 926 sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/etude/param_aster.txt");
396     params_aster->enregistrer((char*)sys_command.c_str());
397    
398     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/lstparam_analyse.txt");
399 couturad 919 f.open(sys_command.c_str(),ios::out);
400 couturad 926 f << "// Fichier de parametres ANALYSE" << std::endl;
401 couturad 919 f << "// Nombre de fichier de parametres" << std::endl;
402 couturad 930 f << "NB_FICHIER_PARAM 21" << std::endl;
403 couturad 926 f << "params_analyse_champ_sigma.txt" << std::endl;
404     f << "params_analyse_champ_sigma_boite.txt" << std::endl;
405 couturad 930 f << "params_analyse_champ_sigma_erosion.txt" << std::endl;
406 couturad 926 f << "params_analyse_champ_epsilon.txt" << std::endl;
407     f << "params_analyse_champ_epsilon_boite.txt" << std::endl;
408 couturad 930 f << "params_analyse_champ_epsilon_erosion.txt" << std::endl;
409 couturad 926 f << "params_analyse_orientation.txt" << std::endl;
410     f << "params_analyse_orientation_boite.txt" << std::endl;
411 couturad 930 f << "params_analyse_orientation_erosion.txt" << std::endl;
412 couturad 926 f << "params_analyse_cao_all.txt" << std::endl;
413     f << "params_analyse_cao_cylindre.txt" << std::endl;
414     f << "params_analyse_cao_matrice.txt" << std::endl;
415     f << "params_analyse_mg_maillage_all.txt" << std::endl;
416     f << "params_analyse_mg_maillage_cylindre.txt" << std::endl;
417     f << "params_analyse_mg_maillage_matrice.txt" << std::endl;
418     f << "params_analyse_fem_maillage_all.txt" << std::endl;
419     f << "params_analyse_fem_maillage_all_boite.txt" << std::endl;
420 couturad 930 f << "params_analyse_fem_maillage_all_erosion.txt" << std::endl;
421 couturad 926 f << "params_analyse_fem_maillage_cylindre.txt" << std::endl;
422 couturad 930 f << "params_analyse_fem_maillage_cylindre_erosion.txt" << std::endl;
423 couturad 926 f << "params_analyse_fem_maillage_matrice.txt" << std::endl;
424 couturad 919 f.close();
425 couturad 926 OT_PARAMETRES params_analyse_champ_sigma;
426     ini_param_analyse(&params_analyse_champ_sigma,MSTRUCT_ANALYSE::TYPE_ANALYSE::CHAMP);
427     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_champ_sigma.txt");
428     params_analyse_champ_sigma.change_valeur("Identifiant","Champ_SIGMA");
429     params_analyse_champ_sigma.enregistrer((char*)sys_command.c_str());
430     OT_PARAMETRES params_analyse_champ_sigma_boite;
431     ini_param_analyse(&params_analyse_champ_sigma_boite,MSTRUCT_ANALYSE::TYPE_ANALYSE::CHAMP);
432     params_analyse_champ_sigma_boite.change_valeur("Boite_analyse",1);
433     params_analyse_champ_sigma_boite.change_valeur("Identifiant","Champ_SIGMA_boite");
434     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_champ_sigma_boite.txt");
435     params_analyse_champ_sigma_boite.enregistrer((char*)sys_command.c_str());
436 couturad 930 OT_PARAMETRES params_analyse_champ_sigma_erosion;
437     ini_param_analyse(&params_analyse_champ_sigma_erosion,MSTRUCT_ANALYSE::TYPE_ANALYSE::CHAMP,true);
438     params_analyse_champ_sigma_erosion.change_valeur("Identifiant","Champ_SIGMA_erosion");
439     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_champ_sigma_erosion.txt");
440     params_analyse_champ_sigma_erosion.enregistrer((char*)sys_command.c_str());
441 couturad 919
442 couturad 930
443 couturad 926 OT_PARAMETRES params_analyse_champ_epsilon;
444     ini_param_analyse(&params_analyse_champ_epsilon,MSTRUCT_ANALYSE::TYPE_ANALYSE::CHAMP);
445     params_analyse_champ_epsilon.change_valeur("Identifiant","Champ_EPSILON");
446     params_analyse_champ_epsilon.change_valeur("Num_solution",1);
447     params_analyse_champ_epsilon.change_valeur("Largeur_colonne_distribution",0.00000000001);
448     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_champ_epsilon.txt");
449     params_analyse_champ_epsilon.enregistrer((char*)sys_command.c_str());
450     OT_PARAMETRES params_analyse_champ_epsilon_boite;
451     ini_param_analyse(&params_analyse_champ_epsilon_boite,MSTRUCT_ANALYSE::TYPE_ANALYSE::CHAMP);
452     params_analyse_champ_epsilon_boite.change_valeur("Identifiant","Champ_EPSILON_boite");
453     params_analyse_champ_epsilon_boite.change_valeur("Num_solution",1);
454     params_analyse_champ_epsilon_boite.change_valeur("Boite_analyse",1);
455     params_analyse_champ_epsilon_boite.change_valeur("Largeur_colonne_distribution",0.00000000001);
456     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_champ_epsilon_boite.txt");
457     params_analyse_champ_epsilon_boite.enregistrer((char*)sys_command.c_str());
458 couturad 930 OT_PARAMETRES params_analyse_champ_epsilon_erosion;
459     ini_param_analyse(&params_analyse_champ_epsilon_erosion,MSTRUCT_ANALYSE::TYPE_ANALYSE::CHAMP,true);
460     params_analyse_champ_epsilon_erosion.change_valeur("Identifiant","Champ_EPSILON_erosion");
461     params_analyse_champ_epsilon_erosion.change_valeur("Num_solution",1);
462     params_analyse_champ_epsilon_erosion.change_valeur("Largeur_colonne_distribution",0.00000000001);
463     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_champ_epsilon_erosion.txt");
464     params_analyse_champ_epsilon_erosion.enregistrer((char*)sys_command.c_str());
465 couturad 926
466     OT_PARAMETRES params_analyse_orientation;
467     ini_param_analyse(&params_analyse_orientation,MSTRUCT_ANALYSE::TYPE_ANALYSE::ORIENTATION);
468     params_analyse_orientation.change_valeur("Identifiant","Orientation_Cylindre");
469     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_orientation.txt");
470     params_analyse_orientation.enregistrer((char*)sys_command.c_str());
471     OT_PARAMETRES params_analyse_orientation_boite;
472     ini_param_analyse(&params_analyse_orientation_boite,MSTRUCT_ANALYSE::TYPE_ANALYSE::ORIENTATION);
473     params_analyse_orientation_boite.change_valeur("Identifiant","Orientation_Cylindre_boite");
474     params_analyse_orientation_boite.change_valeur("Boite_analyse",1);
475     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_orientation_boite.txt");
476     params_analyse_orientation_boite.enregistrer((char*)sys_command.c_str());
477 couturad 930 OT_PARAMETRES params_analyse_orientation_erosion;
478     ini_param_analyse(&params_analyse_orientation_erosion,MSTRUCT_ANALYSE::TYPE_ANALYSE::ORIENTATION,true);
479     params_analyse_orientation_erosion.change_valeur("Identifiant","Orientation_Cylindre_erosion");
480     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_orientation_erosion.txt");
481     params_analyse_orientation_erosion.enregistrer((char*)sys_command.c_str());
482 couturad 926
483    
484     OT_PARAMETRES params_analyse_cao_all;
485     ini_param_analyse(&params_analyse_cao_all,MSTRUCT_ANALYSE::TYPE_ANALYSE::CAO);
486     params_analyse_cao_all.change_valeur("Identifiant","Cao_ALL");
487     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_cao_all.txt");
488     params_analyse_cao_all.enregistrer((char*)sys_command.c_str());
489     OT_PARAMETRES params_analyse_cao_cylindre;
490     ini_param_analyse(&params_analyse_cao_cylindre,MSTRUCT_ANALYSE::TYPE_ANALYSE::CAO);
491     params_analyse_cao_cylindre.change_valeur("Identifiant","Cao_Cylindre");
492     params_analyse_cao_cylindre.change_valeur("Nom_groupe_forme","Cylindre");
493     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_cao_cylindre.txt");
494     params_analyse_cao_cylindre.enregistrer((char*)sys_command.c_str());
495     OT_PARAMETRES params_analyse_cao_matrice;
496     ini_param_analyse(&params_analyse_cao_matrice,MSTRUCT_ANALYSE::TYPE_ANALYSE::CAO);
497     params_analyse_cao_matrice.change_valeur("Identifiant","Cao_Matrice");
498     params_analyse_cao_matrice.change_valeur("Nom_groupe_forme","Matrice");
499     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_cao_matrice.txt");
500     params_analyse_cao_matrice.enregistrer((char*)sys_command.c_str());
501    
502     OT_PARAMETRES params_analyse_mg_maillage_all;
503     ini_param_analyse(&params_analyse_mg_maillage_all,MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_MG);
504     params_analyse_mg_maillage_all.change_valeur("Identifiant","MG_maillage_ALL");
505     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_mg_maillage_all.txt");
506     params_analyse_mg_maillage_all.enregistrer((char*)sys_command.c_str());
507     OT_PARAMETRES params_analyse_mg_maillage_cylindre;
508     ini_param_analyse(&params_analyse_mg_maillage_cylindre,MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_MG);
509     params_analyse_mg_maillage_cylindre.change_valeur("Identifiant","MG_maillage_Cylindre");
510     params_analyse_mg_maillage_cylindre.change_valeur("Nom_groupe_forme","Cylindre");
511     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_mg_maillage_cylindre.txt");
512     params_analyse_mg_maillage_cylindre.enregistrer((char*)sys_command.c_str());
513     OT_PARAMETRES params_analyse_mg_maillage_matrice;
514     ini_param_analyse(&params_analyse_mg_maillage_matrice,MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_MG);
515     params_analyse_mg_maillage_matrice.change_valeur("Identifiant","MG_maillage_Matrice");
516     params_analyse_mg_maillage_matrice.change_valeur("Nom_groupe_forme","Matrice");
517     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_mg_maillage_matrice.txt");
518     params_analyse_mg_maillage_matrice.enregistrer((char*)sys_command.c_str());
519    
520     OT_PARAMETRES params_analyse_fem_maillage_all;
521     ini_param_analyse(&params_analyse_fem_maillage_all,MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_FEM);
522     params_analyse_fem_maillage_all.change_valeur("Identifiant","FEM_maillage_ALL");
523     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_fem_maillage_all.txt");
524     params_analyse_fem_maillage_all.enregistrer((char*)sys_command.c_str());
525     OT_PARAMETRES params_analyse_fem_maillage_all_boite;
526     ini_param_analyse(&params_analyse_fem_maillage_all_boite,MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_FEM);
527     params_analyse_fem_maillage_all_boite.change_valeur("Identifiant","FEM_maillage_ALL_boite");
528     params_analyse_fem_maillage_all_boite.change_valeur("Boite_analyse",1);
529     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_fem_maillage_all_boite.txt");
530 couturad 930 params_analyse_fem_maillage_all_boite.enregistrer((char*)sys_command.c_str());
531     OT_PARAMETRES params_analyse_fem_maillage_all_erosion;
532     ini_param_analyse(&params_analyse_fem_maillage_all_erosion,MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_FEM,true);
533     params_analyse_fem_maillage_all_erosion.change_valeur("Identifiant","FEM_maillage_ALL_erosion");
534     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_fem_maillage_all_erosion.txt");
535     params_analyse_fem_maillage_all_erosion.enregistrer((char*)sys_command.c_str());
536 couturad 926 OT_PARAMETRES params_analyse_fem_maillage_cylindre;
537     ini_param_analyse(&params_analyse_fem_maillage_cylindre,MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_FEM);
538     params_analyse_fem_maillage_cylindre.change_valeur("Identifiant","FEM_maillage_Cylindre");
539     params_analyse_fem_maillage_cylindre.change_valeur("Nom_groupe_forme","Cylindre");
540     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_fem_maillage_cylindre.txt");
541 couturad 930 params_analyse_fem_maillage_cylindre.enregistrer((char*)sys_command.c_str());
542     OT_PARAMETRES params_analyse_fem_maillage_cylindre_erosion;
543     ini_param_analyse(&params_analyse_fem_maillage_cylindre_erosion,MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_FEM,true);
544     params_analyse_fem_maillage_cylindre_erosion.change_valeur("Identifiant","FEM_maillage_Cylindre_erosion");
545     params_analyse_fem_maillage_cylindre_erosion.change_valeur("Nom_groupe_forme","Cylindre");
546     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_fem_maillage_cylindre_erosion.txt");
547     params_analyse_fem_maillage_cylindre_erosion.enregistrer((char*)sys_command.c_str());
548 couturad 926 OT_PARAMETRES params_analyse_fem_maillage_matrice;
549     ini_param_analyse(&params_analyse_fem_maillage_matrice,MSTRUCT_ANALYSE::TYPE_ANALYSE::MAILLAGE_FEM);
550     params_analyse_fem_maillage_matrice.change_valeur("Identifiant","FEM_maillage_Matrice");
551     params_analyse_fem_maillage_matrice.change_valeur("Nom_groupe_forme","Matrice");
552     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_analyse_fem_maillage_matrice.txt");
553     params_analyse_fem_maillage_matrice.enregistrer((char*)sys_command.c_str());
554    
555     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/lstparam_ver.txt");
556     f.open(sys_command.c_str(),ios::out);
557     f << "// Fichier de parametres VER" << std::endl;
558     f << "// Nombre de fichier de parametres" << std::endl;
559 couturad 930 f << "NB_FICHIER_PARAM 3" << std::endl;
560 couturad 926 f << "params_ver.txt" << std::endl;
561     f << "params_ver_boite.txt" << std::endl;
562 couturad 930 f << "params_ver_erosion.txt" << std::endl;
563 couturad 926 f.close();
564    
565     OT_PARAMETRES params_ver;
566     ini_param_ver(&params_ver);
567     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_ver.txt");
568     params_ver.enregistrer((char*)sys_command.c_str());
569    
570     OT_PARAMETRES params_ver_boite;
571     ini_param_ver(&params_ver_boite);
572     params_ver_boite.change_valeur("Identifiant_epsilon","Champ_EPSILON_boite");
573     params_ver_boite.change_valeur("Identifiant_sigma","Champ_SIGMA_boite");
574     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_ver_boite.txt");
575     params_ver_boite.enregistrer((char*)sys_command.c_str());
576 couturad 930
577     OT_PARAMETRES params_ver_erosion;
578     ini_param_ver(&params_ver_erosion);
579     params_ver_erosion.change_valeur("Identifiant_epsilon","Champ_EPSILON_erosion");
580     params_ver_erosion.change_valeur("Identifiant_sigma","Champ_SIGMA_erosion");
581     params_ver_erosion.change_valeur("Analyse_erosion","1");
582     sys_command.clear(); sys_command.append(nom_dossier_param); sys_command.append("/analyse/params_ver_erosion.txt");
583     params_ver_erosion.enregistrer((char*)sys_command.c_str());
584 couturad 919 return OK;
585     }
586    
587 couturad 929 int PARAMETRES::generer_script_ves(int type_etude)
588     {
589     if(type_etude==MSTRUCT_VES::TYPE_ETUDE::MECANIQUE)
590     {
591     ofstream ofstrm((char*)"generer_nb_ves",std::ios::out|std::ios::trunc);
592     ofstrm.precision(16);
593     ofstrm.setf(std::ios::showpoint);
594     ofstrm << "#!/bin/bash" << std::endl;
595     ofstrm << "INDEX_DEPART=$1" << std::endl;
596     ofstrm << "INDEX_FIN=$2" << std::endl;
597     ofstrm << "DOSSIER_PARAM=$3" << std::endl;
598     ofstrm << "CHEMINCOMPILE=$4" << std::endl;
599     ofstrm << "if [ -z ${INDEX_DEPART// } ]" << std::endl;
600     ofstrm << "then" << std::endl;
601     ofstrm << " echo \"Erreur : Numero du ves de depart manquant\"" << std::endl;
602     ofstrm << " exit" << std::endl;
603     ofstrm << "fi" << std::endl;
604     ofstrm << "if [ -z ${INDEX_FIN// } ]" << std::endl;
605     ofstrm << "then" << std::endl;
606     ofstrm << " echo \"Erreur : Numero du ves de fin manquant\"" << std::endl;
607     ofstrm << " exit" << std::endl;
608     ofstrm << "fi" << std::endl;
609     ofstrm << "if [ -z ${DOSSIER_PARAM// } ]" << std::endl;
610     ofstrm << "then" << std::endl;
611     ofstrm << " echo \"Erreur : Dossier de parametres manquant\"" << std::endl;
612     ofstrm << " exit" << std::endl;
613     ofstrm << "fi" << std::endl;
614     ofstrm << "if [ -z ${CHEMINCOMPILE// } ]" << std::endl;
615     ofstrm << "then" << std::endl;
616     ofstrm << " echo \"Erreur : Chemin vers exe/app manquant\"" << std::endl;
617     ofstrm << " exit" << std::endl;
618     ofstrm << "fi" << std::endl;
619     ofstrm << "function executer" << std::endl;
620     ofstrm << "{" << std::endl;
621     ofstrm << "OPERATION_REUSSIE=1;" << std::endl;
622     ofstrm << "$*" << std::endl;
623     ofstrm << "VAR=$?" << std::endl;
624     ofstrm << "if (( $VAR )) ; then " << std::endl;
625     ofstrm << " OPERATION_REUSSIE=0; " << std::endl;
626     ofstrm << " STAT=$VAR;" << std::endl;
627     ofstrm << "fi" << std::endl;
628     ofstrm << "}" << std::endl;
629     ofstrm << "" << std::endl;
630     ofstrm << "echo \"*********************************************************\"" << std::endl;
631     ofstrm << "echo \"Generation des microstructures $INDEX_DEPART - $INDEX_FIN\"" << std::endl;
632     ofstrm << "echo \"*********************************************************\"" << std::endl;
633     ofstrm << "STAT=$?" << std::endl;
634     ofstrm << "" << std::endl;
635     ofstrm << "if [ $INDEX_DEPART -eq 1 ]" << std::endl;
636     ofstrm << "then" << std::endl;
637     ofstrm << "mkdir etude" << std::endl;
638     ofstrm << "fi" << std::endl;
639     ofstrm << "cd etude" << std::endl;
640     ofstrm << "for (( k=$INDEX_DEPART; k<=$INDEX_FIN; k++ ))" << std::endl;
641     ofstrm << "do" << std::endl;
642     ofstrm << " echo \"Microstructure #$k\"" << std::endl;
643     ofstrm << " mkdir e$k" << std::endl;
644     ofstrm << " cd e$k" << std::endl;
645     ofstrm << " OPERATION_REUSSIE=0" << std::endl;
646     ofstrm << " until [ $OPERATION_REUSSIE -eq 1 ]; do" << std::endl;
647     ofstrm << " executer ../../ves $DOSSIER_PARAM $CHEMINCOMPILE" << std::endl;
648     ofstrm << " done" << std::endl;
649     ofstrm << " DOSSIER_ACTUEL=$(pwd)" << std::endl;
650     ofstrm << " echo \"$DOSSIER_ACTUEL/ves_analyse_CLDH_sph.ves\" >> ../../lst_ves_analyse_CLDH_sph.txt" << std::endl;
651     ofstrm << " echo \"$DOSSIER_ACTUEL/ves_analyse_CLCH_sph.ves\" >> ../../lst_ves_analyse_CLCH_sph.txt" << std::endl;
652     ofstrm << " echo \"$DOSSIER_ACTUEL/ves_analyse_CLDH_dev.ves\" >> ../../lst_ves_analyse_CLDH_dev.txt" << std::endl;
653     ofstrm << " echo \"$DOSSIER_ACTUEL/ves_analyse_CLCH_dev.ves\" >> ../../lst_ves_analyse_CLCH_dev.txt" << std::endl;
654     ofstrm << " cd .." << std::endl;
655     ofstrm << "done" << std::endl;
656     ofstrm << "cd .." << std::endl;
657     ofstrm << "echo \"**********************\"" << std::endl;
658     ofstrm << "echo \" FIN \"" << std::endl;
659     ofstrm << "echo \"**********************\"" << std::endl;
660     ofstrm << "exit $STAT" << std::endl;
661     ofstrm.close();
662     system((char*)"chmod 744 generer_nb_ves");
663    
664     ofstrm.open((char*)"ves",std::ios::out|std::ios::trunc);
665     ofstrm.precision(16);
666     ofstrm.setf(std::ios::showpoint);
667     ofstrm << "#!/bin/bash" << std::endl;
668     ofstrm << "DOSSIER_PARAM=$1" << std::endl;
669     ofstrm << "CHEMINCOMPILE=$2" << std::endl;
670     ofstrm << "function executer" << std::endl;
671     ofstrm << "{" << std::endl;
672     ofstrm << "$*" << std::endl;
673     ofstrm << "VAR=$?" << std::endl;
674     ofstrm << "if (( $VAR )) ; then " << std::endl;
675     ofstrm << " STATUSTXT=\"ERREUR dans les tests\"; " << std::endl;
676     ofstrm << " STAT=$VAR;" << std::endl;
677     ofstrm << " exit $STAT" << std::endl;
678     ofstrm << "fi" << std::endl;
679     ofstrm << "}" << std::endl;
680     ofstrm << "STATUSTXT=\"Test sans ERREUR\"" << std::endl;
681     ofstrm << "echo \"*******************************\"" << std::endl;
682     ofstrm << "echo \"Generation d'une microstructure\"" << std::endl;
683     ofstrm << "echo \"*******************************\"" << std::endl;
684     ofstrm << "STAT=$?" << std::endl;
685     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_geo -ves_out ves_geo.ves -lstparam $DOSSIER_PARAM/geometrie/lstparam_geo.txt" << std::endl;
686     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_mat -ves_in ves_geo.ves -ves_out ves_mat.ves -lstparam $DOSSIER_PARAM/materiau/lstparam_mat.txt" << std::endl;
687     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_carte -ves_in ves_mat.ves -ves_out ves_carte.ves -param $DOSSIER_PARAM/carte/param_carte.txt" << std::endl;
688     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_mg_maill -ves_in ves_carte.ves -ves_out ves_maill.ves -param $DOSSIER_PARAM/maillage/param_maill.txt -carte carte.ctt" << std::endl;
689     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_fem_maill -ves_in ves_maill.ves -ves_out ves_fem_maill.ves -param $DOSSIER_PARAM/fem_maillage/param_fem_maill.txt" << std::endl;
690     ofstrm << "" << std::endl;
691     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill.ves -ves_out ves_etude_CLDH_sph.ves -param $DOSSIER_PARAM/etude/param_etude_CLDH_sph.txt" << std::endl;
692     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill.ves -ves_out ves_etude_CLCH_sph.ves -param $DOSSIER_PARAM/etude/param_etude_CLCH_sph.txt" << std::endl;
693     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill.ves -ves_out ves_etude_CLDH_dev.ves -param $DOSSIER_PARAM/etude/param_etude_CLDH_dev.txt" << std::endl;
694     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_etude -ves_in ves_fem_maill.ves -ves_out ves_etude_CLCH_dev.ves -param $DOSSIER_PARAM/etude/param_etude_CLCH_dev.txt" << std::endl;
695     ofstrm << "" << std::endl;
696     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -calcul -ves_in ves_etude_CLDH_sph.ves -ves_out ves_calcul_CLDH_sph.ves -param $DOSSIER_PARAM/etude/param_etude_CLDH_sph.txt -paramaster $DOSSIER_PARAM/etude/param_aster.txt" << std::endl;
697     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -calcul -ves_in ves_etude_CLCH_sph.ves -ves_out ves_calcul_CLCH_sph.ves -param $DOSSIER_PARAM/etude/param_etude_CLCH_sph.txt -paramaster $DOSSIER_PARAM/etude/param_aster.txt" << std::endl;
698     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -calcul -ves_in ves_etude_CLDH_dev.ves -ves_out ves_calcul_CLDH_dev.ves -param $DOSSIER_PARAM/etude/param_etude_CLDH_dev.txt -paramaster $DOSSIER_PARAM/etude/param_aster.txt" << std::endl;
699     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -calcul -ves_in ves_etude_CLCH_dev.ves -ves_out ves_calcul_CLCH_dev.ves -param $DOSSIER_PARAM/etude/param_etude_CLCH_dev.txt -paramaster $DOSSIER_PARAM/etude/param_aster.txt" << std::endl;
700     ofstrm << "" << std::endl;
701     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_analyse -ves_in ves_calcul_CLDH_sph.ves -ves_out ves_analyse_CLDH_sph.ves -lstparam $DOSSIER_PARAM/analyse/lstparam_analyse.txt" << std::endl;
702     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_analyse -ves_in ves_calcul_CLCH_sph.ves -ves_out ves_analyse_CLCH_sph.ves -lstparam $DOSSIER_PARAM/analyse/lstparam_analyse.txt" << std::endl;
703     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_analyse -ves_in ves_calcul_CLDH_dev.ves -ves_out ves_analyse_CLDH_dev.ves -lstparam $DOSSIER_PARAM/analyse/lstparam_analyse.txt" << std::endl;
704     ofstrm << "executer $CHEMINCOMPILE/microstructure/microstructure.exe -gen_analyse -ves_in ves_calcul_CLCH_dev.ves -ves_out ves_analyse_CLCH_dev.ves -lstparam $DOSSIER_PARAM/analyse/lstparam_analyse.txt" << std::endl;
705     ofstrm << "" << std::endl;
706     ofstrm << "#executer $CHEMINCOMPILE/microstructure/microstructure.exe -affiche_contenu -ves_in ves_analyse_CLDH_sph.ves" << std::endl;
707     ofstrm << "#executer $CHEMINCOMPILE/microstructure/microstructure.exe -affiche_contenu -ves_in ves_analyse_CLCH_sph.ves" << std::endl;
708     ofstrm << "#executer $CHEMINCOMPILE/microstructure/microstructure.exe -affiche_contenu -ves_in ves_analyse_CLDH_dev.ves" << std::endl;
709     ofstrm << "#executer $CHEMINCOMPILE/microstructure/microstructure.exe -affiche_contenu -ves_in ves_analyse_CLCH_dev.ves" << std::endl;
710 couturad 930 ofstrm << "#find . -type f ! -name '*.ves' -delete" << std::endl;
711 couturad 929 ofstrm << "cd .." << std::endl;
712     ofstrm << "echo \"**********************\"" << std::endl;
713     ofstrm << "echo $STATUSTXT" << std::endl;
714     ofstrm << "echo \"**********************\"" << std::endl;
715     ofstrm << "exit $STAT" << std::endl;
716     ofstrm.close();
717     system((char*)"chmod 744 ves");
718     }
719     }
720    
721    
722    
723 couturad 919 int PARAMETRES::importer_fichier_parametres(char* fichier, OT_PARAMETRES* params)
724     {
725     params->lire(fichier);
726     }
727    
728     int PARAMETRES::importer_fichier_liste_parametres(char* fichier, std::vector< OT_PARAMETRES* >& vector_params)
729     {
730     std::string str_fichier = fichier;
731     std::size_t found = str_fichier.rfind((char*)"/");
732     std::string chemin;
733     if(found!=-1)
734     {
735     chemin = str_fichier.substr(0,found+1);
736     }
737     else chemin = "";
738     FILE* in=fopen(fichier,"rt");
739     char ligne[4000];
740     while(!feof(in))
741     {
742     fgets(ligne,4000,in);
743     char *aide = strstr(ligne,"//");
744     if(aide==NULL)
745     {
746     char message[500];
747     int nb_fichier_param=0;
748     sscanf(ligne,"%s %i ",message,&nb_fichier_param);
749     if(strcmp(message,"NB_FICHIER_PARAM")==0)
750     {
751     for(int i=0;i<nb_fichier_param;i++)
752     {
753     fgets(ligne,4000,in);
754     if(feof(in)) std::cerr << "*** PARAMETRES::PARAMETRES::importer_fichier_liste_parametres : Erreur de lecture ! ***" << std::endl;
755     char sous_fichier[500];
756     sscanf(ligne,"%s",sous_fichier);
757     OT_PARAMETRES *param = new OT_PARAMETRES;
758     std::string str_sous_fichier = chemin;
759     str_sous_fichier.append(sous_fichier);
760     param->lire((char*)str_sous_fichier.c_str());
761     vector_params.push_back(param);
762     }
763     }
764     }
765     }
766     }