Télécharger la plaquette Cast3M


fr en

Télécharger tirilo.eso

Retour à la liste

Numérotation des lignes : 

tirilo

  1. C TIRILO    SOURCE    PV        18/06/18    21:15:41     9860           
  2.       SUBROUTINE TIRILO(RP,XL,YOUN,CISAIL,SD,SREDY,SREDZ,
  3.      .                  XINX,XINY,XINZ)
  4. C-----------------------------------------------------------------------
  5. C
  6. C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
  7. C    POUTRE DROITE DE TIMOSHENKO EN AXES LOCAUX
  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    SREDZ   SECTION REDUITE A L EFFORT TRANCHANT TZ
  15. C    XINX    INERTIE IX
  16. C    XINY    INERTIE IY
  17. C    XINZ    INERTIE IZ
  18. C    YOUN    MODULE D YOUNG
  19. C    CISAIL  MODULE DE CISAILLEMENT
  20. C  SORTIE
  21. C    RP      MATRICE DE RIGIDITE ELEMENTAIRE 12X12
  22. C
  23. C     PP 3/9/92
  24. C-----------------------------------------------------------------------
  25.       IMPLICIT INTEGER(I-N)
  26.       IMPLICIT REAL*8(A-H,O-Z)
  27.       DIMENSION RP(12,12)
  28.       DATA P5,P25/.5D0,.25D0/
  29. C
  30. C  INITIALISATIONS
  31. C
  32.       DO 120 J=1,12
  33.       DO 120 I=1,12
  34.   120 RP(I,J)=0.D00
  35.       RL=1.D00/XL
  36.       RL2=RL*RL
  37. C
  38. C  SECTION "REDUITE"
  39. C
  40.       IF ( SREDZ. NE.0. ) THEN
  41.            SRZ=SREDZ
  42.       ELSE
  43.            SRZ=SD
  44.       ENDIF
  45.       IF ( SREDY. NE.0. ) THEN
  46.            SRY=SREDY
  47.       ELSE
  48.            SRY=SD
  49.       ENDIF
  50. C
  51.       ESL=YOUN*SD*RL
  52.       EIZ=YOUN*XINZ
  53.       GAZ=CISAIL*SRZ
  54.       EIY=YOUN*XINY
  55.       GAY=CISAIL*SRY
  56.       GJL=CISAIL*XINX*RL
  57. C
  58. C  REMPLISSAGE DE LA MATRICE
  59. C  FX1
  60.        RP(1,1)=ESL
  61.       RP(1,7)=-ESL
  62. C  FY1
  63.        RP(2,2)= GAY*RL
  64.       RP(2,6)=  P5*GAY
  65.        RP(2,8)= -GAY*RL
  66.       RP(2,12)= P5*GAY
  67. C  FZ1
  68.        RP(3,3)= GAZ*RL
  69.       RP(3,5)=  -P5*GAZ
  70.        RP(3,9)= -GAZ*RL
  71.       RP(3,11)= -P5*GAZ
  72. C  MX1
  73.        RP(4,4)= GJL
  74.       RP(4,10)=-GJL
  75. C  MY1
  76.       RP(5,5)=  P25*GAZ*XL + EIY*RL
  77.        RP(5,9)= P5*GAZ
  78.       RP(5,11)= P25*GAZ*XL - EIY*RL
  79. C  MZ1
  80.       RP(6,6)=  P25*GAY*XL + EIZ*RL
  81.        RP(6,8)= -P5*GAY
  82.       RP(6,12)= P25*GAY*XL - EIZ*RL
  83. C  FX2
  84.       RP(7,7)=  ESL
  85. C  FY2
  86.        RP(8,8)= GAY*RL
  87.       RP(8,12)= -P5*GAY
  88. C  FZ2
  89.        RP(9,9)= GAZ*RL
  90.       RP(9,11)= P5*GAZ
  91. C  MX2
  92.       RP(10,10)= GJL
  93. C  MY2
  94.       RP(11,11)= P25*GAZ*XL + EIY*RL
  95. C  MZ2
  96.       RP(12,12)= P25*GAY*XL + EIZ*RL
  97. C
  98.       DO 9 I=1,12
  99.       DO 9 J=I,12
  100.    9  RP(J,I)=RP(I,J)
  101.       RETURN
  102.       END
  103.  
  104.  
  105.  

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