Télécharger la plaquette Cast3M


fr en

Télécharger pret_scal1.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  pret_scal1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **********************************************************************
  5. **** APPROCHE VF                                                  ****
  6. ****                                                              ****
  7. **** OPERATEUR PRET                                               ****
  8. **** Operateur qui reconstruit les variables primitives aux faces ****
  9. ****                                                              ****
  10. **** A. BECCANTINI SFME/LTMF                       NOVEMBRE 01    ****
  11. **********************************************************************
  12.  
  13.  'OPTION'  'DIME' 2  'ELEM' QUA4  'ECHO'  0  'TRAC' 'X' ;
  14.  
  15. *
  16. *** GRAPH
  17. *
  18.  
  19. * GRAPH = VRAI ;
  20.  GRAPH = FAUX ;
  21.  
  22.  
  23. *************************** 
  24. ***** DOMAINE SPATIAL  ****
  25. ***************************
  26.  
  27. *
  28. **** Deux carre
  29. *
  30.  
  31.  A1 = 0.0D0 0.0D0;
  32.  A2 = 1.0D0 0.0D0;
  33.  A3 = 2.0D0 0.0D0;
  34.  A4 = 2.0D0 1.0D0;
  35.  A5 = 1.0D0 1.0D0;
  36.  A6 = 0.0D0 1.0D0;
  37.  
  38.  L12 = A1 'DROIT' 1 A2;
  39.  L23 = A2 'DROIT' 1 A3;
  40.  L34 = A3 'DROIT' 1 A4;
  41.  L45 = A4 'DROIT' 1 A5;
  42.  L56 = A5 'DROIT' 1 A6;
  43.  L61 = A6 'DROIT' 1 A1;
  44.  L25 = A2 'DROIT' 1 A5;
  45.  
  46.  
  47.  DOM10  = 'DALL' L12 L25 L56 L61
  48.           'PLANE';
  49.  DOM20  = 'DALL' L23 L34 L45 ('INVERSE' L25)
  50.           'PLANE';
  51.  
  52. *
  53. *** Point face entre le deux carre, ou on fait les controlles
  54. *
  55.  
  56.  P10 = 1.0 0.5;
  57.  
  58.  
  59.  DOM1 =  DOM10 ;
  60.  DOM2 =  DOM20 ;
  61.  
  62.  
  63. 'ELIMINATION' (DOM1 ET DOM2) 1D-6;
  64. DOMTOT = DOM1 ET DOM2;
  65.  
  66.  $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  67.  
  68.  $DOM1 = 'MODELISER'  DOM1 'EULER';
  69.  $DOM2 = 'MODELISER'  DOM2 'EULER';
  70.  
  71.  TDOMTOT = 'DOMA'  $DOMTOT 'VF';
  72.  
  73.  TDOM1 = 'DOMA'  $DOM1 'VF';
  74.  TDOM2 = 'DOMA'  $DOM2 'VF';
  75.  
  76.  MDOM1 = TDOM1 . 'QUAF' ;
  77.  MDOM2 = TDOM2 . 'QUAF' ;
  78.  MDOMTOT = TDOMTOT . 'QUAF' ;
  79.  
  80.  'ELIMINATION'  (MDOMTOT ET MDOM1) 0.0001 ;
  81.  'ELIMINATION'  (MDOMTOT ET MDOM2) 0.0001 ;
  82.  
  83.  P1 = ('DOMA' $DOMTOT 'FACE') 'POIN' 'PROC' P10 ;
  84.  
  85.  'SI' GRAPH;
  86.    'TRACER' (('DOMA' $DOMTOT 'MAILLAGE') 'ET' ('DOMA' $DOMTOT 'FACEL')
  87.    'ET' P1) 'TITRE' 'Domaine et FACEL';
  88.  'FINSI' ;
  89.  
  90.  ETATG = 'PROG' 1.0 2.0 3.0 4.0 ;
  91.  ETATD = ETATG '*' 10 ;
  92.  
  93.  SN =('MANUEL' 'CHPO' ('DOMA' $DOM1 'CENTRE') 4 'C1' 1.0 'C2' 2.0 
  94.      'C3' 3.0 'C4' 4.0 'NATU' 'DISCRET') 'ET'
  95.      ('MANUEL' 'CHPO' ('DOMA' $DOM2 'CENTRE') 4 'C1' 10.0 'C2' 20.0 
  96.      'C3' 30.0 'C4' 40.0 'NATU' 'DISCRET') ;
  97.  
  98.  SF =  'PRET' 'CLAUDEIS' 'FACE' 1  $DOMTOT SN ;
  99.  
  100.  
  101. *********************************************************
  102. *** Control des etats sur la surface qui contient P1 ****
  103. *********************************************************
  104.  
  105.  
  106.  GEOP1    = ('DOMA' $DOMTOT 'FACEL') 'ELEM' 'APPUYE' 'LARGEMENT' P1;
  107.  GEOP2    = ('DOMA' $DOMTOT 'FACE')  'ELEM' 'APPUYE' 'LARGEMENT' P1;
  108.  
  109.  
  110.  
  111.  SGEOP1  = 'REDU' SF GEOP1;
  112.  
  113.  C1g = 'EXTRAIRE' SGEOP1 'C1' 1 1 1;
  114.  C1d = 'EXTRAIRE' SGEOP1 'C1' 1 1 3;
  115.  
  116.  C2g = 'EXTRAIRE' SGEOP1 'C2' 1 1 1;
  117.  C2d = 'EXTRAIRE' SGEOP1 'C2' 1 1 3;
  118.  
  119.  C3g = 'EXTRAIRE' SGEOP1 'C3' 1 1 1;
  120.  C3d = 'EXTRAIRE' SGEOP1 'C3' 1 1 3;
  121.  
  122.  C4g = 'EXTRAIRE' SGEOP1 'C4' 1 1 1;
  123.  C4d = 'EXTRAIRE' SGEOP1 'C4' 1 1 3;
  124.  
  125.  CHPNOR = 'DOMA' $DOMTOT 'NORMALE' ;
  126.  
  127. *
  128. **** Orientation de la normal n de castem par raport a la
  129. *    notre; t est par consequence
  130. *
  131.  
  132.  ORIENT = 'EXTRAIRE' CHPNOR P1 'UX' ;
  133.  
  134. *
  135. **** ORIENT = -1 -> Mon etat gauche est son etat droite
  136. *
  137.  
  138.  'SI' (ORIENT > 0);
  139.  
  140.     ERRLIG = 'PROG' C1g C2g C3g C4g;
  141.  
  142.     ERRLID = 'PROG' C1d C2d C3d C4d; 
  143.  
  144.  'SINON' ;                
  145.  
  146.     ERRLID = 'PROG' C1g C2g C3g C4g;
  147.  
  148.     ERRLIG = 'PROG' C1d C2d C3d C4d;
  149.  
  150.  'FINSI' ;
  151.  
  152.  ERRO = 'MAXIMUM' ('PROG' 
  153.    ('MAXIMUM' (ETATG '-' ERRLIG) 'ABS')
  154.    ('MAXIMUM' (ETATD '-' ERRLID) 'ABS')
  155.  );
  156.  
  157.  
  158.  'SI' (ERRO > 1.0D-14)
  159.     'MESSAGE' 'Ordre en espace = 1';
  160.     'MESSAGE' 'PRET = ???'
  161.     'ERREUR' 5 ;
  162.  'FINSI' ;
  163.  
  164. 'FIN' ;
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  

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