Télécharger la plaquette Cast3M


fr en

Télécharger muliq.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : muliq.dgibi
  2. *
  3. *
  4. * Viscosité dynamique de l'eau
  5. * Valeurs de référence : tables VDI
  6. * Comparaisons de différentes formules et tests de 'VARI' 'MULIQ'
  7. *
  8. *
  9. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' 'TRAC' 'PSC' 'EPTR' 10 'DENS' 1. ;
  10. graph = faux ;
  11. *
  12. *- Calcul de la viscosité dynamique de l'eau (formule Cathare)
  13. *
  14. 'DEBP' muliqk TW0*'LISTREEL' ;
  15. *
  16. UN  = TW0 / TW0 ;
  17. *TW1 = 273.15 * UN + TW0 ;
  18. *
  19. a0 = 570.58058 * UN ;
  20. b0 = -140.D0 * UN ;
  21. L1 = 2.414D-5 * ('EXP' (a0 / (b0 + TW0))) ;
  22. *
  23. 'FINP' L1 ;
  24. *
  25. *- Calcul de la viscosité dynamique de l'eau (formule VDI)
  26. *
  27. 'DEBP' muliqv TW0*'LISTREEL' PW0*'FLOTTANT' ;
  28. *
  29. UN  = TW0 / TW0 ;
  30. *TW1 = 273.15 * UN + TW0 ;
  31. PW1 = PW0 * UN ;
  32. *
  33. a1 = 241.4 ;
  34. a2 = 0.3828209486 ;
  35. a3 = 0.2162830218 ;
  36. a4 = 0.1498693949 ;
  37. a5 = 0.4711880117 ;
  38. *
  39. PC0   = 221.2D5 * UN ;
  40. TC0   = (374.15 + 273.15) * UN ;
  41. BETA  = PW1 / PC0 ;
  42. TETA  = TW0 / TC0 ;
  43. PS1   = 'VARI' 'PSATT' TW0 ;
  44. BETAS = PS1 / PC0 ;
  45. *
  46. L1 = (a4 * (BETA - BETAS) * (TETA - (a5 * UN)) + UN) * a1 * 1.D-7 ;
  47. L2 = ('EXP' ( ((a2*UN) / (teta - (a3*un))) * ('LOG' 10.) ) ) * L1 ; 
  48. *
  49. 'FINP' L2 ;
  50. *
  51. *
  52. *----------------------------------------------
  53. *- Température et pression pour comparaison
  54. *----------------------------------------------
  55. *
  56. *
  57. TW1   = ('PROG' 20. 40. 60. 80. 90.) + 273.15 ;
  58. PW1   = 1.D5 ;
  59. LREF1 = 1.D-6 * ('PROG' 1000. 651. 463. 351. 311.) ;
  60. L10 = MULIQK TW1 ;
  61. L11 = MULIQV TW1 PW1 ;
  62. *
  63. TW2   = ('PROG' 10. 30. 50. 100. 130. 150. 170.) + 273.15 ;
  64. PW2   = 10.D5 ;
  65. LREF2 = 1.D-6 * ('PROG' 1300. 797. 544. 279. 211. 181. 159.) ;
  66. L20 = MULIQK TW2 ;
  67. L21 = MULIQV TW2 PW2 ;
  68. *
  69. TW3   = ('PROG' 10. 50. 100. 150. 200. 220.) + 273.15 ;
  70. PW3   = 25.D5 ;
  71. LREF3 = 1.D-6 * ('PROG' 1300. 544. 279. 182. 134. 122.) ;
  72. L30 = MULIQK TW3 ;
  73. L31 = MULIQV TW3 PW3 ;
  74. *
  75. * Test de VARI avec T sous forme de FLOTTANT
  76. 'DEBP' testf t*'LISTREEL' ;
  77.    nb     = 'DIME' t ;
  78.    listmu = 'PROG' ;
  79.    'REPE' blo1 nb ;
  80.       valt   = 'EXTR' t &blo1 ;
  81.       valmu  = 'VARI' 'MULIQ' valt ; 
  82.       listmu = listmu 'ET' ('PROG' valmu) ;
  83.    'FIN' blo1 ;
  84. 'FINP' listmu ;
  85. test_f1 = testf tw1 ;
  86. test_f2 = testf tw2 ;
  87. test_f3 = testf tw3 ;
  88. *
  89. * Test de VARI avec T sous forme de LISTREEL
  90. test_l1 = vari 'MULIQ' tw1 ;
  91. test_l2 = vari 'MULIQ' tw2 ;
  92. test_l3 = vari 'MULIQ' tw3 ;
  93. *
  94. * Test de VARI avec T sous forme de CHPO
  95. 'DEBP' testc t*'CHPOINT' mail1*'MAILLAGE' ;
  96.    test1  = 'VARI' 'MULIQ' t ;
  97.    listmu = 'EXTR' test1 'VALE' mail1 ;
  98. 'FINP' listmu ;
  99. P1 = 1. 0. ; P2 = 2. 0. ;
  100. GEO1 = P1 'DROI' 4 P2 ; PT1 = 'CHAN' 'POI1' GEO1 ;
  101. GEO2 = P1 'DROI' 6 P2 ; PT2 = 'CHAN' 'POI1' GEO2 ;
  102. GEO3 = P1 'DROI' 5 P2 ; PT3 = 'CHAN' 'POI1' GEO3 ;
  103. T1 = 'MANU' 'CHPO' PT1 'SCAL' TW1 ; 
  104. T2 = 'MANU' 'CHPO' PT2 'SCAL' TW2 ; 
  105. T3 = 'MANU' 'CHPO' PT3 'SCAL' TW3 ; 
  106. test_c1 = testc t1 pt1 ;
  107. test_c2 = testc t2 pt2 ;
  108. test_c3 = testc t3 pt3 ;
  109. *
  110. * Comparaisons mise en oeuvre Esope/Gibi pour Cathare 
  111. d10 = 'EVOL'        'MANU' tw1 (l10 - test_L1 / l10 * 100.) ;
  112. d20 = 'EVOL' 'VERT' 'MANU' tw2 (l20 - test_l2 / l20 * 100.) ;
  113. d30 = 'EVOL' 'ROUG' 'MANU' tw3 (l30 - test_l3 / l30 * 100.) ;
  114. *
  115. *-------------------------------------
  116. *- Sorties Graphiques
  117. *-------------------------------------
  118. *
  119. 'SI' graph ;
  120. *
  121. ev10   = 'EVOL'        'MANU' tw1 l10 ;
  122. ev20   = 'EVOL' 'VERT' 'MANU' tw2 l20 ;
  123. ev30   = 'EVOL' 'ROUG' 'MANU' tw3 l30 ;
  124. ev11   = 'EVOL'        'MANU' tw1 l11 ;
  125. ev21   = 'EVOL' 'VERT' 'MANU' tw2 l21 ;
  126. ev31   = 'EVOL' 'ROUG' 'MANU' tw3 l31 ;
  127. evr1   = 'EVOL'        'MANU' tw1 lref1 ;
  128. evr2   = 'EVOL' 'VERT' 'MANU' tw2 lref2 ;
  129. evr3   = 'EVOL' 'ROUG' 'MANU' tw3 lref3 ;
  130. rela10 = 'EVOL'        'MANU' tw1 (l10 - lref1 / lref1 * 100.) ;
  131. rela20 = 'EVOL' 'VERT' 'MANU' tw2 (l20 - lref2 / lref2 * 100.) ;
  132. rela30 = 'EVOL' 'ROUG' 'MANU' tw3 (l30 - lref3 / lref3 * 100.) ;
  133. rela11 = 'EVOL'        'MANU' tw1 (l11 - lref1 / lref1 * 100.) ;
  134. rela21 = 'EVOL' 'VERT' 'MANU' tw2 (l21 - lref2 / lref2 * 100.) ;
  135. rela31 = 'EVOL' 'ROUG' 'MANU' tw3 (l31 - lref3 / lref3 * 100.) ;
  136. *
  137. TAB1 = TABLE ;
  138. TAB1 . 'TITRE' = TABLE ;
  139. TAB1 .  1 = ' MARQ PLEIN ROND NOLI' ;
  140. TAB1 .  2 = ' MARQ PLEIN CARR NOLI' ;
  141. TAB1 .  3 = ' MARQ PLEIN LOSA NOLI' ;
  142. TAB1 .  4 = ' MARQ ROND NOLI' ;
  143. TAB1 .  5 = ' MARQ CARR NOLI' ;
  144. TAB1 .  6 = ' MARQ LOSA NOLI' ;
  145. *
  146. TIT1 = 'Dynamic viscosity of water (Cathare)' ;
  147. TIT2 = ' C' ;
  148. TIT3 = 'kg/m/s' ;
  149. TAB1 . 'TITRE' .  1 = 'MOT' ' 1 bar' ;
  150. TAB1 . 'TITRE' .  2 = 'MOT' '10 bar' ;
  151. TAB1 . 'TITRE' .  3 = 'MOT' '25 bar' ;
  152. TAB1 . 'TITRE' .  4 = 'MOT' 'ref  1' ;
  153. TAB1 . 'TITRE' .  5 = 'MOT' 'ref 10' ;
  154. TAB1 . 'TITRE' .  6 = 'MOT' 'ref 25' ;
  155. 'DESS' (ev10 'ET' ev20 'ET' ev30 'ET' evr1 'ET' evr2 'ET' evr3)
  156. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  157. *
  158. TIT1 = 'Dynamic viscosity of water (VDI)' ;
  159. TIT2 = ' C' ;
  160. TIT3 = 'kg/m/s' ;
  161. TAB1 . 'TITRE' .  1 = 'MOT' ' 1 bar' ;
  162. TAB1 . 'TITRE' .  2 = 'MOT' '10 bar' ;
  163. TAB1 . 'TITRE' .  3 = 'MOT' '25 bar' ;
  164. TAB1 . 'TITRE' .  4 = 'MOT' 'ref  1' ;
  165. TAB1 . 'TITRE' .  5 = 'MOT' 'ref 10' ;
  166. TAB1 . 'TITRE' .  6 = 'MOT' 'ref 25' ;
  167. 'DESS' (ev11 'ET' ev21 'ET' ev31 'ET' evr1 'ET' evr2 'ET' evr3)
  168. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  169. *
  170. TIT1 = 'Dynamic viscosity of water. Relative error of Cathare (x100).' ;
  171. TIT2 = ' C' ;
  172. TIT3 = ' %' ;
  173. TAB1 . 'TITRE' .  1 = 'MOT' 'Cat.  (1 bar)' ;
  174. TAB1 . 'TITRE' .  2 = 'MOT' 'Cat. (10 bar)' ;
  175. TAB1 . 'TITRE' .  3 = 'MOT' 'Cat. (25 bar)' ;
  176. 'DESS' (rela10 'ET' rela20 'ET' rela30)
  177. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  178. *
  179. TIT1 = 'Dynamic viscosity of water. Relative error of VDI (x100).' ;
  180. TIT2 = ' C' ;
  181. TIT3 = ' %' ;
  182. TAB1 . 'TITRE' .  1 = 'MOT' 'VDI.  (1 bar)' ;
  183. TAB1 . 'TITRE' .  2 = 'MOT' 'VDI. (10 bar)' ;
  184. TAB1 . 'TITRE' .  3 = 'MOT' 'VDI. (25 bar)' ;
  185. 'DESS' (rela11 'ET' rela21 'ET' rela31)
  186. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  187. *
  188. TIT1 = 'Dynamic viscosity of water. Relative error (x100).' ;
  189. TIT2 = ' C' ;
  190. TIT3 = ' %' ;
  191. TAB1 . 'TITRE' .  1 = 'MOT' 'Cat.  (1 bar)' ;
  192. TAB1 . 'TITRE' .  2 = 'MOT' 'Cat. (10 bar)' ;
  193. TAB1 . 'TITRE' .  3 = 'MOT' 'Cat. (25 bar)' ;
  194. TAB1 . 'TITRE' .  4 = 'MOT' 'VDI.  (1 bar)' ;
  195. TAB1 . 'TITRE' .  5 = 'MOT' 'VDI. (10 bar)' ;
  196. TAB1 . 'TITRE' .  6 = 'MOT' 'VDI. (25 bar)' ;
  197. 'DESS' (rela10 'ET' rela20 'ET' rela30 'ET'
  198.         rela11 'ET' rela21 'ET' rela31)
  199. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  200. *
  201. TIT1 =
  202. 'Esope/Gibi relative error for Cathare formula (x100)' ;
  203. TIT2 = ' C' ;
  204. TIT3 = ' %' ;
  205. TAB1 . 'TITRE' .  1 = 'MOT' 'Esope/Gibi (1 bar)' ;
  206. TAB1 . 'TITRE' .  2 = 'MOT' 'Esope/Gibi (10 bar)' ;
  207. TAB1 . 'TITRE' .  3 = 'MOT' 'Esope/Gibi (25 bar)' ;
  208. 'DESS' (d10 'ET' d20 'ET' d30)
  209. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  210. 'FINSI' ;
  211. *
  212. *----------------------------------------------------------------------
  213. *- Tests de bon fonctionnement
  214. *- (Comparaisons des différents types et de la mise en oeuvre en gibi)
  215. *----------------------------------------------------------------------
  216. *
  217. diff1 = 'MAXI' (test_f1 - test_l1 / l10) 'ABS' ;
  218. diff2 = 'MAXI' (test_f2 - test_l2 / l20) 'ABS' ;
  219. diff3 = 'MAXI' (test_f3 - test_l3 / l30) 'ABS' ;
  220. diff4 = 'MAXI' (test_f1 - test_c1 / l10) 'ABS' ;
  221. diff5 = 'MAXI' (test_f2 - test_c2 / l20) 'ABS' ;
  222. diff6 = 'MAXI' (test_f3 - test_c3 / l30) 'ABS' ;
  223. diff7 = 'MAXI' (test_f1 - l10 / l10)     'ABS' ;
  224. diff8 = 'MAXI' (test_f2 - l20 / l20)     'ABS' ;
  225. diff9 = 'MAXI' (test_f3 - l30 / l30)     'ABS' ;
  226. didi  = diff1 + diff2 + diff3 + diff4 + diff5
  227.       + diff6 + diff7 + diff8 + diff9 ;
  228. 'MESS' 'Sommes du max des différences relatives ' didi ;
  229. 'SI' (didi > 1.D-14) ;
  230.    'ERRE' 5 ;
  231. 'FINSI' ;
  232. *
  233. 'FIN' ;
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  

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