Télécharger elpdk1.eso

Retour à la liste

Numérotation des lignes :

  1. C ELPDK1 SOURCE CHAT 05/01/12 23:36:51 5004
  2. SUBROUTINE ELPDK1(XLP0
  3. & ,XD,XNU
  4. & ,NTRAP2,CS1,CB)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-B,D-H,O-Z)
  7. IMPLICIT COMPLEX*16(C)
  8. ************************************************************************
  9. *
  10. * PLAQUES PAR EQUATION INTEGRALE:
  11. *
  12. * ON INTEGRE PAR LA METHODE DES TRAPEZES EN STATIQUE
  13. * ET EN DYNAMIQUE
  14. * op (w-ws)
  15. ************************************************************************
  16. *
  17. *
  18. DIMENSION P0(2)
  19. DIMENSION QI(2)
  20. DIMENSION XOP(8)
  21. DIMENSION COPD(8)
  22. *
  23. DIMENSION XN (2)
  24. *
  25. DIMENSION XT (2)
  26. *
  27. *
  28. DIMENSION CS1 (8)
  29. *
  30. * 1- TAILLE DES TRAPEZES ( DEMI SEGMENT SEULEMENT)
  31. *
  32. X0 = XLP0/(4*NTRAP2)
  33. X1 = XLP0 /(2*NTRAP2)
  34. *
  35. * 2- POSITIONS DE REFERENCE
  36. *
  37. P0(1) = 0D0
  38. P0(2) = 0D0
  39. XN(1) = 1D0
  40. XN(2) = 0D0
  41. XT(1) = 0D0
  42. XT(2) = 1D0
  43. QI(1) = 0D0
  44. QI(2) = 0D0
  45. *
  46. * 2- INTEGRATION
  47. *
  48. DO 10 J=1,8
  49. CS1(J) = 0D0
  50. COPD(J) = 0D0
  51. 10 CONTINUE
  52. DO 100 I=1,NTRAP2
  53. QI(2) = X0 + (I-1)*X1
  54.  
  55. CALL ELPDKP(P0,XN,QI,XN,XT,XD,XNU,COPD,CB)
  56.  
  57. DO 110 J = 1,8
  58. CS1(J) = CS1(J) + 2*( X1 * COPD(J))
  59. 110 CONTINUE
  60. 100 CONTINUE
  61. *
  62. RETURN
  63. END
  64.  
  65.  

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