Esope source files

   1 : C INCREE    SOURCE    CHAT      05/01/13    00:35:37     5004           
   2 :       SUBROUTINE INCREE (TAU,SIG,P,QTLD,DECRA,DESWA,EPSVPT,EC0PT,ESW0PT,
   3 :      &                  NSTRS1,IVIEX,MFR1,IFOURB)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z)
   6 : C-----------------------------------------------------------------------
   7 : C
   8 : C     DESCRIPTION FONCTIONNELLE :
   9 : C     -------------------------
  10 : C
  11 : C     Lois 'VISCO_EXTERNE' : exploitation des sorties de CREEP pour le
  12 : C     calcul du tenseur des vitesses de deformation inelastique
  13 : C
  14 : C     Entrees : TAU,SIG,P,QTLD,DECRA,DESWA
  15 : C     Sorties : EPSVPT,EC0PT,ESW0PT
  16 : C
  17 : C     Parametres : NSTRS1,IVIEX,MFR1,IFOURB
  18 : C                  Les indicateurs IVIEX,MFR1,IFOURB ne servent pas
  19 : C                  pour l'instant mais seront utiles si extensions :
  20 : C                  autres modeles que 'GENERAL', autres formulations EF,
  21 : C                  autres options de calcul
  22 : C
  23 : C-----------------------------------------------------------------------
  24 : C     Arguments de l'interface
  25 : C
  26 :       REAL*8     SIG(*), DECRA(*), DESWA(*), EPSVPT(*)
  27 : C
  28 : C     Variables locales
  29 : C
  30 : C     Tableaux de travail
  31 : C
  32 :       REAL*8     SS(8)
  33 : C
  34 : C------------------- Debut du code executable --------------------------
  35 : C
  36 : C 1 - Initialisation du tenseur des vitesses de deformation inelastique
  37 : C
  38 :       DO 10 I=1,NSTRS1
  39 :          EPSVPT(I)=0.0D0
  40 :   10  CONTINUE
  41 : C
  42 : C 2 - Calcul des vitesses de deformation inelastique
  43 : C     (grandeurs scalaires)
  44 : C
  45 :       EC0PT = DECRA(1)/TAU
  46 :       ESW0PT = DESWA(1)/TAU
  47 : C
  48 : C 3 - Premiere contribution au tenseur des vitesses de deformation
  49 : C     inelastique : 1/3 * ESW0PT * [Id]
  50 : C
  51 :       CONS=0.333333333333333D0*ESW0PT
  52 :       DO 30 I=1,3
  53 :          EPSVPT(I)=EPSVPT(I)+CONS
  54 :   30  CONTINUE
  55 : C
  56 : C 4 - Seconde contribution au tenseur des vitesses de deformation
  57 : C     inelastique : (3/2)/QTLD * [SS]
  58 : C
  59 :       IF ( QTLD.LE.1.D-10 ) GOTO 1000
  60 : C
  61 : C     Calcul du deviateur des contraintes
  62 : C
  63 :       DO 40 I=1,NSTRS1
  64 :          A=0.0D0
  65 :          IF (I.LE.3) A=1.0D0
  66 :          SS(I)=SIG(I)-A*P
  67 :   40  CONTINUE
  68 : C
  69 :       CONS=(1.5D0*EC0PT)/QTLD
  70 :       DO 41 I=1,NSTRS1
  71 :          EPSVPT(I)=EPSVPT(I)+CONS*SS(I)
  72 :   41  CONTINUE
  73 : C
  74 : 1000  CONTINUE
  75 :       RETURN
  76 :       END
  77 :  
  78 :