Télécharger la plaquette Cast3M


fr en

Télécharger vecoul2D.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : vecoul2D.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ********************************************************************
  5. 'OPTION' 'ECHO' 1 ;
  6. 'SAUTER' 'PAGE';
  7. *
  8. 'TITRE' 'Vecteurs couleur'   ;
  9. OPTI DIME 2 ELEM QUA4   ;
  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 = 12.0D0 ;
  24. Y1 = 15.0D0 ;
  25. *
  26. A1 = 2.D0 0.0D0 ;
  27. B1 = X1 0.0D0  ;
  28. C1 = X1 Y1 ;
  29. D1 = 0.0D0 Y1 ;
  30. *
  31. *--------------------- Lignes du plan de base --------------
  32. *                          -- Densites --
  33. NIV1 = 1.0 ;
  34. NBX1 = ENTIER (NIV1 * 40) ;
  35. NBY1 = ENTIER (NIV1 * 40) ;
  36. LG1 = DROIT NBX1 A1 B1 ;
  37. LG2 = DROIT NBY1 B1 C1 ;
  38. LG3 = DROIT NBX1 C1 D1 ;
  39. LG4 = DROIT NBY1 D1 A1 ;
  40. *
  41. *----------------------- Surface plan de base --------------
  42. *
  43. MASSIF0 = COUL BLEU (DALLER LG1 LG2 LG3 LG4) ;
  44. *
  45. *=====================================================================
  46. *                             -- MAILLAGES QUAF --
  47. QFTOT = CHANGE MASSIF0 QUAF ;
  48. *
  49. *=====================================================================
  50. *                                 -- MODELE --
  51. MODHYB = MODE QFTOT 'DARCY' 'ANISOTROPE'     ;
  52. XXC YYC = 'COOR' (DOMA MODHYB 'CENTRE') ;
  53. *
  54. *=====================================================================
  55. *                   -- Vecteur de base --
  56. CH1 = (NOMC 'VX' (2.0*XXC)) ET (NOMC 'VY' YYC) ;
  57. CH2 = 1.0D-12 * CH1 ;
  58. VECT1 = VECT CH2 VX VY ROUG ;
  59. SI GRAPH1 ;
  60.    TRAC VECT1 QFTOT ;
  61. FINSI ;
  62. AMPL2 = 0.61677E+11 ;
  63. *
  64. VCTOT1 = @VECOUL CH2 AMPL2 ;
  65. SI GRAPH1 ;
  66.    TRAC VCTOT1 QFTOT ;
  67. FINSI ;
  68. *
  69. VCTOT1 = @VECOUL CH2 AMPL2 0.5 ;
  70. SI GRAPH1 ;
  71.    TRAC VCTOT1 QFTOT ;
  72. FINSI ;
  73. *
  74. VCTOT1 = @VECOUL AMPL2 CH2 'ALTR' 0.5 ;
  75. SI GRAPH1 ;
  76.    TRAC VCTOT1 QFTOT ;
  77. FINSI ;
  78. *
  79. VCTOT1 = @VECOUL CH2 AMPL2 1. ;
  80. SI GRAPH1 ;
  81.    TRAC VCTOT1 QFTOT ;
  82. FINSI ;
  83. *
  84. VCTOT1 = @VECOUL AMPL2 CH2 'ALTR' 1. ;
  85. SI GRAPH1 ;
  86.    TRAC VCTOT1 QFTOT ;
  87. FINSI ;
  88. *
  89. VCTOT1 = @VECOUL CH2 AMPL2 150 ;
  90. SI GRAPH1 ;
  91.    TRAC VCTOT1 QFTOT ;
  92. FINSI ;
  93. *
  94. VCTOT1 = @VECOUL AMPL2 CH2 'ALTR' 150 ;
  95. SI GRAPH1 ;
  96.    TRAC VCTOT1 QFTOT ;
  97. FINSI ;
  98. *
  99. VCTOT1 = @VECOUL CH2 AMPL2 2 ;
  100. SI GRAPH1 ;
  101.    TRAC VCTOT1 QFTOT ;
  102. FINSI ;
  103. *
  104. VCTOT1 = @VECOUL AMPL2 CH2 'ALTR' 2 ;
  105. SI GRAPH1 ;
  106.    TRAC VCTOT1 QFTOT ;
  107. FINSI ;
  108. *
  109. FIN ;
  110. *
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

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