fini
/clear
/prep7
et,1,solid70
mptemp,1,0,227,727,927
mpdata,kxx,1,1,5
mpdata,kyy,1,1,5
mpdata,kzz,1,1,5
blc4,,,0.01,0.01,0.005
wpoffs,0.0046
wprota,,,90
vsbw,all
wpoffs,,,0.0008
vsbw,all
lsel,s,length,,0.0046
cm,l1,line
lesize,l1,,,20,0.25
lsel,s,length,,0.005
lsel,u,,,9,12,3
cm,l2,line
lesize,l2,,,20,0.25
lsel,s,,,9,12,3
cm,l3,line
lesize,l3,,,20,4
lsel,s,length,,0.01
cm,l4,line
lesize,l4,0.0002
lsel,s,length,,0.0008
cm,l4,line
lesize,l4,,,1
vmesh,all
physics,write,thermal
physics,clear
etchg,tts
mptemp,1,0,227,727,927
mpdata,ex,1,1,1e11,0.9815e11,0.6878e11,0.5711e11
mp,prxy,1,0.403
mp,dens,1,4520
mp,alpx,1,7.06e-6
tb,bkin,1,4
tbtemp,0,1
tbdata,1, 8e6,1e10
tbtemp,227,2
tbdata,1, 7.97e6,0.9815e10
tbtemp,727,3
tbdata,1, 5.80e6,0.6878e10
tbtemp,927,4
tbdata,1, 5.78e6,0.5711e10
physics,write,struct
physics,clear
fini
/solu
physics,read,thermal
tunif,20
nsel,s,loc,x,0.0046,0.0054
nsel,r,loc,z,0.0048,0.005
esln,s,1
a=elmiqr(0,13)
*get,b,elem,0,num,min
*dim,c,array,a
*do,i,1,a
c(i)=elnext(b+i-2)
*enddo
antype,trans
trnopt,full
autots,on
timint,on
allsel
*do,i,1,a
sfdel,all,all
time,10/6/a*i
deltim,0.01,0.005,0.01
kbc,1
esel,s,elem,,c(i)
nsle
nsel,r,loc,z,0.005
sf,all,hflux,3e7
sfdel,all,all
*enddo
nsel,s,loc,x,0
nsel,a,loc,x,0.01
sf,all,conv,20,20
nsel,s,loc,y,0
nsel,a,loc,y,0.01
sf,all,conv,20,20
nsel,s,loc,z,0
sf,all,conv,20,20
allsel
outres,all,all
solve
fini
/post26
nsol,2,26,temp
plvar,2
帮我看看,问题出在哪里?加热后,温度一直不变,为什么?问题在哪?? |