Télécharger la plaquette Cast3M


fr en

Télécharger verfdg.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  verfdg.dgibi
  2. *
  3. * Ce test vérifie les matrices masses diagonales
  4. * modèle NAVIER_STOKES
  5. *
  6. opti dime 2 elem QUA8 ;
  7. p1= 0. 0. ; p2 = 1.13 0. ;
  8. l1= p1 d dini 0.2 dfin 0.1 p2;
  9.  
  10. *DISCR='LINE' ;
  11. *DISCR='QUAF' ;
  12. n1= 10 ;
  13. mt=l1 trans n1 (0. 1.) ;
  14. Mmt= chan mt QUAF ;
  15.  
  16. DISCR='MACRO';
  17. $mt=mode Mmt 'NAVIER_STOKES' DISCR ;
  18. Dg= doma $mt 'XXCTREP1' ;
  19. list (somt dg) ;
  20. er = abs ( (somt dg) - 1.13 ) ;
  21. mess ' er= ' er ;
  22. Si (er > 1.e-5) ; erreur 5 ; Finsi ;
  23.  
  24. Dg= doma $mt 'XXCTREP0' ;
  25. list (somt dg) ;
  26. er = abs ( (somt dg) - 1.13 ) ;
  27. mess ' er= ' er ;
  28. Si (er > 1.e-5) ; erreur 5 ; Finsi ;
  29.  
  30. Dg= doma $mt 'XXMSOMME' ;
  31. list (somt dg) ;
  32. er = abs ( (somt dg) - 1.13 ) ;
  33. mess ' er= ' er ;
  34. Si (er > 1.e-5) ; erreur 5 ; Finsi ;
  35.  
  36. DISCR='QUAF';
  37. $mt=mode Mmt 'NAVIER_STOKES' DISCR ;
  38. Dg= doma $mt 'XXCTREP1' ;
  39. list (somt dg) ;
  40. er = abs ( (somt dg) - 1.13 ) ;
  41. mess ' er= ' er ;
  42. Si (er > 1.e-5) ; erreur 5 ; Finsi ;
  43.  
  44. Dg= doma $mt 'XXCTREP0' ;
  45. list (somt dg) ;
  46. er = abs ( (somt dg) - 1.13 ) ;
  47. mess ' er= ' er ;
  48. Si (er > 1.e-5) ; erreur 5 ; Finsi ;
  49.  
  50. Dg= doma $mt 'XXMSOMME' ;
  51. list (somt dg) ;
  52. er = abs ( (somt dg) - 1.13 ) ;
  53. mess ' er= ' er ;
  54. Si (er > 1.e-5) ; erreur 5 ; Finsi ;
  55.  
  56. aa= extr (doma $mt 'XXMSOMME') maillage ;
  57. aa = aa coul rouge ;
  58. mt= doma $mt maillage ;
  59. *trace (MT et aa) ;
  60.  
  61. DISCR='LINE';
  62. $mt=mode Mmt 'NAVIER_STOKES' DISCR ;
  63. Dg= doma $mt 'XXMSOMME' ;
  64. list (somt dg) ;
  65. er = abs ( (somt dg) - 1.13 ) ;
  66. mess ' er= ' er ;
  67. Si (er > 1.e-5) ; erreur 5 ; Finsi ;
  68.  
  69. FIN ;
  70.  
  71.  
  72.  

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