24小时热门版块排行榜

返回列表

121100220

捐助贵宾 (初入文坛)

应助: 0 (幼儿园)
金币: 83.2
帖子: 9
在线: 3.1小时
虫号: 15855106
注册: 2019-08-06
专业: 高分子材料的加工与成型

[求助] ANSYS 生死单元加载时间问题已有1人参与

求助论坛大神，小弟最近在做增材制造温度场的仿真，用到ANSYS生死单元功能，在运行过程中遇到一个关于加载时间的问题。主要问题就是在执行*DO循环时一直报错“present time is equal or less than previous time”。自己实在找不到问题在哪，恳请论坛大神出手相助，帮忙看一下问题出在哪，小弟不胜感激。APDL代码如下:

/CONFIG,NRES,40000
LTH_MD=0.06  ！模具长度
LTH_E=0.001 ！单元长度
NA=6             ！层数
THICC=0.000125 ！单层厚
THICK=0.01       ！模具厚度
THICA=THICC*NA  ！铺层厚度
T_I=0.5             ！时间步
YXRZ=16          ！有效热载长度
DERTAX=5       ！铺放步
hflux_cons=88541  ！热流密度

/PREP7
ET,1,13,4
ET,1,13,4
！-----------------------------材料热参数---------------------------------
MPTEMP,1,0,50,100,150,200,250$MPTEMP,7,300,350,400,600
MPDATA,DENS,1,1,1601,1598,1593,1586,1575,1563
MPDATA,DENS,1,7,1551,1537,1524,1472
MPDATA,C,1,1,800,930,1040,1260,1300,1400
MPDATA,C,1,7,1550,1650,1700,1900
MPDATA,KXX,1,1,3.5,4.6,5.1,5.9,5.9,6.1
MPDATA,KXX,1,7,6.7,6.8,7,7.8
MPDATA,KYY,1,1,0.42,0.52,0.6,0.7,0.7,0.7
MPDATA,KYY,1,7,0.75,0.68,0.65,0.53

MPTEMP,1,0,50,100,150,200,250,$MPTEMP,7,300,330,343,350,400,600
MPDATA,EX,1,1,128e9,128e9,128e9,128e9,128e9,128e9
MPDATA,EX,1,7,128e9,128e9,128e9,128e9,128e9,128e9
MPDATA,EY,1,1,11E9,11E9,11E9,11E9,11E9,11E9
MPDATA,EY,1,7,11E9,11E9,1E-6,1E-6,1E-6,1E-6
MPDATA,EZ,1,1,11E9,11E9,11E9,11E9,11E9,11E9
MPDATA,EZ,1,7,11E9,11E9,11E9,11E9,11E9,11E9
MPDATA,NUXY,1,1,0.28,0.28,0.28,0.28,0.28,0.28
MPDATA,NUXY,1,7,0.28,0.28,1E-9,1E-9,1E-9,1E-9,1E-9
MPDATA,NUXZ,1,1,1E-9,1E-9,1E-9,1E-9,1E-9,1E-9
MPDATA,NUXZ,1,7,1E-9,1E-9,1E-9,1E-9,1E-9,1E-9
MPDATA,NUYZ,1,1,1E-9,1E-9,1E-9,1E-9,1E-9,1E-9
MPDATA,NUYZ,1,7,1E-9,1E-9,1E-9,1E-9,1E-9,1E-9
MPDATA,GXY,1,1,5.1E9,5.1E9,5.1E9,5.1E9,5.1E9,5.1E9
MPDATA,GXY,1,7,5.1E9,5.1E9,5.1E9,5.1E9,5.1E9,5.1E9

