三维棱边单元法静磁场分析结果误差太大

hbbdyzf

我用三维棱边单元法求解静磁场，求救结果与实际值相差很大，命令改了好几天，都是这样，不知道什么原因。求高手解答，下面是命令流

/prep7
emunit,mks
h1=0.123
h2=0.047
h3=0.027
h4=0.010
t=0.060
l1=0.057
l2=0.129
i=85
n=435
jx=0
jy=0
jz=1435365
*afun,deg
et,1,solid117
mp,murx,1,1
mp,murx,2,1
mp,rsvx,1,1e-8
mp,murx,3,4000
local,11,0,0,0,-0.833
wpcsys,11
k,1,-0.075,0.415,0
k,2,-0.365,0.235,0
k,3,-0.405,0.010,0
k,4,-0.570,0.010,0
k,5,-0.570,-0.770,0
k,6,-0.345,-1.020,0
k,7,0.345,-1.020,0
k,8,0.570,-0.770,0
k,9,0.570,0.010,0
k,10,0.405,0.010,0
k,11,0.365,0.235,0
k,12,0.075,0.415,0
lstr,1,2
lstr,2,3
lstr,3,4
lstr,4,5
lstr,5,6
lstr,6,7
lstr,7,8
lstr,8,9
lstr,9,10
lstr,10,11
lstr,11,12
lstr,12,1
al,all
voffst,1,0.260
wpoff,0,0,0.260
cyl4,0,0,0.288,360,,,0.165
local,12,0,0,0,-0.408
wpcsys,12
k,33,0,0,0
k,34,0.4275,0,0
k,35,0.4275,0,0.012
k,36,0.4315,0,0.012
k,37,0.480,0,0.067
k,38,0,0,0.067
lstr,33,34
lstr,34,35
lstr,35,36
lstr,36,37
lstr,37,38
lstr,38,33
al,47,48,49,50,51,52
vrotat,19,,,,,,33,38,360
vadd,3,4,5,6
numcmp,volu
local,13,0,0,0,-0.341
wpcsys,13
k,51,0,0,0
k,52,0,-0.480,0
k,53,0,-0.525,0.056
k,54,0,-0.525,h1
k,55,0,-0.525+t,h1
k,56,0,-0.525+t,h1-h2
k,57,0,-0.525+t+l1,h1-h2-h3
k,58,0,-0.525+t+l1+l2,h1-h2-h3
k,59,0,-0.525+t+l1+l2,h1-h2-h3+h4
k,60,0,0,h1-h2-h3+h4
lstr,51,52
lstr,52,53
lstr,53,54
lstr,54,55
lstr,55,56
lstr,56,57
lstr,57,58
lstr,58,59
lstr,59,60
lstr,60,51
al,52,84,85,86,87,88,89,90,91,92
vrotat,19,,,,,,51,60,360
vadd,4,5,6,7
numcmp,volu
local,14,1,0,0,-0.5605
wpcsys,14
cylind,0.3,0.484,0,0.14,360,360
csys,0
wpcsys,0
allsel,all
vsymm,z,all,,,,0,0
numcmp,volu
local,15,0,0,0,-0.573
wpcsys,15
k,201,0.005,-0.89826,0
k,202,0.46472,-0.77818,0  
k,203,0.52974,-1.02995,0
k,204,0.06499,-1.14995,0
lstr,201,202
lstr,202,203
lstr,203,204
lstr,204,201
al,341,342,343,344
voffst,161,-1.146
csys,15
wpcsys,15
k,209,-0.005,-0.89829,0
k,210,-0.46472,-0.77818,0
k,211,-0.52974,-1.02995,0
k,212,-0.06499,-1.14995,0
lstr,209,210
lstr,210,211
lstr,211,212
lstr,212,209
al,353,354,355,356
voffst,167,1.146
csys,0
wpcsys,0
block,-1.14,1.14,-2.58495,1.85,-2.499,2.499
allsel,all
vovlap,all
numcmp,volu
vsel,s,volu,,13
vatt,2,0,1
vsel,s,volu,,3,12,1
vatt,3,0,1
local,16,1,0,0,0.5605
vsel,s,volu,,2
vatt,1,0,1
csys,14
vsel,s,volu,,1
vatt,1,0,1
csys,0
allsel,all
smrtsize,3
mshape,1,3d
mshkey,0
vmesh,all
eplot
csys,16
vsel,s,volu,,2
eslv
emod,all,esys,16
bfe,all,js,1,jx,jy,jz
csys,14
vsel,s,volu,,1
eslv
emod,all,esys,14
bfe,all,js,1,jx,jy,jz
csys,0
allsel
finish
/solu
antype,static
eqslv,sparse
csys,0
asel,s,area,,173,178,1
da,all,az,0
allsel,all
magsolv
save
finish
/post1
plnsol,b,sum

path,pat1,2
ppath,1,,0,0,-0.18
ppath,2,,0,0,0.18
pdef,line1,bz,sum
plpath,line1