Télécharger asegtr.eso

Retour à la liste

Numérotation des lignes :

  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.  
  53.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales