Esope source files

   1 : C PARAM1    SOURCE    CHAT      05/01/13    02:09:00     5004           
   2 :       SUBROUTINE PARAM1(CMAT,NMAT,AMBDA1,AMBDA2,X1,X2)
   3 : C
   4 : C ROUTINE APPELE DANS ESTIM
   5 : C_______________________________________________________________________
   6 : C
   7 : C CALCUL DE PARAMETRE DE LA LOI DE FLUAGE PROPRE SPHERIQUE
   8 : C_______________________________________________________________________
   9 : C
  10 :       IMPLICIT INTEGER(I-N)
  11 :       IMPLICIT REAL*8(A-H,O-Z)
  12 :       DIMENSION CMAT(NMAT)
  13 : C
  14 : C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
  15 : C
  16 :       RRSP = CMAT(24)
  17 :       VRSP = CMAT(25)
  18 :       RISP = CMAT(26)
  19 :       VISP = CMAT(27)
  20 : C
  21 : C PARAMETRES INTERMEDIAIRES => EQUATION (2.12-1)
  22 : C
  23 :       URR = RRSP/VRSP
  24 :       UII = RISP/VISP
  25 :       URI = RRSP/VISP
  26 : C
  27 : C PARAMETRE INTERMEDIAIRE => EQUATION (2.12-2)
  28 : C
  29 :       DELTA = (URR - UII)**2 + 8*URI*(URR + UII) + 16*URI**2
  30 : C
  31 : C INVERSE DES TEMPS CARACTERISTIQUES => EQUATION (2.12-4)
  32 : C
  33 :       AMBDA1 = - (URR + 4*URI + UII + DSQRT(DELTA))/2
  34 :       AMBDA2 = - (URR + 4*URI + UII - DSQRT(DELTA))/2
  35 : C
  36 : C PARAMETRES INTERMEDIAIRES => EQUATION (2.12-3)
  37 : C
  38 :       X1 = (AMBDA1 + UII)/(2*URI)
  39 :       X2 = 2*URI/(AMBDA2 + UII)
  40 : C
  41 :       END
  42 :  
  43 :