Télécharger gonfl2Dex.dgibi

Retour à la liste

Numérotation des lignes :

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

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