Esope source files

   1 : C TUFADV    SOURCE    CHAT      05/01/13    03:54:58     5004           
   2 : C TUFADV    SOUCE    BROC      94/04/09    21:19:18
   3 : C----------------------------------------------------------------
   4 :       SUBROUTINE TUFADV(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,
   5 :      &     SIGPT,EPSVPT,VARPT,XMAT,XCAR,
   6 :      &     NSTRS,NVARI,IFOUR,INPLAS,NCOMAT     )
   7 : C -----------------------------------------------------------------
   8 :       IMPLICIT INTEGER(I-N)
   9 :       IMPLICIT REAL*8(A-H,O-Z)
  10 : C    Include contenant quelques constantes dont XPI :
  11 : -INC CCREEL
  12 :       DIMENSION SIG(*),EPSV(*),VAR(*),SIG1(*),EPSV1(*),VAR1(*)
  13 :       DIMENSION SIGPT(*),EPSVPT(*),VARPT(*),XMAT(*),XCAR(*)
  14 :       DIMENSION DTUFI(2,2),CTUFI(2,2)
  15 : C----------------------------------------------------------------
  16 : C  CETTE ROUTINE CALCULE LES VALEURS FINALES DE SIG EPSV VAR A
  17 : C  PARTIR DES VALEURS INITIALES DE L'INCREMENT DE TEMPS ET DES
  18 : C  DERIVEES SIGPT,EPSVPT,VARPT
  19 : C----------------------------------------------------------------
  20 : 
  21 :       THETA = XCAR(9)/2
  22 :       EPAIS = XCAR(2)
  23 :       RAYOM = XCAR(1)-(EPAIS/2)
  24 :       YOUN = XMAT(1)
  25 :       ANG=XPI*THETA/180.0D0
  26 : 
  27 : C  EVALUATION DE L INCREMENT D ELONGATION DE FLUAGE
  28 : 
  28 :       CALL TUFIHO(ANG,DTUFI,CTUFI,YOUN,RAYOM,EPAIS)
  29 :       DF = DTUFI(1,1)
  30 :       DEF= DTUFI(1,2)
  31 :       DE = DTUFI(2,2)
  32 :       CM = CTUFI(1,1)
  33 :       CMP= CTUFI(1,2)
  34 :       CP = CTUFI(2,2)
  35 :       PHIF = CM*SIG(6) + CMP*SIG(1)
  36 :       ELF = CMP*SIG(6) + CP*SIG(1)
  37 :       RAPF=ELF/PHIF
  38 : 
  39 :       DO 10 I=1,6
  40 :          SIG1(I)=SIG(I)
  41 :  10   CONTINUE
  42 : 
  43 : C________________________________________________________________
  44 : C   CALCUL DES DEFORMATIONS VISQUEUSES ELONGATION & ROTATION
  45 : C                DUES AU FLUAGE
  46 : C________________________________________________________________
  47 :       DO 20 I =2,6
  48 :          EPSV1(I) =EPSV(I)+TAU*EPSVPT(I)
  49 :  20   CONTINUE
  50 :       EPSV1(1)=RAPF*EPSV1(6)
  51 :       EPSVPT(1)=RAPF*EPSVPT(6)
  52 : C________________________________________________________________
  53 : C       MISE A JOUR DES VARIABLES INTERNES
  54 : C________________________________________________________________
  55 :       DO 30 I=1,5
  56 :          VAR1(I)=VARPT(I)
  57 :  30   CONTINUE
  58 :       VAR1(6)=VAR(6)+TAU*VARPT(6)
  59 :       VAR1(7)=TAU*VARPT(7)
  60 :       SIG1(1) = SIG(1) + TAU*(SIGPT(1) - DEF*EPSVPT(6)-DE*EPSVPT(1))
  61 :       SIG1(6) = SIG(6) + TAU*(SIGPT(6) - DF*EPSVPT(6)-DEF*EPSVPT(1))
  62 :  99   CONTINUE
  63 :       RETURN
  64 :       END
  65 : 
  66 : 
  67 : 
  68 : 
  69 :  
  70 :  
  71 :