Télécharger la plaquette Cast3M


fr en

Télécharger tifdl3.eso

Retour à la liste

Numérotation des lignes : 

tifdl3

  1. C TIFDL3    SOURCE    CHAT      05/01/13    03:38:51     5004
  2.        SUBROUTINE TIFDL3(RP,XL,LHOOK,DDHOOK,OMEG)
  3. C-----------------------------------------------------------------------
  4. C
  5. C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
  6. C    POUTRE DROITE DE TIMOSHENKO 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 RIGIDITE 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.       RL=1.D00/XL
  31. C
  32.       GA = OMEG*0.5D0*(DDHOOK( 2, 2) + DDHOOK( 3, 3))
  33.       EI = OMEG*0.5D0*(DDHOOK( 5, 5) + DDHOOK( 6, 6))
  34. C
  35. C  REMPLISSAGE DE LA MATRICE
  36. C  FX1
  37. C      RP( 1, 1)=DDHOOK( 1, 1)*RL
  38. C      RP( 1, 7)=-RP( 1, 1)
  39. C      RP( 7, 1)=-RP( 1, 1)
  40. C  FY1
  41.       RP( 2, 3)=GA*RL
  42.       RP( 3, 2)=-GA*RL
  43.       RP( 2, 5)=-GA/2.0D0
  44.       RP( 5, 2)=GA/2.0D0
  45.       RP( 2, 9)=-GA*RL
  46.       RP( 9, 2)=GA*RL
  47.       RP( 2,11)=-GA/2.0D0
  48.       RP(11, 2)=GA/2.0D0
  49. C  FZ1
  50.       RP( 3, 6)=-GA/2.0D0
  51.       RP( 3, 8)=GA*RL
  52.       RP( 3, 12)=-GA/2.0D0
  53.       RP( 6, 3)=GA/2.0D0
  54.       RP( 8, 3)=-GA*RL
  55.       RP( 12, 3)=GA/2.0D0
  56. C  MX1
  57. C      RP( 4, 4)=DDHOOK( 4, 4)*RL
  58. C      RP( 4,10)=-RP( 4, 4)
  59. C      RP( 10,4)=-RP( 4, 4)
  60. C  MY1
  61.       RP( 5, 6)=EI*RL+GA/(4.0D0*RL)
  62.       RP( 5, 8)=-GA/2.0D0
  63.       RP( 5,12)=-EI*RL+GA/(4.0D0*RL)
  64.       RP( 6, 5)=-EI*RL-GA/(4.0D0*RL)
  65.       RP( 8, 5)=GA/2.0D0
  66.       RP(12, 5)=EI*RL-GA/(4.0D0*RL)
  67. C  MZ1
  68.       RP( 6, 9)=-GA/2.0D0
  69.       RP( 6,11)=EI*RL-GA/(4.0D0*RL)
  70.       RP( 9, 6)=GA/2.0D0
  71.       RP( 11, 6)=-EI*RL+GA/(4.0D0*RL)
  72. C  FX2
  73. C      RP( 7, 7)=RP( 1, 1)
  74. C  FY2
  75.       RP( 8, 9)=GA*RL
  76.       RP( 9, 8)=-GA*RL
  77.       RP( 8, 11)=GA/2.0D0
  78.       RP( 11, 8)=-GA/2.0D0
  79. C  FZ2
  80.       RP( 9, 12)=GA/2.0D0
  81.       RP( 12, 9)=-GA/2.0D0
  82. C  MX2
  83. C      RP(10,10)=RP( 4, 4)
  84. C  MY2
  85.       RP(11,12)=EI*RL+GA/(4.0D0*RL)
  86. C  MZ2
  87.       RP(12,11)=-EI*RL-GA/(4.0D0*RL)
  88. C
  89.       RETURN
  90.       END
  91.  
  92.  
  93.  
  94.  

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