Sources Esope | Cast3M

dphook
C DPHOOK    SOURCE    CHAT      05/01/12    22:54:48     5004
      SUBROUTINE DPHOOK(EPSX,EPSY,EPSZ,EPSXY,SIGX,SIGY,SIGZ,SIGXY,
     1 YNU,YOUN)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C------------------------------------------------------------------
C       PERMET DE CALCULER LE TENSEUR SIGMA EN FONCTION DU TENSEUR
C       EPSILON EN DEFORMATION PLANE
C------------------------------------------------------------------
C
      HOOK=YOUN/(1.D0+YNU)/(1.D0-2.D0*YNU)
      SIGX=HOOK*((1.D0-YNU)*EPSX+YNU*EPSY)
      SIGY=HOOK*(YNU*EPSX+(1.D0-YNU)*EPSY)
      SIGZ=HOOK*YNU*(EPSX+EPSY)
      SIGXY=HOOK*(1.D0-2.D0*YNU)*EPSXY
      RETURN
      END