24小时热门版块排行榜

返回列表

wudongquan

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 87.7
帖子: 54
在线: 83.9小时
虫号: 2675168
注册: 2013-09-23
专业: 固体力学

[交流] abaqus蠕变子程序编译一直不通过，请教各位大神，子程序中有什么问题吗？已有2人参与

最近编了个简单的蠕变子程序，基于韧性耗散模型以及与应力有关的蠕变韧性模型。但是abaqus编译一直不通过，请教各位大神，子程序中有什么问题吗？
   SUBROUTINE CREEP(DECRA,DESWA,STATEV,SERD,EC,ESW,P,QTILD,
   1 TEMP,DTEMP,PREDEF,DPRED,TIME,DTIME,CMNAME,LEXIMP,LEND,
   2 COORDS,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)
C
   INCLUDE 'ABA_PARAM.INC'
C
   CHARACTER*80 CMNAME
C
   DIMENSION DECRA(5),DESWA(5),STATEV(*),PREDEF(*),DPRED(*),TIME(2),
   # COORDS(*),EC(2),ESW(2)
   DIMENSION  PROPS(4),PS(3)
   REAL ECRI
   DOUBLE PRECISION DCRD,DCR,DAMAGE,EVOL
   PARAMETER (ONE=1.0D0,TWO=2.0D0,THREE=3.0D0,SIX=6.0D0, HALF=0.5D0,
   # PAI=3.1415926D0, TOLER=1.0D-6,NEWTON=40D0,THETA=0.8D0,
   # LLMT=6.3D1,HLMT=1.05D2)
   COMMON /NEWTV/STORD(800)
SAVE /NEWTV/

   PROPS(3)=2.6353E-16
   PROPS(4)=5.11005

   DO K1=1,5
      DECRA(K1)=0.D0
      DESWA(K1)=0.D0
   END DO

   IF(QTILD.LE.LLMT) THEN
         ECRI=0.02
   ELSE IF(QTILD.GE.HLMT) THEN
         ECRI=0.85
   ELSE
         ECRI=1.14E-14*(QTILD**5.8)+0.00537
   END IF
      CRSTRANF1=SINH(2*(PROPS(4)-0.5)/((PROPS(4)+0.5)*3))

      CRSTRANF2=SINH(2*((PROPS(4)-0.5)/(PROPS(4)+0.5))*ABS(STATEV(2)
   # /STATEV(1)))

      CRSTRANF=CRSTRANF1/CRSTRANF2

      DECRA(1)=DTIME*PROPS(3)*(QTILD**PROPS(4))

      DCRD=DECRA(1)/(CRSTRANF*ECRI)
      DAMAGE=STATEV(8)+DCRD
      IF(DAMAGE.GE.ONE) DAMAGE=ONE

      IF((ABS(DAMAGE-0.999)).LT.0.0001) DAMAGE=0.99901

      STATEV(8)=DAMAGE
      STATEV(9)=DECRA(1)

         DECRA(5)=DTIME*PROPS(3)*(QTILD**(PROPS(4)-1))*PROPS(4)

      IF(NOEL.EQ.1610.AND.NPT.EQ.1)  STORD(1)=STATEV(8)
      IF(NOEL.EQ.1610.AND.NPT.EQ.3)  STORD(2)=STATEV(8)
      IF(NOEL.EQ.1201.AND.NPT.EQ.1)  STORD(3)=STATEV(8)
      IF(NOEL.EQ.1201.AND.NPT.EQ.3)  STORD(4)=STATEV(8)
      IF(NOEL.EQ.796.AND.NPT.EQ.1)  STORD(5)=STATEV(8)
      IF(NOEL.EQ.796.AND.NPT.EQ.3)  STORD(6)=STATEV(8)



      DO K2=1,396
            I2=397-K2
            IF(NOEL.EQ.I2.AND.NPT.EQ.1)  STORD(5+2*K2)=STATEV(8)
            IF(NOEL.EQ.I2.AND.NPT.EQ.3)  STORD(6+2*K2)=STATEV(8)
         END DO



      IF(NOEL.EQ.397.AND.NPT.EQ.1) STORD(799)=STATEV(8)
      IF(NOEL.EQ.397.AND.NPT.EQ.3) STORD(800)=STATEV(8)

C
   RETURN
   END

回复此楼