C SAVCR3    SOURCE    FD218221  26/02/13    21:15:30     12462          
      SUBROUTINE SAVCR3(VARF,NVARI,CMAT,NMAT,SIGF,
     $           SIGEF,NSTRS,IFOUR,EPSFLU,AN,BN,CN,SIGFI)
C
C_______________________________________________________________________
C
C  SOUS PROGRAMME APPELE PAR SAVCR QUI CALCUL LES MATRICES DE 
C  DEFORMATION DE FLUAGE DE DESSICCATION A LA FIN DU PAS DE TEMPS
C  D APRES LE MODELE DE F. BENBOUDJEMA (similaire au modele de Bazant) 
C_______________________________________________________________________
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION VARF(NVARI),CMAT(NMAT),SIGF(NSTRS),SIGEF(NSTRS)
      DIMENSION EPSFLU(6),SIGFI(NSTRS)
      DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
C
C RECUPERATION DES VARIABLES INTERNES INITIALES
C
C   FLUAGE PROPRE
C
      EPSFLU(1) = VARF(21)
      EPSFLU(2) = VARF(22)
      EPSFLU(3) = VARF(23)
      EPSFLU(4) = VARF(24)
      EPSFLU(5) = VARF(25)
      EPSFLU(6) = VARF(26)
C
C  TEST SI LE TEMPS EST NUL
C
      TDT = CMAT(28)
      IF (TDT.EQ.0.D0) THEN
        GOTO 30
      ENDIF
C
C INITIALISATION DES MATRICES
C
      CALL ZERO(AN,NSTRS,1)
      CALL ZERO(BN,NSTRS,NSTRS)
      CALL ZERO(CN,NSTRS,NSTRS)
      CALL ZERO(SIGFI,NSTRS,1)
C
      DO 10 I=1,NSTRS
          SIGFI(I) = SIGF(I)
  10  CONTINUE
C
C_______________________________________________________________________
C
C  FLUAGE PROPRE AMORTISSEUR
C
      NDI=NSTRS
      CALL DEFCR3(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS,
     $                 AN,BN,CN,NDI)
C
      DO 20 I=1,NSTRS
        EPSFLU(I) = EPSFLU(I) + AN(I)
        DO 21 J=1,NSTRS
          EPSFLU(I) = EPSFLU(I) + BN(I,J)*SIGEF(J) + CN(I,J)*SIGFI(J)
  21    CONTINUE
  20  CONTINUE
C
  30  CONTINUE
C
C  SAUVEGARDE DES DEFORMATIONS
C
C    FLUAGE DESSICCATION
C
      VARF(21) = EPSFLU(1)
      VARF(22) = EPSFLU(2)
      VARF(23) = EPSFLU(3)
      VARF(24) = EPSFLU(4)
      VARF(25) = EPSFLU(5)
      VARF(26) = EPSFLU(6)
C
      RETURN
      END