Télécharger pokslo.eso

Retour à la liste

Numérotation des lignes :

  1. C POKSLO SOURCE PV 18/06/18 21:15:23 9860
  2. SUBROUTINE POKSLO(RP,XL,SIGM)
  3. C-----------------------------------------------------------------------
  4. C
  5. C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE GEOMETRIQUE
  6. C ELEMENTAIRE D UNE POUTRE DROITE EN AXES LOCAUX
  7. C D APRES TEDEL
  8. C
  9. C ENTREE
  10. C XL LONGUEUR
  11. C SIGM EFFORT NORMAL MOYEN
  12. C SORTIE
  13. C RP MATRICE DE RIGIDITE GEOMETRIQUE ELEMENTAIRE 12X12
  14. C
  15. C MAI 85 MILLARD
  16. C-----------------------------------------------------------------------
  17. IMPLICIT INTEGER(I-N)
  18. IMPLICIT REAL*8(A-H,O-Z)
  19. DIMENSION RP(12,12)
  20. C
  21. C INITIALISATIONS
  22. C
  23. DO 120 J=1,12
  24. DO 120 I=1,12
  25. 120 RP(I,J)=0.D00
  26. RL=1.D00/XL
  27. RN=SIGM*RL
  28. U1=1.2D00*RN
  29. U2=0.1D00*SIGM
  30. U3=0.13333333333333333D00*SIGM*XL
  31. U4=0.03333333333333333D00*SIGM*XL
  32. RP(1 ,1)=RN
  33. RP(1,7)=-RN
  34. RP(7,7)=RN
  35. RP(2,2)=U1
  36. RP(2,6)=U2
  37. RP(2,8)=-U1
  38. RP(2,12)=U2
  39. RP(3,3)=U1
  40. RP(3,5)=-U2
  41. RP(3,9)=-U1
  42. RP(3,11)=-U2
  43. RP(5,5)=U3
  44. RP(5,9)=U2
  45. RP(5,11)=-U4
  46. RP(6,6)=U3
  47. RP(6,8)=-U2
  48. RP(6,12)=-U4
  49. RP(8,8)=U1
  50. RP(8,12)=-U2
  51. RP(9,9)=U1
  52. RP(9,11)=U2
  53. RP(11,11)=U3
  54. RP(12,12)=U3
  55. C
  56. C REMPLISSAGE DE LA MATRICE
  57. C
  58. DO 9 I=1,12
  59. DO 9 J=I,12
  60. 9 RP(J,I)=RP(I,J)
  61. RETURN
  62. END
  63.  
  64.  
  65.  

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