Esope source files

   1 : C AVANXX    SOURCE    CHAT      05/01/12    21:30:06     5004           
   2 :       SUBROUTINE AVANXX(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,
   3 :      &SIGD,EPSVD,VARD,NSTRS,NVARI)
   4 : C
   5 : C----------------------------------------------------------------
   6 : C Objet: Mise a jour des variables internes p,r,D
   7 : C                    des deformations plastiques 
   8 : C                    des contraintes
   9 : C----------------------------------------------------------------
  10 : C
  11 : C----------------------------------------------------------------
  12 : C Entree: TAU pas de temps
  13 : C         SIG(NSTRS,NSTRS) contraintes au debut du pas
  14 : C         EPSV(NSTRS,NSTRS) deformations plastiques
  15 : C           au debut du pas
  16 : C         VAR(NVARI) variables internes au debut du pas
  17 : C         EPSVD(NSTRS,NSTRS) vitesse de deformations
  18 : C           visco-plastiques calculees au debut du pas
  19 : C         VARD(NVARI) vitesse des variables internes au debut du pas
  20 : C         SIGD(NSTRS) vitesse des contraintes
  21 : C-------------------------------------------------------------
  22 : C
  23 : C-------------------------------------------------------------
  24 : C Sortie: SIG1(NSTRS,NSTRS) etat de contraintes a la fin du pas
  25 : C         EPSV1(NSTRS,NSTRS) deformations plastiques la fin
  26 : C           du pas
  27 : C         VAR1(NVARI) variables internes a la fin du pas        
  28 : C--------------------------------------------------------------                 
  29 : C
  30 :       IMPLICIT INTEGER(I-N)
  31 :       IMPLICIT REAL*8(A-H,O-Z)     
  32 :       logical dtlibr                                             
  33 :       DIMENSION SIG(*),EPSV(*),VAR(*)
  34 :       DIMENSION SIG1(*),EPSV1(*),VAR1(*)                  
  35 :       DIMENSION EPSVD(*),VARD(*),SIGD(*)                                                
  36 : C
  37 : C------------------------------------
  38 : C CALCULATE   VARIABIES :  P , R , D  
  39 : C------------------------------------                            
  40 : 20    DO 50 I = 1,NVARI                                                         
  41 :  50      VAR1(I) = VAR(I)+TAU*VARD(I)
  42 :         IF (VAR1(3).LE.0.99)  GOTO  60 
  43 : c modif antoine
  44 : c             TAU = (0.989D0-VAR(3))/VARD(3)
  45 : c             GOTO   20                                                              
  46 : C                                                                                                                                                              
  47 : C-------------------------------            
  48 : C CALCULATE  PLASTIC  STRAINS                                     
  49 : C-------------------------------                                                                              
  50 : 60      DO  70 I = 1,NSTRS                                                        
  51 : 70    EPSV1(I) = EPSV(I)+TAU*EPSVD(I)                                          
  52 : C
  53 : C -------------------------------
  54 : C CALCULATE STRESS
  55 : C--------------------------------
  56 :       DO 80 I=1,NSTRS
  57 : 80      SIG1(I)=SIG(I)+TAU*SIGD(I)
  58 :       RETURN
  59 :       END
  60 :  
  61 :  
  62 :  
  63 :