Télécharger tifilo.eso

Retour à la liste

Numérotation des lignes :

  1. C TIFILO SOURCE PV 16/03/23 21:15:09 8864
  2. SUBROUTINE TIFILO(RP,XL,LHOOK,DDHOOK)
  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,*)
  22. DIMENSION DDHOOK(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. C REMPLISSAGE DE LA MATRICE
  33. C FX1
  34. RP( 1, 1)=DDHOOK( 1, 1)*RL
  35. RP( 1, 5)=DDHOOK( 1, 5)*RL
  36. RP( 1, 6)=DDHOOK( 1, 6)*RL
  37. RP( 1, 7)=-RP( 1, 1)
  38. RP( 1,11)=-RP( 1, 5)
  39. RP( 1,12)=-RP( 1, 6)
  40. C FY1
  41. RP( 2, 2)=DDHOOK( 2, 2)*RL
  42. RP( 2, 4)=DDHOOK( 2, 4)*RL
  43. RP( 2, 6)=DDHOOK( 2, 2)/2.0D0
  44. RP( 2, 8)=-RP( 2, 2)
  45. RP( 2,10)=-RP( 2, 4)
  46. RP( 2,12)= RP( 2, 6)
  47. C FZ1
  48. RP( 3, 3)=DDHOOK( 3, 3)*RL
  49. RP( 3, 4)=DDHOOK( 3, 4)*RL
  50. RP( 3, 5)=-DDHOOK( 3, 3)/2.0D0
  51. RP( 3, 9)=-RP( 3, 3)
  52. RP( 3,10)=-RP( 3, 4)
  53. RP( 3,11)= RP( 3, 5)
  54. C MX1
  55. RP( 4, 4)=DDHOOK( 4, 4)*RL
  56. RP( 4, 5)=-DDHOOK( 4, 3)/2.0D0
  57. RP( 4, 6)=DDHOOK( 4, 2)/2.0D0
  58. RP( 4, 8)=-DDHOOK( 4, 2)*RL
  59. RP( 4, 9)=-DDHOOK( 4, 3)*RL
  60. RP( 4,10)=-RP( 4, 4)
  61. RP( 4,11)=RP( 4, 5)
  62. RP( 4,12)=RP( 4, 6)
  63. C MY1
  64. RP( 5, 5)=DDHOOK( 5, 5)*RL+DDHOOK( 3, 3)/(4.0D0*RL)
  65. RP( 5, 6)=DDHOOK( 5, 6)*RL
  66. RP( 5, 7)=-DDHOOK( 5, 1)*RL
  67. RP( 5, 9)=DDHOOK( 3, 3)/2.0D0
  68. RP( 5,10)=DDHOOK( 3, 4)/2.0D0
  69. RP( 5,11)=-DDHOOK( 5, 5)*RL+DDHOOK( 3, 3)/(4.0D0*RL)
  70. RP( 5,12)=-RP( 5, 6)
  71. C MZ1
  72. RP( 6, 6)=DDHOOK( 6, 6)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  73. RP( 6, 7)=-DDHOOK( 6, 1)*RL
  74. RP( 6, 8)=-DDHOOK( 2, 2)/2.0D0
  75. RP( 6,10)=-DDHOOK( 2, 4)/2.0D0
  76. RP( 6,11)=-DDHOOK( 6, 5)*RL
  77. RP( 6,12)=-DDHOOK( 6, 6)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  78. C FX2
  79. RP( 7, 7)=RP( 1, 1)
  80. RP( 7,11)=RP( 1, 5)
  81. RP( 7,12)=RP( 1, 6)
  82. C FY2
  83. RP( 8, 8)=RP( 2, 2)
  84. RP( 8,10)=RP( 2, 4)
  85. RP( 8,12)=-RP( 2, 6)
  86. C FZ2
  87. RP( 9, 9)=RP( 3, 3)
  88. RP( 9,10)=RP( 3, 4)
  89. RP( 9,11)=-RP( 3, 5)
  90. C MX2
  91. RP(10,10)=RP( 4, 4)
  92. RP(10,11)=-RP( 4, 5)
  93. RP(10,12)=-RP( 4, 6)
  94. C MY2
  95. RP(11,11)=RP( 5, 5)
  96. RP(11,12)=RP( 5, 6)
  97. C MZ2
  98. RP(12,12)=RP( 6, 6)
  99. C
  100. DO 9 I=1,11
  101. DO 9 J=I+1,12
  102. 9 RP(J,I)=RP(I,J)
  103. RETURN
  104. END
  105.  
  106.  
  107.  

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