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    CB215821  26/08/24    21:18:37     12622          
  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 121 J=1,12
  32.         DO 120 I=1,12
  33.           RP(I,J)=0.D00
  34.   120   CONTINUE
  35.   121 CONTINUE
  36. C
  37.       RHOSD = RHO*SD
  38.       FAC=RHOSD*XL/6.D0
  39.       FAC2=FAC*2.D0
  40.       TOR1= RHO*XINX*XL/3.D0
  41.       TOR2= RHO*XINX*XL/6.D0
  42.       FLZ1= RHO*XINY*XL/3.D0
  43.       FLZ2= RHO*XINY*XL/6.D0
  44.       FLY1= RHO*XINZ*XL/3.D0
  45.       FLY2= RHO*XINZ*XL/6.D0
  46. C
  47. C
  48.       ROME(1,1) = (-1.)*((VROTL(2)**2) + (VROTL(3)**2))
  49.       ROME(2,2) = (-1.)*((VROTL(1)**2) + (VROTL(3)**2))
  50.       ROME(3,3) = (-1.)*((VROTL(1)**2) + (VROTL(2)**2))
  51.       ROME(1,2) = VROTL(1)*VROTL(2)
  52.       ROME(1,3) = VROTL(1)*VROTL(3)
  53.       ROME(2,3) = VROTL(2)*VROTL(3)
  54.       ROME(2,1) = ROME(1,2)
  55.       ROME(3,1) = ROME(1,3)
  56.       ROME(3,2) = ROME(2,3)
  57. C
  58. C  REMPLISSAGE DE LA MATRICE
  59. C  FX1
  60.       RP(1,1)=FAC2*ROME(1,1)
  61.       RP(1,2)=FAC2*ROME(1,2)
  62.       RP(1,3)=FAC2*ROME(1,3)
  63.       RP(1,7)=FAC*ROME(1,1)
  64.       RP(1,8)=FAC*ROME(1,2)
  65.       RP(1,9)=FAC*ROME(1,3)
  66. C  FY1
  67.       RP(2,1)=FAC2*ROME(2,1)
  68.       RP(2,2)=FAC2*ROME(2,2)
  69.       RP(2,3)=FAC2*ROME(2,3)
  70.       RP(2,7)=FAC*ROME(2,1)
  71.       RP(2,8)=FAC*ROME(2,2)
  72.       RP(2,9)=FAC*ROME(2,3)
  73. C  FZ1
  74.       RP(3,1)=FAC2*ROME(3,1)
  75.       RP(3,2)=FAC2*ROME(3,2)
  76.       RP(3,3)=FAC2*ROME(3,3)
  77.       RP(3,7)=FAC*ROME(3,1)
  78.       RP(3,5)=FAC*ROME(3,2)
  79.       RP(3,9)=FAC*ROME(3,3)
  80. C  FX2
  81.       RP(7,1)=FAC*ROME(1,1)
  82.       RP(7,2)=FAC*ROME(1,2)
  83.       RP(7,3)=FAC*ROME(1,3)
  84.       RP(7,7)=FAC2*ROME(1,1)
  85.       RP(7,8)=FAC2*ROME(1,2)
  86.       RP(7,9)=FAC2*ROME(1,3)
  87. C  FY2
  88.       RP(8,1)=FAC*ROME(2,1)
  89.       RP(8,2)=FAC*ROME(2,2)
  90.       RP(8,3)=FAC*ROME(2,3)
  91.       RP(8,7)=FAC2*ROME(2,1)
  92.       RP(8,8)=FAC2*ROME(2,2)
  93.       RP(8,9)=FAC2*ROME(2,3)
  94. C  FZ2
  95.       RP(9,1)=FAC2*ROME(3,1)
  96.       RP(9,2)=FAC2*ROME(3,2)
  97.       RP(9,3)=FAC2*ROME(3,3)
  98.       RP(9,7)=FAC*ROME(3,1)
  99.       RP(9,8)=FAC*ROME(3,2)
  100.       RP(9,9)=FAC*ROME(3,3)
  101. C
  102. C On rajoute le moment
  103. C
  104. C  MX1
  105. C      RP(4,4)=TOR1*ROME(1,1)
  106. C      RP(4,5)=TOR1*ROME(1,2)
  107. C      RP(4,6)=TOR1*ROME(1,3)
  108. C      RP(4,10)=TOR2*ROME(1,1)
  109. C      RP(4,11)=TOR2*ROME(1,2)
  110. C      RP(4,12)=TOR2*ROME(1,3)
  111. C  MY1
  112. C      RP(5,4)=FLY1*ROME(2,1)
  113. C      RP(5,5)=FLY1*ROME(2,2)
  114. C      RP(5,6)=FLY1*ROME(2,3)
  115. C      RP(5,10)=FLY2*ROME(2,1)
  116. C      RP(5,11)=FLY2*ROME(2,2)
  117. C      RP(5,12)=FLY2*ROME(2,3)
  118. C  MZ1
  119. C      RP(6,4)=FLZ1*ROME(3,1)
  120. C      RP(6,5)=FLZ1*ROME(3,2)
  121. C      RP(6,6)=FLZ1*ROME(3,3)
  122. C      RP(6,10)=FLZ2*ROME(3,1)
  123. C      RP(6,11)=FLZ2*ROME(3,2)
  124. C      RP(6,12)=FLZ2*ROME(3,3)
  125. C  MX2
  126. C      RP(10,10)=TOR1*ROME(1,1)
  127. C      RP(10,11)=TOR1*ROME(1,2)
  128. C      RP(10,12)=TOR1*ROME(1,3)
  129. C      RP(10,4)=TOR2*ROME(1,1)
  130. C      RP(10,5)=TOR2*ROME(1,2)
  131. C      RP(10,6)=TOR2*ROME(1,3)
  132. C  MY2
  133. C      RP(11,4)=FLY2*ROME(2,1)
  134. C      RP(11,5)=FLY2*ROME(2,2)
  135. C      RP(11,6)=FLY2*ROME(2,3)
  136. C      RP(11,10)=FLY1*ROME(2,1)
  137. C      RP(11,11)=FLY1*ROME(2,2)
  138. C      RP(11,12)=FLY1*ROME(2,3)
  139. C  MZ2
  140. C      RP(12,4)=FLZ2*ROME(3,1)
  141. C      RP(12,5)=FLZ2*ROME(3,2)
  142. C      RP(12,6)=FLZ2*ROME(3,3)
  143. C      RP(12,10)=FLZ1*ROME(3,1)
  144. C      RP(12,11)=FLZ1*ROME(3,2)
  145. C      RP(12,12)=FLZ1*ROME(3,3)
  146. C
  147. C
  148.       RETURN
  149.       END
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  

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