我用LMATRIX命令算电感到底错在哪？能显示磁力线

xianzhang

/title,keti
/filname,keti,1
keyw,magnod,1
/prep7
ET,1,PLANE53           !空气
KEYOPT,1,3,1
ET,2,PLANE53           !线圈
KEYOPT,2,3,1
ET,3,INFIN110
KEYOPT,3,3,1            ! 远场
MP,MURX,1,1             !空气材料
MP,MURX,2,1               !线圈材料
MP,RSVX,2,1/(5.8*10E7)
!!!!!!!!!!!!!!!!!!!!
S=0.2
N1=500
N2=450
N3=400
R=1.5*S

*dim,wind,char,3


RECTNG,S,2*S,0,S
RECTNG,2.5*S,3.5*S,0,S
RECTNG,4*S,5*S,0,S
CYL4,0,0,0,0,18*S,90
CYL4,0,0,0,0,12*S,90
ASEL,ALL
AOVLAP,ALL

!!!!!!!!!!!!!!!!!!!!!!
/pnum,area,1                       !!!!进行剖分
aplot,all
numcmp,area
aplot,all

r,1,2*s*s,500,,1,1
r,2,2*s*s,450,,1,1
r,3,2*s*s,400,,1,1


asel,s,area,,1
aatt,2,1,2
smrtsize,2
mshape,0,2d
mshkey,0
amesh,1

asel,s,area,,2
aatt,2,2,2
smrtsize,2
mshape,0,2d
mshkey,0
amesh,2

asel,s,area,,3
aatt,2,3,2
smrtsize,2
mshape,0,2d
mshkey,0
amesh,3
allsel,all
aplot,all
SMRT,OFF
ASEL,S,AREA,,4
AATT,1,,3
LOCAL,11,1
WPCSYS,,11
LSEL,S,LOC,Y,45
LESIZE,ALL,,,10
ESIZE,,1
AMESH,4
ALLSEL,ALL
SAVE
APLOT,ALL
ASEL,S,AREA,,5
smrt,on
SMRTSIZE,5
MSHKEY,0
MSHAPE,1,2D
AMESH,5
allsel,all
finish
!!!!!!!!!!!!!!!!!!!!!!!!!!
                  !!!!!!!!!!!!!!!!!进行求解
/SOLU

cnvtol,csg,,1.0e-4!!!
eqslv,front

CSYS,0
NSEL,S,LOC,X,0
D,ALL,AZ,0
CSYS,1
NSEL,S,LOC,X,18*S
SF,ALL,INF
CSYS,0
ALLSEL,ALL
SAVE

asel,s,area,,1
esla,s
cm,wind1,elem
bfe,all,js,,,,1.2*500/(s*s)

asel,s,area,,2
esla,s
cm,wind2,elem
bfe,all,js,,,,1.5*450/(s*s)
!bfe,all,volt,,12

asel,s,area,,3
esla,s
cm,wind3,elem
bfe,all,js,,,,1.7*400/(s*s)

allsel,all
antype,static
magsolv
finish
!!!!!!!!!!!!!!!!!!!!!!!!!
/post1
allsel,all
plf2d
finish
/solu
*dim,cur,array,3

cur(1)=1.2
cur(2)=1.5
cur(3)=1.7
lmatrix,1,'wind','cur','ind'

*stat,ind

zhendfliuzp

你说说是什么错误啊？