Télécharger la plaquette Cast3M


fr en

Télécharger carre3D_therper.dgibi

Retour à la liste

Numérotation des lignes : 

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

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