Télécharger porilo.eso

Retour à la liste

Numérotation des lignes :

  1. C PORILO SOURCE AM 08/05/19 21:15:08 6108
  2. SUBROUTINE PORILO(RP,XL,YOUN,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 YOUN MODULE D YOUNG
  19. C CISAIL MODULE DE 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,144
  32. 120 RP(I,1)=0.D00
  33. RL=1.D00/XL
  34. RL2=RL*RL
  35. RL3=RL2*RL
  36. X2=0.D00
  37. X3=0.D00
  38.  
  39. IF( SREDZ. NE.0.) X2= YOUN/CISAIL*XINY*RL2/ SREDZ
  40. IF( SREDY. NE.0.) X3= YOUN/CISAIL*XINZ*RL2/ SREDY
  41. ESL=YOUN*SD*RL
  42. EIZ=YOUN*XINZ
  43. EIY=YOUN*XINY
  44. GJL=CISAIL*XINX*RL
  45.  
  46. X12= 1./(1.+12.*X2)
  47. X13= 1./(1.+12.*X3)
  48.  
  49. C
  50. C REMPLISSAGE DE LA MATRICE
  51. C
  52. RP(1,1)=ESL
  53. RP(1,7)=-ESL
  54. RP(2,2)= 12.*EIZ*RL3*X13
  55. RP(2,6)=6.*EIZ*RL2*X13
  56. RP(2,8)= -12. *EIZ*RL3*X13
  57. RP(2,12)= 6.*EIZ*RL2*X13
  58. RP(3,3)=12.*EIY*RL3*X12
  59. RP(3,5)= -6.*EIY*RL2*X12
  60. RP(3,9)= -12.*EIY*RL3*X12
  61. RP(3,11)= -6.*EIY*RL2*X12
  62. RP(4,4)= GJL
  63. RP(4,10)=-GJL
  64. RP(5,5)= 4.*EIY*RL
  65. . *(1.+3.*X2)*X12
  66. RP(5,9)= 6.*EIY*RL2*X12
  67. RP(5,11)=2.*EIY*RL
  68. .*(1.-6.*X2)*X12
  69. RP(6,6)= 4.*EIZ*RL
  70. . *(1.+3.*X3)*X13
  71. RP(6,8)=-6.*EIZ*RL2*X13
  72. RP(6,12)=2.*EIZ*RL
  73. . *(1.-6.*X3)*X13
  74. RP(7,7)=ESL
  75. RP(8,8)= 12.*EIZ*RL3*X13
  76. RP(8,12)=-6.*EIZ*RL2*X13
  77. RP(9,9)= 12.*EIY*RL3*X12
  78. RP(9,11)= 6.*EIY*RL2*X12
  79. RP(10,10)= GJL
  80. RP(11,11)= 4.*EIY*RL
  81. . *(1.+3.*X2)*X12
  82. RP(12,12)= 4.*EIZ*RL
  83. . *(1.+3.*X3)*X13
  84. DO 9 I=1,12
  85. DO 9 J=I,12
  86. 9 RP(J,I)=RP(I,J)
  87. RETURN
  88. END
  89.  
  90.  
  91.  
  92.  

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