Télécharger barkce.eso

Retour à la liste

Numérotation des lignes :

  1. C BARKCE SOURCE KK2000 05/05/25 21:15:06 5088
  2. SUBROUTINE BARKCE(REL,LRE,CARAC,COORD,VROT,WORK,KERRE)
  3. C-----------------------------------------------------------------------
  4. C ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR CENTRIFUGE ELEMENTAIRE
  5. C DE L ELEMENT DE BARRE (REPERE TOURNANT)
  6. C
  7. C ENTREE
  8. C LRE DIMENSION DU TABLEAU REL
  9. C CARAC CARACTERISTIQUE MATERIAU
  10. C COORD COORDONNEES DES 2 NOEUDS
  11. C VROT VECTEUR VITESSE DE ROTATION
  12. C
  13. C SORTIE
  14. C REL MATRICE DE COUPLAGE ELEMENTAIRE
  15. C
  16. C Didier COMBESCURE mars 2003
  17. C
  18. C-----------------------------------------------------------------------
  19. IMPLICIT INTEGER(I-N)
  20. IMPLICIT REAL*8(A-H,O-Z)
  21. -INC CCREEL
  22. DIMENSION REL(LRE,*),CARAC(*),COORD(3,*),V(3),WORK(*),VROT(*)
  23. DIMENSION ROME(3,3)
  24. C
  25. SEC0=CARAC(1)
  26. RHO=CARAC(2)
  27. C
  28. DO 1 I=1,36
  29. 1 REL(I,1)=0.D00
  30. C
  31. V(3)=0.
  32. DO 2 I=1,3
  33. 2 V(I)=COORD(I,2)-COORD(I,1)
  34. XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3))
  35. C WRITE(*,*) XL
  36. XL=SQRT(XL)
  37. C
  38. FAC=SEC0*RHO*XL/6.
  39. FAC2=FAC*2.
  40. C
  41. ROME(1,1) = (-1.)*((VROT(2)**2) + (VROT(3)**2))
  42. ROME(2,2) = (-1.)*((VROT(1)**2) + (VROT(3)**2))
  43. ROME(3,3) = (-1.)*((VROT(1)**2) + (VROT(2)**2))
  44. ROME(1,2) = VROT(1)*VROT(2)
  45. ROME(1,3) = VROT(1)*VROT(3)
  46. ROME(2,3) = VROT(2)*VROT(3)
  47. ROME(2,1) = ROME(1,2)
  48. ROME(3,1) = ROME(1,3)
  49. ROME(3,2) = ROME(2,3)
  50. C
  51. C REMPLISSAGE DE LA MATRICE
  52. C FX1
  53. REL(1,1)=FAC2*ROME(1,1)
  54. REL(1,2)=FAC2*ROME(1,2)
  55. REL(1,3)=FAC2*ROME(1,3)
  56. REL(1,4)=FAC*ROME(1,1)
  57. REL(1,5)=FAC*ROME(1,2)
  58. REL(1,6)=FAC*ROME(1,3)
  59. C FY1
  60. REL(2,1)=FAC2*ROME(2,1)
  61. REL(2,2)=FAC2*ROME(2,2)
  62. REL(2,3)=FAC2*ROME(2,3)
  63. REL(2,4)=FAC*ROME(2,1)
  64. REL(2,5)=FAC*ROME(2,2)
  65. REL(2,6)=FAC*ROME(2,3)
  66. C FZ1
  67. REL(3,1)=FAC2*ROME(3,1)
  68. REL(3,2)=FAC2*ROME(3,2)
  69. REL(3,3)=FAC2*ROME(3,3)
  70. REL(3,4)=FAC*ROME(3,1)
  71. REL(3,5)=FAC*ROME(3,2)
  72. REL(3,6)=FAC*ROME(3,3)
  73. C FX2
  74. REL(4,1)=FAC*ROME(1,1)
  75. REL(4,2)=FAC*ROME(1,2)
  76. REL(4,3)=FAC*ROME(1,3)
  77. REL(4,4)=FAC2*ROME(1,1)
  78. REL(4,5)=FAC2*ROME(1,2)
  79. REL(4,6)=FAC2*ROME(1,3)
  80. C FY2
  81. REL(5,1)=FAC*ROME(2,1)
  82. REL(5,2)=FAC*ROME(2,2)
  83. REL(5,3)=FAC*ROME(2,3)
  84. REL(5,4)=FAC2*ROME(2,1)
  85. REL(5,5)=FAC2*ROME(2,2)
  86. REL(5,6)=FAC2*ROME(2,3)
  87. C FZ2
  88. REL(6,1)=FAC2*ROME(3,1)
  89. REL(6,2)=FAC2*ROME(3,2)
  90. REL(6,3)=FAC2*ROME(3,3)
  91. REL(6,4)=FAC*ROME(3,1)
  92. REL(6,5)=FAC*ROME(3,2)
  93. REL(6,6)=FAC*ROME(3,3)
  94. C
  95. KERRE = 0
  96. C
  97. RETURN
  98. END
  99.  
  100.  
  101.  
  102.  

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