Procedures | Cast3M

   1 : * INDUCTIO  PROCEDUR  CHAT      05/02/04    21:15:05     5023           
   2 :  'DEBPROC' INDUCTIO GEO*'MAILLAGE'  SOL1*'CHPOINT'  AXI*'LOGIQUE';
   3 : ******************************************************************
   4 : *  estimation de B aux points de Gauss et puis chan chpo
   5 : *                 2D  UNIQUEMENT                                        
   6 : *          calcul; de l induction  en potentiel vecteur        *        
   7 : *          GEO  maillage sur lequel on recherche B             *        
   8 : *          SOL1   solution en potentiel vecteur                *        
   9 : *          AXI    logique  vrai si axi                         *        
  10 : ****************************************************************        
  11 :  OBJ0 = GEO 'MODL' 'THERMIQUE' 'ISOTROPE' ;                             
  12 :  DERIV   = 'GRAD' OBJ0 SOL1 ;                                           
  13 :  BXX= ('EXCO' 'T,Y' DERIV   'SCAL' )    ;                               
  14 :  BYY= ('EXCO' 'T,X' DERIV   'SCAL' )    ; 
  15 :  BXX = BXX 'CHAN' 'TYPE' 'SCALAIRE';
  16 :  BYY = BYY 'CHAN' 'TYPE' 'SCALAIRE';                              
  17 :               'SI'   AXI   ;                                            
  18 :  'MESS' '*  axisymetrique   '   ;                                       
  19 :  1SRAY = (COOR 1 DERIV ) ** -1. ;                                       
  20 :  BXX=EXCO ( BXX  * 1SRAY * -1. ) 'SCAL' 'BX' ;                          
  21 :  BYY=EXCO ( BYY  * 1SRAY       ) 'SCAL' 'BY' ; 
  22 :  BXX = BXX 'CHAN' 'TYPE' 'SCALAIRE';
  23 :  BYY = BYY 'CHAN' 'TYPE' 'SCALAIRE';                               
  24 :               'SINON' ;                                                 
  25 : 'MESS' '*  probleme plan   ';                                           
  26 :  BXX =EXCO  BXX 'SCAL'    'BX' ;                                        
  27 :  BYY =EXCO (BYY * -1. ) 'SCAL' 'BY'  ;  
  28 :  BXX = BXX 'CHAN' 'TYPE' 'SCALAIRE';
  29 :  BYY = BYY 'CHAN' 'TYPE' 'SCALAIRE';                                      
  30 :               'FINSI' ;                                                 
  31 :  BTOT = BXX + BYY  ;                                                    
  32 :  BTOT = 'CHAN' 'CHPO' BTOT OBJ0 ;                                       
  33 :  'FINPROC' BTOT  ;                                                      
  34 : 
  35 :  
  36 :