Télécharger la plaquette Cast3M


fr en

Télécharger forbloc.procedur

Retour à la liste

Numérotation des lignes : 

  1. * FORBLOC   PROCEDUR  FANDEUR   14/10/10    21:15:15     8178           
  2. 'DEBPROC'  FORBLOC  BLOC*'MAILLAGE' SOL1*'CHPOINT' COUR*'FLOTTANT'  
  3.                                             AXI/'LOGIQUE'  ;            
  4. ************************************************************************
  5. *      MAGNETOSTATIQUE  2D POTENTIEL VECTEUR                            
  6. *   BLOC  maillage  non complexe                                        
  7. *   SOL1   SOLTUTION EN POTENTIEL  DEFINIE AU MOINS SUR BLOC            
  8. *  integration de j vectoriel b  sur les elements                       
  9. *   RPT     champ par point ( FX FY ) aux noeuds                        
  10. *   RFB     resultante  des forces au barycentre de la GEOMETRIE        
  11. *   AXI     logique valant vrai si pb axisym ( defaut  faux pb plan) 
  12. *     en plan  forces/unite de longueur
  13. *     en axisymetrique  forces/radian attention signification 
  14. ************************************************************************
  15.  'SI'( 'NON' ('EXISTE' AXI) ) ; AXI = FAUX ; 'FINSI' ;                  
  16.  BT = INDUCTIO BLOC  SOL1 AXI  ;                                        
  17.  OBJ0 = BLOC  'MODE' 'THERMIQUE' 'ISOTROPE' ;                           
  18.  BBY = 'EXCO' BT  BY    ;                                               
  19.  BBX = 'EXCO' BT  BX    ;                                               
  20. *  DENSITE DE FORCES AUX NOEUDS                                         
  21.  DFX = COUR * BBY  * -1. ;                                              
  22.  DFY = COUR * BBX  ;                                                    
  23.   SI AXI ;
  24.  RX = COOR 1 BLOC ;
  25.  DFX = DFX * RX   ;
  26.  DFY = DFY * RX   ;
  27.   FINSI ;
  28. *  ON INTEGRE  SUR LES ELEMENTS                                         
  29.  RPTX = ('SOURCE' OBJ0 DFX  ) 'NOMC' 'FX' ;                             
  30.  RPTY = ('SOURCE' OBJ0 DFY  )'NOMC' 'FY' ;                              
  31.  RPT = RPTX + RPTY ;                                                    
  32.  CDG = 'BARY' BLOC ;                                                    
  33.  RFB = 'RESU' RPT ;                                                     
  34.  IP1 =  ( 'EXTR' RFB 'MAILLAGE'  ) 'POINT' 1 ;                          
  35.  RFB = RFB 'PLUS' ( CDG 'MOINS' IP1 ) ;                                 
  36.   'SI' axi  ;                                                           
  37.  rpt = rpt * -1. ;                                                      
  38.  rfb = rfb * -1. ;                                                      
  39.    'FINSI' ;                                                            
  40. *                                                                       
  41.  'FINPROC'  RPT RFB ;                                     
  42.  
  43.  
  44.  

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