《ansys电磁学有限元分析从入门到精通》三侧向测井仪的问题

fly名不好取

！这是该书上的一个例子，程序运行了很多遍出现了两个大问题（出现出有注释），程序如下：
/fillname,logwell_K
keyw,magelc,1
/prep7
r0=0.0445
k,1,0,-35
k,2,0,-1.8
k,3,r0,-1.8
k,4,r0,-0.7
k,5,r0,-0.5
k,6,r0,-0.1
k,7,r0,-0.075
k,8,r0,0.075
k,9,r0,0.1
k,10,r0,0.5
k,11,r0,0.7
k,12,r0,1.8
k,13,0,1.8
k,14,0,35
k,15,0.1,35
k,16,0.1,-35
a,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16
rectng,0.1,35,1,35
rectng,0.1,35,-1,1
rectng,0.1,35,-1,-35
rectng,0.1,50,1,51
rectng,35,50,-1,1
rectng,0.1,50,-1,-51
rectng,0,0.1,35,51
rectng,0,0.1,-35,-51
aovlap,all
aglue,all
nummrg,all
numcmp,all
/pnum,area,1
/replot
et,1,67,,,1
mp,rsvx,1,1
mp,rsvx,2,1
asel,s,area,,1,3,2
asel,a,area,,5,9
aatt,1,,1
allsel
asel,s,area,,2,4,2
aatt,2,,1
allsel
lesize,1,,,80,0.05,,,,1
lesize,2,,,4
lesize,3,,,80
lesize,4,,,15
lesize,5,,,35
lesize,6,,,2
lesize,7,,,12
lesize,8,,,2
lesize,9,,,35
lesize,10,,,15
lesize,11,,,80
lesize,12,,,4
lesize,13,,,80,20,,,,1
lesize,14,,,1
lesize,15,,,1
lesize,24,,,1
lesize,27,,,1
lesize,16,,,85,20,,,,0
lesize,28,,,85,0.05,,,,0
lesize,29,,,14,,,,,0
lesize,17,,,14,,,,,0
lesize,18,,,85,20,,,,0
lesize,32,,,85,0.05,,,,0
lesize,33,,,14,,,,,0
lesize,34,,,14,,,,,0
lesize,30,,,35,,,,,1
lesize,20,,,14,,,,,0
lesize,19,,,14,,,,,0
lesize,21,,,14,,,,,0
lesize,38,,,14,,,,,0
lesize,23,,,7,,,,,0
lesize,35,,,7,,,,,0
lesize,36,,,7,,,,,0
lesize,26,,,7,,,,,0
lccat,19,20
lccat,17,29
lccat,33,34
lccat,21,38

mshape,1,2d
mshkey,2
amesh,3,5,2

mshape,1,2d
mshkey,0
amesh,4

mshape,1,2d
mshkey,1
amesh,7,8

mshape,1,2d
mshkey,1
amesh,6

mshape,1,2d
mshkey,1
amesh,1,2

smrtsize,1
mshape,1
mshkey,0
amesh,9

lsel,,,,7
nsll,s,1
cp,1,volt,all
allsel

lsel,,,,9
lsel,a,,,5
nsll,s,1
cp,2,volt,all
allsel

lsel,,,,11
lsel,a,,,3
nsll,s,1
cp,3,volt,all
allsel

*dim,aa,,3
*vget,aa(1),kp,7,attr,node
*vget,aa(2),kp,9,attr,node
*vget,aa(3),kp,11,attr,node
finish

/solu
lsel,s,loc,x,50
lsel,a,loc,y,-51
lsel,a,loc,y,51
dl,all,,volt,0
allsel,all
fk,7,amps,1
fk,9,amps,0
fk,11,amps,0
sbctran
solve           !1、在加载荷后求解时，会出现327个警告，在错误报告中说 Finer mesh required near center line for element即在线单元划分时有问题
save         
finish

/post1
set,last
*dim,u1,,3
*vget,u1(1),node,aa(1),volt
*vget,u1(2),node,aa(2),volt
*vget,u1(3),node,aa(3),volt
*cfoen,u10,txt,d:\ansys\result,
*vwrite,u1(1)
(f10.8)
*cfclos
finish

/fillname,logwell_K
resum
fk,7,amps,0
fk,9,amps,1
fk,11,amps,0
solve
save
finish

/post1
set,last
*dim,u2,,3
*vget,u2(1),node,aa(1),volt
*vget,u2(2),node,aa(2),volt
*vget,u2(3),node,aa(3),volt
*cfoen,u20,txt,d:\ansys\result,
*vwrite,u2(1)
(f10.8)
*cfclos
finish

/fillname,logwell_K
resum
fk,7,amps,0
fk,9,amps,0
fk,11,amps,1
solve
save
finish

/post1
set,last
*dim,u3,,3
*vget,u3(1),node,aa(1),volt
*vget,u3(2),node,aa(2),volt
*vget,u3(3),node,aa(3),volt
*cfoen,u30,txt,d:\ansys\result,
*vwrite,u3(1)
(f10.8)
*cfclos
*dim,u1,,3
*dim,u2,,3
*vread,u1(1),u10,txt,d:\ansys\result,
(f10.8)
*vread,u2(1),u20,txt,d:\ansys\result,
(f10.8)
*dim,a,,2,2
*dim,b,,2
*dim,c,,2
a(1,1)=u2(1)-u2(2)
a(1,2)=u3(1)-u3(2)
a(2,1)=u2(1)-u2(3)
a(2,2)=u3(1)-u3(3)
b(1)=u1(2)-u1(1)
b(2)=u1(3)-u1(1)
*moper,c(1),a(1,1),solve,b(1)  !2、矩阵求解时显示Matrix is singular for SOLV command

fly名不好取

注释的内容在程序运行时没有输入，程序的内容也是在美式输入法下写出来的