Télécharger kbick1.eso

Retour à la liste

Numérotation des lignes :

  1. C KBICK1 SOURCE CB215821 16/04/21 21:17:29 8920
  2. SUBROUTINE KBICK1(X,XKI1)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. C Include contenant quelques constantes dont XPI :
  6. -INC CCREEL
  7. C*********************************************************************
  8. C Calcul des facteurs de forme en milieu absorbant en 2D-plan
  9. C SP appele par KRCOA2
  10. C Calcul de la fonction de Bickley-Naylor à l'ordre 3 en X
  11. C normalisée
  12. C (integrale de o à pi/2)
  13. C entree:
  14. C X : abscisse
  15. C NT : discrétisation de l'intervalle (0,pi/2)
  16. C resultat:
  17. C XKI1 : valeur de la fonction
  18. C
  19. C*********************************************************************
  20. C
  21. C méthode des trapèzes
  22. C découpage en NT segments identiques
  23. NT=100
  24.  
  25. DTETA = XPI/2./NT
  26.  
  27. G = 0.D0
  28. DO K=1,NT
  29. TETA = DTETA*(0.5 + (K-1))
  30. SINTETA = SIN(TETA)
  31. G1 = EXP(-X/SINTETA)
  32. G1 = G1 * SINTETA *SINTETA
  33. G = G + G1
  34. ENDDO
  35. XKI1= G*2.D0/NT
  36. C WRITE(6,*) 'BICKLEY: X KI1: ',X, XKI1
  37. RETURN
  38. END
  39.  
  40.  
  41.  
  42.  
  43.  

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