MPTEMP,1,0,50,100,140,143,250
MPTEMP,7,300,330,343,350,400,600
MPDATA,ALPX,1,1,0.6E-6,0.6E-6,0.6E-6,0.6E-6,0.6E-6,0.6E-6
MPDATA,ALPX,1,7,0.6E-6,0.6E-6,0.6E-6,0.6E-6,0.6E-6,0.6E-6
MPDATA,ALPY,1,1,34.3E-6,34.3E-6,34.3E-6,34.3E-6,66.0E-6,66.0E-6
MPDATA,ALPY,1,7,66.0E-6,66.0E-6,1E-30,1E-30,1E-30,1E-30

MPTEMP,1,0,50,100,150,200,250
MPTEMP,7,300,330,343,350,400,600
MPDATA,ENTH,1,1,0.00E+00,7.43E+07,1.57E+08,2.57E+08,3.59E+08,4.69E+08
MPDATA,ENTH,1,7,5.89E+08,6.63E+08,7.31E+08,7.48E+08,8.78E+08,1.44E+09

MP,EX,2.73E9
MP,NUXY,2,0.3
MP,ALPX,2,23E-6
MP,KXX,2,237
MP,C,2,905
MP,DENS,2,2700
！---------------几何模型--------------------------
BLC4,0,THICK,LTH_MD,THICA
BLC4,0,0,LTH_MD,0.01
AGLUE,ALL
  ！---------------------------划分网格-------------------------
LSEL,S,LINE,,1,5,2
LESIZE,ALL,,,60,,,,,1
LSEL,S,LINE,,2,4,2
LESIZE,ALL,,,6,,,,,1
LSEL,S,LINE,,9,10,1
LESIZE,ALL,,,10,5,,,,1
ASEL,S,AREA,,1
AATT,1,,1,0
MAT,1
AMESH,ALL
ASEL,S,AREA,,3
AATT,2,,1,0
MAT,2
AMESH,ALL$FINISH

！----------------------加载-------------------------------------
/sol
ANTYPE,TRANS
AUTOTS,1
TINTP,0.005,,,1,,,
LNSRCH,1
TRNOPT,FULL
TIMINT,ON
TUNIF,20
NROPT,FULL
ESEL,S,ELEM,,1,(LTH_MD/LTH_E)*NA,1
EKILL,ALL
NSLE,S
D,ALL,,20,,,,TEMP
TIME,T_I
NSUBST,10
LSEL,S,LINE,,5
SFL,ALL,CONV,13,,25
LSEL,S,LINE,,9,10,1
SFL,ALL,CONV,13,,25
ESEL,S,ELEM,,1,DERTAX,1
EALIVE,ALL
NSLE,S
DDELE,ALL
SFE,ALL,,CONV,2,hflux_cons
ALLSEL
EQSLV,ICCG
SOLVE
FINISH

/sol
ANTYPE,,REST
T=T_I*2
TIME,T
NSUBST,10
AUTOTS,1
TINTP,0.005,,,1,,,
LNSRCH,1
ESEL,S,ELEM,,1,DERTAX*2,1
EALIVE,ALL
NSLE,S
DDELE,ALL
SFE,ALL,,CONV,2,hflux_cons
ALLSEL
OUTRES,ALL,ALL
EQSLV,ICCG
SOLVE
T=T+T_I
FINISH

/sol
ANTYPE,,REST
*DO,I,YXRZ+DERTAX,(LTH_MD/LTH_E),DERTAX
TIME,T
NSUBST,10
AUTOTS,1
TINTP,0.005,,,1,,,
LNSRCH,1
ESEL,S,ELEM,,I-DERTAX+1,I,1
EALIVE,ALL
NSLE,S
DDELE,ALL
SFE,ALL,,CONV,2,hflux_cons
ESEL,S,ELEM,,I-(YXRZ+DERTAX)+1,I-YXRZ,1
SFE,ALL,,CONV,,13
SFE,ALL,,CONV,2,25
ALLSEL
OUTRES,ALL,ALL
EQSLV,ICCG
SOLVE
T=T+T_I
*ENDDO
FINISH@wuming524

回复此楼