板材淬火，应力结果跟理论相差很大，求高手查看代码

干吃泡面

FINISH
/CLEAR
/PREP7
ET,1,70                        !选用单元70为研究单元
MP,KXX,1,113             !设置的热导率
MP,DENS,1,2820                  !设置密度
MP,C,1,888                      !设置比热容
mp,alpx,1,27.6e-6
MPTEMP,1,20,200
mpdata,ex,1,,68E9,56E9
MPDATA,prxy,1,,0.3,0.3
BLC4,0,0,0.5,0.25,0.1
ESIZE,0.05
MSHKEY,1
VMESH,ALL
ASEL,S,,,2,6,2
CM,AA,AREA
ALLSEL

/SOLU                        !进入求解器
TUNIF,540                !设置初始温度
SFA,AA,1,CONV,3000,20       !设置表面换热系数及末了温度
ANTYPE,TRANS        !定义瞬态分析
KBC,1             !设置载荷类型为阶跃载荷
OUTRES,,ALL                !定义结果输出，即选择的是写出每个子步的结果
TOFFST,273            !定义温度偏移量

t=0
j=10
*do,i,1,60
DELTIM,10,10,10!定义初始时间子步步长，最小及最大时间步长
t=t+j
TIME,t            !设置淬火时间
SOLVE
*enddo
save

!*************************温度场转换成应力场显示***************************
/prep7
/com,structural
alls
eplot
etchg,tts
lsclear,all
/solu
ASEL,S,,,1,5,2
CM,KD,AREA
ALLSEL
DA,KD,SYMM
antype,4
trnopt,full            !完全瞬态积分法
lumpm,0                !使用与单元相关的质量矩阵公式
t=0
j=10
*do,i,1,60
/solu
ldread,temp,i,,,,,rth!
/pbf,temp,,1  
t=t+j
time,t
deltim,10,10,10
solve
*enddo
求解出来的应力，心部应力应该是先压后拉，但时间历程显示出来的不对，全是拉应力或者压应力。求高手看一下，是不是代码错了，或者其他方面指点迷津。先谢谢了

雷雷Cruz

楼主，单位统一没？ m,k......

干吃泡面

2# 雷雷Cruz
单位都是统一的，就是没有设置那个命令

bxx

LZ，我也是做这方面的，方便的下加下QQ 1035626813