Esope source files

   1 : C ADVAR     SOURCE    PV        07/11/23    21:15:16     5978           
   2 : C ADVAR   SOURCE    ROUX G-M      03/02/24    21:15:01     4587           
   3 :        SUBROUTINE ADVAR(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,DSPT,EPSVPT,                 
   4 :      & DEPST,VARPT,XMAT,NSTRS,NVARI,IFOUR,INPLAS,MFR)
   5 :             
   6 :        IMPLICIT REAL*8(A-H,O-Z)                                                 
   7 :        DIMENSION SIG(*),EPSV(*),VAR(*),SIG1(*),EPSV1(*),VAR1(*)                 
   8 :        DIMENSION DSPT(*),EPSVPT(*),VARPT(*),XMAT(*),DEPST(*)                     
   9 :        real*8 lamb,mu
  10 :        DATA UN/1.D0/
  11 :        DATA DEUX/2.D0/
  12 : 
  13 : C*******************************************************
  14 : C                               ENTREE DE ADVAC
  15 : C*******************************************************
  16 : 
  17 : 
  18 : C----------------------------------------                                                                             
  19 : C mise à jour des déformations plastiques                                             
  20 : C----------------------------------------                                            
  21 :        YOU=XMAT(1)                                                            
  22 :        XNU=XMAT(2)                                                                                                                  
  23 : 
  24 :        DO 60 I=1,NSTRS                                                        
  25 :          EPSV1(I)=EPSV(I)+TAU*EPSVPT(I)                         
  26 :   60   CONTINUE 
  27 : 
  28 : C----------------------------
  29 : C mise à jour des contraintes
  30 : C----------------------------                                                              
  31 :          LAMB = YOU*XNU/((UN+XNU)*(UN-DEUX*XNU))
  32 :          MU = YOU/(DEUX*(UN+XNU))
  33 :          TRACE=EPSVPT(1)+EPSVPT(2)+EPSVPT(3)                              
  34 :            DO 70 I=1,3                            
  35 :            SIG1(I) = SIG(I) + TAU*(DSPT(I)-LAMB*TRACE-DEUX*MU*EPSVPT(I))  
  36 :  70        CONTINUE
  37 :            DO 80 I=4,NSTRS
  38 :            SIG1(I) = SIG(I) + TAU*(DSPT(I)-DEUX*MU*EPSVPT(I))
  39 :  80        CONTINUE     
  40 :                                                            
  41 : C-----------------------------------
  42 : C mise à jour des variables internes                                    
  43 : C-----------------------------------
  44 :           DO 90 I=1,NVARI                                                            
  45 :           VAR1(I)=VAR(I)+TAU*VARPT(I)   
  46 :           VAR1(I)=VAR(I)+TAU*VARPT(I)   
  47 : 90        CONTINUE                     
  48 : 
  49 :        RETURN                                                                   
  50 :        END                                                                      
  51 :  
  52 :  
  53 :  
  54 : 
  55 :  
  56 :