Télécharger la plaquette Cast3M


fr en

Télécharger tliqu.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  tliqu.dgibi
  2. *****************  CAS TEST : tliqu.dgibi   **************************
  3. *
  4. * Test de l'opérateur VARI TLIQUID(P,H)
  5. * Les données sont un FLOTTANT, un LISTREEL ou un CHPO
  6. *
  7. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  8. 'OPTI' 'ECHO' 0 ;
  9. *
  10. *------------------------------------------------------> FLOTTANT
  11. *
  12. X1 = 1.0133D5 ;
  13. X2 = 4.1906D5 ;
  14. X3 = VARI TLIQUID X1 X2 ;
  15. XSTO1 = X3 ;
  16. *
  17. X1 = 5.021D5  ;
  18. X2 = 6.4079D5 ;
  19. X3 = VARI TLIQUID X1 X2 ;
  20. XSTO2 = X3 ;
  21. *
  22. X1 = 10.027D5 ;
  23. X2 = 7.6312D5 ;
  24. X3 = VARI TLIQUID x1 x2 ;
  25. XSTO3 = X3 ;
  26. *------------------------------------------------------> LISTREEL
  27. Y1 = 'PROG' 1.0133D5 5.021D5  10.027D5 ;
  28. Y2 = 'PROG' 4.1906D5 6.4079D5 7.6312D5 ;
  29. Y3 = VARI TLIQUID Y1 Y2 ;
  30. YSTO1 = EXTR Y3 1 ;
  31. YSTO2 = EXTR Y3 2 ;
  32. YSTO3 = EXTR Y3 3 ;
  33. *------------------------------------------------------> CHPO
  34. P1   = 1.0133 4.1906 ;
  35. P2   = 5.021 6.4079 ;
  36. P3   = 10.027 7.6312 ;
  37. P1P3 = P1 'DROI' 1 P2 'DROI' 1 P3 ;
  38. Z1   = 1.D5 * ('COOR' 1 P1P3) ;
  39. Z2   = 1.D5 * ('COOR' 2 P1P3) ;
  40. Z3   = VARI TLIQUID Z1 Z2 ;
  41. ZSTO1 = 'EXTR' Z3 'SCAL' P1  ;
  42. ZSTO2 = 'EXTR' Z3 'SCAL' P2  ;
  43. ZSTO3 = 'EXTR' Z3 'SCAL' P3  ;
  44. *-------------------------------------------> Controle
  45. CTRL1 = XSTO1 + YSTO1 + ZSTO1 / 3.d0 - XSTO1 / XSTO1 ;
  46. CTRL2 = XSTO2 + YSTO2 + ZSTO2 / 3.d0 - XSTO2 / XSTO2 ;
  47. CTRL3 = XSTO3 + YSTO3 + ZSTO3 / 3.d0 - XSTO3 / XSTO3 ;
  48. CTRL4 = CTRL1 + CTRL2 + CTRL3 ;
  49. XREF1 = 373.15 ;
  50. XREF2 = 425.15 ;
  51. XREF3 = 453.15 ;
  52. CTRL5 = XSTO1 + YSTO1 + ZSTO1 / 3.d0 - XREF1 / XREF1 ;
  53. CTRL6 = XSTO2 + YSTO2 + ZSTO2 / 3.d0 - XREF2 / XREF2 ;
  54. CTRL7 = XSTO3 + YSTO3 + ZSTO3 / 3.d0 - XREF3 / XREF3 ;
  55. CTRL8 = CTRL5 + CTRL6 + CTRL7 ;
  56. *-------------------------------------------> Affichage
  57. 'MESS' '     ' ;
  58. 'MESS' '     ' ;
  59. 'MESS' 'Point 1  VARI TLIQUID ----->' XSTO1 YSTO1 ZSTO1 XREF1 ;
  60. 'MESS' 'Point 2  VARI TLIQUID ----->' XSTO2 YSTO2 ZSTO2 XREF2 ;
  61. 'MESS' 'Point 3  VARI TLIQUID ----->' XSTO3 YSTO3 ZSTO3 XREF3 ;
  62. 'MESS' '     ' ;
  63. 'MESS' 'Comparaison calculs rel.----->' CTRL4 ;
  64. 'MESS' 'Erreur rel. VDI         ----->' CTRL8 ;
  65. 'MESS' '     ' ;
  66. 'MESS' '     ' ;
  67. 'MESS' '     ' ;
  68. 'MESS' '     ' ;
  69. 'MESS' '     ' ;
  70. *-------------------------------------------> Compte-rendu et sortie
  71. EPS4  = 1.E-14 ;
  72. EPS8  = 1.2E-3  ;
  73. CTRL4 = ABS CTRL4 ;
  74. CTRL8 = ABS CTRL8 ;
  75. LOG4 = CTRL4 > EPS4 ;
  76. LOG8 = CTRL8 > EPS8 ;
  77. L0   = LOG4 'OU' LOG8 ;
  78. 'SI' L0        ;
  79.    'ERREUR' 5  ;
  80. 'SINON'        ;
  81.    'ERREUR' 0  ;
  82. 'FINSI'        ;
  83. 'FIN' ;
  84.  
  85.  
  86.  

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