Télécharger la plaquette Cast3M


fr en

Télécharger tidhl3.eso

Retour à la liste

Numérotation des lignes : 

tidhl3

  1. C TIDHL3    SOURCE    CHAT      05/01/13    03:38:25     5004
  2.       SUBROUTINE TIDHL3(RP,XL,YOUN,CISAIL,SD,SREDY,SREDZ,
  3.      .                  XINX,XINY,XINZ,KERRE)
  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 I=1,12
  33.       DO 120 J=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.       IF (ABS (SRY - SRZ).GT.(MAX(SRY,SRZ)*1.D-8)) GOTO 999
  51.       IF (ABS (XINY - XINZ).GT.(MAX(XINY,XINZ)*1.D-8)) GOTO 999
  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
  61. C
  62. C
  63. C  REMPLISSAGE DE LA MATRICE
  64. C  FX1
  65. C      RP(1,1)=ESL
  66. C      RP(1,7)=-ESL
  67. C      RP(7,1)=-ESL
  68. C  FY1
  69.       RP(2,3)= GAY*RL
  70.       RP(3,2)= -GAY*RL
  71.       RP(2,5)=  -P5*GAY
  72.       RP(5,2)=  P5*GAY
  73.       RP(2,9)= -GAY*RL
  74.       RP(9,2)= GAY*RL
  75.       RP(2,11)= -P5*GAY
  76.       RP(11,2)= P5*GAY
  77. C  FZ1
  78.       RP(3,6)=  -P5*GAZ
  79.       RP(6,3)=  P5*GAZ
  80.       RP(3,8)= GAZ*RL
  81.       RP(8,3)= -GAZ*RL
  82.       RP(3,12)= -P5*GAZ
  83.       RP(12,3)= P5*GAZ
  84. C  MX1
  85. C      RP(4,4)= -GJL
  86. C      RP(4,10)=GJL
  87. C      RP(10,4)=-GJL
  88. C  MY1
  89.       RP(5,6)=  P25*GAZ*XL + EIY*RL
  90.       RP(6,5)=  -P25*GAZ*XL - EIY*RL
  91.       RP(5,8)= -P5*GAZ
  92.       RP(8,5)= P5*GAZ
  93.       RP(5,12)= P25*GAZ*XL - EIY*RL
  94.       RP(12,5)=  -P25*GAZ*XL + EIY*RL
  95. C  MZ1
  96.       RP(6,9)= -P5*GAY
  97.       RP(9,6)= P5*GAY
  98.       RP(6,11)= -P25*GAY*XL + EIZ*RL
  99.       RP(11,6)= P25*GAY*XL - EIZ*RL
  100. C  FX2
  101. C      RP(7,7)=  ESL
  102. C  FY2
  103.       RP(8,9)= GAY*RL
  104.       RP(9,8)= -GAY*RL
  105.       RP(8,11)= P5*GAY
  106.       RP(11,8)= -P5*GAY
  107. C  FZ2
  108.       RP(9,12)= P5*GAZ
  109.       RP(12,9)= -P5*GAZ
  110. C  MX2
  111. C      RP(10,10)= -GJL
  112. C  MY2
  113.       RP(11,12)= P25*GAZ*XL + EIY*RL
  114. C  MZ2
  115.       RP(12,11)= -P25*GAY*XL - EIZ*RL
  116. C
  117.       RETURN
  118. 999   KERRE=1
  119.       RETURN
  120.       END
  121.  
  122.  
  123.  
  124.  

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