Télécharger la plaquette Cast3M


fr en

Télécharger eautemp.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  eautemp.dgibi
  2. * repertoire des fichiers "divers"
  3. DIVERS = VENV 'CASTEM_DIVERS';
  4. *
  5.  SAUT PAGE ;
  6. ********************************************************************
  7. *               Test de validation des operateurs CHI1 et CHI2
  8. *               eau + gradient thermique + pression CO2
  9. *********************************************************************
  10. OPTION DIME 2 ;                                                         
  11. A= 0. 0. ;                                                              
  12. B= 5. 0. ;
  13. OPTION  ELEM QUA4 ;                                                     
  14.  AB= A DROIT 1 B ;                                                      
  15.  VEC1= 0. 40. ;
  16.  VEC2= 0. 60. ;
  17.  M1 = AB TRANS DINI 10. DFIN 20. VEC1 ;
  18.  CD= INVE ( COTE 3 M1 );
  19.  M2= M1 TRANS DINI 30. DFIN 50. VEC2 ;
  20.  EF= COTE 3 M2 ;
  21.  TASS M2 ;
  22.  MT = CHANGE M2 POI1  ;
  23.   XX YY = COOR MT ;
  24. *
  25. *                table de données pour CHI1
  26. *
  27. TABDON=TABLE ;                                                          
  28. TABDON.BDD= 'STRASBG' ;
  29. TABDON.IDEN= LECT 1 50 60 61 101 165 ;
  30.  TABDON.CHXMX= LECT 2144 2148 2157 2166 2168 2192 2198 2200
  31.     2224 2231 2249 2278 2281 2372 ;
  32.  TABCLIM=TABLE ;
  33.  TABCLIM.TYP6= LECT 2002 ;
  34.  TABCLIM.TYP3= LECT 2400 ;
  35.  TABCLIM.COMP3= LECT 101 ;
  36.  TABDON.CLIM=TABCLIM ;
  37.  TABDON.TEMPERATURE = 'OUI' ;
  38. *
  39.    TB1=CHI1 TABDON    
  40.        LOGK ('CHAINE' DIVERS '/COMPSM')
  41.        ENTH  ('CHAINE' DIVERS '/COMPSM')
  42.        COMP  ('CHAINE' DIVERS '/COMPSM') ;
  43. *
  44. *
  45. *                   table de données pour CHI2
  46. *
  47.  TBPARM= TABLE ;
  48.  TBPARM.ITMAX = 95;
  49.  TBPARM.ITERSOLI  = 10 ;
  50.  TBPARM.NFI  = 8;
  51.  TBPARM.EPS= 1.D-4 ;
  52.  TBPARM.SORTIE= MOTS 'PREC' 'FION' 'TYP5' ;
  53.  TBPAR2= TABLE ;
  54.  TBPAR2.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  55.  TBPAR2.LOGC= MANU CHPO MT 6 X001 -5.6 X050 -7.0 X060 -5.0
  56.               X061 -3. X101 -7. X165 -5.  ;
  57.   TOT001= MANU CHPO MT 1 X001  6.D-4 ;
  58.   TOT050= MANU CHPO MT 1 X050  -1.15D-3  ;
  59.   TOT060= MANU CHPO MT 1 X060  0.D0 ;
  60.   TOT061= MANU CHPO MT 1 X061  2.8D-2 ;
  61.   TOT101= MANU CHPO MT 1 X101  0.D0 ;
  62.   TOT165= MANU CHPO MT 1 X165  5.5D-5 ;
  63.  TBPAR2.TOT= TOT001+ TOT050 + TOT060 + TOT061 + TOT101 + TOT165 ;
  64. *
  65. * définition du champ de temperature
  66. * gradient constant entre l'entrée (180°C) et la sortie (25°C)
  67. *
  68.   B= 180. ;
  69.   A= -1.55 ;
  70.   XX YY = COOR MT ;
  71.   TMPE= A*YY + B ;
  72.   TBPAR2.TEMPE= TMPE ;
  73. *
  74. *
  75. *  calcul de la contrainte KCO2 pour avoir une pression  partielle
  76. *  de CO2 = 10**(-2.32)
  77. *  KCO2 = log(k(T)) - log ( pco2)
  78. *
  79.   COO = 2.32 ;
  80. *
  81.   TSTA1 = 0.d0  ;
  82.   TSTA2 = 25.d0 ;
  83.   TSTA3 = 60.d0 ;
  84.   TSTA4 = 100.d0;
  85.   TSTA5 = 150.d0;
  86.   TSTA6 = 200.d0;
  87.   TSTA7 = 250.d0;
  88.   TSTA8 = 300.d0;
  89. *
  90.   K11 = -18.14  ;
  91.   K21 = -18.13  ;
  92.   K31 = -18.25  ;
  93.   K41 = -18.58  ;
  94.   K51 = -19.09  ;
  95.   K61 = -19.82  ;
  96.   K71 = -21.04  ;
  97.   K81 = -24.07  ;
  98.   K12 = 1.001232E-03    ;
  99.   K22 = -8.024372E-04   ;
  100.   K32 = -6.155775E-03   ;
  101.   K42 = -9.206214E-03   ;
  102.   K52 = -1.241482E-02   ;
  103.   K62 = -1.553448E-02   ;
  104.   K72 = -4.244730E-02   ;
  105.   K82 = 0.000000E+00    ;
  106.   K13 = 0.000000E+00    ;
  107.   K23 = -7.214677E-05   ;
  108.   K33 = -8.080574E-05   ;
  109.   K43 = 4.544751E-06    ;
  110.   K53 = -6.871693E-05   ;
  111.   K63 = 6.323791E-06    ;
  112.   K73 = -5.445802E-04   ;
  113.   K83 = 0.000000E+00    ;
  114.   K14 = -9.619569E-07   ;
  115.   K24 = -8.246646E-08   ;
  116.   K34 = 7.112541E-07    ;
  117.   K44 = -4.884112E-07   ;
  118.   K54 = 5.002715E-07    ;
  119.   K64 = -3.672694E-06   ;
  120.   K74 = 3.630535E-06    ;
  121.   K84 = 0.000000E+00    ;
  122. *
  123.  MASQ1= TMPE MASQ 'EGALE'  TSTA1 ;
  124.   KCO21= (-1.* K11) * MASQ1   ;
  125.  MASQ1= TMPE MASQ 'EGALE'  TSTA2 ;
  126.   KCO22= (-1.* K21) * MASQ1   ;
  127.  MASQ1= TMPE MASQ 'EGALE'  TSTA3 ;
  128.   KCO23= (-1.* K31) * MASQ1   ;
  129.  MASQ1= TMPE MASQ 'EGALE'  TSTA4 ;
  130.   KCO24= (-1.* K41) * MASQ1   ;
  131.  MASQ1= TMPE MASQ 'EGALE'  TSTA5 ;
  132.   KCO25= (-1.* K51) * MASQ1   ;
  133.  MASQ1= TMPE MASQ 'EGALE'  TSTA6 ;
  134.   KCO26= (-1.* K61) * MASQ1   ;
  135.  MASQ1= TMPE MASQ 'EGALE'  TSTA7 ;
  136.   KCO27= (-1.* K71) * MASQ1   ;
  137.  MASQ1= TMPE MASQ 'EGALE'  TSTA8 ;
  138.   KCO28= (-1.* K81) * MASQ1   ;
  139.   kco2= kco21 + kco22 +kco23 +kco24+kco25+kco26+kco27+kco28 ;
  140.    TMPE2= TMPE - TSTA7 ;
  141.    TSTAD = TSTA8 - TSTA7 ;
  142.    KCOCC= MANU CHPO MT 1 SCAL (COO - K71) ;
  143.    MASQ1= TMPE2 MASQ 'INFERIEUR'  TSTAD ;
  144.    TMPE1= MASQ1* TMPE2 ;
  145.    KCOC1= MASQ1 * KCOCC ;
  146.    MASQ1= TMPE1 MASQ 'SUPERIEUR' 0.    ;
  147.    TMPE2= MASQ1* TMPE1 ;
  148.    KCOCC= MASQ1* KCOC1 ;
  149.    KCOC1=   ((((K74*TMPE2)+K73)*TMPE2)+K72)*TMPE2 ;
  150.    KCOI7= KCOCC - KCOC1 ;
  151.    TMPE2= TMPE - TSTA6 ;
  152.    TSTAD = TSTA7 - TSTA6 ;
  153.    KCOCC= MANU CHPO MT 1 SCAL (COO - K61) ;
  154.    MASQ1= TMPE2 MASQ 'INFERIEUR'  TSTAD ;
  155.    TMPE1= MASQ1* TMPE2 ;
  156.    KCOC1= MASQ1 * KCOCC ;
  157.    MASQ1= TMPE1 MASQ 'SUPERIEUR' 0.    ;
  158.    TMPE2= MASQ1* TMPE1 ;
  159.    KCOCC= MASQ1* KCOC1 ;
  160.    KCOC1=   ((((K64*TMPE2)+K63)*TMPE2)+K62)*TMPE2 ;
  161.    KCOI6= KCOCC - KCOC1 ;
  162.    TMPE2= TMPE - TSTA5 ;
  163.    TSTAD = TSTA6 - TSTA5 ;
  164.    KCOCC= MANU CHPO MT 1 SCAL (COO - K51) ;
  165.    MASQ1= TMPE2 MASQ 'INFERIEUR'  TSTAD ;
  166.    TMPE1= MASQ1* TMPE2 ;
  167.    KCOC1= MASQ1 * KCOCC ;
  168.    MASQ1= TMPE1 MASQ 'SUPERIEUR' 0.    ;
  169.    TMPE2= MASQ1* TMPE1 ;
  170.    KCOCC= MASQ1* KCOC1 ;
  171.    KCOC1=   ((((K54*TMPE2)+K53)*TMPE2)+K52)*TMPE2 ;
  172.    KCOI5= KCOCC - KCOC1 ;
  173.    TMPE2= TMPE - TSTA4 ;
  174.    TSTAD = TSTA5 - TSTA4 ;
  175.    KCOCC= MANU CHPO MT 1 SCAL (COO - K41) ;
  176.    MASQ1= TMPE2 MASQ 'INFERIEUR'  TSTAD ;
  177.    TMPE1= MASQ1* TMPE2 ;
  178.    KCOC1= MASQ1 * KCOCC ;
  179.    MASQ1= TMPE1 MASQ 'SUPERIEUR' 0.    ;
  180.    TMPE2= MASQ1* TMPE1 ;
  181.    KCOCC= MASQ1* KCOC1 ;
  182.    KCOC1=   ((((K44*TMPE2)+K43)*TMPE2)+K42)*TMPE2 ;
  183.    KCOI4= KCOCC - KCOC1 ;
  184.    TMPE2= TMPE - TSTA3 ;
  185.    TSTAD = TSTA4 - TSTA3 ;
  186.    KCOCC= MANU CHPO MT 1 SCAL (COO - K31) ;
  187.    MASQ1= TMPE2 MASQ 'INFERIEUR'  TSTAD ;
  188.    TMPE1= MASQ1* TMPE2 ;
  189.    KCOC1= MASQ1 * KCOCC ;
  190.    MASQ1= TMPE1 MASQ 'SUPERIEUR' 0.    ;
  191.    TMPE2= MASQ1* TMPE1 ;
  192.    KCOCC= MASQ1* KCOC1 ;
  193.    KCOC1=   ((((K34*TMPE2)+K33)*TMPE2)+K32)*TMPE2 ;
  194.    KCOI3= KCOCC - KCOC1 ;
  195.    TMPE2= TMPE - TSTA2 ;
  196.    TSTAD = TSTA3 - TSTA2 ;
  197.    KCOCC= MANU CHPO MT 1 SCAL (COO - K21) ;
  198.    MASQ1= TMPE2 MASQ 'INFERIEUR'  TSTAD ;
  199.    TMPE1= MASQ1* TMPE2 ;
  200.    KCOC1= MASQ1 * KCOCC ;
  201.    MASQ1= TMPE1 MASQ 'SUPERIEUR' 0.    ;
  202.    TMPE2= MASQ1* TMPE1 ;
  203.    KCOCC= MASQ1* KCOC1 ;
  204.    KCOC1=   ((((K24*TMPE2)+K23)*TMPE2)+K22)*TMPE2 ;
  205.    KCOI2= KCOCC - KCOC1 ;
  206.    TMPE2= TMPE - TSTA1 ;
  207.    TSTAD = TSTA2 - TSTA1 ;
  208.    KCOCC= MANU CHPO MT 1 SCAL (COO - K11) ;
  209.    MASQ1= TMPE2 MASQ 'INFERIEUR'  TSTAD ;
  210.    TMPE1= MASQ1* TMPE2 ;
  211.    KCOC1= MASQ1 * KCOCC ;
  212.    MASQ1= TMPE1 MASQ 'SUPERIEUR' 0.    ;
  213.    TMPE2= MASQ1* TMPE1 ;
  214.    KCOCC= MASQ1* KCOC1 ;
  215.    KCOC1=   ((((K14*TMPE2)+K13)*TMPE2)+K12)*TMPE2 ;
  216.    KCOI1= KCOCC - KCOC1 ;
  217.    KCOI= KCO2 + KCOI1 +KCOI2 +KCOI3+ KCOI4 +KCOI5+KCOI6+KCOI7 ;
  218.    KCOI= KCOI- 18.13 ;
  219.    KCO2= NOMC W016 KCOI ;
  220.   TBPAR2.CLIM= KCO2 ;
  221.  TB3= CHI2 TB1 TBPARM TBPAR2 ;
  222. *
  223. *
  224. *              controle des résultats
  225. *
  226.  PRECAB=REDU AB TB3.PREC ;
  227.  PRECCD=REDU CD TB3.PREC ;
  228.  PRECEF=REDU EF TB3.PREC ;
  229.  PREC1= MANU CHPO AB 14 W018 0. W019 2.83058E-04 W020  2.49697E-02
  230.   W021 0. W022 0. W023 0. W024 0. W025 0. W026 0. W027 0. W028 0.
  231.   W029 0. W030 0. W031 0. 'NATURE' DISCRET ;
  232.  PREC2= MANU CHPO CD 14 W018 0. W019 2.84979E-04 W020  2.67385E-02
  233.   W021 0. W022 0. W023 0. W024 0. W025 1.84781E-05 W026 0. W027 0.
  234.  W028 0.  W029 0. W030 0. W031 0. 'NATURE' DISCRET ;
  235.  PREC3= MANU CHPO EF 14 W018 0. W019 0. W020  2.78445E-02
  236.   W021 2.74960E-05 W022 0. W023 0. W024 0. W025 0. W026 0. W027 0.
  237.  W028 0.  W029 0. W030 0. W031 0. 'NATURE' DISCRET ;
  238. *PRECMT= PRECAB+PRECCD+PRECEF ;
  239.  PRECMT= PRECAB+PRECCD ;
  240. *PRECTES=PREC1+PREC2+PREC3 ;
  241.  PRECTES=PREC1+PREC2 ;
  242.  PRECD= PRECTES/50. ;
  243.  VERR1= (ABS ( PRECTES - PRECMT)) MASQUE SUPERIEUR PRECD SOMME ;
  244.  FIONAB= REDU AB TB3.FION ;
  245.  FIONCD= REDU CD TB3.FION ;
  246.  FIONEF= REDU EF TB3.FION ;
  247. *FIONMT= FIONAB+ FIONCD+FIONEF ;
  248.  FIONMT= FIONAB+ FIONCD ;
  249.  FIONTE1= MANU CHPO AB 1 SCAL 5.57739E-04 'NATURE' DISCRET ;
  250.  FIONTE2= MANU CHPO CD 1 SCAL 8.16893E-04 'NATURE' DISCRET ;
  251.  FIONTE3= MANU CHPO EF 1 SCAL 1.78045E-03 'NATURE' DISCRET ;
  252. *FIONTES= FIONTE1+FIONTE2+FIONTE3;
  253.  FIONTES= FIONTE1+FIONTE2;
  254.  FIOND= FIONTES/2.E5 ;
  255.  VERR2= (ABS ( FIONTES - FIONMT)) MASQUE SUPERIEUR FIOND SOMME ;
  256.  TYP5AB= REDU AB TB3.TYP5 ;
  257.  TYP5CD= REDU CD TB3.TYP5 ;
  258.  TYP5EF= REDU EF TB3.TYP5 ;
  259.  TYP51=MANU CHPO AB 14 W018 2.04302E-02 W019 0. W020 0.
  260.  W021  4.07418E-03 W022 9.08803E-04 W023 1.41667E-01 W024 6.49819E-04
  261.  W025  2.92431E-02 W026  4.43300E-02 W027  2.49896E-01
  262.  W028 2.15059E-02  W029 4.66326E-07 W030 2.92533E-07 W031 3.37606E-04
  263.  'NATURE' DISCRET ;
  264.  TYP52=MANU CHPO CD 14 W018 1.42911E-01 W019 0. W020 0.
  265. W021 3.55161E-01 W022 2.35987E-02 W023 5.19600E-01 W024 4.09491E-02
  266.  W025 0. W026  5.85450E-01  W027 1.97493E-01  W028 1.39410E-01
  267.  W029 9.40715E-04 W030 2.15057E-04 W031 1.57713E-03 'NATURE' DISCRET ;
  268.  TYP53=MANU CHPO EF 14 W018 2.97852E-01 W019 9.53362E-02 W020 0.
  269.  W021 0. W022 2.24386E-02 W023 2.09538E-04  W024 6.91829E-02
  270.  W025 1.06065E-01 W026 7.06324E-02 W027 1.04954E-01 W028 1.33045E-01
  271.  W029 1.80148E-01 W030 2.17771E-04 W031 2.50476E-05 'NATURE' DISCRET ;
  272. *TYP5TES= TYP51+TYP52+TYP53 ;
  273.  TYP5TES= TYP51+TYP52 ;
  274. *TYP5MT=TYP5AB+TYP5CD+TYP5EF ;
  275.  TYP5MT=TYP5AB+TYP5CD ;
  276.  TYP5D= TYP5TES/50. ;
  277.  VERR3= (ABS ( TYP5TES - TYP5MT)) MASQUE SUPERIEUR TYP5D SOMME ;
  278.  VERR= VERR1+VERR2+VERR3 ;
  279.  SI (VERR EGA 0 ) ;
  280.       ERREUR 0 ;
  281.  SINON ;
  282.       ERREUR 5 ;
  283.  FINSI ;
  284. *
  285. fin ;                                                                   
  286.  
  287.  
  288.  
  289.  
  290.  
  291.  
  292.  
  293.  
  294.  
  295.  
  296.  
  297.  
  298.  

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