Accueil | Cast3M

   1 : $$$$ INDI     NOTICE  GOUNAND   26/01/09    21:15:36     12441          
   2 :                                              DATE     26/01/09
   3 : 
   4 :   Operateur INDI                           Voir aussi : TRIA REMA
   5 :     --------------  
   6 :     ----------------------  
   7 :     ----------------------
   8 : 
   9 :     CHAM1 = INDI GEO1 NOMi ...   ;
  10 : 
  11 :     Objet :
  12 :     -------
  13 : 
  14 :     L'operateur INDI fournit des indicateurs permettant d'apprecier
  15 :     la qualite d'un maillage.
  16 :     Les criteres testes sont specifies par des mots-cles.
  17 :     L'operateur cree un champ par elements (type MCHAML) dont les
  18 :     composantes sont ces criteres.
  19 : 
  20 :     Commentaires :
  21 :     --------------
  22 : 
  23 :     GEO1   : objet maillage (type MELEME)
  24 : 
  25 :     NOMi   : mot-cle (type mot). Les valeurs possibles sont :
  26 : 
  27 :    'PLAN'  : Critere de planearite pour les elements a 4 noeuds.
  28 :              La valeur calculee est directement proportionnelle a
  29 :              l'angle entre les normales N1 et N2 definies par :
  30 : 
  31 :              - N1 = produit vectoriel (S1 S2,S1 S4)
  32 :              - N2 = produit vectoriel (S3 S4,S3 S2)
  33 :                (S1,S2,S3,S4 sont les sommets de l'element)
  34 : 
  35 :              Les valeurs sont comprises entre 0 et 100.
  36 :              Un mchaml resultat de valeur constante egale a 100
  37 :              indiquera que tous les elements sont plans.
  38 : 
  39 :    'ASPE'  : Pour les elements TRI3 et TET4 seulement
  40 :              Rapport d'aspect (aspect ratio), il s'agit du rapport de la
  41 :              plus grande arete sur la plus petite hauteur, normalise
  42 : 
  43 :              - pour les TRI3 = SQRT(3)/2 * grande arete / petite hauteur
  44 :              - pour les TET4 = SQRT(2/3) * grande arete / petite hauteur
  45 : 
  46 :              Les valeurs sont comprises entre 1 (element equilateral)
  47 :              et + infini
  48 : 
  49 :    'SKEW'  : Pour les elements TRI3 et TET4 seulement
  50 :              Pointicite (skewness), il s'agit de l'ecart relatif entre la
  51 :              mesure de l'element et celle de l'element "optimal"
  52 : 
  53 :              SKEW = (mes_optimal - mes) / mes_optimal
  54 : 
  55 :              La mesure est l'aire (pour les TRI3) ou le volume (pour les TET4)
  56 : 
  57 :              L'element "optimal" est l'element equilateral place dans la
  58 :              meme sphere circonscrite que l'element etudie
  59 : 
  60 :              Les valeurs sont comprises entre 0 (element equilateral)
  61 :              et 1 (element degenere)
  62 : 
  63 : 
  64 :     ----------------------
  65 :     |  2eme possibilite  |
  66 :     ----------------------
  67 : 
  68 :      CHAM1 = INDI |'COHE'| MAIL1 |(METR1)      | ;
  69 :                   |'EQLT'|       |(METR1 MOD1) |
  70 :                   |'ISOD'|
  71 :                   |'DENS'|
  72 :                   |'TOPO'|
  73 : 
  74 :      Objet :
  75 :      _______
  76 : 
  77 :      L'operateur INDI fournit un indicateur de qualite d'un maillage
  78 :      MAIL1 eventuellement dans une metrique METR1 donnee.
  79 : 
  80 :      Commentaire :
  81 :      _____________
  82 : 
  83 :      MAIL1   : maillage (type MAILLAGE)
  84 :      METR1   : objet de type FLOTTANT, CHPOINT ou MCHAML
  85 :                Si METR1 est de type FLOTTANT, il s'agit de la taille
  86 :                d'arete voulue (densite)
  87 :                Si METR1 est de type CHPOINT ou MCHAML, il s'agit de
  88 :                l'inverse de la metrique voulue (unite : longueur^-2)
  89 :      MOD1    : modele associe a METR1 de type MCHAML
  90 : 
  91 :      CHAM1   : champ (type MCHAML) de qualite
  92 : 
  93 :      Remarques :
  94 :      ___________
  95 : 
  96 :      1) Si la metrique voulue est isotrope, le nom de composante est G.
  97 :         Si la metrique voulue est anisotrope, les noms des composantes
  98 :         sont :  G11, G21, G22, (G31, G32, G33 en 3D)
  99 :         Si la metrique n'est pas donnee, la metrique unite est supposee.
 100 : 
 101 :      2) Les indicateurs suivant sont disponibles :
 102 :         - 'COHE' : coherence des formes d'elements avec la metrique
 103 :                    METR1 (indicateur local) ;
 104 :         - 'EQLT' : equilateralite des elements (indicateur local)  ;
 105 :         - 'ISOD' : isodistribution des tailles d'elements (volume)
 106 :                    selon METR1 (indicateur global) ;
 107 :         - 'DENS' : coherence des tailles d'elements avec la metrique
 108 :                    METR1 (indicateur local) ;
 109 :         - 'TOPO' : combinaison de 'COHE' et 'DENS' utilisee par TRIA
 110 :                    TOPO et REMA.
 111 :         Tous les indicateurs sont homogènes à des rapports de longueur
 112 :         et varient entre 0 et 1 sauf 'ISOD' (rapports de volume de
 113 :         moyenne 1).
 114 : 
 115 :      3) L'indicateur 'TOPO' ne fonctionne que pour un maillage
 116 :         MAIL1 de TRI3 ou TET4.
 117 : 
 118 :      4) On peut demander a l'operateur plusieurs indicateurs en
 119 :         precisant plusieurs mots-cles :
 120 :         CHAM1 = INDI 'COHE' 'TOPO' MAIL1 ... ;
 121 : 
 122 :      5) CHAM1 a pour noms de composantes les noms des indicateurs.
 123 :