Procedures | Cast3M

   1 : * @COUPLE   PROCEDUR  CHAT      96/01/03    21:17:53     1971           
   2 : ************************************************************************        
   3 : *                                                                               
   4 : *   PROCEDURE  @COUPLE                  MLN/CERN JANVIER 91                     
   5 : *                                                                               
   6 : *                                                                               
   7 : *   LISSAGE D'UNE LISSE PAR LA METHODE D'ENERGIE MINIMUM DE                     
   8 : *   DEFORMATION                                                                 
   9 : *                                                                               
  10 : *-----------------------------------------------------------------------        
  11 : *                                                                               
  12 : *   ENTREE                                                                      
  13 : *                                                                               
  14 : *                                                                               
  15 : *    EVO_1            Evolution dans le plan YOZ                                
  16 : *                                                                               
  17 : *    EVO_2            Evolution dans le plan XOZ                                
  18 : *                                                                               
  19 : *    ENT_1            Entier = 1 ,le couple est dans le plan XOZ                
  20 : *                              2 ,le couple est dans le plan YOZ                
  21 : *                                                                               
  22 : *                                                                               
  23 : *-----------------------------------------------------------------------        
  24 : *                                                                               
  25 : *   SORTIE                                                                      
  26 : *                                                                               
  27 : *    COU_1            Couple 3D                                                 
  28 : *                                                                               
  29 : *                                                                               
  30 : *-----------------------------------------------------------------------        
  31 : 'DEBP' @COUPLE  EVO_1*'EVOLUTION' EVO_2*'EVOLUTION' ENT_1*'ENTIER' ;            
  32 : 'TITRE' ' COUPLE' ;                                                             
  33 : 'OPTI' 'DIME' 3 'ELEM' 'SEG2' 'MODE' 'TRIDIM' ;                                 
  34 : 'SI' ('EGA' ENT_1 1) ;                                                          
  35 : EVO_1X = 'EXTR' EVO_2 'ABSC' ;                                                  
  36 : EVO_1Y = 'EXTR' EVO_1 'ABSC' ;                                                  
  37 : EVO_1Z = 'EXTR' EVO_2 'ORDO' ;                                                  
  38 : 'FINSI' ;                                                                       
  39 : 'SI' ('EGA' ENT_1 2) ;                                                          
  40 : EVO_1X = 'EXTR' EVO_2 'ABSC' ;                                                  
  41 : EVO_1Y = 'EXTR' EVO_1 'ABSC' ;                                                  
  42 : EVO_1Z = 'EXTR' EVO_1 'ORDO' ;                                                  
  43 : 'FINSI' ;                                                                       
  44 : NBP = 'DIME' EVO_1X ;                                                           
  45 : I = 0 ;                                                                         
  46 : PX = 'TABLE' ;                                                                  
  47 : PY = 'TABLE' ;                                                                  
  48 : PZ = 'TABLE' ;                                                                  
  49 : PP = 'TABLE' ;                                                                  
  50 : 'REPE' BOU_1 NBP ;                                                              
  51 : I =  I + 1 ;                                                                    
  52 : PX.I =  ('EXTR' EVO_1X I) ;                                                     
  53 : PY.I =  ('EXTR' EVO_1Y I) ;                                                     
  54 : PZ.I =  ('EXTR' EVO_1Z I) ;                                                     
  55 : PP.I = (PX.I PY.I PZ.I) ;                                                       
  56 : 'SI' ('EGA' 2 I) ;                                                              
  57 :     MAI_1 = (PP.(I - 1) 'DROI' 1 PP.I) ;                                        
  58 : 'FINSI' ;                                                                       
  59 : 'SI' (I >  2) ;                                                                 
  60 :     MAI_1 = MAI_1 'ET' (PP.(I - 1) 'DROI' 1 PP.I) ;                             
  61 : 'FINSI' ;                                                                       
  62 : 'FIN' BOU_1 ;                                                                   
  63 : CON_1 = 'MAXI' ('ABS' EVO_1Z) ;                                                 
  64 : CON_1 = CON_1 / 10. ;                                                           
  65 : O = 0. 0. 0. ;                                                                  
  66 : O1 = 1E5 1E5 1E5 ;                                                              
  67 : X = O 'DROI' 1 (CON_1 0. 0.) ;                                                  
  68 : Y = O 'DROI' 1 (0. CON_1 0.) ;                                                  
  69 : Z = O 'DROI' 1 (0. 0. CON_1) ;                                                  
  70 : X = X 'COUL' 'VERT' ;                                                           
  71 : Y = Y 'COUL' 'VERT' ;                                                           
  72 : Z = Z 'COUL' 'VERT' ;                                                           
  73 : MAI_1 = MAI_1 'COUL' 'ROUG' ;                                                   
  74 : 'TRAC' O1 'QUAL' (X ET Y ET Z ET MAI_1) ;                                       
  75 : *                                                                               
  76 : 'FINPROC' MAI_1 ;       
  77 :