关于abaqus文档给的UMIXMODEFATIGUE子程序为什么报错？

littlehuster

2018年后abaqus提供了自定义裂纹扩展速率的接口UMIXMODEFATIGUE

abaqus文档给了一个简单的例子：

           SUBROUTINE UMIXMODEFATIGUE (DADN, GI_MAX, GII_MAX, GIII_MAX,     1 GI_MIN, GII_MIN, GIII_MIN, TEMP, DTEMP, PREDEF, DPRED,     2 NFIELD, NPROPS, PROPS(NPROPS), NSTATV, STATEV,     3 NIARRAY, I_ARRAY, NRARRAY, R_ARRAY, NCARRAY, C_ARRAY, DA)      INCLUDE 'ABA_PARAM.INC'C      DIMENSION PROPS(NPROPS),STATEV(NSTATV),PREDEF(NFIELD),DPRED(NFIELD),     1 I_ARRAY(NIARRAY), R_ARRAY(NRARRAY), C_ARRAY(NCARRAY)      A=PROPS(1)      B=PROPS(2)      C=PROPS(3)      D=PROPS(4)      G_SHEAR=GII_MAX+GIII_MAX      G_TOTAL=GI_MAX+GII_MAX+GIII_MAX      RMODEMIX=G_SHEAR/G_TOTAL      IF(RMODEMIX.LE.0.5)THEN                DADN=A*G_TOTAL**B      ELSE         DADN=C*G_TOTAl**D      ENDIFC      RETURN      END

但是，这个例子却报了语法错误：PROPS(NPROPS)不可以为形参。

       能否提供解决这个问题的思路；或者给个能在这个子程序下，运行成功的inp文件也行。