Sources Esope | Cast3M

hujcri
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