Télécharger tifclo.eso

Retour à la liste

Numérotation des lignes :

tifclo
  1. C TIFCLO SOURCE CB215821 26/08/24 21:18:33 12622
  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 202 I=1,12
  28. DO 120 J=1,12
  29. RP(I,J)=0.D00
  30. 120 CONTINUE
  31. 202 CONTINUE
  32. C
  33. RHOSD = DDHOOK( 1, 1)
  34. FAC=RHOSD*XL/6.D0
  35. FAC2=FAC*2.D0
  36. C
  37. DO 203 I=1,6
  38. DO 121 J=1,6
  39. ROME(I,J)=0.D00
  40. 121 CONTINUE
  41. 203 CONTINUE
  42. C
  43. ROME(1,1) = (-1.)*((VROTL(2)**2) + (VROTL(3)**2))
  44. ROME(2,2) = (-1.)*((VROTL(1)**2) + (VROTL(3)**2))
  45. ROME(3,3) = (-1.)*((VROTL(1)**2) + (VROTL(2)**2))
  46. ROME(1,2) = VROTL(1)*VROTL(2)
  47. ROME(1,3) = VROTL(1)*VROTL(3)
  48. ROME(2,3) = VROTL(2)*VROTL(3)
  49. ROME(2,1) = ROME(1,2)
  50. ROME(3,1) = ROME(1,3)
  51. ROME(3,2) = ROME(2,3)
  52. C
  53. ROME(4,4) = ROME(1,1)
  54. ROME(4,5) = ROME(1,2)
  55. ROME(4,6) = ROME(1,3)
  56. ROME(5,4) = ROME(2,1)
  57. ROME(5,5) = ROME(2,2)
  58. ROME(5,6) = ROME(2,3)
  59. ROME(6,4) = ROME(3,1)
  60. ROME(6,5) = ROME(3,2)
  61. ROME(6,6) = ROME(3,3)
  62. C
  63. C
  64. C REMPLISSAGE DE LA MATRICE
  65. C
  66. DO 204 I=1,3
  67. DO 201 J=1,3
  68. RP(I,J) = DDHOOK( I, J)*ROME(I,J)*XL/3.D0
  69. RP(I+6,J) =DDHOOK( I, J)*ROME(I,J)*XL/6.D0
  70. RP(I,J+6) =DDHOOK( I, J)*ROME(I,J)*XL/6.D0
  71. RP(I+6,J+6) =DDHOOK( I, J)*ROME(I,J)*XL/3.D0
  72. 201 CONTINUE
  73. 204 CONTINUE
  74. C
  75. RETURN
  76. END
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  

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