在降雨过程中同时实现边坡强度折减

小鸿

各位大神好，我已经获得土体弹塑性参数与含水率的关系，并应用usdfld子程序编了代码，那在边坡降雨过程中如何同时实现强度折减呢？是当做两个场变量处理还是重新编代码呢？非常感谢。
以下是我的子程序代码，请各位指导！
      SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,
     1 TIME,DTIME,CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,
     2 KSPT,KSTEP,KINC,NDI,NSHR,COORD,JMAC,JMATYP,MATLAYO,LACCFLA)
C
      INCLUDE 'ABA_PARAM.INC'
C
      CHARACTER*80 CMNAME,ORNAME
      CHARACTER*3  FLGRAY(15)
      DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),
     1 T(3,3),TIME(2)
      DIMENSION ARRAY(15),JARRAY(15),JMAC(*),JMATYP(*),COORD(*)
C     
      real W
      
C     取出饱和度数值
      CALL GETVRM('SAT',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,
     1 MATLAYO,LACCFLA)
      STATEV(1) = ARRAY(1)
C     取出孔隙比数值
      CALL GETVRM('VOIDR',ARRAY,JARRAY,FLGRAY,JRCD,JMAC,JMATYP,
     1 MATLAYO,LACCFLA)
      STATEV(2) = ARRAY(1)
C
      W = (STATEV(1) * STATEV(2))
      FIELD(1) = W
C
      STATEV(3) = （27907.97 - 122383.04 * W + 217548.04 * W * W）*100
      STATEV(4) = 0.19 + 0.3*W
      
      RETURN
      END

bizhonghui

老哥，可以加QQ互相讨论以下。1158877067