Esope source files

   1 : C POKCLO    SOURCE    CHAT      05/01/13    02:16:54     5004           
   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 I=1,144
  32 :   120 RP(I,1)=0.D00
  33 : C
  34 :       RHOSD = RHO*SD
  35 :       FAC=RHOSD*XL/6.D0                                                             
  36 :       FAC2=FAC*2.D0 
  37 :       TOR1= RHO*XINX*XL/3.D0   
  38 :       TOR2= RHO*XINX*XL/6.D0  
  39 :       FLZ1= RHO*XINY*XL/3.D0   
  40 :       FLZ2= RHO*XINY*XL/6.D0   
  41 :       FLY1= RHO*XINZ*XL/3.D0   
  42 :       FLY2= RHO*XINZ*XL/6.D0
  43 : C   
  44 : C
  45 :       ROME(1,1) = (-1.)*((VROTL(2)**2) + (VROTL(3)**2))
  46 :       ROME(2,2) = (-1.)*((VROTL(1)**2) + (VROTL(3)**2))
  47 :       ROME(3,3) = (-1.)*((VROTL(1)**2) + (VROTL(2)**2))
  48 :       ROME(1,2) = VROTL(1)*VROTL(2)
  49 :       ROME(1,3) = VROTL(1)*VROTL(3)
  50 :       ROME(2,3) = VROTL(2)*VROTL(3)
  51 :       ROME(2,1) = ROME(1,2) 
  52 :       ROME(3,1) = ROME(1,3) 
  53 :       ROME(3,2) = ROME(2,3) 
  54 : C
  55 : C  REMPLISSAGE DE LA MATRICE
  56 : C  FX1
  57 :       RP(1,1)=FAC2*ROME(1,1)
  58 :       RP(1,2)=FAC2*ROME(1,2)
  59 :       RP(1,3)=FAC2*ROME(1,3)
  60 :       RP(1,7)=FAC*ROME(1,1)
  61 :       RP(1,8)=FAC*ROME(1,2)
  62 :       RP(1,9)=FAC*ROME(1,3)
  63 : C  FY1
  64 :       RP(2,1)=FAC2*ROME(2,1)
  65 :       RP(2,2)=FAC2*ROME(2,2)
  66 :       RP(2,3)=FAC2*ROME(2,3)
  67 :       RP(2,7)=FAC*ROME(2,1)
  68 :       RP(2,8)=FAC*ROME(2,2)
  69 :       RP(2,9)=FAC*ROME(2,3)
  70 : C  FZ1
  71 :       RP(3,1)=FAC2*ROME(3,1)
  72 :       RP(3,2)=FAC2*ROME(3,2)
  73 :       RP(3,3)=FAC2*ROME(3,3)
  74 :       RP(3,7)=FAC*ROME(3,1)
  75 :       RP(3,5)=FAC*ROME(3,2)
  76 :       RP(3,9)=FAC*ROME(3,3)
  77 : C  FX2
  78 :       RP(7,1)=FAC*ROME(1,1)
  79 :       RP(7,2)=FAC*ROME(1,2)
  80 :       RP(7,3)=FAC*ROME(1,3)
  81 :       RP(7,7)=FAC2*ROME(1,1)
  82 :       RP(7,8)=FAC2*ROME(1,2)
  83 :       RP(7,9)=FAC2*ROME(1,3)
  84 : C  FY2
  85 :       RP(8,1)=FAC*ROME(2,1)
  86 :       RP(8,2)=FAC*ROME(2,2)
  87 :       RP(8,3)=FAC*ROME(2,3)
  88 :       RP(8,7)=FAC2*ROME(2,1)
  89 :       RP(8,8)=FAC2*ROME(2,2)
  90 :       RP(8,9)=FAC2*ROME(2,3)
  91 : C  FZ2
  92 :       RP(9,1)=FAC2*ROME(3,1)
  93 :       RP(9,2)=FAC2*ROME(3,2)
  94 :       RP(9,3)=FAC2*ROME(3,3)
  95 :       RP(9,7)=FAC*ROME(3,1)
  96 :       RP(9,8)=FAC*ROME(3,2)
  97 :       RP(9,9)=FAC*ROME(3,3)
  98 : C
  99 : C 
 100 :       RETURN
 101 :       END
 102 :  
 103 :  
 104 :