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