Télécharger la plaquette Cast3M


fr en

Télécharger vecoul3D.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : vecoul3D.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *********************************************************************
  5. 'OPTION' 'ECHO' 1 ;
  6. 'SAUTER' 'PAGE';
  7. *
  8. 'TITRE' 'Vecteurs couleur'   ;
  9. OPTI DIME 3 ELEM CUB8   ;
  10. *OPTI ISOV SURFACE ;
  11. OPTI ISOV SULI ;
  12. *OPTI TRAC PSC ;
  13. *
  14. GRAPH1 = FAUX ;
  15. ************************ CAS TEST : vecoul.dgibi ******************
  16. *******************************************************************
  17. *
  18. *--------------------- Création du maillage 3D ---------------------
  19. *
  20. EPSI1 = 0.0001 ;
  21. *--------------------- Points du plan de base --------------
  22. *
  23. X1 = 10.0D0 ;
  24. Y1 = 9.0D0 ;
  25. Z1 = 13.5D0 ;
  26. *
  27. A1 = 2D0 0.0D0 0.0D0 ;
  28. B1 = X1 0.0D0 0.0D0 ;
  29. C1 = X1 Y1 0.0D0 ;
  30. D1 = 0.0D0 Y1 0.0D0 ;
  31. *
  32. *--------------------- Lignes du plan de base --------------
  33. *                          -- Densites --
  34. NIV1 = 1.0 ;
  35. NBX1 = ENTIER (NIV1 * 5) ;
  36. NBY1 = ENTIER (NIV1 * 8) ;
  37. NBZ1 = ENTIER (NIV1 * 6) ;
  38. LG1 = DROIT NBX1 A1 B1 ;
  39. LG2 = DROIT NBY1 B1 C1 ;
  40. LG3 = DROIT NBX1 C1 D1 ;
  41. LG4 = DROIT NBY1 D1 A1 ;
  42. *
  43. *----------------------- Surface plan de base --------------
  44. *
  45. GP1 = DALLER LG1 LG2 LG3 LG4 ;
  46. *
  47. *------------------------------- Volume --------------------
  48. *
  49. VD1 = 0.0D0 0.0D0 Z1 ;
  50. GP2 = GP1 PLUS VD1 ;
  51. MASSIF0 = COUL BLEU (GP1 VOLU NBZ1 GP2) ;
  52. SI GRAPH1 ;
  53.    TRAC CACH MASSIF0 ;
  54. FINSI ;
  55. *
  56. *
  57. *=====================================================================
  58. *                             -- MAILLAGES QUAF --
  59. QFTOT = CHANGE MASSIF0 QUAF ;
  60. *
  61. *=====================================================================
  62. *                                 -- MODELE --
  63. MODHYB = MODE QFTOT 'DARCY' 'ANISOTROPE'     ;
  64. XXC YYC ZZC = 'COOR' (DOMA MODHYB 'CENTRE') ;
  65. *
  66. *=====================================================================
  67. *                   -- Vecteur de base --
  68. CH1 = (NOMC 'VX' (2.0*XXC)) ET (NOMC 'VY' YYC) ET (NOMC 'VZ' ZZC) ;
  69. CH2 = 1.0D+11 * CH1 ;
  70. VECT1 = VECT CH2 VX VY VZ ROUG ;
  71. AFTOT = 'ARETE' QFTOT ;
  72. SI GRAPH1 ;
  73.    TRAC VECT1 QFTOT AFTOT ;
  74. FINSI ;
  75. AMPL2 = 0.61224E-12 ;
  76. *
  77. VCTOT1 = @VECOUL AMPL2 CH2 ;
  78. SI GRAPH1 ;
  79.    TRAC VCTOT1 QFTOT AFTOT ;
  80. FINSI ;
  81. *
  82. VCTOT1 = @VECOUL AMPL2 CH2 0.34 ;
  83. SI GRAPH1 ;
  84.    TRAC VCTOT1 QFTOT AFTOT ;
  85. FINSI ;
  86. *
  87. VCTOT1 = @VECOUL AMPL2 CH2 'ALTR' 0.34 ;
  88. SI GRAPH1 ;
  89.    TRAC VCTOT1 QFTOT AFTOT ;
  90. FINSI ;
  91. *
  92. LMOT = 'MOTS' 'VX' 'VY' 'VZ' ;
  93. *
  94. VCTOT1 = @VECOUL AMPL2 ('*' CH2 0.) LMOT -1. ;
  95. SI GRAPH1 ;
  96.    TRAC VCTOT1 QFTOT AFTOT ;
  97. FINSI ;
  98. *
  99. FIN ;
  100. *
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  

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