Télécharger elpdkp.eso

Retour à la liste

Numérotation des lignes :

  1. C ELPDKP SOURCE CHAT 05/01/12 23:36:58 5004
  2. SUBROUTINE ELPDKP(P0,XN0,QI,XN,XT,XD,XNU,COP,CB)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-B,D-H,O-Z)
  5. IMPLICIT COMPLEX*16(C)
  6. ************************************************************************
  7. *
  8. * PLAQUES PAR EQUATION INTEGRALE:
  9. *
  10. * CALCUL DES OPERATEURS PLAQUE ET DE LEURS DERIVEES
  11. * SUR WD - WS
  12. *
  13. ************************************************************************
  14. *
  15. *
  16. DIMENSION P0(2)
  17. DIMENSION QI(2)
  18. DIMENSION COP(8)
  19. *
  20. DIMENSION XN (2)
  21. DIMENSION XN0(2)
  22. *
  23. DIMENSION XT (2)
  24. DIMENSION XT0(2)
  25. DIMENSION CG (5)
  26. DIMENSION CG0(5)
  27. DIMENSION XG (5)
  28. DIMENSION XR (2)
  29. *
  30. * 1- CALCUL DE R ET DU VECTEUR UNITAIRE SUIVANT R
  31. *
  32. R= ( ( QI(1) - P0(1) ) ** 2 + (QI(2) - P0(2) ) ** 2 ) **.5
  33. XR(1) = (QI(1) - P0(1)) / R
  34. XR(2) = (QI(2) - P0(2)) / R
  35. *
  36. *
  37. *
  38. *
  39. * 2- ON CALCULE LES FONCTIONS DE GREEN ET LEURS DERIVEES,R
  40. *
  41. CALL ELPDGF(R,CB,CG0)
  42. CALL ELPGFS(R,XG)
  43. DO 10 I=1,5
  44. CG(I) = CG0(I) - XG(I)
  45. 10 CONTINUE
  46. *
  47. * 3- CALCUL DE V
  48. *
  49. COP(1) = CG(1)
  50. *
  51. * 5- CALCUL DE MN= -D * (V,NN + NU* V,TT)
  52. *
  53. CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,2,CAUX1)
  54. CALL ELPDDR(CG,R,XR,XT,XT,XN,XN,2,CAUX2)
  55. COP(3) = -1* XD* ( CAUX1 + XNU*CAUX2)
  56. *
  57. * 6- CALCUL DE KN= -D * (V,NNN + (2 - NU)* V,NTT)
  58. * ON NE LE CALCULE PAS
  59. * CALL ELPDDR(CG,R,XR,XN,XN,XN,XN,3,CAUX1)
  60. * CALL ELPDDR(CG,R,XR,XN,XT,XT,XN,3,CAUX2)
  61. * COP(4) = -1* XD* ( CAUX1 + (2. - XNU) *CAUX2)
  62. *****
  63. * TERMES DERIVES
  64. *****
  65. *
  66. * 8- CALCUL DE V,NN0
  67. *
  68. CALL ELPDDR(CG,R,XR,XN,XN0,XN,XN,2,CAUX1)
  69. COP(6) = - CAUX1
  70. *
  71. * 10- CALCUL DE KN,N0
  72. *
  73. CALL ELPDDR(CG,R,XR,XN,XN,XN,XN0,4,CAUX1)
  74. CALL ELPDDR(CG,R,XR,XN,XT,XT,XN0,4,CAUX2)
  75. COP(8) = -1* XD* (- CAUX1 - (2. - XNU) *CAUX2)
  76.  
  77. RETURN
  78. END
  79.  
  80.  

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