Télécharger hy2.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : hy2.dgibi
  2. **
  3. ** --- 2 JUIN 1993 ---
  4. ** PERTES DE CHARGE et GMV (Groupe Moto Ventilateur)
  5. **
  6. ** CANAL LONGUEUR 10. LARGEUR 1.
  7. **
  8. ** test cas isotherme NS FROT ET GMV
  9. **
  10. ** On considŠre l'‚coulement dans un canal plan vertical
  11. ** compris entre les plans x=0 et x=1. la hauteur est de 10.
  12. ** A proximite de la sortie du canal on a place un faisceau de tube.
  13. ** modelise par une perte de charge du type K U**BETA
  14. ** Dans la zone d'entree on a place un GMV
  15. **
  16. **
  17. **
  18. **
  19. **
  20. **
  21. OPTION DIME 2 ELEM QUA4 ;
  22. GRAPH=VRAI ;
  23. GRAPH=FAUX ;
  24. FIN;
  25. p1=0 0.;
  26. p2=1. 0.;
  27. entree= p1 d 14 p2 ;
  28. c1= 5. 0. ;
  29.  
  30. * ikas= 0 --> DROIT (option par defaut) ikas=1 --> COURBE
  31. ikas=0 ;
  32. mess 'ikas= 0 --> DROIT (option par defaut) ikas=1 --> COURBE ? ';
  33. *obtenir ikas*entier ;
  34.  
  35. si (EGA ikas 1) ;
  36.  
  37. q1=p1 tour c1 -90. ;
  38. q2=p2 tour c1 -90. ;
  39. pp1=p1 c c1 q1 20 ;
  40. pp2=p2 c c1 q2 20 ;
  41. sortie=entree tour c1 -90 ;
  42.  
  43. sinon ;
  44.  
  45. q1=p1 plus (0 10) ;
  46. q2=p2 plus (0 10) ;
  47. pp1=p1 d q1 20 ;
  48. pp2=p2 d q2 20 ;
  49. sortie=entree plus (0 10) ;
  50. finsi ;
  51.  
  52. pp1=inve pp1 ;
  53. sortie=inve sortie;
  54. elim 0.0001 (sortie et pp1 et pp2 );
  55. cnt=entree et pp2 et sortie et pp1 ;
  56. *bell=surf cnt ;
  57.  
  58. bell=daller entree pp2 sortie pp1 ;
  59. angle=0. ;
  60. *obtenir angle*flottant ;
  61.  
  62.  
  63. bell=bell tour c1 angle ;
  64. entree=entree tour c1 angle ;
  65. sortie=sortie tour c1 angle ;
  66. pp1=pp1 tour c1 angle ;
  67. pp2=pp2 tour c1 angle ;
  68.  
  69. *pp1=inve pp1 ;
  70. elim (pp1 et pp2 et entree et sortie et bell) 0.0001 ;
  71. entref=entree elem seg2 (lect 2 pas 1 13) ;
  72. *entree=inve entree ;
  73. *sortie=inve sortie ;
  74.  
  75. $bell=doma bell ;
  76. $bell.titre= 'CHAINE' 'TEST PERTES DE CHARGE ET GMV' ;
  77. $entree = doma entree 'INCL' $bell 1.e-4 ;
  78. $entref = doma entref 'INCL' $bell 1.e-4 ;
  79. $sortie = doma sortie 'INCL' $bell 1.e-4 ;
  80. $pp1 = doma pp1 'INCL' $bell 1.e-4 ;
  81. $pp2 = doma pp2 'INCL' $bell 1.e-4 ;
  82. ech=bell elem qua4 (lect 211 pas 1 252) ;
  83. $ech = doma ech 'INCL' $bell 1.e-4 ;
  84.  
  85. gmvi=bell elem qua4 ( (lect 62 pas 1 65) et (lect 76 pas 1 79)) ;
  86. cgmv=cont gmvi ;
  87. legmv=cgmv elem seg2 (lect 1 pas 1 4) ;
  88. lsgmv=cgmv elem seg2 (lect 7 pas 1 10) ;
  89. $gmv=doma gmvi 'INCL' $bell 1.e-5 ;
  90. egmv=gmvi elem 'APPUYE' 'LARGEMENT' legmv;
  91. sgmv=gmvi elem 'APPUYE' 'LARGEMENT' lsgmv;
  92. $egmv=doma egmv 'INCL' $gmv 1.e-5 ;
  93. $sgmv=doma sgmv 'INCL' $gmv 1.e-5 ;
  94.  
  95.  
  96. CK= 100. 100. ;
  97. CB= 2. 2. ;
  98.  
  99. tabgmv=table ;
  100. tabgmv.'DIR'= 0. 1. ;
  101. tabgmv.'PENTREE'=($egmv.centre);
  102. tabgmv.'PSORTIE'=($sgmv.centre);
  103. tabgmv.'LDEBIT'=legmv ;
  104. tabgmv.'IMPR'=5 ;
  105. *tabgmv.'KIMP'=10. ;
  106. tabgmv.'GMV'= EVOL 'MANU'
  107. 'DEBIT'
  108. (prog 0.3 0.6 0.8 1.)
  109. 'PRESSION'
  110. (prog 180. 130. 80. 0.) ;
  111. Si GRAPH ;dessin tabgmv.'GMV' ;FINSI ;
  112. tabgmv.omega=0.2 ;
  113. list tabgmv ;
  114.  
  115. nu=5.E-1 ;
  116. tpsc=1500 ;
  117. zero=1.e-20 ;
  118.  
  119. rv=eqex $bell 'DUMP' ALFA 0.7 ITMA 400
  120. ZONE $BELL OPER NS NU INCO 'UN'
  121. ZONE $ech OPER 'FROT' CK CB INCO 'UN'
  122. ZONE $gmv OPER 'GMV' tabgmv INCO 'UN'
  123. CLIM 'UN' UIMP GMVi 0. ;
  124.  
  125. rvp= eqpr $bell zone $bell oper PRESSION 0.
  126. zone $PP1 oper VNIMP 0.
  127. zone $PP1 oper VTIMP 0.
  128. zone $PP2 oper VNIMP 0.
  129. zone $PP2 oper VTIMP 0. ;
  130.  
  131. rv.'INCO'.'UN' = kcht $bell vect sommet (1.e-5 1. ) ;
  132. rv.pression=rvp ;
  133. rv.'FIDT'=20;
  134.  
  135. lh= (noeu 10) et (noeu 20) et (noeu 30) et (noeu 40) et
  136. (noeu 50) et (noeu 60) ;
  137. lj= (manu poi1 ($bell.maillage poin proc( 0.5 0.5) ) ) ;
  138. his = khis 'UN' 1 lh 'UN' 2 (lh et lj) ;
  139. rv.'HIST'=his ;
  140. opti veri 1;
  141. exec rv ;
  142.  
  143. PN = rvp.'PN' ;
  144. un= (rv.'INCO'.'UN') ;
  145. ung1 = vect 0.1 un ux uy jaune ;
  146.  
  147. qe=dbit un $entree ;
  148. qs=dbit un $sortie ;
  149. dq=(abs qe )-(abs qs) ;
  150. mess (' BILAN : dq=') dq ;
  151. Si GRAPH ;
  152. dessin his.'TABD' his.'1UN' ;
  153. dessin his.'TABD' his.'2UN' ;
  154.  
  155. trace ung1 bell ;
  156. trace pn bell ;
  157. FINSI ;
  158. FIN ;
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  

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