Esope source files

   1 : C SEGTRA    SOURCE    CHAT      05/01/13    03:13:18     5004           
   2 : C
   3 : C*********************************************************************
   4 : C
   5 :       SUBROUTINE SEGTRA(COUR,HI,HE,RI,RE,QB,QH,YMIN,YMAX,
   6 :      &X,Z,BX,BY,BZ)
   7 :       IMPLICIT INTEGER(I-N)
   8 :       IMPLICIT REAL*8(A-H,O-Z) 
   9 : C
  10 : -INC CCREEL
  11 : C
  12 : C*********************************************************************
  13 : C  Calcul du champ induit par un segment de courant COUR,
  14 : C  d'axe Oy, de section trapezoidale symetrique par rapport a xOy:
  15 : C     - rayons interieur et exterieur RI et RE (RI=-RE)
  16 : C     - hauteurs interieure et exterieure HI et HE
  17 : C     - pentes basse et haute QB et QH
  18 : C     - extension longitudinale YMIN a YMAX par rapport au plan xOz
  19 : C       contenant le point de coordonnees (X,0,Z)
  20 : C       oº on calcule le champ (BX,BY,BZ)
  21 : C*********************************************************************
  22 : C
  23 : C
  24 :       EXTERNAL FONCF
  25 :       EXTERNAL FONCG
  26 : C
  27 : C
  28 : C  calcul des grandeurs geometriques reduites
  29 : C
  30 :       AMU=4.E-7*XPI
  31 :       R=(RI+RE)/2.
  32 :       A=(RE-RI)/2.
  33 :       B=(HI+HE)/4.
  34 :       COEF=COUR*AMU/(16.*XPI*A*B)
  35 :       Q2B=SQRT(QB*QB+1.)
  36 :       Q2H=SQRT(QH*QH+1.)
  37 : C
  38 :          BX=COEF*
  39 :      &   ((FONCF(QH,+A,YMAX,B-Z,X)-FONCF(QH,+A,YMIN,B-Z,X)
  40 :      &   -FONCF(QH,-A,YMAX,B-Z,X)+FONCF(QH,-A,YMIN,B-Z,X))/Q2H
  41 :      &   +(-FONCF(-QB,+A,YMAX,B+Z,X)+FONCF(-QB,+A,YMIN,B+Z,X)
  42 :      &   +FONCF(-QB,-A,YMAX,B+Z,X)-FONCF(-QB,-A,YMIN,B+Z,X))/Q2B)
  43 : C
  44 :          BY=0.
  45 : C
  46 :          BZ=COEF*
  47 :      &   (FONCG(QH,+A,YMAX,B-Z,X)-FONCG(QH,+A,YMIN,B-Z,X)
  48 :      &   -FONCG(QH,-A,YMAX,B-Z,X)+FONCG(QH,-A,YMIN,B-Z,X)
  49 :      &   +FONCG(-QB,+A,YMAX,B+Z,X)-FONCG(-QB,+A,YMIN,B+Z,X)
  50 :      &   -FONCG(-QB,-A,YMAX,B+Z,X)+FONCG(-QB,-A,YMIN,B+Z,X))
  51 :      &        +COEF*QH/Q2H*
  52 :      &   (FONCF(QH,+A,YMAX,B-Z,X)-FONCF(QH,+A,YMIN,B-Z,X)
  53 :      &   -FONCF(QH,-A,YMAX,B-Z,X)+FONCF(QH,-A,YMIN,B-Z,X))
  54 :      &        -COEF*QB/Q2B*
  55 :      &   (FONCF(-QB,+A,YMAX,B+Z,X)-FONCF(-QB,+A,YMIN,B+Z,X)
  56 :      &   -FONCF(-QB,-A,YMAX,B+Z,X)+FONCF(-QB,-A,YMIN,B+Z,X))
  57 :       RETURN
  58 :       END
  59 :  
  60 :  
  61 :  
  62 :