| 查看: 974 | 回复: 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年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有274人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有3人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
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
10