flowvp.for

本程序是D.R.J.OWEN主编的FINITE ELEMENTS IN PLASTICITYG一书中第八章例题的有限元程序

C$DEBUG
C$LARGE
	SUBROUTINE FLOWVP(AVECT,PROPS,LPROP,STEFF,NSTR1,MTOTG,VIVEL,
     .                  YIELD,KGAUS,MMATS,NCRIT,FNORM,ALLOW)  !(Page 294)
C******************************************************************
C
C****THIS SUBROUTINE EVALUATES THE VISCOPLASTIC STRAIN RATE
C
C******************************************************************
	DIMENSION AVECT(4),PROPS(MMATS,10),VIVEL(5,MTOTG)
	ALLOW=0.01
	IF(STEFF.EQ.0.0) GO TO 90
	YOUNG=PROPS(LPROP,1)
	POISS=PROPS(LPROP,2)
	HARDS=PROPS(LPROP,6)
	FRICT=PROPS(LPROP,7)
	GAMMA=PROPS(LPROP,8)
	DELTA=PROPS(LPROP,9)
	NFLOW=PROPS(LPROP,10)
      ROOT3=1.73205080757
	FDATM=PROPS(LPROP,5)
	FRICT=FRICT*0.017453292
	IF(NCRIT.EQ.3) FDATM=FDATM*COS(FRICT)
	IF(NCRIT.EQ.4) FDATM=6.0*FDATM*COS(FRICT)/(ROOT3*(3.0-SIN(FRICT)))
	IF(HARDS.GT.0.0) FDATM=FDATM+VIVEL(5,KGAUS)*HARDS
	IF(FDATM.LT.0.001) FDATM=1.0
	FCURR=YIELD-FDATM
	FNORM=FCURR/FDATM
	IF(FNORM.LT.ALLOW) GO TO 90
	IF(NFLOW.EQ.1) GO TO 50
	CMULT=GAMMA*(EXP(DELTA*FNORM)-1.0)
	GO TO 60
   50 CMULT=GAMMA*(FNORM**DELTA)
   60 DO 70 ISTR1=1,NSTR1
   70 AVECT(ISTR1)=CMULT*AVECT(ISTR1)
      DO 80 ISTR1=1,NSTR1
   80 VIVEL(ISTR1,KGAUS)=AVECT(ISTR1)
      RETURN
   90 DO 100 ISTR1=1,NSTR1
  100 VIVEL(ISTR1,KGAUS)=0.0

	RETURN
	END