Procedures | Cast3M

   1 : * BIOVOL    PROCEDUR  CHAT      93/01/08    21:15:08     780            
   2 :  'DEBPROC' BIOVOL  BOBI*'MAILLAGE' PCAL*'MAILLAGE' J*'CHPOINT'  
   3 :                                             MU0* 'FLOTTANT'  ;
   4 : ***********************************************************************
   5 : *    CALCUL DU CHAMP BS  CREE SUR PCAL PAR L INDUCTEUR VOLUMIQUE BOBI *
   6 : *    PARCOURU PAR DES COURANTS DES DENSITES J (COMPOSANTES JX JY JZ)  *
   7 : *    MU0  PERMEABILITE DU VIDE                                        *
   8 : ***********************************************************************    
   9 :    fact = 1. / (4. * pi) ;
  10 :    nni = bobi 'NBNO' ;
  11 :    JIX = 'EXCO' J 'JX' ;
  12 :    JIY = 'EXCO' J 'JY';
  13 :    JIZ  = 'EXCO' J 'JZ' ;
  14 : *   
  15 :    X Y Z = 'COOR' BOBI ;
  16 :    NGEO = PCAL 'NBNO' ;IGEO = 0;
  17 :    'REPETER' BO_GEO  NGEO ;
  18 :      IGEO = IGEO + 1 ;
  19 :      PCAL1 = PCAL 'POINT' IGEO ;
  20 :      XP YP ZP = 'COOR' PCAL1 ;
  21 :      RX = XP - X ;RY = YP - Y ; RZ = ZP  - Z ;
  22 :      RM2 = (RX * RX) + ( RY * RY ) + ( RZ * RZ ) ;
  23 :      R3 = RM2 ** 1.5 ;
  24 :      HX = (((JIY * RZ ) - ( JIZ * RY) )/ R3) ; 
  25 :      HY = (((JIZ * RX ) - ( JIX * RZ) )/ R3) ;
  26 :      HZ  = (((JIX  * RY  ) - ( JIY  * RX ) )/ R3) ; 
  27 :      HPT =(HX 'NOMC' 'HX')+( HY 'NOMC' 'HY' )+( HZ 'NOMC' 'HZ' ) ;
  28 :      OBJM = 'MODL' bobi  'THERMIQUE' 'ISOTROPE' ;
  29 :      HEPT = 'CHAN' CHAM HPT  OBJM  'NOEUDS' ;
  30 :      hheX = 'EXCO' hept 'HX' ;  
  31 :      hhey = 'EXCO'hept 'HY' ;  
  32 :      hhez = 'EXCO'  hept 'HZ' ;  
  33 :      HTX = 'INTG' HHEX OBJM ;
  34 :      HTY = 'INTG' HHEY OBJM ;
  35 :      HTZ = 'INTG' HHEZ OBJM ;
  36 : * 
  37 :      HT ='MANU' 'CHPO' pcal1 3 'HX'  HTX 'HY' HTY  'HZ' HTZ ; 
  38 :      'SI' (IGEO 'EGA' 1) ;
  39 :         HGEO = HT ;
  40 :      'SINON' ;
  41 :         HGEO =HGEO + HT ;
  42 :      'FINSI' ;
  43 :  'FIN' BO_GEO ;
  44 : ******************************************************
  45 :  HGEO  = HGEO * fact ;
  46 :  'FINPROC' HGEO ;