| 查看: 948 | 回复: 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 |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有20人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有166人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
-
桂林理工大学物理学专业招收调剂
已经有7人回复
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
