Télécharger la plaquette Cast3M


fr en

Télécharger tifclo.eso

Retour à la liste

Numérotation des lignes : 

tifclo

  1. C TIFCLO    SOURCE    KK2000    05/05/25    21:15:27     5088
  2.       SUBROUTINE TIFCLO(RP,XL,LHOOK,DDHOOK,VROTL)
  3. C-----------------------------------------------------------------------
  4. C
  5. C    ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR CENTRIFUGE ELEMENTAIRE
  6. C    D UNE POUTRE DE TIMOSHENKO DROITE MULTI-FIBRE EN AXES LOCAUX
  7. C
  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 RAIDEUR ELEMENTAIRE 12X12
  15. C
  16. C-----------------------------------------------------------------------
  17. C     Didier Combescure 2003
  18. C-----------------------------------------------------------------------
  19.       IMPLICIT INTEGER(I-N)
  20.       IMPLICIT REAL*8(A-H,O-Z)
  21. -INC CCREEL
  22.       DIMENSION RP(12,12),VROTL(3),ROME(6,6)
  23.       DIMENSION DDHOOK(LHOOK,LHOOK)
  24. C
  25. C  INITIALISATIONS
  26. C
  27.       DO 120 I=1,12
  28.       DO 120 J=1,12
  29.   120 RP(I,J)=0.D00
  30. C
  31.       RHOSD = DDHOOK( 1, 1)
  32.       FAC=RHOSD*XL/6.D0
  33.       FAC2=FAC*2.D0
  34. C
  35.       DO 121 I=1,6
  36.       DO 121 J=1,6
  37.   121 ROME(I,J)=0.D00
  38. C
  39.       ROME(1,1) = (-1.)*((VROTL(2)**2) + (VROTL(3)**2))
  40.       ROME(2,2) = (-1.)*((VROTL(1)**2) + (VROTL(3)**2))
  41.       ROME(3,3) = (-1.)*((VROTL(1)**2) + (VROTL(2)**2))
  42.       ROME(1,2) = VROTL(1)*VROTL(2)
  43.       ROME(1,3) = VROTL(1)*VROTL(3)
  44.       ROME(2,3) = VROTL(2)*VROTL(3)
  45.       ROME(2,1) = ROME(1,2)
  46.       ROME(3,1) = ROME(1,3)
  47.       ROME(3,2) = ROME(2,3)
  48. C
  49.       ROME(4,4) = ROME(1,1)
  50.       ROME(4,5) = ROME(1,2)
  51.       ROME(4,6) = ROME(1,3)
  52.       ROME(5,4) = ROME(2,1)
  53.       ROME(5,5) = ROME(2,2)
  54.       ROME(5,6) = ROME(2,3)
  55.       ROME(6,4) = ROME(3,1)
  56.       ROME(6,5) = ROME(3,2)
  57.       ROME(6,6) = ROME(3,3)
  58. C
  59. C
  60. C  REMPLISSAGE DE LA MATRICE
  61. C
  62.       DO 201 I=1,3
  63.         DO 201 J=1,3
  64.         RP(I,J) =  DDHOOK( I, J)*ROME(I,J)*XL/3.D0
  65.         RP(I+6,J) =DDHOOK( I, J)*ROME(I,J)*XL/6.D0
  66.         RP(I,J+6) =DDHOOK( I, J)*ROME(I,J)*XL/6.D0
  67.         RP(I+6,J+6) =DDHOOK( I, J)*ROME(I,J)*XL/3.D0
  68.  201  CONTINUE
  69. C
  70.       RETURN
  71.       END
  72.  
  73.  
  74.  
  75.  
  76.  

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