Esope source files

   1 : C TIMGF1    SOURCE    AM        10/09/20    21:15:24     6756           
   2 :       SUBROUTINE TIMGF1(COORD,DEPLA,VECT,GRAD,WORK)
   3 : C-----------------------------------------------------------------------
   4 : C  ROUTINE DE CALCUL DES GRADIENTS DE FLEXION POUR LA POUTRE DROITE
   5 : C  DE TIMOSHENKO SIMPLE OU MULTI-FIBRE
   6 : C
   7 : C  ENTREE
   8 : C          DEPLA   DEPLACEMENTS
   9 : C          VECT    ORIENTATION DE L'AXE OY LOCAL
  10 : C          COORD   COORDONNEES DES 2 NOEUDS
  11 : C
  12 : C  SORTIE
  13 : C          GRAD    GRADIENTS DE FLEXION  AU CENTRE 3 PAR NOEUD
  14 : C  TRAVAIL
  15 : C          WORK(300)
  16 : C-----------------------------------------------------------------------
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)
  19 :       DIMENSION WORK(*)
  20 :       DIMENSION GRAD(*),COORD(3,*),VECT(3),P(3,3),
  21 :      .          XE(2),YE(2),ZE(2),DEPLA(*)
  22 : C
  23 :       DO I=1,2
  24 :         XE(I)=COORD(1,I)
  25 :         YE(I)=COORD(2,I)
  26 :         ZE(I)=COORD(3,I)
  27 :       ENDDO
  28 :       XL=(XE(2)-XE(1))**2+(YE(2)-YE(1))**2+(ZE(2)-ZE(1))**2
  29 :       XL=SQRT(XL)
  30 : C
  31 : C  CALCUL DE LA MATRICE DE PASSAGE
  32 : C
  32 :       CALL POUPAS(XE,YE,ZE,VECT,P,KERRE)
  33 : C
  34 : C  PASSAGE DES DEPLACEMENTS EN AXES LOCAUX
  35 : C
  35 :       CALL POUVEC(DEPLA,WORK,P,1)
  36 : C
  37 : C  CALCUL DES GRADIENTS DE FLEXION AU CENTRE
  38 : C
  39 :       GRAD(1)=(WORK(10)-WORK(4))/XL
  40 :       GRAD(2)=(WORK(11)-WORK(5))/XL
  41 :       GRAD(3)=(WORK(12)-WORK(6))/XL
  42 :       RETURN
  43 :       END
  44 :  
  45 :