感应加热出现无.rst文件可用的问题

killna

我作的是电磁感应加热，可是在计算时出现了这样的error.
No results file ( xx.rst ) available. The LDREAD command is ignored.
这是为什么，请大侠帮忙解决一下。

/BATCH  
/CPLANE,1  
/CLEAR,NOSTART
/VIEW,1,1,1,1  
WPSTYLE,,,,,,,,0
/config,nres,3000000
/Nerr,1,1500000
/filenam,induction1
  
/prep7
/title,Induction Heating
et,1,117                   !non-steel region   
et,2,117,1                 !steel region
et,3,117,1              !coil

*SET,freq,40000               ! frequency (Hz.)   
*SET,pi,4*atan(1)              ! pi   
*SET,cond,.392e7               ! maximum conductivity
*SET,muzero,4e-7*pi            ! free-space permeability  
*SET,mur,200                   ! maximum relative permeability
*SET,skind,sqrt(1/(pi*freq*cond*muzero*mur)) ! skin depth   
*SET,ftime,20                  ! final time   
*SET,tinc,1                    ! time increment for harmonic analysis
*SET,time,0                    ! initialize time  
*SET,delt,1                    ! maximum delta time step  

emunit,mks               ! set magnetic units   
mp,murx,1,1                ! air relative permeability
mptemp  
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,20,100,200,400,600             ! temps for resistivity
mpdata,rsvx,2,1,2.23e-7,2.71e-7,3.42e-7,5.29e-7,7.85e-7
mptemp  
mptemp,1,0,100,300,500,1200             ! temps for conductivity
mpdata,kxx,2,1,42.7,42.7,40.6,37.3,30.1
mp,dens,2,7850   
mp,c,2,465  
mp,dens,3, 8960   
mp,RSVX,3,1.7e-8   
mp,murx,3,1

BLOCK,-10.5/1000,-0.5/1000,-35/1000,35/1000,0/1000,80/1000   
BLOCK,0.5/1000,10.5/1000,-35/1000,35/1000,0/1000,80/1000  
WPAVE,10.5/1000,50/1000,7.5/1000
wpro,,,90
CYL4, , ,4/1000
CYL4, , ,3/1000
ASBA,      13,      14
wpro,,,-90
K,38,25/1000,50/1000,7.5/1000
K,39,25/1000,-50/1000,7.5/1000
K,40,-25/1000,-50/1000,7.5/1000
K,41,-25/1000,50/1000,7.5/1000
L,38,39
L,39,40
L,40,41
L,41,38
LFILLT,36,33,15/1000,
LFILLT,33,34,15/1000,
LFILLT,34,35,15/1000,
LFILLT,35,36,15/1000,
FLST,8,8,4  
FITEM,8,37  
FITEM,8,33  
FITEM,8,38  
FITEM,8,34  
FITEM,8,39  
FITEM,8,35  
FITEM,8,40  
FITEM,8,36  
VDRAG,      15, , , , , ,P51X
alls
WPAVE,0,0,0
WPSTYLE,,,,,,,,0
FLST,3,8,6,ORDE,2   
FITEM,3,3   
FITEM,3,-10
VGEN,5,P51X, , , , ,16/1000, ,0
BLOCK,0.2,-0.2,0.16,-0.16,0.18,-0.1
alls
VOVLAP,ALL  
VGLUE,ALL
NUMCMP,AREA
NUMCMP,VOLUME

