请教：感应加热分析时出现的错误

twt716

我做一个感应加热的分析，是参照ansys帮助里面的那个例子做的，只把模型做了稍微的修改，别的都和例子一样，但在算到热环境的求解的时候出现了下面的错误，“ Element 1063 has an undefined node number 0.  ”一共有38个，我将这些element画出来看了，就些结点就是在定义的热分析表面（SURF151单元）上的结点，请教下大家这问题是怎么回事？

下面是我的命令流，大家帮助看看，谢谢！

/batch,list
/filenam, induc-heating
/prep7
shpp,off
/title,induction heating
/com,
et,1,13,,,1     !plane13,axisymmetric, AZ dof,gang
et,2,13,,,1     !kongqi
et,3,151,,,1,1,1        ! SURF151, thermal, radiation
r,3,0
r1=0.1
r2=0.115
r3=0.165
r4=0.245
r5=0.7
h1=.05
h2=.15
freq=50
pi=4*atan(1)
cond=.392e7
muzero=4e-7*pi       ! free-space permeability
mur=200              ! maximum relative permeability
skind=sqrt(1/(pi*freq*cond*muzero*mur))    ! skin depth
ftime=3
tinc=.05
time=0
delt=.01

emunit,mks    ! set magnetic units
mp,murx,1,1             ! air relative permeability
mp,murx,3,1             ! coil relative permeability
mptemp,1,25.5,160,291.5,477.6,635,698   ! temps for relative permeability
mptemp,7,709,720.3,742,761,1000
mpdata,murx,2,1,200,190,182,161,135,104 ! steel relative permeability
mpdata,murx,2,7,84,35,17,1,1
mptemp
mptemp,1,0,125,250,375,500,625          ! temps for resistivity
mptemp,7,750,875,1000
mpdata,rsvx,2,1,.184e-6,.272e-6,.384e-6,.512e-6,.656e-6,.824e-6
mpdata,rsvx,2,7,1.032e-6,1.152e-6,1.2e-6            ! steel resistivity

rectng,0,r1,0,h2        ! billet
rectng,r1,r2,0,h1      ! air-gap
rectng,r3,r4,0,h1      ! coil
rectng,r1,r5,0,h2       ! outer air
aovlap,all
numcmp,area
asel,s,area,,1      ! set attributes for coil region
aatt,3,1,2
asel,s,area,,2       ! set attributes for billet region
aatt,2,1,1
asel,s,area,,3
aatt,2,1,1
asel,s,area,,4        ! set attributes for air region
aatt,1,1,2
asel,all
lsel,s,loc,x,r3,r4
lesize,all,,,8
lsel,all
mshape,0,2d
mshk,1
amesh,1
lsel,s,loc,x,r2
lesize,all,skind/2
lsel,s,loc,x,0
lesize,all,40*skind
lsel,all
lsel,s,loc,x,r5
lesize,all,.12
lsel,all
mshape,1,2d
mshk,0
amesh,all
n                       ! create space node for SURF151   
*get,nmax,node,,num,max
lsel,s,loc,x,r2
type,3
real,3
mat,2
lmesh,all            ! mesh billet outer radius with SURF151      
*get,emax,elem,,num,max
emodif,emax,3,nmax   ! modify element to add space node for radiation
et,3,0               ! reset type 3 to null element

nsel,s,loc,x,0
d,all,az,0              ! apply flux-normal b.c.
nsel,s,loc,x,r5
d,all,az,0              ! apply flux-normal b.c.
nsel,all
esel,s,mat,,3
bfe,all,js,,,,15e6      ! apply current density to coil
esel,all
finish
/solu
antyp,harm
harfrq,50
physics,write,emag      ! write emag physics file
finish

/prep7
lsclear,all             ! clear all b.c.'s and options
et,1,55,,,1             ! PLANE55 thermal element, axisymmetric
et,2,0                  ! null element type for coil and air region
et,3,151,,,1,1,1        ! SURF151 element for radiation
keyopt,3,9,1
r,3,1,5.67e-8           ! form factor, Stefan-Boltzmann constant
mptemp
mptemp,1,0,730,930,1000         ! temps for conductivity
mpdata,kxx,2,1,60.64,29.5,28,28
mptemp                          ! temps for enthalpy
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.548547e9,3.548556e9,4.3520e9
mp,emis,2,.68                   ! emissivity
finish
/solu
antype,trans
toffst,273
tunif,30                       ! initial uniform temperature
d,nmax,temp,25.5                  ! ambient temperature
cnvtol,heat,1                   ! convergence tolerance
kbc,1                           ! step loads
trnopt,full
autos,on                        ! auto time-stepping
deltim,1e-5,1e-6,delt,on        ! time step control
outres,basic,all                ! save all load step information
physics,write,thermal           ! write thermal physics file
finish
*do,i,1,ftime/tinc   ! solution *do loop
time=time+tinc               ! increment time
physics,read,emag    ! read emag physics file
/solu
*if,i,eq,1,then
tunif,25.5                    ! initial temperature
*else
ldread,temp,last,,,,,rth     ! read thermal analysis temperatures
*endif
solve                         ! solve harmonic analysis
finish
physics,read,thermal          ! read thermal physics file
/assign,esav,therm,esav       ! redirect files for use in thermal restart
/assign,emat,therm,emat
/solu
*if,i,gt,1,then
antype,trans,rest            ! thermal restart
*endif
time,time                     ! time at end of thermal run
esel,s,mat,,2                 ! select billet region
ldread,hgen,,,,2,,rmg         ! apply coupled joule heating load from emag
esel,all
solve

twt716

没有人回答吗，求求大家，我刚接触ansys不久，很急，希望大家能帮帮忙，谢谢