Télécharger la plaquette Cast3M


fr en

Télécharger elements_vtk.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : element_vtk.dgibi
  2. *
  3. OPTI 'DIME' 3 'ELEM' 'CUB8' ;
  4. ************************************************************************
  5. ************************************************************************
  6.  
  7. GRAPH = FAUX ;
  8.  
  9. * repertoire des fichiers "divers"
  10. DIVERS = VENV 'CASTEM_DIVERS';
  11. *
  12.  
  13. * ****************************************
  14. * CREATION DES DIFFERENTS TYPES D'ELEMENTS
  15. * ****************************************
  16.  
  17. P0   = 0. 0. 0. ;
  18. PX   = 1. 0. 0. ;
  19. PY   = 0. 1. 0. ;
  20. PZ   = 0. 0. 1. ;
  21. PXY  = 1. 1. 0. ;
  22. PXZ  = 1. 0. 1. ;
  23. PYZ  = 0. 1. 1. ;
  24. PXYZ = 1. 1. 1. ;
  25.  
  26. * Elements lineaires
  27. * ==================
  28. _POI1 = P0 ;
  29. _SEG2 = MANU 'SEG2' P0 PX ;
  30. _TRI3 = MANU 'TRI3' P0 PX PY ;
  31. _QUA4 = MANU 'QUA4' P0 PX PXY PY ;
  32. _CUB8 = _QUA4 VOLU 'TRAN' 1 PZ ;
  33. _TET4 = MANU 'TET4' P0 PX PY PZ ;
  34. _PYR5 = MANU 'PYR5' (0.5*PX) (0.5*PY) (-0.5*PX) (-0.5*PY) PZ ;
  35. _PRI6 = _TRI3 VOLU 'TRAN' 1 PZ ;
  36.  
  37. * Elements quadratiques
  38. * =====================
  39. _SEG3  = CHAN 'QUAD' _SEG2 ;
  40. _TRI6  = CHAN 'QUAD' _TRI3 ;
  41. _QUA8  = CHAN 'QUAD' _QUA4 ;
  42. _CUB20 = CHAN 'QUAD' _CUB8 ;
  43. _TET10 = CHAN 'QUAD' _TET4 ;
  44. _PYR13 = CHAN 'QUAD' _PYR5 ;
  45. _PRI15 = CHAN 'QUAD' _PRI6 ;
  46.  
  47. * Elements a bulles
  48. * =================
  49. _TRI7  = CHAN 'QUAF' _TRI3 ;
  50. _QUA9  = CHAN 'QUAF' _QUA4 ;
  51. _CUB27 = CHAN 'QUAF' _CUB8 ;
  52.  
  53. * Concatenation de tous les maillages
  54. * ===================================
  55. _SEG2  = _SEG2  PLUS (2. *PX) ;
  56. _TRI3  = _TRI3  PLUS (4. *PX) ;
  57. _QUA4  = _QUA4  PLUS (6. *PX) ;
  58. _CUB8  = _CUB8  PLUS (8. *PX) ;
  59. _TET4  = _TET4  PLUS (10.*PX) ;
  60. _PYR5  = _PYR5  PLUS (12.*PX) ;
  61. _PRI6  = _PRI6  PLUS (14.*PX) ;
  62. _SEG3  = _SEG3  PLUS (2. *PX) PLUS (2.*PY) ;
  63. _TRI6  = _TRI6  PLUS (4. *PX) PLUS (2.*PY) ;
  64. _QUA8  = _QUA8  PLUS (6. *PX) PLUS (2.*PY) ;
  65. _CUB20 = _CUB20 PLUS (8. *PX) PLUS (2.*PY) ;
  66. _TET10 = _TET10 PLUS (10.*PX) PLUS (2.*PY) ;
  67. _PYR13 = _PYR13 PLUS (12.*PX) PLUS (2.*PY) ;
  68. _PRI15 = _PRI15 PLUS (14.*PX) PLUS (2.*PY) ;
  69. _TRI7  = _TRI7  PLUS (4. *PX) PLUS (4.*PY) ;
  70. _QUA9  = _QUA9  PLUS (6. *PX) PLUS (4.*PY) ;
  71. _CUB27 = _CUB27 PLUS (8. *PX) PLUS (4.*PY) ;
  72.  
  73. MTOT = _POI1 ET
  74.        _SEG2 ET _TRI3 ET _QUA4 ET _CUB8 ET _TET4 ET _PYR5 ET _PRI6 ET
  75.        _SEG3 ET _TRI6 ET _QUA8 ET _CUB20 ET _TET10 ET _PYR13 ET
  76.         _PRI15 ET _TRI7 ET _QUA9 ET _CUB27 ;
  77.  
  78. MTOT1 = CHAN 'POI1' MTOT ;
  79.  
  80.  
  81. * ****************************
  82. * ORDRE DES NOEUDS DANS CAST3M
  83. * ****************************
  84.  
  85. SI GRAPH ;
  86.     REPE BLOC1 (NBNO MTOT1) ;
  87.         TRAC (MTOT ET (COUL 'ROUG' (MANU 'POI1' (MTOT1 POIN &BLOC1)))) ;
  88.     FIN BLOC1 ;
  89. FINS ;
  90.  
  91.  
  92. * *************************
  93. * EXPORTATION VERS PARAVIEW
  94. * *************************
  95.  
  96. REPE BLOC1 (NBNO MTOT1) ;
  97.  
  98.     CHPO1 = MANU 'CHPO' (MTOT1 POIN &BLOC1) 1 'SCAL' 1. ;
  99.  
  100. *    OPTI 'SORT' ('CHAINE' DIVERS '/elements_vtk') ;
  101.      OPTI 'SORT' 'elements_vtk' ;
  102.      SORT 'VTK' MTOT CHPO1 'ID_NOEUD' 'TEMPS' (&BLOC1*0.001) ;
  103.  
  104. FIN BLOC1 ;
  105.  
  106.  
  107. FIN ;
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales