关于creep子程序问题
本帖最后由 zhoulan 于 2019-11-12 17:33 编辑写了一个蠕变子程序,不能收敛,望大佬解惑。 SUBROUTINE CREEP(DECRA,DESWA,STATEV,SERD,EC,ESW,P,QTILD,
1TEMP,DTEMP,PREDEF,DPRED,TIME,DTIME,CMNAME,LEXIMP,LEND,
2COORDS,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(*),
1TIME(3),EC(2),ESW(2),COORDS(*)
C
A=6.2886E-6
AN=6.0629
C1=-13.68442+0.19957*TEMP-6.9875*0.0001*TEMP**2
C2=-13.35202+0.11178*TEMP
C3=37.23182-0.46776*TEMP+0.001488*TEMP**2
C4=-6.07565+0.04292*TEMP
C5=-49.85311+0.6351*TEMP-0.00198*TEMP**2
Q=17439.84337-191.83477*TEMP+0.70838*TEMP**2
C
DECRA(1)=A*QTILD**AN*(TIME(1)**C1+C2*TIME(1)**C3+C4*TIME(1)**C5-
1 (TIME(1)-DTIME)**C1-C2*(TIME(1)-DTIME)**C3-
2 C4*(TIME(1)-DTIME)**C5)*EXP(1.0)**(-Q/TEMP)
write(*,*)DECRA(1)
IF(LEXIMP.EQ.1) THEN
DECRA(5)=A*AN*QTILD**(AN-1)*(TIME(1)**C1+C2*TIME(1)**C3+C4*TIME(1)
1 **C5-(TIME(1)-DTIME)**C1-C2*(TIME(1)-DTIME)**C3-
2 C4*(TIME(1)-DTIME)**C5)*EXP(1.0)**(-Q/TEMP)
write(*,*)DECRA(5)
END IF
C
RETURN
END
你好,目前我也在做CREEP的工作,能请教几个问题吗
页:
[1]