C AVANXX    SOURCE    CHAT      05/01/12    21:30:06     5004           
      SUBROUTINE AVANXX(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,
     &SIGD,EPSVD,VARD,NSTRS,NVARI)
C
C----------------------------------------------------------------
C Objet: Mise a jour des variables internes p,r,D
C                    des deformations plastiques 
C                    des contraintes
C----------------------------------------------------------------
C
C----------------------------------------------------------------
C Entree: TAU pas de temps
C         SIG(NSTRS,NSTRS) contraintes au debut du pas
C         EPSV(NSTRS,NSTRS) deformations plastiques
C           au debut du pas
C         VAR(NVARI) variables internes au debut du pas
C         EPSVD(NSTRS,NSTRS) vitesse de deformations
C           visco-plastiques calculees au debut du pas
C         VARD(NVARI) vitesse des variables internes au debut du pas
C         SIGD(NSTRS) vitesse des contraintes
C-------------------------------------------------------------
C
C-------------------------------------------------------------
C Sortie: SIG1(NSTRS,NSTRS) etat de contraintes a la fin du pas
C         EPSV1(NSTRS,NSTRS) deformations plastiques la fin
C           du pas
C         VAR1(NVARI) variables internes a la fin du pas        
C--------------------------------------------------------------                 
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)     
      logical dtlibr                                             
      DIMENSION SIG(*),EPSV(*),VAR(*)
      DIMENSION SIG1(*),EPSV1(*),VAR1(*)                  
      DIMENSION EPSVD(*),VARD(*),SIGD(*)                                                
C
C------------------------------------
C CALCULATE   VARIABIES :  P , R , D  
C------------------------------------                            
20    DO 50 I = 1,NVARI                                                         
 50      VAR1(I) = VAR(I)+TAU*VARD(I)
        IF (VAR1(3).LE.0.99)  GOTO  60 
c modif antoine
c             TAU = (0.989D0-VAR(3))/VARD(3)
c             GOTO   20                                                              
C                                                                                                                                                              
C-------------------------------            
C CALCULATE  PLASTIC  STRAINS                                     
C-------------------------------                                                                              
60      DO  70 I = 1,NSTRS                                                        
70    EPSV1(I) = EPSV(I)+TAU*EPSVD(I)                                          
C
C -------------------------------
C CALCULATE STRESS
C--------------------------------
      DO 80 I=1,NSTRS
80      SIG1(I)=SIG(I)+TAU*SIGD(I)
      RETURN
      END