Procedures | Cast3M

   1 : * H_B       PROCEDUR  BAZE      97/01/30    21:15:07     2468           
   2 :  'DEBPROC' H_B MU0*'FLOTTANT'  EVO1/'EVOLUTION'  F/'FLOTTANT'  ;
   3 : ************************************************************************* 
   4 : *     definition de la courbe  de permeabilite magnetostatique 2d 3d    *
   5 : *                                                                       *
   6 : *   ENTREES :
   7 : *           MU0   permeabilité du vide       ( S.I.  4 pi 10-7 )        *       
   8 : *           EVO1  évolution    abscisse  H  ordonne  B                  *       
   9 : *           F   Packing factor     1 par defaut 
  10 : *************************************************************************       
  11 : *   SORTIES :
  12 : *                                                                       *
  13 : *    MU(H)   3d methode a 2 potentiels scalaires                        *  
  14 : * 1./MU(B)   2d plan ou axi potentiel vecteur                           *
  15 : *                                                                       *       
  16 : *                                                                       * 
  17 : *************************************************************************
  18 :   
  19 :  'SI' ('NON'  ('EXISTE'  EVO1 )) ; 
  20 :   TTTT = 'B(H)__STANDARD'  ;
  21 :   MESS '***********************************************************';
  22 :   MESS ' ON UTILISE UNE COURBE B(H) STANDARD  TESLA-AMPERE/METRE ' ;
  23 :   MESS '***********************************************************';  
  24 :  B= 'PROG' 0. .8 1.2 1.4 1.5  'PAS' .05 2.15 2.175 2.2 2.25 2.28 2.3
  25 :  2.3443
  26 :  2.3996 2.4905 2.5627 2.6706 2.8498 3.2074 3.5644 4.2782 4.8134 5.7052
  27 :  6.4186 7.4887 17.48 27.48   ;
  28 :  H= 'PROG' 0. 159.2 294.4 501.3 795.8 1154. 1795. 2862. 4383. 6044.
  29 :  8122.
  30 :  10590. 13610. 17220. 21170. 26750. 33760. 43800. 52440. 66000. 99470.
  31 :  120960. 141210. 169600. 212170. 283130. 339890. 425040. 566950.
  32 :  850760.
  33 :  1134600. 1702300. 2128000. 2837700. 3405100. 4256700. 12215578.
  34 :   20174457. ;  
  35 :  'SINON' ;
  36 :   TTTT = 'B(H)__UTILISATEUR'  ;
  37 :  H = 'EXTR' EVO1 'ABSC'  ; B = 'EXTR' EVO1 'ORDO'  ;                    
  38 :  'FINSI' ;    
  39 :   'SI' ('NON' ('EXISTE' F)) ; F= 1. ;'FINSI' ;
  40 :   
  41 :   TUTU = CHAIN TTTT '_______PACKING   '  F  ;      
  42 :   'MESS' '********************************************' ;
  43 :   'MESS'   TUTU ;
  44 :   'MESS' '********************************************' ;
  45 :    NN = DIME H ;                 
  46 :  H_MAX = MAXI H ;
  47 :  B_MAX = MAXI B ;
  48 :  B_EXT = 5.e4  ;  H_EXT = H_MAX + (( B_EXT - B_MAX)/MU0 ) ;
  49 :  H =  H ET ( 'PROG' H_EXT) ;
  50 :  B =  B ET ( 'PROG' B_EXT) ;  
  51 :    'SI' (f < 1.) ;
  52 :      ip = 1 ;
  53 :    BRED = PROG 0. ;
  54 :        'REPETER' BUB  nn   ;
  55 :      ip = ip + 1 ;
  56 :      Bi = EXTR B IP ; Hi = EXTR H IP ; MUI = Bi/Hi ;
  57 :      MURED =( MU0 * ( 1. - f )) +  (MUI *  f )  ;
  58 :      BRED =  BRED et (PROG ( Hi * MURED)) ;
  59 :      'FIN' BUB ;
  60 :   B = BRED ;
  61 :   'FINSI' ;
  62 :  
  63 :   BHEVO=EVOL MANU  'H'   H   'B' B ;
  64 : *  DESS BHEVO;
  65 :  H_1=  'ENLEVER' H  1 ;B_1= 'ENLEVER' B 1  ;                                    
  66 :  MUV= H_1 / B_1;TU= 'EXTRA' MUV 1;MUV= ( 'PROG' TU )  'ET' MUV ;
  67 :   IFOR = 'VALE' DIME ;         
  68 :  'SI' ('EGA' IFOR   3 );                                           
  69 :  BOBO= 'TEXTE'   ' MU(H)______POTENTIEL_SCALAIRE' ;                             
  70 :  'TITRE' BOBO ;                                                         
  71 :  REVOL = 'EVOL' MANU 'TEMPERATURE ' H   'CONDUCTIVITE ' ( MUV ** -1.) ; 
  72 :   'SINON' ;                                                               
  73 :  BOBO= 'TEXTE'   ' MU(B)______POTENTIEL_VECTEUR ';                              
  74 :  'TITRE' BOBO ;                                                         
  75 :  REVOL='EVOL' 'MANU' 'TEMPERATURE' B 'CONDUCTIVITE'  MUV   ;            
  76 :  'FINSI' ;
  77 : * 
  78 :  'SAUTER' 1 'LIGNE' ;
  79 :  'MESS' '********************************************' ;                        
  80 :  'MESS' BOBO ;
  81 :  'MESS' '********************************************' ;                        
  82 :  'SAUTER' 1 'LIGNE' ;                                                  
  83 :  'FINPROC' REVOL ; 
  84 : 
  85 :