FLST,5,40,6,ORDE,2  
FITEM,5,3   
FITEM,5,-42
VSEL,S, , ,P51X
VATT,3,1,3
FLST,5,2,6,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
VSEL,S, , ,P51X
VATT,2,1,2
FLST,5,6,6,ORDE,2   
FITEM,5,43  
FITEM,5,-48
VSEL,S, , ,P51X
VATT,1,1,1
alls
FLST,5,2,6,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
VSEL,S, , ,P51X
FLST,5,4,4,ORDE,4   
FITEM,5,2   
FITEM,5,7   
FITEM,5,16  
FITEM,5,-17
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,20,10, , , ,1
FLST,5,4,4,ORDE,4   
FITEM,5,4   
FITEM,5,-5  
FITEM,5,14  
FITEM,5,19  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,20,0.1, , , ,1   
FLST,5,16,4,ORDE,9  
FITEM,5,1   
FITEM,5,3   
FITEM,5,6   
FITEM,5,8   
FITEM,5,-13
FITEM,5,15  
FITEM,5,18  
FITEM,5,20  
FITEM,5,-24
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,30, , , , ,1
FLST,5,2,6,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
CM,_Y,VOLU  
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VSWEEP,_Y1  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2
FLST,5,40,6,ORDE,2  
FITEM,5,3   
FITEM,5,-42
VSEL,S, , ,P51X
vsweep,all
alls
FLST,5,5,6,ORDE,2   
FITEM,5,43  
FITEM,5,-47
VSEL,S, , ,P51X
SMRT,6
MSHAPE,1,3D
MSHKEY,0
VMESH,all
alls
CM,_Y,VOLU  
VSEL, , , ,      48
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2
alls
VPLOT
SAVE
FINISH
/SOLU            
ASEL,S, EXT
DA,ALL,AZ,0
ASEL,ALL
FLST,5,5,6,ORDE,2   
FITEM,5,38  
FITEM,5,-42
VSEL,S, , ,P51X
eslv
esln,r,1
bfe,all,js,1,4.55e6,,
FLST,5,5,6,ORDE,5   
FITEM,5,6   
FITEM,5,13  
FITEM,5,20  
FITEM,5,27  
FITEM,5,34  
VSEL,S, , ,P51X
eslv
esln,r,1
bfe,all,js,1,-4.55e6,,
FLST,5,15,6,ORDE,10
FITEM,5,3   
FITEM,5,-5  
FITEM,5,10  
FITEM,5,-12
FITEM,5,17  
FITEM,5,-19
FITEM,5,24  
FITEM,5,-26
FITEM,5,31  
FITEM,5,-33
VSEL,S, , ,P51X
eslv
esln,r,1
BFE,ALL,JS, 1, ,-4.55e6,
FLST,5,15,6,ORDE,10
FITEM,5,7   
FITEM,5,-9  
FITEM,5,14  
FITEM,5,-16
FITEM,5,21  
FITEM,5,-23
FITEM,5,28  
FITEM,5,-30
FITEM,5,35  
FITEM,5,-37
VSEL,S, , ,P51X
eslv
esln,r,1
BFE,ALL,JS, 1, ,4.55e6,
ALLS
/solu   
antyp,harm  
harfrq,40000
HROPT,FULL  
HROUT,ON
EQSLV, ,1e-008,
physics,write,emag      ! write emag physics file
solve
SAVE  
finish

/prep7  
lsclear,all             ! clear all b.c.'s and options  
et,1,0                  ! non-steel
et,2,SOLID90                  ! steel
et,3,0
finish  
/solu   
antype,trans
toffst,273  
tunif,50                       ! initial uniform temperature
cnvtol,heat,1                   ! convergence tolerance
kbc,1                           ! step loads
trnopt,full
autos,on                        ! auto time-stepping
deltim,1e-1,1e-2,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
*SET,time,time+tinc             ! increment time
physics,read,emag               ! read emag physics file
/solu   
*if,i,eq,1,then
tunif,50                   ! 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,,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
save
finish

killna

做磁热耦合分析的都没有遇过这个错误吗?
该如何解决，请大家帮忙解救
谢谢

xbzhang

写成rth，试试
以前好像有人问过，直接用rth

killna

磁分析结果为rst,再把此结果带入分析热,副档名为rth,您的意思是磁分析直接写成rth,那我热分析不就要写成rst,印象中记得不同的副档名有限定储存分析的领域范围,热分析就要写成rth,这样好像就会与您提出的想法出现冲突,或是我误会您的意思
我有找寻之前相同问题发问的文章,有人说是单元体的问题,但我用相同的单元体在别的case上是可以分析,偏偏在我这case是不行,因为那篇是前年的问题,想问问板上大侠,此问题有新的解决发法或是有新的见解,能提供给我参考,谢谢

xbzhang

哦，看错了，后一个磁场的结果应该是rmg文件吧，你看你的盘里应该有rmg。
ansys理论上，磁场分析的结果是rmg，热场是rth文件。
以前有人问过没有rmg文件，在论坛上讨论，后来也没有准确的解释。没有rmg就用rst了。

killna

我只有rst,但是它档案为0MB,这个错误要怎解决阿
很奇妙的错误,有时候可以但改变模型就会有这问题

xbzhang

我觉得你还是重新建模吧，换一种建模方法，仔细检查各项参数，大体的磁热耦合过程就是这样的。有限元分析急不得。
你的命令流好多都是交互的，需要正规的命令流。

quhaiduan

楼主
你的问题解决了吗？
如果解决了的话
能把你的命令流发到我邮箱里吗》
qu717@163.com