Télécharger INTG_test_integration_reduite.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. ************************************************************************
  3. * Cas-test de calculs d'integrales avec les elements a integration
  4. * reduite C20R et P15R
  5. ************************************************************************
  6.  
  7. OPTI 'DIME' 3 'ELEM' 'CU20';
  8. XZPREC = (VALE 'PREC') * 100.;
  9.  
  10. **************
  11. * CUBE C20R
  12. **************
  13.  
  14. P1 = 0. 0. 0.;
  15. P2 = 1. 0. 0.;
  16. D1 = DROI 1 P1 P2;
  17. S1 = D1 TRAN 1 (0. 1. 0.);
  18. V1 = S1 VOLU 'TRAN' 1 (0. 0. 1.);
  19. * TRAC V1;FIN;
  20.  
  21. MO = MODE V1 'MECANIQUE' 'ELASTIQUE' 'C20R';
  22.  
  23. * VOLUME A PARTIR D’UN CHAMP CONSTANT
  24. UN = MANU 'CHML' MO 'SCAL' 1.;
  25. VOL = INTG MO UN;
  26. ECART = ABS (VOL - 1.);
  27. MESS ECART XZPREC;
  28. SI (ECART > XZPREC);
  29. MESS 'ERREUR DANS INTG POUR UN CHAMP CONSTANT';
  30. ERRE 5;
  31. FINSI;
  32.  
  33. * CHAMP LINEAIRE
  34. Z = COOR 3 UN;
  35. RES = INTG MO Z;
  36. ECART = ABS (RES - 0.5);
  37. MESS ECART XZPREC;
  38. SI (ECART > XZPREC);
  39. MESS 'ERREUR DANS INTG POUR UN CHAMP LINEAIRE';
  40. ERRE 5;
  41. FINSI;
  42.  
  43. * CHAMP QUADRATIQUE
  44. RES = INTG MO (Z**2);
  45. ECART = ABS (RES - (1./3.));
  46. MESS ECART XZPREC;
  47. SI (ECART > XZPREC);
  48. MESS 'ERREUR DANS INTG POUR UN CHAMP QUADRATIQUE';
  49. ERRE 5;
  50. FINSI;
  51.  
  52. **************
  53. * PRISME P15R
  54. **************
  55.  
  56. P1 = 0. 0. 0.;
  57. P2 = 1. 0. 0.;
  58. P3 = 0. 1. 0.;
  59. D1 = DROI 1 (DROI 1 (DROI 1 P1 P2) P3) P1;
  60. S1 = SURF 'PLAN' D1;
  61. V1 = S1 VOLU 'TRAN' 1 (0. 0. 1.);
  62. * TRAC V1;FIN;
  63.  
  64. MO = MODE V1 'MECANIQUE' 'ELASTIQUE' 'P15R';
  65.  
  66. * VOLUME A PARTIR D’UN CHAMP CONSTANT
  67. UN = MANU 'CHML' MO 'SCAL' 1.;
  68. VOL = INTG MO UN;
  69. ECART = ABS (VOL - 0.5);
  70. MESS ECART XZPREC;
  71. SI (ECART > XZPREC);
  72. MESS 'ERREUR DANS INTG POUR UN CHAMP CONSTANT';
  73. ERRE 5;
  74. FINSI;
  75.  
  76. * CHAMP LINEAIRE
  77. Z = COOR 3 UN;
  78. RES = INTG MO Z;
  79. ECART = ABS (RES - 0.25);
  80. MESS ECART XZPREC;
  81. SI (ECART > XZPREC);
  82. MESS 'ERREUR DANS INTG POUR UN CHAMP LINEAIRE';
  83. ERRE 5;
  84. FINSI;
  85.  
  86. * CHAMP QUADRATIQUE
  87. RES = INTG MO (Z**2);
  88. ECART = ABS (RES - (1./6.));
  89. MESS ECART XZPREC;
  90. SI (ECART > XZPREC);
  91. MESS 'ERREUR DANS INTG POUR UN CHAMP QUADRATIQUE';
  92. ERRE 5;
  93. FINSI;
  94.  
  95. FIN;
  96.  
  97.  
  98.  

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