Télécharger gonfl2Dex.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : gonfl2Dex.dgibi
  2. *****************************************************************
  3. ** Jeu de données MISTRA pour le maillage 2D de l enceinte
  4. **
  5. ** Il a été choisi de mailler tous les volumes à charge
  6. ** ensuite à l utilisateur de ne retenir que ce qui
  7. ** l intéresse pour le calcul
  8. **
  9. *
  10. * Essai de gonflement de l"enceinte en AIR isotherme
  11. *
  12. ** Auteur E. Studer
  13. ** Date : 13/03/02
  14. *****************************************************************
  15. COMPLET = FAUX ;
  16. GRAPH = FAUX ;
  17.  
  18. 'SI' (COMPLET) ;
  19. NPTS = 897 ;
  20. 'SINON' ;
  21. NPTS = 10 ;
  22. 'FINSI' ;
  23.  
  24. 'OPTI' 'DIME' 2 ;
  25. 'OPTI' 'TRAC' 'PSC' ;
  26. 'OPTI' 'MODE' 'AXIS' ;
  27. 'OPTI' 'ELEM' 'QUA4' ;
  28.  
  29. * on modifie CALCP pour avoir un Cp qui ne depend pas de la
  30. * temperature pour valider la montée en pression et en temperature
  31.  
  32. DEBPROC CALCP ;
  33. 'ARGU' T/'CHPOINT'
  34. TF/'FLOTTANT'
  35. TL/'LISTREEL' ;
  36. *-------------------------------------------------*
  37. * calcul de la capacite calorifique a pression *
  38. * constante pour les incondensables en fonction *
  39. * de la temperature (J/kg/K) *
  40. * *
  41. * Entree T en Celsius *
  42. * Sortie Cp en J/kg/K *
  43. * *
  44. * PS : le resultat est du meme type que l'entree *
  45. *-------------------------------------------------*
  46. A0 = 14.346 ;
  47. A1 = 8.16D-4 ;
  48. A2 = 4.35D-7 ;
  49. B0 = 5.193 ;
  50. B1 = 0.0 ;
  51. B2 = 0.0 ;
  52. C0 = 0.917 ;
  53. C1 = 1.404D-4 ;
  54. C2 = 3.0D-8 ;
  55. D0 = 1.0309 ;
  56. D1 = 2.044D-4 ;
  57. D2 = -3.3D-8 ;
  58. E0 = 0.836 ;
  59. E1 = 7.432D-4 ;
  60. E2 = -2.511D-7 ;
  61. F0 = 1.0302 ;
  62. F1 = 2.396D-4 ;
  63. F2 = -4.950D-8 ;
  64.  
  65. 'SI' ( 'EXIST' T ) ;
  66. Cph2 = ( ((A2 * T) * T ) + ( T * A1 ) + A0 ) * 1.D3 ;
  67. Cphe = ( ((B2 * T) * T ) + ( T * B1 ) + B0 ) * 1.D3 ;
  68. Cpo2 = ( ((C2 * T) * T ) + ( T * C1 ) + C0 ) * 1.D3 ;
  69. Cpn2 = ( ((D2 * T) * T ) + ( T * D1 ) + D0 ) * 1.D3 ;
  70. Cpco2 = ( ((E2 * T) * T ) + ( T * E1 ) + E0 ) * 1.D3 ;
  71. Cpco = ( ((F2 * T) * T ) + ( T * F1 ) + F0 ) * 1.D3 ;
  72. * Cpair = (0.8*Cpn2) + (0.2*Cpo2) ;
  73. Cpair = 1000.0 ;
  74. 'RESPROC' Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair ;
  75. 'QUITTER' CALCP ;
  76. 'FINSI' ;
  77. 'SI' ( 'EXIST' TF) ;
  78. Cph2 = ( ((A2 * TF) * TF ) + ( TF * A1 ) + A0 ) * 1.D3 ;
  79. Cphe = ( ((B2 * TF) * TF ) + ( TF * B1 ) + B0 ) * 1.D3 ;
  80. Cpo2 = ( ((C2 * TF) * TF ) + ( TF * C1 ) + C0 ) * 1.D3 ;
  81. Cpn2 = ( ((D2 * TF) * TF ) + ( TF * D1 ) + D0 ) * 1.D3 ;
  82. Cpco2 = ( ((E2 * TF) * TF ) + ( TF * E1 ) + E0 ) * 1.D3 ;
  83. Cpco = ( ((F2 * TF) * TF ) + ( TF * F1 ) + F0 ) * 1.D3 ;
  84. * Cpair = (0.8*Cpn2) + (0.2*Cpo2) ;
  85. Cpair = (0.0*Cpn2) + 1000.0 ;
  86. 'RESPROC' Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair ;
  87. 'QUITTER' CALCP ;
  88. 'FINSI' ;
  89. 'SI' ( 'EXIST' TL) ;
  90. T00 = TL / TL ;
  91. A00 = A0 * T00 ;
  92. A10 = A1 * T00 ;
  93. A20 = A2 * T00 ;
  94. B00 = B0 * T00 ;
  95. B10 = B1 * T00 ;
  96. B20 = B2 * T00 ;
  97. C00 = C0 * T00 ;
  98. C10 = C1 * T00 ;
  99. C20 = C2 * T00 ;
  100. D00 = D0 * T00 ;
  101. D10 = D1 * T00 ;
  102. D20 = D2 * T00 ;
  103. E00 = E0 * T00 ;
  104. E10 = E1 * T00 ;
  105. E20 = E2 * T00 ;
  106. F00 = F0 * T00 ;
  107. F10 = F1 * T00 ;
  108. F20 = F2 * T00 ;
  109.  
  110. Cph2 = ( ((A20* TL) * TL ) + ( TL * A10) + A00) * 1.D3 ;
  111. Cphe = ( ((B20* TL) * TL ) + ( TL * B10) + B00) * 1.D3 ;
  112. Cpo2 = ( ((C20* TL) * TL ) + ( TL * C10) + C00) * 1.D3 ;
  113. Cpn2 = ( ((D20* TL) * TL ) + ( TL * D10) + D00) * 1.D3 ;
  114. Cpco2 = ( ((E20* TL) * TL ) + ( TL * E10) + E00) * 1.D3 ;
  115. Cpco = ( ((F20* TL) * TL ) + ( TL * F10) + F00) * 1.D3 ;
  116. * Cpair = (0.8*Cpn2) + (0.2*Cpo2) ;
  117. Cpair = (0.0*Cpn2) + 1000.0 ;
  118. 'RESPROC' Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair ;
  119. 'QUITTER' CALCP ;
  120.  
  121. 'FINSI' ;
  122. 'ERRE' 5 ;
  123. FINPROC ;
  124.  
  125. CTRV = FAUX ;
  126.  
  127. DMIN = 0.05 ;
  128. DMIN1 = 0.02 ;
  129. DMAX = 0.2 ;
  130.  
  131. * Dimensions sur plan
  132. R0 = 2.125 ;
  133. R0S2 = R0 '/' 2.0 ;
  134.  
  135. EPB = 0.025 ;
  136. EPF = 0.015 ;
  137. EPH = 0.119 ;
  138. EPSI = 1.E-5 ;
  139.  
  140. HCOND1 = 2.187 ;
  141. HCOND2 = 1.784 ;
  142. HCOND3 = 1.784 ;
  143.  
  144. REQC = 1.9079 ;
  145. EPC = 0.040 ;
  146. EPCC = 0.120 ;
  147. HCC = 0.150 ;
  148. HEC = 0.120 ;
  149.  
  150.  
  151. EPBR = 0.110 ;
  152. HBR1 = 0.190 ;
  153. HBR2 = 0.195 ;
  154.  
  155. RINJ = 0.1 ;
  156. RS = 3.2529 ;
  157.  
  158. ZCOND1 = 1.285 ;
  159. ZCOND2 = 3.592 ;
  160. ZCOND3 = 5.4965 ;
  161.  
  162. ZHB1 = 1.938 ;
  163. ZHB2 = 4.341 ;
  164. ZHB3 = 6.559 ;
  165.  
  166. * Entiers pour le nombre de mailles
  167. NINJ = 2 ;
  168. NX1 = 6 ;
  169. NX2 = 2 ;
  170. NX3 = 3 ;
  171. NX4 = 3 ;
  172.  
  173. NF5 = 2 ;
  174. NM1 = 2 ;
  175. NM2 = 4 ;
  176. NZ1 = 2 ;
  177. NZ2 = 6 ;
  178. NZ3 = 2 ;
  179. *NZ4 = 5 ;
  180. *NZ5 = 2 ;
  181. NZ4 = 12 ;
  182. NZ5 = 4 ;
  183.  
  184. * Injection
  185. A = 0.0 ZCOND1 ;
  186. B = RINJ ZCOND1 ;
  187. AB = A 'DROI' NINJ B ;
  188. AB = AB 'COUL' 'BLAN' ;
  189.  
  190. * jusqu"a R/2
  191. C = R0S2 ZCOND1 ;
  192. BC = B 'DROI' C 'DINI' DMIN 'DFIN' DMAX ;
  193. BC = BC 'COUL' 'BLAN' ;
  194.  
  195. * jusqu"au condenseurs
  196. D = REQC ZCOND1 ;
  197. CD = C 'DROI' D 'DINI' DMAX 'DFIN' DMIN1 ;
  198. CD = CD 'COUL' 'BLAN' ;
  199.  
  200. * condenseurs plus isolant
  201. E = (REQC '+' EPC) ZCOND1 ;
  202. DE = D 'DROI' NX2 E ;
  203. DE = DE 'COUL' 'BLEU' ;
  204.  
  205. * collecteurs des condenseurs
  206. F = (REQC '+' EPC '+' EPCC) ZCOND1 ;
  207. EF = E 'DROI' NX3 F ;
  208. EF = EF 'COUL' 'JAUN' ;
  209.  
  210. * volume mort derriere les condenseurs
  211. G = R0 ZCOND1 ;
  212. FG = F 'DROI' NX4 G ;
  213. FG = FG 'COUL' 'BLAN' ;
  214.  
  215. * fut de l"enceinte
  216. H = (R0 '+' EPF) ZCOND1 ;
  217. GH = G 'DROI' NM1 H ;
  218. GH = GH 'COUL' 'ROUG' ;
  219.  
  220. * Brides
  221. I = (R0 '+' EPF '+' EPBR) ZCOND1 ;
  222. HI = H 'DROI' NM2 I ;
  223. HI = HI 'COUL' 'VERT' ;
  224.  
  225. * Passage au maillage 2D
  226. BASE = BC 'ET' CD 'ET' (DE 'COUL' 'BLAN') 'ET'
  227. (EF 'COUL' 'BLAN') 'ET' FG 'ET' GH 'ET' HI ;
  228. BASE_1 = BASE 'PLUS' (0. (HBR1 '*' -1.0)) ;
  229. VB_1 = 'TRAN' NZ1 BASE (0. (HBR1 '*' -1.0)) ;
  230. 'ELIM' EPSI (BASE_1 'ET' VB_1) ;
  231.  
  232. 'SI' GRAPH ;
  233. 'TRACE' 'CACH' (VB_1)
  234. 'TITR' 'Volume VB_1 ' ;
  235. 'FINSI' ;
  236.  
  237. V1 = 0. ZCOND1 ;
  238. BASE_21 = (BASE_1 'ELEM' 'BLAN') ;
  239. BASE_22 = (BASE_1 'ELEM' 'ROUG') ;
  240. BASE_31 = BASE_21 'PROJ' 'CYLI' (-1.0 '*' V1) 'CERC' (0. RS)
  241. ( 0. 0. ) ;
  242. BASE_32 = BASE_22 'PROJ' 'CYLI' (-1.0 '*' V1) 'CERC' (0. RS)
  243. ( 0. 0. ) ;
  244. BASE_3 = BASE_31 'ET' BASE_32 ;
  245. 'ELIM' EPSI BASE_3 ;
  246.  
  247. VB_2_1 = 'REGLER' BASE_21 BASE_31 NZ2 ;
  248. VB_2_2 = 'REGLER' BASE_22 BASE_32 NZ2 ;
  249. VB_2 = VB_2_1 'ET' VB_2_2 ;
  250. 'ELIM' EPSI VB_2 ;
  251.  
  252. 'SI' GRAPH ;
  253. 'TRACE' 'CACH' (VB_1 'ET' VB_2)
  254. 'TITR' 'Volume VB_1,VB_2 ' ;
  255. 'FINSI' ;
  256.  
  257. BASE_4 = BASE_3 'HOMO' (0. RS) (RS '+' EPB '/' RS) ;
  258. VB_3 = 'REGLER' BASE_3 BASE_4 NM1 ;
  259. VB_3 = VB_3 'COUL' 'ROUG' ;
  260. 'SI' GRAPH ;
  261. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3)
  262. 'TITR' 'Volume VB_1,VB_2,VB_3 ' ;
  263. 'FINSI' ;
  264.  
  265. * Maintenant on monte en Z
  266. *
  267. * Collecteur
  268. HAUT = AB 'ET' BC 'ET' CD 'ET' DE 'ET' EF 'ET' FG 'ET'
  269. GH 'ET' HI ;
  270. HAUT_2 = HAUT 'PLUS' (0. HCC) ;
  271. VB_4 = 'TRAN' NZ3 HAUT (0. HCC) ;
  272. 'ELIM' EPSI (HAUT_2 'ET' VB_4) ;
  273. H_V = VB_4 'ELEM' 'VERT' ;
  274. VB_4 = 'DIFF' VB_4 H_V ;
  275. SC1_1 = INTER ('CONT' (VB_4 'ELEM' 'BLAN'))
  276. ('CONT' (VB_4 'ELEM' 'BLEU'));
  277. * il faut construire l"axe
  278. CNT = 'CONT' VB_4 ;
  279. CNT = 'DIFF' CNT HAUT ;
  280. CNT = 'DIFF' CNT HAUT_2 ;
  281. CNT = 'CCON' CNT ;
  282. AXE = CNT.2 ;
  283. 'ELIM' EPSI (AXE 'ET' VB_4) ;
  284.  
  285. 'SI' GRAPH ;
  286. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4)
  287. 'TITR' 'Volume VB_1,VB_2,VB_3,VB_4 ' ;
  288. 'FINSI' ;
  289.  
  290. *Bride et le collecteur devient fluide (BLAN)
  291. HAUT_3 = HAUT_2 'PLUS' (0. HBR2) ;
  292. VB_5 = 'TRAN' NZ1 HAUT_2 (0. HBR2) ;
  293. 'ELIM' EPSI (HAUT_3 'ET' VB_5) ;
  294. SC1_2 = INTER ('CONT' (VB_5 'ELEM' 'BLAN'))
  295. ('CONT' (VB_5 'ELEM' 'BLEU'));
  296. H1_J = VB_5 'ELEM' 'JAUN' ;
  297. H1_J = H1_J 'COUL' 'BLAN' ;
  298. VB_5 = 'DIFF' VB_5 H1_J ;
  299. VB_5 = VB_5 'ET' H1_J ;
  300.  
  301. CNT = 'CONT' VB_5 ;
  302. CNT = 'DIFF' CNT HAUT_2 ;
  303. CNT = 'DIFF' CNT HAUT_3 ;
  304. CNT = 'CCON' CNT ;
  305. AXE1 = CNT.2 ;
  306. 'ELIM' EPSI (AXE1 'ET' VB_5) ;
  307. AXE = AXE 'ET' AXE1 ;
  308. 'SI' GRAPH ;
  309. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  310. 'ET' VB_5 )
  311. 'TITR' 'Volume VB_1,VB_2,VB_3,VB_4,VB_5 ' ;
  312. 'FINSI' ;
  313.  
  314. * entre les deux brides on ajoute un plan pour le hublot
  315. Z3 = ( 'MAXI' ('COOR' HAUT_3 2)) ;
  316. HAUT_3M = HAUT_3 'PLUS' (0. (ZHB1 '-' Z3)) ;
  317. HAUT_4 = HAUT_3 'PLUS'
  318. (0. (HCOND1 '-' HBR2 '-' HBR2 '-' HCC '-' HCC )) ;
  319. VB_6_1 = 'TRAN' NZ5 HAUT_3 (0. (ZHB1 '-' Z3)) ;
  320. VB_6_2 = 'TRAN' NZ4 HAUT_3M
  321. (0. (HCOND1 '-' HBR2 '-' HBR2 '-' HCC '-' HCC '-' ZHB1 '+' Z3)) ;
  322. VB_6 = VB_6_1 'ET' VB_6_2 ;
  323. 'ELIM' EPSI (HAUT_3M 'ET' HAUT_4 'ET' VB_6) ;
  324.  
  325. SC1_3 = INTER ('CONT' (VB_6 'ELEM' 'BLAN'))
  326. ('CONT' (VB_6 'ELEM' 'BLEU'));
  327. H2_V = VB_6 'ELEM' 'VERT' ;
  328. VB_6 = 'DIFF' VB_6 H2_V ;
  329. H1_J = VB_6 'ELEM' 'JAUN' ;
  330. H1_J = H1_J 'COUL' 'BLAN' ;
  331. VB_6 = 'DIFF' VB_6 H1_J ;
  332. VB_6 = VB_6 'ET' H1_J ;
  333.  
  334. CNT = 'CONT' VB_6 ;
  335. CNT = 'DIFF' CNT HAUT_3 ;
  336. CNT = 'DIFF' CNT HAUT_4 ;
  337. CNT = 'CCON' CNT ;
  338. AXE1 = CNT.2 ;
  339. 'ELIM' EPSI (AXE1 'ET' VB_6) ;
  340. AXE = AXE 'ET' AXE1 ;
  341. 'SI' GRAPH ;
  342. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  343. 'ET' VB_5 'ET' VB_6 )
  344. 'TITR' 'Volume VB_1,VB_2,VB_3,VB_4,VB_5,VB_6 ' ;
  345. 'FINSI' ;
  346.  
  347. * une nouvelle bride
  348. HAUT_5 = HAUT_4 'PLUS' (0. HBR2) ;
  349. VB_7 = 'TRAN' NZ1 HAUT_4 (0. HBR2) ;
  350. 'ELIM' EPSI (HAUT_5 'ET' VB_7) ;
  351. SC1_4 = INTER ('CONT' (VB_7 'ELEM' 'BLAN'))
  352. ('CONT' (VB_7 'ELEM' 'BLEU'));
  353. H1_J = VB_7 'ELEM' 'JAUN' ;
  354. H1_J = H1_J 'COUL' 'BLAN' ;
  355. VB_7 = 'DIFF' VB_7 H1_J ;
  356. VB_7 = VB_7 'ET' H1_J ;
  357. CNT = 'CONT' VB_7 ;
  358. CNT = 'DIFF' CNT HAUT_4 ;
  359. CNT = 'DIFF' CNT HAUT_5 ;
  360. CNT = 'CCON' CNT ;
  361. AXE1 = CNT.2 ;
  362. 'ELIM' EPSI (AXE1 'ET' VB_7) ;
  363. AXE = AXE 'ET' AXE1 ;
  364.  
  365. 'SI' GRAPH ;
  366. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  367. 'ET' VB_5 'ET' VB_6 'ET' VB_7 )
  368. 'TITR' 'VB_1,VB_2,VB_3,VB_4,VB_5,VB_6, VB_7 ' ;
  369. 'FINSI' ;
  370.  
  371. * un collecteur haut du condenseur bas
  372. HAUT_6 = HAUT_5 'PLUS' (0. HCC) ;
  373. VB_8 = 'TRAN' NZ1 HAUT_5 (0. HCC) ;
  374. 'ELIM' EPSI (HAUT_6 'ET' VB_8) ;
  375. SC1_5 = INTER ('CONT' (VB_8 'ELEM' 'BLAN'))
  376. ('CONT' (VB_8 'ELEM' 'BLEU'));
  377. H2_V = VB_8 'ELEM' 'VERT' ;
  378. VB_8 = 'DIFF' VB_8 H2_V ;
  379. CNT = 'CONT' VB_8 ;
  380. CNT = 'DIFF' CNT HAUT_5 ;
  381. CNT = 'DIFF' CNT HAUT_6 ;
  382. CNT = 'CCON' CNT ;
  383. AXE1 = CNT.2 ;
  384. 'ELIM' EPSI (AXE1 'ET' VB_8) ;
  385. AXE = AXE 'ET' AXE1 ;
  386.  
  387. 'SI' GRAPH ;
  388. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  389. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 )
  390. 'TITR' 'VB1,VB2,VB3,VB4,VB5,VB6,VB7,VB8 ' ;
  391. 'FINSI' ;
  392.  
  393. * entre le condenseur 1 et le condenseur 2
  394. HAUT_7 = HAUT_6 'PLUS' (0. HEC) ;
  395. VB_9 = 'TRAN' NZ1 HAUT_6 (0. HEC) ;
  396. 'ELIM' EPSI (HAUT_7 'ET' VB_9) ;
  397. H2_V = VB_9 'ELEM' 'VERT' ;
  398. VB_9 = 'DIFF' VB_9 H2_V ;
  399. H1_J = VB_9 'ELEM' 'JAUN' ;
  400. H1_J = H1_J 'COUL' 'BLAN' ;
  401. VB_9 = 'DIFF' VB_9 H1_J ;
  402. VB_9 = VB_9 'ET' H1_J ;
  403. H1_B = VB_9 'ELEM' 'BLEU' ;
  404. H1_B = H1_B 'COUL' 'BLAN' ;
  405. VB_9 = 'DIFF' VB_9 H1_B ;
  406. VB_9 = VB_9 'ET' H1_B ;
  407. CNT = 'CONT' VB_9 ;
  408. CNT = 'DIFF' CNT HAUT_6 ;
  409. CNT = 'DIFF' CNT HAUT_7 ;
  410. CNT = 'CCON' CNT ;
  411. AXE1 = CNT.2 ;
  412. 'ELIM' EPSI (AXE1 'ET' VB_9) ;
  413. AXE = AXE 'ET' AXE1 ;
  414. 'SI' GRAPH ;
  415. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  416. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9 )
  417. 'TITR' 'VB1,VB2,VB3,VB4,VB5,VB6,VB7,VB8,VB9 ' ;
  418. 'FINSI' ;
  419.  
  420. * collecteur bas du condenseur moyen
  421. HAUT_8 = HAUT_7 'PLUS' (0. HCC) ;
  422. VB_10 = 'TRAN' NZ1 HAUT_7 (0. HCC) ;
  423. 'ELIM' EPSI (HAUT_8 'ET' VB_10) ;
  424. SC2_1 = INTER ('CONT' (VB_10 'ELEM' 'BLAN'))
  425. ('CONT' (VB_10 'ELEM' 'BLEU'));
  426. H2_V = VB_10 'ELEM' 'VERT' ;
  427. VB_10 = 'DIFF' VB_10 H2_V ;
  428. CNT = 'CONT' VB_10 ;
  429. CNT = 'DIFF' CNT HAUT_7 ;
  430. CNT = 'DIFF' CNT HAUT_8 ;
  431. CNT = 'CCON' CNT ;
  432. AXE1 = CNT.2 ;
  433. 'ELIM' EPSI (AXE1 'ET' VB_10) ;
  434. AXE = AXE 'ET' AXE1 ;
  435. 'SI' GRAPH ;
  436. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  437. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  438. 'ET' VB_10 )
  439. 'TITR' 'VB1,VB2,VB3,VB4,VB5,VB6,VB7,VB8,VB9,VB10 ' ;
  440. 'FINSI' ;
  441.  
  442. * Bride du condenseur 2
  443. HAUT_9 = HAUT_8 'PLUS' (0. HBR2) ;
  444. VB_11 = 'TRAN' NZ1 HAUT_8 (0. HBR2) ;
  445. 'ELIM' EPSI (HAUT_9 'ET' VB_11) ;
  446. SC2_2 = INTER ('CONT' (VB_11 'ELEM' 'BLAN'))
  447. ('CONT' (VB_11 'ELEM' 'BLEU'));
  448. H1_J = VB_11 'ELEM' 'JAUN' ;
  449. H1_J = H1_J 'COUL' 'BLAN' ;
  450. VB_11 = 'DIFF' VB_11 H1_J ;
  451. VB_11 = VB_11 'ET' H1_J ;
  452. CNT = 'CONT' VB_11 ;
  453. CNT = 'DIFF' CNT HAUT_8 ;
  454. CNT = 'DIFF' CNT HAUT_9 ;
  455. CNT = 'CCON' CNT ;
  456. AXE1 = CNT.2 ;
  457. 'ELIM' EPSI (AXE1 'ET' VB_11) ;
  458. AXE = AXE 'ET' AXE1 ;
  459.  
  460. 'SI' GRAPH ;
  461. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  462. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  463. 'ET' VB_10 'ET' VB_11 )
  464. 'TITR' 'VB1 A VB11 ' ;
  465. 'FINSI' ;
  466.  
  467. * entre bride et collecteur haut pour condenseur moyen
  468. Z3 = ( 'MAXI' ('COOR' HAUT_9 2)) ;
  469. HAUT_9M = HAUT_9 'PLUS' (0. (ZHB2 '-' Z3)) ;
  470. HAUT_10 = HAUT_9 'PLUS'
  471. (0. (HCOND2 '-' HBR2 '-' HCC '-' HCC )) ;
  472. VB_12_1 = 'TRAN' NZ5 HAUT_9 (0. (ZHB2 '-' Z3)) ;
  473. VB_12_2 = 'TRAN' (NZ4 '-' NZ5) HAUT_9M
  474. (0. (HCOND2 '-' HBR2 '-' HCC '-' HCC '-' ZHB2 '+' Z3)) ;
  475. VB_12 = VB_12_1 'ET' VB_12_2 ;
  476. 'ELIM' EPSI (HAUT_9M 'ET' HAUT_10 'ET' VB_12) ;
  477. SC2_3 = INTER ('CONT' (VB_12 'ELEM' 'BLAN'))
  478. ('CONT' (VB_12 'ELEM' 'BLEU'));
  479. H2_V = VB_12 'ELEM' 'VERT' ;
  480. VB_12 = 'DIFF' VB_12 H2_V ;
  481. H1_J = VB_12 'ELEM' 'JAUN' ;
  482. H1_J = H1_J 'COUL' 'BLAN' ;
  483. VB_12 = 'DIFF' VB_12 H1_J ;
  484. VB_12 = VB_12 'ET' H1_J ;
  485. CNT = 'CONT' VB_12 ;
  486. CNT = 'DIFF' CNT HAUT_9 ;
  487. CNT = 'DIFF' CNT HAUT_10 ;
  488. CNT = 'CCON' CNT ;
  489. AXE1 = CNT.2 ;
  490. 'ELIM' EPSI (AXE1 'ET' VB_12) ;
  491. AXE = AXE 'ET' AXE1 ;
  492. 'SI' GRAPH ;
  493. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  494. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  495. 'ET' VB_10 'ET' VB_11 'ET' VB_12 )
  496. 'TITR' 'VB1 A VB12 ' ;
  497. 'FINSI' ;
  498.  
  499. * collecteur haut du condenseur moyen
  500. HAUT_11 = HAUT_10 'PLUS' (0. HCC) ;
  501. VB_13 = 'TRAN' NZ1 HAUT_10 (0. HCC) ;
  502. 'ELIM' EPSI (HAUT_11 'ET' VB_13) ;
  503. SC2_4 = INTER ('CONT' (VB_13 'ELEM' 'BLAN'))
  504. ('CONT' (VB_13 'ELEM' 'BLEU'));
  505. H2_V = VB_13 'ELEM' 'VERT' ;
  506. VB_13 = 'DIFF' VB_13 H2_V ;
  507. CNT = 'CONT' VB_13 ;
  508. CNT = 'DIFF' CNT HAUT_10 ;
  509. CNT = 'DIFF' CNT HAUT_11 ;
  510. CNT = 'CCON' CNT ;
  511. AXE1 = CNT.2 ;
  512. 'ELIM' EPSI (AXE1 'ET' VB_13) ;
  513. AXE = AXE 'ET' AXE1 ;
  514. 'SI' GRAPH ;
  515. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  516. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  517. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 )
  518. 'TITR' 'VB1 A VB13 ' ;
  519. 'FINSI' ;
  520.  
  521. * entre le condenseur moyen et le condenseur haut
  522. HAUT_12 = HAUT_11 'PLUS' (0. HEC) ;
  523. VB_14 = 'TRAN' NZ1 HAUT_11 (0. HEC) ;
  524. 'ELIM' EPSI (HAUT_12 'ET' VB_14) ;
  525. H2_V = VB_14 'ELEM' 'VERT' ;
  526. VB_14 = 'DIFF' VB_14 H2_V ;
  527. H1_J = VB_14 'ELEM' 'JAUN' ;
  528. H1_J = H1_J 'COUL' 'BLAN' ;
  529. VB_14 = 'DIFF' VB_14 H1_J ;
  530. VB_14 = VB_14 'ET' H1_J ;
  531. H1_B = VB_14 'ELEM' 'BLEU' ;
  532. H1_B = H1_B 'COUL' 'BLAN' ;
  533. VB_14 = 'DIFF' VB_14 H1_B ;
  534. VB_14 = VB_14 'ET' H1_B ;
  535. CNT = 'CONT' VB_14 ;
  536. CNT = 'DIFF' CNT HAUT_11 ;
  537. CNT = 'DIFF' CNT HAUT_12 ;
  538. CNT = 'CCON' CNT ;
  539. AXE1 = CNT.2 ;
  540. 'ELIM' EPSI (AXE1 'ET' VB_14) ;
  541. AXE = AXE 'ET' AXE1 ;
  542. 'SI' GRAPH ;
  543. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  544. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  545. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  546. VB_14 )
  547. 'TITR' 'VB1 A VB14 ' ;
  548. 'FINSI' ;
  549.  
  550. * collecteur bas du condenseur haut
  551. HAUT_13 = HAUT_12 'PLUS' (0. HCC) ;
  552. VB_15 = 'TRAN' NZ1 HAUT_12 (0. HCC) ;
  553. 'ELIM' EPSI (HAUT_13 'ET' VB_15) ;
  554. SC3_1 = INTER ('CONT' (VB_15 'ELEM' 'BLAN'))
  555. ('CONT' (VB_15 'ELEM' 'BLEU'));
  556. H2_V = VB_15 'ELEM' 'VERT' ;
  557. VB_15 = 'DIFF' VB_15 H2_V ;
  558. CNT = 'CONT' VB_15 ;
  559. CNT = 'DIFF' CNT HAUT_12 ;
  560. CNT = 'DIFF' CNT HAUT_13 ;
  561. CNT = 'CCON' CNT ;
  562. AXE1 = CNT.2 ;
  563. 'ELIM' EPSI (AXE1 'ET' VB_15) ;
  564. AXE = AXE 'ET' AXE1 ;
  565. 'SI' GRAPH ;
  566. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  567. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  568. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  569. VB_14 'ET' VB_15 )
  570. 'TITR' 'VB1 A VB15 ' ;
  571. 'FINSI' ;
  572.  
  573. * entre collecteur bas et bas de la bride sup
  574. Z3 = ( 'MAXI' ('COOR' HAUT_13 2)) ;
  575. HAUT_13M = HAUT_13 'PLUS' (0. (ZHB3 '-' Z3)) ;
  576. HAUT_14 = HAUT_13 'PLUS'
  577. (0. (HCOND3 '-' HBR2 '-' HCC '-' HCC )) ;
  578. VB_16_1 = 'TRAN' (NZ4 '-' NZ5) HAUT_13 (0. (ZHB3 '-' Z3)) ;
  579. VB_16_2 = 'TRAN' NZ5 HAUT_13M
  580. (0. (HCOND3 '-' HBR2 '-' HCC '-' HCC '-' ZHB3 '+' Z3)) ;
  581. VB_16 = VB_16_1 'ET' VB_16_2 ;
  582. 'ELIM' EPSI (HAUT_13M 'ET' HAUT_14 'ET' VB_16) ;
  583. SC3_2 = INTER ('CONT' (VB_16 'ELEM' 'BLAN'))
  584. ('CONT' (VB_16 'ELEM' 'BLEU'));
  585. H2_V = VB_16 'ELEM' 'VERT' ;
  586. VB_16 = 'DIFF' VB_16 H2_V ;
  587. H1_J = VB_16 'ELEM' 'JAUN' ;
  588. H1_J = H1_J 'COUL' 'BLAN' ;
  589. VB_16 = 'DIFF' VB_16 H1_J ;
  590. VB_16 = VB_16 'ET' H1_J ;
  591. CNT = 'CONT' VB_16 ;
  592. CNT = 'DIFF' CNT HAUT_13 ;
  593. CNT = 'DIFF' CNT HAUT_14 ;
  594. CNT = 'CCON' CNT ;
  595. AXE1 = CNT.2 ;
  596. 'ELIM' EPSI (AXE1 'ET' VB_16) ;
  597. AXE = AXE 'ET' AXE1 ;
  598. 'SI' GRAPH ;
  599. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  600. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  601. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  602. VB_14 'ET' VB_15 'ET' VB_16 )
  603. 'TITR' 'VB1 A VB16 ' ;
  604. 'FINSI' ;
  605.  
  606. * une nouvelle bride
  607. HAUT_15 = HAUT_14 'PLUS' (0. HBR2) ;
  608. VB_17 = 'TRAN' NZ1 HAUT_14 (0. HBR2) ;
  609. 'ELIM' EPSI (HAUT_15 'ET' VB_17) ;
  610. SC3_3 = INTER ('CONT' (VB_17 'ELEM' 'BLAN'))
  611. ('CONT' (VB_17 'ELEM' 'BLEU'));
  612. H1_J = VB_17 'ELEM' 'JAUN' ;
  613. H1_J = H1_J 'COUL' 'BLAN' ;
  614. VB_17 = 'DIFF' VB_17 H1_J ;
  615. VB_17 = VB_17 'ET' H1_J ;
  616. CNT = 'CONT' VB_17 ;
  617. CNT = 'DIFF' CNT HAUT_14 ;
  618. CNT = 'DIFF' CNT HAUT_15 ;
  619. CNT = 'CCON' CNT ;
  620. AXE1 = CNT.2 ;
  621. 'ELIM' EPSI (AXE1 'ET' VB_17) ;
  622. AXE = AXE 'ET' AXE1 ;
  623. 'SI' GRAPH ;
  624. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  625. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  626. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  627. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 )
  628. 'TITR' 'VB1 A VB17 ' ;
  629. 'FINSI' ;
  630.  
  631. * un collecteur haut du condenseur haut
  632. HAUT_16 = HAUT_15 'PLUS' (0. HCC) ;
  633. VB_18 = 'TRAN' NZ1 HAUT_15 (0. HCC) ;
  634. 'ELIM' EPSI (HAUT_16 'ET' VB_18) ;
  635. SC3_4 = INTER ('CONT' (VB_18 'ELEM' 'BLAN'))
  636. ('CONT' (VB_18 'ELEM' 'BLEU'));
  637. H2_V = VB_18 'ELEM' 'VERT' ;
  638. VB_18 = 'DIFF' VB_18 H2_V ;
  639. CNT = 'CONT' VB_18 ;
  640. CNT = 'DIFF' CNT HAUT_15 ;
  641. CNT = 'DIFF' CNT HAUT_16 ;
  642. CNT = 'CCON' CNT ;
  643. AXE1 = CNT.2 ;
  644. 'ELIM' EPSI (AXE1 'ET' VB_18) ;
  645. AXE = AXE 'ET' AXE1 ;
  646. 'SI' GRAPH ;
  647. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  648. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  649. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  650. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 'ET' VB_18 )
  651. 'TITR' 'VB1 A VB18 ' ;
  652. 'FINSI' ;
  653.  
  654. * entre le condenseur haut et le plafond
  655. HAUT_17 = HAUT_16 'PLUS' (0. 0.100) ;
  656. VB_19 = 'TRAN' NZ1 HAUT_16 (0. 0.100) ;
  657. 'ELIM' EPSI (HAUT_17 'ET' VB_19) ;
  658. H2_V = VB_19 'ELEM' 'VERT' ;
  659. VB_19 = 'DIFF' VB_19 H2_V ;
  660. H1_J = VB_19 'ELEM' 'JAUN' ;
  661. H1_J = H1_J 'COUL' 'BLAN' ;
  662. VB_19 = 'DIFF' VB_19 H1_J ;
  663. VB_19 = VB_19 'ET' H1_J ;
  664. H1_B = VB_19 'ELEM' 'BLEU' ;
  665. H1_B = H1_B 'COUL' 'BLAN' ;
  666. VB_19 = 'DIFF' VB_19 H1_B ;
  667. VB_19 = VB_19 'ET' H1_B ;
  668. CNT = 'CONT' VB_19 ;
  669. CNT = 'DIFF' CNT HAUT_16 ;
  670. CNT = 'DIFF' CNT HAUT_17 ;
  671. CNT = 'CCON' CNT ;
  672. AXE1 = CNT.2 ;
  673. 'ELIM' EPSI (AXE1 'ET' VB_19) ;
  674. AXE = AXE 'ET' AXE1 ;
  675. 'SI' GRAPH ;
  676. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  677. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  678. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  679. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 'ET' VB_18
  680. 'ET' VB_19 )
  681. 'TITR' 'VB1 A VB19' ;
  682. 'FINSI' ;
  683.  
  684. * le plafond
  685. HAUT_18 = HAUT_17 'PLUS' (0. EPH) ;
  686. VB_20 = 'TRAN' NM1 HAUT_17 (0. EPH) ;
  687. 'ELIM' EPSI (HAUT_18 'ET' VB_20) ;
  688. H2_V = VB_20 'ELEM' 'VERT' ;
  689. VB_20 = 'DIFF' VB_20 H2_V ;
  690. H1_J = VB_20 'ELEM' 'JAUN' ;
  691. H1_J = H1_J 'COUL' 'ROUG' ;
  692. VB_20 = 'DIFF' VB_20 H1_J ;
  693. VB_20 = VB_20 'ET' H1_J ;
  694. H1_B = VB_20 'ELEM' 'BLEU' ;
  695. H1_B = H1_B 'COUL' 'ROUG' ;
  696. VB_20 = 'DIFF' VB_20 H1_B ;
  697. VB_20 = VB_20 'ET' H1_B ;
  698. H1_BL = VB_20 'ELEM' 'BLAN' ;
  699. H1_BL = H1_BL 'COUL' 'ROUG' ;
  700. VB_20 = 'DIFF' VB_20 H1_BL ;
  701. VB_20 = VB_20 'ET' H1_BL ;
  702. 'SI' GRAPH ;
  703. 'TRACE' 'CACH' (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  704. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  705. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  706. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 'ET' VB_18
  707. 'ET' VB_19 'ET' VB_20 )
  708. 'TITR' 'VB1 A VB20' ;
  709. 'FINSI' ;
  710.  
  711. ***************************************************************
  712. * On effectue des controles par rapport aux plans usines
  713. ***************************************************************
  714. MISTRA = (VB_1 'ET' VB_2 'ET' VB_3 'ET' VB_4
  715. 'ET' VB_5 'ET' VB_6 'ET' VB_7 'ET' VB_8 'ET' VB_9
  716. 'ET' VB_10 'ET' VB_11 'ET' VB_12 'ET' VB_13 'ET'
  717. VB_14 'ET' VB_15 'ET' VB_16 'ET' VB_17 'ET' VB_18
  718. 'ET' VB_19 'ET' VB_20 )
  719. ;
  720.  
  721. 'SI' CTRV ;
  722.  
  723. * volume du fluide
  724. VFF = MISTRA 'ELEM' 'BLAN' ;
  725. 'TRACE' 'CACH' VFF 'TITR' 'Volume fluide' ;
  726. MVFF = 'CHAN' VFF 'QUAF' ;
  727. $VFF = 'MODE' MVFF 'NAVIER_STOKES' 'QUAF' ;
  728. VVFF = 'DOMA' $VFF 'XXDIAGSI' ;
  729. VFF = 'DOMA' $VFF 'MAILLAGE' ;
  730. 'MESS' '==> Volume du fluide : '
  731. ('SOMT' VVFF) ;
  732.  
  733. * Elevation
  734. ZZ = 'COOR' VFF 2 ;
  735. 'MESS' '==> Zfluide min =' ('MINI' ZZ) ;
  736. 'MESS' '==> Zfluide max =' ('MAXI' ZZ) ;
  737.  
  738. * Surfaces de condensation
  739. SC1 = SC1_1 'ET' SC1_2 'ET' SC1_3 'ET' SC1_4 'ET'
  740. SC1_5 ;
  741. SC2 = SC2_1 'ET' SC2_2 'ET' SC2_3 'ET' SC2_4 ;
  742. SC3 = SC3_1 'ET' SC3_2 'ET' SC3_3 'ET' SC3_4 ;
  743. 'TRACE' 'CACH' (SC1 'ET' SC2 'ET' SC3)
  744. 'TITR' 'Surface des condenseurs' ;
  745. MS1 = 'CHAN' SC1 'QUAF' ;
  746. $S1 = 'MODE' MS1 'NAVIER_STOKES' 'QUAF' ;
  747. VS1 = 'DOMA' $S1 'XXDIAGSI' ;
  748. MS2 = 'CHAN' SC2 'QUAF' ;
  749. $S2 = 'MODE' MS2 'NAVIER_STOKES' 'QUAF' ;
  750. VS2 = 'DOMA' $S2 'XXDIAGSI' ;
  751. MS3 = 'CHAN' SC3 'QUAF' ;
  752. $S3 = 'MODE' MS3 'NAVIER_STOKES' 'QUAF' ;
  753. VS3 = 'DOMA' $S3 'XXDIAGSI' ;
  754. 'MESS' '==> Surface condenseur BAS : '
  755. ('SOMT' VS1) ;
  756. 'MESS' '==> Surface condenseur MOY : '
  757. ('SOMT' VS2) ;
  758. 'MESS' '==> Surface condenseur HAU : '
  759. ('SOMT' VS3) ;
  760.  
  761. * Volume des condenseurs
  762. CS = (MISTRA 'ELEM' 'BLEU') 'ET' (MISTRA 'ELEM' 'JAUN') ;
  763. 'TRACE' 'CACH' CS 'TITR' 'Volume des condenseurs' ;
  764. MCS = 'CHAN' CS 'QUAF' ;
  765. $CS = 'MODE' MCS 'NAVIER_STOKES' 'QUAF' ;
  766. VCS = 'DOMA' $CS 'XXDIAGSI' ;
  767. 'MESS' '==> Volume des condenseurs : '
  768. ('SOMT' VCS);
  769.  
  770. * Volume des brides
  771. BRIDE = MISTRA 'ELEM' 'VERT' ;
  772. 'TRACE' 'CACH' BRIDE 'TITR' 'Volume des brides' ;
  773. MBR = 'CHAN' BRIDE 'QUAF' ;
  774. $BR = 'MODE' MBR 'NAVIER_STOKES' 'QUAF' ;
  775. VBR = 'DOMA' $BR 'XXDIAGSI' ;
  776. 'MESS' '==> Masse des Brides : '
  777. (('SOMT' VBR) '*' 7800.0) ;
  778.  
  779. * Volume acier enceinte
  780. ACI = MISTRA 'ELEM' 'ROUG' ;
  781. 'TRACE' 'CACH' ACI 'TITR' 'Volume Acier Enceinte' ;
  782. MAC = 'CHAN' ACI 'QUAF' ;
  783. $AC = 'MODE' MAC 'NAVIER_STOKES' 'QUAF' ;
  784. VAC = 'DOMA' $AC 'XXDIAGSI' ;
  785. 'MESS' '==> Masse des Aciers : '
  786. (('SOMT' VAC) '*' 7800.0) ;
  787. 'FINSI' ;
  788.  
  789. *************************************************************
  790. * On s"occupe maintenant du maillage 2D Complet
  791. *
  792. *************************************************************
  793. SC1 = SC1_1 'ET' SC1_2 'ET' SC1_3 'ET' SC1_4 'ET'
  794. SC1_5 ;
  795. SC2 = SC2_1 'ET' SC2_2 'ET' SC2_3 'ET' SC2_4 ;
  796. SC3 = SC3_1 'ET' SC3_2 'ET' SC3_3 'ET' SC3_4 ;
  797. *************************************************************
  798. ** FIN DU MAILLAGE **
  799. *************************************************************
  800. VFF = MISTRA 'ELEM' 'BLAN' ;
  801. PAR = (MISTRA 'ELEM' 'ROUG') 'ET' (MISTRA 'ELEM' 'VERT') ;
  802.  
  803. * Recherche de la paroi externe
  804.  
  805. CPAR = 'CONT' PAR ;
  806. CVTF = 'CONT' VFF ;
  807. CINT = 'INTE' CPAR CVTF ;
  808. CEXT = 'DIFF' CPAR CINT ;
  809. NPT = 'NBEL' CEXT ;
  810. B1 = CEXT 'ELEM' ('LECT' 1 'PAS' 1 NM1) ;
  811. B2 = CEXT 'ELEM' ('LECT' (NPT '-' NM1 '+' 1) 'PAS' 1 NPT) ;
  812. CEXT = 'DIFF' CEXT B1 ;
  813. CEXT = 'DIFF' CEXT B2 ;
  814.  
  815.  
  816. * Mise en place de la procedure
  817.  
  818. RXT = 'TABLE' ;
  819. RXT.'vtf' = VFF ;
  820. RXT.'breche' = AB ;
  821. RXT.'axe' = AXE ;
  822. RXT.'diru1' = (0 1) ;
  823. RXT.'pi' = ( VFF 'POIN' 'PROC' (0.0 2.0)) ;
  824. RXT.'IMPR' = 2 ;
  825.  
  826. *RXT.'THERMP' = VRAI ;
  827. *RXT.'vtp' = PAR ;
  828. *RXT.'LAMBDA' = 17.0 ;
  829. *RXT.'ROCP' = (7800.0 '*' 500.0) ;
  830. *RXT.'Tp0' = 20.0 ;
  831. *RXT.'ECHAN' = 5.0 ;
  832.  
  833. *RXT.'ECHEXT' = VRAI ;
  834. *RXT.'parext' = CEXT ;
  835. *RXT.'HEXT' = 4.0 ;
  836. *RXT.'TPEXT' = 20.0 ;
  837.  
  838.  
  839. *RXT.'TIMP1' = 'TABLE' ;
  840. *RXT.'TIMP1'.'MAILLAGE' = SC1 ;
  841. *RXT.'TIMP1'.'t' = 'PROG' 0.0 25.5 345.5 100000.0 ;
  842. *RXT.'TIMP1'.'TIMP' = 'PROG' 135.0 135.0 120.0 120.0 ;
  843. *RXT.'TIMP1'.'ECHAN' = 30.0 ;
  844. *RXT.'TIMP2' = 'TABLE' ;
  845. *RXT.'TIMP2'.'MAILLAGE' = SC2 ;
  846. *RXT.'TIMP2'.'t' = 'PROG' 0.0 25.5 345.5 100000.0 ;
  847. *RXT.'TIMP2'.'TIMP' = 'PROG' 135.0 135.0 120.0 120.0 ;
  848. *RXT.'TIMP2'.'ECHAN' = 30.0 ;
  849. *RXT.'TIMP3' = 'TABLE' ;
  850. *RXT.'TIMP3'.'MAILLAGE' = SC3 ;
  851. *RXT.'TIMP3'.'t' = 'PROG' 0.0 25.5 345.5 100000.0 ;
  852. *RXT.'TIMP3'.'TIMP' = 'PROG' 135.0 135.0 120.0 120.0;
  853. *RXT.'TIMP3'.'ECHAN' = 30.0 ;
  854.  
  855. RXT.'VAPEUR' = FAUX ;
  856. RXT.'TF0' = 20.0 ;
  857. RXT.'PT0' = 1.E5 ;
  858. * RXT.'Yvap0' = 0.5 ;
  859.  
  860. *RXT.'MODTURB' = 'NUTURB' ;
  861. *RXT.'NUT' = 5.E-3 ;
  862. RXT.'MODTURB' = 'LMEL' ;
  863. RXT.'LMEL' = 0.08 ;
  864.  
  865. RXT.'epsi' = 1.E-5 ;
  866. RXT.'DISCR' = 'LINE' ;
  867. RXT.'KPRE' = 'MSOMMET' ;
  868. RXT.'GRAPH' = GRAPH ;
  869.  
  870. RXT.'scenario' = 'TABLE' ;
  871. RXT.'scenario'.'t' = 'PROG' 0.0 10000.0 ;
  872. RXT.'scenario'.'qair' = 'PROG' 0.1302 0.1302 ;
  873. * RXT.'scenario'.'qeau' = 'PROG' 0.150 0.150 0.0 0.0 0.12 0.12 ;
  874. RXT.'scenario'.'tinj' = 'PROG' 20.0 20.0 ;
  875.  
  876. *-- Definition des capteurs de mesures dans l"enceinte
  877. R00 = 0. 1.575 ;
  878. R01 = 0. 2.025 ;
  879. R02 = 0. 2.625 ;
  880. R03 = 0. 3.625 ;
  881. R04 = 0. 4.625 ;
  882. R05 = 0. 5.625 ;
  883. R06 = 0. 6.425 ;
  884. R07 = 0. 7.125 ;
  885. R10 = 0.475 1.575 ;
  886. R11 = 0.475 2.025 ;
  887. R12 = 0.475 2.625 ;
  888. R13 = 0.475 3.625 ;
  889. R14 = 0.475 4.625 ;
  890. R15 = 0.475 5.625 ;
  891. R16 = 0.475 6.425 ;
  892. R17 = 0.475 7.125 ;
  893. R20 = 0.950 1.575 ;
  894. R21 = 0.950 2.025 ;
  895. R22 = 0.950 2.625 ;
  896. R23 = 0.950 3.625 ;
  897. R24 = 0.950 4.625 ;
  898. R25 = 0.950 5.625 ;
  899. R26 = 0.950 6.425 ;
  900. R27 = 0.950 7.125 ;
  901. R28 = 0.950 0.600 ;
  902. R29 = 0.950 1.100 ;
  903. R30 = 1.425 1.575 ;
  904. R31 = 1.425 2.025 ;
  905. R32 = 1.425 2.625 ;
  906. R33 = 1.425 3.625 ;
  907. R34 = 1.425 4.625 ;
  908. R35 = 1.425 5.625 ;
  909. R36 = 1.425 6.425 ;
  910. R37 = 1.425 7.125 ;
  911. R40 = 1.814 1.575 ;
  912. R41 = 1.814 2.025 ;
  913. R42 = 1.814 2.625 ;
  914. R43 = 1.814 3.625 ;
  915. R44 = 1.814 4.625 ;
  916. R45 = 1.814 5.625 ;
  917. R46 = 1.814 6.425 ;
  918. R47 = 1.814 7.125 ;
  919. R51 = 2.000 2.625 ;
  920. R52 = 2.000 4.625 ;
  921. R53 = 2.000 6.425 ;
  922.  
  923. PERSO=TABLE;
  924. PERSO.CAPT0 = VRAI ;
  925. PERSO.TITRE = 'Tableperso' ;
  926. PERSO.TEMPS = 'PROG' ;
  927. PERSO.FREQ0 = 10 ;
  928. PERSO.LPOINTS= R00 et R01 et R02 et R03 et R04 et R05
  929. et R06 et R07 et R10 et R11 et R12 et R13 et R14 et
  930. R15 et R16 et R17 et R20 et R21 et R22 et R23 et R24
  931. et R25 et R26 et R27 et R28 et R29 et R30 et R31 et R32
  932. et R33 et R34 et R35 et R36 et R37 et R40 et R41 et
  933. R42 et R43 et R44 et R45 et R46 et R47 et R51 et R52
  934. et R53 ;
  935. PERSO.NPOINTS=45 ;
  936. PERSO.CAPT='MOTS' C00 C01 C02 C03 C04 C05 C06 C07
  937. C10 C11 C12 C13 C14 C15 C16 C17 C20 C21 C22 C23 C24
  938. C25 C26 C27 C28 C29
  939. C30 C31 C32 C33 C34 C35 C36 C37 C40 C41
  940. C42 C43 C44 C45 C46 C47 C51 C52 C53 ;
  941. PERSO.COMP= MOTS 'TF' ;
  942. * Trace de champs 3D
  943. PERSO.CH3D0 = FAUX ;
  944. PERSO.COMP3D = MOTS 'TF' ;
  945. PERSO.TPS3D = PROG ;
  946. PERSO.FREQ3D = 4 ;
  947.  
  948. *rxt.'PERSO' = VRAI ;
  949. *rxt.'PRCPERSO'= 'VISUP' ;
  950. *rxt.'TABPERSO'= PERSO ;
  951.  
  952. *************************************
  953. ** EXECUTION **
  954. *************************************
  955. RXT.'DT0' = 0.05 ;
  956. NSAV = 0 ;
  957. NPAS = 10 ;
  958. EXECRXT NPAS RXT ;
  959. NSAV = NSAV '+' NPAS ;
  960.  
  961. RXT.'DT0' = 0.5 ;
  962. NPAS = 11 ;
  963. EXECRXT NPAS RXT ;
  964. NSAV = NSAV '+' NPAS ;
  965.  
  966. RXT.'DT0' = 2.0 ;
  967. NPAS = NPTS ;
  968. EXECRXT NPAS RXT ;
  969. NSAV = NSAV '+' NPAS ;
  970.  
  971. *'OPTI' 'SAUVE' 'gonflairiso.save' ;
  972. *'SAUVE' ;
  973.  
  974. Si (COMPLET) ;
  975.  
  976. NPTS = 'DIME' RXT.'TIC'.'LTPS' ;
  977. PT = 'EXTR' RXT.'TIC'.'PT' NPTS ;
  978. TF = 'EXTR' RXT.'TIC'.'Tfm' NPTS ;
  979. PSOL = 3.90785E+05 ;
  980. TSOL = 99.64 ;
  981. erp=abs( PT - PSOL )/PSOL ;
  982. ert=abs( TF - TSOL )/TSOL;
  983. mess ' Erreur Pression = ' erp
  984. ' Erreur Temperature = ' ert ;
  985. Si (erp > 0.002); erreur 5 ; Finsi ;
  986. Si (ert > 0.002); erreur 5 ; Finsi ;
  987.  
  988. Finsi ;
  989.  
  990.  
  991. Si ('NON' COMPLET) ;
  992.  
  993. NPTS = 'DIME' RXT.'TIC'.'LTPS' ;
  994. PT = 'EXTR' RXT.'TIC'.'PT' NPTS ;
  995. TF = 'EXTR' RXT.'TIC'.'Tfm' NPTS ;
  996. PSOL = 1.04200E+05 ;
  997. TSOL = 23.432 ;
  998. erp=abs( PT - PSOL )/PSOL ;
  999. ert=abs( TF - TSOL )/TSOL;
  1000. mess ' Erreur Pression = ' erp
  1001. ' Erreur Temperature = ' ert ;
  1002. Si (erp > 0.002); erreur 5 ; Finsi ;
  1003. Si (ert > 0.002); erreur 5 ; Finsi ;
  1004.  
  1005. Finsi ;
  1006.  
  1007. 'FIN' ;
  1008. *
  1009.  
  1010.  
  1011.  
  1012.  
  1013.  

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