Esope source files

   1 : C INCRE8    SOURCE    CHAT      05/01/13    00:35:34     5004           
   2 :         
   3 :        SUBROUTINE INCRE8(SIG,VAR,T0,TF,EPSVPT,VARPT,XMAT,NSTRS0,                      
   4 :      &     MFR,NVARI,NCOMAT)  
   5 :       IMPLICIT INTEGER(I-N)
   6 :        IMPLICIT REAL*8(A-H,O-Z)                                                 
   7 :        DIMENSION SIG(*),VAR(*),EPSVPT(*)
   8 :        DIMENSION VARPT(*),XMAT(*)                       
   9 :        DIMENSION XX(6)              
  10 :        DIMENSION SIG0(6),EPS0(6)  
  11 :        PARAMETER (AMAX = 1.0D20 , AMIN = 1.D-10 ,ZER=0.D0 )                                
  12 :        DETIER = 2.0D0/3.0D0                                                     
  13 :        ROOT = SQRT(3.0D0/2.0D0)                                                      
  13 :        CALL ZDANUL(VARPT,NVARI)
  14 :        IF (MFR.EQ.5) THEN
  15 :          NSTRS=6
  16 :          SIG0(1)=SIG(1)
  17 :          SIG0(2)=SIG(2)
  18 :          SIG0(3)=0.D0
  19 :          SIG0(4)=SIG(3)
  20 :          SIG0(5)=SIG(4)
  21 :          SIG0(6)=SIG(5)
  22 :        ELSE
  23 :          NSTRS=NSTRS0
  24 :          DO 10 I=1,NSTRS
  25 :           SIG0(I)=SIG(I)
  26 :   10     CONTINUE
  27 :        ENDIF 
  28 :           DO 20 I=1,NSTRS
  29 :           EPS0(I) = ZER  
  30 :   20     CONTINUE                                                                      
  31 : C-------------------------------------------------------------------|
  32 : C                      COEFFICIENTS MATERIAU                        |
  33 : C-------------------------------------------------------------------|
  34 :        YOU =XMAT( 1)                                                            
  35 :        XNU =XMAT( 2)                                                            
  36 :        AA=XMAT( 5)
  37 :        B  =XMAT( 6)
  38 :        RM  =XMAT( 7)
  39 :        Q =XMAT( 8)
  40 :        R =XMAT( 9)
  41 :        H0 =XMAT(10)
  42 :        AP =XMAT(11)
  43 :        SB =XMAT(12)
  44 :        RN =XMAT(13)
  45 :        S0 =XMAT(14) 
  46 : C
  47 : C ------------------------------------------------------------------|           
  48 : C                        PARTIE VISCOPLASTIQUE                      |           
  49 : C-------------------------------------------------------------------I
  50 : C
  51 : C
  52 :        TRACE =(SIG0(1)+SIG0(2)+SIG0(3))/3.0D0  
  53 :        DO 71 I=1,NSTRS                                                          
  54 :        A = 0.0D0                                                                
  55 :        IF (I.LE.3) A=1.0D0     
  56 :        XX(I) = SIG0(I)-A*TRACE                   
  57 :   71   CONTINUE                                                                 
  58 : C
  59 : C *** CRITERE VISCOPLASTIQUE
  60 : C
  61 :        AJ2 = PROCON(XX,XX,NSTRS)                               
  62 :        AJ2 = SQRT(1.5D0*AJ2)  
  63 :        IF (VAR(2).LT.AMIN) VAR(2)= S0        
  64 :        RR =  SINH (B * AJ2/VAR(2))
  65 : C       
  66 : C   temperature au debut du pas 
  67 : C       
  68 :        QSRT = Q/(R*T0)
  69 :        Z =  AA*XPUISS(RR,1.D0/RM,AMAX) 
  70 :        VARPT(1) = Z*EXP(-QSRT)
  71 :        SSAT = SB * XPUISS(Z/AA,RN,AMAX) 
  72 :        DSAT =  SSAT - S0      
  73 :        SSIG =  ABS(DSAT )/DSAT
  74 :        RK =   ABS((SSAT - VAR(2))/ DSAT)
  75 :        VARPT(2)= H0*XPUISS(RK,AP,AMAX)*SSIG*VARPT(1)       
  76 :        IF (AJ2.LT.AMIN)  GOTO 100  
  77 :        DO  70 I=1,NSTRS
  78 :        EPS0(I) =XX(I)/AJ2*VARPT(1)/DETIER
  79 : 70     CONTINUE 
  80 : C
  81 : 100    CONTINUE      
  82 : C    
  83 :         IF (MFR.EQ.5) THEN  
  84 :          EPSVPT(1)=EPS0(1)  
  85 :          EPSVPT(2)=EPS0(2) 
  86 :          EPSVPT(3)=EPS0(4) 
  87 :          EPSVPT(4)=EPS0(5) 
  88 :          EPSVPT(5)=EPS0(6) 
  89 :        ELSE 
  90 :          DO 11 I=1,NSTRS 
  91 :          EPSVPT(I)=EPS0(I) 
  92 :  11      CONTINUE 
  93 :        ENDIF
  94 :        RETURN                                                                   
  95 :        END                                                                      
  96 :  
  97 :  
  98 :  
  99 :  
 100 :  
 101 :