Télécharger aarctr.eso

Retour à la liste

Numérotation des lignes :

  1. C AARCTR SOURCE CB215821 17/11/30 21:15:01 9639
  2. SUBROUTINE AARCTR(SANGLE,COUR,HI,HE,RI,RE,QB,QH,TMIN,TMAX,
  3. &XP,ZP,AX,AY,AZ)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. C
  7. -INC CCREEL
  8. C*******************************************************************
  9. C Calcul du potentiel vecteur induit par un arc de cercle de courant COUR,
  10. C d'axe Oz, de section trapezoidale asymetrique par rapport a xOy:
  11. C - rayons interieur et exterieur RI et RE
  12. C - hauteurs interieure et exterieure HI et HE
  13. C - pentes basse et haute QB et QH
  14. C - extension angulaire TMIN a TMAX par rapport au plan meridien
  15. C xOz contenant le point de coordonnees (X,0,Z)
  16. C ou on calcule le champ (AX,AY,0)
  17. C*******************************************************************
  18. C
  19. SEGMENT SANGLE
  20. REAL*8 TETM(NT4)
  21. REAL*8 TETI(NT4)
  22. REAL*8 DTE(NT4)
  23. ENDSEGMENT
  24. C
  25. C champ magnetique selon Ox, Oy, Or et Oz cree par l'element d'arc de
  26. C longueur DTE situe en TETA, et grandeurs dependant de TETA pour
  27. C l'integration
  28. C
  29. C
  30. EXTERNAL FARCA
  31. C
  32. C calcul des grandeurs geometriques independantes de TETA
  33. C
  34. X=XP
  35. Z=ZP
  36. AMU=4.E-7*XPI
  37. R=(RI+RE)/2.
  38. A=(RE-RI)/2.
  39. B=(HI+HE)/4.
  40. COEF=COUR*AMU/(16.*XPI*A*B)
  41. RE2=RE**2
  42. RI2=RI**2
  43. C
  44. QH2=QH**2
  45. AH=1.+QH2
  46. QHR=QH*R
  47. QHRPZ=QHR+Z
  48. QHZ=QH*Z
  49. CH=X**2+B**2+QHR**2+Z**2-2.*B*QHRPZ+2.*QHR*Z
  50. DH=QH
  51. EH=B-QHRPZ
  52. C
  53. QB2=QB**2
  54. AB=1.+QB2
  55. QBR=QB*R
  56. QBRPZ=QBR+Z
  57. QBZ=QB*Z
  58. CB=X**2+B**2+QBR**2+Z**2+2.*B*QBRPZ+2.*QBR*Z
  59. DB=QB
  60. EB=-B-QBRPZ
  61. C
  62. C recherche des bornes d'integration
  63. C
  64. IDEB=0
  65. IFIN=0
  66. TTAI=-1.
  67. TTAF=-1.
  68. NT4=TETM(/1)
  69. C
  70. DO 1 I=1,NT4
  71. DMIN=TMIN-TETI(I)
  72. DMAX=TMAX-TETI(I)
  73. TTBI=SIGN(1.D0,DMIN)
  74. TTBF=SIGN(1.D0,DMAX)
  75. TESTI=0.5*(1.-TTAI*TTBI)
  76. TESTF=0.5*(1.-TTAF*TTBF)
  77. IDEB=IDEB+I*INT(TESTI)
  78. IFIN=IFIN+I*INT(TESTF)
  79. TTAI=TTBI
  80. TTAF=TTBF
  81. 1 CONTINUE
  82. C
  83. DTINI=DTE(IDEB)
  84. DTFIN=DTE(IFIN)
  85. TETINI=TETM(IDEB)
  86. TETFIN=TETM(IFIN)
  87. DTE(IDEB)=TETI(IDEB+1)-TMIN
  88. DTE(IFIN)=TMAX-TETI(IFIN)
  89. TETM(IDEB)=0.5*(TETI(IDEB+1)+TMIN)
  90. TETM(IFIN)=0.5*(TETI(IFIN)+TMAX)
  91. C
  92. AX=0.
  93. AY=0.
  94. AZ=0.
  95. C
  96. DO 21 I=IDEB,IFIN
  97. COST=COS(TETM(I))
  98. SINT=SIN(TETM(I))
  99. BH=-2.*(-B*QH+QH*QHRPZ+X*COST)
  100. BB=-2.*(B*QB+QB*QBRPZ+X*COST)
  101. C
  102. HR11=RE2*LOG(SQRT(AH*RE2+BH*RE+CH)+DH*RE+EH)
  103. & -RI2*LOG(SQRT(AH*RI2+BH*RI+CH)+DH*RI+EH)
  104. HR12=RE2*LOG(SQRT(AB*RE2+BB*RE+CB)+DB*RE+EB)
  105. & -RI2*LOG(SQRT(AB*RI2+BB*RI+CB)+DB*RI+EB)
  106. HR1=HR11-HR12
  107. HR1=0.5*HR1
  108. HR211=FARCA(RE,AH,BH,CH,DH,EH,X,COST,SINT)
  109. HR212=FARCA(RI,AH,BH,CH,DH,EH,X,COST,SINT)
  110. HR21=HR211-HR212
  111. HR221=FARCA(RE,AB,BB,CB,DB,EB,X,COST,SINT)
  112. HR222=FARCA(RI,AB,BB,CB,DB,EB,X,COST,SINT)
  113. HR22=HR221-HR222
  114. HR2=HR21-HR22
  115. HR=HR1-HR2
  116. HR=HR*DTE(I)
  117. C
  118. HX=HR*SINT
  119. HY=HR*COST
  120. C
  121. AX=AX-HX
  122. AY=AY+HY
  123. C
  124. 21 CONTINUE
  125. C
  126. AX=AX*COEF
  127. AY=AY*COEF
  128. C
  129. DTE(IDEB)=DTINI
  130. DTE(IFIN)=DTFIN
  131. TETM(IDEB)=TETINI
  132. TETM(IFIN)=TETFIN
  133. C
  134. RETURN
  135. END
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  

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