Télécharger poksl2.eso

Retour à la liste

Numérotation des lignes :

  1. C POKSL2 SOURCE CHAT 05/01/13 02:16:58 5004
  2. C POKSL2 SOURCE CASTEM 85/06/13 12:33:17
  3. SUBROUTINE POKSL2(RP,XL,SIGM)
  4. C-----------------------------------------------------------------------
  5. C
  6. C ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE GEOMETRIQUE
  7. C ELEMENTAIRE D UNE POUTRE DROITE EN AXES LOCAUX
  8. C D APRES TEDEL
  9. C
  10. C ENTREE
  11. C XL LONGUEUR
  12. C SIGM EFFORT NORMAL MOYEN
  13. C SORTIE
  14. C RP MATRICE DE RIGIDITE GEOMETRIQUE ELEMENTAIRE 12X12
  15. C
  16. C MAI 85 MILLARD
  17. C-----------------------------------------------------------------------
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8(A-H,O-Z)
  20. DIMENSION RP(6,6)
  21. C
  22. C INITIALISATIONS
  23. C
  24. DO 120 I=1,36
  25. 120 RP(I,1)=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,4)=-RN
  34. RP(4,4)=RN
  35. RP(2,2)=U1
  36. RP(2,3)=U2
  37. RP(2,5)=-U1
  38. RP(2,6)=U2
  39. RP(3,3)=U3
  40. RP(3,5)=-U2
  41. RP(3,6)=-U4
  42. RP(5,5)=U1
  43. RP(5,6)=-U2
  44. RP(6,6)=U3
  45. C
  46. C REMPLISSAGE DE LA MATRICE
  47. C
  48. DO 9 I=1,6
  49. DO 9 J=I,6
  50. 9 RP(J,I)=RP(I,J)
  51. RETURN
  52. END
  53.  
  54.  
  55.  

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