Télécharger test_norm_env.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : test_norm_env.dgibi
  2.  
  3. GRAPH=VRAI ;
  4. GRAPH=FAUX ;
  5. OPTI ECHO 0 ;
  6.  
  7. opti dime 2 elem tri6 ;
  8. p1=1 0; p2=0 0 ;
  9. l1= p1 d 1 p2 ;
  10. m1= l1 trans 1 (0 1) ;
  11. opti dime 2 elem qua8 ;
  12. m2=(cote 3 m1) trans 1 (0 1) ;
  13. m1 = m1 et m2 ;
  14. Si GRAPH ;trace m1 ;Finsi;
  15.  
  16. Mm1= chan m1 QUAF ;
  17. ERROR = 0 ;
  18. * ** 2D ***
  19.  
  20. ******************** LINE *******************************
  21. DISCR='LINE' ;
  22. $m1=mode Mm1 'NAVIER_STOKES' DISCR ;
  23. an= doma $m1 'NORMALEV' ;
  24. an1= vect an 0.1 ux uy jaune;
  25. Si GRAPH ; trace an1 (DOMA $m1 'MAILLAGE') ; Finsi;
  26. *list an ;
  27.  
  28. 'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
  29. DB1 = 'DBIT' ('MODE' ('CONT' Mm1) 'NAVIER_STOKES' DISCR) an 'IMPR' ;
  30. DB2 = 4.82842712 ;
  31. err1=1.E-8 ;
  32. Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;
  33.  
  34. ******************** QUAF *******************************
  35. DISCR='QUAF' ;
  36. $m1=mode Mm1 'NAVIER_STOKES' DISCR ;
  37. an= doma $m1 'NORMALEV' ;
  38. an1= vect an 0.1 ux uy jaune;
  39. Si GRAPH ; trace an1 (DOMA $m1 'MAILLAGE') ; Finsi;
  40. *list an ;
  41.  
  42. 'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
  43. DB1 = 'DBIT' ('MODE' ('CONT' Mm1) 'NAVIER_STOKES' DISCR) an 'IMPR' ;
  44. DB2 = 5.60947571 ;
  45. err1=1.E-8 ;
  46. Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;
  47.  
  48. ******************** MACRO ******************************
  49. DISCR='MACRO';
  50. $m1=mode Mm1 'NAVIER_STOKES' DISCR ;
  51. an= doma $m1 'NORMALEV' ;
  52. an1= vect an 0.1 ux uy jaune;
  53. Si GRAPH ; trace an1 (DOMA $m1 'MAILLAGE') ; Finsi;
  54. *list an ;
  55.  
  56. 'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
  57. DB1 = 'DBIT' ('MODE' ('CONT' Mm1) 'NAVIER_STOKES' DISCR) an 'IMPR' ;
  58. DB2 = 5.41421356 ;
  59. err1=1.E-8 ;
  60. Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;
  61.  
  62.  
  63. **********************************
  64. * ** 3D ***
  65. **********************************
  66.  
  67. opti dime 3 elem cu20 ;
  68. v1= m1 volu 1 trans (0 0 1) ;
  69. Si GRAPH; trace v1 ; Finsi;
  70. Mv1= chan QUAF v1 ;
  71.  
  72. ******************** LINE *******************************
  73. DISCR='LINE' ;
  74. $v1= mode Mv1 'NAVIER_STOKES' DISCR ;
  75. na= doma $v1 'NORMALEV' ;
  76. na1= vect na 0.1 ux uy uz jaune ;
  77. Si GRAPH ; trace na1 (DOMA $v1 'MAILLAGE') ; Finsi;
  78. *list an ;
  79.  
  80. 'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
  81. DB1 = 'DBIT' ('MODE' ('DOMA' $v1 'ENVELOPP') 'NAVIER_STOKES' DISCR)
  82. na 'IMPR' ;
  83. DB2 = 6.32391111 ;
  84. err1=1.E-8 ;
  85. Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;
  86.  
  87. ******************** QUAF *******************************
  88. DISCR='QUAF' ;
  89. $v1= mode Mv1 'NAVIER_STOKES' DISCR ;
  90. na= doma $v1 'NORMALEV' ;
  91. na1= vect na 0.1 ux uy uz jaune ;
  92. Si GRAPH ; trace na1 (DOMA $v1 'MAILLAGE') ; Finsi;
  93. *list an ;
  94.  
  95. 'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
  96. DB1 = 'DBIT' ('MODE' ('DOMA' $v1 'ENVELOPP') 'NAVIER_STOKES' DISCR)
  97. na 'IMPR' ;
  98. DB2 = 8.62820114 ;
  99. err1=1.E-8 ;
  100. Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;
  101.  
  102. ******************** MACRO ******************************
  103. DISCR='MACRO' ;
  104. $v1= mode Mv1 'NAVIER_STOKES' DISCR ;
  105. na= doma $v1 'NORMALEV' ;
  106. na1= vect na 0.1 ux uy uz jaune ;
  107. Si GRAPH ; trace na1 (DOMA $v1 'MAILLAGE') ; Finsi;
  108. *list an ;
  109.  
  110. 'MESS' 'Dimension ' ('VALE' 'DIME') ' Type ' DISCR ;
  111. DB1 = 'DBIT' ('MODE' ('DOMA' $v1 'ENVELOPP') 'NAVIER_STOKES' DISCR)
  112. na 'IMPR' ;
  113. DB2 = 7.87130633 ;
  114. err1=1.E-8 ;
  115. Si ( (abs (DB1 - DB2)) > err1 ) ; ERROR = ERROR + 1 ; finsi ;
  116.  
  117. 'SI' (ERROR > 0) ;
  118. 'MESS' 'Il y a ' ERROR ' erreur(s)' ;
  119. 'ERRE' 5 ;
  120. 'SINON' ;
  121. 'MESS' 'Tout s est bien passé !' ;
  122. 'FINSI' ;
  123.  
  124.  
  125. FIN ;
  126.  
  127.  
  128.  
  129.  
  130.  

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