Télécharger la plaquette Cast3M


fr en

Télécharger asegtr.eso

Retour à la liste

Numérotation des lignes : 

asegtr

  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