Télécharger la plaquette Cast3M


fr en

Télécharger testlgQUAF.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  testlgQUAF.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. * On teste le bon fonctionnement de l'option ARETE de DOMA.
  6. * Ceci necessite la bonne description des lignes et des faces
  7. * des elements QUAF
  8. *
  9. * EN 2D  QUA9 et TRI7
  10. opti dime 2 elem QUA8;
  11.  
  12. Graph = FAUX ;
  13.  
  14. DISCR= QUAF ;
  15. p0 = 0. 0. ;
  16. p1 = 1. 0. ;
  17. AB = p0 d 4  p1 ;
  18. mt = ab trans 4  (0. 2.) ;
  19. cd = cote 3 mt ;
  20. opti elem tri6 ;
  21. mtt = cd trans 2 (1. 1.);
  22. mt = mt et mtt  ;
  23. Mmt = chan mt QUAF ;
  24. $mt = mode Mmt 'NAVIER_STOKES' DISCR;
  25.  
  26. arm = doma $mt 'ARETE';
  27. lt = elem arm type;
  28. Si (NON (EGA (dime lt) 1));
  29. ERREUR 5;
  30. FINSI;
  31. A= extr lt 1 ;
  32. Si (NON (EGA a 'SEG3'));
  33. ERREUR 5;
  34. FINSI;
  35.  
  36. Si Graph;
  37. xv= coor 1 Mmt ;
  38. trace xv Mmt ;
  39. Finsi ;
  40.  
  41. * EN 3D  CU27 et PR21
  42. opti dime 3 ELEM cu20 ;
  43.  
  44. vt = mt volu 8 trans (0. 0. 2.);
  45. Mvt = chan QUAF vt ;
  46. $vt = mode Mvt 'NAVIER_STOKES' DISCR;
  47. arm = doma $vt 'ARETE';
  48.  
  49. lt = elem arm type;
  50. Si (NON (EGA (dime lt) 1));
  51. ERREUR 5;
  52. FINSI;
  53. A= extr lt 1 ;
  54. Si (NON (EGA a 'SEG3'));
  55. ERREUR 5;
  56. FINSI;
  57.  
  58. Si Graph;
  59. xv= (coor 1 Mvt) ;
  60. trace xv Mvt ;
  61. Finsi ;
  62.  
  63. * EN 3D  les 4 elements CU27 PR21 TE15 et PY19
  64.  ae = enve vt ;
  65.  u  = volu ae ;
  66.  
  67.  Mu = chan u QUAF ;
  68. le = elem Mu type ;
  69.  
  70. Te=( (EXIST le 'CU27') et (EXIST le 'PR21') et
  71.      (EXIST le 'TE15') et (EXIST le 'PY19'));
  72. Si (NON TE);
  73. ERREUR 5;
  74. FINSI;
  75.  
  76.  ul = chan ligne Mu ;
  77.  
  78. lt = elem ul  type;
  79. Si (NON (EGA (dime lt) 1));
  80. ERREUR 5;
  81. FINSI;
  82. A= extr lt 1 ;
  83. Si (NON (EGA a 'SEG3'));
  84. ERREUR 5;
  85. FINSI;
  86.  
  87.  
  88. Si Graph;
  89. xv =  (coor 1 Mu);
  90. trace xv Mu;
  91. Finsi ;
  92. Fin ;
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  

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