Télécharger la plaquette Cast3M


fr en

Télécharger fornod.procedur

Retour à la liste

Numérotation des lignes : 

  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.  
  84.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales