请教：感应加热（间接耦合）问题

banana_wly

描述：单相交流电源的导线，穿过空心杆，并与杆在末端相连，从而形成一闭合串联电路。当回路中有交变电流通过时，导线周围将有交变磁场产生。由于空心杆属铁磁性物质(可视为一铁芯)，当其所处的磁场强度发生变化时，通过空心杆轴截面的磁通也随之变化，此时便有感生电动势产生，其大小和方向可由麦克斯韦电磁感应定律来确定。在空心杆的外壁附近，由于整个空心杆是闭合的，所以其中的感生电动势在垂直于导线的磁力线的平面上就会产生感生电流，即电涡流，它存在于整个抽油杆的轴截面内。空心杆本身的电阻增大了，他就成为发热体。
我想知道抽油杆被加热时的温度变化。
我是仿照ansys,help中2.10. Example Induction-heating Analysis Using Physics Environments做的。结果没有得出结果，老是死循环。请高手帮忙解决，万分感激！

banana_wly

/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Thermal  
/COM,  Magnetic-Nodal   
!*  
/PREP7  
!*  
ET,1,PLANE13
!*  
KEYOPT,1,1,6
KEYOPT,1,2,0
KEYOPT,1,3,0
KEYOPT,1,4,0
KEYOPT,1,5,0
!*  
ET,2,PLANE13
!*  
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,0
KEYOPT,2,4,0
KEYOPT,2,5,0
!*  
ET,3,PLANE13
!*  
KEYOPT,3,1,6
KEYOPT,3,2,0
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,5,0
!*  
ET,4,SURF151
!*  
KEYOPT,4,1,0
KEYOPT,4,2,0
KEYOPT,4,3,0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,6,0
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,0
!*  
!*  
!*  
R,4, , , , , ,  
RMORE, , , , , ,
RMORE, , ,  
!*  
!*  
!*  
EMUNIT,MKS  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,1,,1                  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,RSVX,1,,1.678e-8
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MURX,2,,1
MPTEMP,,,,,,,,  
MPTEMP,1,0
MPTEMP,2,25.5   
MPTEMP,3,160
MPTEMP,4,291.5  
MPTEMP,5,477.6  
MPTEMP,6,635
MPTEMP,7,698
MPTEMP,8,709
MPTEMP,9,720.3  
MPTEMP,10,742
MPTEMP,11,761   
MPTEMP,12,1000
MPDATA,MURX,3,,205
MPDATA,MURX,3,,200  
MPDATA,MURX,3,,190  
MPDATA,MURX,3,,182  
MPDATA,MURX,3,,161  
MPDATA,MURX,3,,135  
MPDATA,MURX,3,,104  
MPDATA,MURX,3,,84   
MPDATA,MURX,3,,35   
MPDATA,MURX,3,,17   
MPDATA,MURX,3,,1
MPDATA,MURX,3,,1
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPTEMP,2,125
MPTEMP,3,250
MPTEMP,4,375
MPTEMP,5,500
MPTEMP,6,625
MPTEMP,7,750
MPTEMP,8,875
MPTEMP,9,1000   
MPDATA,RSVX,3,,0.184e-6
MPDATA,RSVX,3,,0.272e-6
MPDATA,RSVX,3,,0.384e-6
MPDATA,RSVX,3,,0.512e-6
MPDATA,RSVX,3,,0.656e-6
MPDATA,RSVX,3,,0.824e-6
MPDATA,RSVX,3,,1.032e-6
MPDATA,RSVX,3,,1.152e-6
MPDATA,RSVX,3,,1.2e-6   
SAVE
PCIRC,0.003, ,0,360,            !cable
PCIRC,0.012, ,0,360,           !inner radius of guan   
PCIRC,0.018, ,0,360,         !outer radius  of guan  
PCIRC,0.060, ,0,360,         !outer radius of model
FLST,2,4,5,ORDE,2   
FITEM,2,1   
FITEM,2,-4  
AOVLAP,P51X
SAVE
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1,       4,   1,       0,   !copper  cable
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,5   
FITEM,5,7   
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2,       4,   2,       0,   !air
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       6
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       3,       4,   3,       0,   !guan
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
SMRT,6                      !mesh all
SMRT,2  
MSHAPE,0,2D
MSHKEY,0
!*  
FLST,5,4,5,ORDE,3   
FITEM,5,1   
FITEM,5,5   
FITEM,5,-7  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FLST,2,4,5,ORDE,3   
FITEM,2,1   
FITEM,2,5   
FITEM,2,-7  
ACLEAR,P51X
APLOT   
MSHAPE,1,2D
!*  
FLST,5,4,5,ORDE,3   
FITEM,5,1   
FITEM,5,5   
FITEM,5,-7  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   '
CMDELE,_Y1  
CMDELE,_Y2  
!*  
APLOT   
EPLOT   
SAVE
n                       ! create space node for SURF151   
*get,nmax,node,,num,max

