Esope source files

   1 : C ASEGTR    SOURCE    CHAT      05/01/12    21:25:35     5004           
   2 :       SUBROUTINE ASEGTR(COUR,HI,HE,RI,RE,QB,QH,YMIN,YMAX,
   3 :      &XP,ZP,AX,AY,AZ)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)
   6 :       REAL*8 MBMZP
   7 : C
   8 : -INC CCREEL
   9 : C*********************************************************************
  10 : C  Calcul du potentiel vecteur induit par un segment de courant COUR,
  11 : C  d'axe Oy, de section trapezoidale symetrique par rapport a xOy:
  12 : C     - rayons interieur et exterieur RI et RE (RI=-RE)
  13 : C     - hauteurs interieure et exterieure HI et HE
  14 : C     - pentes basse et haute QB et QH
  15 : C     - extension longitudinale YMIN a YMAX par rapport au plan xOz
  16 : C       contenant le point de coordonnees (XP,0,ZP)
  17 : C       et on calcule le champ (0,AY,0)
  18 : C*********************************************************************
  19 : C
  20 :       EXTERNAL FSEGA
  21 : C
  22 : C  calcul des grandeurs geometriques reduites
  23 : C
  24 :       AMU=4.E-7*XPI
  25 :       R=(RI+RE)/2.
  26 :       A=(RE-RI)/2.
  27 :       B=(HI+HE)/4.
  28 :       COEF=COUR*AMU/(16.*XPI*A*B)
  29 : C
  30 :       QH2P1=1.+QH**2
  31 :       BMZP=B-ZP
  32 :       QB2P1=1.+QB**2
  33 :       MBMZP=-B-ZP
  34 : C
  35 :       RES1=    -FSEGA(XP,YMAX,RE,QH2P1,QH,BMZP)
  36 :       RES1=RES1+FSEGA(XP,YMAX,RI,QH2P1,QH,BMZP)
  37 :       RES1=RES1+FSEGA(XP,YMAX,RE,QB2P1,QB,MBMZP)
  38 :       RES1=RES1-FSEGA(XP,YMAX,RI,QB2P1,QB,MBMZP)
  39 :       RES1=RES1+FSEGA(XP,YMIN,RE,QH2P1,QH,BMZP)
  40 :       RES1=RES1-FSEGA(XP,YMIN,RI,QH2P1,QH,BMZP)
  41 :       RES1=RES1-FSEGA(XP,YMIN,RE,QB2P1,QB,MBMZP)
  42 :       RES1=RES1+FSEGA(XP,YMIN,RI,QB2P1,QB,MBMZP)
  43 : C
  44 :       AX=0.
  45 :       AY=COEF*RES1
  46 :       AZ=0.
  47 :       RETURN
  48 :       END
  49 :  
  50 :  
  51 :  
  52 :