帮忙看看写的如何 关于温度场和应力场的

/prep7
ETCHG,TTS ！将热单元转化为结构单元
ex2=
alpx2=
nuxy2=
et,1,solid45 ！选择单元模型
mp,ex,2,ex2 ！定义杨氏模量
mp,alpx,2,alpx2 ！定义热膨胀系数
mp,nuxy,2,nuxy2 ！定义泊松比
allsel,all
dof,ux,uy,uz,rotx,roty,rotz
cmsel,s,Aduichena (垂直于Z的平面）
nsla,s
d,all,uz,0
d,all,rotx,0
d,all,roty,0
allsel,all
dof,ux,uy,uz,rotx,roty,rotz
cmsel,s,Aduichenb （垂直于x的平面） //选名字为Aduichenb 的组
nsla,s
d,all,ux,0
d,all,roty,0
d,all,rotz,0
allsel,all
nsel,s,loc,x，0.888 （边界）
d,all,uy,0
d,all,uz,0
d,all,rotx,0
d,all,roty,0
d,all,rotz,0
allsel,all
finish
/solu
tref,70 //参考温度
OUTRES,,ALL
allsel,all
antype,static 定义为稳态分析
ACEL, acelx, acely, acelz //对物体设立加速度
ACEL, ,-9.8
SFGRAD, PRES,0,y,0.222,76440（
CMSEL,s,down
SF,all,PRES,0 //对上述节点群施加面荷载
allsel,all
TM_START=1e-6
TM_END=2000
TM_INCR=200
*DO,TM,TM_START,TM_END,TM_INCR
time,TM
ldread,temp,,,TM,,bbbb,rth
SOLVE
*ENDDO