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    CB215821  26/08/24    21:18:31     12622          
  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 1000 I=1,12
  33.       DO 120 J=1,12
  34.       RP(I,J)=0.D00
  35.   120 CONTINUE
  36.  1000 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.       IF (ABS (SRY - SRZ).GT.(MAX(SRY,SRZ)*1.D-8)) GOTO 999
  53.       IF (ABS (XINY - XINZ).GT.(MAX(XINY,XINZ)*1.D-8)) GOTO 999
  54. C
  55.       ESL=YOUN*SD*RL
  56.       EIZ=YOUN*XINZ
  57.       GAZ=CISAIL*SRZ
  58.       EIY=YOUN*XINY
  59.       GAY=CISAIL*SRY
  60.       GJL=CISAIL*XINX*RL
  61. C
  62. C
  63. C
  64. C
  65. C  REMPLISSAGE DE LA MATRICE
  66. C  FX1
  67. C      RP(1,1)=ESL
  68. C      RP(1,7)=-ESL
  69. C      RP(7,1)=-ESL
  70. C  FY1
  71.       RP(2,3)= GAY*RL
  72.       RP(3,2)= -GAY*RL
  73.       RP(2,5)=  -P5*GAY
  74.       RP(5,2)=  P5*GAY
  75.       RP(2,9)= -GAY*RL
  76.       RP(9,2)= GAY*RL
  77.       RP(2,11)= -P5*GAY
  78.       RP(11,2)= P5*GAY
  79. C  FZ1
  80.       RP(3,6)=  -P5*GAZ
  81.       RP(6,3)=  P5*GAZ
  82.       RP(3,8)= GAZ*RL
  83.       RP(8,3)= -GAZ*RL
  84.       RP(3,12)= -P5*GAZ
  85.       RP(12,3)= P5*GAZ
  86. C  MX1
  87. C      RP(4,4)= -GJL
  88. C      RP(4,10)=GJL
  89. C      RP(10,4)=-GJL
  90. C  MY1
  91.       RP(5,6)=  P25*GAZ*XL + EIY*RL
  92.       RP(6,5)=  -P25*GAZ*XL - EIY*RL
  93.       RP(5,8)= -P5*GAZ
  94.       RP(8,5)= P5*GAZ
  95.       RP(5,12)= P25*GAZ*XL - EIY*RL
  96.       RP(12,5)=  -P25*GAZ*XL + EIY*RL
  97. C  MZ1
  98.       RP(6,9)= -P5*GAY
  99.       RP(9,6)= P5*GAY
  100.       RP(6,11)= -P25*GAY*XL + EIZ*RL
  101.       RP(11,6)= P25*GAY*XL - EIZ*RL
  102. C  FX2
  103. C      RP(7,7)=  ESL
  104. C  FY2
  105.       RP(8,9)= GAY*RL
  106.       RP(9,8)= -GAY*RL
  107.       RP(8,11)= P5*GAY
  108.       RP(11,8)= -P5*GAY
  109. C  FZ2
  110.       RP(9,12)= P5*GAZ
  111.       RP(12,9)= -P5*GAZ
  112. C  MX2
  113. C      RP(10,10)= -GJL
  114. C  MY2
  115.       RP(11,12)= P25*GAZ*XL + EIY*RL
  116. C  MZ2
  117.       RP(12,11)= -P25*GAY*XL - EIZ*RL
  118. C
  119.       RETURN
  120. 999   KERRE=1
  121.       RETURN
  122.       END
  123.  
  124.  
  125.  
  126.  
  127.  

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