Télécharger eautemp.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : eautemp.dgibi
  2. SAUT PAGE ;
  3. ********************************************************************
  4. * Test de validation des operateurs CHI1 et CHI2
  5. * eau + gradient thermique + pression CO2
  6. *********************************************************************
  7. A= 0. 0. ;
  8. B= 5. 0. ;
  9. OPTION ELEM QUA4 ;
  10. AB= A DROIT 1 B ;
  11. VEC1= 0. 40. ;
  12. VEC2= 0. 60. ;
  13. M1 = AB TRANS DINI 10. DFIN 20. VEC1 ;
  14. CD= INVE ( COTE 3 M1 );
  15. M2= M1 TRANS DINI 30. DFIN 50. VEC2 ;
  16. EF= COTE 3 M2 ;
  17. TASS M2 ;
  18. MT = CHANGE M2 POI1 ;
  19. XX YY = COOR MT ;
  20. *
  21. * table de données pour CHI1
  22. *
  23. TABDON=TABLE ;
  24. TABDON.BDD= 'STRASBG' ;
  25. TABDON.IDEN= LECT 1 50 60 61 101 165 ;
  26. TABDON.CHXMX= LECT 2144 2148 2157 2166 2168 2192 2198 2200
  27. 2224 2231 2249 2278 2281 2372 ;
  28. TABCLIM=TABLE ;
  29. TABCLIM.TYP6= LECT 2002 ;
  30. TABCLIM.TYP3= LECT 2400 ;
  31. TABCLIM.COMP3= LECT 101 ;
  32. TABDON.CLIM=TABCLIM ;
  33. TABDON.TEMPERATURE = 'OUI' ;
  34. *
  35. TB1=CHI1 TABDON
  36. LOGK '/u/castem/divers/COMPSM'
  37. ENTH '/u/castem/divers/COMPSM'
  38. COMP '/u/castem/divers/COMPSM' ;
  39. *
  40. *
  41. * table de données pour CHI2
  42. *
  43. TBPARM= TABLE ;
  44. TBPARM.ITMAX = 95;
  45. TBPARM.ITERSOLI = 10 ;
  46. TBPARM.NFI = 8;
  47. TBPARM.EPS= 1.D-4 ;
  48. TBPARM.SORTIE= MOTS 'PREC' 'FION' 'TYP5' ;
  49. TBPAR2= TABLE ;
  50. TBPAR2.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  51. TBPAR2.LOGC= MANU CHPO MT 6 X001 -5.6 X050 -7.0 X060 -5.0
  52. X061 -3. X101 -7. X165 -5. ;
  53. TOT001= MANU CHPO MT 1 X001 6.D-4 ;
  54. TOT050= MANU CHPO MT 1 X050 -1.15D-3 ;
  55. TOT060= MANU CHPO MT 1 X060 0.D0 ;
  56. TOT061= MANU CHPO MT 1 X061 2.8D-2 ;
  57. TOT101= MANU CHPO MT 1 X101 0.D0 ;
  58. TOT165= MANU CHPO MT 1 X165 5.5D-5 ;
  59. TBPAR2.TOT= TOT001+ TOT050 + TOT060 + TOT061 + TOT101 + TOT165 ;
  60. *
  61. * définition du champ de temperature
  62. * gradient constant entre l'entrée (180°C) et la sortie (25°C)
  63. *
  64. B= 180. ;
  65. A= -1.55 ;
  66. XX YY = COOR MT ;
  67. TMPE= A*YY + B ;
  68. TBPAR2.TEMPE= TMPE ;
  69. *
  70. *
  71. * calcul de la contrainte KCO2 pour avoir une pression partielle
  72. * de CO2 = 10**(-2.32)
  73. * KCO2 = log(k(T)) - log ( pco2)
  74. *
  75. COO = 2.32 ;
  76. *
  77. TSTA1 = 0.d0 ;
  78. TSTA2 = 25.d0 ;
  79. TSTA3 = 60.d0 ;
  80. TSTA4 = 100.d0;
  81. TSTA5 = 150.d0;
  82. TSTA6 = 200.d0;
  83. TSTA7 = 250.d0;
  84. TSTA8 = 300.d0;
  85. *
  86. K11 = -18.14 ;
  87. K21 = -18.13 ;
  88. K31 = -18.25 ;
  89. K41 = -18.58 ;
  90. K51 = -19.09 ;
  91. K61 = -19.82 ;
  92. K71 = -21.04 ;
  93. K81 = -24.07 ;
  94. K12 = 1.001232E-03 ;
  95. K22 = -8.024372E-04 ;
  96. K32 = -6.155775E-03 ;
  97. K42 = -9.206214E-03 ;
  98. K52 = -1.241482E-02 ;
  99. K62 = -1.553448E-02 ;
  100. K72 = -4.244730E-02 ;
  101. K82 = 0.000000E+00 ;
  102. K13 = 0.000000E+00 ;
  103. K23 = -7.214677E-05 ;
  104. K33 = -8.080574E-05 ;
  105. K43 = 4.544751E-06 ;
  106. K53 = -6.871693E-05 ;
  107. K63 = 6.323791E-06 ;
  108. K73 = -5.445802E-04 ;
  109. K83 = 0.000000E+00 ;
  110. K14 = -9.619569E-07 ;
  111. K24 = -8.246646E-08 ;
  112. K34 = 7.112541E-07 ;
  113. K44 = -4.884112E-07 ;
  114. K54 = 5.002715E-07 ;
  115. K64 = -3.672694E-06 ;
  116. K74 = 3.630535E-06 ;
  117. K84 = 0.000000E+00 ;
  118. *
  119. MASQ1= TMPE MASQ 'EGALE' TSTA1 ;
  120. KCO21= (-1.* K11) * MASQ1 ;
  121. MASQ1= TMPE MASQ 'EGALE' TSTA2 ;
  122. KCO22= (-1.* K21) * MASQ1 ;
  123. MASQ1= TMPE MASQ 'EGALE' TSTA3 ;
  124. KCO23= (-1.* K31) * MASQ1 ;
  125. MASQ1= TMPE MASQ 'EGALE' TSTA4 ;
  126. KCO24= (-1.* K41) * MASQ1 ;
  127. MASQ1= TMPE MASQ 'EGALE' TSTA5 ;
  128. KCO25= (-1.* K51) * MASQ1 ;
  129. MASQ1= TMPE MASQ 'EGALE' TSTA6 ;
  130. KCO26= (-1.* K61) * MASQ1 ;
  131. MASQ1= TMPE MASQ 'EGALE' TSTA7 ;
  132. KCO27= (-1.* K71) * MASQ1 ;
  133. MASQ1= TMPE MASQ 'EGALE' TSTA8 ;
  134. KCO28= (-1.* K81) * MASQ1 ;
  135. kco2= kco21 + kco22 +kco23 +kco24+kco25+kco26+kco27+kco28 ;
  136. TMPE2= TMPE - TSTA7 ;
  137. TSTAD = TSTA8 - TSTA7 ;
  138. KCOCC= MANU CHPO MT 1 SCAL (COO - K71) ;
  139. MASQ1= TMPE2 MASQ 'INFERIEUR' TSTAD ;
  140. TMPE1= MASQ1* TMPE2 ;
  141. KCOC1= MASQ1 * KCOCC ;
  142. MASQ1= TMPE1 MASQ 'SUPERIEUR' 0. ;
  143. TMPE2= MASQ1* TMPE1 ;
  144. KCOCC= MASQ1* KCOC1 ;
  145. KCOC1= ((((K74*TMPE2)+K73)*TMPE2)+K72)*TMPE2 ;
  146. KCOI7= KCOCC - KCOC1 ;
  147. TMPE2= TMPE - TSTA6 ;
  148. TSTAD = TSTA7 - TSTA6 ;
  149. KCOCC= MANU CHPO MT 1 SCAL (COO - K61) ;
  150. MASQ1= TMPE2 MASQ 'INFERIEUR' TSTAD ;
  151. TMPE1= MASQ1* TMPE2 ;
  152. KCOC1= MASQ1 * KCOCC ;
  153. MASQ1= TMPE1 MASQ 'SUPERIEUR' 0. ;
  154. TMPE2= MASQ1* TMPE1 ;
  155. KCOCC= MASQ1* KCOC1 ;
  156. KCOC1= ((((K64*TMPE2)+K63)*TMPE2)+K62)*TMPE2 ;
  157. KCOI6= KCOCC - KCOC1 ;
  158. TMPE2= TMPE - TSTA5 ;
  159. TSTAD = TSTA6 - TSTA5 ;
  160. KCOCC= MANU CHPO MT 1 SCAL (COO - K51) ;
  161. MASQ1= TMPE2 MASQ 'INFERIEUR' TSTAD ;
  162. TMPE1= MASQ1* TMPE2 ;
  163. KCOC1= MASQ1 * KCOCC ;
  164. MASQ1= TMPE1 MASQ 'SUPERIEUR' 0. ;
  165. TMPE2= MASQ1* TMPE1 ;
  166. KCOCC= MASQ1* KCOC1 ;
  167. KCOC1= ((((K54*TMPE2)+K53)*TMPE2)+K52)*TMPE2 ;
  168. KCOI5= KCOCC - KCOC1 ;
  169. TMPE2= TMPE - TSTA4 ;
  170. TSTAD = TSTA5 - TSTA4 ;
  171. KCOCC= MANU CHPO MT 1 SCAL (COO - K41) ;
  172. MASQ1= TMPE2 MASQ 'INFERIEUR' TSTAD ;
  173. TMPE1= MASQ1* TMPE2 ;
  174. KCOC1= MASQ1 * KCOCC ;
  175. MASQ1= TMPE1 MASQ 'SUPERIEUR' 0. ;
  176. TMPE2= MASQ1* TMPE1 ;
  177. KCOCC= MASQ1* KCOC1 ;
  178. KCOC1= ((((K44*TMPE2)+K43)*TMPE2)+K42)*TMPE2 ;
  179. KCOI4= KCOCC - KCOC1 ;
  180. TMPE2= TMPE - TSTA3 ;
  181. TSTAD = TSTA4 - TSTA3 ;
  182. KCOCC= MANU CHPO MT 1 SCAL (COO - K31) ;
  183. MASQ1= TMPE2 MASQ 'INFERIEUR' TSTAD ;
  184. TMPE1= MASQ1* TMPE2 ;
  185. KCOC1= MASQ1 * KCOCC ;
  186. MASQ1= TMPE1 MASQ 'SUPERIEUR' 0. ;
  187. TMPE2= MASQ1* TMPE1 ;
  188. KCOCC= MASQ1* KCOC1 ;
  189. KCOC1= ((((K34*TMPE2)+K33)*TMPE2)+K32)*TMPE2 ;
  190. KCOI3= KCOCC - KCOC1 ;
  191. TMPE2= TMPE - TSTA2 ;
  192. TSTAD = TSTA3 - TSTA2 ;
  193. KCOCC= MANU CHPO MT 1 SCAL (COO - K21) ;
  194. MASQ1= TMPE2 MASQ 'INFERIEUR' TSTAD ;
  195. TMPE1= MASQ1* TMPE2 ;
  196. KCOC1= MASQ1 * KCOCC ;
  197. MASQ1= TMPE1 MASQ 'SUPERIEUR' 0. ;
  198. TMPE2= MASQ1* TMPE1 ;
  199. KCOCC= MASQ1* KCOC1 ;
  200. KCOC1= ((((K24*TMPE2)+K23)*TMPE2)+K22)*TMPE2 ;
  201. KCOI2= KCOCC - KCOC1 ;
  202. TMPE2= TMPE - TSTA1 ;
  203. TSTAD = TSTA2 - TSTA1 ;
  204. KCOCC= MANU CHPO MT 1 SCAL (COO - K11) ;
  205. MASQ1= TMPE2 MASQ 'INFERIEUR' TSTAD ;
  206. TMPE1= MASQ1* TMPE2 ;
  207. KCOC1= MASQ1 * KCOCC ;
  208. MASQ1= TMPE1 MASQ 'SUPERIEUR' 0. ;
  209. TMPE2= MASQ1* TMPE1 ;
  210. KCOCC= MASQ1* KCOC1 ;
  211. KCOC1= ((((K14*TMPE2)+K13)*TMPE2)+K12)*TMPE2 ;
  212. KCOI1= KCOCC - KCOC1 ;
  213. KCOI= KCO2 + KCOI1 +KCOI2 +KCOI3+ KCOI4 +KCOI5+KCOI6+KCOI7 ;
  214. KCOI= KCOI- 18.13 ;
  215. KCO2= NOMC W016 KCOI ;
  216. TBPAR2.CLIM= KCO2 ;
  217. TB3= CHI2 TB1 TBPARM TBPAR2 ;
  218. *
  219. *
  220. * controle des résultats
  221. *
  222. PRECAB=REDU AB TB3.PREC ;
  223. PRECCD=REDU CD TB3.PREC ;
  224. PRECEF=REDU EF TB3.PREC ;
  225. PREC1= MANU CHPO AB 14 W018 0. W019 2.83058E-04 W020 2.49697E-02
  226. W021 0. W022 0. W023 0. W024 0. W025 0. W026 0. W027 0. W028 0.
  227. W029 0. W030 0. W031 0. 'NATURE' DISCRET ;
  228. PREC2= MANU CHPO CD 14 W018 0. W019 2.84979E-04 W020 2.67385E-02
  229. W021 0. W022 0. W023 0. W024 0. W025 1.84781E-05 W026 0. W027 0.
  230. W028 0. W029 0. W030 0. W031 0. 'NATURE' DISCRET ;
  231. PREC3= MANU CHPO EF 14 W018 0. W019 0. W020 2.78445E-02
  232. W021 2.74960E-05 W022 0. W023 0. W024 0. W025 0. W026 0. W027 0.
  233. W028 0. W029 0. W030 0. W031 0. 'NATURE' DISCRET ;
  234. *PRECMT= PRECAB+PRECCD+PRECEF ;
  235. PRECMT= PRECAB+PRECCD ;
  236. *PRECTES=PREC1+PREC2+PREC3 ;
  237. PRECTES=PREC1+PREC2 ;
  238. PRECD= PRECTES/50. ;
  239. VERR1= (ABS ( PRECTES - PRECMT)) MASQUE SUPERIEUR PRECD SOMME ;
  240. FIONAB= REDU AB TB3.FION ;
  241. FIONCD= REDU CD TB3.FION ;
  242. FIONEF= REDU EF TB3.FION ;
  243. *FIONMT= FIONAB+ FIONCD+FIONEF ;
  244. FIONMT= FIONAB+ FIONCD ;
  245. FIONTE1= MANU CHPO AB 1 SCAL 5.57739E-04 'NATURE' DISCRET ;
  246. FIONTE2= MANU CHPO CD 1 SCAL 8.16893E-04 'NATURE' DISCRET ;
  247. FIONTE3= MANU CHPO EF 1 SCAL 1.78045E-03 'NATURE' DISCRET ;
  248. *FIONTES= FIONTE1+FIONTE2+FIONTE3;
  249. FIONTES= FIONTE1+FIONTE2;
  250. FIOND= FIONTES/2.E5 ;
  251. VERR2= (ABS ( FIONTES - FIONMT)) MASQUE SUPERIEUR FIOND SOMME ;
  252. TYP5AB= REDU AB TB3.TYP5 ;
  253. TYP5CD= REDU CD TB3.TYP5 ;
  254. TYP5EF= REDU EF TB3.TYP5 ;
  255. TYP51=MANU CHPO AB 14 W018 2.04302E-02 W019 0. W020 0.
  256. W021 4.07418E-03 W022 9.08803E-04 W023 1.41667E-01 W024 6.49819E-04
  257. W025 2.92431E-02 W026 4.43300E-02 W027 2.49896E-01
  258. W028 2.15059E-02 W029 4.66326E-07 W030 2.92533E-07 W031 3.37606E-04
  259. 'NATURE' DISCRET ;
  260. TYP52=MANU CHPO CD 14 W018 1.42911E-01 W019 0. W020 0.
  261. W021 3.55161E-01 W022 2.35987E-02 W023 5.19600E-01 W024 4.09491E-02
  262. W025 0. W026 5.85450E-01 W027 1.97493E-01 W028 1.39410E-01
  263. W029 9.40715E-04 W030 2.15057E-04 W031 1.57713E-03 'NATURE' DISCRET ;
  264. TYP53=MANU CHPO EF 14 W018 2.97852E-01 W019 9.53362E-02 W020 0.
  265. W021 0. W022 2.24386E-02 W023 2.09538E-04 W024 6.91829E-02
  266. W025 1.06065E-01 W026 7.06324E-02 W027 1.04954E-01 W028 1.33045E-01
  267. W029 1.80148E-01 W030 2.17771E-04 W031 2.50476E-05 'NATURE' DISCRET ;
  268. *TYP5TES= TYP51+TYP52+TYP53 ;
  269. TYP5TES= TYP51+TYP52 ;
  270. *TYP5MT=TYP5AB+TYP5CD+TYP5EF ;
  271. TYP5MT=TYP5AB+TYP5CD ;
  272. TYP5D= TYP5TES/50. ;
  273. VERR3= (ABS ( TYP5TES - TYP5MT)) MASQUE SUPERIEUR TYP5D SOMME ;
  274. VERR= VERR1+VERR2+VERR3 ;
  275. SI (VERR EGA 0 ) ;
  276. ERREUR 0 ;
  277. SINON ;
  278. ERREUR 5 ;
  279. FINSI ;
  280. *
  281. fin ;
  282.  
  283.  
  284.  
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  

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