Esope source files

   1 : C AVANX1    SOURCE    AM        08/12/19    21:15:03     6237           
   2 :       SUBROUTINE AVANX1(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
   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 :       DIMENSION D1(3,3),Di(3),ROT(3,3)
  37 : C     write(6,*) 'Entrée dans AVANX1'
  38 : C
  39 : C----------------------
  40 : C CALCULATE   VARIABLES   
  41 : C----------------------                          
  42 : 20    DO 50 I = 1,7                                                         
  43 : 50      VAR1(I) = VAR(I)+TAU*VARD(I)                                                  
  44 :       XD11=VAR1(2)
  45 :       XD22=VAR1(3)
  46 :       XD33=VAR1(4)
  47 :       XD12=VAR1(5)
  48 :       XD13=VAR1(6)
  49 :       XD23=VAR1(7)
  49 :       CALL ZERO(D1,3,3)                                                                                         
  50 :       D1(1,1) = XD11 
  51 :       D1(2,2) = XD22
  52 :       D1(3,3) = XD33                                                                                                                
  53 :       D1(1,2) = XD12
  54 :       D1(1,3) = XD13
  55 :       D1(2,3) = XD23
  56 :       D1(2,1) = D1(1,2)
  57 :       D1(3,1) = D1(1,3)
  58 :       D1(3,2) = D1(2,3)
  59 : C
  59 :       CALL JACOB3(D1,3,Di,ROT)
  60 : C
  61 :       ADMAX = DMAX1(Di(1),Di(2),Di(3))
  62 :       VAR1(8)=ADMAX
  63 : C                                                                                                                                                              
  64 : C-------------------------------            
  65 : C CALCULATE  PLASTIC  STRAINS                                     
  66 : C-------------------------------                                                                              
  67 :       DO  60 I = 1,NSTRS                                                        
  68 : 60    EPSV1(I) = EPSV(I)+TAU*EPSVD(I)                                          
  69 : C
  70 : C -------------------------------
  71 : C CALCULATE STRESS
  72 : C--------------------------------
  73 :       DO 70 I=1,NSTRS
  74 : 70      SIG1(I)=SIG(I)+TAU*SIGD(I)
  75 :       RETURN
  76 :       END
  77 :  
  78 :  
  79 :  
  80 :  
  81 :