Télécharger la plaquette Cast3M


fr en

Télécharger afilci.eso

Retour à la liste

Numérotation des lignes : 

afilci

  1. C AFILCI    SOURCE    CB215821  25/03/25    21:15:01     12217          
  2.       SUBROUTINE AFILCI(SANGLE,COUR,R,TMIN,TMAX,X,Z,AX,AY,AZ)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5. C
  6. -INC CCREEL
  7. C
  8. C*******************************************************************
  9. C  Calcul du potentiel vecteur induit par un arc de cercle de courant COUR,
  10. C  d'axe Oz, de section filaire :
  11. C     - rayon R
  12. C     - extension angulaire TMIN a TMAX par rapport au plan meridien
  13. C       xOz contenant le point de coordonnees (X,0,Z)
  14. C       ou on calcule le champ (AX,AY,AZ)
  15. C*******************************************************************
  16. C
  17.       SEGMENT SANGLE
  18.          REAL*8 TETM(NT4)
  19.          REAL*8 TETI(NT4)
  20.          REAL*8 DTE(NT4)
  21.       ENDSEGMENT
  22. C
  23. C  champ magnetique selon Ox, Oy, Or et Oz cree par l'element d'arc de
  24. C  longueur DTE situe en TETA, et grandeurs dependant de TETA pour
  25. C  l'integration
  26. C
  27. C
  28.       XMUN  =-1.D0
  29.       XUN   = 1.D0
  30.       XDEUX = 2.D0
  31.       XQUAT = 4.D0
  32.  
  33.       XDEMI = 0.5D0
  34.  
  35.       XPREC = 4.D-7
  36.  
  37.       NT4=TETM(/1)
  38. C
  39. C  calcul des grandeurs geometriques independantes de TETA
  40. C
  41.       AMU=XPREC*XPI
  42.       COEF=COUR*AMU*R/(XQUAT*XPI)
  43.       X2=X**2
  44.       R2=R**2
  45.       Z2=Z**2
  46.       XR=X*R
  47. C
  48.       AX=XZERO
  49.       AY=XZERO
  50.       AZ=XZERO
  51. C
  52. C  cas de l'axe
  53. C
  54.       IF(A_EGALE_B(X,XZERO)) THEN
  55.         DIST= SQRT(R2+Z2)
  56.         AX  = COEF/DIST*(COS(TMAX)-COS(TMIN))
  57.         AY  = COEF/DIST*(SIN(TMAX)-SIN(TMIN))
  58.         RETURN
  59.       ENDIF
  60. C
  61. C  recherche des bornes d'integration
  62. C
  63.       IDEB=0
  64.       IFIN=0
  65.       TTAI=XMUN
  66.       TTAF=XMUN
  67. C
  68.       DO 1 I=1,NT4
  69.          DMIN=TMIN-TETI(I)
  70.          DMAX=TMAX-TETI(I)
  71.          TTBI=SIGN(XUN,DMIN)
  72.          TTBF=SIGN(XUN,DMAX)
  73.          TESTI=XDEMI*(XUN-TTAI*TTBI)
  74.          TESTF=XDEMI*(XUN-TTAF*TTBF)
  75.          IDEB=IDEB+I*INT(TESTI)
  76.          IFIN=IFIN+I*INT(TESTF)
  77.          TTAI=TTBI
  78.          TTAF=TTBF
  79.  1    CONTINUE
  80. C
  81.       DTINI=DTE(IDEB)
  82.       DTFIN=DTE(IFIN)
  83.       TETINI=TETM(IDEB)
  84.       TETFIN=TETM(IFIN)
  85.       DTE(IDEB)=TETI(IDEB+1)-TMIN
  86.       DTE(IFIN)=TMAX-TETI(IFIN)
  87.       TETM(IDEB)=XDEMI*(TETI(IDEB+1)+TMIN)
  88.       TETM(IFIN)=XDEMI*(TETI(IFIN)+TMAX)
  89. C
  90.       DO 21 I=IDEB,IFIN
  91.          COST=COS(TETM(I))
  92. C         SINT=SIN(TETM(I))
  93.          DIST=SQRT(R2+X2+Z2-XDEUX*XR*COST)
  94.          AY=AY+COST/DIST*DTE(I)
  95.  21   CONTINUE
  96.       AY=COEF*AY
  97. C
  98. C     INTEGRATION ANALYTIQUE DE AX
  99.       DIST1=SQRT(R2+X2+Z2-XDEUX*XR*COS(TMIN))
  100.       DIST2=SQRT(R2+X2+Z2-XDEUX*XR*COS(TMAX))
  101.       AX=-COEF/XR*(DIST2-DIST1)
  102. C
  103.       DTE(IDEB)=DTINI
  104.       DTE(IFIN)=DTFIN
  105.       TETM(IDEB)=TETINI
  106.       TETM(IFIN)=TETFIN
  107. C
  108.       RETURN
  109.       END
  110.  
  111.  

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