Esope source files

   1 : C TIMAL2    SOURCE    CHAT      05/01/13    03:39:51     5004           
   2 : C TIMAL2    SOURCE    AM        93/01/12    21:35:51     781            
   3 :       SUBROUTINE TIMAL2(RP,XL,RHO,SD,XINZ)
   4 : C-----------------------------------------------------------------------
   5 : C
   6 : C    ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D UNE
   7 : C    POUTRE DROITE EN AXES LOCAUX
   8 : C    D APRES TEDEL
   9 : C    DANS CETTE VERSION,ON PREND EN COMPTE LES INERTIES DE ROTATION
  10 : C    MAIS PAS LA DEFORMABILITE A L EFFORT TRANCHANT
  11 : C
  12 : C  ENTREE
  13 : C    XL      LONGUEUR
  14 : C    SD      SECTION DROITE
  15 : C    XINZ    INERTIE IZ
  16 : C    RHO     MASSE VOLUMIQUE
  17 : C  SORTIE
  18 : C    RP      MATRICE DE MASSE ELEMENTAIRE 6X6
  19 : C
  20 : C    D'APRES TIMALO DC 98
  21 : C
  22 : C-----------------------------------------------------------------------
  23 :       IMPLICIT INTEGER(I-N)
  24 :       IMPLICIT REAL*8(A-H,O-Z)
  25 :       DIMENSION RP(6,6)
  26 : C
  27 : C  INITIALISATIONS
  28 : C
  29 :       DO 120 I=1,36
  30 :   120 RP(I,1)=0.D00
  31 : C
  32 :       RHOSL=RHO*SD*XL
  33 :       RSP1= RHO*SD*  XL*0.333333333333333333
  34 :       RSP2= RHO*SD*  XL*0.166666666666666666
  35 :       FLZ1= RHO*XINZ*XL*0.333333333333333333
  36 :       FLZ2= RHO*XINZ*XL*0.166666666666666666
  37 : C
  38 : C  REMPLISSAGE DE LA MATRICE
  39 : C  FX1
  40 :       RP(1,1)=RSP1
  41 :       RP(1,4)=RSP2
  42 : C  FY1
  43 :       RP(2,2)=RSP1
  44 :       RP(2,5)=RSP2
  45 : C  MZ1
  46 :       RP(3,3)=FLZ1
  47 :       RP(3,6)=FLZ2
  48 : C  FX2
  49 :       RP(4,4)=RSP1
  50 : C  FY2
  51 :       RP(5,5)=RSP1
  52 : C  MZ2
  53 :       RP(6,6)=FLZ1
  54 : C
  55 :       DO 9 I=1,6
  56 :       DO 9 J=I,6
  57 :   9   RP(J,I)= RP(I,J)
  58 :       RETURN
  59 :       END
  60 :  
  61 :