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    CB215821  26/08/24    21:18:32     12622          
  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 121 I=1,12
  27.       DO 120 J=1,12
  28.       RP(I,J)=0.D00
  29.   120 CONTINUE
  30.   121 CONTINUE
  31. C
  32. C  REMPLISSAGE DE LA MATRICE
  33. C  FX1
  34.       RP( 1, 1)=DDHOOK( 1, 1)*XL/3.0D0
  35.       RP( 1, 5)=DDHOOK( 1, 5)*XL/3.0D0
  36.       RP( 1, 6)=DDHOOK( 1, 6)*XL/3.0D0
  37.       RP( 1, 7)=RP( 1, 1)/2.0D0
  38.       RP( 1,11)=RP( 1, 5)/2.0D0
  39.       RP( 1,12)=RP( 1, 6)/2.0D0
  40. C  FY1
  41.       RP( 2, 2)=DDHOOK( 2, 2)*XL/3.0D0
  42.       RP( 2, 4)=DDHOOK( 2, 4)*XL/3.0D0
  43.       RP( 2, 8)=RP( 2, 2)/2.0D0
  44.       RP( 2,10)=RP( 2, 4)/2.0D0
  45. C  FZ1
  46.       RP( 3, 3)=DDHOOK( 3, 3)*XL/3.0D0
  47.       RP( 3, 4)=DDHOOK( 3, 4)*XL/3.0D0
  48.       RP( 3, 9)=RP( 3, 3)/2.0D0
  49.       RP( 3,10)=RP( 3, 4)/2.0D0
  50. C  MX1
  51.       RP( 4, 4)=DDHOOK( 4, 4)*XL/3.0D0
  52.       RP( 4, 8)=DDHOOK( 4, 2)*XL/6.0D0
  53.       RP( 4, 9)=DDHOOK( 4, 3)*XL/6.0D0
  54.       RP( 4,10)=RP( 4, 4)/2.0D0
  55. C  MY1
  56.       RP( 5, 5)=DDHOOK( 5, 5)*XL/3.0D0
  57.       RP( 5, 6)=DDHOOK( 5, 6)*XL/3.0D0
  58.       RP( 5, 7)=DDHOOK( 5, 1)*XL/6.0D0
  59.       RP( 5,11)=RP( 5, 5)/2.0D0
  60.       RP( 5,12)=RP( 5, 6)/2.0D0
  61. C  MZ1
  62.       RP( 6, 6)=DDHOOK( 6, 6)*XL/3.0D0
  63.       RP( 6, 7)=DDHOOK( 6, 1)*XL/6.0D0
  64.       RP( 6,11)=DDHOOK( 6, 5)*XL/6.0D0
  65.       RP( 6,12)=RP( 6, 6)/2.0D0
  66. C  FX2
  67.       RP( 7, 7)=RP( 1, 1)
  68.       RP( 7,11)=RP( 1, 5)
  69.       RP( 7,12)=RP( 1, 6)
  70. C  FY2
  71.       RP( 8, 8)=RP( 2, 2)
  72.       RP( 8,10)=RP( 2, 4)
  73. C  FZ2
  74.       RP( 9, 9)=RP( 3, 3)
  75.       RP( 9,10)=RP( 3, 4)
  76. C  MX2
  77.       RP(10,10)=RP( 4, 4)
  78. C  MY2
  79.       RP(11,11)=RP( 5, 5)
  80.       RP(11,12)=RP( 5, 6)
  81. C  MZ2
  82.       RP(12,12)=RP( 6, 6)
  83. C
  84.       DO 122 I=1,11
  85.       DO 9 J=I+1,12
  86.       RP(J,I)= RP(I,J)
  87.     9 CONTINUE
  88.   122 CONTINUE
  89.       RETURN
  90.       END
  91.  
  92.  
  93.  
  94.  

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