Télécharger la plaquette Cast3M


fr en

Télécharger pokclo.eso

Retour à la liste

Numérotation des lignes : 

pokclo

  1. C POKCLO    SOURCE    PV        18/06/18    21:15:22     9860           
  2.       SUBROUTINE POKCLO(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 DROITE 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
  101.       RETURN
  102.       END
  103.  
  104.  
  105.  
  106.  
  107.  

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