计算直流线圈电感（求大虾指点）

lqtading

各位大虾，你们好  在下初学ansys，今天参考资料试着编了一个计算直流线圈电感的程序，但结果不太对，而且线圈中的磁力线分布也不正确，哪位大虾有时间可以指点一下，不胜感激！！（施加边界条件那部分自己不是很懂）
FINISH
/CLEAR

/TITLE,3-D STATIC COID
/FILNAME,SELF_INDUCTANCE,1
KEYW,MAGNOD,1

/PREP7
EMUNIT,MKS  ! MKS单位制

N=2450      ! 线圈匝数
I=170       ! 线圈电流（A）   
HC=0.1        ! 线圈的长度（米）  
R=0.25      ! 线圈内径（米）
WC=0.25    ! 线圈厚度（米）


ET,1,36          !电流源单元
R,1,1,N*I,WC,HC  !1号实常数：线圈类型，电流，厚度，长度

N,1,R,0,0.5
N,2,0,R,0.5
N,3,-R,0,0.5


TYPE,1
REAL,1

E,1,2,3

/ESHAPE,1
/VIEW,1,1,1,1
/VUP,1,Z
/TRIAD,LBOT
/TYPE,1,HIDP
EPLOT

ESEL,S,TYPE,,1
CM,COIL1,ELEM


ET,2,96
MP,MURX,1,1

/PNUM,VOLU,1
CYL4,,,0,0,1,360,1
VSEL,S,VOLU,,1
VATT,1,,2

ALLSEL,ALL
SMRT,6
MSHAPE,1,3D
MSHKEY,0
VMESH,1
/PNUM,MAT,1
/NUMBER,1
EPLOT


/SOLU
ALLSEL,ALL

NSEL,EXT
D,all,MAG,0
ALLSEL,ALL
MAGSOLV,3,,,,,1
FINISH
SAVE

/SOLU
*DIM,CUR,ARRAY,1
CUR(1)=I
LMATRIX,1,'COIL','CUR'
FINISH
SAVE



/post1
PLVECT,B

lqtading

补充一下，我建模的思路是：先建立线圈模型 然后外加圆柱形的空气场并对空气场进行剖分，从而建立有限元模型