Télécharger la plaquette Cast3M


fr en

Télécharger tikclo.eso

Retour à la liste

Numérotation des lignes : 

tikclo

  1. C TIKCLO    SOURCE    PV        18/06/18    21:15:38     9860           
  2.       SUBROUTINE TIKCLO(RP,XL,RHO,SD,XINX,XINY,XINZ,VROTL)
  3. C-----------------------------------------------------------------------
  4. C
  5. C    ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR
  6. C   CENTRIFUGE ELEMENTAIRE D UNE
  7. C    POUTRE DE TIMOSHENKO EN AXES LOCAUX
  8. C
  9. C
  10. C  ENTREE
  11. C    XL      LONGUEUR
  12. C    SD      SECTION DROITE
  13. C    XINX    INERTIE IX
  14. C    XINY    INERTIE IY
  15. C    XINZ    INERTIE IZ
  16. C    RHO     MASSE VOLUMIQUE
  17. C  SORTIE
  18. C    RP      MATRICE DE RAIDEUR CENTRIFUGE ELEMENTAIRE 12X12
  19. C
  20. C                         Didier COMBESCURE mars 2003
  21. C
  22. C-----------------------------------------------------------------------
  23.       IMPLICIT INTEGER(I-N)
  24.       IMPLICIT REAL*8(A-H,O-Z)
  25. -INC CCREEL
  26.       DIMENSION RP(12,12),VROTL(3),ROME(3,3)
  27. C
  28. C  INITIALISATIONS
  29. C
  30. C
  31.       DO 120 J=1,12
  32.       DO 120 I=1,12
  33.   120 RP(I,J)=0.D00
  34. C
  35.       RHOSD = RHO*SD
  36.       FAC=RHOSD*XL/6.D0
  37.       FAC2=FAC*2.D0
  38.       TOR1= RHO*XINX*XL/3.D0
  39.       TOR2= RHO*XINX*XL/6.D0
  40.       FLZ1= RHO*XINY*XL/3.D0
  41.       FLZ2= RHO*XINY*XL/6.D0
  42.       FLY1= RHO*XINZ*XL/3.D0
  43.       FLY2= RHO*XINZ*XL/6.D0
  44. C
  45. C
  46.       ROME(1,1) = (-1.)*((VROTL(2)**2) + (VROTL(3)**2))
  47.       ROME(2,2) = (-1.)*((VROTL(1)**2) + (VROTL(3)**2))
  48.       ROME(3,3) = (-1.)*((VROTL(1)**2) + (VROTL(2)**2))
  49.       ROME(1,2) = VROTL(1)*VROTL(2)
  50.       ROME(1,3) = VROTL(1)*VROTL(3)
  51.       ROME(2,3) = VROTL(2)*VROTL(3)
  52.       ROME(2,1) = ROME(1,2)
  53.       ROME(3,1) = ROME(1,3)
  54.       ROME(3,2) = ROME(2,3)
  55. C
  56. C  REMPLISSAGE DE LA MATRICE
  57. C  FX1
  58.       RP(1,1)=FAC2*ROME(1,1)
  59.       RP(1,2)=FAC2*ROME(1,2)
  60.       RP(1,3)=FAC2*ROME(1,3)
  61.       RP(1,7)=FAC*ROME(1,1)
  62.       RP(1,8)=FAC*ROME(1,2)
  63.       RP(1,9)=FAC*ROME(1,3)
  64. C  FY1
  65.       RP(2,1)=FAC2*ROME(2,1)
  66.       RP(2,2)=FAC2*ROME(2,2)
  67.       RP(2,3)=FAC2*ROME(2,3)
  68.       RP(2,7)=FAC*ROME(2,1)
  69.       RP(2,8)=FAC*ROME(2,2)
  70.       RP(2,9)=FAC*ROME(2,3)
  71. C  FZ1
  72.       RP(3,1)=FAC2*ROME(3,1)
  73.       RP(3,2)=FAC2*ROME(3,2)
  74.       RP(3,3)=FAC2*ROME(3,3)
  75.       RP(3,7)=FAC*ROME(3,1)
  76.       RP(3,5)=FAC*ROME(3,2)
  77.       RP(3,9)=FAC*ROME(3,3)
  78. C  FX2
  79.       RP(7,1)=FAC*ROME(1,1)
  80.       RP(7,2)=FAC*ROME(1,2)
  81.       RP(7,3)=FAC*ROME(1,3)
  82.       RP(7,7)=FAC2*ROME(1,1)
  83.       RP(7,8)=FAC2*ROME(1,2)
  84.       RP(7,9)=FAC2*ROME(1,3)
  85. C  FY2
  86.       RP(8,1)=FAC*ROME(2,1)
  87.       RP(8,2)=FAC*ROME(2,2)
  88.       RP(8,3)=FAC*ROME(2,3)
  89.       RP(8,7)=FAC2*ROME(2,1)
  90.       RP(8,8)=FAC2*ROME(2,2)
  91.       RP(8,9)=FAC2*ROME(2,3)
  92. C  FZ2
  93.       RP(9,1)=FAC2*ROME(3,1)
  94.       RP(9,2)=FAC2*ROME(3,2)
  95.       RP(9,3)=FAC2*ROME(3,3)
  96.       RP(9,7)=FAC*ROME(3,1)
  97.       RP(9,8)=FAC*ROME(3,2)
  98.       RP(9,9)=FAC*ROME(3,3)
  99. C
  100. C On rajoute le moment
  101. C
  102. C  MX1
  103. C      RP(4,4)=TOR1*ROME(1,1)
  104. C      RP(4,5)=TOR1*ROME(1,2)
  105. C      RP(4,6)=TOR1*ROME(1,3)
  106. C      RP(4,10)=TOR2*ROME(1,1)
  107. C      RP(4,11)=TOR2*ROME(1,2)
  108. C      RP(4,12)=TOR2*ROME(1,3)
  109. C  MY1
  110. C      RP(5,4)=FLY1*ROME(2,1)
  111. C      RP(5,5)=FLY1*ROME(2,2)
  112. C      RP(5,6)=FLY1*ROME(2,3)
  113. C      RP(5,10)=FLY2*ROME(2,1)
  114. C      RP(5,11)=FLY2*ROME(2,2)
  115. C      RP(5,12)=FLY2*ROME(2,3)
  116. C  MZ1
  117. C      RP(6,4)=FLZ1*ROME(3,1)
  118. C      RP(6,5)=FLZ1*ROME(3,2)
  119. C      RP(6,6)=FLZ1*ROME(3,3)
  120. C      RP(6,10)=FLZ2*ROME(3,1)
  121. C      RP(6,11)=FLZ2*ROME(3,2)
  122. C      RP(6,12)=FLZ2*ROME(3,3)
  123. C  MX2
  124. C      RP(10,10)=TOR1*ROME(1,1)
  125. C      RP(10,11)=TOR1*ROME(1,2)
  126. C      RP(10,12)=TOR1*ROME(1,3)
  127. C      RP(10,4)=TOR2*ROME(1,1)
  128. C      RP(10,5)=TOR2*ROME(1,2)
  129. C      RP(10,6)=TOR2*ROME(1,3)
  130. C  MY2
  131. C      RP(11,4)=FLY2*ROME(2,1)
  132. C      RP(11,5)=FLY2*ROME(2,2)
  133. C      RP(11,6)=FLY2*ROME(2,3)
  134. C      RP(11,10)=FLY1*ROME(2,1)
  135. C      RP(11,11)=FLY1*ROME(2,2)
  136. C      RP(11,12)=FLY1*ROME(2,3)
  137. C  MZ2
  138. C      RP(12,4)=FLZ2*ROME(3,1)
  139. C      RP(12,5)=FLZ2*ROME(3,2)
  140. C      RP(12,6)=FLZ2*ROME(3,3)
  141. C      RP(12,10)=FLZ1*ROME(3,1)
  142. C      RP(12,11)=FLZ1*ROME(3,2)
  143. C      RP(12,12)=FLZ1*ROME(3,3)
  144. C
  145. C
  146.       RETURN
  147.       END
  148.  
  149.  
  150.  
  151.  
  152.  

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