Télécharger la plaquette Cast3M


fr en

Télécharger tifil2.eso

Retour à la liste

Numérotation des lignes : 

tifil2

  1. C TIFIL2    SOURCE    CB215821  26/08/24    21:18:34     12622          
  2. C TIFIL2    SOURCE    AM        94/02/25    22:55:47     1061
  3.       SUBROUTINE TIFIL2(RP,XL,LHOOK,DDHOOK)
  4. C-----------------------------------------------------------------------
  5. C
  6. C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
  7. C    POUTRE DROITE DE TIMOSHENKO MULTI-FIBRE EN AXES LOCAUX
  8. C    "D APRES TEDEL"
  9. C
  10. C  ENTREE
  11. C    XL      LONGUEUR
  12. C    DDHOOK(LHOOK,LHOOK) MATRICE DE HOOK
  13. C
  14. C  SORTIE
  15. C    RP      MATRICE DE RIGIDITE ELEMENTAIRE 12X12
  16. C
  17. C D'APRES TIFILO DC 98
  18. C-----------------------------------------------------------------------
  19.       IMPLICIT INTEGER(I-N)
  20.       IMPLICIT REAL*8(A-H,O-Z)
  21.       DIMENSION RP(6,*)
  22.       DIMENSION DDHOOK(LHOOK,*)
  23. C
  24. C  INITIALISATIONS
  25. C
  26.       DO 121 I=1,6
  27.         DO 120 J=1,6
  28.           RP(I,J)=0.D0
  29.   120   CONTINUE
  30.   121 CONTINUE
  31. C
  32.       RL=1.D00/XL
  33. C
  34. C  REMPLISSAGE DE LA MATRICE
  35. C  FX1
  36.       RP( 1, 1)=DDHOOK( 1, 1)*RL
  37.       RP( 1, 3)=DDHOOK( 1, 3)*RL
  38.       RP( 1, 4)=-RP( 1, 1)
  39.       RP( 1, 6)=-RP( 1, 3)
  40. C  FY1
  41.       RP( 2, 2)=DDHOOK( 2, 2)*RL
  42.       RP( 2, 3)=DDHOOK( 2, 2)/2.0D0
  43.       RP( 2, 5)=-RP( 2, 2)
  44.       RP( 2, 6)= RP( 2, 3)
  45. C  MZ1
  46.       RP( 3, 3)=DDHOOK( 3, 3)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  47.       RP( 3, 4)=-DDHOOK( 3, 1)*RL
  48.       RP( 3, 5)=-DDHOOK( 2, 2)/2.0D0
  49.       RP( 3, 6)=-DDHOOK( 3, 3)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  50. C  FX2
  51.       RP( 4, 4)=RP( 1, 1)
  52.       RP( 4, 6)=-RP( 1, 6)
  53. C  FY2
  54.       RP( 5, 5)=RP( 2, 2)
  55.       RP( 5, 6)=-RP( 2, 3)
  56. C  MZ2
  57.       RP( 6, 6)=RP( 3, 3)
  58. C
  59.       DO 122 I=1,5
  60.         DO 9 J=I+1,6
  61.           RP(J,I)=RP(I,J)
  62.     9   CONTINUE
  63.   122 CONTINUE
  64.       RETURN
  65.  
  66.       END
  67.  
  68.  
  69.  
  70.  
  71.  

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