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    CB215821  26/08/24    21:18:40     12622          
  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 121 J=1,12
  33.         DO 120 I=1,12
  34.           RP(I,J)=0.D00
  35.   120   CONTINUE
  36.   121 CONTINUE
  37.       RL=1.D00/XL
  38.       RL2=RL*RL
  39. C
  40. C  SECTION "REDUITE"
  41. C
  42.       IF ( SREDZ. NE.0. ) THEN
  43.            SRZ=SREDZ
  44.       ELSE
  45.            SRZ=SD
  46.       ENDIF
  47.       IF ( SREDY. NE.0. ) THEN
  48.            SRY=SREDY
  49.       ELSE
  50.            SRY=SD
  51.       ENDIF
  52. C
  53.       ESL=YOUN*SD*RL
  54.       EIZ=YOUN*XINZ
  55.       GAZ=CISAIL*SRZ
  56.       EIY=YOUN*XINY
  57.       GAY=CISAIL*SRY
  58.       GJL=CISAIL*XINX*RL
  59. C
  60. C  REMPLISSAGE DE LA MATRICE
  61. C  FX1
  62.        RP(1,1)=ESL
  63.       RP(1,7)=-ESL
  64. C  FY1
  65.        RP(2,2)= GAY*RL
  66.       RP(2,6)=  P5*GAY
  67.        RP(2,8)= -GAY*RL
  68.       RP(2,12)= P5*GAY
  69. C  FZ1
  70.        RP(3,3)= GAZ*RL
  71.       RP(3,5)=  -P5*GAZ
  72.        RP(3,9)= -GAZ*RL
  73.       RP(3,11)= -P5*GAZ
  74. C  MX1
  75.        RP(4,4)= GJL
  76.       RP(4,10)=-GJL
  77. C  MY1
  78.       RP(5,5)=  P25*GAZ*XL + EIY*RL
  79.        RP(5,9)= P5*GAZ
  80.       RP(5,11)= P25*GAZ*XL - EIY*RL
  81. C  MZ1
  82.       RP(6,6)=  P25*GAY*XL + EIZ*RL
  83.        RP(6,8)= -P5*GAY
  84.       RP(6,12)= P25*GAY*XL - EIZ*RL
  85. C  FX2
  86.       RP(7,7)=  ESL
  87. C  FY2
  88.        RP(8,8)= GAY*RL
  89.       RP(8,12)= -P5*GAY
  90. C  FZ2
  91.        RP(9,9)= GAZ*RL
  92.       RP(9,11)= P5*GAZ
  93. C  MX2
  94.       RP(10,10)= GJL
  95. C  MY2
  96.       RP(11,11)= P25*GAZ*XL + EIY*RL
  97. C  MZ2
  98.       RP(12,12)= P25*GAY*XL + EIZ*RL
  99. C
  100.       DO 122 I=1,12
  101.         DO 9 J=I,12
  102.           RP(J,I)=RP(I,J)
  103.     9   CONTINUE
  104.   122 CONTINUE
  105.       RETURN
  106.       END
  107.  
  108.  
  109.  
  110.  

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