Télécharger la plaquette Cast3M


fr en

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' 'DIRE'  (-1.0 '*' V1) 'CERC' (0. RS)  
  140.  ( 0. 0. ) ; 
  141. BASE_32 = BASE_22 'PROJ' 'DIRE'  (-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. * Mise en place de la procedure 
  715. RXT = 'TABLE' ;
  716. rxt . 'VERSION' = 'V0' ;
  717. rxt . 'vtf'     = VFF ;
  718. rxt . 'axe'     = AXE ;
  719. rxt . 'epsi'    = 1.E-5 ;
  720. *
  721. rxt . 'pi' = ( VFF 'POIN' 'PROC' (0.0 2.0)) ;
  722. *
  723. rxt . 'DISCR' = 'MACRO' ;
  724. rxt . 'KPRE'  = 'CENTREP1' ;
  725. *
  726. rxt . 'MODTURB' = 'LMEL' ; 
  727. rxt . 'LMEL'    = 0.080  ;     
  728. rxt . 'MODCOND' = 'CHIL1' ;  
  729. *
  730. rxt . 'THERMP' = VRAI ;
  731. rxt . 'vtp'    = PAR ;
  732. rxt . 'LAMBDA' = 17.0 ;
  733. *ext . 'ROCP' = (7800.0 '*' 500.0) ;   
  734. rxt . 'ROCP'   = (78.0 '*' 500.0) ;   
  735. rxt . 'Tp0'    = 135.0 ;  
  736. rxt . 'ECHAN'  = 10.0 ;
  737.  
  738. rxt . 'ECHEXT' = VRAI ;
  739. rxt . 'parext' = CEXT ;
  740. rxt . 'HEXT'   = 4.0 ;
  741. rxt . 'TPEXT'  = 20.0 ;
  742.  
  743. rxt . 'TIMP' = 'TABLE' ;  
  744. rxt . 'TIMP' . 'TIMP1' = 'TABLE' ;  
  745. rxt . 'TIMP' . 'TIMP1'.'MAILLAGE' = SC1 ;  
  746. rxt . 'TIMP' . 'TIMP1'.'t' = 'PROG' 0.0 25.5 345.5 100000.0 ;  
  747. rxt . 'TIMP' . 'TIMP1'.'TIMP' = 'PROG' 135.0 135.0 115.0 115.0 ;   
  748. rxt . 'TIMP' . 'TIMP1'.'ECHAN' = 25.0 ; 
  749. rxt . 'TIMP' . 'TIMP2' = 'TABLE' ;  
  750. rxt . 'TIMP' . 'TIMP2'.'MAILLAGE' = SC2 ;  
  751. rxt . 'TIMP' . 'TIMP2'.'t' = 'PROG' 0.0 25.5 345.5 100000.0 ;  
  752. rxt . 'TIMP' . 'TIMP2'.'TIMP' = 'PROG' 135.0 135.0 115.0 115.0 ;   
  753. rxt . 'TIMP' . 'TIMP2'.'ECHAN' = 25.0 ; 
  754. rxt . 'TIMP' . 'TIMP3' = 'TABLE' ;  
  755. rxt . 'TIMP' . 'TIMP3'.'MAILLAGE' = SC3 ;  
  756. rxt . 'TIMP' . 'TIMP3'.'t' = 'PROG' 0.0 25.5 345.5 100000.0 ;  
  757. rxt . 'TIMP' . 'TIMP3'.'TIMP' = 'PROG' 135.0 135.0 115.0 115.0;
  758. rxt . 'TIMP' . 'TIMP3'.'ECHAN' = 25.0 ; 
  759.  
  760. rxt . 'VAPEUR' = VRAI ;
  761. rxt . 'HE'     = VRAI ;
  762. RXT . 'Yvap0'  = 0.6138 ;
  763. RXT . 'Yhe0'   = 0.0 ;   
  764.  
  765. RXT . 'TF0'    = 140.0 ;  
  766. RXT . 'PT0'    = 5.E5 ;   
  767. *
  768. rxt . 'Breches' = 'TABLE' ;
  769. rxt . 'Breches' . 'A' = 'TABLE' ;
  770. rxt . 'Breches' . 'A' . 'Maillage' = AB ;
  771. rxt . 'Breches' . 'A' . 'diru'  = (0. 1.) ;
  772. rxt . 'Breches' . 'A' . 'scenario' = 'TABLE' ;
  773. rxt . 'Breches' . 'A' . 'scenario' . 't'    = 'PROG' 0.0    100000.0 ;
  774. rxt . 'Breches' . 'A' . 'scenario' . 'qair' = 'PROG' 0.     0.    ;
  775. rxt . 'Breches' . 'A' . 'scenario' . 'qhe'  = 'PROG' 0.     0.    ;
  776. rxt . 'Breches' . 'A' . 'scenario' . 'qeau' = 'PROG' 0.130  0.130 ;
  777. rxt . 'Breches' . 'A' . 'scenario' . 'tinj' = 'PROG' 200.0  200.0 ;
  778. *
  779. rxt . 'IMPR' = 2 ;
  780. rxt . 'GRAPH' = FAUX ; 
  781. rxt . 'FRPREC' = 2 ;
  782.  
  783.  
  784. *************************************   
  785. **  EXECUTION  **   
  786. *************************************   
  787. RXT.'DT0' = 0.05 ;
  788. NSAV = 0 ;
  789. NPAS = 10 ;
  790. EXECRXT  NPAS RXT ;
  791. NSAV = NSAV '+' NPAS ;
  792.  
  793.  
  794. RXT.'DT0' = 0.5 ;
  795. NPAS = 10 ;
  796. EXECRXT  NPAS RXT ;
  797. NSAV = NSAV '+' NPAS ;
  798.  
  799. RXT.'DT0' = 2.0 ;
  800. NPAS = 10 ;
  801. EXECRXT  NPAS RXT ;
  802. NSAV = NSAV '+' NPAS ;
  803.  
  804. RXT.'DT0' = 0.05 ;
  805. NPAS = 10 ;
  806. EXECRXT  NPAS RXT ;
  807. NSAV = NSAV '+' NPAS ;
  808.  
  809. RXT.'DT0' = 0.5 ;
  810. NPAS = 11 ;
  811. EXECRXT  NPAS RXT ;
  812. NSAV = NSAV '+' NPAS ;
  813.  
  814. RXT.'DT0' = 1. ;
  815. NPAS = 50 ;
  816. EXECRXT  NPAS RXT ;
  817. NSAV = NSAV '+' NPAS ;
  818.  
  819. *'OPTI' 'SAUVE' 'isp472D.sauv'  ;
  820. *'SAUVE' ;
  821.  
  822. ***----------------------------------------***
  823. ***TESTS de NON-REGRESSION ***
  824. ***----------------------------------------***
  825.  
  826. *** Test sur la valeur moyenne sur la paroi
  827. *** du coefficient de transfert de masse
  828.  
  829. IERR  = 0 ;
  830. EPSNR = 9.E-4 ;
  831. KKC2mref = 9.57867E-03;
  832. KKC2mref = 9.15854E-03;
  833. KKC2mref = 4.21656E-03;
  834. KKC2mref = 4.27938E-03;
  835. KKC2mref = 4.33798E-03 ;
  836. Ntps   = 'DIME' rxt.'TIC'.'LTPS';
  837. *
  838. 'OPTI' 'ECHO' 0 ; 'SAUT' 2 'LIGN' ;
  839. $mtp2    = rxt . 'TIMP' . 'TIMP2' . '$mtpi' ;
  840. diag2    = 'DOMA' $mtp2 'XXDIAGSI';
  841. kkc2     = rxt . 'TIMP' . 'TIMP2' . 'KKC1' ;
  842. KKC2     = 'ELNO' rxt . 'TIC' . kkc2  $mtp2 ;
  843. KKC2mcal = ('SOMT' (diag2 '*' KKC2)) '/' ('SOMT' diag2) ;
  844. A1       = 'ABS' (KKC2mref '-' KKC2mcal '/' KKC2mref) ;
  845. 'MESS' 'Test valeur moyenne du coeff d echange de masse kc ' ;
  846. 'MESS' 'Err. rel. sur kc :' ' ' A1 ' '  '(tol.' ' '  epsnr')' ;
  847. 'MESS' 'Ref. :' ' ' KKC2mref ' ' 'm/s |' ' '
  848.        'Calc. :' ' ' KKC2mcal ' ' 'm/s' ;
  849. 'SI' (A1 '>EG' EPSNR) ;
  850.   IERR  = IERR '+' 1 ;
  851. 'FINSI';
  852.  
  853. 'SI' (IERR '>' 0) ;
  854.   'MESS' 'Il y a des problemes !!!' ;
  855.   'ERRE' 5 ;
  856. 'SINO';
  857.   'MESS' 'Tout s est bien passe!' ;
  858. 'FINS' ;
  859.  
  860. 'FIN' ;
  861. * 
  862.  
  863.  
  864.  
  865.  
  866.  
  867.  
  868.  
  869.  
  870.  
  871.  
  872.  
  873.  
  874.  
  875.  
  876.  
  877.  
  878.  
  879.  
  880.  
  881.  
  882.  

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