Procedures | Cast3M

   1 : * FORBLOC   PROCEDUR  BAZE      97/01/30    21:15:01     2468           
   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  'MODL' '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 :