Télécharger INTG_test.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : INTG_test.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5.  
  6. ************************************************************************
  7. * Ce Cas-Test permet de tester l'operateur INTG dans differentes
  8. * configurations d'options
  9. ************************************************************************
  10.  
  11. OPTI DIME 3 ELEM 'CUB8';
  12. XZPREC = (VALE 'PREC') * 100.D0;
  13.  
  14. ************************************************************************
  15. * MAILLAGE d'un CUBE de cote L
  16. ************************************************************************
  17. L = 2.D0 ;
  18. NBE1 = 10 ;
  19.  
  20. P1 = 0. 0. 0. ;
  21. P2 = L 0. 0. ;
  22.  
  23. L1 = DROI NBE1 P1 P2;
  24. L2 = L1 PLUS (0. L 0.);
  25.  
  26. S1 = L1 REGL NBE1 L2;
  27. S2 = S1 PLUS (0. 0. L);
  28.  
  29. V1 = S1 VOLU NBE1 S2;
  30.  
  31. *TRAC CACH V1;
  32.  
  33.  
  34. ************************************************************************
  35. * MODELE et MCHAML
  36. ************************************************************************
  37. MOD1 = MODE V1 'MECANIQUE' 'ELASTIQUE';
  38. MCH1 = MANU 'CHML' MOD1 'VOL1' 1.D0 ; COMM 'MCHAML CONSTANT';
  39.  
  40. MOD2 = REDU MOD1 (VIDE 'MAILLAGE');
  41.  
  42.  
  43. ************************************************************************
  44. * INTEGRALES
  45. ************************************************************************
  46. INTG11 = INTG MOD1 MCH1;
  47.  
  48. TEST = ABS (INTG11 - (L ** 3)) ;
  49. MESS TEST XZPREC;
  50. SI (TEST > XZPREC);
  51. MESS 'ERREUR 1 DANS INTG';
  52. ERRE 5;
  53. FINS;
  54.  
  55.  
  56. INTG12 = INTG MOD1 MCH1 'ELEM';
  57. MIN12 = MINI INTG12;
  58. MAX12 = MAXI INTG12;
  59. DELTA = MAX12 - MIN12;
  60. TEST = MAX12 - (INTG11 / (NBE1 **3));
  61. MESS MIN12 MAX12 DELTA TEST;
  62. SI ((DELTA > XZPREC) OU (TEST > XZPREC));
  63. MESS 'ERREUR 2 DANS INTG';
  64. ERRE 5;
  65. FINS;
  66.  
  67.  
  68. INTG21 = INTG MOD2 MCH1;
  69. MESS INTG21;
  70. SI (INTG21 > 0.D0);
  71. MESS 'ERREUR 3 DANS INTG';
  72. ERRE 5;
  73. FINS;
  74.  
  75.  
  76. INTG22 = INTG MOD2 MCH1 'ELEM';
  77. MAIL22 = EXTR INTG22 'MAIL';
  78. TEST = NBEL MAIL22;
  79. SI (TEST NEG 0);
  80. MESS 'ERREUR 4 DANS INTG';
  81. ERRE 5;
  82. FINS;
  83.  
  84. FIN;
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  

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