ansys电磁感应加热，有偿

hailinli

命令流神马的基本上完成了，电磁分析这块结果没有问题，主要是温度场分析结果不正确
先将命令流贴上来，帮我看看问题出在哪里了，谢谢

finish
/clear
freq=150000             ! frequency (Hz.)
pi=4*atan(1)            ! pi
cond=.392e7             ! maximum conductivity
muzero=4e-7*pi          ! free-space permeability
mur=200                 ! maximum relative permeability
skind=sqrt(1/(pi*freq*cond*muzero*mur)) ! skin depth
ftime=3            ! final time
tinc=.05                ! time increment for harmonic analysis
time=0                  ! initialize time
delt=.01
/prep7
emunit,mks
et,1,233,1,,1   
et,2,233,1,,1  
et,3,233,,,1  
et,4,151,,,1,,1        ! SURF151, thermal, radiation

mp,murx,1,1                ! air relative permeability
mp,murx,2,1   
mp,rsvx,2,1.2e-6
mp,murx,3,1   
mp,rsvx,3,1.2e-6

RECTNG,0,0.1,0,0.2
RECTNG,0.15,0.2,0,0.2
RECTNG,0,1.0,0,1.0
aovlap,ALL
NUMCMP,ALL
/pnum,area,1

r,4,0                   ! Real constant set for SURF151
ASEL,S,,,1
AATT,2,,1
ALLS
ASEL,S,,,2
AATT,3,,2
ALLS
ASEL,S,,,3
AATT,1,,3
ALLS


LSEL,S,LOC,X,0.05,0.08
LSEL,R,LOC,Y,0.2
LESIZE,ALL,,,20,5
ALLS
LSEL,S,LOC,X,0.05,0.08
LSEL,R,LOC,Y,0
LESIZE,ALL,,,20,0.2
ALLS
ASEL,S,,,1
LSLA,R
LSEL,R,LOC,X,0
LESIZE,ALL,,,20
ALLS
LSEL,S,LOC,X,0.1
LESIZE,ALL,,,20
ALLS
AMESH,1

LSEL,S,LOC,X,0.16,0.18
LSEL,R,LOC,Y,0.2
LESIZE,ALL,,,20,0.2

LSEL,S,LOC,X,0.16,0.18
LSEL,R,LOC,Y,0
LESIZE,ALL,,,20,5


LSEL,S,LOC,X,0.15
LESIZE,ALL,,,20
LSEL,S,LOC,X,0.2
LESIZE,ALL,,,20
AMESH,2

SMRTSIZE,2
AMESH,3
alls
*get,nmax,node,,num,max
lsel,s,loc,x,0.1
type,4

real,4
mat,2
lmesh,all
*get,emax,elem,,num,max
emodif,emax,4,nmax   ! modify element to add space node for radiation
et,4,0   
alls
ASEL,S,,,2
NSLA,R
CP,1,VOLT,ALL
*get,NODE_AMPS,NODE,,num,max
F,NODE_AMPS,AMPS,10000

ALLS
LSEL,S,LOC,X,0
LSEL,A,LOC,X,1.0
LSEL,A,LOC,Y,1.0
NSLL,R
D,ALL,AZ,0
ALLS

ASEL,S,,,1
NSLA,R
D,ALL,VOLT,0!短路导体涡流计算
ALLS
/SOLU
ANTYPE,HARMIC
HARFRQ,150000
KBC,1
SOLVE
alls
physics,write,emg
physics,clear
finish

/prep7
lsclear,all
et,1,55,,,1  
et,2,0
et,3,0
et,4,151,,,1,1,1        ! SURF151 element for radiation
keyopt,4,9,1
r,4,1,5.67e-8           ! form factor, Stefan-Boltzmann constant
mp,kxx,2,120
mp,c,2,1280
mp,dens,2,7850
mp,emis,2,.68                   ! emissivity

/solu
antype,trans
toffst,273
tunif,100
d,nmax,temp,25                  ! ambient temperature
cnvtol,heat,1                   ! convergence tolerance
kbc,1                           ! step loads
trnopt,full
autos,on                        ! auto time-stepping
deltim,1e-1,1e-1,delt,on        ! time step control
outres,basic,all   
physics,write,thermal
finish

*do,i,1,30   ! solution *do loop
time=time+tinc  
physics,read,emg
/solu

*if,i,eq,1,then
tunif,100   
*else
ldread,temp,last,,,,,rth
*endif
solve
finish
physics,read,thermal
/assign,esav,therm,esav
/assign,emat,therm,emat
/solu
parsav,scalar,parameter,sav
antype,trans
*if,i,gt,1,then
antype,trans,rest            ! thermal restart
*endif
parres,new,parameter,sav
time,time
esel,s,mat,,2
ldread,hgen,,,,2,,rst
esel,all
solve
finish
/assign,esav
/assign,emat
*enddo
finish
save