Télécharger tiril2.eso

Retour à la liste

Numérotation des lignes :

tiril2
  1. C TIRIL2 SOURCE CB215821 26/08/24 21:18:39 12622
  2. C TIRILO SOURCE AM 93/01/12 21:36:57 781
  3. SUBROUTINE TIRIL2(RP,XL,YOUN,CISAIL,SD,SREDY,XINZ)
  4. C-----------------------------------------------------------------------
  5. C
  6. C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
  7. C POUTRE DROITE DE TIMOSHENKO EN AXES LOCAUX
  8. C "D APRES TEDEL"
  9. C
  10. C ENTREE
  11. C XL LONGUEUR
  12. C SD SECTION DROITE
  13. C SREDY SECTION REDUITE A L EFFORT TRANCHANT TY
  14. C XINZ INERTIE IZ
  15. C YOUN MODULE D YOUNG
  16. C CISAIL MODULE DE CISAILLEMENT
  17. C SORTIE
  18. C RP MATRICE DE RIGIDITE ELEMENTAIRE 6X6
  19. C
  20. C D'PRES TIRILO DC 98
  21. C-----------------------------------------------------------------------
  22. IMPLICIT INTEGER(I-N)
  23. IMPLICIT REAL*8(A-H,O-Z)
  24. DIMENSION RP(6,6)
  25. DATA P5,P25/.5D0,.25D0/
  26. C
  27. C INITIALISATIONS
  28. C
  29. DO 121 J=1,6
  30. DO 120 I=1,6
  31. RP(I,J)=0.D00
  32. 120 CONTINUE
  33. 121 CONTINUE
  34. RL=1.D00/XL
  35. RL2=RL*RL
  36. C
  37. C SECTION "REDUITE"
  38. C
  39. IF ( SREDY. NE.0. ) THEN
  40. SRY=SREDY
  41. ELSE
  42. SRY=SD
  43. ENDIF
  44. C
  45. ESL=YOUN*SD*RL
  46. EIZ=YOUN*XINZ
  47. GAY=CISAIL*SRY
  48. C
  49. C REMPLISSAGE DE LA MATRICE
  50. C FX1
  51. RP(1,1)=ESL
  52. RP(1,4)=-ESL
  53. C FY1
  54. RP(2,2)= GAY*RL
  55. RP(2,3)= P5*GAY
  56. RP(2,5)= -GAY*RL
  57. RP(2,6)= P5*GAY
  58. C MZ1
  59. RP(3,3)= P25*GAY*XL + EIZ*RL
  60. RP(3,5)= -P5*GAY
  61. RP(3,6)= P25*GAY*XL - EIZ*RL
  62. C FX2
  63. RP(4,4)= ESL
  64. C FY2
  65. RP(5,5)= GAY*RL
  66. RP(5,6)= -P5*GAY
  67. C MZ2
  68. RP(6,6)= P25*GAY*XL + EIZ*RL
  69. C
  70. DO 122 I=1,6
  71. DO 9 J=I,6
  72. RP(J,I)=RP(I,J)
  73. 9 CONTINUE
  74. 122 CONTINUE
  75. RETURN
  76. END
  77.  
  78.  
  79.  
  80.  
  81.  

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