Esope source files

   1 : C TIFOLO    SOURCE    CHAT      05/01/13    03:39:23     5004           
   2 :       SUBROUTINE TIFOLO(XL,V,CONT,LHOOK,DDHOOK,IDEUX)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C    ROUTINE DE CALCUL DES CONTRAINTES DE LA POUTRE DROITE
   6 : C    DE TIMOSHENKO MULTI-FIBRE D APRES TEDEL
   7 : C
   8 : C  ENTREE
   9 : C    XL      LONGUEUR DE LA POUTRE DROITE
  10 : C    V       VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE
  11 : C    DDHOOK(LHOOK,LHOOK) MATRICE DE HOOK
  12 : C    IDEUX   INDICE POUR DEFORMATIONS DU SECOND ORDRE
  13 : C  SORTIE
  14 : C    CONT    CONTRAINTES AU MILIEU 6 PAR NOEUD
  15 : C
  16 : C-----------------------------------------------------------------------
  17 : C     Pierre Pegon (ISPRA) Juillet/Aout 1993
  18 : C-----------------------------------------------------------------------
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)
  21 :       DIMENSION CONT(*),V(*),DDHOOK(LHOOK,LHOOK),EPSI(6)
  22 :       DATA P5/0.5D0/
  23 : C
  24 :       RL=1.D00/XL
  25 : C
  26 : C  DEFORMATIONS GENERALISEES
  27 : C
  28 :       EPSI(1) = (V(7)-V(1))*RL
  29 :       EPSI(2) = (V(8)-V(2))*RL - (V(6)+V(12))*P5
  30 :       EPSI(3) = (V(9)-V(3))*RL + (V(5)+V(11))*P5
  31 :       EPSI(4) = (V(10)-V(4))*RL
  32 :       EPSI(5) =-(V(5)-V(11))*RL
  33 :       EPSI(6) =-(V(6)-V(12))*RL
  34 : C
  35 : C  TERMES DU SECOND ORDRE
  36 : C
  37 :       IF(IDEUX.NE.0) THEN
  38 :       ENDIF
  39 : C
  40 : C  CONTRAINTES GENERALISEES
  41 : C
  42 :       DO I=1,6
  43 :         CONT(I)=0.D0
  44 :         DO J=1,6
  45 :           CONT(I)=CONT(I)+DDHOOK(I,J)*EPSI(J)
  46 :         END DO
  47 :       END DO
  48 : C
  49 :       RETURN
  50 :       END
  51 :