Procedures | Cast3M

   1 : * @MESU     PROCEDUR  AM        92/06/12    21:15:10     595            
   2 : *-------------------------------------------------                              
   3 : ******          PROCEDURE @MESU             ******
   4 : *-------------------------------------------------                              
   5 : *
   6 : *
   7 : *    CETTE PROCEDURE A ETE MISE GRACIEUSEMENT
   8 : *   A DISPOSOTION DE LA COMMUNAUTE  CASTEM2000
   9 : *    PAR M. LIBEYRE ( CEA/DSM/DRFC )
  10 : *
  11 : *     TEL : ( 33 1 ) 42 25 46 03
  12 : *
  13 : *-------------------------------------------------                              
  14 : *
  15 : * MESURE D'UN MAILLAGE (LONGUEUR, AIRE OU VOLUME)                               
  16 : *
  17 : *--------------------------------------------------------                       
  18 : DEBPROC @MESU GEO*MAILLAGE ;
  19 :  LTY = GEO ELEM TYPE ;                                                          
  20 :  ITCAR=FAUX;
  21 : *  mess 'liste des types d elements :' ;                                        
  22 : *  list LTY ;                                                                   
  23 :  NTY = DIME LTY ;                                                               
  24 :  I = 0 ;                                                                        
  25 :  E1 = FAUX ;                                                                    
  26 :  E2 = FAUX ;                                                                    
  27 :  E3 = FAUX ;                                                                    
  28 :  REPETER BOUC NTY ;                                                             
  29 :  I = I + 1 ;                                                                    
  30 :  ITY = EXTR LTY I ;                                                             
  31 :  EE1 = (EGA ITY SEG2) OU (EGA ITY SEG3) ;                                       
  32 :  EE2 = (EGA ITY TRI3) OU (EGA ITY TRI6) OU (EGA ITY QUA4)                       
  33 :       OU (EGA ITY QUA8) ;                                                       
  34 :  EE3 = (EGA ITY CU20) OU (EGA ITY PRI6) OU (EGA ITY PR15)                       
  35 :       OU (EGA ITY PYR5) OU (EGA ITY PY13) OU (EGA ITY TET4)                     
  36 :       OU (EGA ITY TE10) ;                                                       
  37 :  E1 = E1 OU EE1 ;                                                               
  38 :  E2 = E2 OU EE2 ;                                                               
  39 :  E3 = E3 OU EE3 ;                                                               
  40 :  FIN BOUC ;                                                                     
  41 : SI E1;                                                                          
  42 : * MESS 'MAILLAGE DE TYPE LIGNE' ;                                               
  43 : FINSI ;                                                                         
  44 : SI E2;                                                                          
  45 : * MESS 'MAILLAGE DE TYPE SURFACE' ;                                             
  46 : FINSI ;                                                                         
  47 : SI E3;                                                                          
  48 : * MESS 'MAILLAGE DE TYPE MASSIF' ;                                              
  49 : FINSI ;                                                                         
  50 : SI ((E1 ET E2) OU (E1 ET E3) OU (E2 ET E3)) ;                                   
  51 :  MESS 'LE MAILLAGE DOIT ETRE D UN TYPE UNIQUE : LIGNE,                          
  52 :        SURFACE OU MASSIF' ;                                                     
  53 :  P1 =0. ;                                                                       
  54 : SINON ;                                                                         
  55 :  SI E1 ;                                                                        
  56 :   IMA = 1 ;                                                                     
  57 :  FINSI ;                                                                        
  58 :  SI E2 ;                                                                        
  59 :   IMA = 2 ;                                                                     
  60 :  FINSI ;                                                                        
  61 :  SI E3 ;                                                                        
  62 :   IMA = 3 ;                                                                     
  63 :  FINSI ;                                                                        
  64 : * MESS 'DIMENSION DE GEO :';                                                    
  65 : * LIST IMA ;                                                                    
  66 :  SI (EGA IMA 1) ;                                                               
  67 : *  MESS 'CALCUL DE LA LONGUEUR' ;                                               
  68 :   GEO1 = CHANGER GEO SEG2 ;                                                     
  69 :   MODL1 = MODL GEO1 MECANIQUE  ELASTIQUE BARR ;                                 
  70 :   MAT1 = MATR MODL1 'RHO' 1. ;                                                  
  71 :   CAR1 = CARB MODL1 'SECT' 1.  ;                                                
  72 :   ITCAR=VRAI;
  73 :  FINSI ;                                                                        
  74 :  SI (EGA IMA 2) ;                                                               
  75 : *  MESS 'CALCUL DE L AIRE' ;                                                    
  76 :   I = 0 ;                                                                       
  77 :   REPETER BOUC2 NTY ;                                                           
  78 :   I = I + 1 ;                                                                   
  79 :   GEOT =  GEO ELEM (EXTR LTY I) ;                                               
  80 :   SI (EGA I 1) ;                                                                
  81 :    GEO1 = CHANGER GEOT TRI3 ;                                                   
  82 :   SINON ;                                                                       
  83 :    GEO1 = GEO1 ET (CHANGER GEOT TRI3) ;                                         
  84 :   FINSI ;                                                                       
  85 :   FIN BOUC2 ;                                                                   
  86 :   IDI = VALEUR DIME ;                                                           
  87 :   SI (EGA IDI 2) ;                                                              
  88 :   MODL1 = MODL GEO1 MECANIQUE  ELASTIQUE TRI3 ;                                 
  89 :   SINON ;                                                                       
  90 :   MODL1 = MODL GEO1 MECANIQUE  ELASTIQUE COQ3 ;
  91 :   CAR1 = CARB MODL1 'EPAI' 1.  ;                                                
  92 :   ITCAR=VRAI;
  93 :   FINSI ;                                                                       
  94 :   MAT1 = MATR MODL1 'RHO' 1. ;                                                  
  95 :  FINSI ;                                                                        
  96 :  SI (EGA IMA 3) ;                                                               
  97 : *  MESS 'CALCUL DU VOLUME' ;                                                    
  98 :   I = 0 ;                                                                       
  99 :   REPETER BOUC3 NTY ;                                                           
 100 :   I = I + 1 ;                                                                   
 101 :   GEOT =  GEO ELEM (EXTR LTY I) ;                                               
 102 :   SI (EGA I 1) ;                                                                
 103 :    GEO1 = CHANGER GEOT TET4 ;                                                   
 104 :   SINON ;                                                                       
 105 :    GEO1 = GEO1 ET (CHANGER GEOT TET4) ;                                         
 106 :   FINSI ;                                                                       
 107 :   FIN BOUC3 ;                                                                   
 108 :   MODL1 = MODL GEO1 MECANIQUE  ELASTIQUE TET4 ;                                 
 109 :   MAT1 = MATR MODL1 'RHO' 1. ;                                                  
 110 :  FINSI ;                                                                        
 111 : SI (IMA < 3) ;                                                                  
 112 :    SI ITCAR ;
 113 :      MASS1 = MASSE MODL1 (MAT1 ET CAR1) ;
 114 :    SINON ;
 115 :     MASS1 = MASSE MODL1 MAT1  ;
 116 :    FINSI ;
 117 : SINON ;                                                                         
 118 :  MASS1 = MASSE MODL1 MAT1  ;                                                    
 119 : FINSI ;                                                                         
 120 : GEO2 = CHANGER GEO POI1 ;                                                       
 121 : GAM1 = MANU CHPO GEO2 1 UX 1. ;                                                 
 122 : POIP = MASS1 * GAM1 ;                                                           
 123 : P1 = @TOTAL POIP GEO2 FX ;
 124 : FINSI;
 125 : FINPROC P1 ;