Télécharger la plaquette Cast3M


fr en

Télécharger filtreke.procedur

Retour à la liste

Numérotation des lignes : 

  1. * FILTREKE  PROCEDUR  MAGN      15/07/27    21:15:04     8592           
  2. *
  3. * VERSION DU 16-12-97
  4. *
  5.  
  6.  DEBP FILTREKE ;
  7.  ARGU RX*TABLE ;
  8. * Filtre sur K et Epsilon
  9. * - Echelle de vitesse (K**0.5) inférieure à une fraction (alfk)
  10. *   de Uref (vitesse caractéristique) (alfk=1 pour l'instant)
  11. *   Uref=max(UN,U0) 
  12. *  - K > K0 
  13. *  - Epsilon tel que l'echelle de longueur reste inférieure
  14. *   à (L0/a) où L0 = diamètre enceinte et a=f(Re)
  15. *   => Nut < Uref*L0/a
  16.  
  17.  rv=rx.'EQEX' ;
  18.  iarg=rx.'IARG' ;
  19.  NASTOK = rv.'NAVISTOK' ;
  20.  
  21.  'SI' (('EGA' NASTOK 0) ou ('EGA' IARG 2)) ;
  22.  
  23. * Ancien Filtre sur K et Epsilon
  24. *- Echelle de vitesse (K**0.5) inférieure à une fraction (alfk)
  25. * de U0 (vitesse caractéristique) (alfk=1 pour l'instant)
  26. * K > 0  Epsilon tel que l'echelle de longueur reste inférieure
  27. * à L0 diamètre enceinte
  28. *   => Nut < U0*L0
  29.  
  30.  si( non ( ega iarg 2)) ;
  31.  mess 'Procedure FILTREKE : nombre d arguments incorrect ' iarg ;
  32.  quitter FILTREKE  ;
  33.  finsi ;
  34.  si ( ega ('TYPE' rx.'ARG1') 'MOT     ') ;
  35.  U1=rv.'INCO'.(rx.'ARG1') ;
  36.  sinon ;
  37.  si ( ega ('TYPE' (rx.'ARG1')) 'FLOTTANT') ;
  38.  U1=rx.'ARG1' ;
  39.  sinon ;
  40.  mess 'Procedure FILTREKE : type argument 1 invalide ' ;
  41.  quitter FILTREKE ;
  42.  finsi ;
  43.  finsi ;
  44.  
  45.  si ( ega ('TYPE' rx.'ARG2') 'MOT     ') ;
  46.  L0=rv.'INCO'.(rx.'ARG2') ;
  47.  sinon ;
  48.  si ( ega ('TYPE' (rx.'ARG2')) 'FLOTTANT') ;
  49.  L0=rx.'ARG2' ;
  50.  sinon ;
  51.  mess 'Procedure FILTREKE : type argument 2 invalide ' ;
  52.  quitter FILTREKE ;
  53.  finsi ;
  54.  finsi ;
  55.  
  56.  
  57.  nic=dime (rx.'LISTINCO') ;
  58.  si( non ( ega nic 2)) ;
  59.  mess 'Procedure FILTREKE : nombre d inconnues incorrect ' nic ;
  60.  quitter FILTREKE ;
  61.  finsi ;
  62.  nomi1=extr 1 (rx.'LISTINCO');
  63.  nomi2=extr 2 (rx.'LISTINCO');
  64.  nom1= mot (text (chai nomi1)) ;
  65.  nom2= mot (text (chai nomi2)) ;
  66.  
  67.  en=rv.'INCO'.nom2 ;
  68.  kn=rv.'INCO'.nom1 ;
  69. *Rec=100.;
  70.  k0 = 1.e-10 ;
  71.  
  72. *lcu1=extr u1 'COMP'         ;
  73. *mdu2=u1 lcu1 'PSCA' u1 lcu1 ;
  74. *mdu=mdu2 ** 0.5 ;
  75. *Re=kops (kops (kops mdu '*' L0) '/' nu) '+' (Rec / 10.) ;
  76. *a= exp (kops Rec '/' Re ) ;
  77.  
  78.  mdu2=u1*u1 ;
  79.  a=0.09 ;
  80.  
  81.  kn=kops kn '|<' k0 ;
  82.  kn=kops kn '>|' mdu2 ;
  83.  
  84.  E0= kops (kops kn '**' 1.5) '*' (a / L0) ;
  85.  en=kops en '|<' E0 ;
  86.  
  87.  rv.'INCO'.'EN'=en ;
  88.  rv.'INCO'.'KN'=kn ;
  89.  
  90.  si( non ( ega NASTOK 0 ) ) ;
  91.  as2 ama1  = 'KOPS' 'MATRIK' ;
  92.  RESPRO as2 ama1 ;
  93.  finsi ;
  94.  
  95. 'QUITTER' FILTREKE ;
  96. 'FINSI'  ;
  97.  
  98.  si( non ( ega iarg 4)) ;
  99.  mess 'Procedure FILTREKE : nombre d arguments incorrect ' iarg ;
  100.  quitter FILTREKE  ;
  101.  finsi ;
  102.  si ( ega ('TYPE' rx.'ARG1') 'MOT     ') ;
  103.  U1=rv.'INCO'.(rx.'ARG1') ;
  104.  sinon ;
  105.  si ( ega ('TYPE' (rx.'ARG1')) 'FLOTTANT') ;
  106.  U1=rx.'ARG1' ;
  107.  sinon ;
  108.  mess 'Procedure FILTREKE : type argument 1 invalide ' ;
  109.  quitter FILTREKE ;
  110.  finsi ;
  111.  finsi ;
  112.  
  113.  si ( ega ('TYPE' rx.'ARG2') 'MOT     ') ;
  114.  L0=rv.'INCO'.(rx.'ARG2') ;
  115.  sinon ;
  116.  si ( ega ('TYPE' (rx.'ARG2')) 'FLOTTANT') ;
  117.  L0=rx.'ARG2' ;
  118.  sinon ;
  119.  mess 'Procedure FILTREKE : type argument 2 invalide ' ;
  120.  quitter FILTREKE ;
  121.  finsi ;
  122.  finsi ;
  123.  
  124.  si ( ega ('TYPE' rx.'ARG3') 'MOT     ') ;
  125.  NU=rv.'INCO'.(rx.'ARG3') ;
  126.  sinon ;
  127.  si ( ega ('TYPE' (rx.'ARG3')) 'FLOTTANT') ;
  128.  NU=rx.'ARG3' ;
  129.  sinon ;
  130.  mess 'Procedure FILTREKE : type argument 3 invalide ' ;
  131.  quitter FILTREKE ;
  132.  finsi ;
  133.  finsi ;
  134.  
  135.  si ( ega ('TYPE' rx.'ARG4') 'MOT     ') ;
  136.  UN=rv.'INCO'.(rx.'ARG4') ;
  137.  sinon ;
  138.  si ( ega ('TYPE' (rx.'ARG4')) 'CHPOINT') ;
  139.  UN=rx.'ARG4' ;
  140.  sinon ;
  141.  mess 'Procedure FILTREKE : type argument 4 invalide ' ;
  142.  quitter FILTREKE ;
  143.  finsi ;
  144.  finsi ;
  145.  
  146.  
  147.  nic=dime (rx.'LISTINCO') ;
  148.  si( non ( ega nic 2)) ;
  149.  mess 'Procedure FILTREKE : nombre d inconnues incorrect ' nic ;
  150.  quitter FILTREKE ;
  151.  finsi ;
  152.  
  153.  nomi1=extr 1 (rx.'LISTINCO');
  154.  nomi2=extr 2 (rx.'LISTINCO');
  155.  nom1= mot (text (chai nomi1)) ;
  156.  nom2= mot (text (chai nomi2)) ;
  157.  
  158.  en=rv.'INCO'.nom2 ;
  159.  kn=rv.'INCO'.nom1 ;
  160.  
  161.  Rec=100.;
  162.  k0 = 1.e-10 ;
  163.  cnu=0.09;
  164.  
  165.  lcu=extr un 'COMP'       ;
  166.  mdu=un lcu 'PSCA' un lcu ;
  167.  mdu=mdu ** 0.5 ;
  168.  
  169.  Re=kops (kops (kops mdu '*' L0) '/' nu) '+' (Rec / 10.) ;
  170.  a= exp (kops Rec '/' Re ) ;
  171.  
  172.  mdu = kops mdu '|<' u1 ;
  173.  mdu2= kops mdu '*' mdu ;
  174.  
  175.  kn=kops kn '|<' k0 ;
  176.  kn=kops kn '>|' mdu2 ;
  177.  
  178.  E0= kops (kops kn '**' 1.5) '*' (a / L0) ;
  179.  en=kops en '|<' E0 ;
  180.  
  181.  rv.'INCO'.nom2=en ;
  182.  rv.'INCO'.nom1=kn ;
  183.  
  184.  as2 ama1  = 'KOPS' 'MATRIK' ;
  185.  RESPRO as2 ama1 ;
  186.  
  187.  FINPROC ;
  188.  
  189.  
  190.  

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