板材淬火,应力结果跟理论相差很大,求高手查看代码
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
求解出来的应力,心部应力应该是先压后拉,但时间历程显示出来的不对,全是拉应力或者压应力。求高手看一下,是不是代码错了,或者其他方面指点迷津。先谢谢了 楼主,单位统一没? m,k...... 2# 雷雷Cruz
单位都是统一的,就是没有设置那个命令 LZ,我也是做这方面的,方便的下加下QQ 1035626813
页:
[1]