Télécharger la plaquette Cast3M


fr en

Télécharger kbick1.eso

Retour à la liste

Numérotation des lignes : 

kbick1

  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