HARDSUB (SYIEL0,HARD,EQPLAS,PROP(3),NVALUE) 如何解释

liliu328

请问哪位用过这样的一个子程序，在help里找了半天也没有找到。只知道通过这个子程序能将等效塑性应变——转换成——屈服应力。哪位能具体解释一下？
SYIEL0,HARD 又是什么意思？
谢谢！

denghuakun

喔，我知道了，你肯定看的Abaqus材料用户子程序UMAT基础知识及手册例子完整解释的PDF吧？那个不是Abaqus自带的用户子程序，你把他那个程序往后翻，那其实是作者自己写的一个子程序，然后自己调用

liliu328

一个问题挂了五天了，求助呀！！！！！

liliu328

半个月过去了，哪位能提点一下呀？

denghuakun

我也想问，同求帮顶...

pl02707822

顶楼主  同样想知道HARDINI子程序的详细使用方法

liuyingtao

:victory::victory:

TerlisDL

那个是作者自己写的子程序，不是abaqus自带的。具体如下：
SUBROUTINE HARDSUB(SYIELD,HARD,EQPLAS,TABLE,NVALUE )
C
      INCLUDE 'ABA_PARAM.INC'
C
      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