Procedures | Cast3M

   1 : * @TOLE2    PROCEDUR  CHAT      96/01/03    21:31:52     1971           
   2 : **********************************************************************          
   3 : 'DEBP' @TOLE2  L_1*'MAILLAGE' RA_1*'FLOTTANT' N_1*'ENTIER' ;                    
   4 : EL_1 = L_1 'ELEM' 'TYPE' ;                                                      
   5 : DEL_1 = 'DIME' EL_1 ;                                                           
   6 : DI_1 = 'VALE' 'DIME' ;                                                          
   7 : MS = 'TEXT' 'MECANIQUE ELASTIQUE' ;                                             
   8 : 'SI' ('EGA' DI_1  2 ) ;                                                         
   9 :     OP_1 = 'MOT' 'REGL' ;                                                       
  10 :     'OPTI' 'ELEM' 'QUA4' ;                                                      
  11 :     AF_1 = 'MODL' L_1 MS 'COQ2' ;                                               
  12 :     VE_1 = 'PRES' 'COQU' AF_1 1. 'NORM' ;                                       
  13 :     VE_1 = 'REDU' VE_1 L_1 ;                                                    
  14 :     V1NX = 'EXCO' 'FX' VE_1 'SCAL' ;                                            
  15 :     V1NY = 'EXCO' 'FY' VE_1 'SCAL' ;                                            
  16 :     N1VN = ((V1NX ** 2) + (V1NY ** 2)) ** .5 ;                                  
  17 :     V1UNX = V1NX / N1VN ;                                                       
  18 :     V1UNY = V1NY / N1VN ;                                                       
  19 :     RA_1 = RA_1 / 2. ;                                                          
  20 :     R1X = V1UNX * RA_1 ;                                                        
  21 :     R1Y = V1UNY * RA_1 ;                                                        
  22 :     R2X = V1UNX * (RA_1 * -1.) ;                                                
  23 :     R2Y = V1UNY * (RA_1 * -1.) ;                                                
  24 :     TR_1 = ('NOMC' UX R1X) + ('NOMC' UY R1Y) ;                                  
  25 :     TR_2 = ('NOMC' UX R2X) + ('NOMC' UY R2Y) ;                                  
  26 :     L_3 = (L_1 'PLUS' TR_1) 'COUL' 'ROUG'  ;                                    
  27 :     L_4 = (L_1 PLUS TR_2) COUL ROUG  ;                                          
  28 :     RE_1 = L_3 OP_1 N_1 L_4 ;                                                   
  29 : 'SINON' ;                                                                       
  30 :   OP_1 = 'MOT' 'VOLU' ;                                                         
  31 :   'SI' ('EGA' DEL_1 1) ;                                                        
  32 :     MT1 = 'EXTR' EL_1 1 ;                                                       
  33 :     'SI' ('EGA' MT1 'QUA4') ;                                                   
  34 :       'OPTI' 'DIME' 3 'ELEM' 'CUB8' 'MODE' 'TRIDIM' ;                           
  35 :       AF_1 = 'MODL' L_1 MS 'COQ4' ;                                             
  36 :     'FINSI' ;                                                                   
  37 :     'SI' ('EGA' MT1 'TRI3') ;                                                   
  38 :       'OPTI' 'DIME' 3 'ELEM' 'PRI6' 'MODE' 'TRIDIM' ;                           
  39 :       AF_1 = 'MODL' L_1 MS 'COQ3' ;                                             
  40 :     'FINSI' ;                                                                   
  41 :   'SINON' ;                                                                     
  42 :       'OPTI' 'ELEM' 'CUB8' ;                                                    
  43 :       AF_1 = 'MODL' L_1 MS 'COQ4' 'COQ3' ;                                      
  44 :   'FINSI' ;                                                                     
  45 :     VE_1 = 'PRES' 'COQU' AF_1 1. 'NORM' ;                                       
  46 :     VE_1 = 'REDU' VE_1 L_1 ;                                                    
  47 :     V1NX = 'EXCO' 'FX' VE_1 'SCAL' ;                                            
  48 :     V1NY = 'EXCO' 'FY' VE_1 'SCAL' ;                                            
  49 :     V1NZ = 'EXCO' 'FZ' VE_1 'SCAL' ;                                            
  50 :     N1VN = ((V1NX ** 2) + (V1NY ** 2) + (V1NZ ** 2)) ** .5 ;                    
  51 :     V1UNX = V1NX / N1VN ;                                                       
  52 :     V1UNY = V1NY / N1VN ;                                                       
  53 :     V1UNZ = V1NZ / N1VN ;                                                       
  54 :     R1X = V1UNX * RA_1 ;                                                        
  55 :     R1Y = V1UNY * RA_1 ;                                                        
  56 :     R1Z = V1UNZ * RA_1 ;                                                        
  57 :     TR_1 = ('NOMC' UX R1X) + ('NOMC' UY R1Y) + ('NOMC' UZ R1Z) ;                
  58 :     L_2 = (L_1 'PLUS' TR_1) 'COUL' 'ROUG'  ;                                    
  59 :     RE_1 = L_1 OP_1 N_1 L_2 ;                                                   
  60 : 'FINSI' ;                                                                       
  61 : 'FINPROC' RE_1 ;  
  62 :