3d运动钢坯感应加热 Error occurred in tree gauging for edge formulation

fangdong116

运动的钢坯感应加热例子，出现如题错误，望大牛指教
/title 3-d billet induction heating
/filname 3-d coupled field analysis
/prep7
blen=3.0        !钢坯总长
coil=1.0        !线圈长度
tot=blen+coil   !模型总长
div=100         !划分段数
bsidlen=0.06   !钢坯边长
clen=blen/div   !每段钢坯长度
cinsidlen=0.115          !线圈内边长/2
cthic=0.025        !线圈厚度
coutsidlen=cinsidlen+cthic  !线圈外边长
ra=2*coutsidlen            !空气范围
f=1000
!pi=3.1415926
!cond=0.392e7
!muzero=4e-7*pi
!mur=200
!skind=sqrt(1/(pi*f*cond*muzero*mur))

vol=5                  !钢坯运动速度

ftime=60*tot/vol      !结束时间
tinc=60*clen/vol      !载荷步长
tim=0                 !开始时间
/view,,1,1,1
et,1,117,1              !钢坯
et,2,117,1              !线圈
et,3,117              !空气
emunit,mks
mptemp,1,25.5,477.6,635,698,709,720.3    !钢随温度的相对磁导率
mptemp,7,742,761,1000,1100,1200
mpdata,murx,1,1,200,161,135,104,84,35
mpdata,murx,1,7,17,1,1,1,1
mptemp
mptemp,1,20,400,500,600,700,800          !钢的电阻率
mptemp,7,900,1000,1100,1200,1300
mpdata,rsvx,1,1,0.16e-6,0.483e-6,0.610e-6,0.755e-6,0.922e-6,1.095e-6
mpdata,rsvx,1,7,1.135e-6,1.168e-6,1.195e-6,1.220e-6,1.240e-6 ! steel resistivity
mp,murx,2,1           !线圈磁导率
mp,rsvx,2,2e-8        !线圈电阻率
mp,murx,3,1           !空气磁导率

cyl4,0,0,0,0,ra,90,blen
block,0,coutsidlen,0,coutsidlen,0,blen
vsbv,1,2
numcmp,all
block,0,coutsidlen,0,coutsidlen,0,clen
block,0,cinsidlen,0,cinsidlen,0,clen
vsbv,2,3
numcmp,all
vgen,100,2,,,0,0,clen
vglue,all
numcmp,all
block,0,cinsidlen,0,cinsidlen,0,blen
block,0,bsidlen,0,bsidlen,0,blen
vovlap,102,103
vglue,all
numcmp,all             !线圈编号1-100，外围空气编号101，内空气编号102，钢坯编号103
vsel,s,,,103       !选择钢坯体
vatt,1,1,1
vsel,s,,,101,102     !选择空气体
vatt,3,1,3
vsel,s,,,1,100        !选择线圈体
vatt,2,1,2      
vsel,all
mshape,0,3d
mshkey,1
!划分钢坯
lesize,1228,,,10,1/4         !宽划分
lesize,1222,,,10,1/4         !厚度划分
lesize,1223,,,10,4            !宽划分
lesize,1230,,,10,4            !厚度划分
lesize,1225,,,200
vsweep,103
mshape,0,3d                  !单元形状为六面体
mshk,1                       !映射网格划分
vsel,all
vsel,u,,,103
esize,0.05
vsweep,all
alls
! set ground volt=0
*do,k,1,100,1
        vsel,s,,,k
        aslv
        asel,r,loc,x,0
        nsla,s,1
        d,all,volt,0
        nsel,all
        !cp source node
        vsel,s,,,k
        aslv
        asel,r,loc,y,0
        nsla,s,1
        cp,1,volt,all
*enddo
! short curcuit conductor: set volt=0 on symmetry plane
vsel,s,,,103
aslv
asel,r,loc,x,0
nsla,s,1
d,all,volt,0
nsel,all
vsel,s,,,103
aslv
asel,r,loc,y,0
nsla,s,1
d,all,volt,0
nsel,all

