Télécharger porimo.eso

Retour à la liste

Numérotation des lignes :

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

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