Télécharger elpdi3.eso

Retour à la liste

Numérotation des lignes :

  1. C ELPDI3 SOURCE CHAT 05/01/12 23:36:39 5004
  2. SUBROUTINE ELPDI3(P0,XN0
  3. & ,A,XN1,XN2,XT1,XT2,XLL
  4. & ,XD,XNU
  5. & ,CS2,CB,ISTAT)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-B,D-H,O-Z)
  8. IMPLICIT COMPLEX*16(C)
  9. ************************************************************************
  10. *
  11. * PLAQUES EN DYNAMIQUE PAR EQUATION INTEGRALE:
  12. *
  13. * CALCUL DU TERME DE COIN : 12 COEFFICIENTS A PLACER DANS LES
  14. * RELATIONS
  15. *
  16. ************************************************************************
  17. *
  18. *
  19. DIMENSION P0(2)
  20. DIMENSION A (2)
  21. *
  22. DIMENSION XR (2)
  23. DIMENSION XN0(2)
  24. DIMENSION XN1(2)
  25. DIMENSION XN2(2)
  26. *
  27. DIMENSION XT1(2)
  28. DIMENSION XT2(2)
  29. *
  30. DIMENSION XLL(4)
  31. *
  32. DIMENSION XG (5 )
  33. DIMENSION CG (5 )
  34. DIMENSION CS2 (12)
  35. *
  36. * 1- CALCUL DE GREEN AU COIN
  37. *
  38. R = ( (A(1)- P0(1))**2 + ( A(2) - P0(2)) ** 2 ) ** .5
  39. XR(1)=(A(1)- P0(1) )/ R
  40. XR(2)=(A(2)- P0(2) )/ R
  41. IF ( ISTAT .EQ. 1)THEN
  42. CALL ELPGFS(R,XG)
  43. DO 10 I=1,5
  44. CG(I) = XG(I)
  45. 10 CONTINUE
  46. ELSE
  47. CALL ELPDGF(R,CB,CG)
  48. ENDIF
  49. *
  50.  
  51. * 2- COEFFICIENTS SUR LA PREMIERE LIGNE D'EQUATION
  52. *
  53. CAUX1 = -1* XD* (1. - XNU ) * CG(1)
  54. CS2(1) = CAUX1 / ( ( XLL(3) + XLL(4) ) / 2)
  55. CS2(2) = - CS2(1)
  56. CS2(3) = - CAUX1 / ( ( XLL(1) + XLL(2) ) / 2)
  57. CS2(4) = - CS2(3)
  58. *
  59. * CALCUL DE MNT EN A
  60. *
  61. CALL ELPDDR(CG,R,XR,XN1,XT1,XT1,XT1,2,CVA1)
  62. CALL ELPDDR(CG,R,XR,XN2,XT2,XT1,XT1,2,CVA2)
  63. CMNT1 = - 1* ( XD* (1. - XNU ) ) * CVA1
  64. CMNT2 = - 1* ( XD* (1. - XNU ) ) * CVA2
  65. AX = XLL(2) / (XLL(2) + XLL(3))
  66. CS2(5) = AX * ( - CMNT2 + CMNT1)
  67. CS2(6) = (1 - AX) * ( - CMNT2 + CMNT1)
  68.  
  69. *
  70. * 3- COEFFICIENTS SUR LA DEUXIEME LIGNE
  71. *
  72. CALL ELPDDR(CG,R,XR,XN0,XN0,XN0,XN0,1,CVA )
  73. CAUX1 = -1* XD* (1. - XNU ) * (- CVA)
  74. CS2(7) = CAUX1 / ( ( XLL(3) + XLL(4) ) / 2)
  75. CS2(8) = - CS2(7)
  76. CS2(9) = - CAUX1 / ( ( XLL(1) + XLL(2) ) / 2)
  77. CS2(10) = - CS2(9)
  78. *
  79. * CALCUL DE MNT,N0 EN A
  80. *
  81. CALL ELPDDR(CG,R,XR,XN1,XT1,XN0,XN0,3,CVA1)
  82. CALL ELPDDR(CG,R,XR,XN2,XT2,XN0,XN0,3,CVA2)
  83. CMNT1 = - 1* ( XD* (1. - XNU ) ) * (- CVA1)
  84. CMNT2 = - 1* ( XD* (1. - XNU ) ) * (- CVA2 )
  85. CS2(11) = AX * ( - CMNT2 + CMNT1)
  86. CS2(12) = (1. - AX) * ( - CMNT2 + CMNT1)
  87. *
  88. RETURN
  89. END
  90.  
  91.  

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