Télécharger la plaquette Cast3M


fr en

Télécharger bplqdr.eso

Retour à la liste

Numérotation des lignes : 

bplqdr

  1. C BPLQDR    SOURCE    CHAT      05/01/12    21:42:27     5004
  2.       SUBROUTINE BPLQDR(COUR,RI,RE,Y1,Y2,X,Z,BX,BY,BZ)
  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 champ induit par une plaque rectangulaire de courant COUR,
  10. C  d'axe Oy:
  11. C     - rayons interieur et exterieur RI et RE (RI=-RE)
  12. C     - extension longitudinale Y1 a Y2 par rapport au plan xOz
  13. C       contenant le point de coordonnees (X,0,Z)
  14. C       et on calcule le champ (BX,BY,BZ)
  15. C*********************************************************************
  16. C
  17. C
  18. C  calcul des grandeurs geometriques reduites
  19. C
  20.       AMU=4.E-7*XPI
  21.       R=(RI+RE)/2.
  22.       A=(RE-RI)/2.
  23.       COEF=COUR*AMU/(8.*XPI*A)
  24. C
  25.       XI1=Y1**2+Z**2
  26.       XI2=Y2**2+Z**2
  27.       XMA2=(X-A)**2
  28.       XPA2=(X+A)**2
  29.       Z2=Z**2
  30.       DEXT1=XMA2+Z2
  31.       DEXT2=XPA2+Z2
  32.       DEXT3=XMA2-Z2
  33.       DEXT4=XPA2-Z2
  34.       IF(Z.EQ.0.) THEN
  35.          BX=0.
  36.       ELSE
  37.          IF(DEXT1.EQ.0.) THEN
  38.             FAXI2=-XPI/2.
  39.             FAXI1=-XPI/2.
  40.          ELSE
  41.             FAXI2=ASIN((DEXT3*XI2-2.*Z2*XMA2)/(DEXT1*XI2))
  42.             FAXI1=ASIN((DEXT3*XI1-2.*Z2*XMA2)/(DEXT1*XI1))
  43.          ENDIF
  44.          IF(DEXT2.EQ.0.) THEN
  45.             FMAXI2=-XPI/2.
  46.             FMAXI1=-XPI/2.
  47.          ELSE
  48.             FMAXI2=ASIN((DEXT4*XI2-2.*Z2*XPA2)/(DEXT2*XI2))
  49.             FMAXI1=ASIN((DEXT4*XI1-2.*Z2*XPA2)/(DEXT2*XI1))
  50.          ENDIF
  51. C
  52.          BX=0.5*COEF*(FAXI2-FAXI1-FMAXI2+FMAXI1)
  53.       ENDIF
  54.       BY=0.
  55.       DCOIN1=SQRT(DEXT1+Y1**2)+Y1
  56.       DCOIN2=SQRT(DEXT1+Y2**2)+Y2
  57.       DCOIN3=SQRT(DEXT2+Y1**2)+Y1
  58.       DCOIN4=SQRT(DEXT2+Y2**2)+Y2
  59.       IF(DCOIN1*DCOIN2*DCOIN3*DCOIN4.NE.0) THEN
  60.          BZ=COEF*LOG((DCOIN4*DCOIN1)/(DCOIN2*DCOIN3))
  61.       ELSE
  62.          BZ=XGRAND
  63.       ENDIF
  64.       RETURN
  65.       END
  66.  
  67.  
  68.  
  69.  
  70.  

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