Esope source files

   1 : C TIMALO    SOURCE    CHAT      05/01/13    03:39:56     5004           
   2 :       SUBROUTINE TIMALO(RP,XL,RHO,SD,XINX,XINY,XINZ)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C    ROUTINE DE CALCUL DE LA MATRICE DE MASSE ELEMENTAIRE D UNE
   6 : C    POUTRE DROITE EN AXES LOCAUX
   7 : C    D APRES TEDEL
   8 : C    DANS CETTE VERSION,ON PREND EN COMPTE LES INERTIES DE ROTATION
   9 : C    MAIS PAS LA DEFORMABILITE A L EFFORT TRANCHANT
  10 : C
  11 : C  ENTREE
  12 : C    XL      LONGUEUR
  13 : C    SD      SECTION DROITE
  14 : C    XINX    INERTIE IX
  15 : C    XINY    INERTIE IY
  16 : C    XINZ    INERTIE IZ
  17 : C    RHO     MASSE VOLUMIQUE
  18 : C  SORTIE
  19 : C    RP      MATRICE DE MASSE ELEMENTAIRE 12X12
  20 : C
  21 : C    PP 8/9/92
  22 : C
  23 : C-----------------------------------------------------------------------
  24 :       IMPLICIT INTEGER(I-N)
  25 :       IMPLICIT REAL*8(A-H,O-Z)
  26 :       DIMENSION RP(12,12)
  27 : C
  28 : C  INITIALISATIONS
  29 : C
  30 :       DO 120 I=1,144
  31 :   120 RP(I,1)=0.D00
  32 : C
  33 :       RHOSL=RHO*SD*XL
  34 :       RSP1= RHO*SD*  XL*0.333333333333333333
  35 :       RSP2= RHO*SD*  XL*0.166666666666666666
  36 :       TOR1= RHO*XINX*XL*0.333333333333333333
  37 :       TOR2= RHO*XINX*XL*0.166666666666666666
  38 :       FLY1= RHO*XINY*XL*0.333333333333333333
  39 :       FLY2= RHO*XINY*XL*0.166666666666666666
  40 :       FLZ1= RHO*XINZ*XL*0.333333333333333333
  41 :       FLZ2= RHO*XINZ*XL*0.166666666666666666
  42 : C
  43 : C  REMPLISSAGE DE LA MATRICE
  44 : C  FX1
  45 :       RP(1,1)=RSP1
  46 :       RP(1,7)=RSP2
  47 : C  FY1
  48 :       RP(2,2)=RSP1
  49 :       RP(2,8)=RSP2
  50 : C  FZ1
  51 :       RP(3,3)=RSP1
  52 :       RP(3,9)=RSP2
  53 : C  MX1
  54 :       RP(4,4)= TOR1
  55 :       RP(4,10)= TOR2
  56 : C  MY1
  57 :       RP(5,5)=FLY1
  58 :       RP(5,11)=FLY2
  59 : C  MZ1
  60 :       RP(6,6)=FLZ1
  61 :       RP(6,12)=FLZ2
  62 : C  FX2
  63 :       RP(7,7)=RSP1
  64 : C  FY2
  65 :       RP(8,8)=RSP1
  66 : C  FZ2
  67 :       RP(9,9)=RSP1
  68 : C  MX2
  69 :       RP(10,10)=TOR1
  70 : C  MY2
  71 :       RP(11,11)=FLY1
  72 : C  MZ2
  73 :       RP(12,12)=FLZ1
  74 : C
  75 :       DO 9 I=1,12
  76 :       DO 9 J=I,12
  77 :   9   RP(J,I)= RP(I,J)
  78 :       RETURN
  79 :       END
  80 :