lgd5401 发表于 2011-8-25 18:29:00

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

本帖最后由 lgd5401 于 2011-8-25 21:30 编辑

小弟想计算空间两块磁铁之间的吸力,永磁铁尺寸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 发表于 2011-8-29 09:27:47

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

lgd5401 发表于 2011-8-30 20:49:56

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

majianzzu 发表于 2011-9-2 15:46:36

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

frankygtd 发表于 2012-3-17 15:21:18

三维的应该不是点边界把
页: [1]
查看完整版本: 关于永磁铁之间吸力的三维静磁场分析,帮帮忙,怎么算