Télécharger isp472d_cond_Fick.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : isp472d_cond_Fick.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *****************************************************************
  5. * Cas test de non-régression pour la bas de jdd de CAST3M
  6. * Test de non régression du modèle de Fick
  7. *****************************************************************
  8. ** Jeu de données MISTRA pour le maillage 2D de l enceinte
  9. **
  10. ** Il a été choisi de mailler tous les volumes à charge
  11. ** ensuite à l utilisateur de ne retenir que ce qui
  12. ** l intéresse pour le calcul
  13. **
  14. ** Auteur E. Studer
  15. ** Date : 13/03/02
  16. *****************************************************************
  17.  
  18. 'OPTI' 'DIME' 2 ;
  19. 'OPTI' 'TRAC' 'PSC' ;
  20. 'OPTI' 'MODE' 'AXIS' ;
  21. 'OPTI' 'ELEM' 'QUA8' ;
  22.  
  23. RAF = 1.0 ;
  24. GRAPH = FAUX ;
  25. CTRV = FAUX ;
  26.  
  27. DMIN = 0.1 '/' RAF ;
  28. DMIN1 = 0.1 '/' RAF ;
  29. DMAX = 0.2 '/' RAF;
  30.  
  31. * Dimensions sur plan
  32. R0 = 2.125 ;
  33. R0S2 = R0 '/' 2.0 ;
  34.  
  35. EPB = 0.025 ;
  36. EPF = 0.015 ;
  37. EPH = 0.119 ;
  38. EPSI = 1.E-5 ;
  39.  
  40. HCOND1 = 2.187 ;
  41. HCOND2 = 1.784 ;
  42. HCOND3 = 1.784 ;
  43.  
  44. REQC = 1.9079 ;
  45. EPC = 0.040 ;
  46. EPCC = 0.120 ;
  47. HCC = 0.150 ;
  48. HEC = 0.120 ;
  49.  
  50.  
  51. EPBR = 0.110 ;
  52. HBR1 = 0.190 ;
  53. HBR2 = 0.195 ;
  54.  
  55. RINJ = 0.1 ;
  56. RS = 3.2529 ;
  57.  
  58. ZCOND1 = 1.285 ;
  59. ZCOND2 = 3.592 ;
  60. ZCOND3 = 5.4965 ;
  61.  
  62. ZHB1 = 1.938 ;
  63. ZHB2 = 4.341 ;
  64. ZHB3 = 6.559 ;
  65.  
  66. * Entiers pour le nombre de mailles
  67. NINJ = 2 '*' ('ENTIER' RAF) '-' 1 ;
  68. NX1 = 6 '*' ('ENTIER' RAF);
  69. NX2 = 2 '*' ('ENTIER' RAF) '-' 1 ;
  70. NX3 = 2 '*' ('ENTIER' RAF) '-' 1 ;
  71. NX4 = 2 '*' ('ENTIER' RAF) '-' 1 ;
  72.  
  73. NF5 = 2 '*' ('ENTIER' RAF) '-' 1 ;
  74. NM1 = 2 '*' ('ENTIER' RAF) '-' 1 ;
  75. NM2 = 2 '*' ('ENTIER' RAF) '-' 1 ;
  76. NZ1 = 2 '*' ('ENTIER' RAF) '-' 1 ;
  77. NZ2 = 6 '*' ('ENTIER' RAF);
  78. NZ3 = 2 '*' ('ENTIER' RAF) '-' 1 ;
  79. NZ4 = 6 '*' ('ENTIER' RAF);
  80. NZ5 = 3 '*' ('ENTIER' RAF);
  81. NZ5B = 3 '*' ('ENTIER' RAF) '-' 1 ;
  82.  
  83. * Injection
  84. A = 0.0001 ZCOND1 ;
  85. B = RINJ ZCOND1 ;
  86. AB = A 'DROI' NINJ B ;
  87. AB = AB 'COUL' 'BLAN' ;
  88.  
  89. * jusqu"a R/2
  90. C = R0S2 ZCOND1 ;
  91. BC = B 'DROI' C 'DINI' DMIN 'DFIN' DMAX ;
  92. BC = BC 'COUL' 'BLAN' ;
  93.  
  94. * jusqu"au condenseurs
  95. D = REQC ZCOND1 ;
  96. CD = C 'DROI' D 'DINI' DMAX 'DFIN' DMIN1 ;
  97. CD = CD 'COUL' 'BLAN' ;
  98.  
  99. * condenseurs plus isolant
  100. E = (REQC '+' EPC) ZCOND1 ;
  101. DE = D 'DROI' NX2 E ;
  102. DE = DE 'COUL' 'BLEU' ;
  103.  
  104. * collecteurs des condenseurs
  105. F = (REQC '+' EPC '+' EPCC) ZCOND1 ;
  106. EF = E 'DROI' NX3 F ;
  107. EF = EF 'COUL' 'JAUN' ;
  108.  
  109. * volume mort derriere les condenseurs
  110. G = R0 ZCOND1 ;
  111. FG = F 'DROI' NX4 G ;
  112. FG = FG 'COUL' 'BLAN' ;
  113.  
  114. * fut de l"enceinte
  115. H = (R0 '+' EPF) ZCOND1 ;
  116. GH = G 'DROI' NM1 H ;
  117. GH = GH 'COUL' 'ROUG' ;
  118.  
  119. * Brides
  120. I = (R0 '+' EPF '+' EPBR) ZCOND1 ;
  121. HI = H 'DROI' NM2 I ;
  122. HI = HI 'COUL' 'VERT' ;
  123.  
  124. * Passage au maillage 2D
  125. BASE = BC 'ET' CD 'ET' (DE 'COUL' 'BLAN') 'ET'
  126. (EF 'COUL' 'BLAN') 'ET' FG 'ET' GH 'ET' HI ;
  127. BASE_1 = BASE 'PLUS' (0. (HBR1 '*' -1.0)) ;
  128. VB_1 = 'TRAN' NZ1 BASE (0. (HBR1 '*' -1.0)) ;
  129. 'ELIM' EPSI (BASE_1 'ET' VB_1) ;
  130.  
  131. 'SI' GRAPH ;
  132. 'TRACE' 'CACH' (VB_1)
  133. 'TITR' 'Volume VB_1 ' ;
  134. 'FINSI' ;
  135.  
  136. V1 = 0. ZCOND1 ;
  137. BASE_21 = (BASE_1 'ELEM' 'BLAN') ;
  138. BASE_22 = (BASE_1 'ELEM' 'ROUG') ;
  139. BASE_31 = BASE_21 'PROJ' 'CYLI' (-1.0 '*' V1) 'CERC' (0. RS)
  140. ( 0. 0. ) ;
  141. BASE_32 = BASE_22 'PROJ' 'CYLI' (-1.0 '*' V1) 'CERC' (0. RS)
  142. ( 0. 0. ) ;
  143. BASE_3 = BASE_31 'ET' BASE_32 ;
  144. 'ELIM' EPSI BASE_3 ;
  145.  
  146. VB_2_1 = 'REGLER' BASE_21 BASE_31 NZ2 ;
  147. VB_2_2 = 'REGLER' BASE_22 BASE_32 NZ2 ;
  148. VB_2 = VB_2_1 'ET' VB_2_2 ;
  149. 'ELIM' EPSI VB_2 ;
  150.  
  151. 'SI' GRAPH ;
  152. 'TRACE' 'CACH' (VB_1 'ET' VB_2)
  153. 'TITR' 'Volume VB_1,VB_2 ' ;
  154. 'FINSI' ;
  155.  
  156. BASE_4 = BASE_3 'HOMO' (0. RS) (RS '+' EPB '/' RS) ;
  157. VB_3 = 'REGLER' BASE_3 BASE_4 NM1 ;
  158. VB_3 = VB_3 'COUL' 'ROUG' ;
  159. 'SI' GRAPH ;
  160. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3)
  161. 'TITR' 'Volume VB_1,VB_2,VB_3 ' ;
  162. 'FINSI' ;
  163.  
  164. * Maintenant on monte en Z
  165. *
  166. * Collecteur
  167. HAUT = AB 'ET' BC 'ET' CD 'ET' DE 'ET' EF 'ET' FG 'ET'
  168. GH 'ET' HI ;
  169. HAUT_2 = HAUT 'PLUS' (0. HCC) ;
  170. VB_4 = 'TRAN' NZ3 HAUT (0. HCC) ;
  171. 'ELIM' EPSI (HAUT_2 'ET' VB_4) ;
  172. H_V = VB_4 'ELEM' 'VERT' ;
  173. VB_4 = 'DIFF' VB_4 H_V ;
  174. SC1_1 = INTER ('CONT' (VB_4 'ELEM' 'BLAN'))
  175. ('CONT' (VB_4 'ELEM' 'BLEU'));
  176. * il faut construire l"axe
  177. CNT = 'CONT' VB_4 ;
  178. CNT = 'DIFF' CNT HAUT ;
  179. CNT = 'DIFF' CNT HAUT_2 ;
  180. CNT = 'CCON' CNT ;
  181. AXE = CNT.2 ;
  182. 'ELIM' EPSI (AXE 'ET' VB_4) ;
  183.  
  184. 'SI' GRAPH ;
  185. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4)
  186. 'TITR' 'Volume VB_1,VB_2,VB_3,VB_4 ' ;
  187. 'FINSI' ;
  188.  
  189. *Bride et le collecteur devient fluide (BLAN)
  190. HAUT_3 = HAUT_2 'PLUS' (0. HBR2) ;
  191. VB_5 = 'TRAN' NZ1 HAUT_2 (0. HBR2) ;
  192. 'ELIM' EPSI (HAUT_3 'ET' VB_5) ;
  193. SC1_2 = INTER ('CONT' (VB_5 'ELEM' 'BLAN'))
  194. ('CONT' (VB_5 'ELEM' 'BLEU'));
  195. H1_J = VB_5 'ELEM' 'JAUN' ;
  196. H1_J = H1_J 'COUL' 'BLAN' ;
  197. VB_5 = 'DIFF' VB_5 H1_J ;
  198. VB_5 = VB_5 'ET' H1_J ;
  199.  
  200. CNT = 'CONT' VB_5 ;
  201. CNT = 'DIFF' CNT HAUT_2 ;
  202. CNT = 'DIFF' CNT HAUT_3 ;
  203. CNT = 'CCON' CNT ;
  204. AXE1 = CNT.2 ;
  205. 'ELIM' EPSI (AXE1 'ET' VB_5) ;
  206. AXE = AXE 'ET' AXE1 ;
  207. 'SI' GRAPH ;
  208. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  209. 'ET' VB_5 )
  210. 'TITR' 'Volume VB_1,VB_2,VB_3,VB_4,VB_5 ' ;
  211. 'FINSI' ;
  212.  
  213. * entre les deux brides on ajoute un plan pour le hublot
  214. Z3 = ( 'MAXI' ('COOR' HAUT_3 2)) ;
  215. HAUT_3M = HAUT_3 'PLUS' (0. (ZHB1 '-' Z3)) ;
  216. HAUT_4 = HAUT_3 'PLUS'
  217. (0. (HCOND1 '-' HBR2 '-' HBR2 '-' HCC '-' HCC )) ;
  218. VB_6_1 = 'TRAN' NZ5B HAUT_3 (0. (ZHB1 '-' Z3)) ;
  219. VB_6_2 = 'TRAN' NZ4 HAUT_3M
  220. (0. (HCOND1 '-' HBR2 '-' HBR2 '-' HCC '-' HCC '-' ZHB1 '+' Z3)) ;
  221. VB_6 = VB_6_1 'ET' VB_6_2 ;
  222. 'ELIM' EPSI (HAUT_3M 'ET' HAUT_4 'ET' VB_6) ;
  223.  
  224. SC1_3 = INTER ('CONT' (VB_6 'ELEM' 'BLAN'))
  225. ('CONT' (VB_6 'ELEM' 'BLEU'));
  226. H2_V = VB_6 'ELEM' 'VERT' ;
  227. VB_6 = 'DIFF' VB_6 H2_V ;
  228. H1_J = VB_6 'ELEM' 'JAUN' ;
  229. H1_J = H1_J 'COUL' 'BLAN' ;
  230. VB_6 = 'DIFF' VB_6 H1_J ;
  231. VB_6 = VB_6 'ET' H1_J ;
  232.  
  233. CNT = 'CONT' VB_6 ;
  234. CNT = 'DIFF' CNT HAUT_3 ;
  235. CNT = 'DIFF' CNT HAUT_4 ;
  236. CNT = 'CCON' CNT ;
  237. AXE1 = CNT.2 ;
  238. 'ELIM' EPSI (AXE1 'ET' VB_6) ;
  239. AXE = AXE 'ET' AXE1 ;
  240. 'SI' GRAPH ;
  241. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  242. 'ET' VB_5 'ET' VB_6 )
  243. 'TITR' 'Volume VB_1,VB_2,VB_3,VB_4,VB_5,VB_6 ' ;
  244. 'FINSI' ;
  245.  
  246. * une nouvelle bride
  247. HAUT_5 = HAUT_4 'PLUS' (0. HBR2) ;
  248. VB_7 = 'TRAN' NZ1 HAUT_4 (0. HBR2) ;
  249. 'ELIM' EPSI (HAUT_5 'ET' VB_7) ;
  250. SC1_4 = INTER ('CONT' (VB_7 'ELEM' 'BLAN'))
  251. ('CONT' (VB_7 'ELEM' 'BLEU'));
  252. H1_J = VB_7 'ELEM' 'JAUN' ;
  253. H1_J = H1_J 'COUL' 'BLAN' ;
  254. VB_7 = 'DIFF' VB_7 H1_J ;
  255. VB_7 = VB_7 'ET' H1_J ;
  256. CNT = 'CONT' VB_7 ;
  257. CNT = 'DIFF' CNT HAUT_4 ;
  258. CNT = 'DIFF' CNT HAUT_5 ;
  259. CNT = 'CCON' CNT ;
  260. AXE1 = CNT.2 ;
  261. 'ELIM' EPSI (AXE1 'ET' VB_7) ;
  262. AXE = AXE 'ET' AXE1 ;
  263.  
  264. 'SI' GRAPH ;
  265. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  266. 'ET' VB_5 'ET' VB_6 'ET' VB_7 )
  267. 'TITR' 'VB_1,VB_2,VB_3,VB_4,VB_5,VB_6, VB_7 ' ;
  268. 'FINSI' ;
  269.  
  270. * un collecteur haut du condenseur bas
  271. HAUT_6 = HAUT_5 'PLUS' (0. HCC) ;
  272. VB_8 = 'TRAN' NZ1 HAUT_5 (0. HCC) ;
  273. 'ELIM' EPSI (HAUT_6 'ET' VB_8) ;
  274. SC1_5 = INTER ('CONT' (VB_8 'ELEM' 'BLAN'))
  275. ('CONT' (VB_8 'ELEM' 'BLEU'));
  276. H2_V = VB_8 'ELEM' 'VERT' ;
  277. VB_8 = 'DIFF' VB_8 H2_V ;
  278. CNT = 'CONT' VB_8 ;
  279. CNT = 'DIFF' CNT HAUT_5 ;
  280. CNT = 'DIFF' CNT HAUT_6 ;
  281. CNT = 'CCON' CNT ;
  282. AXE1 = CNT.2 ;
  283. 'ELIM' EPSI (AXE1 'ET' VB_8) ;
  284. AXE = AXE 'ET' AXE1 ;
  285.  
  286. 'SI' GRAPH ;
  287. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  288. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 )
  289. 'TITR' 'VB1,VB2,VB3,VB4,VB5,VB6,VB7,VB8 ' ;
  290. 'FINSI' ;
  291.  
  292. * entre le condenseur 1 et le condenseur 2
  293. HAUT_7 = HAUT_6 'PLUS' (0. HEC) ;
  294. VB_9 = 'TRAN' NZ1 HAUT_6 (0. HEC) ;
  295. 'ELIM' EPSI (HAUT_7 'ET' VB_9) ;
  296. H2_V = VB_9 'ELEM' 'VERT' ;
  297. VB_9 = 'DIFF' VB_9 H2_V ;
  298. H1_J = VB_9 'ELEM' 'JAUN' ;
  299. H1_J = H1_J 'COUL' 'BLAN' ;
  300. VB_9 = 'DIFF' VB_9 H1_J ;
  301. VB_9 = VB_9 'ET' H1_J ;
  302. H1_B = VB_9 'ELEM' 'BLEU' ;
  303. H1_B = H1_B 'COUL' 'BLAN' ;
  304. VB_9 = 'DIFF' VB_9 H1_B ;
  305. VB_9 = VB_9 'ET' H1_B ;
  306. CNT = 'CONT' VB_9 ;
  307. CNT = 'DIFF' CNT HAUT_6 ;
  308. CNT = 'DIFF' CNT HAUT_7 ;
  309. CNT = 'CCON' CNT ;
  310. AXE1 = CNT.2 ;
  311. 'ELIM' EPSI (AXE1 'ET' VB_9) ;
  312. AXE = AXE 'ET' AXE1 ;
  313. 'SI' GRAPH ;
  314. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  315. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9 )
  316. 'TITR' 'VB1,VB2,VB3,VB4,VB5,VB6,VB7,VB8,VB9 ' ;
  317. 'FINSI' ;
  318.  
  319. * collecteur bas du condenseur moyen
  320. HAUT_8 = HAUT_7 'PLUS' (0. HCC) ;
  321. VB_10 = 'TRAN' NZ1 HAUT_7 (0. HCC) ;
  322. 'ELIM' EPSI (HAUT_8 'ET' VB_10) ;
  323. SC2_1 = INTER ('CONT' (VB_10 'ELEM' 'BLAN'))
  324. ('CONT' (VB_10 'ELEM' 'BLEU'));
  325. H2_V = VB_10 'ELEM' 'VERT' ;
  326. VB_10 = 'DIFF' VB_10 H2_V ;
  327. CNT = 'CONT' VB_10 ;
  328. CNT = 'DIFF' CNT HAUT_7 ;
  329. CNT = 'DIFF' CNT HAUT_8 ;
  330. CNT = 'CCON' CNT ;
  331. AXE1 = CNT.2 ;
  332. 'ELIM' EPSI (AXE1 'ET' VB_10) ;
  333. AXE = AXE 'ET' AXE1 ;
  334. 'SI' GRAPH ;
  335. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  336. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  337. 'ET' VB_10 )
  338. 'TITR' 'VB1,VB2,VB3,VB4,VB5,VB6,VB7,VB8,VB9,VB10 ' ;
  339. 'FINSI' ;
  340.  
  341. * Bride du condenseur 2
  342. HAUT_9 = HAUT_8 'PLUS' (0. HBR2) ;
  343. VB_11 = 'TRAN' NZ1 HAUT_8 (0. HBR2) ;
  344. 'ELIM' EPSI (HAUT_9 'ET' VB_11) ;
  345. SC2_2 = INTER ('CONT' (VB_11 'ELEM' 'BLAN'))
  346. ('CONT' (VB_11 'ELEM' 'BLEU'));
  347. H1_J = VB_11 'ELEM' 'JAUN' ;
  348. H1_J = H1_J 'COUL' 'BLAN' ;
  349. VB_11 = 'DIFF' VB_11 H1_J ;
  350. VB_11 = VB_11 'ET' H1_J ;
  351. CNT = 'CONT' VB_11 ;
  352. CNT = 'DIFF' CNT HAUT_8 ;
  353. CNT = 'DIFF' CNT HAUT_9 ;
  354. CNT = 'CCON' CNT ;
  355. AXE1 = CNT.2 ;
  356. 'ELIM' EPSI (AXE1 'ET' VB_11) ;
  357. AXE = AXE 'ET' AXE1 ;
  358.  
  359. 'SI' GRAPH ;
  360. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  361. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  362. 'ET' VB_10 'ET' VB_11 )
  363. 'TITR' 'VB1 A VB11 ' ;
  364. 'FINSI' ;
  365.  
  366. * entre bride et collecteur haut pour condenseur moyen
  367. Z3 = ( 'MAXI' ('COOR' HAUT_9 2)) ;
  368. HAUT_9M = HAUT_9 'PLUS' (0. (ZHB2 '-' Z3)) ;
  369. HAUT_10 = HAUT_9 'PLUS'
  370. (0. (HCOND2 '-' HBR2 '-' HCC '-' HCC )) ;
  371. VB_12_1 = 'TRAN' NZ5 HAUT_9 (0. (ZHB2 '-' Z3)) ;
  372. VB_12_2 = 'TRAN' NZ4 HAUT_9M
  373. (0. (HCOND2 '-' HBR2 '-' HCC '-' HCC '-' ZHB2 '+' Z3)) ;
  374. VB_12 = VB_12_1 'ET' VB_12_2 ;
  375. 'ELIM' EPSI (HAUT_9M 'ET' HAUT_10 'ET' VB_12) ;
  376. SC2_3 = INTER ('CONT' (VB_12 'ELEM' 'BLAN'))
  377. ('CONT' (VB_12 'ELEM' 'BLEU'));
  378. H2_V = VB_12 'ELEM' 'VERT' ;
  379. VB_12 = 'DIFF' VB_12 H2_V ;
  380. H1_J = VB_12 'ELEM' 'JAUN' ;
  381. H1_J = H1_J 'COUL' 'BLAN' ;
  382. VB_12 = 'DIFF' VB_12 H1_J ;
  383. VB_12 = VB_12 'ET' H1_J ;
  384. CNT = 'CONT' VB_12 ;
  385. CNT = 'DIFF' CNT HAUT_9 ;
  386. CNT = 'DIFF' CNT HAUT_10 ;
  387. CNT = 'CCON' CNT ;
  388. AXE1 = CNT.2 ;
  389. 'ELIM' EPSI (AXE1 'ET' VB_12) ;
  390. AXE = AXE 'ET' AXE1 ;
  391. 'SI' GRAPH ;
  392. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  393. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  394. 'ET' VB_10 'ET' VB_11 'ET' VB_12 )
  395. 'TITR' 'VB1 A VB12 ' ;
  396. 'FINSI' ;
  397.  
  398. * collecteur haut du condenseur moyen
  399. HAUT_11 = HAUT_10 'PLUS' (0. HCC) ;
  400. VB_13 = 'TRAN' NZ1 HAUT_10 (0. HCC) ;
  401. 'ELIM' EPSI (HAUT_11 'ET' VB_13) ;
  402. SC2_4 = INTER ('CONT' (VB_13 'ELEM' 'BLAN'))
  403. ('CONT' (VB_13 'ELEM' 'BLEU'));
  404. H2_V = VB_13 'ELEM' 'VERT' ;
  405. VB_13 = 'DIFF' VB_13 H2_V ;
  406. CNT = 'CONT' VB_13 ;
  407. CNT = 'DIFF' CNT HAUT_10 ;
  408. CNT = 'DIFF' CNT HAUT_11 ;
  409. CNT = 'CCON' CNT ;
  410. AXE1 = CNT.2 ;
  411. 'ELIM' EPSI (AXE1 'ET' VB_13) ;
  412. AXE = AXE 'ET' AXE1 ;
  413. 'SI' GRAPH ;
  414. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  415. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  416. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 )
  417. 'TITR' 'VB1 A VB13 ' ;
  418. 'FINSI' ;
  419.  
  420. * entre le condenseur moyen et le condenseur haut
  421. HAUT_12 = HAUT_11 'PLUS' (0. HEC) ;
  422. VB_14 = 'TRAN' NZ1 HAUT_11 (0. HEC) ;
  423. 'ELIM' EPSI (HAUT_12 'ET' VB_14) ;
  424. H2_V = VB_14 'ELEM' 'VERT' ;
  425. VB_14 = 'DIFF' VB_14 H2_V ;
  426. H1_J = VB_14 'ELEM' 'JAUN' ;
  427. H1_J = H1_J 'COUL' 'BLAN' ;
  428. VB_14 = 'DIFF' VB_14 H1_J ;
  429. VB_14 = VB_14 'ET' H1_J ;
  430. H1_B = VB_14 'ELEM' 'BLEU' ;
  431. H1_B = H1_B 'COUL' 'BLAN' ;
  432. VB_14 = 'DIFF' VB_14 H1_B ;
  433. VB_14 = VB_14 'ET' H1_B ;
  434. CNT = 'CONT' VB_14 ;
  435. CNT = 'DIFF' CNT HAUT_11 ;
  436. CNT = 'DIFF' CNT HAUT_12 ;
  437. CNT = 'CCON' CNT ;
  438. AXE1 = CNT.2 ;
  439. 'ELIM' EPSI (AXE1 'ET' VB_14) ;
  440. AXE = AXE 'ET' AXE1 ;
  441. 'SI' GRAPH ;
  442. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  443. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  444. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  445. VB_14 )
  446. 'TITR' 'VB1 A VB14 ' ;
  447. 'FINSI' ;
  448.  
  449. * collecteur bas du condenseur haut
  450. HAUT_13 = HAUT_12 'PLUS' (0. HCC) ;
  451. VB_15 = 'TRAN' NZ1 HAUT_12 (0. HCC) ;
  452. 'ELIM' EPSI (HAUT_13 'ET' VB_15) ;
  453. SC3_1 = INTER ('CONT' (VB_15 'ELEM' 'BLAN'))
  454. ('CONT' (VB_15 'ELEM' 'BLEU'));
  455. H2_V = VB_15 'ELEM' 'VERT' ;
  456. VB_15 = 'DIFF' VB_15 H2_V ;
  457. CNT = 'CONT' VB_15 ;
  458. CNT = 'DIFF' CNT HAUT_12 ;
  459. CNT = 'DIFF' CNT HAUT_13 ;
  460. CNT = 'CCON' CNT ;
  461. AXE1 = CNT.2 ;
  462. 'ELIM' EPSI (AXE1 'ET' VB_15) ;
  463. AXE = AXE 'ET' AXE1 ;
  464. 'SI' GRAPH ;
  465. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  466. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  467. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  468. VB_14 'ET' VB_15 )
  469. 'TITR' 'VB1 A VB15 ' ;
  470. 'FINSI' ;
  471.  
  472. * entre collecteur bas et bas de la bride sup
  473. Z3 = ( 'MAXI' ('COOR' HAUT_13 2)) ;
  474. HAUT_13M = HAUT_13 'PLUS' (0. (ZHB3 '-' Z3)) ;
  475. HAUT_14 = HAUT_13 'PLUS'
  476. (0. (HCOND3 '-' HBR2 '-' HCC '-' HCC )) ;
  477. VB_16_1 = 'TRAN' NZ4 HAUT_13 (0. (ZHB3 '-' Z3)) ;
  478. VB_16_2 = 'TRAN' NZ5 HAUT_13M
  479. (0. (HCOND3 '-' HBR2 '-' HCC '-' HCC '-' ZHB3 '+' Z3)) ;
  480. VB_16 = VB_16_1 'ET' VB_16_2 ;
  481. 'ELIM' EPSI (HAUT_13M 'ET' HAUT_14 'ET' VB_16) ;
  482. SC3_2 = INTER ('CONT' (VB_16 'ELEM' 'BLAN'))
  483. ('CONT' (VB_16 'ELEM' 'BLEU'));
  484. H2_V = VB_16 'ELEM' 'VERT' ;
  485. VB_16 = 'DIFF' VB_16 H2_V ;
  486. H1_J = VB_16 'ELEM' 'JAUN' ;
  487. H1_J = H1_J 'COUL' 'BLAN' ;
  488. VB_16 = 'DIFF' VB_16 H1_J ;
  489. VB_16 = VB_16 'ET' H1_J ;
  490. CNT = 'CONT' VB_16 ;
  491. CNT = 'DIFF' CNT HAUT_13 ;
  492. CNT = 'DIFF' CNT HAUT_14 ;
  493. CNT = 'CCON' CNT ;
  494. AXE1 = CNT.2 ;
  495. 'ELIM' EPSI (AXE1 'ET' VB_16) ;
  496. AXE = AXE 'ET' AXE1 ;
  497. 'SI' GRAPH ;
  498. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  499. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  500. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  501. VB_14 'ET' VB_15 'ET' VB_16 )
  502. 'TITR' 'VB1 A VB16 ' ;
  503. 'FINSI' ;
  504.  
  505. * une nouvelle bride
  506. HAUT_15 = HAUT_14 'PLUS' (0. HBR2) ;
  507. VB_17 = 'TRAN' NZ1 HAUT_14 (0. HBR2) ;
  508. 'ELIM' EPSI (HAUT_15 'ET' VB_17) ;
  509. SC3_3 = INTER ('CONT' (VB_17 'ELEM' 'BLAN'))
  510. ('CONT' (VB_17 'ELEM' 'BLEU'));
  511. H1_J = VB_17 'ELEM' 'JAUN' ;
  512. H1_J = H1_J 'COUL' 'BLAN' ;
  513. VB_17 = 'DIFF' VB_17 H1_J ;
  514. VB_17 = VB_17 'ET' H1_J ;
  515. CNT = 'CONT' VB_17 ;
  516. CNT = 'DIFF' CNT HAUT_14 ;
  517. CNT = 'DIFF' CNT HAUT_15 ;
  518. CNT = 'CCON' CNT ;
  519. AXE1 = CNT.2 ;
  520. 'ELIM' EPSI (AXE1 'ET' VB_17) ;
  521. AXE = AXE 'ET' AXE1 ;
  522. 'SI' GRAPH ;
  523. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  524. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  525. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  526. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 )
  527. 'TITR' 'VB1 A VB17 ' ;
  528. 'FINSI' ;
  529.  
  530. * un collecteur haut du condenseur haut
  531. HAUT_16 = HAUT_15 'PLUS' (0. HCC) ;
  532. VB_18 = 'TRAN' NZ1 HAUT_15 (0. HCC) ;
  533. 'ELIM' EPSI (HAUT_16 'ET' VB_18) ;
  534. SC3_4 = INTER ('CONT' (VB_18 'ELEM' 'BLAN'))
  535. ('CONT' (VB_18 'ELEM' 'BLEU'));
  536. H2_V = VB_18 'ELEM' 'VERT' ;
  537. VB_18 = 'DIFF' VB_18 H2_V ;
  538. CNT = 'CONT' VB_18 ;
  539. CNT = 'DIFF' CNT HAUT_15 ;
  540. CNT = 'DIFF' CNT HAUT_16 ;
  541. CNT = 'CCON' CNT ;
  542. AXE1 = CNT.2 ;
  543. 'ELIM' EPSI (AXE1 'ET' VB_18) ;
  544. AXE = AXE 'ET' AXE1 ;
  545. 'SI' GRAPH ;
  546. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  547. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  548. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  549. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 'ET' VB_18 )
  550. 'TITR' 'VB1 A VB18 ' ;
  551. 'FINSI' ;
  552.  
  553. * entre le condenseur haut et le plafond
  554. HAUT_17 = HAUT_16 'PLUS' (0. 0.100) ;
  555. VB_19 = 'TRAN' NZ1 HAUT_16 (0. 0.100) ;
  556. 'ELIM' EPSI (HAUT_17 'ET' VB_19) ;
  557. H2_V = VB_19 'ELEM' 'VERT' ;
  558. VB_19 = 'DIFF' VB_19 H2_V ;
  559. H1_J = VB_19 'ELEM' 'JAUN' ;
  560. H1_J = H1_J 'COUL' 'BLAN' ;
  561. VB_19 = 'DIFF' VB_19 H1_J ;
  562. VB_19 = VB_19 'ET' H1_J ;
  563. H1_B = VB_19 'ELEM' 'BLEU' ;
  564. H1_B = H1_B 'COUL' 'BLAN' ;
  565. VB_19 = 'DIFF' VB_19 H1_B ;
  566. VB_19 = VB_19 'ET' H1_B ;
  567. CNT = 'CONT' VB_19 ;
  568. CNT = 'DIFF' CNT HAUT_16 ;
  569. CNT = 'DIFF' CNT HAUT_17 ;
  570. CNT = 'CCON' CNT ;
  571. AXE1 = CNT.2 ;
  572. 'ELIM' EPSI (AXE1 'ET' VB_19) ;
  573. AXE = AXE 'ET' AXE1 ;
  574. 'SI' GRAPH ;
  575. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  576. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  577. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  578. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 'ET' VB_18
  579. 'ET' VB_19 )
  580. 'TITR' 'VB1 A VB19' ;
  581. 'FINSI' ;
  582.  
  583. * le plafond
  584. HAUT_18 = HAUT_17 'PLUS' (0. EPH) ;
  585. VB_20 = 'TRAN' NM1 HAUT_17 (0. EPH) ;
  586. 'ELIM' EPSI (HAUT_18 'ET' VB_20) ;
  587. H2_V = VB_20 'ELEM' 'VERT' ;
  588. VB_20 = 'DIFF' VB_20 H2_V ;
  589. H1_J = VB_20 'ELEM' 'JAUN' ;
  590. H1_J = H1_J 'COUL' 'ROUG' ;
  591. VB_20 = 'DIFF' VB_20 H1_J ;
  592. VB_20 = VB_20 'ET' H1_J ;
  593. H1_B = VB_20 'ELEM' 'BLEU' ;
  594. H1_B = H1_B 'COUL' 'ROUG' ;
  595. VB_20 = 'DIFF' VB_20 H1_B ;
  596. VB_20 = VB_20 'ET' H1_B ;
  597. H1_BL = VB_20 'ELEM' 'BLAN' ;
  598. H1_BL = H1_BL 'COUL' 'ROUG' ;
  599. VB_20 = 'DIFF' VB_20 H1_BL ;
  600. VB_20 = VB_20 'ET' H1_BL ;
  601. 'SI' GRAPH ;
  602. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  603. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  604. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  605. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 'ET' VB_18
  606. 'ET' VB_19 'ET' VB_20 )
  607. 'TITR' 'VB1 A VB20' ;
  608. 'FINSI' ;
  609.  
  610. ***************************************************************
  611. * On effectue des controles par rapport aux plans usines
  612. ***************************************************************
  613. MISTRA = (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  614. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  615. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  616. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 'ET' VB_18
  617. 'ET' VB_19 'ET' VB_20 )
  618. ;
  619.  
  620. 'SI' CTRV ;
  621.  
  622. * volume du fluide
  623. VFF = MISTRA 'ELEM' 'BLAN' ;
  624. *'TRACE' 'CACH' VFF 'TITR' 'Volume fluide' ;
  625. MVFF = 'CHAN' VFF 'QUAF' ;
  626. $VFF = 'MODE' MVFF 'NAVIER_STOKES' 'QUAF' ;
  627. VVFF = 'DOMA' $VFF 'XXDIAGSI' ;
  628. VFF = 'DOMA' $VFF 'MAILLAGE' ;
  629. 'MESS' '==> Volume du fluide : '
  630. ('SOMT' VVFF) ;
  631.  
  632. * Elevation
  633. ZZ = 'COOR' VFF 2 ;
  634. 'MESS' '==> Zfluide min =' ('MINI' ZZ) ;
  635. 'MESS' '==> Zfluide max =' ('MAXI' ZZ) ;
  636.  
  637. * Surfaces de condensation
  638. SC1 = SC1_1 'ET' SC1_2 'ET' SC1_3 'ET' SC1_4 'ET'
  639. SC1_5 ;
  640. SC2 = SC2_1 'ET' SC2_2 'ET' SC2_3 'ET' SC2_4 ;
  641. SC3 = SC3_1 'ET' SC3_2 'ET' SC3_3 'ET' SC3_4 ;
  642. *'TRACE' 'CACH' (SC1 'ET' SC2 'ET' SC3)
  643. * 'TITR' 'Surface des condenseurs' ;
  644. MS1 = 'CHAN' SC1 'QUAF' ;
  645. $S1 = 'MODE' MS1 'NAVIER_STOKES' 'QUAF' ;
  646. VS1 = 'DOMA' $S1 'XXDIAGSI' ;
  647. MS2 = 'CHAN' SC2 'QUAF' ;
  648. $S2 = 'MODE' MS2 'NAVIER_STOKES' 'QUAF' ;
  649. VS2 = 'DOMA' $S2 'XXDIAGSI' ;
  650. MS3 = 'CHAN' SC3 'QUAF' ;
  651. $S3 = 'MODE' MS3 'NAVIER_STOKES' 'QUAF' ;
  652. VS3 = 'DOMA' $S3 'XXDIAGSI' ;
  653. 'MESS' '==> Surface condenseur BAS : '
  654. ('SOMT' VS1) ;
  655. 'MESS' '==> Surface condenseur MOY : '
  656. ('SOMT' VS2) ;
  657. 'MESS' '==> Surface condenseur HAU : '
  658. ('SOMT' VS3) ;
  659.  
  660. * Volume des condenseurs
  661. CS = (MISTRA 'ELEM' 'BLEU') 'ET' (MISTRA 'ELEM' 'JAUN') ;
  662. *'TRACE' 'CACH' CS 'TITR' 'Volume des condenseurs' ;
  663. MCS = 'CHAN' CS 'QUAF' ;
  664. $CS = 'MODE' MCS 'NAVIER_STOKES' 'QUAF' ;
  665. VCS = 'DOMA' $CS 'XXDIAGSI' ;
  666. 'MESS' '==> Volume des condenseurs : '
  667. ('SOMT' VCS);
  668.  
  669. * Volume des brides
  670. BRIDE = MISTRA 'ELEM' 'VERT' ;
  671. *'TRACE' 'CACH' BRIDE 'TITR' 'Volume des brides' ;
  672. MBR = 'CHAN' BRIDE 'QUAF' ;
  673. $BR = 'MODE' MBR 'NAVIER_STOKES' 'QUAF' ;
  674. VBR = 'DOMA' $BR 'XXDIAGSI' ;
  675. 'MESS' '==> Masse des Brides : '
  676. (('SOMT' VBR) '*' 7800.0) ;
  677.  
  678. * Volume acier enceinte
  679. ACI = MISTRA 'ELEM' 'ROUG' ;
  680. *'TRACE' 'CACH' ACI 'TITR' 'Volume Acier Enceinte' ;
  681. MAC = 'CHAN' ACI 'QUAF' ;
  682. $AC = 'MODE' MAC 'NAVIER_STOKES' 'QUAF' ;
  683. VAC = 'DOMA' $AC 'XXDIAGSI' ;
  684. 'MESS' '==> Masse des Aciers : '
  685. (('SOMT' VAC) '*' 7800.0) ;
  686. 'FINSI' ;
  687.  
  688. *************************************************************
  689. * On s"occupe maintenant du maillage 2D Complet
  690. *
  691. *************************************************************
  692. SC1 = SC1_1 'ET' SC1_2 'ET' SC1_3 'ET' SC1_4 'ET'
  693. SC1_5 ;
  694. SC2 = SC2_1 'ET' SC2_2 'ET' SC2_3 'ET' SC2_4 ;
  695. SC3 = SC3_1 'ET' SC3_2 'ET' SC3_3 'ET' SC3_4 ;
  696. *************************************************************
  697. ** FIN DU MAILLAGE **
  698. *************************************************************
  699. VFF = MISTRA 'ELEM' 'BLAN' ;
  700. PAR = (MISTRA 'ELEM' 'ROUG') 'ET' (MISTRA 'ELEM' 'VERT') ;
  701.  
  702. * Recherche de la paroi externe
  703.  
  704. CPAR = 'CONT' PAR ;
  705. CVTF = 'CONT' VFF ;
  706. CINT = 'INTE' CPAR CVTF ;
  707. CEXT = 'DIFF' CPAR CINT ;
  708. NPT = 'NBEL' CEXT ;
  709. B1 = CEXT 'ELEM' ('LECT' 1 'PAS' 1 NM1) ;
  710. B2 = CEXT 'ELEM' ('LECT' (NPT '-' NM1 '+' 1) 'PAS' 1 NPT) ;
  711. CEXT = 'DIFF' CEXT B1 ;
  712. CEXT = 'DIFF' CEXT B2 ;
  713.  
  714.  
  715. * Mise en place de la procedure
  716.  
  717. RXT = 'TABLE' ;
  718. RXT.'vtf' = VFF ;
  719. RXT.'breche' = AB ;
  720. RXT.'axe' = AXE ;
  721. RXT.'diru1' = (0 1) ;
  722. RXT.'pi' = ( VFF 'POIN' 'PROC' (0.0 2.0)) ;
  723. RXT.'IMPR' = 2 ;
  724.  
  725. RXT.'THERMP' = VRAI ;
  726. RXT.'vtp' = PAR ;
    4/li>8lI style="font-weight: normal; vertical-align:top;">
    RXT.'LAMBDA' = 17.0 ;
  727. *RXT.'ROCP' = (7800.0 '*' 500.0) ;
  728. RXT.'ROCP' = (78.0 '*' 500.0) ;
  729. RXT.'Tp0' = 135.0 ;
  730. RXT.'ECHAN' = 10.0 ;
  731.  
  732.  
  733. RXT.'ECHEXT' = VRAI ;
  734. RXT.'parext' = CEXT ;
  735. RXT.'HEXT' = 4.0 ;
  736. RXT.'TPEXT' = 20.0 ;
  737.  
  738.  
  739. RXT.'TIMP1' = 'TABLE' ;
  740. RXT.'TIMP1'.'MAILLAGE' = SC1 ;
  741. RXT.'TIMP1'.'t' = 'PROG' 0.0 25.5 345.5 100000.0 ;
  742. RXT.'TIMP1'.'TIMP' = 'PROG' 135.0 135.0 115.0 115.0 ;
  743. RXT.'TIMP1'.'ECHAN' = 25.0 ;
  744. RXT.'TIMP2' = 'TABLE' ;
  745. RXT.'TIMP2'.'MAILLAGE' = SC2 ;
  746. RXT.'TIMP2'.'t' = 'PROG' 0.0 25.5 345.5 100000.0 ;
  747. RXT.'TIMP2'.'TIMP' = 'PROG' 135.0 135.0 115.0 115.0 ;
  748. RXT.'TIMP2'.'ECHAN' = 25.0 ;
  749. RXT.'TIMP3' = 'TABLE' ;
  750. RXT.'TIMP3'.'MAILLAGE' = SC3 ;
  751. RXT.'TIMP3'.'t' = 'PROG' 0.0 25.5 345.5 100000.0 ;
  752. RXT.'TIMP3'.'TIMP' = 'PROG' 135.0 135.0 115.0 115.0;
  753. RXT.'TIMP3'.'ECHAN' = 25.0 ;
  754.  
  755. RXT.'VAPEUR' = VRAI ;
  756. RXT.'AIR' = VRAI ;
  757. RXT.'HE' = VRAI ;
  758. RXT.'TF0' = 140.0 ;
  759. RXT.'PT0' = 5.E5 ;
  760. RXT.'Yvap0' = 0.6138 ;
  761. RXT.'Yhe0' = 0.0 ;
  762.  
  763. RXT.'MODTURB' = 'LMEL' ;
  764. RXT.'LMEL' = 0.08 ;
  765.  
  766. RXT.'MODCOND' = 'CHIL1';
  767.  
  768. RXT.'epsi' = 1.E-5 ;
  769. RXT.'DISCR' = 'MACRO' ;
  770. RXT.'KPRE' = 'CENTREP1' ;
  771. RXT.'GRAPH' = FAUX ;
  772. RXT.'FRPREC' = 2 ;
  773.  
  774. RXT.'scenario' = 'TABLE' ;
  775. RXT.'scenario'.'t' = 'PROG' 0.0 100000.0 ;
  776. RXT.'scenario'.'qair' = 'PROG' 0.0 0.0 ;
  777. RXT.'scenario'.'qhe' = 'PROG' 0.0 0.0 ;
  778. RXT.'scenario'.'qeau' = 'PROG' 0.13 0.13 ;
  779. RXT.'scenario'.'tinj' = 'PROG' 200.0 200.0 ;
  780.  
  781.  
  782. *************************************
  783. ** EXECUTION **
  784. *************************************
  785. RXT.'DT0' = 0.05 ;
  786. NSAV = 0 ;
  787. NPAS = 10 ;
  788. ENCEINTE NPAS RXT ;
  789. NSAV = NSAV '+' NPAS ;
  790.  
  791.  
  792. RXT.'DT0' = 0.5 ;
  793. NPAS = 10 ;
  794. ENCEINTE NPAS RXT ;
  795. NSAV = NSAV '+' NPAS ;
  796.  
  797. RXT.'DT0' = 2.0 ;
  798. NPAS = 10 ;
  799. ENCEINTE NPAS RXT ;
  800. NSAV = NSAV '+' NPAS ;
  801.  
  802. RXT.'DT0' = 0.05 ;
  803. NPAS = 10 ;
  804. ENCEINTE NPAS RXT ;
  805. NSAV = NSAV '+' NPAS ;
  806.  
  807. RXT.'DT0' = 0.5 ;
  808. NPAS = 11 ;
  809. ENCEINTE NPAS RXT ;
  810. NSAV = NSAV '+' NPAS ;
  811.  
  812. RXT.'DT0' = 1. ;
  813. NPAS = 50 ;
  814. ENCEINTE NPAS RXT ;
  815. NSAV = NSAV '+' NPAS ;
  816.  
  817. *'OPTI' 'SAUVE' 'isp472D.sauv' ;
  818. *'SAUVE' ;
  819.  
  820. ***----------------------------------------***
  821. *** TESTS de NON-REGRESSION ***
  822. ***----------------------------------------***
  823.  
  824. *** Test sur les valeurs finales (regime permanant)
  825. *** de la température et de la pression
  826.  
  827. IERR = 0 ;
  828. EPSNR = 9.E-4 ;
  829. KKC2mref = 9.57867E-03;
  830. KKC2mref = 9.15854E-03;
  831. Ntps = 'DIME' rxt.'TIC'.'LTPS';
  832. ********
  833.  
  834. diag2 = 'DOMA' rxt.'GEO'.'$mtp2' 'XXDIAGSI';
  835. KKC2 = 'ELNO' rxt.'TIC'.'KKC2' rxt.'GEO'.'$mtp2';
  836. KKC2mcal= ('SOMT' (diag2 * KKC2)) / ('SOMT' diag2);
  837. 'MESS' 'KKC2mcal = ' KKC2mcal ' m/s' ;
  838. A1 = ('ABS' (KKC2mref '-' KKC2mcal)) / KKC2mref ;
  839. 'MESS' 'Erreur relative sur le coeff d echange de masse: ' A1 ;
  840. 'SI' (A1 '>EG' EPSNR);
  841. 'MESS' 'Probleme avec le coeff d echange de masse !' ;
  842. IERR = IERR '+' 1 ;
  843. 'FINSI';
  844.  
  845.  
  846. 'SI' (IERR '>' 0) ;
  847. 'MESS' 'Il y a des problemes !!!';
  848. 'ERRE' 5;
  849. 'SINON';
  850. 'MESS' 'Tout s est bien passe!' ;
  851. 'FINSI';
  852.  
  853.  
  854. 'FIN' ;
  855. *
  856.  
  857.  
  858.  
  859.  
  860.  
  861.  
  862.  
  863.  
  864.  

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