C POKCLO    SOURCE    CB215821  26/08/24    21:17:43     12622          
      SUBROUTINE POKCLO(RP,XL,RHO,SD,XINX,XINY,XINZ,VROTL)
C-----------------------------------------------------------------------
C
C    ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR
C   CENTRIFUGE ELEMENTAIRE D UNE
C    POUTRE DROITE EN AXES LOCAUX
C
C
C  ENTREE
C    XL      LONGUEUR
C    SD      SECTION DROITE
C    XINX    INERTIE IX
C    XINY    INERTIE IY
C    XINZ    INERTIE IZ
C    RHO     MASSE VOLUMIQUE
C  SORTIE
C    RP      MATRICE DE RAIDEUR CENTRIFUGE ELEMENTAIRE 12X12
C
C                         Didier COMBESCURE mars 2003
C
C-----------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
-INC CCREEL
      DIMENSION RP(12,12),VROTL(3),ROME(3,3)
C
C  INITIALISATIONS
C
C
      DO 121 J=1,12
      DO 120 I=1,12
      RP(I,J)=0.D00
  120 CONTINUE
  121 CONTINUE
C
      RHOSD = RHO*SD
      FAC=RHOSD*XL/6.D0
      FAC2=FAC*2.D0
      TOR1= RHO*XINX*XL/3.D0
      TOR2= RHO*XINX*XL/6.D0
      FLZ1= RHO*XINY*XL/3.D0
      FLZ2= RHO*XINY*XL/6.D0
      FLY1= RHO*XINZ*XL/3.D0
      FLY2= RHO*XINZ*XL/6.D0
C
C
      ROME(1,1) = (-1.)*((VROTL(2)**2) + (VROTL(3)**2))
      ROME(2,2) = (-1.)*((VROTL(1)**2) + (VROTL(3)**2))
      ROME(3,3) = (-1.)*((VROTL(1)**2) + (VROTL(2)**2))
      ROME(1,2) = VROTL(1)*VROTL(2)
      ROME(1,3) = VROTL(1)*VROTL(3)
      ROME(2,3) = VROTL(2)*VROTL(3)
      ROME(2,1) = ROME(1,2)
      ROME(3,1) = ROME(1,3)
      ROME(3,2) = ROME(2,3)
C
C  REMPLISSAGE DE LA MATRICE
C  FX1
      RP(1,1)=FAC2*ROME(1,1)
      RP(1,2)=FAC2*ROME(1,2)
      RP(1,3)=FAC2*ROME(1,3)
      RP(1,7)=FAC*ROME(1,1)
      RP(1,8)=FAC*ROME(1,2)
      RP(1,9)=FAC*ROME(1,3)
C  FY1
      RP(2,1)=FAC2*ROME(2,1)
      RP(2,2)=FAC2*ROME(2,2)
      RP(2,3)=FAC2*ROME(2,3)
      RP(2,7)=FAC*ROME(2,1)
      RP(2,8)=FAC*ROME(2,2)
      RP(2,9)=FAC*ROME(2,3)
C  FZ1
      RP(3,1)=FAC2*ROME(3,1)
      RP(3,2)=FAC2*ROME(3,2)
      RP(3,3)=FAC2*ROME(3,3)
      RP(3,7)=FAC*ROME(3,1)
      RP(3,5)=FAC*ROME(3,2)
      RP(3,9)=FAC*ROME(3,3)
C  FX2
      RP(7,1)=FAC*ROME(1,1)
      RP(7,2)=FAC*ROME(1,2)
      RP(7,3)=FAC*ROME(1,3)
      RP(7,7)=FAC2*ROME(1,1)
      RP(7,8)=FAC2*ROME(1,2)
      RP(7,9)=FAC2*ROME(1,3)
C  FY2
      RP(8,1)=FAC*ROME(2,1)
      RP(8,2)=FAC*ROME(2,2)
      RP(8,3)=FAC*ROME(2,3)
      RP(8,7)=FAC2*ROME(2,1)
      RP(8,8)=FAC2*ROME(2,2)
      RP(8,9)=FAC2*ROME(2,3)
C  FZ2
      RP(9,1)=FAC2*ROME(3,1)
      RP(9,2)=FAC2*ROME(3,2)
      RP(9,3)=FAC2*ROME(3,3)
      RP(9,7)=FAC*ROME(3,1)
      RP(9,8)=FAC*ROME(3,2)
      RP(9,9)=FAC*ROME(3,3)
C
C
      RETURN
      END