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

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