Télécharger la plaquette Cast3M


fr en

Télécharger elpdk1.eso

Retour à la liste

Numérotation des lignes : 

elpdk1

  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