Télécharger la plaquette Cast3M


fr en

Télécharger flamcrebcom.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  flamcrebcom.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **********************************************************************
  5. **** OPERATEUR FLAM                                               ****
  6. **** Critere CREBCOM                                              ****
  7. ****                                                              ****
  8. **** A. BECCANTINI DEN/DM2S/SFME/LTMF        JUIN 2001            ****
  9. **********************************************************************
  10.  
  11.  'OPTION'  'DIME' 2 ;
  12.  'OPTION'  'ELEM' QUA4 ;
  13.  'OPTION'  'TRAC' 'X' ;
  14.  'OPTION'  'ECHO'  0  ;
  15.  
  16.  
  17.  GRAPH = FAUX ;
  18.  
  19.  
  20.  LIG1 = (0.0 0.0) 'DROIT' 10 (0.0 1.0) ;
  21.  DOM1 = 'TRANSLATION' LIG1 10 (1.0 0.0) ;
  22.  DOM2 = 'TRANSLATION' LIG1 10 (-1.0 0.0) ;
  23.  DOMTOT = DOM1 'ET' DOM2 ;
  24.  
  25.  $DOMTOT = 'MODE' DOMTOT 'EULER'  ;
  26.  $DOM1   = 'MODE' DOM1   'EULER'  ;
  27.  $DOM2   = 'MODE' DOM2   'EULER'  ;
  28.  
  29.  TDOMTOT = 'DOMA' $DOMTOT 'VF'  ;
  30.  TDOM1   = 'DOMA' $DOM1   'VF'  ;
  31.  TDOM2   = 'DOMA' $DOM2   'VF'  ;
  32.  
  33.  MDOMTOT = TDOMTOT . 'QUAF';
  34.  MDOM1   = TDOM1   . 'QUAF';
  35.  MDOM2   = TDOM2   . 'QUAF';
  36. 'ELIM' (MDOMTOT 'ET' MDOM1 'ET' MDOM2) 1D-5;
  37.  
  38.  EPSILON = 1.0D-2 ;
  39.  CSIMAX  = 9.0D-1 ;
  40.  
  41.  
  42. *
  43. **** CSIMAX: maximum value of CSI
  44. *
  45. *    CSIN = CSIMAX if x   <  0
  46. *           0.0   x  =>  0
  47. *
  48. *
  49.  
  50.  
  51.  CSIN = ('MANUEL' 'CHPO' ('DOMA' $DOM2 'CENTRE') 1 'H2O2'
  52.          (1.0001 '*' CSIMAX) 'NATURE' 'DISCRET')
  53.         'ET'
  54.         ('MANUEL' 'CHPO' ('DOMA' $DOM1 'CENTRE') 1 'H2O2' 0.0
  55.          'NATURE' 'DISCRET') ;
  56.  
  57.  
  58.  MOD1  =  'MODELISER'  ('DOMA' $DOMTOT 'MAILLAGE') 'THERMIQUE' ;
  59.  
  60.  
  61.  'SI' GRAPH ;
  62.     CHM_CSI  =  'KCHA' $DOMTOT 'CHAM' CSIN ;
  63.     'TRAC' CHM_CSI  MOD1 'TITR' ('CHAINE' 'csi');
  64.  'FINSI' ;
  65.  
  66.  
  67.  CHP1 = 'FLAM' 'CREBCOM' $DOMTOT EPSILON CSIMAX CSIN ;
  68.  
  69.  
  70. *
  71. **** CHP1 =  1.0 in x = 0.05
  72. *            0   ailleurs
  73. *
  74.  
  75.  'SI' GRAPH ;
  76.     CHM_C1   =  'KCHA' $DOMTOT 'CHAM' CHP1 ;
  77.     'TRAC' CHM_C1   MOD1 'TITR' ('CHAINE' 'C1');
  78.  'FINSI' ;
  79.  
  80.  GEO1 = ('COORDONNEE' 1 ('DOMA' $DOMTOT 'CENTRE')) 'POIN' 'COMPRIS'
  81.    0.049  0.051 ;
  82.  CHPTEST = 'MANUEL' 'CHPO' GEO1 1 'SCAL' 1.0 ;
  83.  
  84.  ERRO = 'MAXIMUM' (CHP1 '-' CHPTEST) 'ABS' ;
  85.  
  86.  'SI' (ERRO > 1.0D-6) ;
  87.      'ERREUR' 5 ;
  88.  'FINSI' ;
  89.  
  90.  'FIN' ;
  91.  
  92.  
  93.  
  94.  
  95.  
  96. *
  97.  
  98.  
  99.  
  100.  
  101.  

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