Télécharger la plaquette Cast3M


fr en

Télécharger carre3D_therper_scalpass.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  carre3D_therper_scalpass.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ******************************************************************
  5. * CALCUL DU CUBE   A CHOC                                        *
  6. * Modele: gaz mono-espece "thermally perfect"                    *
  7. * Gaz mono-espece "thermally perfect"                            *
  8. * Splitting de scalaire passifs                                  *
  9. * Controle de la symetrie                                        *
  10. *                                                                *
  11. * A. BECCANTINI,    DRN/DMT/SEMT/LTMF,     FEVRIER 2000          *
  12. *Modif, 10/07/01, syntaxe de PENT changée                        *
  13. *Modif, 25/09/02, syntaxe de KONV changée                        *
  14. ******************************************************************
  15.  
  16. 'OPTION'  'DIME' 3 'ELEM' 'CUB8' 'ISOV' 'SULI' 'ECHO' 1
  17.   'TRAC' 'X';
  18.  
  19. GRAPH   = FAUX ;
  20. * GRAPH   = VRAI ;
  21.  
  22. *
  23. *** Methodes possibles :
  24. *
  25. *   'VLH', 'SS'
  26. *
  27.  
  28. METO =  'VLH' ;
  29. * SAFFAC = safety facton on the CFL
  30. SAFFAC = 0.5 ;
  31. * LOGSO = second order reconstruction
  32. LOGSO = VRAI ;
  33. TFINAL = 10. ;
  34. NITER = 1 ;
  35. ZERO = 1.0D-12 ;
  36.  
  37.  
  38. ************ 
  39. * MAILLAGE *
  40. ************
  41.  
  42.  A1 = 0.0 0.0 0.0 ;
  43.  A2 = 0.0 0.0 5.0 ;
  44.  LAT12 = A1 'DROIT' 5 A2 ;
  45.  SUR1 = LAT12 'TRANSLATION' 5 (0.0 5.0 0.0) ;
  46.  DOMTOT  = 'VOLUME' SUR1 'TRANSLATION' 5 (5.0 0.0 0.0) ;
  47.  
  48.  MDNS = 'EULER' ;
  49.  $DOMTOT = 'MODE' DOMTOT MDNS ;
  50.  TABM = 'DOMA' $DOMTOT 'VF' ;
  51.  MDOMTOT = TABM . 'QUAF' ;
  52.  'ELIM' MDOMTOT  1.E-5 ;
  53.  
  54. * Element au centre
  55.  
  56.  PCEN = ('DOMA' $DOMTOT 'CENTRE') 'POIN' 'PROC' (2.5 2.5 2.5);
  57.  DCEN = DOMTOT 'ELEM' 'CONTENANT' PCEN ;
  58.  'ELIMINATION' 1.D-4 DCEN DOMTOT ;
  59.  
  60.  $DCEN  = 'MODE' DCEN MDNS  ;
  61.  TABCEN = 'DOMA' $DCEN 'VF' ;
  62.  MDCEN  = TABCEN . 'QUAF' ;
  63.  'ELIM' (MDOMTOT 'ET' MDCEN) 1.E-5 ;
  64.  
  65. 'SI' GRAPH ;
  66.  'TRACER' ('DOMA' $DOMTOT 'MAILLAGE') 'TITRE' ('CHAINE' 'Maillage ');
  67. 'FINSI' ;   
  68.  
  69. ***********************
  70. **** LA TABLE PGAZ ****
  71. ***********************
  72.  
  73. PGAZ = 'TABLE' ;
  74.  
  75. *
  76. **** Ordre des polynoms
  77. *
  78.  
  79. PGAZ . 'NORD' = 4 ;
  80.  
  81. *
  82. **** Especes qui sont dans les equations d'Euler
  83. *
  84.  
  85. PGAZ . 'ESPEULE' = 'MOTS' 'H2  ' 'O2  ' 'H2O ' ;
  86.  
  87. *
  88. **** Espece qui n'y est pas
  89. *
  90.  
  91.  
  92. PGAZ . 'ESPNEULE' = 'N2  ';
  93.  
  94. *
  95.  
  96. PGAZ .  'H2  ' = 'TABLE'  ;
  97. PGAZ .  'H2O ' = 'TABLE'  ;
  98. PGAZ .  'N2  ' = 'TABLE'  ;
  99. PGAZ .  'O2  ' = 'TABLE'  ;
  100.  
  101. *
  102. **** R (J/Kg/K)
  103. *
  104.  
  105. PGAZ .  'H2  ' . 'R' = 4130.0 ;
  106. PGAZ .  'H2O ' . 'R' = 461.4 ;
  107. PGAZ .  'N2  ' . 'R' = 296.8 ;
  108. PGAZ .  'O2  ' . 'R' = 259.8 ;
  109.  
  110.  
  111. *
  112. **** Regressions polynomials
  113. *
  114.  
  115. PGAZ .  'H2  ' . 'A' = 'PROG'  9834.91866 0.54273926 0.000862203836
  116.                                -2.37281455E-07 1.84701105E-11 ;
  117. PGAZ .  'H2O ' . 'A' = 'PROG' 1155.95625 0.768331151 -5.73129958E-05
  118.                               -1.82753232E-08 2.44485692E-12 ;
  119. PGAZ .  'N2  ' . 'A' = 'PROG' 652.940766 0.288239099 -7.80442298E-05
  120.                               8.78233606E-09 -3.05514485E-13 ;
  121. PGAZ .  'O2  ' . 'A' = 'PROG' 575.012333  0.350522002 -0.000128294865
  122.                               2.33636971E-08 -1.53304905E-12;
  123.  
  124. *
  125. **** "Enthalpies" (ou energies) de formations a OK (J/Kg)
  126. *
  127.  
  128. PGAZ .  'H2  ' . 'H0K' = -4.195D6 ;
  129. PGAZ .  'H2O ' . 'H0K' = -1.395D7 ;
  130. PGAZ .  'N2  ' . 'H0K' = -2.953D5 ;
  131. PGAZ .  'O2  ' . 'H0K' = -2.634D5 ;
  132.  
  133. * Scalaires passives
  134.  
  135. PGAZ . 'SCALPASS' = 'MOTS' 'TAU ' 'EPS ' ;
  136.  
  137. *
  138. *** Fin PGAZ
  139. *
  140. *
  141. **** ETAT CENTRE
  142. *
  143.  
  144. rog   = 1.0  ;
  145. retg  = .4291145555695540D+04 ;
  146. yh2g = 0.01 ;
  147. yo2g = 0.2 ;
  148. yh2og = 0.3 ;
  149. taug = 1.0 ;
  150. epsg = 2.0 ;
  151.  
  152.  
  153. rouxg  = 0.0 ;
  154. rouyg  = 0.0 ;
  155. rouzg  = 0.0 ;
  156.  
  157.  
  158. ROC1 = KCHT $DCEN 'SCAL' 'CENTRE' rog ;
  159. ROVC1 = KCHT $DCEN 'VECT' 'CENTRE' (rouxg rouyg rouzg);
  160. RETC1 = KCHT $DCEN 'SCAL' 'CENTRE' retg ;
  161. RYH2C1 = KCHT $DCEN 'SCAL' 'CENTRE' 'COMP' 'H2  '  (rog*yh2g) ;
  162. RYO2C1 = KCHT $DCEN 'SCAL' 'CENTRE' 'COMP' 'O2  '  (rog*yo2g) ;
  163. RYH2OC1 = KCHT $DCEN 'SCAL' 'CENTRE' 'COMP' 'H2O ' (rog*yh2og) ;
  164. RSCAC11 = 'KCHT' $DCEN 'SCAL' 'CENTRE' 'COMP' 'TAU' (rog*taug) ;
  165. RSCAC12 = 'KCHT' $DCEN 'SCAL' 'CENTRE' 'COMP' 'EPS' (rog*epsg) ;
  166.  
  167.  
  168. *
  169. *** Le reste
  170. *
  171.  
  172. rod   =  1.250D-1 ;
  173. retd  = .3598345082089522D+01 ;
  174. yh2d = 0.2 ;
  175. yo2d = 0.4 ;
  176. yh2od = 0.1 ;
  177. taud = 11.0 ;
  178. epsd = 21.01 ;
  179.  
  180. rouxd  = 0.0 ;
  181. rouyd  = 0.0 ;
  182. rouzd  = 0.0 ;
  183.  
  184. ROC2 = KCHT $DOMTOT 'SCAL' 'CENTRE' rod ;
  185. ROVC2 = KCHT $DOMTOT 'VECT' 'CENTRE' (rouxd rouyd rouzd) ;
  186. RETC2 = KCHT $DOMTOT 'SCAL' 'CENTRE' retd ;
  187. RYH2C2 = KCHT $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'H2  '  (rod * yh2d) ;
  188. RYO2C2 = KCHT $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'O2  '  (rod * yo2d) ;
  189. RYH2OC2 = KCHT $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'H2O ' (rod * yh2od) ;
  190. RSCAC21 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'TAU' (rod*taud) ;
  191. RSCAC22 = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'EPS' (rod*epsd) ;
  192.  
  193. RN  = 'KCHT'  $DOMTOT 'SCAL' 'CENTRE' ROC2 ROC1 ;
  194. GN  = 'KCHT'  $DOMTOT 'VECT' 'CENTRE' ROVC2 ROVC1 ;
  195. RETN = 'KCHT'  $DOMTOT 'SCAL' 'CENTRE' RETC2 RETC1 ;
  196. RSCAN = ('KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'TAU' RSCAC21 RSCAC11)
  197.   'ET'  ('KCHT' $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'EPS' RSCAC22 RSCAC12) ;
  198.  
  199. RYH2  =  KCHT $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'H2  ' RYH2C2 RYH2C1 ;
  200. RYO2  =  KCHT $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'O2  ' RYO2C2 RYO2C1 ;
  201. RYH2O =  KCHT $DOMTOT 'SCAL' 'CENTRE' 'COMP' 'H2O ' RYH2OC2 RYH2OC1 ;
  202.  
  203. RYN =  RYH2 'ET' RYO2 'ET' RYH2O ;
  204.  
  205. ********************************************************
  206. *** CREATION DE 'MODE' POUR GRAPHIQUER LE CHAMELEM   ***
  207. ********************************************************
  208.  
  209. MOD1     =  'MODELISER'  ('DOMA' $DOMTOT 'MAILLAGE') 'THERMIQUE' ;
  210.  
  211. *
  212. *** GRAPHIQUE DES C.I.
  213. *
  214.  
  215. 'SI' GRAPH ;
  216. * 'SI' FAUX ;
  217. *
  218. *** CREATION DE CHAMELEM
  219. *
  220.  
  221.   CHM_RN   =  'KCHA' $DOMTOT 'CHAM' RN ;
  222.   CHM_GN  =  'KCHA' $DOMTOT 'CHAM' GN ;
  223.   CHM_RETN  =  'KCHA' $DOMTOT 'CHAM' RETN ;
  224.   CHM_RYN   =  'KCHA' $DOMTOT 'CHAM' RYN ;
  225.   CHM_RSN  =  'KCHA' $DOMTOT 'CHAM' RSCAN ;
  226.   TRAC CHM_RN MOD1 'TITR'  ('CHAINE'  'RN at t=' 0.0);
  227.   TRAC CHM_GN MOD1 'TITR'  ('CHAINE'  'GN at t=' 0.0);
  228.   TRAC CHM_RETN MOD1 'TITR'  ('CHAINE'  'ROET at t=' 0.0);
  229.   TRAC CHM_RYN MOD1 'TITR'  ('CHAINE'  'RYN  at t=' 0.0);
  230.   TRAC CHM_RSN MOD1 'TITR'  ('CHAINE'  'RSCAN  at t=' 0.0);
  231.  
  232.  
  233. 'FINSI' ;
  234.  
  235.  RN0 = 'COPIER' RN ;
  236.  GN0 = 'COPIER' GN ;
  237.  RETN0 = 'COPIER' RETN ;
  238.  RYN0 = 'COPIER' RYN ;
  239.  RSCAN0 = 'COPIER' RSCAN ;
  240.  
  241. *
  242. **** Parameter for the time loop 
  243. *
  244.  
  245. * Names of the residuum components
  246.  
  247.  LISTINCO = 'MOTS' 'RN' 'RUX' 'RUY' 'RUZ' 'RETN' 'H2' 'O2' 'H2O'
  248.                    'TAU' 'EPS' ;
  249.  NOMRN = 'MOTS' 'SCAL' ;
  250.  NOMVN = 'MOTS' 'UX' 'UY' 'UZ' ;
  251.  NOMY  = 'MOTS' 'H2' 'O2' 'H2O' ;
  252.  NOMT  = 'MOTS' 'TAU' 'EPS' ;
  253.  
  254. *
  255. **** Geometrical coefficient to compute gradients
  256. *
  257.  
  258.  GRADR CACCA COEFSCAL =  'PENT' $DOMTOT 'CENTRE' 'EULESCAL' 'NOLIMITE'
  259.           NOMRN RN ;
  260.  
  261.  GRADV CACCA COEFVECT =  'PENT' $DOMTOT 'CENTRE' 'EULEVECT' 'NOLIMITE'
  262.           NOMVN GN ;
  263.  
  264.  TPS = 0.0 ;
  265.  
  266.  
  267. *
  268. **** Temporal loop
  269. *
  270.  
  271.  'MESSAGE' ;
  272.  'MESSAGE' ('CHAINE' 'Methode = ' METO) ;
  273.  'MESSAGE' ('CHAINE' 'SAFFAC      = ' SAFFAC) ;
  274.  'MESSAGE' ;
  275.  
  276.   VN PN TN YN SN GAMMAN = 'PRIM' 'PERFTEMP' PGAZ RN GN RETN RYN RSCAN ;
  277.  
  278.   TNM1 = 'COPIER' TN ;
  279.  
  280.  'TEMPS' 'ZERO' ;
  281.  
  282.  'REPETER' BL1 NITER ;
  283.  
  284. *
  285. **** Primitive variables
  286. *
  287.  
  288.     VN PN TN YN SN GAMMAN = 'PRIM' 'PERFTEMP' PGAZ RN GN RETN RYN RSCAN
  289.       TNM1  ;
  290.  
  291.     TNM1 = 'COPIER' TN ;
  292.  
  293.     'SI' LOGSO ;
  294.  
  295.        GRADR ALR = 'PENT' $DOMTOT 'CENTRE' 'EULESCAL' 'LIMITEUR'
  296.           NOMRN RN  'GRADGEO' COEFSCAL ;
  297.  
  298.        GRADP ALP = 'PENT' $DOMTOT 'CENTRE' 'EULESCAL' 'LIMITEUR'
  299.           NOMRN PN  'GRADGEO' COEFSCAL ;
  300.  
  301.        GRADV ALV = 'PENT' $DOMTOT 'CENTRE' 'EULEVECT' 'LIMITEUR'
  302.           NOMVN VN   'GRADGEO'  COEFVECT ;
  303.  
  304.        GRADY ALY = 'PENT' $DOMTOT 'CENTRE' 'EULESCAL' 'LIMITEUR'
  305.           NOMY YN  'GRADGEO' COEFSCAL ;
  306.  
  307.        GRADS ALS = 'PENT' $DOMTOT 'CENTRE' 'EULESCAL' 'LIMITEUR'
  308.           NOMT SN  'GRADGEO' COEFSCAL ;
  309.  
  310.        ROF VITF PF YF SF = 'PRET' 'PERFTEMP'  2 1
  311.                             $DOMTOT  PGAZ
  312.                             RN    GRADR  ALR
  313.                             VN    GRADV  ALV
  314.                             PN    GRADP  ALP
  315.                             YN    GRADY  ALY
  316.                             SN    GRADS  ALS
  317.                             ;
  318.  
  319.     'SINON' ;
  320.  
  321.        ROF VITF PF YF SF = 'PRET' 'PERFTEMP'  1 1 
  322.                               $DOMTOT  PGAZ
  323.                               RN   
  324.                               VN   
  325.                               PN   
  326.                               YN
  327.                               SN ;
  328.  
  329.     'FINSI' ;
  330.  
  331.     RESIDU DELTAT =  'KONV' 'VF' 'PERFTEMP' 'RESI' METO
  332.          $DOMTOT PGAZ LISTINCO  ROF VITF PF YF SF ;
  333.  
  334.     DT_CON = SAFFAC '*' DELTAT ;     
  335. *
  336. **** The time step linked to tps
  337. *
  338.  
  339.     DTTPS = (TFINAL '-' TPS) * (1. '+' ZERO) ;
  340.  
  341. *
  342. **** Total time step
  343. *
  344.  
  345.     DTMIN = 'MINIMUM' ('PROG' DT_CON DTTPS) ;
  346.  
  347. *
  348. **** Increment of the variables (convection)
  349. *
  350.  
  351.     RESIDU = DTMIN '*' RESIDU ;
  352.  
  353.     DRN = 'EXCO' 'RN' RESIDU 'SCAL' ;
  354.     DGN = 'EXCO' ('MOTS' 'RUX' 'RUY' 'RUZ') RESIDU
  355.                  ('MOTS' 'UX' 'UY' 'UZ') ;
  356.     DRETN = 'EXCO' 'RETN' RESIDU 'SCAL' ;
  357.     DRYN  = 'EXCO' ('MOTS' 'H2' 'O2' 'H2O') RESIDU
  358.             ('MOTS' 'H2' 'O2' 'H2O') ;
  359.     DRSN  = 'EXCO' ('MOTS' 'TAU' 'EPS') RESIDU
  360.             ('MOTS' 'TAU' 'EPS') ;
  361.  
  362.     TPS  = TPS '+' DTMIN ;
  363.     RN   = RN '+' DRN ;
  364.     GN   = GN '+' DGN ;
  365.     RETN = RETN '+' DRETN ;
  366.     RYN = RYN '+' DRYN ;
  367.     RSCAN = RSCAN '+' DRSN ;
  368.  
  369.     'SI' (((&BL1 '/' 20) '*' 20) 'EGA' &BL1) ;
  370.         'MESSAGE' ('CHAINE' 'ITER =' &BL1 '  TPS =' TPS) ;
  371.     'FINSI' ;
  372.  
  373.     'SI' (TPS > TFINAL) ;
  374.        'QUITTER' BL1 ;
  375.     'FINSI' ;
  376.  
  377.  'FIN' BL1 ;
  378.  'TEMPS' ;
  379.  
  380. *
  381. **** Les variables primitives
  382. *
  383.  
  384. VN PN TN YN SN GAMN =   'PRIM' 'PERFTEMP' PGAZ
  385.            RN GN RETN RYN RSCAN ;
  386.  
  387. CN = (GAMN * (PN / RN)) '**' 0.5 ;
  388.  
  389. *
  390. *** GRAPHIQUE DES SOLUTIONS
  391. *
  392.  
  393. 'SI' GRAPH ;
  394. *
  395. *** CREATION DE CHAMELEM
  396. *
  397.   CHM_RN   =  'KCHA' $DOMTOT 'CHAM' RN ;
  398.   CHM_GN   =  'KCHA' $DOMTOT 'CHAM' GN ;
  399.   CHM_RETN  =  'KCHA' $DOMTOT 'CHAM' RETN ;
  400.   CHM_RYN  =  'KCHA' $DOMTOT 'CHAM' RYN ;  
  401.   TRAC CHM_RN  MOD1 'TITR'  ('CHAINE'  'RO at t=' TFIN);
  402.   TRAC CHM_GN  MOD1 'TITR'  ('CHAINE'  'GN at t=' TFIN);
  403.   TRAC CHM_RETN MOD1 'TITR'  ('CHAINE'  'ROET at t=' TFIN);
  404.   TRAC CHM_RYN MOD1 'TITR'  ('CHAINE'  'RYN at t=' TFIN);
  405. 'FINSI' ;
  406.  
  407.  
  408.  
  409. * Tests: invariance par rotation
  410.  
  411.  RN1 = 'COPIER'  RN ;
  412.  GN1 = 'COPIER'  GN ;
  413.  RETN1 = 'COPIER'  RETN ;
  414.  RYN1 = 'COPIER'  RYN ;
  415.  RSCAN1 = 'COPIER'  RSCAN ;
  416.  
  417.  MAIL0 = 'EXTRAIRE'  RN 'MAILLAGE' ;
  418.  
  419. * Axe de rotation: PCEN = 2.5 2.5 2.5 '+' VECX ;
  420.  
  421.  VECX = 1.0 0.0 0.0 ;
  422.  
  423.  'REPETER'  BL1 3 ;
  424.  
  425.     RN1 = RN1 'TOURNER'  90 PCEN (PCEN 'PLUS' VECX) ;
  426.     GN1 = GN1 'TOURNER'  90  PCEN (PCEN 'PLUS' VECX)  ;
  427.     RETN1 = RETN1 'TOURNER'  90   PCEN (PCEN 'PLUS' VECX)  ;
  428.     RYN1 = RYN1 'TOURNER'  90  PCEN (PCEN 'PLUS' VECX)  ;
  429.     RSCAN1 = RSCAN1 'TOURNER'  90  PCEN (PCEN 'PLUS' VECX)  ;
  430.  
  431.     MAIL1 = 'EXTRAIRE'  RN1 'MAILLAGE' ;
  432.     'ELIMINATION'  0.001 MAIL0  MAIL1 ;
  433.     ERR1 = 'MAXIMUM'  ((RN1 '-' RN) / RN) 'ABS' ;
  434.     'SI'  (ERR1 > 1.0D-10);
  435.        'ERREUR'  5;
  436.     'FINSI'  ;
  437.  
  438.     MAIL1 = 'EXTRAIRE'  GN1 'MAILLAGE' ;
  439.     'ELIMINATION'  0.001 MAIL0 MAIL1 ;
  440.     ERR1 = 'MAXIMUM'  ((GN1 '-' GN) / (RN * CN)) 'ABS' ;
  441.     'SI'  (ERR1 > 1.0D-10);
  442.        'ERREUR'  5;
  443.     'FINSI'  ;
  444.  
  445.     MAIL1 = 'EXTRAIRE'  RETN1 'MAILLAGE' ;
  446.     'ELIMINATION'  0.001 MAIL0 MAIL1 ;
  447.     ERR1 = 'MAXIMUM'  ((RETN1 '-' RETN) / RETN) 'ABS' ;
  448.     'SI'  (ERR1 > 1.0D-10);
  449.        'ERREUR'  5;
  450.     'FINSI'  ;
  451.  
  452.     MAIL1 = 'EXTRAIRE'  RYN1 'MAILLAGE' ;
  453.     'ELIMINATION'  0.001 MAIL0 MAIL1 ;
  454.     ERR1 = 'MAXIMUM'  ((RYN1 '-' RYN) / RN) 'ABS' ;
  455.     'SI'  (ERR1 > 1.0D-10);
  456.        'ERREUR'  5;
  457.     'FINSI'  ;
  458.  
  459.     MAIL1 = 'EXTRAIRE'  RSCAN1 'MAILLAGE' ;
  460.     'ELIMINATION'  0.001 MAIL0 MAIL1 ;
  461.     ERR1 = 'MAXIMUM'  ((RSCAN1 '-' RSCAN) / RN) 'ABS' ;
  462.  
  463.     'SI'  (ERR1 > 1.0D-10);
  464.       'ERREUR'  5;
  465.     'FINSI'  ;
  466.  
  467.  'FIN'  BL1 ;
  468.  
  469.  VECX = 0.0 0.0 1.0 ;
  470.  
  471.  'REPETER'  BL1 3 ;
  472.  
  473.     RN1 = RN1 'TOURNER'  90 PCEN (PCEN 'PLUS' VECX) ;
  474.     GN1 = GN1 'TOURNER'  90  PCEN (PCEN 'PLUS' VECX)  ;
  475.     RETN1 = RETN1 'TOURNER'  90   PCEN (PCEN 'PLUS' VECX)  ;
  476.     RYN1 = RYN1 'TOURNER'  90  PCEN (PCEN 'PLUS' VECX)  ;
  477.     RSCAN1 = RSCAN1 'TOURNER'  90  PCEN (PCEN 'PLUS' VECX)  ;
  478.  
  479.     MAIL1 = 'EXTRAIRE'  RN1 'MAILLAGE' ;
  480.     'ELIMINATION'  0.001 MAIL0  MAIL1 ;
  481.     ERR1 = 'MAXIMUM'  ((RN1 '-' RN) / RN) 'ABS' ;
  482.     'SI'  (ERR1 > 1.0D-10);
  483.        'ERREUR'  5;
  484.     'FINSI'  ;
  485.  
  486.     MAIL1 = 'EXTRAIRE'  GN1 'MAILLAGE' ;
  487.     'ELIMINATION'  0.001 MAIL0 MAIL1 ;
  488.     ERR1 = 'MAXIMUM'  ((GN1 '-' GN) / (RN * CN)) 'ABS' ;
  489.     'SI'  (ERR1 > 1.0D-10);
  490.        'ERREUR'  5;
  491.     'FINSI'  ;
  492.  
  493.     MAIL1 = 'EXTRAIRE'  RETN1 'MAILLAGE' ;
  494.     'ELIMINATION'  0.001 MAIL0 MAIL1 ;
  495.     ERR1 = 'MAXIMUM'  ((RETN1 '-' RETN) / RETN) 'ABS' ;
  496.     'SI'  (ERR1 > 1.0D-10);
  497.        'ERREUR'  5;
  498.     'FINSI'  ;
  499.  
  500.     MAIL1 = 'EXTRAIRE'  RYN1 'MAILLAGE' ;
  501.     'ELIMINATION'  0.001 MAIL0 MAIL1 ;
  502.     ERR1 = 'MAXIMUM'  ((RYN1 '-' RYN) / RN) 'ABS' ;
  503.     'SI'  (ERR1 > 1.0D-10);
  504.        'ERREUR'  5;
  505.     'FINSI'  ;
  506.  
  507.     MAIL1 = 'EXTRAIRE'  RSCAN1 'MAILLAGE' ;
  508.     'ELIMINATION'  0.001 MAIL0 MAIL1 ;
  509.     ERR1 = 'MAXIMUM'  ((RSCAN1 '-' RSCAN) / RN) 'ABS' ;
  510.  
  511.     'SI'  (ERR1 > 1.0D-10);
  512.        'ERREUR'  5;
  513.     'FINSI'  ;
  514.  
  515.  'FIN'  BL1 ;
  516.  
  517. * Symetrie
  518.  
  519.  MAIL1 RN1 = MAIL0 RN 'SYMETRIE'
  520.           'PLAN' (2.5 2.5 2.5) ((2.5 2.5 2.5) 'PLUS' (0.0 1.0 0.0))
  521.                  ((2.5 2.5 2.5) 'PLUS' (0.0 0.0 1.0)) ;
  522.  'ELIMINATION'  0.001  MAIL0  MAIL1 ;
  523.  ERR1 = 'MAXIMUM'  ((RN1 - RN) / RN) 'ABS' ;
  524.  
  525.  'SI'  (ERR1 > 1.0D-10);
  526.     'ERREUR'  5;
  527.  'FINSI'  ;
  528.  
  529.  'FIN'  ;
  530.  
  531.  
  532.  
  533.  
  534.  
  535.  
  536.  
  537.  
  538.  
  539.  
  540.  
  541.  
  542.  
  543.  

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