找回密码
 注册
Simdroid-非首页
查看: 762|回复: 1

[子程序] 关于creep子程序问题

[复制链接]
发表于 2019-11-11 17:07:41 | 显示全部楼层 |阅读模式 来自 中国
本帖最后由 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


发表于 2021-11-10 16:25:12 | 显示全部楼层 来自 陕西西安
Simdroid开发平台
你好,目前我也在做CREEP的工作,能请教几个问题吗
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-28 16:43 , Processed in 0.028944 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表