Esope source files

   1 : C HUJCRI    SOURCE    FANDEUR   13/01/29    21:15:53     7683           
   2 :       SUBROUTINE HUJCRI(SIG,PI,XKSI,XMAT,XINV,Y,SELAS)
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 :       DIMENSION SIG(*),XMAT(*)
   6 :       DIMENSION XINV(2),COVNMS(1)
   7 : *
   8 :       ALFAH=0.D0
   9 :       COVNMS(1) = 0.D0
  10 : 
  11 :       XINV(1)=TRACE(SIG)/3.D0
  12 :       XINV(2)=VONMIS(SIG,1,ALFAH,COVNMS)
  13 :       COHE = XMAT(6)
  14 :       GM   = XMAT(7)
  15 :       HUB  = XMAT(10)
  16 :       DENOM =(PI+COHE) * EXP(-1.D0/HUB)
  17 :       RAP=(COHE-XINV(1))/DENOM
  18 :       XLRAP=LOG(RAP)
  19 :       Y=XINV(2)/(GM*(COHE-XINV(1))) + XKSI*(HUB*XLRAP - 1.D0)
  20 :       SELAS = XKSI
  21 :       RETURN                                                                    
  22 :       END                                                                       
  23 :  
  24 :