Notices | Cast3M

   1 : $$$$ INDI     NOTICE  GOUNAND   25/07/24    21:15:03     12334          
   2 :                                              DATE     25/07/24
   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 'TOPO' MAIL1 (METR1) | ('ARIT') | ('LISTREEL') ;
  69 :                                        |  'GEOM'  |
  70 : 
  71 :      Objet :
  72 :      _______
  73 : 
  74 :      L'operateur INDI fournit un indicateur de qualite d'un maillage
  75 :      de TRI3 ou de TET4, eventuellement dans une metrique donnee.
  76 : 
  77 :      Commentaire :
  78 :      _____________
  79 : 
  80 :      MAIL1   : maillage (type MAILLAGE) constitué d'un seul type
  81 :                d'elements TRI3 ou TET4
  82 :      METR1   : objet de type FLOTTANT ou CHPOINT
  83 :                Si METR1 est de type FLOTTANT, il s'agit de la taille
  84 :                d'arete voulue (densite)
  85 :                Si METR1 est de type CHPOINT, il s'agit de l'inverse
  86 :                de la metrique voulue (unite : longueur^-2)
  87 : 
  88 :      CHAM1   : champ (type MCHAML ou LISTREEL) de qualite
  89 : 
  90 :      Remarques :
  91 :      ___________
  92 : 
  93 :      1) Si la metrique voulue est isotrope, le nom de composante est G.
  94 :         Si la metrique voulue est anisotrope, les noms des composantes
  95 :         sont :  G11, G21, G22, (G31, G32, G33 en 3D)
  96 : 
  97 :      2) En presence du mot-cle LISTREEL, l'operateur renvoie un objet
  98 :         de ce type plutot que du type MCHAML.
  99 : 
 100 :      3) L'indicateur varie entre 0 (element degenere) et 1
 101 :         (element equilateral, de cote 1 dans la métrique voulue si
 102 :         METR1 est donnee)
 103 : 
 104 :      4) Si le mot-clef 'ARIT' (par defaut), les tenseurs de metrique
 105 :         sont interpoles lineairement dans l'element.
 106 :         Si le mot-clef 'GEOM' est donne, les tenseurs de metrique
 107 :         sont interpoles geometriquement (Log-euclidien).
 108 :