Procedures | Cast3M

   1 : * FORNOD    PROCEDUR  WANG      96/09/17    21:16:05     2278           
   2 :                                                                                 
   3 : 'DEBPROC' FORNOD LIG1*'MAILLAGE' COM1*'MOT' FO*'FLOTTANT' ;                     
   4 :                                                                                 
   5 : *---------------------------------------------------------------------          
   6 : *                                     
   7 : *     POUR CALCULER LES FORCES NODALES EQUIVALENTES D'UNE FORCE
   8 : *     SURFACIQUE  Fo 
   9 : *     ELEMENT HORIZONTAL DE TYPE SEG3 
  10 : *     MODE DE CALCUL : 'PLANDEFO' 'AXIS' ET 'FOUR 0, 1'
  11 : *                                                                               
  12 : *     LIG1     'MAIILAGE'     ligne horizontale sur laquelle on                 
  13 : *                             applique la force surfacique                      
  14 : *     COM      'MOT'          nom de la composante de la force                  
  15 : *                             PLAN DEFO  : FX, FY    
  16 : *                             FOUR 1     : FR ,FZ ,FT     
  17 : *                                                                               
  18 : *     FO       'FLOTTANT'     valeur de la composante de la force               
  19 : *                             surfacique                                        
  20 : *                                                                               
  21 : *---------------------------------------------------------------------          
  22 :                                                                                 
  23 : TYP_MODE  = 'MOT' ('VALEUR' 'MODE') ;
  24 : 'SI'   ('EGA' TYP_MODE 'PLANDEFO') ;                                         
  25 :   ITYP = 1 ;                                                             
  26 : 'SINON' ; 'SI' ('EGA' TYP_MODE 'FOUR') ;                                
  27 :   ITYP = 2 ;                                                                    
  28 : 'SINON' ;                                                                       
  29 :   'MESSAGE' '           *****  ERREUR  *****' ;                                 
  30 :   'MESSAGE' '   On ne traite pas les modèles tridimensionnelles.'
  31 :   'MESSAGE' '   Arrêt de la procédure.' ;
  32 : 'QUITTER' BOUCERR ;                                                             
  33 : 'FINSI' ; 'FINSI' ; 'FINSI' ;                                                   
  34 : *                                                                               
  35 : *   Boucle sur les éléments
  36 : *                                                                               
  37 : NB_ELEM = 'NBEL' LIG1 ;                                                         
  38 : I_EL = 0 ;                                                                      
  39 : 'REPETER' BOUCK NB_ELEM            ;                                            
  40 :   I_EL = I_EL + 1                  ;                                            
  41 :   ELEM_I = LIG1 'ELEMENT' I_EL     ;                                            
  42 :   ELEM_B = 'CHANGER' ELEM_I 'POI1' ;                                            
  43 :   NOEU1 = ELEM_B 'POINT' 1 ;                                                    
  44 :   NOEU2 = ELEM_B 'POINT' 2 ;                                                    
  45 :   NOEU3 = ELEM_B 'POINT' 3 ;                                                    
  46 :   R1 = 'COOR' NOEU1 1      ;                                    
  47 :   R3 = 'COOR' NOEU3 1      ;                      
  48 :   LONG = 'ABS' (R3 - R1)   ;                                  
  49 : 
  50 : ***   MODE DEFO PLAN    ***        
  51 : 
  52 : 'SI' ('EGA' ITYP 1)      ;                                                      
  53 : F1 = LONG * FO / 6.      ;                                                      
  54 : F2 = F1 * 4.             ;                                                      
  55 : F3 = F1                  ;                                                      
  56 :                                                          
  57 : 'SI' ('EGA' I_EL 1) ;                                                           
  58 : FOR1 = ( 'FORCE' COM1 F1 NOEU1 ) 'ET' ( 'FORCE' COM1 F2 NOEU2 )                 
  59 :   'ET' ( 'FORCE' COM1 F3 NOEU3 ) ;                                              
  60 : 'SINON' ; FOR1 = FOR1 'ET' ( 'FORCE' COM1 F1 NOEU1 ) 'ET'                       
  61 :          ( 'FORCE' COM1 F2 NOEU2 ) ET ( 'FORCE' COM1 F3 NOEU3 ) ;               
  62 : 'FINSI' ;                                                                       
  63 : 'FINSI' ;                                                                       
  64 :                                                                                 
  65 : ***   MODE AXIS ET FOURIER   ***                                            
  66 :                                                                                 
  67 : 'SI' ( 'EGA' ITYP 2 ) ;                                                         
  68 :  F1 = PI * LONG * FO *   R1 / 6. ;                                              
  69 :  F2 = PI * LONG * FO * ( R1 + R3 ) / 3. ;                                       
  70 :  F3 = PI * LONG * FO *   R3 / 6. ;                                              
  71 : 'SI' ('EGA' I_EL 1) ;                                                           
  72 :  FOR1 = ( 'FORCE' COM1 F1 NOEU1 ) 'ET' ( 'FORCE' COM1 F2 NOEU2 )                
  73 :        'ET' ( 'FORCE' COM1 F3 NOEU3 ) ;                                         
  74 : 'SINON' ;                                                                       
  75 :  FOR1 = FOR1 'ET' ( 'FORCE' COM1 F1 NOEU1 ) 'ET'                                
  76 :           ( 'FORCE' COM1 F2 NOEU2 ) 'ET' ( 'FORCE' COM1 F3 NOEU3 ) ;            
  77 : 'FINSI' ;                                                                       
  78 : 'FINSI' ;                                                                       
  79 : 'FIN' BOUCK ;                                                                   
  80 : 'FINPROC' FOR1  ;                                                               
  81 : ***********************************************************************         
  82 : 
  83 :