killna 发表于 2009-3-30 16:30:25

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

我作的是电磁感应加热,可是在计算时出现了这样的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 发表于 2009-3-31 14:02:00

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

xbzhang 发表于 2009-3-31 20:09:19

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

killna 发表于 2009-4-1 16:05:58

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

xbzhang 发表于 2009-4-1 19:43:04

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

killna 发表于 2009-4-3 17:59:50

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

xbzhang 发表于 2009-4-4 07:32:11

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

quhaiduan 发表于 2009-7-23 14:17:23

楼主
你的问题解决了吗?
如果解决了的话
能把你的命令流发到我邮箱里吗》
qu717@163.com
页: [1]
查看完整版本: 感应加热出现无.rst文件可用的问题