Notices | Cast3M

$$$$ MAILTOPO NOTICE  GOUNAND   21/04/09    21:15:04     10944          
                                             DATE     21/04/09

    Procedure MAILTOPO                      Voir aussi : TRIA REMA
    ------------------                                   INDI  
    MAIL1 (METR1) = MAILTOPO | 'TRIA' MAIL2         | (METR2) ...
                             | 'REMA' MAIL3 (MAIL4) | 

                 ... | ('AJNO') | ('IPOL') (TAB1) ;
                     |  'NOAJ'  |
    Objet :
    _______

    Cette procédure n'est pas destinee a etre appelee par l'utilisateur
    (voir operateurs TRIA et REMA).
    Elle implemente un algorithme topologique de génération ('TRIA') ou
    d'optimisation ('REMA') d'un maillage de simplex anisotrope du a T.
    Coupez et al. (voir bibliographie ci-apres)

    Commentaire :
    _____________

      MAIL1 : maillage genere (type MAILLAGE)

      METR1 : si le mot-cle 'IPOL' est donne, METR1 (type CHPOINT) est
              la metrique interpolee sur le nouveau maillage MAIL1

      MAIL2 : bord du maillage a generer (type MAILLAGE)

      MAIL3 : maillage a optimiser (type MAILLAGE)

      MAIL4 : partie du bord de MAIL3 que le mailleur ne doit pas
              modifier (par defaut le mailleur peut retirer ou ajouter
              des noeuds sur les parties planes du bord de MAIL3)

      METR2 : objet de type FLOTTANT ou CHPOINT               
              Si METR2 est de type FLOTTANT, il s'agit de la taille
              d'arete voulue (densite)
              Si METR2 est de type CHPOINT, il s'agit de l'inverse
              de la metrique voulue (unite : longueur^-2)

      TAB1  : objet optionnel de type TABLE dont les indices sont des
              parametres d'entree ou de sortie du mailleur
              Entree :
              TAB1 . 'debug' = 0..2 ; (0 par defaut)
                   Niveau d'information
              TAB1 . 'graph' = faux..vrai ; (faux par defaut)
                   Sortie de graphiques
              TAB1 . 'id_cas' = MOT1 ;
                   Chaine de caracteres decrivant le cas
              TAB1 . 'max_iter' = ENTI1 ; (100 par defaut)
                   Nombre d'iterations maximum de l'algorithme
              TAB1 . 'precrel_volume' = FLOT1 ; (1.D-11 par defaut)
                   Precision relative en-dessous de laquelle un element
                   est considere comme ayant un volume nul
              TAB1 . 'precrel_qualite' = FLOT1 ; (1.D-2 par defaut)
                   Precision relative en-dessous de laquelle deux
                   elements sont consideres comme ayant une qualite
                   identique
              TAB1 . 'verif' = 0..2 ; (0 par defaut)
                   Niveau de verification (debug) dans l'operateur
                   appele (OPTO)
              TAB1 . 'impr_segadj' = 0..1 ; (0 par defaut)
                   Impressions (debug) lors des ajustements de segments
                   dans l'operateur appele (OPTO)
              Sortie :
              TAB1 . 'curtopo' = MAIL1 ;
                   Maillage courant obtenu par l'algorithme
              TAB1 . 'dvol' = FLOT1 ;
                   Difference entre le volume du maillage courant et le
                   volume souhaite
              TAB1 . 'nnul' = ENTI1 ;
                   Nombre d'elements du maillage courant ayant un volume
                   nul
              TAB1 . 'miq' = FLOT1 ;
                   Qualite minimum (cf. operateur INDI) des elements du
                   maillage courant
              TAB1 . 'moq' = FLOT1 ;
                   Qualite moyenne des elements du maillage courant
              TAB1 . 'maq' = FLOT1 ;
                   Qualite maximal des elements du maillage courant
                   

    Remarques :
    ___________

     1) Si la metrique voulue est isotrope, le nom de composante est G.
        Si la metrique voulue est anisotrope, les noms des composantes
        sont :  G11, G21, G22, (G31, G32, G33 en 3D)
        
     2) Si le mot-clef 'AJNO' (par defaut) est donne, le mailleur peut
        generer de nouveaux noeuds.
        Si le mot-clef 'NOAJ' est donne, le mailleur ne genere pas de
        nouveaux noeuds.


    Bibliographie :
    _______________

    @article{author = {Coupez, Thierry},
         title  = {Generation de maillage et adaptation de maillage par
                  optimisation locale},
         journal = {Revue Européenne des Éléments Finis},
         volume  = {9}, number = {4}, pages = {403-423}, year = {2000},
         doi     = {10.1080/12506559.2000.10511454}}
     

    @PhdThesis{author = {Cyril Gruau},
           title =  {Generation de métriques pour adaptation anisotrope
              de maillage, application à la mise en forme des matériaux}, 
              school = {ENSMP}, year = {2004}}
      
    @article{author = "Cyril Gruau and Thierry Coupez",
         title = "3D tetrahedral, unstructured and anisotropic mesh
           generation with adaptation to natural and multidomain metric", 
         journal = "Computer Methods in Applied Mechanics and Engineering",
         volume = "194", number = "48 - 49", pages = "4951 - 4976",
         year = "2005",
         doi = "10.1016/j.cma.2004.11.020"}