Télécharger la plaquette Cast3M


fr en

Télécharger fimpvf.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  fimpvf.dgibi
  2. ***********************************************************
  3. **** APPROCHE VF "Cell-Centred Formulation" pour la    ****
  4. **** solution des                                      **** 
  5. **** Equations d'Euler pour un gaz parfait             ****
  6. **** OPERATEUR  FIMP                                   ****
  7. **** Test de la gravité                                ****
  8. ****                                                   ****
  9. **** A. BECCANTINI DEN/DM2S/SFME/LTMF   JANV  2003     ****
  10. ***********************************************************
  11.  
  12.  
  13.  'OPTION'  'DIME' 2 'ELEM' QUA4 'ECHO'  0 ;
  14.  
  15. ***********************
  16. **** TEST1 ************
  17. ***********************
  18. * 2D
  19.  
  20.  
  21.  A1 = 0.0D0 0.0D0;
  22.  A2 = 2.0D0 0.0D0;
  23.  A3 = 3.0D0 0.0D0;
  24.  A4 = 4.0D0 1.0D0;
  25.  A5 = 1.0D0 1.0D0;
  26.  A6 = 0.0D0 1.0D0;
  27.  
  28.  MAI1 = ('MANUEL' 'POI1' A1) 'ET'
  29.         ('MANUEL' 'POI1' A2) 'ET'
  30.         ('MANUEL' 'POI1' A3) 'ET'
  31.         ('MANUEL' 'POI1' A4) 'ET'
  32.         ('MANUEL' 'POI1' A5) 'ET'
  33.         ('MANUEL' 'POI1' A6) ;
  34.  
  35. *
  36. **** CHPOINTs definis sur la meme rapresentation du maillage
  37. *    
  38.  
  39.  RN    = 'BRUI' 'BLAN' 'UNIF' 1.24 0.5  MAI1 ;
  40.  VNY   = ('NOMC' 'UY'
  41.         ('BRUI' 'BLAN' 'UNIF' 2.24 7.5  MAI1)
  42.         'NATU' 'DISCRET') ; 
  43.  VNX   = ('NOMC' 'UX'
  44.         ('BRUI' 'BLAN' 'UNIF' 7.24 0.5  MAI1)
  45.         'NATU' 'DISCRET') ;
  46.  VN    = VNX 'ET' VNY ;
  47.  GN    = RN '*' VN ;
  48.  
  49.  GY   = ('NOMC' 'UY'
  50.         ('BRUI' 'BLAN' 'UNIF' 1.24 7.5  MAI1)
  51.         'NATU' 'DISCRET') ; 
  52.  GX   = ('NOMC' 'UX'
  53.         ('BRUI' 'BLAN' 'UNIF' 0.24 0.5  MAI1)
  54.         'NATU' 'DISCRET') ;
  55.  
  56.  GRAV  = GX 'ET' GY ;
  57.  
  58.  LISTINCO = 'MOTS' 'RN' 'RUX' 'RUY' 'RET' ;
  59.  
  60.  RES  = ('NOMC' ('MOTS' 'UX' 'UY') (RN * GRAV) ('MOTS' 'RUX' 'RUY')
  61.          'NATU' 'DISCRET') 'ET'
  62.         ('NOMC' 'RET' ('PSCAL' GN GRAV ('MOTS' 'UX' 'UY')
  63.          ('MOTS' 'UX' 'UY')) 'NATU' 'DISCRET') ;
  64.  
  65.  RES1 = 'FIMP' 'VF' 'GRAVMONO' 'RESI' LISTINCO RN GN GRAV ;
  66.  
  67.  ERRES  = ('MAXIMUM' (RES '-' RES1) 'ABS') ;
  68.  
  69.  'SI' (ERRES '>' 1.0D-12);
  70.     'ERREUR'  5;
  71.  'FINSI' ;
  72.  
  73. * Jacobienne
  74.  
  75.  JAC = 'FIMP' 'VF' 'GRAVMONO' 'JACOCONS' LISTINCO RN GN GRAV ;
  76.  CHPO1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RN' 1.0 ;
  77.  
  78.  CHPRE = 'KOPS' 'MULT' JAC CHPO1 ;
  79.  
  80.  CHPRE1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 2 'RUX'
  81.      ('EXTRAIRE' GRAV A4 'UX') 'RUY'
  82.      ('EXTRAIRE' GRAV A4 'UY') ;
  83.  
  84.  ERRES  = ('MAXIMUM' (CHPRE1 '-' CHPRE) 'ABS') ;
  85.  
  86.  'SI' (ERRES '>' 1.0D-12);
  87.     'ERREUR'  5;
  88.  'FINSI' ;
  89.  
  90.  CHPO1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RUX' 1.0 '+'
  91.   ('MANUEL' 'CHPO' MAI1  1 'RUX' 0.0) ;
  92.  
  93.  CHPRE = 'KOPS' 'MULT' JAC CHPO1 ;
  94.  
  95.  CHPRE1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RET'
  96.      ('EXTRAIRE' GRAV A4 'UX') ;
  97.  
  98.  ERRES  = ('MAXIMUM' (CHPRE1 '-' CHPRE) 'ABS') ;
  99.  
  100.  'SI' (ERRES '>' 1.0D-12);
  101.     'ERREUR'  5;
  102.  'FINSI' ;
  103.  
  104.  CHPO1 = ('MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RUY' 1.0) '+'
  105.   ('MANUEL' 'CHPO' MAI1  1 'RUY' 0.0) ;
  106.  
  107.  CHPRE = 'KOPS' 'MULT' JAC CHPO1 ;
  108.  
  109.  CHPRE1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RET'
  110.      ('EXTRAIRE' GRAV A4 'UY') ;
  111.  
  112.  ERRES  = ('MAXIMUM' (CHPRE1 '-' CHPRE) 'ABS') ;
  113.  
  114.  'SI' (ERRES '>' 1.0D-12);
  115.     'ERREUR'  5 ;
  116.  'FINSI' ;
  117.  
  118. ***********************
  119. **** TEST2 ************
  120. ***********************
  121. *
  122. *** 3D
  123. *
  124.  
  125.  'OPTION'  'DIME' 3 'ELEM' 'CUB8' 'ECHO'  1 ;
  126.  
  127. *
  128. **** GRAPH
  129. *
  130.  
  131.  GRAPH = FAUX ;
  132. * GRAPH = VRAI ;
  133.  
  134.  A1 = 0.0D0 0.0D0 0.0 ;
  135.  A2 = 2.0D0 0.0D0 0.0 ;
  136.  A3 = 3.0D0 0.0D0 0.0 ;
  137.  A4 = 4.0D0 1.0D0 0.0 ;
  138.  A5 = 1.0D0 1.0D0 0.0 ;
  139.  A6 = 0.0D0 1.0D0 0.0 ;
  140.  
  141.  MAI1 = ('MANUEL' 'POI1' A1) 'ET'
  142.         ('MANUEL' 'POI1' A2) 'ET'
  143.         ('MANUEL' 'POI1' A3) 'ET'
  144.         ('MANUEL' 'POI1' A4) 'ET'
  145.         ('MANUEL' 'POI1' A5) 'ET'
  146.         ('MANUEL' 'POI1' A6) ;
  147.  
  148. *    
  149.  
  150.  RN    = 'BRUI' 'BLAN' 'UNIF' 1.24 0.5  MAI1 ;
  151.  VNY   = ('NOMC' 'UY'
  152.         ('BRUI' 'BLAN' 'UNIF' 4.24 7.5  MAI1)
  153.         'NATU' 'DISCRET') ; 
  154.  VNX   = ('NOMC' 'UX'
  155.         ('BRUI' 'BLAN' 'UNIF' 0.24 0.5  MAI1)
  156.         'NATU' 'DISCRET') ;
  157.  VNZ   = ('NOMC' 'UZ'
  158.         ('BRUI' 'BLAN' 'UNIF' 1.74 0.5  MAI1)
  159.         'NATU' 'DISCRET') ;
  160.  VN    = VNX 'ET' VNY 'ET' VNZ ;
  161.  GN    = RN '*' VN ;
  162.  
  163.  GY   = ('NOMC' 'UY'
  164.         ('BRUI' 'BLAN' 'UNIF' 3.24 7.5  MAI1)
  165.         'NATU' 'DISCRET') ; 
  166.  GX   = ('NOMC' 'UX'
  167.         ('BRUI' 'BLAN' 'UNIF' 2.24 0.5  MAI1)
  168.         'NATU' 'DISCRET') ;
  169.  GZ   = ('NOMC' 'UZ'
  170.         ('BRUI' 'BLAN' 'UNIF' 0.24 0.5  MAI1)
  171.         'NATU' 'DISCRET') ;
  172.  
  173.  GRAV  = GX 'ET' GY 'ET' GZ ;
  174.  
  175.  LISTINCO = 'MOTS' 'RN' 'RUX' 'RUY' 'RUZ' 'RET' ;
  176.  
  177.  RES  = ('NOMC' ('MOTS' 'UX' 'UY' 'UZ') (RN * GRAV)
  178.         ('MOTS' 'RUX' 'RUY' 'RUZ')
  179.          'NATU' 'DISCRET') 'ET'
  180.         ('NOMC' 'RET' ('PSCAL' GN GRAV ('MOTS' 'UX' 'UY' 'UZ')
  181.          ('MOTS' 'UX' 'UY' 'UZ')) 'NATU' 'DISCRET') ;
  182.  
  183.  RES1 = 'FIMP' 'VF' 'GRAVMONO' 'RESI' LISTINCO RN GN GRAV ;
  184.  
  185.  ERRES  = ('MAXIMUM' (RES '-' RES1) 'ABS') ;
  186.  
  187.  'SI' (ERRES '>' 1.0D-12);
  188.     'ERREUR'  5;
  189.  'FINSI' ;
  190.  
  191. * Jacobienne
  192.  
  193.  JAC = 'FIMP' 'VF' 'GRAVMONO' 'JACOCONS' LISTINCO RN GN GRAV ;
  194.  
  195.  CHPO1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RN' 1.0 '+'
  196.   ('MANUEL' 'CHPO' MAI1  1 'RN' 0.0) ;
  197.  
  198.  CHPRE = 'KOPS' 'MULT' JAC CHPO1 ;
  199.  
  200.  CHPRE1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 3 'RUX'
  201.      ('EXTRAIRE' GRAV A4 'UX') 'RUY'
  202.      ('EXTRAIRE' GRAV A4 'UY') 'RUZ'
  203.      ('EXTRAIRE' GRAV A4 'UZ')  ;
  204.  
  205.  ERRES  = ('MAXIMUM' (CHPRE1 '-' CHPRE) 'ABS') ;
  206.  
  207.  'SI' (ERRES '>' 1.0D-12);
  208.     'ERREUR'  5;
  209.  'FINSI' ;
  210.  
  211.  CHPO1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RUX' 1.0 '+'
  212.   ('MANUEL' 'CHPO' MAI1  1 'RUX' 0.0) ;
  213.  
  214.  CHPRE = 'KOPS' 'MULT' JAC CHPO1 ;
  215.  
  216.  CHPRE1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RET'
  217.      ('EXTRAIRE' GRAV A4 'UX') ;
  218.  
  219.  ERRES  = ('MAXIMUM' (CHPRE1 '-' CHPRE) 'ABS') ;
  220.  
  221.  'SI' (ERRES '>' 1.0D-12);
  222.     'ERREUR'  5;
  223.  'FINSI' ;
  224.  
  225.  CHPO1 = ('MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RUY' 1.0) '+'
  226.   ('MANUEL' 'CHPO' MAI1  1 'RUY' 0.0) ;
  227.  
  228.  CHPRE = 'KOPS' 'MULT' JAC CHPO1 ;
  229.  
  230.  CHPRE1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RET'
  231.      ('EXTRAIRE' GRAV A4 'UY') ;
  232.  
  233.  ERRES  = ('MAXIMUM' (CHPRE1 '-' CHPRE) 'ABS') ;
  234.  
  235.  'SI' (ERRES '>' 1.0D-12);
  236.     'ERREUR'  31 ;
  237.  'FINSI' ;
  238.  
  239.  CHPO1 = ('MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RUZ' 1.0) '+'
  240.   ('MANUEL' 'CHPO' MAI1  1 'RUZ' 0.0) ;
  241.  
  242.  CHPRE = 'KOPS' 'MULT' JAC CHPO1 ;
  243.  
  244.  CHPRE1 = 'MANUEL' 'CHPO' ('MANUEL' 'POI1' A4) 1 'RET'
  245.      ('EXTRAIRE' GRAV A4 'UZ') ;
  246.  
  247.  ERRES  = ('MAXIMUM' (CHPRE1 '-' CHPRE) 'ABS') ;
  248.  
  249.  'SI' (ERRES '>' 1.0D-12);
  250.     'ERREUR'  5 ;
  251.  'FINSI' ;
  252.  
  253.  'FIN' ;
  254.  
  255.  
  256.  

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