做磁热耦合遇到了问题，请大家帮帮忙

我爱北京

这是一个线圈加热钢管的问题。总是提示有错误，我觉得是创建辐射单元时有问题，可是我水平有限，不知道该怎样解决。请大家帮忙看一下吧。程序流如下：
finish
/clear
/prep7

et,1,13,,,1
et,2,13,,,1
et,3,151,,,1,1,1
r,3,0

!参数输入
r0=9e-3  !钢管半径
t0=1e-3   !管厚
d0=5e-3   !管与线圈间的距离
h1=10e-3   !线圈的宽度
l1=10e-3    !线圈的长度
l=30e-3   !外层空气



freq=150000 !频率
PI=4*ATAN(1)
COND=0.392E7 !最大电导率
MUZERO=4E-7*PI !真空磁导率
MUR=200   !最大相对磁导率
SKIND=SQRT(1/(PI*FREQ*COND*MUZERO*MUR)) !趋肤深度
FTIME=3   !结束时间
TINC=0.05 !谐波分析时间增量
TIME=0     !起始时间
DELT=0.01   !最大时间步
emunit,mks   !设置磁场单位
!******材料属性**********
MP,MURX,1,1
MP,MURX,3,1
MPTEMP,,,,,,,,
!钢的相对磁导率与温度的关系
MPTEMP,1,25.5,160,291.5,477.6,635,698
MPTEMP,7,709,720.3,742,761,1000
MPDATA,MURX,2,1,200,190,182,161,135,104
MPDATA,MURX,2,7,84,35,17,1,1
MPTEMP,,,,,,,,
!钢的电阻率与温度的关系
mptemp,1,0,125,250,375,500,625
MPTEMP,7,750,875,1000
MPDATA,RSVX,2,1,0.184E-6,0.272E-6,0.384E-6,0.512E-6,0.656E-6,0.824E-6
MPDATA,RSVX,2,7,1.032E-6,1.152E-6,1.2E-6
!*********建模********
rectng,0,r0,0,h1
rectng,r0,r0+t0,0,h1
rectng,r0+t0,r0+t0+d0,0,h1
rectng,r0+t0+d0,r0+t0+d0+l1,0,h1
rectng,r0+t0+d0+l1,r0+t0+d0+l1+l,0,h1
aglue,all
numcmp,area
!***********赋材料属性**********
asel,s,,,2   !tube
aatt,2,,1
asel,s,,,4   !coils
aatt,3,,2
asel,s,,,1,5,2  !air
aatt,1,,2
allsel
!*********分网*************
mshape,0,2d
mshk,1
amesh,2    !tube
lsel,s,loc,y,0
lsel,a,loc,y,h1
lsel,u,loc,x,r0+t0/2
lesize,all,0.001
allsel
amesh,all
!***********创造辐射单元*************
n
*GET,NMAX,NODE,,NUM,MAX
LSEL,S,,,2,6,4
TYPE,3
REAL,3
MAT,2
LMESH,ALL
*GET,EMAX,ELEM,,NUM,MAX
EMODIF,EMAX,3,NMAX
ET,3,0
!磁力平行
nsel,s,loc,x
d,all,az,0
nsel,all
esel,s,mat,,3
bfe,all,js,,,,15e6 !电流密度
esel,all
allsel
antype,harmnic
harfrq,150000
physics,write,emag
!*******以上为磁分析************
!******转入热分析*************
lsclear,all   !清除所有的边界条件和选项
et,1,55,,,1   !热单元 轴对称
et,2,0    !线圈和空气区为0号单元
ET,3,151,,,1,1,1   !辐射单元
KEYOPT,3,9,1
R,3,1,5.67E-8   !波尔兹曼实常数
!*********材料属性**************
mptemp
!钢的热传导率与温度的关系
MPTEMP,1,0,730,930,1000
MPDATA,KXX,2,1,60.64,29.5,28,28
mptemp
!钢的焓与温度的关系
MPTEMP,1,0,27,127,327,527,727
MPTEMP,7,765,765.001,927
MPDATA,ENTH,2,1,0,91609056,453285756,1.2748E9,2.2519E9,3.3396E9
MPDATA,ENTH,2,7,3.548556E9,4.3520E9
!辐射率
MP,EMIS,2,0.68
ANTYPE,TRANS  !热求解选项
TOFFST,273
TUNIF,100   !初始归一化温度
D,NMAX,TEMP,20  !环境温度
cnvtol,heat,1  !收敛公差
KBC,1           !载荷步
trnopt,FULL
AUTOS,ON   !自动时间一步
deltim,1E-5,1E-6,delt,ON !时间步控制
OUTRES,BASIC,ALL !保存所有载荷步信息
physics,write,thermal !写入THERMAL物理文件

!*****循环**********
*DO,i,1,FTIME/TINC
TIME=TIME+TINC !时间增量
PHYSICS,READ,EMAG  !读入EMAG物理文件
/SOLU
*IF,i,EQ,1,THEN
TUNIF,100   !初始化温度
*ELSE
LDREAD,TEMP,LAST,,,,,RTH !读入热分析温度
*ENDIF
SOLVE !求解
FINISH
PHYSICS,READ,THERMAL !读入THERMAL物理文件
/ASSIGN,ESAV,THERM,ESAV
/ASSIGN,EMAT,THERM,EMAT !REDIRECT FILES FOR USE IN THERMAL RESTART
/SOLU
*IF,I,GT,1,THEN
ANTYPE,TRANS,REST !热重启动
*ENDIF
TIME,TIME   !热分析结束时间
ESEL,S,MAT,,2  !选择钢管区域
LDREAD,HGEN,,,,2,,RMG !从EMAG分析中读入焦耳热载荷
ALLSEL
SOLVE
FINISH
/ASSIGN,ESAV
/ASSIGN,EMAT
*ENDDO
FINISH
SAVE
/POST26
/SHOW
NSOL,2,1,TEMP,,TEMPCL
NSOL,3,2,TEMP,,TEMPSURF
PLVAR,2,3
PRVAR,2,3
FINISH