1 |
francois |
737 |
/*! \page Aide Aide MAGiC |
2 |
francois |
868 |
Informations diverses |
3 |
|
|
|
4 |
|
|
- \subpage infodiv |
5 |
|
|
|
6 |
francois |
737 |
Configuration de MAGiC |
7 |
|
|
|
8 |
|
|
- \subpage Configuration |
9 |
|
|
|
10 |
|
|
|
11 |
|
|
Description du format du fichier script MAGiC\n |
12 |
|
|
- \subpage Script |
13 |
francois |
864 |
|
14 |
francois |
902 |
Description de la procédure de reconstruction d'un modèle structurel après optimisation topologique\n |
15 |
francois |
864 |
- \subpage Reconstruction |
16 |
francois |
902 |
|
17 |
francois |
1038 |
Faire une étude SiDolo |
18 |
|
|
|
19 |
|
|
- \subpage SiDolo |
20 |
|
|
|
21 |
|
|
|
22 |
francois |
737 |
*/ |
23 |
francois |
868 |
//----------------------------------------------------------- |
24 |
francois |
737 |
|
25 |
francois |
868 |
|
26 |
|
|
/*! \page infodiv Informations diverses sur le projet MAGiC |
27 |
|
|
Le projet MAGiC est un projet des professeurs Jean-Christophe Cuillière et Vincent Francois de l'Équipe de Recherche en Intégration CAO-CAlcul du département de génie mécanique de l'UQTR.\n\n |
28 |
|
|
Cette page résume les différents liens pour les outils disponibles : \n\n |
29 |
|
|
<ul> |
30 |
francois |
1040 |
<li><a href="http://ericca.uqtr.ca/cgi-bin/viewvc.cgi"> Web SVN </a></li> |
31 |
francois |
868 |
<li><a href="http://ericca.uqtr.ca/biblio"> Biblio ERICCA </a></li> |
32 |
|
|
<li><a href="http://ericca.uqtr.ca/cas_test"> Cas test d'optimisation de topologie </a></li> |
33 |
francois |
869 |
<li><a href="http://ericca.uqtr.ca/webmagic/index_text.html"> Tests journaliers du code </a></li> |
34 |
|
|
|
35 |
francois |
868 |
<li>Documentation de Code ASTER</li> |
36 |
|
|
<ul> |
37 |
|
|
<li><a href="http://ericca.uqtr.ca/fr11.7"> 11.7 </a></li> |
38 |
|
|
<li><a href="http://ericca.uqtr.ca/fr12.4"> 12.4 </a></li> |
39 |
francois |
963 |
<li><a href="http://ericca.uqtr.ca/fr13.6"> 13.6 </a></li> |
40 |
francois |
868 |
</ul> |
41 |
sattarpa |
901 |
<li>Modification dans Code ASTER</li> |
42 |
francois |
868 |
<ul> |
43 |
sattarpa |
901 |
<li>\subpage ExtraireK </a></li> |
44 |
|
|
<li>\subpage CalcParal </a></li> |
45 |
|
|
</ul> |
46 |
|
|
|
47 |
|
|
<li>Script d'installation de l'environnement ERICCA : Downloader les deux fichiers suivants et faire ./instal_magic.bash Version_Aster [tout]. L'option installe les outils de développements. Il faut s'assurer que les scripts soient exécutables.</li> |
48 |
|
|
<ul> |
49 |
francois |
868 |
<li><a href="http://ericca.uqtr.ca/linux/instal_magic.bash"> Script d'installation de MAGiC </a></li> |
50 |
|
|
<li><a href="http://ericca.uqtr.ca/linux/instal_aster.bash"> Script d'installation de code ASTER</a></li> |
51 |
francois |
955 |
<li><a href="http://ericca.uqtr.ca/linux/instal_chrono.bash"> Script d'installation de Project Chrono</a></li> |
52 |
francois |
868 |
|
53 |
|
|
</ul> |
54 |
|
|
</ul> |
55 |
sattarpa |
901 |
*/ |
56 |
francois |
868 |
|
57 |
sattarpa |
901 |
//----------------------------------------------------------- |
58 |
francois |
868 |
|
59 |
sattarpa |
901 |
/*! \page ExtraireK Extraire matrice de raideur |
60 |
|
|
Pour extraire le matrice de rigidité lors d'un calcul de Code Aster dans un fichier en format MATLAB, Il faut utiliser une routine Fortran pour cela (matimp.f) que l'on vient appeler au bon endroit (par exemple dans "op0012.F90"). \n |
61 |
|
|
Pour la version 12.4 de Code Aster, il faut faire la modification suivant dans "op0012.F90". On le trouve dans le paquet d'installation de Code Aster à cet endroit: "aster-full-src-12.4.0/SRC/aster-12.4.0/bibfor/op/op0012.F90". \n |
62 |
|
|
On édite "op0012.F90" on ajoutant "call matimp(matas,29,'MATLAB')" à la lieu montrée. \n |
63 |
francois |
868 |
|
64 |
sattarpa |
901 |
\verbatim |
65 |
|
|
! -- SI MATAS N'EST PAS MPI_COMPLET, ON LA COMPLETE : |
66 |
|
|
call dismoi('MPI_COMPLET', matas, 'MATR_ASSE', repk=kmpic) |
67 |
|
|
ASSERT((kmpic.eq.'OUI').or.(kmpic.eq.'NON')) |
68 |
|
|
if (kmpic .eq. 'NON') call sdmpic('MATR_ASSE', matas) |
69 |
|
|
! |
70 |
|
|
! |
71 |
|
|
! |
72 |
|
|
call matimp(matas,29,'MATLAB') << On appel matimp ICI << |
73 |
|
|
! |
74 |
|
|
! -- MENAGE : |
75 |
|
|
call jedetr(lchci) |
76 |
|
|
call jedetr(lmatel) |
77 |
|
|
! |
78 |
|
|
call jedema() |
79 |
|
|
end subroutine |
80 |
|
|
\endverbatim \n |
81 |
|
|
Et puis, on compile Code ASTER: sudo python setup.py install\n |
82 |
|
|
Pour etre capable de sauvegarder le matrice de rigidité dans le fichier MATLAB, il faut editer le fichier .export de Code ASTER comme montré ci dessous: \n |
83 |
|
|
\verbatim |
84 |
|
|
F comm /home/..../filename.comm D 1 |
85 |
|
|
F mail /home/..../filename.mail D 20 |
86 |
|
|
F mess /home/..../filename.mess R 6 |
87 |
|
|
F m /home/..../filename.m R 29 << On ajoute cette ligne ICI << |
88 |
|
|
\endverbatim \n |
89 |
|
|
*/ |
90 |
francois |
868 |
|
91 |
sattarpa |
901 |
//----------------------------------------------------------- |
92 |
francois |
868 |
|
93 |
sattarpa |
901 |
/*! \page CalcParal Calcul parallèle dans Code ASTER |
94 |
|
|
Référence: <a href="https://sites.google.com/site/codeastersalomemeca/home/code_asterno-heiretuka/parallel-code_aster-12-4-english"> Parallel Code_Aster 12.4 </a>\n\n |
95 |
|
|
<B>Version</B>\n |
96 |
|
|
OS:Ubuntu 14.04 \n |
97 |
|
|
Code_Aster : ver.12.4\n |
98 |
|
|
-----------------------------------------------------------\n\n |
99 |
|
|
<B>Telecharger des fichiers</B>\n |
100 |
|
|
Les fichier sont enregistrés dans: ~/Install_Files\n |
101 |
|
|
Ils sont installé dans : /opt and /opt/aster\n\n |
102 |
|
|
Télécharger les fichier suivant:\n |
103 |
|
|
aster-full-src-12.4.0-1.noarch.tar.gz (Code_Aster)\n |
104 |
|
|
OpenBLAS-0.2.15.tar.gz (OpenBLAS)\n |
105 |
|
|
scalapack_installer.tgz (ScaLAPACK)\n |
106 |
|
|
petsc-3.4.5.tar.gz (PETSc)\n\n |
107 |
|
|
Modifier le propriétaire d'installation par: $ sudo chown username /opt/\n |
108 |
|
|
Installez les bibliothèques pour Code_Aster par: \n |
109 |
|
|
$ sudo apt-get install gfortran g++ python-dev python-qt4 python-numpy liblapack-dev libblas-dev tcl tk zlib1g-dev bison flex checkinstall openmpi-bin libopenmpi-dev libx11-dev cmake qt4-dev-tools libmotif-dev \n |
110 |
|
|
-----------------------------------------------------------\n\n |
111 |
|
|
<B>Compiler OpenBLAS (Math LIbrary pour Code_Aster)</B>\n |
112 |
|
|
$ cd ~/Install_Files/\n |
113 |
|
|
$ tar xfvz OpenBLAS-0.2.15..tar.gz\n |
114 |
|
|
$ cd OpenBLAS-0.2.15\n |
115 |
|
|
$ make NO_AFFINITY=1 USE_OPENMP=1\n |
116 |
|
|
$ make PREFIX=/opt/OpenBLAS install\n |
117 |
|
|
$ echo /opt/OpenBLAS/lib | sudo tee -a /etc/ld.so.conf.d/openblas.conf\n |
118 |
|
|
$ sudo ldconfig\n |
119 |
|
|
-----------------------------------------------------------\n\n |
120 |
|
|
<B>Compiler Code_Aster (séquentiel) avec OpenBLAS.</B>\n |
121 |
|
|
$ cd ~/Install_Files\n |
122 |
|
|
$ tar xfvz aster-full-src-12.4.0-1.noarch.tar.gz\n |
123 |
|
|
$ cd aster-full-src-12.4.0/\n |
124 |
|
|
$ sed -i "s:PREFER_COMPILER\ =\ 'GNU':PREFER_COMPILER\ =\'GNU_without_MATH'\nMATHLIB=\ '/opt/OpenBLAS/lib/libopenblas.a':g" setup.cfg\n |
125 |
|
|
$ python setup.py install\n |
126 |
|
|
<EM>Faire host-file pour calcul parallèle par:</EM>\n |
127 |
|
|
$ echo "$HOSTNAME cpu=$(cat /proc/cpuinfo | grep processor | wc -l)" > /opt/aster/etc/codeaster/mpi_hostfile\n |
128 |
|
|
-----------------------------------------------------------\n\n |
129 |
|
|
<B>Compiler ScaLAPACK</B>\n |
130 |
|
|
$ cd ~/Install_Files\n |
131 |
|
|
$ tar xfvz scalapack_installer.tgz\n |
132 |
|
|
$ cd scalapack_installer_1.0.2\n |
133 |
|
|
$ ./setup.py --lapacklib=/opt/OpenBLAS/lib/libopenblas.a --mpicc=mpicc --mpif90=mpif90 --mpiincdir=/usr/lib/openmpi/include --ldflags_c=-fopenmp --ldflags_fc=-fopenmp --prefix=/opt/scalapack\n |
134 |
|
|
<EM>Une message d'erreur "BLACS: error running BLACS test routines xCbtest" apparaîtra après la compilation, mais vous réussissez, s'il existe le fichier dans "/opt/scalapack/lib/libscalapack.a"</EM> \n |
135 |
|
|
-----------------------------------------------------------\n\n |
136 |
|
|
<B>Compiler MUMPS</B>\n |
137 |
|
|
Copier mumps-4.10.0 dans fichiers source de Code_Aster dans '/opt.\n |
138 |
|
|
Compiler par mpi-compiler\n |
139 |
|
|
$ cp ~/Install_Files/aster-full-src-12.4.0/SRC/mumps-4.10.0-aster3.tar.gz /opt/\n |
140 |
|
|
$ cd /opt\n |
141 |
|
|
$ tar xfvz mumps-4.10.0-aster3.tar.gz\n |
142 |
|
|
$ mv mumps-4.10.0 mumps-4.10.0_mpi\n |
143 |
|
|
$ cd mumps-4.10.0_mpi/\n |
144 |
|
|
Changer 'Makefile.inc.in' pour MUMPS4.10.0 Makefile.inc\n |
145 |
|
|
$ make all \n |
146 |
|
|
-----------------------------------------------------------\n\n |
147 |
|
|
<B>Compiler PETSc avec HYPRE et ML</B>\n |
148 |
|
|
$ cp ~/Install_Files/petsc-3.4.5.tar.gz /opt\n |
149 |
|
|
$ cd /opt\n |
150 |
|
|
$ tar xfvz petsc-3.4.5.tar.gz\n |
151 |
|
|
$ cd petsc-3.4.5\n |
152 |
|
|
$ ./config/configure.py --with-mpi-dir=/usr/lib/openmpi --with-blas-lapack-lib=/opt/OpenBLAS/lib/libopenblas.a --download-hypre=yes --download-ml=yes --with-debugging=0 COPTFLAGS=-O1 CXXOPTFLAGS=-O1 FOPTFLAGS=-O1 --configModules=PETSc.Configure --optionsModule=PETSc.compilerOptions --with-x=0 --with-shared-libraries=0 \n |
153 |
|
|
$ make PETSC_DIR=/opt/petsc-3.4.5 PETSC_ARCH=arch-linux2-c-opt all\n |
154 |
|
|
$ make PETSC_DIR=/opt/petsc-3.4.5 PETSC_ARCH=arch-linux2-c-opt test\n |
155 |
|
|
-----------------------------------------------------------\n\n |
156 |
|
|
<B>Compiler Code_Aster (parallèle)</B>\n |
157 |
|
|
<EM>Changer un partie de 'mpi_get_procid_cmd' de '/opt/aster/etc/codeaster/asrun' par:\n |
158 |
|
|
mpi_get_procid_cmd : echo $OMPI_COMM_WORLD_RANK\n |
159 |
|
|
<EM>Changer 32 a 64 pour les processeur 64 bit: </B>\n |
160 |
|
|
batch_mpi_nbpmax : 32 (ou 64)\n |
161 |
|
|
interactif_mpi_nbpmax : 32 (ou 64)\n |
162 |
|
|
<EM>Décompresser Code_Aster</B>\n |
163 |
|
|
$ cd ~/Install_Files\n |
164 |
|
|
$ cd aster-full-src-12.4.0/SRC\n |
165 |
|
|
$ tar xfvz aster-12.4.0.tgz\n |
166 |
|
|
$ cd aster-12.4.0\n |
167 |
|
|
<EM>Mettre les fichiers de configuration pour calcul parallèle: <a href="https://sites.google.com/site/codeastersalomemeca/home/code_asterno-heiretuka/code_asterno-heiretuka-12-4/Ubuntu_gnu_mpi.py?attredirects=0&d=1"> Ubuntu_gnu_mpi.py </a> et <a href="https://sites.google.com/site/codeastersalomemeca/home/code_asterno-heiretuka/code_asterno-heiretuka-12-4/Ubuntu_gnu.py?attredirects=0&d=1"> Ubuntu_gnu.py </a> dans cette dossier et compiler Code_Aster.</B>\n |
168 |
|
|
$ export ASTER_ROOT=/opt/aster\n |
169 |
|
|
$ ./waf configure --use-config-dir=$ASTER_ROOT/12.4/share/aster --use-config=Ubuntu_gnu_mpi --prefix=$ASTER_ROOT/PAR12.4\n |
170 |
|
|
$ ./waf install -p\n |
171 |
|
|
<EM>Ajouter 'vers : PAR12.4:/opt/aster/PAR12.4/share/aster' en dessous de 'vers : testing' in '/opt/aster/etc/codeaster/aster' et puis 'PAR12.4' est enregistré dans ASTK.</B>\n |
172 |
|
|
|
173 |
francois |
868 |
*/ |
174 |
|
|
|
175 |
francois |
737 |
//----------------------------------------------------------- |
176 |
francois |
868 |
//----------------------------------------------------------- |
177 |
francois |
737 |
|
178 |
|
|
|
179 |
|
|
/*! \page Script Description du contenu du fichier script |
180 |
|
|
Un fichier script est composé de commandes exécutées en séquences depuis le début du fichier jusqu'à la fin.\n |
181 |
|
|
Toute ligne commence par une variable suivie d'un espace suivi d'un signe egal suivi d'un espace et suivi de la commande et de ses arguments.\n |
182 |
|
|
Si le premier caractère de la ligne est le caractère # alors la ligne est ignorée.\n\n |
183 |
|
|
Les commandes disponibles dans le script MAGiC: |
184 |
couturad |
740 |
- \subpage ajoutass |
185 |
couturad |
906 |
- \subpage brepassemblage |
186 |
|
|
- \subpage creemodele |
187 |
couturad |
740 |
- \subpage creeassemblage |
188 |
francois |
737 |
- \subpage diff |
189 |
|
|
- \subpage enregistrer |
190 |
|
|
- \subpage export_mg_gmsh |
191 |
|
|
- \subpage export_fem_gmsh |
192 |
|
|
- \subpage fichier |
193 |
|
|
- \subpage importstepocc |
194 |
|
|
- \subpage inter |
195 |
|
|
- \subpage p_cylindre |
196 |
|
|
- \subpage p_cone |
197 |
|
|
- \subpage p_boite |
198 |
|
|
- \subpage p_sphere |
199 |
couturad |
906 |
- \subpage p_ellipsoide |
200 |
francois |
737 |
- \subpage p_tore |
201 |
|
|
- \subpage union |
202 |
couturad |
906 |
- \subpage fragment |
203 |
|
|
- \subpage tristlmodele |
204 |
|
|
- \subpage tristlassemblage |
205 |
francois |
737 |
*/ |
206 |
|
|
|
207 |
|
|
//----------------------------------------------------------- |
208 |
|
|
|
209 |
|
|
/*! \page Configuration Configuration de MAGiC |
210 |
|
|
Dans chaque compte utilisateur ou MAGiC est utilisé un fichier .magic dans le home directory permet de regler quelques parametres de l'application :\n\n |
211 |
|
|
Affiche_Code_Aster = 0.000000 //0. Affichage redirigé vers le fichier aster.log 1 affichage dans un terminal\n |
212 |
|
|
Convergence_mailleur3d_frontal = 200.000000 //Nombre de front restant à partir duquel la destruction est obligatoire\n |
213 |
|
|
Echantillon_face = 20.000000 //Nombre d'échantillon pour le calcul de la boite englobante des faces\n |
214 |
|
|
Verification_coin_face_mailleur2D = 0.000000 //0. Ne fais pas la vérification 1. Fait la vérification qu'un sommet de face partage au moins deux triangles\n |
215 |
|
|
Optimisation_2D = 1 //Nombre de bits = nombre de passe d'optimisation. Pour chaque bit 1=bouge de noeud 2=inverse diagonale 3=bouge de noeud inserse diagonale simultane 4=inverse diagonale bouge de noeud simulatane\n |
216 |
francois |
883 |
Quadratisation_jmin = 1.000000 //Garantir des tetras quadratiques avec jaconien minimum positif : 1=oui 2=non |
217 |
|
|
Quadratisation_pas = 100. //Nombre de pas pour recaler les noeuds et obtenir un jmin positif |
218 |
|
|
Quadratisation_dis = 0.100000 //Valeur de la distortion du jacobien minimale à atteindre |
219 |
|
|
Separateur_decimale = , //Seprateur décimale pour l'exportation des fichiers excel |
220 |
francois |
737 |
*/ |
221 |
|
|
|
222 |
francois |
864 |
|
223 |
francois |
737 |
//----------------------------------------------------------- |
224 |
francois |
864 |
|
225 |
|
|
/*! \page Reconstruction Description de la procédure de reconstruction d'un nmodèle structurel après optimisation topologique |
226 |
|
|
|
227 |
|
|
Description de la reconstruction |
228 |
|
|
|
229 |
|
|
- \subpage Description_reconstruction |
230 |
|
|
|
231 |
|
|
Commandes de la reconstruction |
232 |
|
|
|
233 |
|
|
- \subpage Commande_reconstruction |
234 |
|
|
|
235 |
|
|
*/ |
236 |
|
|
|
237 |
|
|
//----------------------------------------------------------- |
238 |
|
|
/*! \page Description_reconstruction Description de la reconstruction |
239 |
|
|
La squelettisation consiste à extraire le squelette d’un modèle. Le modèle est supposé structurel, c’est-à-dire contenant des éléments de type poutre. Pour utiliser l’algorithme de squelettisation, l’on a besoin de 2 fichiers .magic et 2 fichiers .txt. Il s’agit du fichier résultant de l’optimisation topologique, du fichier obtenu après lissage (partie de non_design exclue), et les 2 fichiers .txt provenant de MatLab et contenant les coordonnées des points du squelette et les connectivités entre les points de la courbe-squelette. |
240 |
|
|
L’étape 5 ci-dessus résulte en 6 fichiers qui sont : 3 fichiers de squelette (en .magic) et 3 fichiers de solutions (en .sol) sur le calcul des sections.\n\n\n |
241 |
|
|
<B>nomfichier_squelette1.magic</B> est le fichier contenant la courbe-squelette sous forme de ligne. Il est accompagné du fichier Carte_section1.sol qui est le fichier de solution pour la taille des sections en mètre (m).\n\n\n |
242 |
|
|
<B>nomfichier_squelette2.magic</B> est le fichier contenant le squelette après transformation des segments de chaque branche de la courbe squelette en poutre droite. Il est lui aussi accompagné du fichier Carte_section2.sol qui est le fichier de solution pour la taille des sections des poutres droites en mètre (m). Ce squelette est connecté au non_design à l’aide des segments.\n\n\n |
243 |
|
|
<B>nomfichier_squelette3.magic</B> le fichier final de travail. C’est ce fichier qui est utilisé dans la suite pour les calculs. Ce dernier fichier contient le squelette sous forme de poutres droites, est connecté directement au non_design à l’aide des mini-poutres et contient les conditions aux limites et de chargement. Il est accompagné du fichier Carte_section3.sol qui est le fichier de solution pour la taille des sections en mètre (m). |
244 |
|
|
*/ |
245 |
|
|
|
246 |
|
|
//----------------------------------------------------------- |
247 |
|
|
/*! \page Commande_reconstruction Commandes de la reconstruction |
248 |
|
|
|
249 |
|
|
<span style="text-decoration: underline;"><B>Etape 1 :</B> </span> Optimisation topologique\n <i>topo_optis.exe |
250 |
|
|
-optimise -design nom_fichier.magic -nondesign nom_fichier_bloc.magic -out nom_fichier_resu.magic -param nom_fichier_param.txt -paramaster aster.txt</i>\n\n |
251 |
|
|
Si on a déjà un résultat optimal on peut directement passer à l'étape 2. |
252 |
|
|
<BR><BR> |
253 |
|
|
|
254 |
|
|
<span style="text-decoration: underline;"><B>Etape 2 :</B> </span> Lissage du résultat optimal (sans la partie de non_design) \n |
255 |
|
|
<i>post_optis.exe -creepeau -in nom_fichier_resu.magic -out nom_fichier_lisse.magic -param paramliss.txt</i>\n\n |
256 |
|
|
Dans le fichier paramliss, mettre à 1 le paramètre rmimpose_debut pour supprimer le non_design et donner un nom de fichier pour le résultat de la peau sans non_design avec le paramètre nomfichpeau. |
257 |
|
|
<BR><BR> |
258 |
|
|
|
259 |
|
|
|
260 |
|
|
<span style="text-decoration: underline;"><B>Etape 3</B> </span> : Transformation en .off \n |
261 |
|
|
<i>mgoperation.exe -offout -in nom_fichier_lisse_peau.magic</i>\n\n |
262 |
|
|
(Utiliser le fichier de la peau qui ne contient pas le non_design) |
263 |
|
|
<BR><BR> |
264 |
|
|
|
265 |
|
|
|
266 |
|
|
<span style="text-decoration: underline;"><B>Etape 4 </B> </span> : Squelettisation sous MatLab\n |
267 |
|
|
|
268 |
|
|
Partie 1 : génération du squelette\n |
269 |
|
|
Mettre le fichier .off précédent dans le dossier « data » et aller dans le logiciel MatLab. Choisir le chemin vers le dossier « matlab » qui est dans le même répertoire que « data »\n |
270 |
|
|
Ouvrir eg_skeleton_laplacian_rosa.m\n |
271 |
|
|
Mettre le nom (sans son extension) du fichier .off obtenu à l’Étape 3 dans eg_skeleton_laplacian_rosa.m à la ligne 16.\n |
272 |
|
|
executer eg_skeleton_laplacian_rosa.m\n |
273 |
|
|
Des figures (en particulier la figure 4) sont générées.\n\n |
274 |
|
|
|
275 |
|
|
Partie 2 : extraction des données\n |
276 |
|
|
Ouvrir maintenant Entxt.m qui est également dans le dossier « matlab »\n |
277 |
|
|
Donner les noms de fichier désirés dans Entxt.m aux lignes 4, 5 et 33, 34.\n |
278 |
|
|
executer Entxt.m \n |
279 |
|
|
2 fichiers, contenant les coordonnées des points et les connectivités, sont générés dans le dossier « matlab » ; ce sont ces fichiers qui sont utilisés par la suite dans MAGiC.\n |
280 |
|
|
<BR><BR> |
281 |
|
|
|
282 |
|
|
|
283 |
|
|
<span style="text-decoration: underline;"><B>Etape 5 :</B> </span> Exportation du squelette dans MAGiC\n |
284 |
|
|
<i>skeleton.exe -construire_squelette -in nom_fichier_resu.magic -out nom_fichier -inpoints nom_fichier_points.txt -incorres nom_fichier_corres.txt -infilelisse nom_fichier_lisse.magic –connecte 1 –param_lissage paramliss.txt</i>\n\n |
285 |
|
|
|
286 |
|
|
|
287 |
|
|
Pour obtenir une liste des paramètres pour d’éventuelles modifications, faire simplement :\n |
288 |
|
|
<i>skeleton.exe -construire_squelette </i>\n |
289 |
|
|
|
290 |
|
|
–connecte 1 : permet d’introduire des mini-poutres dans toute la zone (les triangles) de contact entre les tétraèdres optimisés et de non_design\n |
291 |
|
|
|
292 |
|
|
–connecte 0 : permet d’introduire des mini-poutres uniquement au nœud de jonction entre le squelette et le non_design\n\n |
293 |
|
|
|
294 |
|
|
Le résultat à utiliser pour la suite est le nom_fichier_squelette3.magic |
295 |
|
|
<BR><BR> |
296 |
|
|
|
297 |
|
|
<span style="text-decoration: underline;"><B>Etape 6 :</B> </span> Lancer le calcul EF \n |
298 |
|
|
<i>mgoperation.exe -calculaster -in nom_fichier_squelette3.magic -out nom_fichier_sortie -param aster.txt -elastique</i>\n |
299 |
|
|
*/ |
300 |
francois |
1038 |
//----------------------------------------------------------- |
301 |
|
|
/*! \page SiDolo Utilisation de SiDolo couplé à MAGiC |
302 |
|
|
<ul> |
303 |
|
|
<li>Avant étude</li> |
304 |
|
|
<ul> |
305 |
|
|
<li>S'assurer que SiDolo est disponible sur son ordinateur. Si ce n'est pas le cas il est diponible <a href="http://ericca.uqtr.ca/SiDolo"> ici </a>. Mettre à jour le fichier sidolo_user dans le répertoire bin afin de mettre le bon répertoire.</li> |
306 |
|
|
</ul> |
307 |
|
|
<li>Étape 1</li> |
308 |
|
|
<ul> |
309 |
|
|
<li>Créer un répertoire pour une étude d'optimisation. Conseil pratique un répetertoire = une étude</li> |
310 |
|
|
<li>Entrer dans ce répertoire</li> |
311 |
|
|
</ul> |
312 |
|
|
<li>Étape 2</li> |
313 |
|
|
<ul> |
314 |
|
|
<li>Créer un modèle MAGiC.</li> |
315 |
|
|
<li>Créer un script qui fait un calcul d'une itération d'optimisation.</li> |
316 |
|
|
</ul> |
317 |
|
|
<li>Étape 3</li> |
318 |
|
|
<ul> |
319 |
|
|
<li>Créer l'ensemble des fichiers nécessaire à une étude Sidolo. <a href="http://ericca.uqtr.ca/SiDolo"> Exemple de fichiers </a>.</li> |
320 |
|
|
<ul> |
321 |
|
|
<li>fichier .dat. Définition des parametre de l'étude </li> |
322 |
|
|
<li>fichier sidvar.h. Définition des variables </li> |
323 |
|
|
<li>fichier interface.f. Programmation des mises à jour de données entre chaque itération</li> |
324 |
|
|
<li>fichier .exp. Fichier de la courbe expérimentale à approcher </li> |
325 |
|
|
</ul> |
326 |
|
|
<li>Configurer les fichiers en se référant à la documentation SiDolo présente dans le répertoire doc de SiDolo.</li> |
327 |
|
|
</ul> |
328 |
|
|
<li>Étape 4</li> |
329 |
|
|
<ul> |
330 |
|
|
<li><span style="color:magenta">source ../ou est SiDolo/bin/sidolo_user </span>. </li> |
331 |
|
|
<li><span style="color:magenta">sde nomoptim</span> pour creer l'excutable de l'optimisation. </li> |
332 |
|
|
<li><span style="color:magenta">sdi interface</span> pour creer l'executable de l'interface. </li> |
333 |
|
|
<li><span style="color:magenta">ln -s "executable de l'interface" "fichier de la courbe expérimentale .exp"</span></li> |
334 |
|
|
</ul> |
335 |
|
|
<li>Étape 5</li>> |
336 |
|
|
<ul> |
337 |
|
|
<li>Lancer l'executable de l'optimisateur</li> |
338 |
|
|
</ul> |
339 |
|
|
</ul> |
340 |
francois |
864 |
|
341 |
francois |
1038 |
*/ |
342 |
|
|
|
343 |
|
|
|
344 |
francois |
864 |
//----------------------------------------------------------- |
345 |
francois |
737 |
/*! \page fichier fichier |
346 |
|
|
var = fichier nom |
347 |
|
|
@brief Mise en mémoire d'un gestionnaire MAGiC |
348 |
couturad |
906 |
@param nom une chaine de caractére. |
349 |
francois |
737 |
@return var une variable representant le gestionnaire MAGiC |
350 |
|
|
@details Si nom = nouveau il y a création d'un nouveau modèle MAGiC\n |
351 |
|
|
@details Sinon ouverture du fichier MAGiC donné par nom\n |
352 |
|
|
*/ |
353 |
|
|
|
354 |
|
|
//----------------------------------------------------------- |
355 |
|
|
|
356 |
|
|
/*! \page p_boite p_boite |
357 |
couturad |
906 |
var = p_boite modele x1 y1 z1 x2 y2 z2 |
358 |
|
|
@brief Création d'une forme primitive volumique de type boîte [MG_CG_FORME_VOLUME_BOITE] entre les points (x1,y1,z1) et (x2,y2,z2) |
359 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
360 |
francois |
737 |
@param x1 un double. |
361 |
|
|
@param y1 un double. |
362 |
|
|
@param z1 un double. |
363 |
|
|
@param x2 un double. |
364 |
|
|
@param y2 un double. |
365 |
|
|
@param z2 un double. |
366 |
couturad |
906 |
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
367 |
francois |
737 |
*/ |
368 |
|
|
|
369 |
|
|
//----------------------------------------------------------- |
370 |
|
|
|
371 |
|
|
/*! \page p_sphere p_sphere |
372 |
couturad |
906 |
var = p_sphere modele x1 y1 z1 r |
373 |
|
|
@brief Création d'une forme primitive volumique de type sphère [MG_CG_FORME_VOLUME_SPHERE] de centre (x1,y1,z1) et de rayon r |
374 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
375 |
francois |
737 |
@param x1 un double. |
376 |
|
|
@param y1 un double. |
377 |
|
|
@param z1 un double. |
378 |
|
|
@param r un double. |
379 |
couturad |
906 |
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
380 |
francois |
737 |
*/ |
381 |
|
|
|
382 |
|
|
//----------------------------------------------------------- |
383 |
|
|
|
384 |
|
|
/*! \page p_cylindre p_cylindre |
385 |
couturad |
906 |
var = p_cylindre modele x1 y1 z1 axe_x axe_y axe_z h r |
386 |
|
|
@brief Création d'une forme primitive volumique de type cylindre [MG_CG_FORME_VOLUME_CYLINDRE] |
387 |
|
|
basée sur le cercle centré en (x1,y1,z1) orienté selon l'axe (axe_x,axe_y,axe_z) de hauteur h et de rayon r |
388 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
389 |
francois |
737 |
@param x1 un double. |
390 |
|
|
@param y1 un double. |
391 |
|
|
@param z1 un double. |
392 |
couturad |
906 |
@param axe_x un double. |
393 |
|
|
@param axe_y un double. |
394 |
|
|
@param axe_z un double. |
395 |
francois |
737 |
@param h un double. |
396 |
|
|
@param r un double. |
397 |
couturad |
906 |
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
398 |
|
|
*/ |
399 |
francois |
737 |
|
400 |
couturad |
906 |
//----------------------------------------------------------- |
401 |
francois |
902 |
|
402 |
couturad |
906 |
/*! \page p_ellipsoide p_ellipsoide |
403 |
|
|
var = p_ellipsoide modele x1 y1 z1 axe_x axe_y axe_z R r |
404 |
|
|
@brief Création d'une forme primitive volumique de type ellipsoide [MG_CG_FORME_VOLUME_ELLIPSOIDE] |
405 |
|
|
centré en (x1,y1,z1) orienté selon l'axe (axe_x,axe_y,axe_z) de rayon majeur R et de rayon mineur r |
406 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
407 |
|
|
@param x1 un double. |
408 |
|
|
@param y1 un double. |
409 |
|
|
@param z1 un double. |
410 |
|
|
@param axe_x un double. |
411 |
|
|
@param axe_y un double. |
412 |
|
|
@param axe_z un double. |
413 |
|
|
@param R un double. |
414 |
|
|
@param r un double. |
415 |
|
|
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
416 |
|
|
*/ |
417 |
|
|
|
418 |
|
|
|
419 |
francois |
737 |
//----------------------------------------------------------- |
420 |
|
|
|
421 |
|
|
/*! \page p_tore p_tore |
422 |
couturad |
906 |
var = p_tore modele x1 y1 z1 axe_x axe_y axe_z R r |
423 |
|
|
@brief Création d'une forme primitive volumique de type tore [MG_CG_FORME_VOLUME_TORE] basée sur le cercle centré en (x1,y1,z1) orienté selon l'axe (axe_x,axe_y,axe_z) de grand rayon R et de petit rayon r |
424 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
425 |
francois |
737 |
@param x1 un double. |
426 |
|
|
@param y1 un double. |
427 |
|
|
@param z1 un double. |
428 |
couturad |
906 |
@param axe_x un double. |
429 |
|
|
@param axe_y un double. |
430 |
|
|
@param axe_z un double. |
431 |
|
|
@param R un double. |
432 |
francois |
737 |
@param r un double. |
433 |
couturad |
906 |
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
434 |
francois |
737 |
*/ |
435 |
|
|
|
436 |
|
|
//----------------------------------------------------------- |
437 |
|
|
|
438 |
|
|
/*! \page p_cone p_cone |
439 |
couturad |
906 |
var = p_cone modele x1 y1 z1 axe_x axe_y axe_z h r |
440 |
|
|
@brief Création d'une forme primitive volumique de type cone [MG_CG_FORME_VOLUME_CONE] basée sur le cercle centré en (x1,y1,z1) orienté selon l'axe (axe_x,axe_y,axe_z) de hauteur h et de petit rayon r |
441 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
442 |
francois |
737 |
@param x1 un double. |
443 |
|
|
@param y1 un double. |
444 |
|
|
@param z1 un double. |
445 |
couturad |
906 |
@param axe_x un double. |
446 |
|
|
@param axe_y un double. |
447 |
|
|
@param axe_z un double. |
448 |
francois |
737 |
@param h un double. |
449 |
|
|
@param r un double. |
450 |
couturad |
906 |
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
451 |
francois |
737 |
*/ |
452 |
couturad |
740 |
//----------------------------------------------------------- |
453 |
francois |
737 |
|
454 |
couturad |
740 |
/*! \page ajoutass ajoutass |
455 |
couturad |
906 |
var = ajoutass assemblage forme |
456 |
|
|
@brief Ajoute une forme [MG_CG_FORME] dans un assemblage [MG_CG_ASSEMBLAGE] |
457 |
|
|
@param assemblage une variable représentant un assemblage [MG_CG_ASSEMBLAGE]. |
458 |
|
|
@param forme une variable représentant forme utile pour la construction géométrique [MG_CG_FORME] |
459 |
couturad |
740 |
@return var une variable inutilisable |
460 |
|
|
*/ |
461 |
francois |
737 |
//----------------------------------------------------------- |
462 |
|
|
/*! \page union union |
463 |
couturad |
906 |
var = union modele forme1 forme2 nom_semantique |
464 |
|
|
@brief Création d'une opération booléene d'union [MG_CG_OP_BOOL_UNION] entre deux formes [MG_CG_FORME] |
465 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
466 |
|
|
@param forme1 une forme [MG_CG_FORME]. |
467 |
|
|
@param forme2 une forme [MG_CG_FORME]. |
468 |
francois |
737 |
@param nom_semantique OPTIONNEL nom pour donner un sens physique à l'opération. Les choix sont : chamfrein, conge, intersection, percage, soustraction, union, extrusion |
469 |
couturad |
906 |
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
470 |
francois |
737 |
@details Si aucun nom de sémantique n'est donné le qualificatif SANS est attribué |
471 |
|
|
*/ |
472 |
couturad |
906 |
|
473 |
francois |
737 |
//----------------------------------------------------------- |
474 |
|
|
/*! \page diff diff |
475 |
couturad |
906 |
var = diff modele forme1 forme2 nom_semantique |
476 |
|
|
@brief Création d'une opération booléene de soustraction [MG_CG_OP_BOOL_DIFFERENCE] entre deux formes [MG_CG_FORME] |
477 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
478 |
|
|
@param forme1 une forme [MG_CG_FORME]. |
479 |
|
|
@param forme2 une forme [MG_CG_FORME]. |
480 |
francois |
737 |
@param nom_semantique OPTIONNEL nom pour donner un sens physique à l'opération. Les choix sont : chamfrein, conge, intersection, percage, soustraction, union, extrusion |
481 |
couturad |
906 |
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
482 |
francois |
737 |
@details Si aucun nom de sémantique n'est donné le qualificatif SANS est attribué |
483 |
|
|
*/ |
484 |
couturad |
906 |
|
485 |
francois |
737 |
//----------------------------------------------------------- |
486 |
|
|
/*! \page inter inter |
487 |
couturad |
906 |
var = inter modele forme1 forme2 nom_semantique |
488 |
|
|
@brief Création d'une opération booléene d'intersection [MG_CG_OP_BOOL_INTERSECTION] entre deux formes [MG_CG_FORME] |
489 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
490 |
|
|
@param forme1 une forme [MG_CG_FORME]. |
491 |
|
|
@param forme2 une forme [MG_CG_FORME]. |
492 |
francois |
737 |
@param nom_semantique OPTIONNEL nom pour donner un sens physique à l'opération. Les choix sont : chamfrein, conge, intersection, percage, soustraction, union, extrusion |
493 |
couturad |
906 |
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
494 |
francois |
737 |
@details Si aucun nom de sémantique n'est donné le qualificatif SANS est attribué |
495 |
|
|
*/ |
496 |
couturad |
906 |
|
497 |
francois |
737 |
//----------------------------------------------------------- |
498 |
couturad |
906 |
/*! \page fragment fragment |
499 |
|
|
var = fragment modele forme1 forme2 nom_semantique |
500 |
|
|
@brief Création d'une opération booléene de fragmentation [MG_CG_OP_BOOL_FRAGMENT] entre deux formes [MG_CG_FORME] |
501 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
502 |
|
|
@param forme1 une forme [MG_CG_FORME]. |
503 |
|
|
@param forme2 une forme [MG_CG_FORME]. |
504 |
|
|
@param nom_semantique OPTIONNEL nom pour donner un sens physique à l'opération. Les choix sont : chamfrein, conge, intersection, percage, soustraction, union, extrusion |
505 |
|
|
@return var une variable representant une forme utile pour la construction géométrique [MG_CG_FORME] |
506 |
|
|
@details Si aucun nom de sémantique n'est donné le qualificatif SANS est attribué |
507 |
francois |
737 |
*/ |
508 |
|
|
|
509 |
|
|
//----------------------------------------------------------- |
510 |
couturad |
740 |
|
511 |
couturad |
906 |
/*! \page creemodele creemodele |
512 |
|
|
var = creemodele gest nom unite eps fusion triSTL eps_tri |
513 |
|
|
@brief Création d'un modele de construction géométrique [MG_CG_MODELE] |
514 |
|
|
@param gest une variable représentant un modèle MAGiC [MG_GESTIONNAIRE] |
515 |
|
|
@param nom une chaine de caractère. |
516 |
|
|
@param unite un double représentant l'unité de mesure utilisé |
517 |
|
|
@param eps un double représentant la précision utilisé |
518 |
|
|
@param fusion une chaine de caractère. Si la valeur est "fusion" alors les entités similaires seront fusionnés à partir d'un algorithme de balayage. |
519 |
|
|
@param triSTL est une chaine de caractère. Si la valeur est "true" alors le modèle sera disposé à importer la triangulation STL. |
520 |
|
|
@param eps_tri un double représentant l'écart entre un entité topologique et son maillage STL. |
521 |
|
|
@return var une variable representant un modele de construction géométrique [MG_CG_MODELE] |
522 |
couturad |
740 |
*/ |
523 |
|
|
|
524 |
|
|
//----------------------------------------------------------- |
525 |
francois |
737 |
|
526 |
couturad |
906 |
/*! \page creeassemblage creeassemblage |
527 |
|
|
var = creeassemblage modele nom |
528 |
|
|
@brief Création d'un assemblage [MG_CG_ASSEMBLAGE] dans un modele de construction géométrique [MG_CG_MODELE] |
529 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
530 |
|
|
@param nom une chaine de caractère. Ce nom n'est pas utilisé par la suite |
531 |
|
|
@return var une variable representant un assemblage [MG_CG_ASSEMBLAGE] |
532 |
couturad |
740 |
*/ |
533 |
|
|
|
534 |
|
|
//----------------------------------------------------------- |
535 |
francois |
737 |
/*! \page enregistrer enregistrer |
536 |
|
|
var = enregistrer gest nom |
537 |
|
|
@brief Sauvegarde d'un fichier MAGiC |
538 |
couturad |
906 |
@param gest une variable représentant un gestionnaire MAGiC [MG_GESTIONNAIRE] |
539 |
francois |
737 |
@param nom une string représentant un nom de fichier (L'extension s'ajoute automatiquement au nom du fichier). |
540 |
|
|
@return var une variable inutilisable |
541 |
|
|
*/ |
542 |
|
|
|
543 |
|
|
//----------------------------------------------------------- |
544 |
couturad |
906 |
/*! \page brepassemblage brepassemblage |
545 |
|
|
var = brepassemblage assemblage nom |
546 |
|
|
@brief Conversion d'un assemblage [MG_CG_ASSEMBLAGE] au format brep d'OpenCascade |
547 |
|
|
@param gest une variable représentant un gestionnaire MAGiC [MG_GESTIONNAIRE] |
548 |
|
|
@param nom une string représentant un nom de fichier brep |
549 |
|
|
@return var une variable inutilisable |
550 |
|
|
*/ |
551 |
|
|
|
552 |
|
|
|
553 |
|
|
//----------------------------------------------------------- |
554 |
francois |
737 |
/*! \page export_mg_gmsh export_mg_gmsh |
555 |
|
|
var = export_mg_gmsh gest nom num |
556 |
|
|
@brief Exportation d'un MG_MAILLAGE MAGiC en format GMSH |
557 |
couturad |
906 |
@param gest une variable représentant un gestionnaire MAGiC [MG_GESTIONNAIRE] |
558 |
francois |
737 |
@param nom une string représentant un nom de fichier gmsh (L'extension s'ajoute automatiquement au nom du fichier). |
559 |
|
|
@param num un integer qui représente un numero sequantielle de maillage (entre 1 et n) |
560 |
|
|
@return var une variable inutilisable |
561 |
|
|
*/ |
562 |
|
|
|
563 |
|
|
//----------------------------------------------------------- |
564 |
|
|
/*! \page export_fem_gmsh export_fem_gmsh |
565 |
|
|
var = export_fem_gmsh gest nom num |
566 |
|
|
@brief Exportation d'un FEM_MAILLAGE MAGiC en format GMSH |
567 |
couturad |
906 |
@param gest une variable représentant un gestionnaire MAGiC [MG_GESTIONNAIRE] |
568 |
francois |
737 |
@param nom une string représentant un nom de fichier gmsh (L'extension s'ajoute automatiquement au nom du fichier). |
569 |
|
|
@param num un integer qui représente un numero sequantielle de maillage (entre 1 et n) |
570 |
|
|
@return var une variable inutilisable |
571 |
|
|
*/ |
572 |
|
|
//----------------------------------------------------------- |
573 |
|
|
/*! \page importstepocc importstepocc |
574 |
couturad |
906 |
var = importstepocc gest fichier_entree fichier_sortie eps fusion triSTL eps_tri |
575 |
|
|
@brief importation d'un fichier STEP en MAGiC via opencascade V2017 |
576 |
|
|
@param gest une variable représentant un gestionnaire MAGiC [MG_GESTIONNAIRE] |
577 |
|
|
@param fichier_entree une string représentant le nom de fichier STEP. |
578 |
|
|
@param fichier_sortie une string représentant le nom du fichier MAGiC. |
579 |
|
|
@param eps un double représentant la précision utilisé |
580 |
|
|
@param fusion une chaine de caractère. Si la valeur est "fusion" alors les entités similaires seront fusionnés à partir d'un algorithme de balayage. |
581 |
|
|
@param triSTL est une chaine de caractère. Si la valeur est "true" alors le modèle sera disposé à importer la triangulation STL. |
582 |
|
|
@param eps_tri un double représentant l'écart entre un entité topologique et son maillage STL. |
583 |
francois |
737 |
@return var une variable de type geometrie |
584 |
|
|
*/ |
585 |
couturad |
906 |
|
586 |
|
|
//----------------------------------------------------------- |
587 |
|
|
/*! \page tristlmodele tristlmodele |
588 |
|
|
var = tristlmodele modele type |
589 |
|
|
@brief Création de la triangulation STL [MG_MAILLAGE] d'un modele de construction géométrique [MG_CG_MODELE] |
590 |
|
|
@param modele une variable représentant un modele de conscrution géométrique [MG_CG_MODELE]. |
591 |
|
|
@param type une string représentant le type de triangulation. Si type = tristl la triangulation ne tiendra pas en compte de la conformité du maillage. Si type = triangulation la conformité du maillage sera assurée. |
592 |
|
|
@return var une variable inutilisable |
593 |
|
|
*/ |
594 |
|
|
|
595 |
|
|
//----------------------------------------------------------- |
596 |
|
|
/*! \page tristlassemblage tristlassemblage |
597 |
|
|
var = tristlassemblage assemblage type |
598 |
|
|
@brief Création de la triangulation STL [MG_MAILLAGE] d'un assemblage [MG_CG_ASSEMBLAGE] |
599 |
|
|
@param assemblage une variable représentant un assemblage [MG_CG_ASSEMBLAGE]. |
600 |
|
|
@param type une string représentant le type de triangulation. Si type = tristl la triangulation ne tiendra pas en compte de la conformité du maillage. Si type = triangulation la conformité du maillage sera assurée. |
601 |
|
|
@return var une variable inutilisable |
602 |
|
|
*/ |