求助高手帮忙看看ansys二维 永磁铁 励磁线圈复合磁场仿真不收敛的问题，十分感激！

我做的是一个永磁铁和载流线圈共同作用的一个二维静态 电磁场分析，但是总是出不来结果，由于本人论文较急，还望有高手可以指点一下，小弟十分感谢！

这是机械结构图

这是有限元模型图，A4,A5为通入电流相反的线圈区域，A6，A7区域为磁极相反的永磁铁
下面是我的有限元仿真程序

/prep7

ET,1,PLANE13 ! 指定单元

KEYOPT,1,3,1 ! 指定分析类型为轴对称

EMUNIT,MKS !单位制

TB,BH, 1, , 14 ! 磁流变液BH曲线

TBPT,DEFI, 18000, 0.125

TBPT,, 35000, 0.25

TBPT,, 47000, 0.375

TBPT,, 61000, 0.5

TBPT,, 85000, 0.625

TBPT,, 120000, 0.75

TBPT,, 155000, 0.875

TBPT,, 200000, 1

TBPT,, 255000, 1.125

TBPT,, 330000, 1.25

TBPT,, 419000, 1.375

TBPT,, 502000, 1.5

TBPT,, 600000, 1.625

TBPT,, 705000, 1.75

TB,BH, 2, , 19 ! 缸体、活塞杆(45#钢)BH曲线

TBPT,DEFI, 410, 0.05

TBPT,, 790, 0.1

TBPT,, 1150,0.2

TBPT,, 1450,0.3

TBPT,, 1650,0.4

TBPT,, 1850,0.5

TBPT,, 2200,0.6

TBPT,, 2420,0.7

TBPT,, 2762,0.8

TBPT,, 2950,0.9

TBPT,, 3220,1

TBPT,, 3680,1.1

TBPT,, 4380,1.2

TBPT,, 5150,1.3

TBPT,, 6300,1.4

TBPT,, 8400,1.5

TBPT,, 11580, 1.6

TBPT,, 16000, 1.7

TBPT,, 24800, 1.85

TB,BH, 3, , 27 ! 活塞(电工纯铁) BH曲线

TBPT,DEFI,33,0.1

TBPT,, 71.8,0.22

TBPT,, 101.1,0.31

TBPT,, 150,0.46

TBPT,, 179.4,0.55

TBPT,, 212,0.65

TBPT,, 248,0.76

TBPT,, 298.2,0.89

TBPT,, 386.5,1.02

TBPT,, 513.6,1.15

TBPT,, 665.3,1.27

TBPT,, 753.6,1.33

TBPT,, 911.5,1.41

TBPT,, 1119.6,1.45

TBPT,, 1471.3,1.53

TBPT,, 1884,1.57

TBPT,, 2334.8,1.6

TBPT,, 2888.8,1.63

TBPT,, 4039,1.67

TBPT,, 5206.3,1.7

TBPT,, 83353,1.76

TBPT,, 10048,1.79

TBPT,, 13098,1.84

TBPT,, 15850,1.88

TBPT,, 21044,1.94

TBPT,, 35168,2.06

TBPT,, 37680,2.0

MP,MURX,4,1 ! 线圈

MP,MURX,5,1

MP,MURX,6,1.415 ! 永磁铁相对磁导率

MP,MGXX,6,-450e3 ! 永磁铁矫顽力

MP,MURX,7,1.415

MP,MGXX,7,450e3

/PNUM,AREA,1 ! 建模

RECTNG,0,0.06,0,0.009

RECTNG,0,0.06,0.009, 0.009+0.0045

RECTNG,0,0.01,0.009+0.0045, 0.057/2

RECTNG,0.01,0.05, 0.009+0.0045, 0.009+0.0045+0.004

RECTNG,0.01,0.05, 0.009+0.0045+0.004, 0.009+0.0045+0.004+0.0095

RECTNG,0.05,0.06, 0.009+0.0045, 0.057/2

RECTNG,0,0.06, 0.06/2, 0.074/2

ARSYM,x,all ! 镜像

AADD,1,8

AADD,6,2,3,10,9,13

AADD,7,14

AGEN,2, all, , , 0.06, , , , 1, 1 ! 坐标轴移动到端部

RECTNG,0,0.12,0, 0.074/2

AOVLAP,all

NUMCMP,AREA

APLOT

asel,s,,,7 ! 赋予材料
磁流变液

AATT,1, 1, 1

asel,s,,,1 ! 赋予材料
缸体，活塞杆

asel,a,,,6

AATT,2, 1, 1

asel,s,,,8 ! 赋予材料
活塞

AATT,3, 1, 1

asel,s,,,2, ! 赋予材料
线圈

AATT,4, 1, 1

asel,s,,,4,

AATT,5, 1, 1

asel,s,,,3 ! 赋予材料
永磁铁

AATT,6, 1, 1

asel,s,,,5

AATT,7, 1, 1

/PNUM,MAT,1 ! 打开材料编号

allsel,all

APLOT

SMRTSIZE,2 ! 智能网格等级

AMESH,ALL ! 划分网格

FINISH

/solu

n=325 ! 线圈匝数

i1=1 ! 电流

wc=0.04

hc=0.0095

acoil=wc*hc

jdens1=n*i1/acoil ! 电流密度

ESEL,S,MAT,,4 ! 选择线圈单元

BFE,ALL,JS,,,,-jdens1 ! 电流载荷

ESEL,S,MAT,,5 ! 选择线圈单元

BFE,ALL,JS,,,,jdens1 ! 相反电流载荷

allsel,all

NSEL,EXT ! 选择外层所有节点

D,ALL,AZ,0 ! 施加磁力线平行边界条件

Magsolv,0,3,0.001,,25 ! 求解

FINISH

/POST1

PLF2D ! 磁力线图

PLVECT,B,, , ,VECT,ELEM,ON,0 ! 磁流密度矢量图

PLNSOL,B,SUM, 0 ! 磁通密度等值云图

程序有什么错误还希望高手指点。