Esope source files

   1 : C INCREM    SOURCE    STRU      06/12/14    21:15:18     5611           
   2 : c     Sous-Programme du module MISTRAL1
   3 : C     --------------------------------------------------------------------------
   4 :       SUBROUTINE INCREM (TT,FI,FIT,
   5 :      &                   SIG,EPS,EPSTH,EPSEL,EPSP,EPSCR,
   6 :      &                   EPSE,EPSEQ,XX,PSI,
   7 :      &                   TTP,FIP,
   8 :      &                   SIGP,VEPS,VEPSTH,VEPSEL,VEPSP,VEPSCR,
   9 :      &                   VEPSE,EPSEQP,XXP,PSIP,
  10 :      &                   NDPI,NDVP,NXX,
  11 :      &                   DT, KINCR,
  12 :      &                   TTF,FIF,FITF,
  13 :      &                   SIGF,EPSF,EPSTHF,EPSELF,EPSPF,EPSCRF,
  14 :      &                   EPSEF,EPSEQF,XXF,PSIF,
  15 :      &                   DELTT,DELFI,DELFIT, DELSIG,DELEPS,DELETH,
  16 :      &                   DELEEL,DELEP,DELECR,DELEE,DELEEQ,DELXX,DELPSI)
  17 : C     --------------------------------------------------------------------------
  18 : C     Initialisation a 0 des variations DELY des variables Y depuis le
  19 : C     debut de la periode d'integration, si KINCR = 0.
  20 : C     Calcul des valeurs YF a l'instant t+DT en fonction des valeurs Y
  21 : C     a l'instant t et de leurs derivees YP a t ou t+DT/2, si KINCR = 1 ou 2.
  22 : C     Incrementation des variations DELY des variables Y depuis le debut de la
  23 : C     periode d'integration en fonction des YP, si KINCR = 2.
  24 : C     --------------------------------------------------------------------------
  25 :       IMPLICIT INTEGER (I-N)
  26 :       IMPLICIT REAL*8 (A-H, O-Z)
  27 : c     std :
  28 :       PARAMETER ( NDVPMA = 4 , NXXMAX = 3 )
  29 : c     std.
  30 :       DIMENSION SIG(1:*),EPS(1:*),EPSTH(1:*),EPSEL(1:*),EPSCR(1:*),
  31 :      &          EPSP(1:6,0:*),EPSE(0:*),EPSEQ(0:*),XX(1:6,1:*)
  32 :       DIMENSION SIGP(1:*),VEPS(1:*),VEPSTH(1:*),VEPSEL(1:*),VEPSCR(1:*),
  33 :      &          VEPSP(1:6,0:*),VEPSE(0:*),EPSEQP(0:*),
  34 :      &          XXP(1:6,1:*)
  35 : c     std :
  36 :       DIMENSION DSIG(1:6),DEPS(1:6),DEPSTH(1:6),DEPSEL(1:6),DEPSCR(1:6),
  37 :      &          DEPSP(1:6,0:NDVPMA),DEPSE(0:NDVPMA),DEPSEQ(0:NDVPMA),
  38 :      &          DX(1:6,1:NXXMAX)
  39 : c     std.
  40 :       DIMENSION SIGF(1:*),EPSF(1:*),EPSTHF(1:*),EPSELF(1:*),EPSCRF(1:*),
  41 :      &          EPSPF(1:6,0:*),EPSEF(0:*),EPSEQF(0:*),
  42 :      &          XXF(1:6,1:*)
  43 :       DIMENSION DELSIG(1:*),
  44 :      &          DELEPS(1:*),DELETH(1:*),DELEEL(1:*),DELECR(1:*),
  45 :      &          DELEP(1:6,0:*),DELEE(0:*),DELEEQ(0:*),
  46 :      &          DELXX(1:6,1:*)
  47 :         IF (KINCR.EQ.0) THEN
  48 :          DELTT = 0.
  49 :          DELFI = 0.
  50 :          DELFIT = 0.
  51 :          DO IJ = 1,6
  52 :           DELSIG(IJ) = 0.
  53 :           DELEPS(IJ) = 0.
  54 :           DELETH(IJ) = 0.
  55 :           DELEEL(IJ) = 0.
  56 :           DELECR(IJ) = 0.
  57 :          END DO
  58 :          DO IDP = 1-NDPI,NDVP
  59 :           DO IJ = 1,6
  60 :            DELEP(IJ,IDP) = 0.
  61 :           END DO
  62 :           DELEE(IDP) = 0.
  63 :           DELEEQ(IDP) = 0.
  64 :          END DO
  65 :        DO IXX = 1,NXX
  66 :           DO IJ = 1,6
  67 :            DELXX(IJ,IXX) = 0.
  68 :           END DO
  69 :        END DO
  70 :          DELPSI = 0.
  71 :          RETURN
  72 :         END IF
  73 :         DTT = TTP*DT
  74 :         TTF = TT+DTT
  75 :         FID = FI
  76 :         DFI = FIP*DT
  77 :         FIF = FI+DFI
  78 :         DFIT = (FID+FIF)/2.*DT
  79 :         FITF = FIT+DFIT
  80 :         DO IJ = 1,6
  81 :          DSIG(IJ) = SIGP(IJ)*DT
  82 :          SIGF(IJ) = SIG(IJ)+DSIG(IJ)
  83 :          DEPS(IJ) = VEPS(IJ)*DT
  84 :          EPSF(IJ) = EPS(IJ)+DEPS(IJ)
  85 :          DEPSTH(IJ) = VEPSTH(IJ)*DT
  86 :          EPSTHF(IJ) = EPSTH(IJ)+DEPSTH(IJ)
  87 :          DEPSEL(IJ) = VEPSEL(IJ)*DT
  88 :          EPSELF(IJ) = EPSEL(IJ)+DEPSEL(IJ)
  89 :          DEPSCR(IJ) = VEPSCR(IJ)*DT
  90 :          EPSCRF(IJ) = EPSCR(IJ)+DEPSCR(IJ)
  91 :         END DO
  92 :         DO IDP = 1-NDPI,NDVP
  93 :          DO IJ = 1,6
  94 :           DEPSP(IJ,IDP) = VEPSP(IJ,IDP)*DT
  95 :           EPSPF(IJ,IDP) = EPSP(IJ,IDP)+DEPSP(IJ,IDP)
  96 :          END DO
  97 :          DEPSE(IDP) = VEPSE(IDP)*DT
  98 :          EPSEF(IDP) = EPSE(IDP)+DEPSE(IDP)
  99 :          DEPSEQ(IDP) = EPSEQP(IDP)*DT
 100 :          EPSEQF(IDP) = EPSEQ(IDP)+DEPSEQ(IDP)
 101 :         END DO
 102 :       DO IXX = 1,NXX
 103 :          DO IJ = 1,6
 104 :           DX(IJ,IXX) = XXP(IJ,IXX)*DT
 105 :           XXF(IJ,IXX) = XX(IJ,IXX)+DX(IJ,IXX)
 106 :          END DO
 107 :       END DO
 108 :         DPSI = PSIP*DT
 109 :         PSIF = PSI+DPSI
 110 :         IF (KINCR.EQ.2) THEN
 111 :          DELTT = DELTT+DTT
 112 :          DELFI = DELFI+DFI
 113 :          DELFIT = DELFIT+DFIT
 114 :          DO IJ = 1,6
 115 :           DELSIG(IJ) = DELSIG(IJ)+DSIG(IJ)
 116 :           DELEPS(IJ) = DELEPS(IJ)+DEPS(IJ)
 117 :           DELETH(IJ) = DELETH(IJ)+DEPSTH(IJ)
 118 :           DELEEL(IJ) = DELEEL(IJ)+DEPSEL(IJ)
 119 :           DELECR(IJ) = DELECR(IJ)+DEPSCR(IJ)
 120 :          END DO
 121 :          DO IDP = 1-NDPI,NDVP
 122 :           DO IJ = 1,6
 123 :            DELEP(IJ,IDP) = DELEP(IJ,IDP)+DEPSP(IJ,IDP)
 124 :           END DO
 125 :           DELEE(IDP) = DELEE(IDP)+DEPSE(IDP)
 126 :           DELEEQ(IDP) = DELEEQ(IDP)+DEPSEQ(IDP)
 127 :          END DO
 128 :        DO IXX = 1,NXX
 129 :           DO IJ = 1,6
 130 :            DELXX(IJ,IXX) = DELXX(IJ,IXX)+DX(IJ,IXX)
 131 :           END DO
 132 :        END DO
 133 :          DELPSI = DELPSI+DPSI
 134 :         END IF
 135 :         RETURN
 136 :         END
 137 :