Télécharger topoptim.notice

Retour à la liste

Afficher cette notice en

Numérotation des lignes :
   1 : $$$$ TOPOPTIM NOTICE  CB215821  18/03/22    21:15:03     9786           
   2 :                                              DATE     18/03/22
   3 :                  
   4 : 
   5 :              CETTE PROCEDURE A ETE MISE GRACIEUSEMENT
   6 :               A DISPOSITION DE LA COMMUNAUTE CAST3M
   7 :                     PAR Guenhael Le Quilliec
   8 :               Laboratoire de Mecanique Gabriel Lame
   9 :  Universite de Tours, Universite d Orleans, INSA Centre Val de Loire
  10 :     Polytech Tours, 7 avenue Marcel Dassault, 37200 Tours, France
  11 : 
  12 : 
  13 :     Procedure TOPOPTIM                       Voir aussi : TOPOCHAN
  14 :     __________________                                    TOPOCRIT  
  15 :                                                           TOPOFILT  
  16 :                                                           TOPOMATE  
  17 :                                                           TOPORSTR  
  18 :                                                           TOPOSAUV  
  19 :                                                           TOPOSURF  
  20 :     TOPOPTIM TAB1 ;
  21 : 
  22 : 
  23 :     TAB1.  CONVERGENCE                     POIDS_ENERGIE_DEFO
  24 :            CRITERE                         POIDS_MECANISME
  25 :            CYCLE                           POIDS_TEMPERATURE
  26 :            CYCLES_SAUVEGARDES              PRECISION
  27 :            FACTEUR_D                       PROCEDURE_TOPOPERS
  28 :            FACTEUR_P                       RAPPORT_RAIDEURS_MECANIQUES
  29 :            FACTEUR_Q                       RAPPORT_RAIDEURS_THERMIQUES
  30 :            FRACTION_VOLUME                 RESOLUTION
  31 :            MAILLAGE                        RESOLUTION_LINEAIRE
  32 :            MAX_CYCLES                      RESOLUTION_PASAPAS
  33 :            MECANISME                       RESTRICTIONS
  34 :            MECANISME_ZERO_SPRING           SEUIL
  35 :            MES_SAUVEGARDES                 TAUX_FILTRAGE
  36 :            OC_B_MIN                        TOPOLOGIE
  37 :            OC_CRITERE                      TOPOLOGIE_MAX_INC
  38 :            OC_L1                           TOPOLOGIE_MIN
  39 :            OC_L2                           TRAC
  40 :            OC_L2_MIN                       ZERO_DIVISION
  41 :            OC_MAX_IT                       ZONE_FIGEE
  42 : 
  43 : 
  44 :     Objet :
  45 :     _______
  46 : 
  47 : Cette procedure permet d effectuer une optimisation topologique d une
  48 : structure soumise a un chargement mecanique et/ou thermique, en
  49 : considerant un comportement lineaire ou non-lineaire, avec ou sans
  50 : restrictions geometriques sur la topologie de sortie. Elle permet aussi
  51 : d effectuer une synthese d un mecanism souple.
  52 : 
  53 : 
  54 :     Commentaires :
  55 :     ______________
  56 : 
  57 : La premiere version de cette procedure s inspirait directement des
  58 : travaux de O. Sigmund ainsi que ceux de W. Hunter.
  59 : 
  60 : La version actuelle (2.1) a ete adaptee pour traiter les non-linearites 
  61 : (contact, plasticite, grandes deformations, grands deplacements et
  62 : grandes rotations), sous chargements multiphysiques. Il est possible
  63 : d ajouter des restrictions geometriques sur la topologie de sortie
  64 : (e.g. une periodicite). Elle est entierement ecrite en Gibiane afin de
  65 : faciliter les developpements pour les utilisateurs avances.
  66 : 
  67 : Le filtrage de la sensibilite est effectue ici par interpolations
  68 : successives entre les noeuds et les points d integration afin d assurer
  69 : des performances correctes tout en restant en language Gibiane. Cette
  70 : operation est assuree par la procedure TOPOFILT que les utilisateurs
  71 : avances peuvent redefinir avant de faire appel a TOPOPTIM.
  72 : 
  73 : Les restrictions geometriques, les mises a jour du materiau (pour tenir
  74 : compte de la densite) et les mises a jour de la topologie sont elles
  75 : aussi assurees par des procedures externes : TOPORSTR, TOPOMATE et
  76 : TOPOCRIT. La encore, les utilisateurs avances peuvent redefinir ces
  77 : procedures avant de faire appel a TOPOPTIM.
  78 : 
  79 : Enfin, une procedure personnelle TOPOPERS peut etre appelee a chaque
  80 : cycle d optimisation, juste apres la resolution, donnant a
  81 : l utilisateur la possibilite d intergir avec TOPOPTIM au cours du
  82 : processus d optimisation. Il peut par exemple appliquer sa propre
  83 : fonction objectif et calculer son propre champs de sensibilite, ou
  84 : bien encore creer ses propres restrictions geometriques...
  85 : 
  86 : 
  87 :     Remarques :
  88 :     ___________
  89 : 
  90 : Pour debuter, il est preferable de garder a 1 les facteurs de penalite
  91 : et d echelle de gris pour choisir convenablement le taux de filtrage.
  92 : Puis on ajustera dans un second temps le facteur de penalite pour
  93 : terminer par ajuster le facteur d echelle de gris.
  94 : 
  95 : Il est possible de stopper le processus a tout instant avec Ctrl+C et
  96 : analyser/modifier/sauvegarder les donnees de TAB1 (en mode interactif)
  97 : puis eventuellement continuer l optimisation en redonnant TAB1 a
  98 : TOPOPTIM.
  99 : 
 100 : Apres avoir utilise TOPOPTIM, la procedure TOPOSURF peut etre lancee
 101 : afin d obtenir une surface maillee lissee de la topologie de sortie.
 102 : Finallement, ce maillage peut etre sauve au format STL en utilisant la
 103 : directive SORTIR. Ce format de fichier est souvent prefere pour les
 104 : applications d impressions 3D.
 105 : 
 106 : 
 107 :     En entree :
 108 :     ___________
 109 : 
 110 : En entree, TAB1 sert a definir les options et les parametres de
 111 : l optimisation. Les indices de l objet TAB1 sont des mots (a ecrire en
 112 : toutes lettres, et en majuscules s ils sont mis entre cotes) dont voici
 113 : la liste :
 114 : 
 115 :  CRITERE : critere de convergence. Cette donnee est facultative et est
 116 :            egale a 0.01 par defaut.
 117 : 
 118 :  CYCLE : numero du precedent cycle a partir duquel l optimisation doit
 119 :          se poursuivre. Cette donnee est facultative et est egale a 0
 120 :          par defaut.
 121 : 
 122 :  CYCLES_SAUVEGARDES : LISTENTI contenant la liste des cycles pour
 123 :                       lesquels la procedure TOPOSAUV sauvegardera les
 124 :                       resultats specifies par MES_SAUVEGARDES. Cette
 125 :                       donnee est facultative et une LISTENTI vide par
 126 :                       defaut. Le dernier cycle est toujours sauvegarde,
 127 :                       qu il fasse ou non partie de LISTENTI.
 128 : 
 129 :  FACTEUR_D : FLOTTANT si constant ou EVOLITION du facteur
 130 :              d amortissement au cours des cycles (fonction contenant
 131 :              en abscisse les cycles et en ordonnee les valeurs du
 132 :              facteur d amortissement). Ce facteur permet d attenuer
 133 :              les phenomene oscillatoire lors des cycles d optimisation.
 134 :              Cette donnee est facultative et est egale a 0.5 par defaut.
 135 :              Une valeur plus faible est generalement souhaitable dans
 136 :              le cas d une synthese de mecanisme souple.
 137 : 
 138 :  FACTEUR_P : FLOTTANT si constant ou EVOLITION du facteur de penalite
 139 :              au cours des cycles (fonction contenant en abscisse les
 140 :              cycles et en ordonnee les valeurs du facteur de penalite).
 141 :              Cette donnee est facultative et est egale a 3.0 par defaut.
 142 : 
 143 :  FACTEUR_Q : FLOTTANT si constant ou EVOLITION du facteur d echelle de
 144 :              gris au cours des cycles d optimisation (fonction contenant
 145 :              en abscisse les cycles et en ordonnee les valeurs du
 146 :              facteur d echelle de gris). Cette donnee est facultative
 147 :              et est egale a 1.0 par defaut.
 148 : 
 149 :  FRACTION_VOLUME : fraction de volume autorisee. Cette donnee est
 150 :                    facultative et est egale a 0.4 par defaut.
 151 : 
 152 :  MAX_CYCLES : nombre maximum de cycles d optimisation. Cette donnee est
 153 :               facultative et est egale a 100 par defaut.
 154 : 
 155 :  MECANISME : cette donnee permet d effectuer une synthese de mecanisme
 156 :              souple. Elle indique la direction dans laquelle les
 157 :              deplacements doivent etre maximises sous l effet du
 158 :              chargement applique. Pour un cas de charge unique, la
 159 :              direction est donnee sous la forme d un simple objet
 160 :              forces. Pour un cas de charges multiples, les directions
 161 :              sont donnees par une TABLE d objets forces. Chaque objet
 162 :              forces est de type CHPOINT dans le cas d une resolution
 163 :              lineaire, ou sous la forme d un objet de type CHARGEME,
 164 :              cree par l operateur CHAR, dans le cas d une resolution
 165 :              non-lineaire. Notons que la synthese de mecanismes souples
 166 :              est encore en developpement et principalement a destination
 167 :              des utilisateurs avances pour le moment.
 168 : 
 169 :  MECANISME_ZERO_SPRING : LOGIQUE indiquant s il faut aussi effectuer
 170 :              une resolution sans les ressorts appliques dans le cas
 171 :              d une synthese de mecanisme souple. Cette donnee est
 172 :              facultative et est egale a FAUX par defaut.
 173 : 
 174 :  MES_SAUVEGARDES : TABLE permettant de specifier quel resultats doivent
 175 :                    etre sauves.
 176 : 
 177 :                    Les indices de cette table sont :
 178 : 
 179 :                    MAILLAGE : LOGIQUE valant VRAI pour sauver le
 180 :                               maillage au cours de l optimisation.
 181 :                               Cette donnee est facultative et est
 182 :                               egale a VRAI par defaut.
 183 : 
 184 :                    RESOLUTION : LOGIQUE valant VRAI pour sauver les
 185 :                                 resultats issus des resolutions
 186 :                                 effectuees au cours de l optimisation.
 187 :                                 Cette donnee est facultative et est
 188 :                                 egale a FAUX par defaut.
 189 : 
 190 :                    TOPOLOGIE : LOGIQUE valant VRAI pour sauver la
 191 :                                topologie au cours de l optimisation.
 192 :                                Cette donnee est facultative et est
 193 :                                egale a VRAI par defaut.
 194 : 
 195 :  OC_B_MIN : valeur minimale de B appliquee dans le cas d une synthese
 196 :             de mecanisme souple avec B**D correspondant au rapport
 197 :             entre le nouveau champ topologique et celui du cycle
 198 :             precedent. Cette donnee est facultative et est egale a
 199 :             1.0e-10 par defaut.
 200 : 
 201 :  OC_CRITERE : valeur du critere de convergence pour l algorithme
 202 :               dicotomique du critere d optimalite. Cette donnee est
 203 :               facultative et est egale a 1.0e-8 par defaut.
 204 : 
 205 :  OC_L1 : valeur de la borne inferieure intiale de l interval de
 206 :          recherche pour l algorithme dicotomique du critere
 207 :          d optimalite. Cette donnee est facultative et est egale a 0.0
 208 :          par defaut.
 209 : 
 210 :  OC_L2 : valeur de la borne superieure intiale de l interval de
 211 :          recherche pour l algorithme dicotomique du critere
 212 :          d optimalite. Cette donnee est facultative et est egale a
 213 :          1.0e5 par defaut.
 214 : 
 215 :  OC_L2_MIN : valeur minimale de la borne superieure de l interval de
 216 :              recherche pour l algorithme dicotomique du critere
 217 :              d optimalite. Cette donnee est facultative et est egale a
 218 :              1.0e-40 par defaut.
 219 : 
 220 :  OC_MAX_IT : nombre maximum d iterations pour l algorithme dicotomique
 221 :              du critere d optimalite. Cette donnee est facultative et
 222 :              est egale a 1.0e5 par defaut.
 223 : 
 224 :  POIDS_TEMPERATURE : poids permettant d affecter la sensibilite a la
 225 :                      minimisation de souplesse thermique par rapport
 226 :                      aux autres criteres de sensibilite eventuels.
 227 :                      Cette donnee est facultative et est egale par
 228 :                      defaut a 0.0 dans le cas d une synthese de
 229 :                      mecanisme souple ou 1.0 sinon. Pour ne pas etre
 230 :                      sensible a la souplesse thermique, mettre ce
 231 :                      parametre a zero.
 232 : 
 233 :  POIDS_MECANISME : poids permettant d affecter la sensibilite du critere
 234 :                    pour la synthese de mecanisme souple par rapport aux
 235 :                    autres criteres de sensibilite eventuels. Cette
 236 :                    donnee est facultative et est egale a 1.0 par defaut.
 237 :                    Pour ne pas etre sensible a ce critere, mettre ce
 238 :                    parametre a zero.
 239 : 
 240 :  POIDS_ENERGIE_DEFO : poids permettant d affecter la sensibilite a la
 241 :                       minimisation de l energie de deformation
 242 :                       elastique par rapport aux autres criteres de
 243 :                       sensibilite eventuels. Cette donnee est
 244 :                       facultative et est egale par defaut a 0.0 dans le
 245 :                       cas d une synthese de mecanisme souple ou 1.0
 246 :                       sinon. Pour ne pas etre sensible a l energie de
 247 :                       deformation elastique, mettre ce parametre a
 248 :                       zero.
 249 : 
 250 :  PRECISION : seuil de precision en deca duquelle on considere qu un
 251 :              nombre est nul. Cette donnee est facultative et est egale
 252 :              a 1.0e-9 par defaut.
 253 : 
 254 :  PROCEDURE_TOPOPERS : LOGIQUE indiquant si une procedure TOPOPERS a
 255 :                       ete definie par l utilisateur. Auquel cas elle
 256 :                       sera appelee a chaque cycle du processus
 257 :                       d optimisation, juste apres la resolution. Cette
 258 :                       donnee est facultative et est egale a FAUX par
 259 :                       defaut.
 260 : 
 261 :  RAPPORT_RAIDEURS_MECANIQUES : valeur du rapport entre la raideur
 262 :                                correspondant aux vides et celle
 263 :                                de la matiere. Cette donnee est
 264 :                                facultative et est egale par defaut a 
 265 :                                1.0e-8 dans le cas d une synthese de
 266 :                                mecanisme souple ou 0.0 sinon. Une valeur
 267 :                                differente de 0.0 revient a considerer
 268 :                                que les vides ont une certaine raideur.
 269 : 
 270 :  RAPPORT_RAIDEURS_THERMIQUES : valeur du rapport entre la conductivite
 271 :                                correspondant aux vides et celle
 272 :                                de la matiere. Cette donnee est
 273 :                                facultative et est egale a 0.001 par
 274 :                                defaut. Par exemple, si les vident
 275 :                                correspondent a de l air, ce rapport
 276 :                                doit etre celui de la conductivite de
 277 :                                l air sur celle de la matiere.
 278 : 
 279 :  RESOLUTION_LINEAIRE : TABLE contenant les donnees de la simulation
 280 :                        pour une resolution lineaire via la procedure
 281 :                        RESO. Si cette table est donnee, alors
 282 :                        RESOLUTION_PASAPAS ne doit pas etre fourni.
 283 :                        
 284 :        Les indices de cette table sont :
 285 :                        
 286 :        BLOCAGES_MECANIQUES : blocages mecaniques (type RIGIDITE) ou
 287 :                              TABLE contenant les differents cas de
 288 :                              blocages mecaniques.
 289 : 
 290 :        BLOCAGES_THERMIQUES : blocages thermiques (type RIGIDITE) ou
 291 :                              TABLE contenant les differents cas de
 292 :                              blocages thermiques.
 293 : 
 294 :        CARACTERISTIQUES : champ de caracteristiques materielles et
 295 :                           geometriques si necessaire.
 296 : 
 297 :        CHARGEMENT : objet de type CHPOINT pour un cas de charge unique,
 298 :                     ou bien de type TABLE contenant les differents
 299 :                     cas de chargement, chacun de type CHPOINT. Dans le
 300 :                     cas d un chargement thermo-mecanique, la resolution
 301 :                     sera decouplee. Pour une resolution couplee,
 302 :                     utiliser RESOLUTION_PASAPAS.
 303 : 
 304 :        MODELE : objet modele (type MMODEL) de la structure de
 305 :                 comportement lineaire.
 306 : 
 307 :  RESOLUTION_PASAPAS : TABLE contenant les donnees de la simulation pour
 308 :                       une resolution non-lineaire via la procedure
 309 :                       PASAPAS. Si cette table est donnee, alors
 310 :                       RESOLUTION_LINEAIRE ne doit pas etre fourni.
 311 : 
 312 :        Les indices de cette table sont ceux decrits dans la procedure
 313 :        PASAPAS a l exception des indices suivants qui peuvent etre
 314 :        donnes sous forme de TABLE pour les cas de chargements
 315 :        multiples :
 316 :        
 317 :        BLOCAGES_MECANIQUES : blocages mecaniques (type RIGIDITE) ou
 318 :                              TABLE contenant les differents cas de
 319 :                              blocages mecaniques, chacun de type
 320 :                              RIGIDITE.
 321 : 
 322 :        BLOCAGES_THERMIQUES : blocages thermiques (type RIGIDITE) ou
 323 :                              TABLE contenant les differents cas de
 324 :                              blocages thermiques, chacun de type
 325 :                              RIGIDITE.
 326 : 
 327 :        CHARGEMENT : chargement en fonction du parametre d evolution
 328 :                     (type CHARGEME, cree par l operateur CHAR), ou
 329 :                     TABLE contenant les differents cas de chargement,
 330 :                     chacun de type CHARGEME.
 331 : 
 332 :  RESTRICTIONS : TABLE contenant les donnees sur les restrictions
 333 :                 geometriques a appliquer. Cette donnee est facultative
 334 :                 et aucune restriction n est applique par defaut.
 335 :                
 336 :                 Les indices de cette TABLE sont des entiers
 337 :                 correspondant aux numeros, de 1 a N, des N restrictions
 338 :                 a appliquer. Les donnees relatives a chaque restriction
 339 :                 sont elles-meme indiquees dans une TABLE.
 340 : 
 341 :        Les indices de la TABLE d une restriction sont :
 342 :        
 343 :        ZONE : partie du maillage sur lequel s applique la restriction.
 344 :               Ce maillage doit lui meme respecter la restriction a
 345 :               imposer a la topologie de sortie. Cette donnee est
 346 :               facultative et la restriction s applique par defaut sur
 347 :               l ensemble du domaine optimise.
 348 :        
 349 :        TYPE : MOT indiquant le type de restriction a appliquer. Les
 350 :               valeurs possibles sont :
 351 : 
 352 :                 1 - SYME_POIN : pour une symetrie centrale.
 353 :                 2 - SYME_DROI : pour une symetrie axiale.
 354 :                 3 - SYME_PLAN : pour une symetrie plane.
 355 :                 4 - PERIODICITE_AXIALE : pour une periodicite axiale
 356 :                 5 - PERIODICITE_CIRCULAIRE : pour une periodicite
 357 :                                              circulaire.
 358 :        
 359 :        POIN1 : un, deux ou trois points suivant le type de restriction.
 360 :        POIN2   Par exemple, pour SYME_PLAN, en 2D, seuls POIN1 et POIN2
 361 :        POIN3   sont attendus. Pour le cas de la periodicite axiale,
 362 :                POIN1 et POIN2 donnent la direction et leur distance
 363 :                doit correspondre a la taille d une periode.
 364 :        
 365 :        PERIODES : nombre de periodes. Cette donnee est attendue
 366 :                   uniquement pour PERIODICITE_AXIALE et pour
 367 :                   PERIODICITE_CIRCULAIRE.
 368 : 
 369 :  SEUIL : valeur seuil du rapport des raideurs en deca duquelle les
 370 :          elements sont elimines pour accelerer le processus
 371 :          d optimisation. Cette donnee est facultative et est egale a
 372 :          1.0e-9 par defaut.
 373 : 
 374 :  TAUX_FILTRAGE : taux de filtrage de la sensibilite. Cette donnee est
 375 :                  facultative et est egale a 1 par defaut.
 376 : 
 377 :  TOPOLOGIE : TABLE facultative qui peut contenir une topologie de
 378 :              depart imposee dont l indice doit etre egal a CYCLE
 379 :              (0 par defaut). Le maillage et le modele associes a ce
 380 :              champ topologique de depart n ont pas besoin d etre les
 381 :              memes que ceux donnes dans la table de resolution. Cela
 382 :              peut etre utile dans le cas d un raffinement de maillage
 383 :              entre plusieurs appels successifs a TOPOPTIM.
 384 : 
 385 : TABLE that can contain an initial topology field whose
 386 :             index should be equal to CYCLE (0 by default). The
 387 :             mesh and model associated to this prescribed initial
 388 :             topology field do not need to be the same as those given
 389 :             in the resolution table. This can be useful in case of
 390 :             mesh refinement between several calls to TOPOPTIM. This
 391 :             data is optional. 
 392 : 
 393 : 
 394 :  TOPOLOGIE_MAX_INC : increment maximum d evolution du critere
 395 :                      d optimalite entre deux cycles d optimisation
 396 :                      successifs. Cette donnee est facultative et est
 397 :                      egale a 0.2 par defaut.
 398 : 
 399 :  TOPOLOGIE_MIN: valeur minimale du champ topologique. Cette donnee est
 400 :                 facultative et est egale par defaut a 0.001 dans le
 401 :                 cas d une synthese de mecanisme souple ou 0.0 sinon.
 402 : 
 403 :  TRAC : LOGIQUE permettant de dessiner l evolution de la topologie a
 404 :         la fin de chaque cycle du processus d optimisation. Cette
 405 :         donnee est facultative et est egale a VRAI par defaut.
 406 : 
 407 :  ZERO_DIVISION : valeur utilisee pour eviter les divisions par zero
 408 :                  lors des operations sur les champs de sensibilite.
 409 :                  Cette donnee est facultative et est egale a 0.001 par
 410 :                  defaut.
 411 : 
 412 :  ZONE_FIGEE : partie du maillage qui fera obligatoirement partie de la
 413 :               topologie optimisee. Cette donnee est facultative.
 414 : 
 415 : 
 416 :     En sortie :
 417 :     ___________
 418 : 
 419 : En sortie, TAB1 permet de retrouver differentes donnes dont voici
 420 : la liste :
 421 : 
 422 :  CONVERGENCE : LOGIQUE indiquant si le critere de convergence a ete
 423 :                atteint. Dans le cas contraire, le nombre maximum
 424 :                de cycles d optimisation a ete atteint.
 425 : 
 426 :  MAILLAGE : TABLE contenant le maillage au cours de l optimisation.
 427 :             Si la sauvegarde du maillage n a pas ete specifiee dans
 428 :             MES_SAUVEGARDES, alors cette table est vide. Dans le cas
 429 :             contraire, les indices de cette table sont ceux des cycles
 430 :             auquels la sauvegarde a ete demandee ainsi que celui du
 431 :             dernier cycle d optimisation effectue.
 432 : 
 433 :  RESOLUTION : TABLE contenant les resultats issus des resolutions
 434 :               effectuees au cours de l optimisation. Si la sauvegarde
 435 :               des resolutions n a pas ete specifiee dans
 436 :               MES_SAUVEGARDES, alors cette table est vide. Dans le cas
 437 :               contraire, les indices de cette table sont ceux des cycles
 438 :               auquels la sauvegarde a ete demandee ainsi que celui du
 439 :               dernier cycle d optimisation effectue. 
 440 : 
 441 :  TOPOLOGIE : TABLE contenant la topologie au cours de l optimisation.
 442 :              Si la sauvegarde de la topologie n a pas ete specifiee dans
 443 :              MES_SAUVEGARDES, alors cette table contient uniquement la
 444 :              topologie du dernier cycle d optimisation effectue. Dans le
 445 :              cas contraire, les indices de cette table sont ceux des
 446 :              cycles auquels la sauvegarde a ete demandee ainsi que celui
 447 :              du dernier cycle d optimisation effectue. Chaque topologie
 448 :              correspond a un champ scalaire de type MCHAML, exprime aux
 449 :              centres de gravite des elements du maillage initial.
 450 : 
 451 : 
 452 :     Exemples :
 453 :     __________
 454 : 
 455 : topoptim_01.dgibi topoptim_02.dgibi topoptim_03.dgibi topoptim_04.dgibi
 456 : topoptim_05.dgibi topoptim_06.dgibi topoptim_07.dgibi topoptim_08.dgibi
 457 : topoptim_09.dgibi topoptim_10.dgibi topoptim_11.dgibi
 458 : 
 459 : 

© Cast3M 2003 - Tous droits réservés.
Mentions légales