可能有的人已经在某些论坛上见到过这个帖子了,我按照这做了好多遍,还是没做出来。谢谢大家了
代码如下,请高手指点!!!
/PREP7
emunit,mks
sd1=0.025 !内环钢性支架内直径
sd2=0.043 !内环钢性支架外直径
id1=0.029 !内磁环1内直径
id2=0.043 !内磁环1外直径
height1=0.008 !内环磁体高度
height2=0.020 !内环钢性支架高度
sd3=0.044 !外环钢性支架内直径
sd4=0.062 !外环钢性支架外直径
od1=0.044 !外磁环内直径
od2=0.058 !外磁环外直径
height3=0.008 !外环磁体高度
height4=0.020 !外环钢性支架高度
height0=0.002 !壁厚
hc=960000 !永磁体矫顽力
gap=0.0005 !单边气隙
n=2
c_x=0
c_y=0
c_z=0.001
!单元,材料选择
ET,1,SOLID96
MP,MURX,1,1 !空气磁导率
MP,MGXX,2,,hc
MP,MGyy,2,,0
MP,MGzz,2,,0
mp,murx,2,1.05
MP,MGXX,3,,-hc
MP,MGyy,3,,0
MP,MGzz,3,,0
mp,murx,3,1.05
!建立模型
!!!!!!
CYLIND,0.0215,0.0145,-0.004,0.004,0,360,
CYLIND,0.029,0.022,0,0.008,0,360,
!指定属性
wpcsys,-1,0
local,11,1,c_x,c_z,c_y,0,-90,0,1,1
local,12,1,0,0,0,0,-90,0,1,1
local,13,1,c_x,c_z+height1,c_y,0,-90,0,1,1
local,14,1,0,height1,0,0,-90,0,1,1
vsel,s,volu,,1 !动环1
vatt,2
vsel,s,volu,,2 !定环1
vatt,3
allsel,all
/pnum,mat,1
vplot
!划分网格
MSHKEY 0 !自由网格
MSHAPE,1,3D !2D四边形单元
SMRTSIZE,6
vMESH,ALL
!!!!修改单元坐标系
vsel,s,volu,,1 !动环1
eslv,s
emodif,all,esys,11
vsel,s,volu,,2 !定环1
eslv,s
emodif,all,esys,12
!!!!!!!!!!!!!!
vsel,s,volu,,1,4,3
vsel,a,volu,,3,6,3
eslv,s
cm,rotor,elem
fmagbc,'rotor'
allsel,all
FINISH
/SOLU
D,2,MAG,0
ALLSEL,ALL
MAGSOLV,2 !选用RSP法求解
FINISH
/POST1
rsys,0
SET,LAST
FMAGSUM,'rotor'
这是有人做出的最终结果
[ 本帖最后由 80790800 于 2008-7-12 08:16 编辑 ] |