FLST,5,4,4,ORDE,2   
FITEM,5,5   
FITEM,5,-8  
LSEL,S, , ,P51X
NSLL,S,1
/UI,MESH,OFF
ALLSEL,ALL  
FLST,5,4,4,ORDE,2   
FITEM,5,5   
FITEM,5,-8  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,S,_Y  
!*  
!*  
CMSEL,S,_Y1
LATT,3,4,4, , , ,   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
ALLSEL,ALL  
FLST,2,4,4,ORDE,2   
FITEM,2,5   
FITEM,2,-8  
LMESH,P51X  
SAVE
*get,emax,elem,,num,max
emodif,emax,4,nmax   ! modify element to add space node for radiation   
EPLOT   
et,4,0               ! reset type 4 to null element
/UI,MESH,OFF

   

APLOT   
ALLSEL,ALL  
ASEL,S, , ,       1
NSLA,S,1
FLST,4,49,1,ORDE,6  
FITEM,4,1   
FITEM,4,-4  
FITEM,4,17  
FITEM,4,-48
FITEM,4,185
FITEM,4,-197
CP,1,VOLT,P51X                  !couple volt in cupper cable
FLST,2,49,1,ORDE,6  
FITEM,2,1   
FITEM,2,-4  
FITEM,2,17  
FITEM,2,-48
FITEM,2,185
FITEM,2,-197
/GO
!*  
F,P51X,AMPS,160              !apply current in copper cable
ALLSEL,ALL  
ASEL,S, , ,       6
NSLA,S,1
FLST,4,148,1,ORDE,8
FITEM,4,5   
FITEM,4,-12
FITEM,4,49  
FITEM,4,-116
FITEM,4,481
FITEM,4,-516
FITEM,4,1124
FITEM,4,-1159   
CP,2,VOLT,P51X                          !couple volt in guan
FLST,2,148,1,ORDE,8
FITEM,2,5   
FITEM,2,-12
FITEM,2,49  
FITEM,2,-116
FITEM,2,481
FITEM,2,-516
FITEM,2,1124
FITEM,2,-1159   
/GO
!*  
F,P51X,AMPS,-160                     !aplly current in guan
ALLSEL,ALL  
FINISH  
/SOL
!*  
ANTYPE,3
HARFRQ,50,0,
NSUBST,0,   
KBC,0   
!*  
FINISH  
/PREP7  
PHYSICS,WRITE,emage, , ,
FINISH  
/prep7  
lsclear,all             ! clear all b.c.'s and options  
et,1,0  
et,2,0  
et,3,55,,,1
et,4,151,,,1,1,1        ! SURF151 element for radiation
keyopt,4,9,1
r,4,1,5.67e-8           ! form factor, Stefan-Boltzmann constant
SAVE
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPTEMP,2,730
MPTEMP,3,930
MPTEMP,4,1000   
MPDATA,KXX,3,,60.64
MPDATA,KXX,3,,29.5  
MPDATA,KXX,3,,28
MPDATA,KXX,3,,28.   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPTEMP,2,27
MPTEMP,3,127
MPTEMP,4,327
MPTEMP,5,527
MPTEMP,6,727
MPTEMP,7,765
MPTEMP,8,765.001
MPTEMP,9,927
MPDATA,ENTH,3,,0
MPDATA,ENTH,3,,91609056
MPDATA,ENTH,3,,453285756
MPDATA,ENTH,3,,1.2748e9
MPDATA,ENTH,3,,2.2519e9
MPDATA,ENTH,3,,3.3396e9
MPDATA,ENTH,3,,3.548547e9   
MPDATA,ENTH,3,,3.548556e9   
MPDATA,ENTH,3,,4.3520e9
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EMIS,3,,0.68

SAVE
FINISH  
/PREP7  
FINISH  
/SOLU
!*  
antype,trans
!*  

TOFFST,273,
!*  
TUNIF,18,   
cnvtol,heat,1                   ! convergence tolerance
kbc,1                           ! step loads
trnopt,full
autos,on                        ! auto time-stepping
*SET,ftime,12    ! final time   
*SET,tinc,.05                ! time increment for harmonic analysis
*SET,time,0                  ! initialize time  
*SET,delt,.01                ! maximum delta time step  
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
*SET,time,time+tinc               ! increment time  
SAVE
physics,read,emag    ! read emag physics file   
/solu   
*if,i,eq,1,then
tunif,18
*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,,3                 ! select billet region
ldread,hgen,,,,2,,rmg         ! 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                          ! save database
/post26                       ! time-history postprocessor
/show
nsol,2,1132,temp,,tempcl         ! store temperature at billet centerline
nsol,3,85,temp,,tempsurf       ! store temperature at billet outer diameter
plvar,2,3                     ! plot temperature rise over time
prvar,2,3                     ! print temperature rise over time
finish

perch88

别的不说，搜索一下
et,4,151,,,1,1,1
相当于你最后把151改成了轴对称，这个肯定不对，有可能是死循环的原因。

从此而引发了对大家的建议：看例题要小心仔细，尽量完全消化再动手不迟

goabroad

我才刚开始学GUI呢

sunbird04

路过，瞧瞧