3d永磁磁场下的线圈磁链总是计算不出来？

tjuchris

我是利用lmatrix宏计算的单个线圈总磁链。为了计算永磁体与线圈形成的匝链，线圈上加极小电流，线圈为空心线圈。可是计算结果为
磁链 coil linkage psi= —0.56924E+05
电感 inductance L=0.638E+08
气隙中的磁场强度最大才为0.7 T，线圈匝数为150，而且线圈的几何参数全是毫米计算，得到的结果不可能有10^5这么大呀
我看过help，读过别的程序，改来改去却始终没有得到正确的结果
问题到底出在哪儿呢？
请大家帮看看程序，多谢了

程序如下：
FINI
/CLEAR
/PREP7
/TITLE,3d arc coil linkage calculation
/PNUM,VOLUME,1
EMUNIT,MKS

INNERRADIUS=0.058      
PMTHICKNESS=0.022
OUTRADIUS=0.3
         
T1=0
T2=0
FAI=60
AIRGAP=0.001            
A=0.006                 
B=0.012                 
HEIGHTCOIL=0.015      
HC=930000               
N=150                    
I=0.0001   !current per turn                  
ACOIL=(B-A)*HEIGHTCOIL   
JDENS=N*I/ACOIL         

*AFUN,DEG
CSYS,4
WPROTA,,90
CYL4,0,0,INNERRADIUS,-90,INNERRADIUS+PMTHICKNESS,90
K,100,0,(INNERRADIUS+PMTHICKNESS)*SIN(FAI/2)
K,101,0.1,(INNERRADIUS+PMTHICKNESS)*SIN(FAI/2)
L,100,101
K,102,0,-(INNERRADIUS+PMTHICKNESS)*SIN(FAI/2)
K,103,0.1,-(INNERRADIUS+PMTHICKNESS)*SIN(FAI/2)
L,102,103
LSEL,S,LINE,,5,6
ASBL,ALL,ALL,DELETE,DELETE
ADELE,2,3,1,1
KDELE,101,104,1,1
NUMCMP,ALL
CSYS,0
K,1000,0,0,-1
K,1001,0,0,1
VROTAT,1,,,,,,1000,1001,60,1
KDELE,1000,1001,1
NUMCMP,ALL
WPROTA,,-90
CSYS,1
VGEN,6,1,1,0,0,60
NUMCMP,ALL
CSYS,0

WPROTA,T1
WPROTA,,,T2
WPOFFS,INNERRADIUS+PMTHICKNESS+AIRGAP
WPROTA,,,90
CYL4,0,0,A,0,B,360,HEIGHTCOIL
WPROTA,0,0,-90
WPAVE,0,0,0
WPROTA,,,-T2
WPROTA,-T1

SPHERE,0,OUTRADIUS,0,360
VOVLAP,ALL
NUMCMP,VOLU

ET,1,SOLID97
KEYOPT,1,1,0

MP,MURX,1,1
MP,MURX,2,1
MP,MGXX,2,HC
MP,MURX,3,1
MP,MGXX,3,-HC
MP,MURX,4,1
LOCAL,11,0,,,,30
VSEL,S,VOLU,,1
VATT,2,0,1,11
LOCAL,12,0,,,,-30
VSEL,S,VOLU,,3
VATT,3,0,1,12
LOCAL,13,0,,,,90
VSEL,S,VOLU,,4
VATT,3,0,1,13
LOCAL,14,0,,,,150
VSEL,S,VOLU,,5
VATT,2,0,1,14
LOCAL,15,0,,,,-150
VSEL,S,VOLU,,6
VATT,3,0,1,15
LOCAL,16,0,,,,-90
VSEL,S,VOLU,,7
VATT,2,0,1,16
CSYS,0
VSEL,S,VOLU,,2
VATT,4,0,1
VSEL,S,VOLU,,8
VATT,1,0,1
allsel,all
/PNUM,MAT,1
vsel,s,mat,,4
mshape,1,3d
mshkey,0
smrtsize,1
vmesh,all
vsel,s,mat,,2,3,1
mshape,1,3d
mshkey,0
smrtsize,1
vmesh,all
vsel,s,mat,,1
mshape,1,3d
mshkey,0
smrtsize,2
vmesh,all

LOCAL,17,1,,,0,T1,0,90+t2
VSEL,S,VOLU,,2
ESLV
EMOD,ALL,ESYS,17      
BFE,ALL,JS,,,JDENS

allsel,all
csys,2
ASEL,S,LOC,X,OUTRADIUS
NSLA
NROTAT,ALL
D,ALL,AY,0
D,ALL,AZ,0
LSEL,S,LOC,X,OUTRADIUS
NSLL
NROTAT,ALL
D,ALL,AY,0
D,ALL,AZ,0
ALLSEL,ALL
FINISH
allsel
/SOLU
allsel
EQSLV,front
solve
SAVE
FINISH

/post1
*dim,cur,array,1   !定义一个N阶数组，N代表线圈个数，N=1
cur(1)=i  !给每个线圈定义一个名义电流，空载时，名义电流用一个极小的电流代替
esel,s,mat,,4
cm,coil1,elem         
allsel,all
lmatrix,1,'coil','cur','linkage'
finish

xuj

我也有过类似的问题，对于lmatrix这个命令我不了解它在幕后是怎么运算的，所以我也不知道为什么答案是错误的。
我用的方法是，先算它的vectorpotential，因为flux就是线圈内的vectorpotential的差，然后你就可以得到答案了。这种方法我在直流电机里面试过，相对于其他的方法来说还是比较准确的。