Esope source files

   1 : C TIDEFO    SOURCE    CHAT      05/01/13    03:38:17     5004           
   2 :       SUBROUTINE TIDEFO(V,DEFO,XL,IDEUX)
   3 : C-----------------------------------------------------------------------
   4 : C
   5 : C    ROUTINE DE CALCUL DES DEFORMATIONS DE LA POUTRE DROITE
   6 : C
   7 : C  ENTREE
   8 : C    V       VECTEUR DES DEPLACEMENTS LOCAUX AUX 2 NOEUDS DE LA POUTRE
   9 : C    XL      LONGUEUR
  10 : C    IDEUX   INDICE POUR DEFORMATIONS DU 2-EME ORDRE
  11 : C
  12 : C  SORTIE
  13 : C    DEFO    DEFORMATIONS AU MILIEU  6 PAR NOEUD
  14 : C
  15 : C    PP 4/9/92
  16 : C-----------------------------------------------------------------------
  17 :       IMPLICIT INTEGER(I-N)
  18 :       IMPLICIT REAL*8(A-H,O-Z)
  19 :       DIMENSION DEFO(*),V(*)
  20 :       DATA P5/0.5D0/
  21 : C
  22 :       RL=1.D00/XL
  23 :       RL2=RL*RL
  24 : C
  25 : C  DEFORMATIONS GENERALISEES
  26 : C
  27 :       EPS  = (V(7)-V(1))*RL
  28 :       EPTY = (V(8)-V(2))*RL - (V(6)+V(12))*P5
  29 :       EPTZ = (V(9)-V(3))*RL + (V(5)+V(11))*P5
  30 :       XSIX = (V(10)-V(4))*RL
  31 :       XSIY =-(V(5)-V(11))*RL
  32 :       XSIZ =-(V(6)-V(12))*RL
  33 : C
  34 : C  REMPLISSAGE DU TABLEAU DEFO
  35 : C
  36 :       DEFO(1)=EPS
  37 :       DEFO(2)=EPTY
  38 :       DEFO(3)=EPTZ
  39 :       DEFO(4)=XSIX
  40 :       DEFO(5)=XSIY
  41 :       DEFO(6)=XSIZ
  42 : C
  43 :       IF(IDEUX.EQ.0) RETURN
  44 : C
  45 : C  TERMES DU SECOND ORDRE pour la membrane
  46 : C
  47 :       AUX1=0.5D0*(EPS**2+((V(8)-V(2))*RL)**2 +((V(9)-V(3))*RL)**2)
  48 :       DEFO(1) = DEFO(1) + AUX1 
  49 :       RETURN
  50 :       END
  51 :  
  52 :