一个关于磁感应强度计算

这是我做的一个计算磁感应强度样例
!/bach,list
/title,2d solenoid static analysis
/filename,mr_2d,1
keyw,magnod,1
/prep7
et,1,plane13
keyopt,1,3,1
mp,murx,1,1
tb,bh,2,,8
tbpt,,800,1.15
tbpt,,1200,1.35
tbpt,,1600,1.45
tbpt,,2800,1.6
tbpt,,4400,1.7
tbpt,,6800,1.8
tbpt,,16000,2
tbpt,,40000,2.2
mp,murx,3,1
mp,murx,4,1
N=938
I=0.612
D=18
TA=1.5
TB=9
TC=0.5
TD=4
L=38
LA=7
WC=2*LA+L
H=D-TB
acoil=H*L
jdens=N*I/acoil
/pnum,area,1
rectng,0,TD,0,WC
rectng,TD,TD+TC,0,LA
rectng,TD,TD+TC,LA,LA+L
rectng,TD,TD+TC,LA+L,WC
rectng,TD+TC,TD+TC+D,0,WC
rectng,TD+TC,TD+TC+TB,LA,LA+L
rectng,TD+TC+D,TD+TC+D+TA,0,WC
aovlap,all
numcmp,area
APLOT
save,mr_2d.db
asel,s,area,,5
aatt,3,1,1,0
asel,s,area,,6
asel,a,area,,7
aatt,2,1,1,0
asel,s,area,,4
aatt,1,1,1,0
asel,s,area,,1
asel,a,area,,2,3
aatt,4,1,1,0
/pnum,mat,1
allsel,all
aplot
smrtsize,1
amesh,all
arscal,all,,,0.001,0.001,1,,0,1
finish
/solu
asel,s,mat,,3
bfe,all,js,1,,,jdens/(0.001**2)
allsel,all
nsel,ext
D,all,az,0
allsel,all
magsolv,0,3,0.001,,25
save,mr_2d_resu.db
save
finish
/post1
plf2d,27,0,10,1
plvext,b,,,vect,elem,on,0
/graphics,power
avres,2
plnsol,b,sum
finish
向大家请教，有没有哪里做得不对！我是新手看着资料一点点编了这个
做完了一点信心也没有