Télécharger la plaquette Cast3M


fr en

Télécharger tifilo.eso

Retour à la liste

Numérotation des lignes : 

tifilo

  1. C TIFILO    SOURCE    CB215821  26/08/24    21:18:35     12622          
  2.       SUBROUTINE TIFILO(RP,XL,LHOOK,DDHOOK)
  3. C-----------------------------------------------------------------------
  4. C
  5. C    ROUTINE DE CALCUL DE LA MATRICE DE RIGIDITE ELEMENTAIRE D UNE
  6. C    POUTRE DROITE DE TIMOSHENKO MULTI-FIBRE EN AXES LOCAUX
  7. C    "D APRES TEDEL"
  8. C
  9. C  ENTREE
  10. C    XL      LONGUEUR
  11. C    DDHOOK(LHOOK,LHOOK) MATRICE DE HOOK
  12. C
  13. C  SORTIE
  14. C    RP      MATRICE DE RIGIDITE ELEMENTAIRE 12X12
  15. C
  16. C-----------------------------------------------------------------------
  17. C     Pierre Pegon (ISPRA) Juillet/Aout 1993
  18. C-----------------------------------------------------------------------
  19.       IMPLICIT INTEGER(I-N)
  20.       IMPLICIT REAL*8(A-H,O-Z)
  21.       DIMENSION RP(12,*)
  22.       DIMENSION DDHOOK(LHOOK,*)
  23. C
  24. C  INITIALISATIONS
  25. C
  26.       DO 121 I=1,12
  27.         DO 120 J=1,12
  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, 5)=DDHOOK( 1, 5)*RL
  38.       RP( 1, 6)=DDHOOK( 1, 6)*RL
  39.       RP( 1, 7)=-RP( 1, 1)
  40.       RP( 1,11)=-RP( 1, 5)
  41.       RP( 1,12)=-RP( 1, 6)
  42. C  FY1
  43.       RP( 2, 2)=DDHOOK( 2, 2)*RL
  44.       RP( 2, 4)=DDHOOK( 2, 4)*RL
  45.       RP( 2, 6)=DDHOOK( 2, 2)/2.0D0
  46.       RP( 2, 8)=-RP( 2, 2)
  47.       RP( 2,10)=-RP( 2, 4)
  48.       RP( 2,12)= RP( 2, 6)
  49. C  FZ1
  50.       RP( 3, 3)=DDHOOK( 3, 3)*RL
  51.       RP( 3, 4)=DDHOOK( 3, 4)*RL
  52.       RP( 3, 5)=-DDHOOK( 3, 3)/2.0D0
  53.       RP( 3, 9)=-RP( 3, 3)
  54.       RP( 3,10)=-RP( 3, 4)
  55.       RP( 3,11)= RP( 3, 5)
  56. C  MX1
  57.       RP( 4, 4)=DDHOOK( 4, 4)*RL
  58.       RP( 4, 5)=-DDHOOK( 4, 3)/2.0D0
  59.       RP( 4, 6)=DDHOOK( 4, 2)/2.0D0
  60.       RP( 4, 8)=-DDHOOK( 4, 2)*RL
  61.       RP( 4, 9)=-DDHOOK( 4, 3)*RL
  62.       RP( 4,10)=-RP( 4, 4)
  63.       RP( 4,11)=RP( 4, 5)
  64.       RP( 4,12)=RP( 4, 6)
  65. C  MY1
  66.       RP( 5, 5)=DDHOOK( 5, 5)*RL+DDHOOK( 3, 3)/(4.0D0*RL)
  67.       RP( 5, 6)=DDHOOK( 5, 6)*RL
  68.       RP( 5, 7)=-DDHOOK( 5, 1)*RL
  69.       RP( 5, 9)=DDHOOK( 3, 3)/2.0D0
  70.       RP( 5,10)=DDHOOK( 3, 4)/2.0D0
  71.       RP( 5,11)=-DDHOOK( 5, 5)*RL+DDHOOK( 3, 3)/(4.0D0*RL)
  72.       RP( 5,12)=-RP( 5, 6)
  73. C  MZ1
  74.       RP( 6, 6)=DDHOOK( 6, 6)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  75.       RP( 6, 7)=-DDHOOK( 6, 1)*RL
  76.       RP( 6, 8)=-DDHOOK( 2, 2)/2.0D0
  77.       RP( 6,10)=-DDHOOK( 2, 4)/2.0D0
  78.       RP( 6,11)=-DDHOOK( 6, 5)*RL
  79.       RP( 6,12)=-DDHOOK( 6, 6)*RL+DDHOOK( 2, 2)/(4.0D0*RL)
  80. C  FX2
  81.       RP( 7, 7)=RP( 1, 1)
  82.       RP( 7,11)=RP( 1, 5)
  83.       RP( 7,12)=RP( 1, 6)
  84. C  FY2
  85.       RP( 8, 8)=RP( 2, 2)
  86.       RP( 8,10)=RP( 2, 4)
  87.       RP( 8,12)=-RP( 2, 6)
  88. C  FZ2
  89.       RP( 9, 9)=RP( 3, 3)
  90.       RP( 9,10)=RP( 3, 4)
  91.       RP( 9,11)=-RP( 3, 5)
  92. C  MX2
  93.       RP(10,10)=RP( 4, 4)
  94.       RP(10,11)=-RP( 4, 5)
  95.       RP(10,12)=-RP( 4, 6)
  96. C  MY2
  97.       RP(11,11)=RP( 5, 5)
  98.       RP(11,12)=RP( 5, 6)
  99. C  MZ2
  100.       RP(12,12)=RP( 6, 6)
  101. C
  102.       DO 122 I=1,11
  103.         DO 9 J=I+1,12
  104.           RP(J,I)=RP(I,J)
  105.     9   CONTINUE
  106.   122 CONTINUE
  107.       RETURN
  108.       END
  109.  
  110.  
  111.  
  112.  

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