『求助』 [电磁场] 关于衔铁受力的问题，高手们都来看看呀！没有明白！

图片上传不上来，不知道为什么!一个铁心加一个衔铁，最基本的计算,SOLID117单元建立模型!
我的整个程序如下：
finish
/clear
/PREP7
!********************************模型结构参数
!**************EE160A(铁心）
A=96e-3,
B=64e-3,
C=32e-3,
D=32e-3,
E=47e-3,
F=63e-3,
!**************EE160B（衔铁）
A1=A，
B1=B，
C1=C，
D1=D，
E1=15e-3,
F1=(F-E)+E1
!**************气隙
H=5e-3
!**************螺旋线圈
L=E-(4e-3) !螺旋线圈的高
K=(B-C)/2-(5e-3) !螺旋线圈的宽

!******************************************建立模型
!******************************建立铁心
block,0,A,0,D,0,F-E, !建立最下部梁
block,0,(A-B)/2,0,D,F-E,F, !建立边缘触点1
block,A-(A-B)/2,A,0,D,F-E,F, !建立边缘触点2
block,(A-B)/2+(B-C)/2,(A-B)/2+(B-C)/2+C,0,D,F-E,F,
vglue,all
numcmp,all !压缩所有编号
!******************************建立衔铁
kwpave,11,
csys,4,
block,0,(A1-B1)/2,0,D1,H,H+E1
block,A1-(A1-B1)/2,A1,0,D1,H,H+E1,
block,(A1-B1)/2+(B1-C1)/2,(A1-B1)/2+(B1-C1)/2+C1,0,D1,H,H+E1, !衔铁中心触点
block,0,A1,0,D1,H+E1,H+F1, !建立衔铁的部梁
vsel,s,loc,z,0,H+F1,
vglue,all
allsel,all
numcmp,all
!******************************建立线圈,线圈和铁心的间距固定2mm
kwpave,10,
block,(B-C)/2-(2e-3),(B-C)/2-(2e-3)-K,0,D1,(E-L)/2,(E-L)/2+L
kwpave,24,
block,(2e-3),(2e-3)+K,0,D1,(E-L)/2,(E-L)/2+L,
kwpave,21,
CYLIND,2e-3,(2e-3)+K,(E-L)/2,(E-L)/2+L,90,180,
kwpave,23,
CYLIND,2e-3,(2e-3)+K,(E-L)/2,(E-L)/2+L,0,90,
kwpave,21,
block,0,C,2e-3,(2e-3)+K,(E-L)/2,(E-L)/2+L,
vsel,s,volu,,9,
vsel,a,volu,,10,
vsel,a,volu,,11,
vsel,a,volu,,12,
vsel,a,volu,,13,
vglue,all
allsel,all
numcmp,all
!*******************************建立气隙
kwpave,11，
block,0,(A-B)/2,0,D,0,H
kwpave,25,
block,0,C,0,D,0,H
kwpave,17,
block,0,(A-B)/2,0,D,0,H
!*******************************建立空气层
wpcsys,-1,0
csys,0,
block,-(A-B)/2,A+(A-B)/2,0,D+(2e-3)+k+(A-B)/2,-(A-B)/2,F+(5e-3)+F1+(A-B)/2,
VOVLAP,ALL
NUMCMP,ALL
!*******************************模型的刨分
kwpave,10，
csys,4,
wpro,,,-90
vsel,s,volu,,4,
vsel,a,volu,,8,
vsbw,all,
allsel,all
kwpave,22,
vsel,s,volu,,20,
vsel,a,volu,,21,
vsbw,all,
allsel,all
kwpave,24,
vsel,s,volu,,22,
vsel,a,volu,,23,
vsbw,all,
allsel,all
kwpave,16,
vsel,s,volu,,24,
vsel,a,volu,,25,
vsbw,all,
allsel,all
!********************************给定材料参数
!***********************铁心材料参数
!铁心B-H曲线（50W400）
TB,BH,1,1,20,
TBTEMP,0
TBPT,,30,0.095
TBPT,,50,0.14
TBPT,,70,0.21
TBPT,,100,0.34
TBPT,,200,0.82
TBPT,,300,1.02
TBPT,,500,1.11
TBPT,,700,1.24
TBPT,,1000,1.3
TBPT,,2000,1.38
TBPT,,3000,1.42
TBPT,,5000,1.46
TBPT,,7000,1.49
TBPT,,10000,1.51
TBPT,,20000,1.58
TBPT,,30000,1.62
TBPT,,50000,1.68
TBPT,,70000,1.72
TBPT,,100000,1.76
TBPT,,300000,1.91
!衔铁B-H曲线（50W400）
*CSET,1,2, 1,2,
MPCOPY, ,1,2
TBCOPY,ALL,1,2
!空气
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,3,,1
!线圈
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,4,,1
!***********************************给定单元属性

et,1,solid117, !铁心
et,2,solid117, !衔铁
et,3,solid117, !空气
et,4,solid117, !线圈
!***********************************给定材料参数
!铁心
vsel,s,volu,,1,
vsel,a,volu,,2,
vsel,a,volu,,4,
vsel,a,volu,,18,
vsel,a,volu,,20,
vsel,a,volu,,23,
vsel,a,volu,,22,
vsel,a,volu,,3,
vatt,1,,1,
allsel,all
!衔铁
vsel,s,volu,,5,
vsel,a,volu,,19,
vsel,a,volu,,8,
vsel,a,volu,,21,
vsel,a,volu,,27,
vsel,a,volu,,26,
vsel,a,volu,,6,
vsel,a,volu,,7,
vatt,2,,2,
allsel,all
!线圈
vsel,s,volu,,9,
vsel,a,volu,,10,
vsel,a,volu,,11,
vsel,a,volu,,12,
vsel,a,volu,,13,
vatt,4,,4,
allsel,all
!空气
vsel,s,volu,,17，
vsel,a,volu,,14,
vsel,a,volu,,16,
vsel,a,volu,,15,
vatt,3,,3,
allsel,all
!************************************划分网格
!铁心和衔铁
VSEL,S,MAT,,1
vsel,a,mat,,2
ASLV,S
LSLA,S
lesize,all,5e-3,
allsel,all
VSEL,S,MAT,,1
vsel,a,mat,,2
MSHAPE,0,3D
MSHKEY,1
vmesh,all
!线圈
vsel,s,volu,,9,
vsel,a,volu,,10,
vsel,a,volu,,13,
ASLV,S
LSLA,S
lesize,all,5e-3,
allsel,all
vsel,s,volu,,9,
vsel,a,volu,,10,
vsel,a,volu,,13,
MSHAPE,0,3D
MSHKEY,1
vmesh,all
allsel,all
lsel,s,line,,132,
lsel,a,line,,129,
lsel,a,line,,135,
lsel,a,line,,134,
lesize,all,5e-3,
allsel,all
vsel,s,volu,,11,
vsel,a,volu,,12,
MSHAPE,0,3D
MSHKEY,1
vmesh,all
allsel,all
!空气
lsel,s,line,,153,
lsel,a,line,,154,
lsel,a,line,,155,
lsel,a,line,,156,
lsel,a,line,,157,
lsel,a,line,,158,
lsel,a,line,,159,
lsel,a,line,,160,
lsel,a,line,,161,
lsel,a,line,,162,
lsel,a,line,,163,
lsel,a,line,,164,
lesize,all,1e-3,
allsel,all
vsel,s,volu,,16,
vsel,a,volu,,14,
vsel,a,volu,,15,
MSHAPE,0,3D
MSHKEY,1
vmesh,all
allsel,all
vsel,s,volu,,17
MSHKEY,0
MSHAPE,1,3d
vmesh,all,
allsel,all
!***************************************给通量定边界条件
WPCSYS,-1,0
asel,s,loc,y,0,
nsla,s,1,
d,all,az,0,
allsel,all
asel,s,area,,77,
asel,a,area,,78,
asel,a,area,,79,
asel,a,area,,80,
asel,a,area,,81,
nsla,s,1,
d,all,az,0,
allsel,all
!**************************************定义衔铁组件
vsel,s,mat,,2, !衔铁的材料属性为2
eslv,s,
cm,arm,elem, !创建衔铁组件
fmagbc,'arm' !施加力标志
allsel,all
allsel,all
!***********************************施加电流密度
csys,0,
acond=L*K
i=8
n=(acond/0.5e-6)*0.9
vsel,s,volu,,9,
ESLV,S
BFE,all,JS,1,,-i*n/acond,, ,
allsel,all
vsel,s,volu,,10,
ESLV,S
BFE,all,JS,1,,i*n/acond,, ,
allsel,all
vsel,s,volu,,13
ESLV,S
BFE,all,JS,1,-i*n/acond,,, ,
allsel,all
kwpave,21，
CSWPLA,11,1,1,1, !在点21处建立柱坐标
vsel,s,volu,,11
ESLV,S
emodif,all,esys,11,
BFE,all,JS,1,,i*n/acond,,,
allsel,all
kwpave,23，
CSWPLA,12,1,1,1, !在点23处建立柱坐标
vsel,s,volu,,12
ESLV,S
emodif,all,esys,12,
BFE,all,JS,1,,i*n/acond,,,
allsel,all
finish