请问三维圆柱体的感应加热模拟结果为什么不对称

ansys平安

这段时间一直在做三维的圆柱体的感应加热的模拟，使用帮助文件里的二维感应加热的例子的循环语句，但一直无法得到一个堆成的结果，以下是我的命令流，就是一个圆柱钢轴，外套一个感应线圈，其他部分用空气包裹。线圈，钢轴和线圈与钢轴间的空气都是使用的扫描网格划分，其他部分使用自由划分，一下是命令流，求高手指正，先谢谢了。
/TITLE, gangguan INDUC         !
/PREP7  
EMUNIT,MKS                    !
*AFUN,DEG                     !
ET,1,SOLID117                 !
ET,2,SOLID117,1               !
ET,3,SOLID117                 !
MP,MURX,1,1                   !
MP,RSVX,1,0.1724E-8           !
MP,MURX,2,1.01                !
MP,RSVX,2,0.88E-6             !
MP,MURX,3,1                   !
*SET,FTIME,8                  !!!   
*SET,TINC,0.5                                  !
*SET,TIME,0                                  !  
*SET,DELT,0.1                                !   
CYLIND,0.015, ,0.02,0.06,0,360  
CYLIND,0.02,0.0175,0.05,0.06,0,360  
CYLIND,0.015,0.0175,0.05,0.06,0,360
Vglue,all   
NUMCMP,VOLU      
lsel,s, , ,5,8,1                   ！根据编号选择线段
lsel,a, , ,23,26,1  
lsel,a, , ,19,22,1  
lesize,all, , ,30   
lsel,s, , ,51   
lesize,all, , ,30   
lsel,s, , ,53   
lesize,all, , ,10   
CYLIND,0.05,0,0,0.07,0,360  
vovlap,all  
NUMCMP,VOLU
VSEL,S,VOLUM, ,1              !将线圈标为1号材料，3
VATT,1,1,3  
VSEL,S,VOLUM, ,3,4,1             !将空气标为3号材料，1
VATT,3,1,1  
VSEL,S,VOLUM, ,2               !将钢标为2号材料，2
VATT,2,1,2  
VSEL,ALL
vsweep,2
vsweep,3
vsweep,1
ESIZE,0.008,0,   
VSEL, , , ,       4
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
MSHKEY,0
MSHAPE,1,3d
FLST,2,1,6,ORDE,1   
FITEM,2,4   
VIMP,P51X,1
FINISH  
/SOL
CSYS,1                                    !
ASEL,S,LOC,X,0.05                         !
ASEL,A,LOC,Z,0.07   
ASEL,A,LOC,Z,0  
DA,ALL,AZ,0                               !
ASEL,ALL
CSWPLA,11,1,1,1,                       !给单元创建柱坐标系
ESEL,S,MAT, ,1  
TYPE,   1   
MAT,       1
REAL,   
ESYS,      11   
SECNUM,
! TSHAP,LINE
!*  
BFE,ALL,JS,1, ,0.84E8   
ESEL,ALL
ANTYP,HARM                                !
HARFR,80000                               !
PHYSICS,WRITE,EMAG                        !     写入emag  
FINISH  
!   
/PREP7  
LSCLEAR,ALL                               !
ET,1,0                                    !     将1号材料转换成0
ET,2,SOLID90                              !     将2号材料转换为SOLID90   
ET,3,0                                    !     将3号单元转换成0
MP,KXX,2,30                                !!!!!传导率各向同性
MP,C,2,500                                !
MP,DENS,2,7900                            !
FINISH  
/SOLU   
ANTYPE,TRANS                              !
TOFFST,273  
TUNIF,25                                  !
CNVTOL,HEAT,1                             !
KBC,1                                     !
TRNOPT,FULL
AUTOS,ON                                  !     自动时间-   
DELTIM,5E-3,5E-4,DELT,ON                  !
OUTRES,BASIC,ALL                           !
PHYSICS,WRITE,THERMAL                      !     写入THERMAL  
SAVE
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,25                             ! 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   
parsav,scalar,parameter,sav !save parameters before multiframe rest
*IF,I,GT,1,THEN
ANTYPE,TRANS,REST                    ! thermal restart  
*ENDIF  
parres,new,parameter,sav !restore parameters after multiframe re
TIME,TIME                             ! time at end of thermal run  
ESEL,S,MAT,,2                         ! select billet region
LDREAD,HGEN,,,,2,,RST                 ! 从emag
ESEL,ALL
SOLVE   
FINISH  
/ASSIGN,ESAV                          ! reassign files to default   
/ASSIGN,EMAT
*ENDDO  




关于单元坐标轴，改不改的结果都是一样。

ansys平安

各位大虾啊，出个主意吧，实在是困在这里好久了。。拜谢。

ansys平安

自己顶一个。。。。。

zjuily

这模拟是很吻合实际阿，看起来像边界效应，不知你的分析是想达到什么情形？

phyee88

看着这图还不错啊