Accueil | Cast3M

   1 : $$$$ PROI     NOTICE  GOUNAND   19/02/19    21:15:05     10114          
   2 :                                              DATE     19/02/19
   3 : 
   4 :    Operateur PROI                           Voir aussi :
   5 :     --------------  
   6 :    OBJ1 = PROI | GEO2   CHEL1               | (FLOT1)            | ; (1)
   7 :                | MOD2   (CARA) CHEL1  (MOT2)|                    |   (2)
   8 :                | MOD2   (CARA) CHEL1  'MINI' (ENT1)              |   (3)
   9 :                | 'POLY' GEO1 GEO2 CHPO1 ENT2 MOT1 ('POID' X1 X2) |   (4)
  10 : 
  11 : 
  12 : 
  13 :     Objet:
  14 :     ______
  15 : 
  16 :     L'operateur PROI effectue la projection des composants d'un MCHAML
  17 :     defini aux noeuds d'une geometrie donnee soit sur une nouvelle
  18 :     geometrie GEO2, soit sur les noeuds (par defaut) ou les points
  19 :     d'intégration d'un modele MOD2 (si MOT2 est present et different de
  20 :     'NOEUDS').
  21 : 
  22 :     Pour les syntaxes 1 et 2, les nouvelles valeurs sont obtenues par
  23 :     interpolation avec les fonctions de forme a partir des valeurs
  24 :     initiales et definissent un champ par points sur la nouvelle
  25 :     geometrie. Ce champ est de nature diffuse si le champ resultat est
  26 :     un champ par point.
  27 : 
  28 : 
  29 :     Commentaire :
  30 :     _____________
  31 : 
  32 :     OBJ1  :  Champ apres projection
  33 :              type CHPOINT  pour les syntaxes (1) et (4)
  34 :              type MCHAML   pour les syntaxes (2) et (3)
  35 : 
  36 :     GEO2  :  Nouvelle geometrie (type MAILLAGE) sur laquelle on veut
  37 :              projeter le champ CHEL1
  38 : 
  39 :     MOD2  :  Modele (type MMODEL) contenant les nouveaux points supports
  40 :              (noeuds ou points d'intégration) sur lesquels on veut
  41 :              projeter les valeurs du champ CHEL1
  42 : 
  43 :     CARA  :  Caracteristiques associees au modele MOD2 (type MCHAML,
  44 :              sous-type CARACTERISTIQUES) pour les elements coque
  45 :              (attention : on rappelle que ce champ doit etre donne avant
  46 :              CHEL1 si CHEL1 est également de sous-type CARACTERISTIQUES)
  47 : 
  48 :     CHEL1  : Champ par elements (type MCHAML) que l'on veut projeter. Il
  49 :              doit être defini aux NOEUDS du massif
  50 : 
  51 :     MOT2   : Mot precisant la localisation des point ou sera defini le
  52 :              champ. Par defaut aux 'NOEUDS' ou aux points d'integration
  53 :              au choix 'GRAVITE', 'RIGIDITE', 'MASSE' ou 'STRESSES'
  54 : 
  55 :     FLOT1  : Critère de rattrapage pour les points de GEO2 ou MOD2
  56 :              situés légèrement en-dehors du domaine de CHEL1
  57 :              Valeur par défaut : 1.e-5
  58 : 
  59 :     'MINI' : Si ce mot 'MINI' est present, l'evaluation est faite
  60 :              pour chaque composante par minimisation de l integrale
  61 :              (Ue-Us)**2 sur chaque element du support de sortie
  62 :              Ue : valeur de la composante au point courant dans
  63 :              le champ origine
  64 :              Us : valeur de la composante au point courant dans
  65 :              le champ resultat
  66 : 
  67 :     ENT1   : Si le mot-clef 'MINI' est present, on peut donner le nombre
  68 :              de points d'integration utilises par direction pour evaluer
  69 :              l'integrale (3 par defaut, i.e. 9 points pour une surface
  70 :              et 27 pour un volume) en précisant l'entier ENT1
  71 : 
  72 :     'POLY' : Si le mot-clef 'POLY' est present, l'operateur PROI
  73 :              effectue la projection sur GEO2 des composantes
  74 :              d'un CHPOINT CHPO1 defini sur GEO1. Il calcule aussi les
  75 :              derivees premieres et secondes du champ par rapport aux
  76 :              variables X et Y. Cette option marche en axisymmetrique et
  77 :              en dimension 2. (Voir remarque b pour plus de details)
  78 : 
  79 :     ENT2   : Si le mot-clef 'POLY' est present, l'entier ENT2 precise le
  80 :              type de symetrie :
  81 :                   1 : pas de symetrie
  82 :                   2 : symetrie par rapport a l'axe X
  83 :                       etc.....
  84 : 
  85 :     CHPO1  : Champ par point (type CHPOINT) a projeter
  86 : 
  87 :     GEO1   : Geometrie associee a CHPO1
  88 : 
  89 :     MOT1   : Mot 'PLAN' ou 'AXIS'
  90 : 
  91 :     'POID' : Mot-clef facultatif introduisant la donnee de deux poids
  92 :              pour l'interpolation des fonctions. Ils valent 1 par defaut
  93 :         X1 : Poids des premiers voisins
  94 :         X2 : Poids des secondes voisins
  95 : 
  96 : 
  97 : 
  98 :     Remarques :
  99 :     __________
 100 : 
 101 :     a) Le champ par element CHEL1 doit obligatoirement etre supporte par
 102 :        des element massifs
 103 : 
 104 :     b) Avec l'option POLY :
 105 :        L'operateur effectue la projection sur GEO2 des composantes d'un
 106 :        CHPOINT CHPO1 defini sur GEO1. Il calcule aussi les derivees
 107 :        premieres et secondes du champ par rapport aux variables X et Y.
 108 :        Cette option marche en axisymmetrique et en dimension 2.
 109 : 
 110 :        La methode POLY consiste pour chaque point de GEO2 a selectionner
 111 :        les premiers voisins et les seconds puis d'ajuster au mieux des
 112 :        fonctions polynomiales connaissant les valeurs aux points voisins
 113 :        et leur poids.
 114 : 
 115 :        Les fonctions polynomiales sont obtenues a partir du
 116 :        developpement en serie entiere des fonctions complexes.
 117 :        L'ajustement est fait par la methode des moindres carres (on
 118 :        attribue des poids aux premiers et aux seconds voisins).
 119 : 
 120 :        Les derivees sont obtenues en derivant les fonctions polynomiales.
 121 : 
 122 :        Le champ resultat a plusieurs composantes :
 123 :             A:   champ interpole.
 124 :             BX:  derivee suivant X (R en axisymmetrique).
 125 :             BY:  derivee suivant Y (Z en axisymmetrique).
 126 :             BXY: derivee seconde suivant X et Y (R et Z en axisymmetrie)
 127 :             BXX: derivee seconde suivant X.
 128 :             BYY: derivee seconde suivant Y.
 129 : 
 130 :        La methode POLY n'est donc pas l'interpolation naturelle utilisee
 131 :        habituellement en elements finis i.e. via les fonctions de forme.
 132 :        Cet operateur doit donc etre utilise avec precaution. Si un point
 133 :        est commun aux maillages GEO1 et GEO2, la valeur associee a ce
 134 :        point dans le CHPOINT resultat OBJ1 sera differente, en general,
 135 :        de la valeur associee a ce point dans le CHPOINT initial CHPO1 !!
 136 : 
 137 :