Esope source files

   1 : C TIFCLO    SOURCE    KK2000    05/05/25    21:15:27     5088           
   2 :       SUBROUTINE TIFCLO(RP,XL,LHOOK,DDHOOK,VROTL)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C    ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR CENTRIFUGE ELEMENTAIRE 
   6 : C    D UNE POUTRE DE TIMOSHENKO DROITE MULTI-FIBRE EN AXES LOCAUX
   7 : C  
   8 : C
   9 : C  ENTREE
  10 : C    XL      LONGUEUR
  11 : C    DDHOOK(LHOOK,LHOOK) MATRICE DE "HOOK"
  12 : C
  13 : C  SORTIE
  14 : C    RP      MATRICE DE RAIDEUR ELEMENTAIRE 12X12
  15 : C
  16 : C-----------------------------------------------------------------------
  17 : C     Didier Combescure 2003
  18 : C-----------------------------------------------------------------------
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)
  21 : -INC CCREEL
  22 :       DIMENSION RP(12,12),VROTL(3),ROME(6,6)
  23 :       DIMENSION DDHOOK(LHOOK,LHOOK)
  24 : C
  25 : C  INITIALISATIONS
  26 : C
  27 :       DO 120 I=1,12
  28 :       DO 120 J=1,12
  29 :   120 RP(I,J)=0.D00
  30 : C
  31 :       RHOSD = DDHOOK( 1, 1)
  32 :       FAC=RHOSD*XL/6.D0                                                             
  33 :       FAC2=FAC*2.D0 
  34 : C
  35 :       DO 121 I=1,6
  36 :       DO 121 J=1,6
  37 :   121 ROME(I,J)=0.D00
  38 : C
  39 :       ROME(1,1) = (-1.)*((VROTL(2)**2) + (VROTL(3)**2))
  40 :       ROME(2,2) = (-1.)*((VROTL(1)**2) + (VROTL(3)**2))
  41 :       ROME(3,3) = (-1.)*((VROTL(1)**2) + (VROTL(2)**2))
  42 :       ROME(1,2) = VROTL(1)*VROTL(2)
  43 :       ROME(1,3) = VROTL(1)*VROTL(3)
  44 :       ROME(2,3) = VROTL(2)*VROTL(3)
  45 :       ROME(2,1) = ROME(1,2) 
  46 :       ROME(3,1) = ROME(1,3) 
  47 :       ROME(3,2) = ROME(2,3) 
  48 : C
  49 :       ROME(4,4) = ROME(1,1) 
  50 :       ROME(4,5) = ROME(1,2) 
  51 :       ROME(4,6) = ROME(1,3) 
  52 :       ROME(5,4) = ROME(2,1) 
  53 :       ROME(5,5) = ROME(2,2) 
  54 :       ROME(5,6) = ROME(2,3) 
  55 :       ROME(6,4) = ROME(3,1) 
  56 :       ROME(6,5) = ROME(3,2) 
  57 :       ROME(6,6) = ROME(3,3) 
  58 : C
  59 : C
  60 : C  REMPLISSAGE DE LA MATRICE
  61 : C  
  62 :       DO 201 I=1,3
  63 :         DO 201 J=1,3 
  64 :         RP(I,J) =  DDHOOK( I, J)*ROME(I,J)*XL/3.D0
  65 :         RP(I+6,J) =DDHOOK( I, J)*ROME(I,J)*XL/6.D0
  66 :         RP(I,J+6) =DDHOOK( I, J)*ROME(I,J)*XL/6.D0
  67 :         RP(I+6,J+6) =DDHOOK( I, J)*ROME(I,J)*XL/3.D0
  68 :  201  CONTINUE
  69 : C
  70 :       RETURN
  71 :       END
  72 :  
  73 :  
  74 :  
  75 :