计算导线串入磁环的电感

青梅煮酒论英雄

现在做在导线上串上高导磁环，计算其产生的电感。请问边界条件如何施加？用Imatrix计算电感时，导线可以看做是单匝的线圈吗？以下是编写的命令流，最终磁感应强度B应该是环绕导线的，但计算结果明显不对，求大神看看哪地方不对？所建模型中最中间为导线，外面的为磁环，其余为空气。
/prep7
n1=1         
s1=0.0000015                 
symfac=1   
nc=1        ! number of coils                              
*dim,cur,array,nc              
*dim,coils,char,nc            
coils(1)='wind1'  

et,1,117  
mp,murx,1,1300
mp,murx,2,1
mp,murx,3,1
mp,rsvx,3,0.0000000175

CYLIND,0.00138, ,0,0.028,0,360,
CYLIND,0.0174,0.009,0,0.028,0,360,  
CYLIND,0.03,0,0,0.028,0,360,
vovlap,all

VSEL, , , ,       1
VATT,       3, ,   1,       0  
VSEL, , , ,       2
VATT,       1, ,   1,       0  
VSEL, , , ,       4
VATT,       2, ,   1,       0  
VSEL, , , ,       5
VATT,       2, ,   1,       0           
allsel,all

SMRTSIZE,4
MSHAPE,1                          
mshkey,0                        
vmesh,all
allsel

cur(1)=0.001
JS1=cur(1)*n1/s1
VSEL,S,MAT,,3,3
bfv,all,js,JS1,0,0
cm,coils(1),elem         

/prep7
/solu
antype,0                          
eqslv,front                    
magsolv                       
finish


/post1
/out
/com,   
*abbr,*uilist,/replot            
lmatrix,symfac,'wind','cur','ind'
L1=ind(1,1)