ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/REPOS_ERICCA/magic/lib/optimisation/src/mgopt.cpp
Revision: 306
Committed: Wed Jan 11 21:14:52 2012 UTC (13 years, 4 months ago) by francois
File size: 4299 byte(s)
Log Message: 
Adaption de la méthode SIMP à la gestion des matériaux

File Contents

# User Rev Content
1 francois 239 #include "gestionversion.h"
2     #include "mgopt.h"
3     #include "mg_file.h"
4     #include "fct_generateur_frontiere.h"
5     #include "fct_generateur_constante.h"
6     #include "fct_generateur_fichier.h"
7     #include "mailleur_bloc.h"
8 francois 243 #include "mg_export.h"
9 francois 239 #include <string.h>
10    
11    
12    
13    
14    
15    
16    
17    
18    
19    
20    
21    
22    
23     MGOPT::MGOPT():gestd(NULL),gestnd(NULL),carte(NULL)
24     {
25 francois 306 params.ajouter("version_aster",2.,"1. gestion des matériaux par catalogue 2. gestion des matériaux par ccf" );
26 francois 243 params.ajouter("maille",0.,"0. le maillage est effectue ou numero de maillage a utiliser" );
27     params.ajouter("dg",7.,"Ecart nodal constant de la carte de taille");
28     params.ajouter("typecarte",1.,"0. calcul d'une nouvelle carte de taille 1. utilisation d'une carte existante");
29     params.ajouter("fechantillonnage",20.,"Nombre de cellules de la carte de taille dans une direction");
30     params.ajouter("festimation",1.,"nombre d'échantillons dans une direction");
31     params.ajouter("fichiercarte","bielle.ctt","Fichier carte de taille a utiliser quand typecarte vaut 1.");
32     params.ajouter("degre",1.,"1. maillage lineaire 2. maillage quadratique");
33 francois 239 }
34    
35     MGOPT::MGOPT(MGOPT &mdd)
36     {
37     params=mdd.params;
38     }
39    
40    
41     MGOPT::~MGOPT()
42     {
43     if (gestd!=NULL) delete gestd;
44     if (gestd!=NULL) delete gestnd;
45     if (carte!=NULL) delete carte;
46     }
47    
48     void MGOPT::active_affichage(void (*fonc)(char*))
49     {
50     affiche=fonc;
51     affichageactif=1;
52     }
53     void MGOPT::lire_params(char *fichier)
54     {
55     params.vide();
56     params.lire(fichier);
57     }
58    
59     void MGOPT::optimise(char *fichier)
60     {
61 francois 272 affiche((char*)"");
62     affiche((char*)"*************************");
63     affiche((char*)"Optimisation de topologie");
64     affiche((char*)"*************************");
65     affiche((char*)"");
66     affiche((char*)"");
67     affiche((char*)"");
68     affiche((char*)"Écriture d'un fichier de parametres par défaut");
69 francois 239 params.enregistrer(fichier);
70     }
71    
72     void MGOPT::optimise(char *nomgestd,char *nomgestnd,char *nomsortie,char *nomparam)
73     {
74 francois 272 affiche((char*)"");
75     affiche((char*)"*************************");
76     affiche((char*)"Optimisation de topologie");
77     affiche((char*)"*************************");
78     affiche((char*)"");
79     affiche((char*)"");
80     affiche((char*)"");
81 francois 271 time_t heuredeb = time(NULL);
82     tdebut = *localtime(&heuredeb);
83 francois 239 if (nomparam!=NULL) lire_params(nomparam);
84     char *p=strchr(nomgestd,'.');
85     strncpy(nometude,nomgestd,p-nomgestd);
86     nometude[p-nomgestd]=0;
87 francois 243 p=strchr(nomsortie,'.');
88     strncpy(nometudesortie,nomsortie,p-nomsortie);
89     nometudesortie[p-nomsortie]=0;
90 francois 272 affiche((char*)"Preparation du modèle");
91 francois 239 gestd=new MG_FILE(nomgestd);
92     gestnd=new MG_FILE(nomgestnd);
93     unsigned long maille=(unsigned long)params.get_valeur("maille");
94     FEM_MAILLAGE* fem;
95     if (maille==0) fem=prepare_model(gestd,gestnd);
96     else fem=gestd->get_fem_maillageid(maille);
97 francois 272 affiche((char*)"Debut optimisation");
98 francois 239 optimisation(fem);
99 francois 272 affiche((char*)"Enregistrement resultat");
100 francois 239 gestd->enregistrer(nomsortie);
101 francois 272 affiche((char*)"Enregistrement resultat sous GMSH");
102 francois 243 MG_EXPORT exp;
103     char nomfichier[500];
104     sprintf(nomfichier,"%s_mg",nometudesortie);
105     exp.gmsh(fem->get_mg_maillage(),nomfichier);
106     sprintf(nomfichier,"%s_fem",nometudesortie);
107     exp.gmsh(fem,nomfichier);
108 francois 239 }
109    
110     FEM_MAILLAGE* MGOPT::prepare_model(MG_GESTIONNAIRE *gestd,MG_GESTIONNAIRE *gestnd)
111     {
112 francois 272 affiche((char*)"Carte de taille");
113 francois 239 int typecarte=(int)params.get_valeur("typecarte");
114     if (typecarte==0)
115     {
116 francois 272 affiche((char*)" Calcul");
117 francois 239 double eng=params.get_valeur("dg");
118     int fechan=(int)params.get_valeur("fechantillonnage");
119     int festim=(int)params.get_valeur("festimation");
120     carte=new FCT_GENERATEUR_CONSTANTE(*gestd,eng);
121     FCT_GENERATEUR_CONSTANTE *carte2=(FCT_GENERATEUR_CONSTANTE*)carte;
122     carte2->construit(fechan,festim);
123     char nomfichiersortie[300];
124     sprintf(nomfichiersortie,"%s.ctt",nometude);
125     carte2->enregistrer(nomfichiersortie);
126     }
127     else
128     {
129 francois 272 affiche((char*)" Lecture");
130 francois 239 carte=new FCT_GENERATEUR_3D<4>;
131     std::string fichiercarte=params.get_nom("fichiercarte");
132     carte->lire((char *)fichiercarte.c_str());
133     }
134 francois 272 affiche((char*)"Maillage par le mailleur bloc");
135 francois 239 FCT_TAILLE* metrique=carte;
136     MAILLEUR_BLOC m(gestd,gestnd,0,0,metrique);
137     m.active_affichage(affiche);
138     m.maille();
139 francois 272 affiche((char*)"Creation du maillage FEM");
140 francois 239 MG_MAILLAGE* mai=gestd->get_mg_maillage(gestd->get_nb_mg_maillage()-1);
141     int degre=params.get_valeur("degre");
142     FEM_MAILLAGE* fem=new FEM_MAILLAGE(mai->get_mg_geometrie(),mai,degre);
143     gestd->ajouter_fem_maillage(fem);
144     fem->construire(0);
145     return fem;
146     }