求助ansys电磁场分析大神

xuanqi_ldr

用ansys计算箔式绕组变压器的绕组涡流损耗，怎么着也出不来集肤效应和邻近效应，跪求指点，命令流如下
finish
/clear

/prep7  
/title, two-winding transformer
et,1,53
et,2,53,1   
et,3,53,1   
mp,murx,1,1
mp,murx,2,10000
mp,murx,3,1
mp,rsvx,3,2.8e-8
mp,murx,4,1
mp,rsvx,4,2.8e-8
f,5000
a,35   
b,100  
c,200  
d,5
e,2
g,88.5
s,50   
h,176  
rectng,0,b+a+a,0,a+a+c                                ! iron
rectng,a,b+a,a,a+c  
rectng,s,s+d,a+(c-h)/2,a+(c-h)/2+h                         ! primary
rectng,s+d+e,s+d+e+d,a+(c-h)/2,a+(c-h)/2+h  
rectng,s+d+e+d+e,s+d+e+d+e+d,a+(c-h)/2,a+(c-h)/2+h  
rectng,g,g+d,a+(c-h)/2,a+(c-h)/2+h                         ! secondary  
rectng,g+e+d,g+e+d+d,a+(c-h)/2,a+(c-h)/2+h  
rectng,g+d+e+d+e,g+d+e+d+e+d,a+(c-h)/2,a+(c-h)/2+h  
/pnum,area,1
aovlap,all  
numcmp,area
aplot   

asel,s,area,,1,3
aatt,3,1,2  
asel,s,area,,4,6
aatt,4,1,3
asel,s,area,,7  
aatt,2,1,1  
asel,s,area,,8  
aatt,1,1,1  
/pnum,mat,1
allsel,all  
aplot   
FLST,5,6,4,ORDE,6   
FITEM,5,11  
FITEM,5,15  
FITEM,5,19  
FITEM,5,23  
FITEM,5,27  
FITEM,5,31  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,10, , , , ,1
!*  
FLST,5,6,4,ORDE,6   
FITEM,5,12  
FITEM,5,16  
FITEM,5,20  
FITEM,5,24  
FITEM,5,28  
FITEM,5,32  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,40, , , , ,1
!*  
MSHAPE,0,2D
MSHKEY,1
!*  
FLST,5,6,5,ORDE,2   
FITEM,5,1   
FITEM,5,-6  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
asel,s,mat,,3,4
ames,all
FINISH  
/SOL
FINISH  
/PREP7  
SMRT,6  
SMRT,3  
/UI,MESH,OFF
ALLSEL,ALL  
APLOT   
MSHKEY,0
!*  
FLST,5,2,5,ORDE,2   
FITEM,5,7   
FITEM,5,-8  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  

esel,s,mat,,3     !线圈组件
cm,coil1,elem   
esel,s,mat,,4   
cm,coil2,elem   
esel,,ename,,53   !磁力线平行边界条件  
nsle

ASEL,S, , ,       1
NSLA,S,1
FLST,4,1301,1,ORDE,2
FITEM,4,1   
FITEM,4,-1301   
CP,1,VOLT,P51X  
ASEL,S, , ,       2
NSLA,S,1
FLST,4,1301,1,ORDE,2
FITEM,4,1302
FITEM,4,-2602   
CP,2,VOLT,P51X  
ASEL,S, , ,       3
ASEL,S, , ,       3
NSLA,S,1
FLST,4,1301,1,ORDE,2
FITEM,4,2603
FITEM,4,-3903   
CP,3,VOLT,P51X  
ASEL,S, , ,       4
NSLA,S,1
FLST,4,1301,1,ORDE,2
FITEM,4,3904
FITEM,4,-5204   
CP,4,VOLT,P51X  
ASEL,S, , ,       5
NSLA,S,1
NPLOT   
FLST,4,1301,1,ORDE,2
FITEM,4,5205
FITEM,4,-6505   
CP,5,VOLT,P51X  
ASEL,S, , ,       6
NSLA,S,1
FLST,4,1301,1,ORDE,2
FITEM,4,6506
FITEM,4,-7806   
CP,6,VOLT,P51X  
ALLSEL,ALL  
EPLOT   

nsel,r,ext  
d,all,az
alls

FINISH  
/solu   
antype,harmic   
harfrq,0,5000   
NSUBST, ,   
KBC,1   
*SET,nn1amps,node(s+d/2,a+c/2,0)
f,nn1amps,amps,200,0
*SET,nn2amps,node(s+d+e+d/2,a+c/2,0)
f,nn2amps,amps,200,0
*SET,nn3amps,node(s+d+e+d+e+d/2,a+c/2,0)
f,nn3amps,amps,200,0
*SET,nn4amps,node(g+d/2,a+c/2,0)
f,nn4amps,amps,-200,0   
*SET,nn5amps,node(g+d+e+d/2,a+c/2,0)
f,nn5amps,amps,-200,0   
*SET,nn6amps,node(g+d+e+d+e+d/2,a+c/2,0)
f,nn6amps,amps,-200,0   
alls
solve   
fini