看看这个可以吗
/prep7
et,1,55
keyopt,1,3,1
mp,kxx,1,1.25
k,1,5
k,2,6
k,3,12
k,4,12,0.25
k,5,6,0.26
k,6,6,1
k,7,5,1
k,8,5,0.25
a,1,2,5,8
a,2,3,4,5
a,8,5,6,7
esize,0.125
amesh,all
eplot
fini
/solu
nsel,s,loc,x,5
sf,all,conv,1,450
nsel,s,loc,x,6,12
nsel,r,loc,y,0.25,1
nsel,a,loc,x,12
sf,all,conv,0.25,70
nsel,all
/psf,conv,hcoef,1 !显示面载荷标记
nplot
solve
fini
/post1
plnsol,temp !把结果显示为连续向量
fini
/prep7
etchg,tts !转换单元类型
keyopt,1,3,1
mp,ex,1,28e6
mp,nuxy,1,0.3
mp,alpx,1,0.9e-5
fini
/solu
nsel,s,loc,y,0
nsel,a,loc,y,1
dsym,symm,y !注明节点上的自由度限制为对称或非对称
nsel,s,loc,x,5
sf,all,pres,1000
nsel,all
/pbc,all,1
/psf,pres,,1
nplot
tref,70 !为热应变设置参考温度
ldread,temp,,,,,,rth
/pbc,all,,0
/pbf,temp,,1
eplot
solve
fini
/post1
plnsol,s,eqv
finish |