Télécharger poril2.eso

Retour à la liste

Numérotation des lignes :

  1. C PORIL2 SOURCE AM 08/05/19 21:15:07 6108
  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 I=1,36
  28. 120 RP(I,1)=0.D00
  29. RL=1.D00/XL
  30. RL2=RL*RL
  31. RL3=RL2*RL
  32. X2=0.D00
  33. X3=0.D00
  34.  
  35. IF( SREDY. NE.0.) X3= YOUN/CISAIL*XINZ*RL2/ SREDY
  36. ESL=YOUN*SD*RL
  37. EIZ=YOUN*XINZ
  38.  
  39. X12= 1./(1.+12.*X2)
  40. X13= 1./(1.+12.*X3)
  41.  
  42. C
  43. C REMPLISSAGE DE LA MATRICE
  44. C
  45. RP(1,1)=ESL
  46. RP(1,4)=-ESL
  47. RP(2,2)= 12.*EIZ*RL3*X13
  48. RP(2,3)=6.*EIZ*RL2*X13
  49. RP(2,5)= -12. *EIZ*RL3*X13
  50. RP(2,6)= 6.*EIZ*RL2*X13
  51. RP(3,3)= 4.*EIZ*RL
  52. . *(1.+3.*X3)*X13
  53. RP(3,5)=-6.*EIZ*RL2*X13
  54. RP(3,6)=2.*EIZ*RL
  55. . *(1.-6.*X3)*X13
  56. RP(4,4)=ESL
  57. RP(5,5)= 12.*EIZ*RL3*X13
  58. RP(5,6)=-6.*EIZ*RL2*X13
  59. RP(6,6)= 4.*EIZ*RL
  60. . *(1.+3.*X3)*X13
  61. DO 9 I=1,6
  62. DO 9 J=I,6
  63. 9 RP(J,I)=RP(I,J)
  64. RETURN
  65. END
  66.  
  67.  
  68.  
  69.  
  70.  

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