关于永磁铁之间吸力的三维静磁场分析，帮帮忙，怎么算

lgd5401

小弟想计算空间两块磁铁之间的吸力，永磁铁尺寸21*21*7mm，间隙4mm，下面是我写的程序，计算结果是零，请大家看一下，指导指导那里错了
!---------------
finish
/clear,start
/filename,mag3d
/title,mg3d
/prep7
emunit,mks     
!define constant
L=0.021  ! length of PM
H=0.007  ! height of PM
D=0.002  ! air gap
!define ele type
et,1,solid98,10
mp,murx,1,1
mp,murx,2,1.3412
mp,murx,3,1.3412
mp,mgyy,2,890000
mp,mgyy,3,-890000
!define model
blc4,0,0,L,H,L
blc4,0,H+D,L,H,L
blc4,0,H,L,D,L
allsel,all

!apply feature
vsel,s,volu,,1
vatt,2,,1
vsel,s,volu,,2
vatt,3,,1
vsel,s,volu,,3
vatt,1,,1
allsel,all
! mesh model
vmesh,all
esel,s,mat,,2
cm,mgdown,elem
esel,s,mat,,3
cm,mgup,elem
!apply loads
fmagbc,'mgdown'
fmagbc,'mgup'
allsel,all
finish
!solution
/solu
magsolv
finish
!postproc
/post1
fmagsum,'mgdown','mgup'

这上面的程序是不对的，整了一晚上，终于得到个靠谱的，下面是我新写的程序，大大们看看有错误没

!---------------
finish
/clear,start
/filename,mag3d
/title,mg3d
/prep7
emunit,mks     
!define ele type
et,1,solid96
keyopt,1,5,0
mp,murx,1,1
mp,murx,2,1.3412
mp,murx,3,1.3412
mp,mgzz,2,-890000
mp,mgzz,3,-890000
!define model
block,10.5,-10.5,-10.5,10.5,2,9
block,10.5,-10.5,-10.5,10.5,-2,-9
block,11.5,-11.5,-11.5,11.5,1,10
block,11.5,-11.5,-11.5,11.5,-1,-10
vovlap,all
numcmp,volu
cylind,50,0,-50,50,0,360,
vovlap,all
!apply feature
vsel,s,volu,,1
vatt,2
vsel,s,volu,,2
vatt,3
vsel,s,volu,,3,4
vatt,1
! mesh model
allsel,all
smrt,7
mshape,1,3d
mshkey,0
vmesh,all
!apply loads
esel,s,mat,,2
cm,mgdown,elem
esel,s,mat,,3
cm,mgup,elem
fmagbc,'mgdown'
fmagbc,'mgup'

allsel,all
vlscale,all,,,.001,.001,.001,,0,1
!apply boundary conditions
d,1466,mag,0   

finish
!solution
/solu
allsel,all
magsolv,2
finish
!postproc
/post1
fmagsum,'mgdown','mgup'

ahua326

想问一个问题，边界条件是怎么加载的呢 怎么就在一个节点上加载呢  

lgd5401

我看一个例子上做的，应该是远点边界吧

majianzzu

我也遇到了这个问题，希望高人指点一下。

frankygtd

三维的应该不是点边界把