Télécharger poril2.eso

Retour à la liste

Numérotation des lignes :

  1. C PORIL2 SOURCE PV 18/06/18 21:15:25 9860
  2. SUBROUTINE PORIL2(RP,XL,YOUN,CISAIL,SD,SREDY,XINZ)
  3. C-----------------------------------------------------------------------
  4. C
  5. C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
  6. C POUTRE DROITE EN AXES LOCAUX EN 2D
  7. C D APRES TEDEL
  8. C
  9. C ENTREE
  10. C XL LONGUEUR
  11. C SD SECTION DROITE
  12. C SREDY SECTION REDUITE A L EFFORT TRANCHANT TY
  13. C XINZ INERTIE IZ
  14. C YOUN MODULE D YOUNG
  15. C CISAIL MODULE DE CISAILLEMENT
  16. C SORTIE
  17. C RP MATRICE DE RIGIDITE ELEMENTAIRE 12X12
  18. C
  19. C VIENT DE TEDEL AOUT 84
  20. C-----------------------------------------------------------------------
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8(A-H,O-Z)
  23. DIMENSION RP(6,6)
  24. C
  25. C INITIALISATIONS
  26. C
  27. DO 120 J=1,6
  28. DO 120 I=1,6
  29. 120 RP(I,J)=0.D00
  30. RL=1.D00/XL
  31. RL2=RL*RL
  32. RL3=RL2*RL
  33. X2=0.D00
  34. X3=0.D00
  35.  
  36. IF( SREDY. NE.0.) X3= YOUN/CISAIL*XINZ*RL2/ SREDY
  37. ESL=YOUN*SD*RL
  38. EIZ=YOUN*XINZ
  39.  
  40. X12= 1./(1.+12.*X2)
  41. X13= 1./(1.+12.*X3)
  42.  
  43. C
  44. C REMPLISSAGE DE LA MATRICE
  45. C
  46. RP(1,1)=ESL
  47. RP(1,4)=-ESL
  48. RP(2,2)= 12.*EIZ*RL3*X13
  49. RP(2,3)=6.*EIZ*RL2*X13
  50. RP(2,5)= -12. *EIZ*RL3*X13
  51. RP(2,6)= 6.*EIZ*RL2*X13
  52. RP(3,3)= 4.*EIZ*RL
  53. . *(1.+3.*X3)*X13
  54. RP(3,5)=-6.*EIZ*RL2*X13
  55. RP(3,6)=2.*EIZ*RL
  56. . *(1.-6.*X3)*X13
  57. RP(4,4)=ESL
  58. RP(5,5)= 12.*EIZ*RL3*X13
  59. RP(5,6)=-6.*EIZ*RL2*X13
  60. RP(6,6)= 4.*EIZ*RL
  61. . *(1.+3.*X3)*X13
  62. DO 9 I=1,6
  63. DO 9 J=I,6
  64. 9 RP(J,I)=RP(I,J)
  65. RETURN
  66. END
  67.  
  68.  
  69.  
  70.  
  71.  
  72.  

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