根据版中关于永磁体的帖子,钻研一个星期,终于求得结果,但是结构与实际值有出入,差一个数量级,且随时间变化不同。特发帖请教!
问题,如图所示(2D):在铝制导板上有一个旋转的永磁轮,n.s极取halbach列,计算在当轮旋转时,通过和铝板之间的磁场感应在轮上产生的x,y轴的力(旋转速度1000转/分,时间为0~0.1s,每子步0.002s)。命令流如下:/prep7
*afun,deg
ang=90
a=45
an=22.5
r_in=0.1725
r_out=0.3
es=1
hc=10e5
mur=1.08
rho=0.2857e-7
g=0.01
d=0.015
w=1.2
h=1
et,1,53!空气
mp,murx,1,1
et,2,53,,1
mp,murx,2,mur
mp,mgxx,2,hc
et,3,53,,1
mp,murx,3,mur
mp,mgyy,3,hc
et,4,53,,1
mp,murx,4,mur
mp,mgxx,4,-hc
et,5,53,,1
mp,murx,5,mur
mp,mgyy,5,-hc
et,6,53!导板
mp,murx,6,1
mp,rsvx,6,rho
!定义实常数
r,2,,,,,,,
rmore,,100/6,0,0,
r,3,,,,,,,
rmore,,100/6,0,0,
r,4,,,,,,,
rmore,,100/6,0,0,
r,5,,,,,,,
rmore,,100/6,0,0,
!建模
local,11,1
csys
/PNUM,AREA,1
CYL4,,,r_in,ang-an,r_out,ang
CYL4,,,r_in,a,r_out,ang-an
CYL4,,,r_in,an,r_out,a
CYL4,,,r_in,0,r_out,an
!将模型分别沿x,y轴映射
FLST,3,4,5,ORDE,2
FITEM,3,1
FITEM,3,-4
ARSYM,Y,P51X, , , ,0,0
FLST,3,8,5,ORDE,2
FITEM,3,1
FITEM,3,-8
ARSYM,X,P51X, , , ,0,0
aglu,all
rectng,-w/2,w/2,-0.3-g-d,-0.3-g !建立导板模型
rectng,-w/2,w/2,-h/2,h/2 !建立外边界模型
aovlap,all
numcmp,area !将生成的面重新编号
aplot
!赋予材料特性
ASEL,S, , , 3
FLST,5,3,5,ORDE,3
FITEM,5,3
FITEM,5,5
FITEM,5,-6
ASEL,S, , ,P51X
/REPLOT
AATT, 1, , 1, 0,
ALLSEL,ALL
CM,_Y,AREA
ASEL, , , , 2
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT, 6, , 1, 0,
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
csys,1
asel,s,loc,x,r_in,r_out
/REPLOT
AATT, 2, , 1, 11,
alls
aplot
!分网,先将导板单独分网0.005,其他0.02,且自由划分
ASEL,S, , , 2
/REPLOT
FLST,2,1,5,ORDE,1
FITEM,2,2
AESIZE,P51X,0.005,
MSHAPE,0,2D
MSHKEY,1
!*
CM,_Y,AREA
ASEL, , , , 2
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
APLOT
ASEL,U, , , 2
/REPLOT
FLST,2,19,5,ORDE,3
FITEM,2,1
FITEM,2,3
FITEM,2,-20
AESIZE,P51X,0.02,
MSHAPE,1,2D
!*
FLST,5,19,5,ORDE,3
FITEM,5,1
FITEM,5,3
FITEM,5,-20
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
MSHKEY,0
!*
FLST,5,19,5,ORDE,3
FITEM,5,1
FITEM,5,3
FITEM,5,-20
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
SAVE
!修改轮中的永磁极性
asel,s,mat,,2
asel,r,loc,y,a,ang-an
alls,belo,area
mat,3
type,3
real,3
esys,11
emod,all
alls
asel,s,mat,,2
asel,r,loc,y,an,a
alls,belo,area
mat,4
type,4
real,4
esys,11
emod,all
alls
asel,s,mat,,2
asel,r,loc,y,0,an
alls,belo,area
mat,5
type,5
real,5
esys,11
emod,all
alls
asel,s,mat,,2
FLST,5,3,5,ORDE,3
FITEM,5,4
FITEM,5,7
FITEM,5,18
ASEL,S, , ,P51X
/REPLOT
alls,belo,area
mat,5
type,5
real,5
esys,11
emod,all
alls
asel,s,mat,,2
/REPLOT
FLST,5,3,5,ORDE,3
FITEM,5,13
FITEM,5,15
FITEM,5,20
ASEL,S, , ,P51X
/REPLOT
alls,belo,area
mat,4
type,4
real,4
esys,11
emod,all
alls
asel,s,mat,,2
/REPLOT
FLST,5,3,5,ORDE,3
FITEM,5,14
FITEM,5,16
FITEM,5,19
ASEL,S, , ,P51X
/REPLOT
alls,belo,area
mat,3
type,3
real,3
esys,11
emod,all
ALLSEL,ALL
/REPLOT
!加载求解
!定义组件wheel
/SOL
ESEL,S,MAT,,2
ESEL,A,MAT,,3
ESEL,A,MAT,,4
ESEL,A,MAT,,5
/REPLOT
CM,wheel,ELEM
FMAGBC,'WHEEL'
ALLSEL,ALL
!施加边界通量平行条件
FLST,5,8,4,ORDE,7
FITEM,5,6
FITEM,5,9
FITEM,5,-10
FITEM,5,13
FITEM,5,15
FITEM,5,21
FITEM,5,-23
LSEL,S, , ,P51X
NSLL,S,1
FLST,2,452,1,ORDE,10
FITEM,2,1
FITEM,2,-2
FITEM,2,482
FITEM,2,-488
FITEM,2,968
FITEM,2,-972
FITEM,2,2797
FITEM,2,-2951
FITEM,2,6791
FITEM,2,-7073
D,P51X,AZ,0
ALLSEL,ALL
/SOL
ANTYPE,4
TRNOPT,FULL
NROPT,AUTO, ,
EQSLV,SPAR, ,0,
PRECISION,0
MSAVE,0
TIME,0.1
AUTOTS,-1
DELTIM,0.02, , ,1
KBC,0
FINISH
save
最后得到力值为:
SUMMARY OF FORCES BY VIRTUAL WORK
Load Step Number: 1.
Substep Number: 1.
Time: 0.1000E+01
Units of Force: (N/m)
Component Force-X Force-Y
WHEEL -0.54290E+02 -0.68597E+03
_________________________________________
SUMMARY OF FORCES BY MAXWELL STRESS TENSOR
Units of Force: (N/m)
Component Force-X Force-Y
WHEEL -0.10416E+04 -0.25962E+03
________________________________________
且x,y轴上的力随时间的变化曲线如图2(线性变化图)所示,但是实际的曲线图如图3(随时间逐渐达到稳态,且为正弦图)所示,两图之间的差别是为什么啊?请教高手!
[ 本帖最后由 NX220 于 2008-12-8 13:39 编辑 ] |