Procedures | Cast3M

   1 : * @TOLE3    PROCEDUR  CHAT      96/01/03    21:31:56     1971           
   2 : 'DEBP'  @TOLE3   A_1*'MAILLAGE' A_2*'MAILLAGE' RA_1*'FLOTTANT'                  
   3 : N_1*'ENTIER' ;                                                                  
   4 : ME = 'TEXT' 'MECANIQUE ELASTIQUE' ;                                             
   5 : 'OPTI' 'ELEM' 'CUB8' ;                                                          
   6 : S_1 = A_1 'REGL' 1 A_2 ;                                                        
   7 :     'OPTI' 'ELEM' 'PRI6' 'MODE' 'TRIDIM' ;                                      
   8 :     AF_1 = 'MODL' S_1 ME 'COQ4' ;                                               
   9 :     VE_1 = 'PRES' 'COQU' AF_1 1. 'NORM' ;                                       
  10 :     VE_1 = 'REDU' VE_1 S_1 ;                                                    
  11 :     V1NX = 'EXCO' 'FX' VE_1 'SCAL' ;                                            
  12 :     V1NY = 'EXCO' 'FY' VE_1 'SCAL' ;                                            
  13 :     V1NZ = 'EXCO' 'FZ' VE_1 'SCAL' ;                                            
  14 :     N1VN = ((V1NX ** 2) + (V1NY ** 2) + (V1NZ ** 2)) ** .5 ;                    
  15 :     V1UNX = V1NX / N1VN ;                                                       
  16 :     V1UNY = V1NY / N1VN ;                                                       
  17 :     V1UNZ = V1NZ / N1VN ;                                                       
  18 :     R1X = V1UNX * RA_1 ;                                                        
  19 :     R1Y = V1UNY * RA_1 ;                                                        
  20 :     R1Z = V1UNZ * RA_1 ;                                                        
  21 :     TR_1 = ('NOMC' UX R1X) + ('NOMC' UY R1Y) + ('NOMC' UZ R1Z) ;                
  22 :     L_1 = (A_1 'PLUS' TR_1) ;                                                   
  23 :     L_2 = (A_2 'PLUS' TR_1) ;                                                   
  24 : NN1 =  'NBNO' A_1 ;                                                             
  25 : NN2 = NN1 / 2 ;                                                                 
  26 : PA_11 = A_1 'POIN' 1 ;                                                          
  27 : PA_12 = A_1 'POIN' NN2 ;                                                        
  28 : PA_13 = A_1 'POIN' NN1 ;                                                        
  29 : PA_21 = A_2 'POIN' 1 ;                                                          
  30 : PA_22 = A_2 'POIN' NN2 ;                                                        
  31 : PA_23 = A_2 'POIN' NN1 ;                                                        
  32 : PA_11X = PA_11 COOR 1 ;                                                         
  33 : PA_11Y = PA_11 COOR 2 ;                                                         
  34 : PA_31 =  PA_11X PA_11Y 1. ;                                                     
  35 : P_1 = 'TABL' ;                                                                  
  36 : P_2 = 'TABL' ;                                                                  
  37 : *                                                                               
  38 : I = 1 ;                                                                         
  39 :   'REPE' BOU100 NN1 ;                                                           
  40 :   PP1 = L_1 'POIN' I ;                                                          
  41 :   PP2 = L_2 'POIN' I ;                                                          
  42 :   DP1 = PP2 'MOIN' PP1 ;                                                        
  43 :   P_1.I = PP1 'PROJ' 'CYLI' DP1 'PLAN' PA_11 PA_12 PA_13 ;                      
  44 :   P_2.I = PP2 'PROJ' 'CYLI' DP1 'PLAN' PA_21 PA_22 PA_23 ;                      
  45 :   I = I + 1 ;                                                                   
  46 :   'FIN' BOU100 ;                                                                
  47 : *                                                                               
  48 : I = 1 ;                                                                         
  49 :   'REPE' BOU101 NN1 ;                                                           
  50 :   P_1IY = P_1.I 'COOR' 2 ;                                                      
  51 :     'SI' (P_1IY < PA_11Y) ;                                                     
  52 :     DP_1 = P_1.I 'MOIN' P_1.(I + 1) ;                                           
  53 :     P_1.I = P_1.I 'PROJ' 'CYLI' DP_1 'PLAN' PA_31 PA_11 PA_21 ;                 
  54 :     P_1.I = P_1.I 'PLUS' ((0. 0. 1.E-5) * (I - 1)) ;                            
  55 :     'SINON' ;                                                                   
  56 :     'QUIT' BOU101 ;                                                             
  57 :     'FINSI';                                                                    
  58 :   I = I + 1 ;                                                                   
  59 :   'FIN' BOU101 ;                                                                
  60 : *                                                                               
  61 : I = 1 ;                                                                         
  62 :   'REPE' BOU103 NN1 ;                                                           
  63 :   P_2IY = P_2.I 'COOR' 2 ;                                                      
  64 :     'SI' (P_2IY < PA_11Y) ;                                                     
  65 :     DP_1 = P_2.I 'MOIN' P_2.(I + 1) ;                                           
  66 :     P_2.I = P_2.I 'PROJ' 'CYLI' DP_1 'PLAN' PA_31 PA_11 PA_21 ;                 
  67 :     P_2.I = P_2.I 'PLUS' ((0. 0. 1.E-5) * (I - 1)) ;                            
  68 :     'SINON' ;                                                                   
  69 :     'QUIT' BOU103 ;                                                             
  70 :     'FINSI' ;                                                                   
  71 :   I = I + 1 ;                                                                   
  72 :   'FIN' BOU103 ;                                                                
  73 : *                                                                               
  74 : I = 2 ;                                                                         
  75 :   'REPE' BOU102 (NN1 - 1) ;                                                     
  76 :   'SI' ('EGA' I 2) ;                                                            
  77 :   C_1 = P_1.(I - 1) 'DROI' 1 P_1.I ;                                            
  78 :   C_2 = P_2.(I - 1) 'DROI' 1 P_2.I ;                                            
  79 :   'FINSI' ;                                                                     
  80 :   'SI' (I > 2) ;                                                                
  81 :   C_1 = C_1 'ET' (P_1.(I - 1) 'DROI' 1 P_1.I) ;                                 
  82 :   C_2 = C_2 'ET' (P_2.(I - 1) 'DROI' 1 P_2.I) ;                                 
  83 :   'FINSI' ;                                                                     
  84 :   I = I + 1 ;                                                                   
  85 :   'FIN' BOU102 ;                                                                
  86 : A_1 = A_1 'COUL' 'ROUG' ;                                                       
  87 : A_2 = A_2 'COUL' 'ROUG' ;                                                       
  88 : C_1 = C_1 'COUL' 'ROUG' ;                                                       
  89 : C_2 = C_2 'COUL' 'ROUG' ;                                                       
  90 : S_2 = C_1 'REGL' 1 C_2 ;                                                        
  91 : V_1 = S_1 'VOLU' N_1 S_2 ;                                                      
  92 : V_1 = V_1 'COUL' 'ROUG' ;                                                       
  93 :  'DETR' VE_1 ; 'DETR' V1NX ; 'DETR' V1NY ; 'DETR' V1NZ ;                        
  94 :  'DETR' V1UNX ; 'DETR' V1UNY ; 'DETR' V1UNZ ;                                   
  95 :  'DETR' R1X ; 'DETR' R1Y ; 'DETR' R1Z ; 'DETR' TR_1 ;                           
  96 : 'FINP'  C_1 C_2 S_1 S_2 V_1 ;  
  97 :