Télécharger porim2.eso

Retour à la liste

Numérotation des lignes :

  1. C PORIM2 SOURCE AM 08/05/19 21:15:10 6108
  2. SUBROUTINE PORIM2(RP,XL,YOUN,CISAIL,SD,SREDY,XINZ,
  3. & 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 EN 2D
  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 12X12
  19. C
  20. C VIENT DE TEDEL AOUT 84
  21. C-----------------------------------------------------------------------
  22. IMPLICIT INTEGER(I-N)
  23. IMPLICIT REAL*8(A-H,O-Z)
  24. DIMENSION RP(6,6)
  25. DIMENSION DDHOOK(LHOOK,*)
  26. C
  27. C INITIALISATIONS
  28. C
  29. DO 120 I=1,36
  30. 120 RP(I,1)=0.D00
  31. RL=1.D00/XL
  32. RL2=RL*RL
  33. RL3=RL2*RL
  34. X2=0.D00
  35. X3=0.D00
  36.  
  37. IF(IMAT.EQ.1) THEN
  38. IF( SREDY. NE.0.) X3= YOUN/CISAIL*XINZ*RL2/ SREDY
  39. ESL=YOUN*SD*RL
  40. EIZ=YOUN*XINZ
  41. ELSE IF(IMAT.EQ.2) THEN
  42. IF(DDHOOK(2,2).NE.0.) X3= DDHOOK(3,3)*RL2/ DDHOOK(2,2)
  43. ESL=DDHOOK(1,1)*RL
  44. EIZ=DDHOOK(3,3)
  45. ENDIF
  46.  
  47. X12= 1./(1.+12.*X2)
  48. X13= 1./(1.+12.*X3)
  49.  
  50. C
  51. C REMPLISSAGE DE LA MATRICE
  52. C
  53. RP(1,1)=ESL
  54. RP(1,4)=-ESL
  55. RP(2,2)= 12.*EIZ*RL3*X13
  56. RP(2,3)=6.*EIZ*RL2*X13
  57. RP(2,5)= -12. *EIZ*RL3*X13
  58. RP(2,6)= 6.*EIZ*RL2*X13
  59. RP(3,3)= 4.*EIZ*RL
  60. . *(1.+3.*X3)*X13
  61. RP(3,5)=-6.*EIZ*RL2*X13
  62. RP(3,6)=2.*EIZ*RL
  63. . *(1.-6.*X3)*X13
  64. RP(4,4)=ESL
  65. RP(5,5)= 12.*EIZ*RL3*X13
  66. RP(5,6)=-6.*EIZ*RL2*X13
  67. RP(6,6)= 4.*EIZ*RL
  68. . *(1.+3.*X3)*X13
  69. DO 9 I=1,6
  70. DO 9 J=I,6
  71. 9 RP(J,I)=RP(I,J)
  72. RETURN
  73. END
  74.  
  75.  
  76.  
  77.  
  78.  

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