感应加热仿真求助

zxc0228

想通过感应加热加热一块石墨板，程序如下，但是求解过程不收敛，请各位帮忙看看是什么问题，谢谢！
/batch,list
/title, Induction heating
/prep7
shpp,off
/com,
/com,
x1=0.01
x2=0.04
x3=0.05
x4=0.08
x5=0.09
x6=0.12
x7=0.13
x8=0.16
x9=0.17
x10=0.20
x11=0.20
y1=0.005
y2=0.055
y3=0.06
freq=30000                ! frequency (Hz.)
pi=4*atan(1)               ! pi
cond=.392e7                ! maximum conductivity
muzero=4e-7*pi             ! free-space permeability
mur=1                    ! maximum relative permeability
skind=sqrt(1/(pi*freq*cond*muzero*mur))  ! skin depth

ftime=1                    ! final time
tinc=.05                   ! time increment for harmonic analysis
time=0                     ! initialize time
delt=.01                   ! maximum delta time step
! Electromagnetic model
et,1,53,,,1                ! kongqi xianquan
et,2,53,,,1                ! shimo
emunit,mks
mp,murx,1,1
mp,murx,2,1
mp,murx,3,1
mp,rsvx,3,13e-6
mp,enth,3,0.548556e4      !enthalpy
mp,kxx,3,129               !conductivity
mp,emis,3,0.9              ! emissivity
rectng,x1,x2,0,y1
rectng,x3,x4,0,y1
rectng,x5,x6,0,y1
rectng,x7,x8,0,y1
rectng,x9,x10,0,y1
rectng,0,x11,y2,y3
rectng,0,0.2,0,0.19
aovlap,all
numcmp,area
asel,s,area,,7,8
aatt,1,1,1
asel,s,area,,1,5,1
aatt,2,2,1
asel,s,area,,6
aatt,3,3,2
asel,all

mshape,0,2d
asel,s,area,,1,5,1                !   划分网格
aesize,all,0.001
amesh,all
asel,s,area,,7
aesize,all,0.002
amesh,all
asel,s,area,,8
aesize,all,0.002
amesh,all                  ! mesh remaining areas
ksel,s,loc,y,y2         
kesize,all,skind/10         
ksel,s,loc,y,y3            
kesize,all,skind /2      
lsel,s,loc,x,x11/2         
lesize,all,,,1            
lsel,all
asel,all
mshape,0,2d
mshk,1
amesh,6
nsel,s,loc,x
d,all,az,0                 
nsel,all
esel,s,mat,,2
bfe,all,js,,,,2e7
alls
! Thermal model
et,4,55,,,1                ! PLANE55 thermal element, axisymmetric
agen,2,6,,,,,,,1
aatt,3,3,4
ksel,s,loc,y,y2           
kesize,all,skind/10      
ksel,s,loc,y,y3   
kesize,all,skind/2
lsel,s,loc,x,x11/2  
lesize,all,,,1
lsel,all
mshape,0,2d
mshk,1
amesh,9
asel,s,area,,9
allsel,below,area
nsel,r,loc,y,y2
sf,all,rdsf,0.9,1         ! Radiation
alls
spctemp,1,25               ! Space temperature
v2dopt,1
radopt,,0.01
stef,5.67e-8               ! Stefan-Boltzman constant
esel,s,mat,,3              ! select shimo material
bfe,all,fvin,,1            ! define volumetric interface
finish
/solu
mfan,on                    ! Activate ANSYS Multi-field solver analysis
mfel,1,1,2                 ! Field #1 ET;s,  Emag
mfel,2,4                   ! Field #2 ET's, Thermal
mfor,1,2                   ! Field solution order
mfti,ftime                 ! Final time
mfdt,tinc                  ! Stagger time increment
mfco,all,1e-2              ! Convergence criteria
antyp,harm                 ! Emag analysis options
harfrq,30000
outres,all,all
tunif,100
mfcm,1,                    ! Write Emag analysis options
mfclear,solu               ! Clear analysis options
antype,trans               ! Thermal analysis options
toffst,273
tunif,100                  ! initial uniform temperature
kbc,1                      ! step loads
trnopt,full
autos,on                   ! auto time-stepping
deltim,0.01,0.005,0.01,on     ! time step control
mfcm,2,                    ! Write Thermal analysis options
mfvo,1,1,hgen,2            ! Transfer hgen from Emag to Thermal
mfvo,1,2,temp,1            ! Transfer Temp from Thermal to Emag
solve
save
finish
/post1
file,field2,rth
set,last                   ! Solution at 1 seconds
esel,s,type,,4             ! select thermal elements
plns,temp                  ! plot temperature
finish