Accueil | Cast3M

   1 : $$$$ @INCLUSI NOTICE  PASCAL    12/11/14    21:15:01     7568           
   2 :                                              DATE     12/11/14
   3 : 
   4 :   Procedure @INCLUSI                    Voir aussi : @P_VORO, @P_BOIT2
   5 :     -----------------  
   6 :     TAB2 = @INCLUSI TAB1 DPAR1 DEXC1 (DENS1) (ITRA1) ;
   7 : 
   8 :     Objet :
   9 :     -------
  10 : 
  11 :     La procedure @INCLUSI maille un echantilon numerique cubique
  12 :     d'un materiau constitue de particules spheriques de meme taille,
  13 :     en inclusions dans une matrice.
  14 : 
  15 :     Pour cela, elle s'appuie sur la partition de Voronoi des centres des
  16 :     particules, obtenue a l'aide des procedures @P_VORO et @P_BOIT2.
  17 : 
  18 :     Commentaire :
  19 :     -------------
  20 : 
  21 :     TAB1     = TABLE, resultat de la procedure @P_BOIT2 ;
  22 : 
  23 :     DPAR1    = FLOTTANT, diametre des particules ;
  24 : 
  25 :     DEXC1    = FLOTTANT, distance minimum entre centres des particules :
  26 :                DEXC1 doit etre strictement superieure a DPAR1 ;
  27 : 
  28 :     DENS1    = FLOTTANT, densite (taille) "moyenne" des elements du
  29 :                maillage, prise egale au quart de la taille moyenne des
  30 :                cellules de Voronoi par defaut. Toutefois, @INCLUSI
  31 :                raffine automatiquement le maillage pour avoir au moins
  32 :                2 elements finis dans chaque ligament de matrice, dont
  33 :                l'epaisseur minimale est egale a (DEXC1-DPAR1) ;
  34 : 
  35 :     ITRA1    = LOGIQUE, active des traces.
  36 : 
  37 :     TAB2     = TABLE, sous-indicee comme suit :
  38 :     . 'MAIL' = MAILLAGE des particules et de la matrice ;
  39 :     . 'PART' = MAILLAGE des particules ;
  40 :     . 'MATR' = MAILLAGE de la matrice ;
  41 :     . 'MPT'  = MAILLAGE de points, centres des particules ;
  42 : 
  43 :     Si PT1 est un point de TAB2 . 'MPT' alors :
  44 :     . PT1 . 'MAIL' = maillage de la particule de centre PT1 et de la
  45 :                      portion de matrice comprise dans la cellule de
  46 :                      Voronoi indicee par PT1 dans TAB1 ;
  47 :     . PT1 . 'PART' = maillage de la particule de centre PT1 ;
  48 :     . PT1 . 'MART' = maillage de la portion de matrice ;
  49 :     . PT1 . 'MPT'  = MAILLAGE de points, centres des particules voisines
  50 :                      a PT1 ;
  51 :     Si PT2 est un point de TAB2 . PT1 . 'MPT' alors :
  52 :     . PT1 . PT2 . 'MATR' : maillage de la face commune aux sous-maillages
  53 :                      des portions de matrice relatives aux particules de
  54 :                      centres PT1 et PT2 dans la partition de Voronoi
  55 :                      definie par TAB1.
  56 : 
  57 :     Remarques :
  58 :     -----------
  59 :        La procedure @INCLUSI ne reussit pas toujours a generer le mail-
  60 :     -lage demande, notamment lorsque des particules sont tangentes a une
  61 :     face ou a une arete du cube. Pour ameliorer sa robustesse, elle
  62 :     s'autorise a deplacer ou a aplanir legerement les particules qui
  63 :     posent probleme. L'utilisateur en est informe par un message. Il
  64 :     peut donc etre utile de faire une copie des affichages, en les
  65 :     redirigeant, par exemple, dans un fichier.
  66 :     
  67 : 
  68 : 
  69 :