Esope source files

   1 : C GRAIN     SOURCE    CHAT      05/01/13    00:20:57     5004           
   2 : C GRAIN     SOURCE    LLB        
   3 :       SUBROUTINE GRAIN(TAU,EVP1,EVP2,SIG1,SIG2,VAR1,VAR2,XMAT,
   4 :      &     NSTRS,NVARI,KERRE)
   5 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
   6 : C
   7 : C       CALCUL DE LA TAILLE DE GRAIN A L'AIDE D'UN SHEMA DE TYPE
   8 : C    RUNGE KUTTA 2.1. LE MODELE UTILISE EST CELUI DE BESSON.
   9 : C ________________________________________________________________________
  10 : C
  11 : C   EN ENTREE :
  12 : C
  13 : C    DE1     : LE TAUX DE DEFAUTS AU TEMPS i
  14 : C    GR1     : LA TAILLE DE GRAIN AU TEMPS i
  15 : C    TAU     : LE PAS DE TEMPS
  16 : C    EVP1    : LA VITESSE DE DEFORMATION VISCOPLASTIQUE AU TEMPS i
  17 : C    EVP2    : LA VITESSE DE DEFORAMTIOM VISCOPLATIQUE AU TEMPS i+TAU/2
  18 : C    SIG1    : LA CONTRAINTE AU TEMPS i
  19 : C    SIG2    : LA CONTRAINTE AU TEMPS i+TAU
  20 : C    XMAT    : LE TABLEAU DES CARACTERISTIQUES MATERIELLES AU TEMPS i+TAU
  21 : C    NSTRS   : LA DIMENSION DES VECTEURS EVP1...SIG1...
  22 : C
  23 : C   EN SORTIE :
  24 : C
  25 : C    DE2     : LE TAUX DE DEFAUT AU TEMPS i+TAU 
  26 : C    GR2     : LA TAILLE DE GRAIN AU TEMPS i+TAU
  27 : C
  28 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  29 : C
  30 :       IMPLICIT INTEGER(I-N)
  31 :       IMPLICIT REAL*8(A-H,O-Z)
  32 : -INC CCOPTIO 
  33 :       DIMENSION EVP1(*),EVP2(*),SIG1(*),VAR1(*),SIG2(*),VAR2(*) 
  34 :       DIMENSION XMAT(*),SIGS2(6)
  35 : C
  36 : C--------------RAPPEL DES DONNES MATERIAU CALCULEE EN DEBUT DE PAS DE TEMPS
  37 : C--------------ATTENTION, DANS L'ETAT ACTUEL, CE SONT LES SEULES DONNEES
  38 : C--------------DISPONIBLES. DANS L'AVENIR, IL FAUDRA REACTUALISER LES K0, K1
  39 : C--------------K2, K3 EN COURS DE PAS DE TEMPS.
  40 : C 
  41 :       GR=VAR1(NVARI-2)
  42 :       DE=VAR1(NVARI-1)
  43 :       G0 = XMAT(20)
  44 :       D0 = XMAT(21)
  45 :       RM = XMAT(22)
  46 :       RK0 = XMAT(23)
  47 :       RK1 = XMAT(24)
  48 :       RK2 = XMAT(25)
  49 :       RK3 = XMAT(26)      
  50 : C
  51 : C--------------ADIMENSIONNALISATION DES VARIABLES INTERNES
  52 : C
  53 :       DE1=DE/D0
  54 :       GR1=GR/G0
  55 : C
  56 : C--------------CALCUL DE L'ENERGIE DE DEFORMATION AU TEMPS i
  57 : C      
  58 :       ENERI=PROCON(EVP1,SIG1,NSTRS)
  59 : C
  60 : C--------------CALCUL DE L'ENERGIE DE DEFORMATION AU TEMPS i+TAU/2
  61 : C
  62 :       DO 10 I=1,NSTRS
  63 :        SIGS2(I)= 0.5*(SIG1(I)+SIG2(I))
  64 : 10    CONTINUE
  65 :       ENERIS2=PROCON(EVP2,SIGS2,NSTRS)
  66 : C
  67 : C--------------CALCUL DE D AU TEMPS i+TAU/2
  68 : C
  69 :       D1S2=DE1+((TAU/2.D0)*(((RK1/D0)*ENERI)-(RK2*DE1)))
  70 : C
  71 : C--------------CALCUL G AU TEMPS i+TAU/2
  72 : C 
  73 :       G1S2=GR1+(0.5D0*TAU*(((RK0/((G0*GR1)**RM))+(RK3*D0*DE1))/G0))
  74 : C
  75 : C--------------CALCUL DE G AU TEMPS i+1
  76 : C
  77 :       GR2=GR1+TAU*(((RK0/((G0*G1S2)**RM))+RK3*D0*D1S2)/G0)
  78 : C
  79 : C--------------CALCUL DE D AU TEMPS i+1
  80 : C
  81 :       DE2=DE1+(TAU*(((RK1/D0)*ENERIS2)-(RK2*D1S2)))
  82 : C
  83 : C   REMPLISSAGE DU TABLEAU DE VARIABLES INTERNES     
  84 : C
  85 :       VAR2(NVARI-2)=GR2*G0
  86 :       VAR2(NVARI-1)=DE2*D0
  87 : C     
  88 :       RETURN
  89 :       END
  90 : 
  91 : 
  92 : 
  93 :  
  94 :