轴向电磁轴承二维电磁场分析 出了问题

vwxy285

file:///C:/Documents
1为轴承定子，2为线圈，3为转轴，4为推力盘。
主要是建立二分之一二维模型，上图上端
结果：磁力线、电流密度等都有，就是没有力输出
我是初学者，请大家帮忙看看，给给高见
命令流如下：
!/clear
/title, axial magnetic bearing(2-D)

/filnam, Maglev_force
keyw,magnod,1

/prep7

et,1,plane53
!FE Defination
keyopt,1,3,1

mp,murx,1,1
!Define air's permeability
mp,murx,2,1
mp,murx,3,661
!rotor
45# steel
tb,bh,4
!BH of thrust disc,axial
tbpt,,251.2,0.77
tbpt,,502.4,1.14
tbpt,,753.6,1.33
tbpt,,1004.8,1.43
tbpt,,1256,1.47
tbpt,,1884,1.57
tbpt,,2512,1.61
tbpt,,3140,1.64
tbpt,,3768,1.66
tbpt,,4369,1.68
tbpt,,6280,1.72
tbpt,,10048,1.79
tbpt,,12560,1.83
tbpt,,15072,1.87
tbpt,,17584,1.9
tbpt,,20096,1.93
tbpt,,23108,1.96
tbpt,,25120,1.98
tbpt,,27632,2
tbpt,,30144,2.02
tbpt,,32656,2.04
tbpt,,35168,2.06
tbpt,,37680,2.07
tbplot,bh,4
!!!!!!!!!!!!!!!!!!!定义常量!!!
!!线圈
wending_length=15
wending_height=6
area_coil=wending_length*wending_height
turns=150
I=2
!I偏置电流
iz=1 !!ix控制电流,可以根据工况改变

!定子
core_r_1=25
core_r_2=30
core_r_3=37
core_r_4=40
core_g=5
!推力盘
thrust_disc_r=40
thrust_disc_h=9
!转子
rotor_r=20
rotor_length=30

!Modeling
/PNUM,AREA,1
! 打开面区域编号

RECTNG,0,rotor_r,0,2*rotor_length+thrust_disc_h,!转子
RECTNG,rotor_r,thrust_disc_r,rotor_length,rotor_length+thrust_disc_h,!推力盘
RECTNG,core_r_2+0.5,core_r_3-0.5,rotor_length+thrust_disc_h+0.6,rotor_length+thrust_disc_h+0.6+wending_length, !上线圈
RECTNG,core_r_2+0.5,core_r_3-0.5,rotor_length-0.6,rotor_length-0.6-wending_length, !下线圈
!上定子
RECTNG,core_r_1,core_r_2,rotor_length+thrust_disc_h+0.6,rotor_length+thrust_disc_h+0.6+wending_length+0.5,
RECTNG,core_r_3,core_r_4,rotor_length+thrust_disc_h+0.6,rotor_length+thrust_disc_h+0.6+wending_length+0.5,
RECTNG,core_r_1,core_r_4,rotor_length+thrust_disc_h+0.6+wending_length+0.5,rotor_length+thrust_disc_h+0.6+wending_length+0.5+core_g,
!下定子
RECTNG,core_r_1,core_r_2,rotor_length-0.6,rotor_length-0.6-wending_length-0.5,
RECTNG,core_r_3,core_r_4,rotor_length-0.6,rotor_length-0.6-wending_length-0.5,
RECTNG,core_r_1,core_r_4,rotor_length-0.6-wending_length-0.5,rotor_length-0.6-wending_length-0.5-core_g,


AADD,5,6,7
AADD,8,9,10

RECTNG,0,40,0,69,

AOVLAP,ALL,
NUMCMP, AREA !面积号压缩
APLOT
ALLSEL,ALL
ARSCALE,ALL,,,0.001,0.001,0.001,,0,1


ASEL,S,AREA, ,1,
!选择推力盘
AATT,4,1,1,0
ASEL,S,AREA, ,6,7,
!选择空气
AATT,1,1,1,0
ASEL,S,AREA, ,2,3,1,
!选择线圈
AATT,2,1,1,0
ASEL,S,AREA, ,4,5,1
!选择定子
ASEL,A,AREA, ,8,
!选择rotor
AATT,3,1,1,0
/PNUM,MAT,1
! 打开材料编号
ALLSEL,ALL
APLOT
!!!!!!!!!!!!!!MESHING!!!!!!!!!!!
smrtsize,1
mshkey,0
mshape,1,2d
Amesh,all
!定义平行边界条件
ALLSEL,ALL
NSEL,EXT
D,ALL,AZ,0
ALLSEL,ALL
!加力标志
ASEL,S,,,1
ESLA,S

CM, FORCE_B,ELEM
FMAGBC,'FORCE_B'
ALLSEL,ALL
FINISH


!load current
/SOLU
ASEL,S, , ,2
ESLA,S

BFE,ALL,JS,1,,,(I+iz)*turns/( area_coil *0.01**2)
ASEL,S, , ,3
ESLA,S

BFE,ALL,JS,1,,,-(I-iz)*turns/( area_coil *0.01**2)


MAGSOLV,0,3,0.001,,25
FINISH

/POST1

FMAGSUM,'FORCE_B'