vol.for

Go to the documentation of this file.

      SUBROUTINE vol(SSEV,PV,PPV,K,DET)
C>     VOLUMETRIC CONTRIBUTION :STRAIN ENERGY FUNCTION AND DERIVATIVES
      IMPLICIT NONE
      include 'PARAM_UMAT.INC'
C
      DOUBLE PRECISION SSEV,PV,PPV
      DOUBLE PRECISION K,G,DET,AUX
C
      g=(one/four)*(det*det-one-two*log(det))
C
      ssev=k*g
C
      pv=k*(one/two)*(det-one/det)
      aux=k*(one/two)*(one+one/(det*det))
      ppv=pv+det*aux
C
      RETURN
      END SUBROUTINE vol