Télécharger podhl3.eso

Retour à la liste

Numérotation des lignes :

podhl3
  1. C PODHL3 SOURCE CHAT 05/01/13 02:16:14 5004
  2. SUBROUTINE PODHL3(RP,XL,VISC,CISAIL,SD,SREDY,SREDZ,
  3. . XINX,XINY,XINZ)
  4. C-----------------------------------------------------------------------
  5. C
  6. C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
  7. C POUTRE DROITE 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 SREDZ SECTION REDUITE A L EFFORT TRANCHANT TZ
  15. C XINX INERTIE IX
  16. C XINY INERTIE IY
  17. C XINZ INERTIE IZ
  18. C VISC MODULE DE VISCOSITE
  19. C CISAIL MODULE DE VISCOSITE EN CISAILLEMENT
  20. C SORTIE
  21. C RP MATRICE DE RIGIDITE ELEMENTAIRE 12X12
  22. C
  23. C VIENT DE TEDEL AOUT 84
  24. C-----------------------------------------------------------------------
  25. IMPLICIT INTEGER(I-N)
  26. IMPLICIT REAL*8(A-H,O-Z)
  27. DIMENSION RP(12,12)
  28. C
  29. C INITIALISATIONS
  30. C
  31. DO 120 I=1,12
  32. DO 120 J=1,12
  33. 120 RP(I,J)=0.D0
  34. C
  35. RL=1.D00/XL
  36. RL2=RL*RL
  37. X2=0.D00
  38. SRED = 0.5D0*(SREDY + SREDZ)
  39. XINYZ = 0.5D0*(XINY + XINZ)
  40. IF((SREDY.NE.0).AND.(SREDZ.NE.0))
  41. . X2= VISC/CISAIL*XINYZ*RL2/SRED
  42. X12= 1./(1.+12.*X2)
  43. C
  44. RL3=RL2*RL
  45. EIYZ=VISC*XINYZ
  46. C
  47. C REMPLISSAGE DE LA MATRICE
  48. C
  49. C
  50. RP(2,3)= 12.*EIYZ*RL3*X12
  51. RP(3,2)= -12.*EIYZ*RL3*X12
  52. RP(2,5)= -6.*EIYZ*RL2*X12
  53. RP(5,2)= 6.*EIYZ*RL2*X12
  54. RP(2,9)= -12. *EIYZ*RL3*X12
  55. RP(9,2)= 12. *EIYZ*RL3*X12
  56. RP(2,11)= -6.*EIYZ*RL2*X12
  57. RP(11,2)= 6.*EIYZ*RL2*X12
  58. C
  59. RP(3,6)= -6.*EIYZ*RL2*X12
  60. RP(3,8)= 12.*EIYZ*RL3*X12
  61. RP(3,12)= -6.*EIYZ*RL2*X12
  62. RP(6,3)= 6.*EIYZ*RL2*X12
  63. RP(8,3)= -12.*EIYZ*RL3*X12
  64. RP(12,3)= 6.*EIYZ*RL2*X12
  65. C
  66. RP(5,6)= 4.*EIYZ*RL
  67. . *(1.+3.*X2)*X12
  68. RP(5,8)= -6.*EIYZ*RL2*X12
  69. RP(5,12)=2.*EIYZ*RL
  70. .*(1.-6.*X2)*X12
  71. RP(6,5)= -4.*EIYZ*RL
  72. . *(1.+3.*X2)*X12
  73. RP(8,5)= 6.*EIYZ*RL2*X12
  74. RP(12,5)=-2.*EIYZ*RL
  75. .*(1.-6.*X2)*X12
  76. C
  77. RP(6,9)=-6.*EIYZ*RL2*X12
  78. RP(6,11)=-2.*EIYZ*RL
  79. . *(1.-6.*X2)*X12
  80. RP(9,6)=6.*EIYZ*RL2*X12
  81. RP(11,6)=2.*EIYZ*RL
  82. . *(1.-6.*X2)*X12
  83. C
  84. RP(8,9)= 12.*EIYZ*RL3*X12
  85. RP(8,11)= 6.*EIYZ*RL2*X12
  86. RP(9,8)= -12.*EIYZ*RL3*X12
  87. RP(11,8)= -6.*EIYZ*RL2*X12
  88. C
  89. RP(9,12)= 6.*EIYZ*RL2*X12
  90. RP(12,9)= -6.*EIYZ*RL2*X12
  91. C
  92. RP(11,12)= 4.*EIYZ*RL
  93. . *(1.D0+3.*X2)*X12
  94. C
  95. RP(12,11)= -4.*EIYZ*RL
  96. . *(1.D0+3.*X2)*X12
  97. C
  98. RETURN
  99. END
  100.  
  101.  
  102.  
  103.  

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