Procedures | Cast3M

   1 : * COUR3D    PROCEDUR  CHAT      93/01/11    21:17:24     780            
   2 :  'DEBPROC' COUR3D BBBB*'MAILLAGE'  FFAC1*'MAILLAGE' FFAC2*'MAILLAGE'         
   3 :   MMM*'MOT'   J*'FLOTTANT'  NORJ/'LOGIQUE' ;                                    
   4 : ************************************************************************
   5 : *         CALCUL DU CHPO DES DENSITES DE COURANTS SUR UNE BOBINE       *
   6 : *         BBBB MAILLAGEGE DE LA BOBINE 3D                              *
   7 : *         FFAC1 FACE  D ENTREE DU COURANT                              *
   8 : *         FFAC2 GACE DE SORTIE DU COURANT                              *
   9 : *         MMMM MOT  AMP  ( AMPERES/UNITE SURFACE OU AT AMPERES TOTAUX  *
  10 : *         J FLOTTANT   DENSITE OU AMP TOTAUX   SUIVANT MMMM            *
  11 : *   SORTIE CHPO DES DENSITES DE COURANT AUX NOEUDS DE BBBB             *
  12 : *                 COMPOSANTES T,X   T,Y   T,Z
  13 : *         NORJ   LOGIQUE  VRAI SI  L ON RECTIFIE J DANS LES COURBURES  *
  14 : ************************************************************************
  15 :  'SI' ( 'NON' ( 'EXISTE' NORJ)) ; NORJ = FAUX ; 'FINSI' ;
  16 :  OBOB = BBBB 'MODL' 'THERMIQUE' 'ISOTROPE' ;                                  
  17 :  MATO = 'MATR' OBOB 'K' 1. ;                                              
  18 :  RRRR = 'COND' OBOB MATO ;                                                
  19 :  ENC1 = 'BLOQUER'  'T' FFAC1 ;                                            
  20 :  FFFF = 'FLUX'  OBOB  1.  FFAC2 ;                                         
  21 :  U = ('EXTR' ( 'RESU' FFFF ) 'MAILLAGE' )'POINT' 1 ;                            
  22 :  SUR = 'EXTR' ( 'RESU' FFFF) 'Q' U ;                                        
  23 :       'SI' ( 'EGA' MMM 'AT') ;                                              
  24 :  J =  J / SUR  ;                                                        
  25 :  'MESS' ' NOMBRES D AMPERES PAR BLOCS '  (J * SUR) ;                      
  26 :       'FINSI' ;                                                           
  27 :  POTE = 'RESOU' ( RRRR 'ET' ENC1 ) ( FFFF * J )  ;                          
  28 :  JJJ = 'CHAN' 'CHPO' ('GRAD' POTE OBOB) OBOB  ;                               
  29 :    'SI'  NORJ  ;
  30 : *****************************************************************
  31 : * normalisation des courants dans les coudes 
  32 : *****************************************************************
  33 :  RRR =( (('EXCO' JJJ T,X ) ** 2 ) + (('EXCO' JJJ T,Y ) ** 2)
  34 :          + (('EXCO' JJJ T,Z ) ** 2)    ) ** .5 ;
  35 :  JJJ = JJJ / RRR * J ;
  36 :    'FINSI' ;
  37 :  JJX = ('EXCO' JJJ T,X) 'NOMC' JX ;                                         
  38 :  JJY = ('EXCO' JJJ T,Y) 'NOMC' JY ;                                         
  39 :  JJZ = ('EXCO' JJJ T,Z) 'NOMC' JZ ;                                         
  40 :  JJJ = JJX + JJY + JJZ ;
  41 : *
  42 :   'FINPROC' JJJ ;