C ASEGTR SOURCE CHAT 05/01/12 21:25:35 5004 SUBROUTINE ASEGTR(COUR,HI,HE,RI,RE,QB,QH,YMIN,YMAX, &XP,ZP,AX,AY,AZ) IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) REAL*8 MBMZP C -INC CCREEL C********************************************************************* C Calcul du potentiel vecteur induit par un segment de courant COUR, C d'axe Oy, de section trapezoidale symetrique par rapport a xOy: C - rayons interieur et exterieur RI et RE (RI=-RE) C - hauteurs interieure et exterieure HI et HE C - pentes basse et haute QB et QH C - extension longitudinale YMIN a YMAX par rapport au plan xOz C contenant le point de coordonnees (XP,0,ZP) C et on calcule le champ (0,AY,0) C********************************************************************* C EXTERNAL FSEGA C C calcul des grandeurs geometriques reduites C AMU=4.E-7*XPI R=(RI+RE)/2. A=(RE-RI)/2. B=(HI+HE)/4. COEF=COUR*AMU/(16.*XPI*A*B) C QH2P1=1.+QH**2 BMZP=B-ZP QB2P1=1.+QB**2 MBMZP=-B-ZP C RES1= -FSEGA(XP,YMAX,RE,QH2P1,QH,BMZP) RES1=RES1+FSEGA(XP,YMAX,RI,QH2P1,QH,BMZP) RES1=RES1+FSEGA(XP,YMAX,RE,QB2P1,QB,MBMZP) RES1=RES1-FSEGA(XP,YMAX,RI,QB2P1,QB,MBMZP) RES1=RES1+FSEGA(XP,YMIN,RE,QH2P1,QH,BMZP) RES1=RES1-FSEGA(XP,YMIN,RI,QH2P1,QH,BMZP) RES1=RES1-FSEGA(XP,YMIN,RE,QB2P1,QB,MBMZP) RES1=RES1+FSEGA(XP,YMIN,RI,QB2P1,QB,MBMZP) C AX=0. AY=COEF*RES1 AZ=0. RETURN END