Télécharger inj.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : inj.dgibi
  2. *
  3. * INJ.DGIBI (Cataldo Caroli) Mai 99
  4. * Teste la procedure ENCEINTE
  5. *
  6. *
  7. *
  8. *
  9. fin ;
  10. *opti trace psc ;
  11. COMPLET= FAUX ;
  12. debproc visu rx*table ;
  13.  
  14. epsi = 1.0e-3 ;
  15. px = -1000.0 0.0 0.0 ;
  16. py = 0.0 -1000.0 0.0 ;
  17. pz = 0.0 0.0 -1000.0 ;
  18. px1 = -1000.0 0.0 10.0 ;
  19. py1 = 0.0 -1000.0 10.0 ;
  20. pz1 = 0.0 0.0 10.0 ;
  21.  
  22. $vtf = rx.'TIC'.'$vtf' ;
  23. $vtp = rx.'TIC'.'$vtp' ;
  24. vtf = doma $vtf maillage ;
  25. vtp = doma $vtp maillage ;
  26.  
  27. un = rx.tic.'UN' ;
  28. tf = rx.tic.'TF' ;
  29. rvp = rx.tic.'RVP' ;
  30. rho = rx.tic.'RHO' ;
  31. tp = rx.tic.'TP' ;
  32. si (rxt.'ASPERS') ;
  33. vn = rx.tic.'VN' ;
  34. xd = rx.tic.'XD' ;
  35. dd = rx.tic.'DD' ;
  36. td = rx.tic.'TD' ;
  37. kgas = rx.tic.'KGAS' ;
  38. kdis = rx.tic.'KDIS' ;
  39. hgas = rx.tic.'HGAS' ;
  40. hdis = rx.tic.'HDIS' ;
  41. sdia = rx.tic.'SDIA' ;
  42. sxd = rx.tic.'SXD' ;
  43. srho = rx.tic.'SRHO' ;
  44. finsi ;
  45. roi = rx.tic.'ROI' ;
  46. rovi = rx.tic.'ROVI' ;
  47.  
  48. x y z = coor vtf ;
  49. ptx = x poin egsupe 0.0 ;
  50. pty = y poin egsupe 0.0 ;
  51. ptz = z poin inferieur 11.0 ;
  52. vtfx = vtf elem appuie strictement ptx ;
  53. vtfy = vtf elem appuie strictement pty ;
  54. vtfz = vtf elem appuie strictement ptz ;
  55.  
  56. x y z = coor vtp ;
  57. ptx = x poin egsupe 0.0 ;
  58. pty = y poin egsupe 0.0 ;
  59. ptz = z poin eginfe 11.0 ;
  60. vtpx = vtp elem appuie strictement ptx ;
  61. vtpy = vtp elem appuie strictement pty ;
  62. vtpz = vtp elem appuie strictement ptz ;
  63.  
  64. ptxz = poin (chan vtf poi1) plan p0 px pz epsi ;
  65. ptyz = poin (chan vtf poi1) plan p0 py pz epsi ;
  66. ptxy = poin (chan vtf poi1) plan px1 py1 pz1 1e-0 ;
  67. plxz = elem (chan vtf poi1) 'APPUIE' 'STRICTEMENT' ptxz ;
  68. plyz = elem (chan vtf poi1) 'APPUIE' 'STRICTEMENT' ptyz ;
  69. plxy = elem (chan vtf poi1) 'APPUIE' 'STRICTEMENT' ptxy ;
  70.  
  71. un_v = vect un 1.0 ux uy uz jaun ;
  72. uxzv = vect (redu un ptxz ) 1.0 ux uy uz jaun ;
  73. uyzv = vect (redu un ptyz ) 1.0 ux uy uz jaun ;
  74. uxyv = vect (redu un ptxy ) 1.0 ux uy uz jaun ;
  75.  
  76. *titr 'Un ' ; trac un_v vtf ;
  77. *titr 'Un yz' ; trac px uyzv plyz ;
  78. titr 'Un xz' ; trac py uxzv plxz ;
  79. titr 'Un xy' ; trac pz uxyv plxy ;
  80.  
  81. si (ega (rxt.'discr') 'QUAF') ;
  82. vtfx = chan vtfx 'DECL' ;
  83. vtfy = chan vtfy 'DECL' ;
  84. vtfz = chan vtfz 'DECL' ;
  85. vtpx = chan vtpx 'DECL' ;
  86. vtpy = chan vtpy 'DECL' ;
  87. vtpz = chan vtpz 'DECL' ;
  88. finsi ;
  89.  
  90. mess 'min Tf ' (mini tf) 'max Tf ' (maxi tf) ;
  91. * titr 'Tf yz' ; trac px vtfx tf ;
  92. titr 'Tf xz' ; trac py vtfy tf ;
  93. * titr 'Tf xy' ; trac pz vtfz tf ;
  94. mess 'min Tp ' (mini tp) 'max Tp ' (maxi tp) ;
  95. * titr 'Tp yz' ; trac px vtpx tp ;
  96. titr 'Tp xz' ; trac py vtpy tp ;
  97. * titr 'Tp xy' ; trac pz vtpz tp ;
  98. mess 'min Rho ' (mini rho) 'max Rho ' (maxi rho) ;
  99. * titr 'Rho yz' ; trac px vtfx rho ;
  100. titr 'Rho xz' ; trac py vtfy rho ;
  101. * titr 'Rho xy' ; trac pz vtfz rho ;
  102. mess 'min Rvp ' (mini rvp) 'max Rvp ' (maxi rvp) ;
  103. * titr 'Rvp yz' ; trac px vtfx rvp ;
  104. titr 'Rvp xz' ; trac py vtfy rvp ;
  105. * titr 'Rvp xy' ; trac pz vtfz rvp ;
  106.  
  107. si (rxt.'ASPERS') ;
  108. vn = rx.tic.'VN' ;
  109. vxzv = vect (redu vn plxz ) 1.0e-0 ux uy uz roug ;
  110. vyzv = vect (redu vn plyz ) 1.0e-0 ux uy uz roug ;
  111. vxyv = vect (redu vn plxy ) 1.0e-0 ux uy uz roug ;
  112. cmps = extr un comp ;
  113. us = (psca un un cmps cmps) ** 0.5 ;
  114. vs = (psca vn vn cmps cmps) ** 0.5 ;
  115. * titr 'Vn yz' ; trac px vyzv plyz ;
  116. titr 'Vn xz' ; trac py vxzv plxz ;
  117. * titr 'Vn xy' ; trac pz vxyv plxy ;
  118. mess 'min Dd ' (mini dd) 'max Dd ' (maxi dd) ;
  119. * titr 'Dd yz' ; trac px vtfx dd ;
  120. titr 'Dd xz' ; trac py vtfy dd ;
  121. * titr 'Dd xy' ; trac pz vtfz dd ;
  122. mess 'min Xd ' (mini xd) 'max Xd ' (maxi xd) ;
  123. * titr 'Xd yz' ; trac px vtfx xd ;
  124. titr 'Xd xz' ; trac py vtfy xd ;
  125. * titr 'Xd xy' ; trac pz vtfz xd ;
  126. mess 'min Td ' (mini td) 'max Td ' (maxi td) ;
  127. * titr 'Td yz' ; trac px vtfx td ;
  128. titr 'Td xz' ; trac py vtfy td ;
  129. * titr 'Td xy' ; trac pz vtfz td ;
  130. finsi ;
  131. finproc ;
  132.  
  133.  
  134.  
  135.  
  136.  
  137. opti dime 3 elem cu20 ;
  138.  
  139. ri = 10.000 ; sp = 1.000 ; hc = 20.000 ;
  140. 'SI' COMPLET ;
  141. n1 = 1 ; n2 = 2 ; n3 = 1 ;
  142. n4 = 4 ; ;
  143. 'SINON' ;
  144. n1 = 1 ; n2 = 1 ; n3 = 1 ;
  145. n4 = 1 ; ;
  146. 'FINSI' ;
  147. epsi = 1.000e-2 ; ;
  148.  
  149. p0 = 0.000 0.000 0.000 ;
  150. px = -1000.000 0.000 0.000 ;
  151. py = 0.000 -1000.000 0.000 ;
  152. pz = 0.000 0.000 1000.000 ;
  153. cd = 0.000 0.000 -20.000 ;
  154. ph0 = 0.000 0.000 hc ;
  155. phx = ri 0.000 hc ;
  156. phy = 0.000 ri hc ;
  157.  
  158. p1 = (ri*0.4) 0.000 0.000 ;
  159. p2 = (ri*0.3) (ri*0.3) 0.000 ;
  160. p3 = 0.000 (ri*0.4) 0.000 ;
  161. p4 = ri 0.000 0.000 ;
  162. p5 = 0.000 ri 0.000 ;
  163. p6 = (ri+sp) 0.000 0.000 ;
  164. p7 = 0.000 (ri+sp) 0.000 ;
  165.  
  166. l1 = droi p0 p1 n1 ;
  167. l2 = droi p1 p2 n1 ;
  168. l3 = droi p2 p3 n1 ;
  169. l4 = droi p3 p0 n1 ;
  170. l5 = cerc p4 p0 p5 (2*n1) ;
  171. l6 = cerc p6 p0 p7 (2*n1) ;
  172.  
  173. basf = dall l1 l2 l3 l4 ;
  174. basf = basf et (regl (l2 et l3) l5 n2) ; elim basf epsi ;
  175. basf = basf et (syme basf droi p0 p3) ; elim basf epsi ;
  176. basf = basf et (syme basf droi p0 p1) ; elim basf epsi ;
  177.  
  178. basm = regl l5 l6 n3 ;
  179. basm = basm et (syme basm droi p0 p3) ; elim basm epsi ;
  180. basm = basm et (syme basm droi p0 p1) ; elim basm epsi ;
  181.  
  182. domf = basf proj cyli pz sphe cd phx ;
  183. domm = basm proj cyli pz plan ph0 phx phy ;
  184.  
  185. aria = volu basf domf n4 ;
  186. muro = (volu basm domm n4) coul roug ;
  187. elim (aria et muro) epsi ;
  188.  
  189. pjg = poin basf proc (0.000 0.000 0.000) ;
  190. *pjg = poin basf proc ((ri/2.0) 0.000 0.000) ;
  191. jg = (elem basf 'APPUIE' 'LARGEMENT' pjg) coul vert ;
  192.  
  193. pjd = poin domf proc (0.000 0.000 hc) ;
  194. jd = (elem domf 'APPUIE' 'LARGEMENT' pjd) coul vert ;
  195.  
  196.  
  197. DT0 = 1.0e-0 ;
  198.  
  199. rxt = table ;
  200.  
  201. rxt.'vtf' = aria ;
  202. rxt.'vtp' = muro ;
  203. rxt.'breche' = jg ;
  204. *rxt.'aspinj' = jd ;
  205. *rxt.'toitf' = domf ;
  206. rxt.'pi' = (0.0 0.0 10.0) ;
  207. *
  208. *rxt.'vtf' = vtf ;
  209. *rxt.'vtp' = vtp ;
  210. *rxt.'breche' = breche ;
  211. *rxt.'paroif' = paroif ;
  212. *rxt.'aspinj' = breche ;
  213. *rxt.'pi' = (-3.0 -3.0 50.0) ; ;
  214.  
  215. *rxt.'ASPERS' = VRAI ;
  216. rxt.'ASPERS' = FAUX ;
  217. *rxt.'rod' = 1.000e+3 ;
  218. *rxt.'Cpd' = 4.180e+3 ;
  219.  
  220. rxt.'scenario' = table ;
  221. rxt.'scenario'.'TF0' = 50.0 ;
  222. rxt.'scenario'.'PT0' = 1.0e5 ;
  223. rxt.'scenario'.'Yvap0' = 0.1 ;
  224. rxt.'scenario'.'Yhe0' = 0.0 ;
  225. rxt.'scenario'.'TP0' = 20.0 ;
  226. rxt.'scenario'.'t' = prog 0.0 100.0 ;
  227. rxt.'scenario'.'Qvap' = prog 10.0 10.0 ;
  228. rxt.'scenario'.'qair' = prog 0.0 0.0 ;
  229. rxt.'scenario'.'qhe' = prog 0.0 0.0 ;
  230. rxt.'scenario'.'tinj' = prog 100.0 100.0 ;
  231. *
  232. *rxt.'scenario'.'vzinj' = prog -1.000e-1 -1.000e-1 ;
  233. *rxt.'scenario'.'xdinj' = prog 1.000e-3 1.000e-3 ;
  234. *rxt.'scenario'.'tdinj' = prog 5.000e+1 5.000e+1 ;
  235. *rxt.'scenario'.'ddinj' = prog 5.000e-4 5.000e-4 ;
  236.  
  237. *rxt.'discr' = 'QUAF' ;
  238. *rxt.'ALGO' = 'IMPL' ;
  239. *rxt.'KPRESS' = 'CENTREP1' ;
  240.  
  241. rxt.'epsi' = 1.0e-4 ;
  242. rxt.'DT0' = DT0 ;
  243.  
  244. rxt.'VAPEUR' = VRAI ;
  245. rxt.'lamb' = 2.300 ;
  246. rxt.'rocpb' = 2.100e+6 ;
  247. rxt.'muair' = 1.800e-5 ;
  248. rxt.'alf' = 1.800e-5 ;
  249. rxt.'lamair' = 2.580e-2 ;
  250. rxt.'db' = 1.000e-5 ;
  251. rxt.'KH' = 1.000e+1 ;
  252. rxt.'MUT' = 1.000e-2 ;
  253.  
  254. rxt.'Cpair' = 1000.0 ;
  255. rxt.'Cph2o' = 1700.0 ;
  256. rxt.'Cphe' = 14550.0 ;
  257. rxt.'Lv' = 2.3E6 ;
  258. rxt.'Rgair' = ((80. * 296.9) + (20. * 259.8))/100. ;
  259. rxt.'Rgh2o' = 461.0 ;
  260. rxt.'Rghe' = 4157.0 ;
  261.  
  262. rxt.'FROG'=1. ;
  263. rxt.'FROC'=0.2 ;
  264.  
  265. rxt.'alfg'=0.5 ;
  266. rxt.'alfc'=0.5 ;
  267.  
  268. 'SI' COMPLET ;
  269. rxt.'GRAPH'= vrai ;
  270. enceinte 10 rxt ;
  271. visu rxt ;
  272. 'SINON' ;
  273. enceinte 3 rxt ;
  274. un = (rxt.tic.'UN') ;
  275. vmax=maxi un ;
  276. vmin=mini un ;
  277. mess ' Vmax = ' vmax 'Vmin = ' vmin ;
  278. er = (abs (vmax - .70325) ) + (abs (vmin + 4.05197E-02)) ;
  279. mess ' Ecart sur Min Max ' er ;
  280. err1= 5.e-5 ;
  281. si ( er > err1 ) ; erreur 5 ; finsi ;
  282.  
  283. 'FINSI' ;
  284. 'FIN' ;
  285.  
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  
  294.  
  295.  
  296.  
  297.  
  298.  

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