求助：计算出的永磁直线电机推力很小（e-14）

丁卯年代

最近在学习ansys做直线电机的仿真，在论坛上看到前辈“grandjiang”以前发表的一篇帖子：
3D永磁直线电机推力计算（APDL）

http://forum.simwe.com/thread-822180-1-1.html
，
于是照着做了一次，并将原帖中提到的问题做了修改，即“加载电流密度的流入和流出端面加边界平行条件”，但是计算出的Z方向的力依然很小，现在比较疑惑的是：到底是边界加载出了问题，抑或是力的计算方法出了问题，还是其他地方有问题，刚接触ansys以及电机分析，希望各位高手能够多加指点，谢谢！



图1、整体结 构



图2、计算得出的力，小得离谱（component是线圈与外部矽钢，即图1中的绿色与红色部分）

程序源文件：
!/BACH,LIST
/TITLE,LINEAR MOTOR
/FILNAME,L0514,1
KEYW,MAGEDG,1
!
/prep7
*afun,deg
*SET,t,12e-3
!结构参数
*SET,h,4e-3
*SET,m,20e-3

*SET,k,43e-3

*SET,s,36e-3
*SET,k1,7/12
*SET,k2,28/43
*SET,K3,4/7
*SET,K4,1/4
*SET,K5,30/36
*SET,ts,k1*t
*SET,d,k2*k
*SET,n,k3*ts
*SET,delta,k4*h
*SET,tau,k5*s

*SET,ml,6*t+t*(1-k1)

*SET,dl ,6*s+s*(1-k5)
*SET,z,300
*SET,I,3
*SET,tt,0
*SET,w,50
*SET,JSA,z*I*sin(w*tt)/(d-0.001)/ts
!动子三相电流密度
*SET,JSB,z*I*sin(w*tt+120)/(d-0.001)/ts
*SET,JSC,z*I*sin(w*tt+240)/(d-0.001)/ts
et,1,plane53
!定义单元
et,2,solid117
emunits,mks
mp,murx,1,1
!空气

mp,murx,2,1

!定子铁轭
TB,BH,3,1,30,
!矽钢片B-H曲线
TBPT,,355,0.7

TBPT,,405,0.8

TBPT,,470,0.9

TBPT,,555,1.0

TBPT,,673,1.1

TBPT,,836,1.2

TBPT,,1065,1.3

TBPT,,1220,1.35
TBPT,,1420,1.4

TBPT,,1720,1.45
TBPT,,2130,1.5

TBPT,,2670,1.55
TBPT,,3480,1.6

TBPT,,4500,1.65
TBPT,,5950,1.7

TBPT,,7650,1.75
TBPT,,10100,1.8
TBPT,,13000,1.85
TBPT,,15900,1.9
TBPT,,21100,1.95
TBPT,,26300,2.0
TBPT,,32900,2.05
TBPT,,42700,2.1
TBPT,,61700,2.15
TBPT,,84300,2.2
TBPT,,110000,2.25

TBPT,,135000,2.3
TBPT,,200000,2.41

TBPT,,400000,2.69

TBPT,,800000,3.22
mp,murx,4,1
!永磁体N极
mp,mgxx,4,895000

mp,murx,5,1
!永磁体S极
mp,mgxx,5,-895000

mp,murx,6,1
mp,rsvx,6,1.0e-8
！参数建模
blc4,0,0,ml,k
!动子模型
blc4,t*(1-k1),1e-3,ts,d
agen,6,2,,,t

blc4,t-ts+0.5*(ts-n),0,n,1e-3

agen,6,8,,,t
allsel,all

aovlap,all
numcmp,all
numcrg,all,,,,low
adele,1,6,1,1
blc4,-(2*S-tau/2)+(2*t+(t-ts)/2),-(m+h+delta),dl,m
!定子模型
blc4,-(2*S-tau/2)+(2*t+(t-ts)/2)+(s-tau),-(h+delta),tau,h
agen,6,2,,,s

blc4,-57e-3,-25e-3,270e-3,118e-3
allsel,all

aovlap,all

numcrg,all,,,,low

numcmp,all

/PNUM,AREA,1
VROTAT,all, , , , , ,81,82,90, ,
allsel,all
wpoff,0,-0.025,0
wpro,,,-90
cswpla,11,1,1,1
csys,11
!赋值材料属性
vsel,s,volum,,12
vatt,3,1,2,11
vsel,s,volum,,6,11,1

vatt,6,1,2,11
vsel,s,volum,,15
vatt,1,1,2,11
vsel,s,volum,,1,5,2
vatt,4,1,2,11
vsel,s,volum,,2,4,2
vsel,a,volum,,13
vatt,5,1,2,11
vsel,s,volum,,14
vatt,2,1,2,11
allsel,all
!网格划分
!
vsel,s,volum,,6,11,1
esize,0.004
vsweep,all
vsel,s,volum,,12
esize,0.004
vsweep,all
vsel,s,volum,,1,5,1
vsel,a,volum,,13,14,1
esize,0.004
vsweep,all
vsel,s,volum,,15
mshape,1,3d
mshkey,0
esize,0.015
vmesh,all
vsel,s,mat,,3,6,3
eslv,s
cm,'mover',elem
FMAGBC,'mover'
allsel,all
finish
/solu
antype,static
emodify,all,esys,11
!更改单元坐标系
vsel,s,volum,,6

!加载电流密度
eslv,s
bfe,all,js,1,,JSA
vsel,s,volum,,7

eslv,s
bfe,all,js,1,,-JSC
vsel,s,volum,,8

eslv,s
bfe,all,js,1,,JSB
vsel,s,volum,,9

eslv,s
bfe,all,js,1,,-JSA
vsel,s,volum,,10

eslv,s
bfe,all,js,1,,JSC
vsel,s,volum,,11

eslv,s
bfe,all,js,1,,-JSB
!*** Flux parallel boundary condition at all extern area of air ambient ***
allsel,all                                !修改后的边界条件
asel,s,ext
da,all,az,0
allsel,all
! Solution parameters
solve
finish
/post1
vsel,s,mat,,3,6,3
eslv,s
nsle,s,all
emft

丁卯年代

自己顶一个！

丁卯年代

任然没有解决，继续期待高手的指点！