Télécharger elpdj1.eso

Retour à la liste

Numérotation des lignes :

  1. C ELPDJ1 SOURCE CHAT 05/01/12 23:36:43 5004
  2. SUBROUTINE ELPDJ1(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. *
  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,0d0)
  50. 10 CONTINUE
  51. DO 100 I=1,NTRAP2
  52. QI(2) = X0 + (I-1)*X1
  53. istat = 0
  54. CALL ELPDOP(P0,XN,QI,XN,XT,XD,XNU,COPD,CB,istat)
  55.  
  56. CALL ELPOP1(P0,XN,QI,XN,XT,XD,XNU,XOP)
  57.  
  58.  
  59.  
  60.  
  61. DO 110 J=1,8
  62. ccop =xop(j)*(1d0,0d0)
  63. CS1(J) = CS1(J) + 2*( X1 * (COPD(J)- ccop))
  64. 110 CONTINUE
  65. 100 CONTINUE
  66. *
  67. RETURN
  68. END
  69.  
  70.  

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