永磁体分析求教！计算的磁场力是实测的十分之一左右，求大侠们帮指点一下啊！谢谢！

niuhelen

/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Magnetic-Nodal   
!*  
*SET,I,6
*SET,N,500  
/PREP7  
!*  
EMUNIT,MKS  !国际单位
*SET,Br,1.18     !定义永磁轴承各个物理量, 这是个剩磁
*SET,Hc,955000   !这个是磁导  率

*SET,mur,4*acos(-1)*1E-7   !这是个弧度值
*SET,Hc,955000   !这个是磁导率
ET,1,SOLID96
!*
MP,MURX,1,1
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,MGXX,2,,
MPDATA,MGYY,2,,
MPDATA,MGZZ,2,,hc   
MPDATA,MURZ,2,,Br/Hc/mur

TBDE,BH,2,,,,  !这个是永磁体的bh曲线
TB,BH,2,1,25,
TBTEMP,O
!TBPT,,-955000+Hc,0
TBPT,,-926800+Hc,0.47
TBPT,,-905400+Hc,0.865
TBPT,,-875600+Hc,1.025
TBPT,,-835800+Hc,1.07
TBPT,,-796000+Hc,1.11
TBPT,,-756200+Hc,1.13
TBPT,,-716400+Hc,1.135
TBPT,,-676600+Hc,1.14
TBPT,,-636800+Hc,1.1425
!描述永磁材料的B-H曲线

TBPT,,-597000+Hc,1.145
TBPT,,-557200+Hc,1.1475
TBPT,,-517400+Hc,1.15
TBPT,,-477600+Hc,1.1505
TBPT,,-437800+Hc,1.1515
TBPT,,-398000+Hc,1.152
TBPT,,-358200+Hc,1.1545
TBPT,,-318400+Hc,1.157
TBPT,,-278600+Hc,1.1595
TBPT,,-238800+Hc,1.162
TBPT,,-199000+Hc,1.165
TBPT,,-159200+Hc,1.168
TBPT,,-119400+Hc,1.171
TBPT,,-79600+Hc,1.174
TBPTT,,-39800+Hc,1.177
TBPT,,10+Hc,1.18
TBPLOT,BH,2,,,,
CYLIND,0.005/2,0.022/2,0.000,0.003,0,90,
CYLIND,0.005/2,0.022/2,0.008,0.011,0,90,
CYLIND,0.030,0,-0.020,0.031,0,90,
VOVLAP,ALL
NUMCMP,ALL
VSEL,S, , ,1
VATT,2, ,   1,       0
VSEL,S, , ,2
VATT,2, ,   1,       0
VSEL,S, , ,3
VATT,1, ,   1,       0
ALLSEL
ESIZE,0.0025,0,  
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL   
VSEL,S, , ,       1
ESLV,S  
CM,arm,ELEM
FMAGBC,'ARM'
allsel
D,2,MAG,0,
D,10,MAG,0,
ALLSEL,ALL  
FINISH  
/SOL

MAGSOLV,3, , ,0.001,25,1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!