2维轴对称电磁分析，电磁力求解问题
如图，部件1和4都是铁芯，2为空气，3为线圈。理论上，部件4对部件3只有X方向的作用力，而由于是轴对称结构，所以X方向的作用力抵消了。而部件4对部件3应该是没有Y方向的作用力的，但是我现在计算的结果是：部件1受到的Y方向的电磁力很大，不知道是为什么？
命令流如下：
/PREP7
N=930 !线圈匝数为930匝
I=0.5 !电流为0.5A
L=40 !线圈长度为40mm
W=25 !线圈宽度为25mm
ACOIL=W*L*0.000001 !线圈面积，单位m2
JDENS=N*I/ACOIL !线圈电流密度
ET,1,PLANE13
KEYOPT,1,3,1 !设置单元表现为轴对称
TB,BH,1
TBPT,,400,1.05
TBPT,,800,1.47
TBPT,,1200,1.51
TBPT,,1600,1.57
TBPT,,2000,1.63
TBPT,,2400,1.65
TBPT,,2800,1.68
TBPT,,3200,1.7
TBPT,,3600,1.72
TBPT,,4000,1.73
TBPT,,4400,1.75
TBPT,,4800,1.76
TBPT,,5200,1.77
TBPT,,5600,1.78
TBPT,,6000,1.8
TBPT,,6400,1.81
TBPT,,6800,1.82
TBPT,,7200,1.83
TBPT,,7600,1.84
TBPT,,8000,1.85
MP,MURX,2,1 !定义2号材料（线圈）的相对磁导率为1
MP,MURX,3,1 !定义3号材料（空气）的相对磁导率为1
MP,MURX,4,2000
K,1,0,0
K,2,29.5,0
K,3,29.5,5
K,4,24.5,5
K,5,24.5,10
K,6,29.5,10
K,7,29.5,15
K,8,24.5,15
K,9,24.5,20
K,10,29.5,20
K,11,29.5,25
K,12,24.5,25
K,13,24.5,30
K,14,29.5,30
K,15,29.5,35
K,16,24.5,35
K,17,24.5,75
K,18,29.5,75
K,19,29.5,80
K,20,24.5,80
K,21,24.5,85
K,22,29.5,85
K,23,29.5,90
K,24,24.5,90
K,25,24.5,95
K,26,29.5,95
K,27,29.5,100
K,28,24.5,100
K,29,24.5,105
K,30,29.5,105
K,31,29.5,110
K,32,0,110
K,33,29.8,0
K,34,29.8,5
K,35,34.8,5
K,36,34.8,10
K,37,29.8,10
K,38,29.8,15
K,39,34.8,15
K,40,34.8,20
K,41,29.8,20
K,42,29.8,25
K,43,34.8,25
K,44,34.8,30
K,45,29.8,30
K,46,29.8,35
K,47,34.8,35
K,48,59.8,35
K,49,59.8,75
K,50,34.8,75
K,51,29.8,75
K,52,29.8,80
K,53,34.8,80
K,54,34.8,85
K,55,29.8,85
K,56,29.8,90
K,57,34.8,90
K,58,34.8,95
K,59,29.8,95
K,60,29.8,100
K,61,34.8,100
K,62,34.8,105
K,63,29.8,105
K,64,29.8,110
K,65,73,110
K,66,73,0
L,1,2
L,2,3
L,3,4
L,4,5
L,5,6
L,6,7
L,7,8
L,8,9
L,9,10
L,10,11
L,11,12
L,12,13
L,13,14
L,14,15
L,15,16
L,16,17
L,17,18
L,18,19
L,19,20
L,20,21
L,21,22
L,22,23
L,23,24
L,24,25
L,25,26
L,26,27
L,27,28
L,28,29
L,29,30
L,30,31
L,31,32
L,32,1
L,2,33
L,31,64
L,33,34
L,34,35
L,35,36
L,36,37
L,37,38
L,38,39
L,39,40
L,40,41
L,41,42
L,42,43
L,43,44
L,44,45
L,45,46
L,46,47
L,47,48
L,48,49
L,49,50
L,50,51
L,47,50
L,51,52
L,52,53
L,53,54
L,54,55
L,55,56
L,56,57
L,57,58
L,58,59
L,59,60
L,60,61
L,61,62
L,62,63
L,63,64
L,64,65
L,65,66
L,33,66
LSEL,S,LINE,,33,48 !选择直线
LSEL,A,LINE,,2,30
LSEL,A,LINE,,52,66
AL,ALL !生成面
LSEL,S,LINE,,1,32
AL,ALL
LSEL,S,LINE,,35,52
LSEL,A,LINE,,54,69
AL,ALL
LSEL,S,LINE,,49,51
LSEL,A,LINE,,53
AL,ALL
AGLUE,ALL
ASEL,S,AREA,,3
AATT,1
ASEL,S,AREA,,4
AATT,2
ASEL,S,AREA,,1
AATT,3
ASEL,S,AREA,,2
AATT,4
ALLSEL,ALL
SMRTSIZE,4
AMESH,ALL
ESEL,S,MAT,,4
CM,ARM,ELEM
FMAGBC,'ARM'
ALLSEL,ALL
ARSCAL,ALL,,,0.001,0.001,1 !缩放模型，将单位从mm转换为m，原模型被删除
FINISH
/SOLU
ESEL,S,MAT,,2
BFE,ALL,JS,1,,,JDENS
ESEL,ALL
NSEL,EXT
D,ALL,AZ,0 !设置通量平行条件
ALLSEL,ALL
FINISH
/SOLU
MAGSOLV
SAVE
FINISH
/POST1
PLF2D
FMAGSUM
