Accueil | Cast3M

   1 : $$$$ MAILTOPO NOTICE  GOUNAND   24/09/27    21:15:13     12019          
   2 :                                              DATE     24/09/27
   3 : 
   4 :   Procedure MAILTOPO                      Voir aussi : TRIA REMA
   5 :     ------------------                                   INDI  
   6 :     MAIL1 (METR1) = MAILTOPO | 'TRIA' MAIL2         | (METR2) ...
   7 :                              | 'REMA' MAIL3 (MAIL4) |
   8 : 
   9 :                  ... | ('AJNO') | ('IPOL') | ('ARIT') | (TAB1) ;
  10 :                      |  'NOAJ'  |          |  'GEOM'  |
  11 :     Objet :
  12 :     _______
  13 : 
  14 :     Cette procédure n'est pas destinee a etre appelee par l'utilisateur
  15 :     (voir operateurs TRIA et REMA).
  16 :     Elle implemente un algorithme topologique de génération ('TRIA') ou
  17 :     d'optimisation ('REMA') d'un maillage de simplex anisotrope du a T.
  18 :     Coupez et al. (voir bibliographie ci-apres)
  19 : 
  20 :     Commentaire :
  21 :     _____________
  22 : 
  23 :       MAIL1 : maillage genere (type MAILLAGE)
  24 : 
  25 :       METR1 : si le mot-cle 'IPOL' est donne, METR1 (type CHPOINT) est
  26 :               la metrique interpolee sur le nouveau maillage MAIL1
  27 : 
  28 :       MAIL2 : bord du maillage a generer (type MAILLAGE)
  29 : 
  30 :       MAIL3 : maillage a optimiser (type MAILLAGE)
  31 : 
  32 :       MAIL4 : partie du bord de MAIL3 que le mailleur ne doit pas
  33 :               modifier (par defaut le mailleur peut retirer ou ajouter
  34 :               des noeuds sur les parties planes du bord de MAIL3)
  35 : 
  36 :       METR2 : objet de type FLOTTANT ou CHPOINT
  37 :               Si METR2 est de type FLOTTANT, il s'agit de la taille
  38 :               d'arete voulue (densite)
  39 :               Si METR2 est de type CHPOINT, il s'agit de l'inverse
  40 :               de la metrique voulue (unite : longueur^-2)
  41 : 
  42 :       TAB1  : objet optionnel de type TABLE dont les indices sont des
  43 :               parametres d'entree ou de sortie du mailleur
  44 :               Entree :
  45 :               TAB1 . 'debug' = 0..2 ; (0 par defaut)
  46 :                    Niveau d'information
  47 :               TAB1 . 'graph' = faux..vrai ; (faux par defaut)
  48 :                    Sortie de graphiques
  49 :               TAB1 . 'id_cas' = MOT1 ;
  50 :                    Chaine de caracteres decrivant le cas
  51 :               TAB1 . 'max_iter' = ENTI1 ; (100 par defaut)
  52 :                    Nombre d'iterations maximum de l'algorithme
  53 :               TAB1 . 'precrel_volume' = FLOT1 ; (1.D-11 par defaut)
  54 :                    Precision relative en-dessous de laquelle un element
  55 :                    est considere comme ayant un volume nul
  56 :               TAB1 . 'precrel_qualite' = FLOT1 ; (1.D-2 par defaut)
  57 :                    Precision relative en-dessous de laquelle deux
  58 :                    elements sont consideres comme ayant une qualite
  59 :                    identique
  60 :               TAB1 . 'verif' = 0..2 ; (0 par defaut)
  61 :                    Niveau de verification (debug) dans l'operateur
  62 :                    appele (OPTO)
  63 :               TAB1 . 'impr_segadj' = 0..1 ; (0 par defaut)
  64 :                    Impressions (debug) lors des ajustements de segments
  65 :                    dans l'operateur appele (OPTO)
  66 :               Sortie :
  67 :               TAB1 . 'curtopo' = MAIL1 ;
  68 :                    Maillage courant obtenu par l'algorithme
  69 :               TAB1 . 'dvol' = FLOT1 ;
  70 :                    Difference entre le volume du maillage courant et le
  71 :                    volume souhaite
  72 :               TAB1 . 'nnul' = ENTI1 ;
  73 :                    Nombre d'elements du maillage courant ayant un volume
  74 :                    nul
  75 :               TAB1 . 'miq' = FLOT1 ;
  76 :                    Qualite minimum (cf. operateur INDI) des elements du
  77 :                    maillage courant
  78 :               TAB1 . 'moq' = FLOT1 ;
  79 :                    Qualite moyenne des elements du maillage courant
  80 :               TAB1 . 'maq' = FLOT1 ;
  81 :                    Qualite maximal des elements du maillage courant
  82 : 
  83 : 
  84 :     Remarques :
  85 :     ___________
  86 : 
  87 :      1) Si la metrique voulue est isotrope, le nom de composante est G.
  88 :         Si la metrique voulue est anisotrope, les noms des composantes
  89 :         sont :  G11, G21, G22, (G31, G32, G33 en 3D)
  90 : 
  91 :      2) Si le mot-clef 'AJNO' (par defaut) est donne, le mailleur peut
  92 :         generer de nouveaux noeuds.
  93 :         Si le mot-clef 'NOAJ' est donne, le mailleur ne genere pas de
  94 :         nouveaux noeuds.
  95 : 
  96 :      3) Si le mot-clef 'ARIT' (par defaut) est donne, l'interpolation
  97 :         de la metrique METR2 est lineaire.
  98 :         Si le mot-clef 'GEOM' est donne, l'interpolation de la metrique
  99 :         METR2 est geometrique.
 100 : 
 101 :     Bibliographie :
 102 :     _______________
 103 : 
 104 :     @article{author = {Coupez, Thierry},
 105 :          title  = {Generation de maillage et adaptation de maillage par
 106 :                   optimisation locale},
 107 :          journal = {Revue Européenne des Éléments Finis},
 108 :          volume  = {9}, number = {4}, pages = {403-423}, year = {2000},
 109 :          doi     = {10.1080/12506559.2000.10511454}}
 110 : 
 111 : 
 112 :     @PhdThesis{author = {Cyril Gruau},
 113 :            title =  {Generation de métriques pour adaptation anisotrope
 114 :               de maillage, application à la mise en forme des matériaux},
 115 :               school = {ENSMP}, year = {2004}}
 116 : 
 117 :     @article{author = "Cyril Gruau and Thierry Coupez",
 118 :          title = "3D tetrahedral, unstructured and anisotropic mesh
 119 :            generation with adaptation to natural and multidomain metric",
 120 :          journal = "Computer Methods in Applied Mechanics and Engineering",
 121 :          volume = "194", number = "48 - 49", pages = "4951 - 4976",
 122 :          year = "2005",
 123 :          doi = "10.1016/j.cma.2004.11.020"}
 124 :