我的3D电磁分析怎么没有结果?

我是个新手,编了个三维电磁分析程序,激励为线圈电流,频率为60Hz,计算结果显示没有任何显示,请高手帮帮忙,看是哪儿出错了.
单元采用SOLID97,介质模型为1/4圆柱体,内部有不具轴旋转对称的介质.
部分代码如下:
r,1,a,1,4*3.14*r3*a,,1
n1=node(r3,0,0)
esel,s,mat,,3
nsle
cp,1,curr,all
nsel,s,loc,x,r6
sf,all,inf
nsel,s,loc,y,z1
sf,all,inf
nsel,s,loc,x,0
d,all,az,0
allsel,all
finish
/solu
antype,harmic
esel,s,mat,,3
bfe,all,js,1,,,1/(2*a)
esel,all
harfrq,60
solv
finish
/post1
set,1
plvect,h
finish