Télécharger la plaquette Cast3M


fr en

Télécharger trkg2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  trkg.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *  Enceinte 2D Axisymetrique Cas de Torrance et Rockett
  6. *
  7. 'OPTI' 'DIME' 2 'ELEM' 'CU20' ;
  8.  'OPTI' 'MODE' 'AXIS' ;
  9.  
  10.  opti trac 'PSC' ;
  11.  COMPLET= VRAI                                          ;
  12.  COMPLET= FAUX                                          ;
  13. *GRAPH  = VRAI                                          ;
  14.  GRAPH  = FAUX                                          ;
  15.  
  16.  echl=2.e-1;
  17. *echl=5.e-1;
  18. *echl=1.e-1;
  19.  'SI' COMPLET ;
  20.   nbit=100    ;
  21.   nbit=20     ;
  22.   DT0 = 1.    ;
  23.   ds1= 0.02*echl;
  24.   ds2= 0.2*echl ;
  25.   n1=5 ;
  26.   n2=10;
  27.   n3=20;
  28.  'SINON'      ;
  29.   nbit=100   ;
  30. * nbit=2    ;
  31.   DT0 = 40.;
  32.   ds1= 0.02*echl;
  33.   ds2= 0.2*echl ;
  34.   n1=5 ;
  35.   n2=10;
  36.   n3=20;
  37.  'FINSI'      ;
  38.  
  39. r0=0.1*echl ;
  40. r1=1.*echl ;
  41. rm=0.5*echl ;
  42. hm=0.5*echl ;
  43. h1=1.*echl ;
  44.  
  45. eps=1.e-2*echl ;
  46. eps=1.e-4*echl ;
  47. epsi=1.e-4*echl ;
  48.  
  49. p1=eps 0 ;
  50. p11=r0 0 ;
  51. p2=rm 0. ;
  52. p3=r1 0 ;
  53. p4=r1 hm;
  54. p5=r1 h1 ;
  55.  
  56.  
  57. breche = p1  'DROI' n1 p11;
  58. bas= p11 'DROI' dini ds1 dfin ds2 p2 d dini ds2 dfin ds1  p3 ;
  59. paroid=p3 'DROI' dini ds1 dfin ds2 p4 'DROI' dini ds2 dfin ds1 p5 ;
  60.  
  61. plaf=inve ( (breche et bas ) plus (0 h1));
  62. axe =inve ( paroid moins ((r1-eps) 0)) ;
  63. elim (breche et bas et paroid et plaf et axe) epsi ;
  64. mt= 'DALLER' (breche 'ET' bas) paroid plaf axe ;
  65. Si GRAPH ; trace mt ; Finsi ;
  66.  
  67. rxt = 'TABLE' ;
  68. rxt.'VERSION'= 'V0'                                           ;
  69. rxt.'vtf'=mt ;
  70. rxt.'axe'   = axe ;
  71. rxt.'epsi'= epsi ;
  72. rxt.'pi'  = 0.5 0.5 ;
  73. rxt.'DT0'  = DT0     ;
  74. rxt.'DISCR'= 'QUAF';
  75. rxt.'KPRE' = 'MSOMMET';
  76. rxt.'TYPINV'=1 ;
  77. rxt.'GRAPH'  = GRAPH ;
  78. rxt.'MODTURB'='NUTURB'  ;
  79. rxt.'NUT'    =1.e-7  ;
  80. rxt.'CORTEMP'=FAUX ;
  81.  
  82.  rxt.'TF0'    =        20.0                         ;
  83.  rxt.'PT0'    =        1.0e5                       ;
  84.  
  85.  rxt.'TIMP'              = table                              ;
  86.  rxt.'TIMP'.'TIMP1'             = table                       ;
  87.  rxt.'TIMP'.'TIMP1'.'MAILLAGE'  = breche                      ;
  88.  rxt.'TIMP'.'TIMP1'.'t'         = prog   0.0 900000.0         ;
  89.  rxt.'TIMP'.'TIMP1'.'TIMP'      = prog  21.     21.           ;
  90.  rxt.'TIMP'.'TIMP1'.'ECHAN'     = 1.e2                        ;
  91.  
  92.  rxt.'TIMP'.'TIMP2'             = table                       ;
  93.  rxt.'TIMP'.'TIMP2'.'MAILLAGE'  = (bas et paroid et plaf)     ;
  94.  rxt.'TIMP'.'TIMP2'.'t'         = prog   0.0 900000.0         ;
  95.  rxt.'TIMP'.'TIMP2'.'TIMP'      = prog  20.0  20.0            ;
  96.  rxt.'TIMP'.'TIMP2'.'ECHAN'     = 1.e2                        ;
  97.  
  98.  
  99. execrxt nbit rxt ;
  100.  
  101.  
  102. $vtf=rxt.'GEO'.'$vtf' ;
  103. vtf=doma $vtf maillage;
  104. rho=rxt.'TIC'.'RHO' ;
  105. tf =rxt.'TIC'.'TF' ;
  106. un =rxt.'TIC'.'UN';
  107.  
  108. $axe=rxt.'GEO'.'$axe';
  109. axe=doma $axe maillage ;
  110. ixe=inve axe ;
  111. evauz = evol chpo (exco un 'UY') ixe ;
  112. evarh = evol chpo rho ixe ;
  113. evatf = evol chpo tf  ixe ;
  114.  
  115.  
  116. luz= (extr evauz 'ORDO' );
  117. lrh= (extr evarh 'ORDO' );
  118. ltf= (extr evatf 'ORDO' );
  119. list luz ; list lrh ; list ltf ;
  120.  
  121. luzr= prog
  122. *  3.07237E-42 -2.34123E-03  2.56887E-03  6.54417E-03  1.62324E-02
  123. *  2.40888E-02  2.50937E-02  2.65036E-02  2.59309E-02  2.53052E-02
  124. *  2.60378E-02  2.60663E-02  2.63601E-02  2.57449E-02  2.55996E-02
  125. *  2.49561E-02  2.43348E-02  2.32284E-02  2.10423E-02  1.78512E-02
  126. *  1.39740E-02  8.69525E-03  3.46774E-03  7.10408E-04 -5.08355E-43;
  127.   1.51890E-42 -2.37360E-03  2.32262E-03  6.07648E-03  1.56239E-02
  128.   2.33899E-02  2.47391E-02  2.65114E-02  2.61377E-02  2.59326E-02
  129.   2.69223E-02  2.75933E-02  2.85105E-02  2.87585E-02  2.92365E-02
  130.   2.90317E-02  2.83702E-02  2.68461E-02  2.42700E-02  2.02241E-02
  131.   1.62119E-02  1.04393E-02  5.00558E-03  9.75119E-04  5.53846E-42 ;
  132.  
  133. lrhr=prog
  134. *   1.1842       1.1845       1.1848       1.1852       1.1855
  135. *   1.1859       1.1862       1.1865       1.1868       1.1870    
  136. *   1.1872       1.1874       1.1875       1.1876       1.1877
  137. *   1.1877       1.1878       1.1878       1.1878       1.1878    
  138. *   1.1879       1.1879       1.1879       1.1880       1.1882;
  139.    1.1842       1.1845       1.1848       1.1852       1.1855
  140.    1.1859       1.1862       1.1865       1.1868       1.1870
  141.    1.1872       1.1874       1.1875       1.1876       1.1876
  142.    1.1877       1.1877       1.1877       1.1878       1.1878
  143.    1.1878       1.1878       1.1879       1.1880       1.1882  ; 
  144.  
  145. ltfr=prog
  146. *   20.994       20.905       20.829       20.740       20.651
  147. *   20.565       20.482       20.404       20.337       20.275
  148. *   20.223       20.177       20.142       20.115       20.095
  149. *   20.082       20.072       20.065       20.058       20.054
  150. *   20.050       20.045       20.036       20.020       20.001;
  151.    20.994       20.906       20.831       20.742       20.655
  152.    20.570       20.489       20.414       20.349       20.292
  153.    20.244       20.202       20.171       20.148       20.132
  154.    20.121       20.113       20.107       20.102       20.099
  155.    20.095       20.091       20.078       20.048       20.003  ;
  156.  
  157. ERUZ=(SOMM (abs (luz - luzr)))*0.04 ;
  158. MESS 'ERUZ=' ERUZ ;
  159. Si (ERUZ '>' 1.e-3) ; erreur 5 ; Finsi ;
  160.  
  161. ERRH=(SOMM (abs (lrh - lrhr)))*0.04 ;
  162. MESS 'ERRH=' ERRH ;
  163. Si (ERRH '>' 1.e-3) ; erreur 5 ; Finsi ;
  164.  
  165. ERTF=(SOMM (abs (ltf - ltfr)))*0.04/20. ;
  166. MESS 'ERTF=' ERTF ;
  167. Si (ERTF '>' 1.e-3) ; erreur 5 ; Finsi ;
  168.  
  169.  
  170.  
  171.  
  172.  
  173. Si GRAPH ;
  174.  
  175.  dess evauz 'TITRE' 'Vitesse sur l axe';
  176.  dess evarh 'TITRE' 'Densite melange sur l axe';
  177.  dess evatf 'TITRE' 'Temperature sur l axe';
  178.  
  179. ung= vect un 1.  ux uy jaune;
  180. trace ung vtf;
  181. opti isov suli ;
  182. trace rho vtf 'TITRE' ' Rho' ;
  183. trace tf  vtf (prog 20. pas 0.001  20.02);
  184. trace tf  vtf (prog 20. 20.0001 20.001 20.005 20.01 20.1) ;
  185. *trace rh2 vtf 'TITRE' ' Rh2' ;
  186. *trace rco2 vtf 'TITRE' ' Rco2' ;
  187. *trace tf  vtf 'TITRE' ' Tf ' ;
  188. Finsi ;
  189.  
  190.  FIN ;
  191.  
  192.  
  193.  
  194.  
  195.  
  196.  
  197.  
  198.  
  199.  
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  

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