请教大家ansys二维电磁感应加热程序，为什么没有进行感应加热？

AlBertt

/batch,list
/filenam, 2D
/title, 2D
/prep7
/com,
shpp,off


et,1,13                             ! PLANE13, AZ dof
et,2,13

ftime=3                            ! final time
tinc=1                                ! time increment for harmonic analysis
time=0                          ! initialize time
delt=0.01                       ! maximum delta time step
       
emunit,mks                           ! set magnetic units
mp,murx,1,1                     ! air relative permeability
mp,murx,3,1                     ! coil relative permeability
mp,dens,3,8960  
mp,RSVX,3,2e-5                         ! resistivity
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
mptemp
mptemp,1,20,100,200,300,400,500                     ! 定义与磁导率对应的温度
mptemp,7,600,700,720,742,760,800
mptemp,13,900,1000
mpdata,kxx,2,1,46.71,46.06,45.59,43.47,40.7,37.67
mpdata,kxx,2,7,34.63,31.66,30.57,29.05,27.57,23.92
mpdata,kxx,2,13,26.32,27.76                        !定义钢坯料热传导系数
mp,dens,2,7850  
mp,c,2,465
mp,emis,2,.68                     ! emissivity


/AUX15                          !A2
!*  
IOPTN,IGES,SMOOTH   
IOPTN,MERGE,YES
IOPTN,SOLID,YES
IOPTN,SMALL,YES
IOPTN,GTOLER, DEFA  
IGESIN,'special-shaped coils01','IGS','Desktop\'
APLOT   
!*  
/VIEW,1,,,1                         !front
/ANG,1  
/REP,FAST
FINISH


/PREP7  
RECTNG,-20,20,-20,20,                !A4
RECTNG,-18,18,-18,18,     
ASBA,       1,       3  

RECTNG,-50,50,-50,50,           !A3,A5,A6
FLST,3,2,5,ORDE,2   
FITEM,3,2   
FITEM,3,4   
ASBA,       1,P51X  

RECTNG,-20,20,-20,20,                !A4
RECTNG,-18,18,-18,18,     
ASBA,       1,       2  

FINISH


/AUX15                          !A2
!*  
IOPTN,IGES,SMOOTH   
IOPTN,MERGE,YES
IOPTN,SOLID,YES
IOPTN,SMALL,YES
IOPTN,GTOLER, DEFA  
IGESIN,'special-shaped coils01','IGS','Desktop\'
APLOT   
!*  
/VIEW,1,,,1                         !front
/ANG,1  
/REP,FAST
FINISH


/PREP7  
FLST,5,3,5,ORDE,3                   !define mesh attributes
FITEM,5,3   
FITEM,5,5   
FITEM,5,-6  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1, ,   2,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       4
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       3, ,   2,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  

     
SMRT,1  
MSHAPE,0,2D
MSHKEY,0
!*  
FLST,5,5,5,ORDE,2   
FITEM,5,2   
FITEM,5,-6  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  


NSEL,S,LOC,X,-50
NSEL,A,LOC,X,50
NSEL,A,LOC,Y,-50
NSEL,A,LOC,Y,50
D,all, ,0, , , ,AZ

ESEL,S,MAT,,3   
FINISH  
/SOL
BFE,all,JS,1, , ,114270, ,

ALLSEL,ALL                         !be necessary


/solu   
antyp,harm  
harfrq,9500
HROPT,FULL  
HROUT,ON
EQSLV, ,1e-008,
physics,write,emag      
solve
SAVE  
finish

/prep7  
lsclear,all                     ! clear all b.c.'s and options  
et,1,55                             ! PLANE55 thermal element
et,2,0  
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,0.01,0.005,delt,on          ! time step control
outres,basic,all                ! save all load step information
physics,write,thermal               ! write thermal physics file
finish  

*SET,time,0
*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,,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

AlBertt

我这个程序时从圆钢坯的感应加热分析中修改过来的，就是把原来的一维模型改成现在的二维模型，电磁热耦合部分基本上没有修改，但是进行运算的时候，程序好像没有对坯料感应加热，不知道哪里出现问题，请教大家，谢谢！

AlBertt

知道什么原因了，是因为我施加载荷是垂直xy平面上的，而不是沿着环形线圈施加电流载荷，现在正找办法如何施加环形电流

madlea

AlBertt 发表于 2014-12-12 14:45
知道什么原因了，是因为我施加载荷是垂直xy平面上的，而不是沿着环形线圈施加电流载荷，现在正找办法如何施 ...

楼主在吗？有问题交流qq272511937，谢谢！