Télécharger la plaquette Cast3M


fr en

Télécharger tifalo.eso

Retour à la liste

Numérotation des lignes : 

tifalo

  1. C TIFALO    SOURCE    CHAT      05/01/13    03:38:33     5004
  2.       SUBROUTINE TIFALO(RP,XL,LHOOK,DDHOOK)
  3. C-----------------------------------------------------------------------
  4. C
  5. C    ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D UNE
  6. C    POUTRE DE TIMOSHENKO DROITE MULTI-FIBRE EN AXES LOCAUX
  7. C    D APRES TEDEL
  8. C
  9. C  ENTREE
  10. C    XL      LONGUEUR
  11. C    DDHOOK(LHOOK,LHOOK) MATRICE DE "HOOK"
  12. C
  13. C  SORTIE
  14. C    RP      MATRICE DE MASSE ELEMENTAIRE 12X12
  15. C
  16. C-----------------------------------------------------------------------
  17. C     Pierre Pegon (ISPRA) Juillet/Aout 1993
  18. C-----------------------------------------------------------------------
  19.       IMPLICIT INTEGER(I-N)
  20.       IMPLICIT REAL*8(A-H,O-Z)
  21.       DIMENSION RP(12,12)
  22.       DIMENSION DDHOOK(LHOOK,LHOOK)
  23. C
  24. C  INITIALISATIONS
  25. C
  26.       DO 120 I=1,12
  27.       DO 120 J=1,12
  28.   120 RP(I,J)=0.D00
  29. C
  30. C  REMPLISSAGE DE LA MATRICE
  31. C  FX1
  32.       RP( 1, 1)=DDHOOK( 1, 1)*XL/3.0D0
  33.       RP( 1, 5)=DDHOOK( 1, 5)*XL/3.0D0
  34.       RP( 1, 6)=DDHOOK( 1, 6)*XL/3.0D0
  35.       RP( 1, 7)=RP( 1, 1)/2.0D0
  36.       RP( 1,11)=RP( 1, 5)/2.0D0
  37.       RP( 1,12)=RP( 1, 6)/2.0D0
  38. C  FY1
  39.       RP( 2, 2)=DDHOOK( 2, 2)*XL/3.0D0
  40.       RP( 2, 4)=DDHOOK( 2, 4)*XL/3.0D0
  41.       RP( 2, 8)=RP( 2, 2)/2.0D0
  42.       RP( 2,10)=RP( 2, 4)/2.0D0
  43. C  FZ1
  44.       RP( 3, 3)=DDHOOK( 3, 3)*XL/3.0D0
  45.       RP( 3, 4)=DDHOOK( 3, 4)*XL/3.0D0
  46.       RP( 3, 9)=RP( 3, 3)/2.0D0
  47.       RP( 3,10)=RP( 3, 4)/2.0D0
  48. C  MX1
  49.       RP( 4, 4)=DDHOOK( 4, 4)*XL/3.0D0
  50.       RP( 4, 8)=DDHOOK( 4, 2)*XL/6.0D0
  51.       RP( 4, 9)=DDHOOK( 4, 3)*XL/6.0D0
  52.       RP( 4,10)=RP( 4, 4)/2.0D0
  53. C  MY1
  54.       RP( 5, 5)=DDHOOK( 5, 5)*XL/3.0D0
  55.       RP( 5, 6)=DDHOOK( 5, 6)*XL/3.0D0
  56.       RP( 5, 7)=DDHOOK( 5, 1)*XL/6.0D0
  57.       RP( 5,11)=RP( 5, 5)/2.0D0
  58.       RP( 5,12)=RP( 5, 6)/2.0D0
  59. C  MZ1
  60.       RP( 6, 6)=DDHOOK( 6, 6)*XL/3.0D0
  61.       RP( 6, 7)=DDHOOK( 6, 1)*XL/6.0D0
  62.       RP( 6,11)=DDHOOK( 6, 5)*XL/6.0D0
  63.       RP( 6,12)=RP( 6, 6)/2.0D0
  64. C  FX2
  65.       RP( 7, 7)=RP( 1, 1)
  66.       RP( 7,11)=RP( 1, 5)
  67.       RP( 7,12)=RP( 1, 6)
  68. C  FY2
  69.       RP( 8, 8)=RP( 2, 2)
  70.       RP( 8,10)=RP( 2, 4)
  71. C  FZ2
  72.       RP( 9, 9)=RP( 3, 3)
  73.       RP( 9,10)=RP( 3, 4)
  74. C  MX2
  75.       RP(10,10)=RP( 4, 4)
  76. C  MY2
  77.       RP(11,11)=RP( 5, 5)
  78.       RP(11,12)=RP( 5, 6)
  79. C  MZ2
  80.       RP(12,12)=RP( 6, 6)
  81. C
  82.       DO 9 I=1,11
  83.       DO 9 J=I+1,12
  84.   9   RP(J,I)= RP(I,J)
  85.       RETURN
  86.       END
  87.  
  88.  
  89.  

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