Télécharger segtra.eso

Retour à la liste

Numérotation des lignes :

segtra
  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.  
  63.  

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