| 查看: 947 | 回复: 2 | |||
[交流]
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 |
回复此楼» 猜你喜欢
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有6人回复
-
物理学I论文润色/翻译怎么收费?
已经有193人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有30人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有62人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
2楼2018-01-22 21:41:55
321caihong
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2783.5
- 帖子: 306
- 在线: 8.1小时
- 虫号: 2684487
- 注册: 2013-09-27
- 性别: MM
- 专业: 生物医用高分子材料
3楼2018-01-23 20:50:02
