我用nmisc 9 求得电感，为什么求出来的值和实际值不符合？求斑竹斑竹！！

damao5277

下面是我的命令流，我求出来的初级线圈的电感值是0.508e-3,而实际值是1.584e-3，这是为什么呢，求斑竹解答，求高手解答！！！！
fini
/clear
/prep7  
/title,2-d elc-c static annalysis   
!/units,CGS
et,1,plane53,,,1 !air   
et,2,plane53,2,,1  !初级coil  
et,3,plane53,,,1  !iron
et,4,plane53,2,,1  !次级coil  
R,1,0.000018,110,23.1418,1,0.217,
R,2,3.5e-6,200,0.942,1,0.0268,
mp,murx,1,1       !air泊松比(磁导率)  
mp,murx,2,10758    !磁铁泊松比(磁导率)  
mp,murx,3,1        !线圈的泊松比(磁导率)   
mp,rsvx,3,1.74e-8    !  
mp,murx,4,1        !线圈的泊松比(磁导率)   
mp,rsvx,4,1.74e-8    !  
!   
*SET,n1,110 !   
*SET,n2,200
*SET,v,160 ! MAX
*SET,h,175000  !   
!   
*SET,r,0.0335   
*SET,s,0.003
*SET,d1,0.032   
*SET,d2,0.006   
*SET,l1,0.007   
*SET,l2,0.001   
*SET,l3,0.0005  
!   
rectng,d1,d1+s,0,d2
rectng,0,l3,1.5*r,1.5*r+l1  
rectng,l3,l2,1.5*r,1.5*r+l1
PCIRC,0,2*R,0,90,   
aovlap,all  
numcmp,area
!   
ASEL,s,,,4  
AATT,1,1,1,0,   
ASEL,s,,,1  
AATT,3,1,2,0,   
ASEL,,,,3   
AATT,4,2,4,0,   
ASEL,s,,,2  
AATT,2,1,3,0,   
allsel,all  
aplot   
MSHAPE,0,2D
MSHKEY,0
asel,s,,,1  
aesize,0.0005   
amesh,1
allsel,all  
asel,s,,,2,3
aesize,0.05
amesh,all   
MSHAPE,1,2D
MSHKEY,0
allsel,all  
asel,s,,,4  
aesize,0.005
amesh,all   
wpstyl,0.01,0.01,,,,,1,1
/view,all,wp
/dist,all,0.2   
eplot   
R,3,1,   
*GET,_zz2,active,,csys  
*GET,_zz3,active,,dsys  
CSYS,0  
DSYS,0  
N,847,-0.2E-01,0.6E-01,0
N,848,-0.2E-01,0.5E-01,0
CSYS,_zz2   
DSYS,_zz3   
RMOD,3,15,-0.1E-01,1
ET,5,CIRCU124,0,0   
TYPE,5  
REAL,3  
MAT,1   
E,847,848   
R,4,1,  
*GET,_zz2,active,,csys  
*GET,_zz3,active,,dsys  
CSYS,0  
DSYS,0  
CSYS,_zz2   
DSYS,_zz3   
RMOD,4,15,0.2E-01,1
ET,6,CIRCU124,5,0   
TYPE,6  
REAL,4  
MAT,1   
E,847,848,236   
FLST,2,1,1,ORDE,1   
FITEM,2,848
/GO
D,P51X,VOLT,0   
eplot   
  
/REPLOT,RESIZE  
R,5,160, ,  
*GET,_zz2,active,,csys  
*GET,_zz3,active,,dsys  
CSYS,0  
DSYS,0  
N,849,0.4E-01,-0.1E-01,0
N,850,0.3E-01,-0.1E-01,0
CSYS,_zz2   
DSYS,_zz3   
RMOD,5,15,0.1E-01,2
ET,7,CIRCU124,4,0   
TYPE,7  
REAL,5  
MAT,1   
N,851,0.35E-01,-0.2E-01,0   
E,849,850,851   
R,6,1,  
*GET,_zz2,active,,csys  
*GET,_zz3,active,,dsys  
CSYS,0  
DSYS,0  
CSYS,_zz2   
DSYS,_zz3   
RMOD,6,15,0,3   
ET,8,CIRCU124,5,0   
TYPE,8  
REAL,6  
MAT,1   
E,850,849,30
FLST,2,1,1,ORDE,1   
FITEM,2,850
/GO

D,P51X,VOLT,0   
esel,s,mat,,3! get coil elements
nsle,s  
cp,1,curr,all   
!cp,2,emf,all   
allsel,all  
esel,s,mat,,4! get coil elements
nsle,s  
cp,3,curr,all   
!cp,4,emf,all   
allsel,all  
!csys,1  
!nsel,s,loc,x,2*r
!d,all,az,0 ! set maxswel surface flag   
!csys,0  
!nsel,s,loc,x,0  
!nsel,a,,,20,21,1
!d,all,az,0! flux-parallel condition

lsel,s,,,18,21,1
lsel,a,,,13
lsel,a,,,1
lsel,a,,,8
!lsel,s,,,5,8,1  
dl,all,,asym
finish  
/solu   
antype,harm
esel,all
harfrq,175000 !175000Hz.
nsubst,1,
kbc,1,
allsel,all  
solve   
   
FINISH  
/POST1  
PRESOL,ELEM
/REPLOT,RESIZE  
/REPLOT,RESIZE  
/POST1  
SET,1,LAST,1,0, , ,
!*  
PRESOL,ELEM
PRESOL,ELEM
ETABLE, ,NMISC, 8   
!*  
ETABLE, ,NMISC, 9   
!*  
ESEL,S,MAT,,3   
SSUM