Télécharger porimo.eso

Retour à la liste

Numérotation des lignes :

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

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