asel,s,loc,x,0
asel,a,loc,y,0
nsla,s,1
d,all,az,0
nsel,all

csys,1
asel,s,loc,x,ra
nsla,s,1
d,all,az,0
csys,0

nsel,all
esel,all
finish
/solu
antype,harm
harfrq,1000
physics,write,emag      ! write emag physics file
finish

/prep7
lsclear,all             ! clear all b.c.'s and options
et,1,90                 !solid90 thermal element,axisymmetic
et,2,0
et,3,0

mptemp
mptemp,1,0,730,930,1000,1200         !钢的导热系数
mpdata,kxx,1,1,60.64,29.5,28,28,28
mptemp                               !钢的焓值
mptemp,1,0,27,127,327,527,727
mptemp,7,765,765.001,927,1099,1225
mpdata,enth,1,1,0,91609056,453285756,1.2748e9,2.2519e9,3.3396e9
mpdata,enth,1,7,3.548547e9,3.548556e9,4.3520e9,5.0984e9,5.6990e9
!mp,emis,1,0.68                      !辐射率
asel,s,loc,y,bsidlen           !在钢坯侧表面辐射
asel,a,loc,x,bsidlen
nsla,s,1
sf,all,rdsf,0.68,1

vsel,s,,,1,100               !在钢坯后表面加辐射系数
aslv
asel,r,loc,z,0
sf,all,rdsf,0.68,1

vsel,s,,,1,100                !在钢坯前表面加辐射系数
aslv
asel,r,loc,z,blen
sf,all,rdsf,0.68,1
alls
finish

/solu
antype,trans              
toffst,273               !绝对温度偏移量
tunif,800              !初始温度
spctemp,1,25        !环境温度
radopt,,0.01               !塞舌尔系数
stef,5.67e-8               ! Stefan-Boltzman constant
cnvtol,heat,1                   ! convergence tolerance
kbc,1                           ! step loads
trnopt,full
autos,on                        ! auto time-stepping
nsubst,10
outres,basic,all                ! save all load step information
physics,write,thermal           ! write thermal physics file
finish


*do,i,0,ftime,tinc  ! solution *do loop
tim=tim+tinc                             ! increment time
physics,read,emag    ! read emag physics file
/solu
*if,i,eq,0,then
tunif,800                    ! initial temperature
*else
ldread,temp,last,,,,,rth     ! read thermal analysis temperatures
*endif
!加载电流
*if,i,le,60*coil/vol,then          !钢坯正在进入线圈，载荷逐渐施加
        *do,k,1,tim/tinc
        vsel,s,,,k
        aslv
        asel,r,loc,y,0
        nsla,s,1
        f,ndnext(0),amps,10000
        *enddo
*elseif,i,gt,60*blen/vol,then !钢坯正在退出线圈，载荷逐渐删除
        *do,k,tim/tinc+1,100
        vsel,s,,,k
        aslv
        asel,r,loc,y,0
        nsla,s,1
        f,ndnext(0),amps,10000
        *enddo
*else                      !钢坯在线圈当中，载荷交替
       
        *do,k,tim/tinc,tim/tinc+coil/clen
        vsel,s,,,k
        aslv
        asel,r,loc,y,0
        nsla,s,1
        f,ndnext(0),amps,10000
        *enddo
*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
parsave,all
*if,i,gt,0,then
antype,trans,rest            ! thermal restart
*endif
parresume
time,tim                    ! time at end of thermal run
esel,s,mat,,1                 ! select billet region
ldread,hgen,,,,2,,rst        ! apply coupled joule heating load from emag
esel,all
solve
finish
/assign,esav                  ! reassign files to default
/assign,emat
*enddo                        ! end of solution looping
finish
save

/post1
set,last                   ! Solution at last seconds
esel,s,type,,1             ! select thermal elements
/expand,4,polar,half,,90
/AUTO,1
/REP,FAST
plns,temp                  ! plot temperature