发布时间 : 星期三 文章一起学习更新完毕开始阅读
DIMENSION TABLE(2,NVALUE) C
PARAMETER(ZERO=0.D0) C
C SET YIELD STRESS TO LAST VALUE OF TABLE, HARDENING TO ZERO C
SYIELD=TABLE(1,NVALUE) HARD=ZERO
C IF MORE THAN ONE ENTRY, SEARCH TABLE C
IF(NVALUE.GT.1) THEN DO K1=1,NVALUE-1
EQPL1=TABLE(2,K1+1)
IF(EQPLAS.LT.EQPL1) THEN EQPL0=TABLE(2,K1)
IF(EQPL1.LE.EQPL0) THEN WRITE(7,1)
1 FORMAT(//,30X,'***ERROR - PLASTIC STRAIN MUST BE `, 1 `ENTERED IN ASCENDING ORDER') CALL XIT ENDIF C
C CURRENT YIELD STRESS AND HARDENING C
DEQPL=EQPL1-EQPL0 SYIEL0=TABLE(1,K1) SYIEL1=TABLE(1,K1+1) DSYIEL=SYIEL1-SYIEL0 HARD=DSYIEL/DEQPL
SYIELD=SYIEL0+(EQPLAS-EQPL0)*HARD GOTO 10 ENDIF END DO 10 CONTINUE ENDIF RETURN END