Télécharger elements_vtk.dgibi

Retour à la liste

Numérotation des lignes :

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

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