ANSYS 中实体表面场强的绘制

原先一直都是通过截面图来间接显示实体表面的场强分布，确定最大场强的位置。由于只可能做有限数量的截面，所以不能保证得到值就一定是最大场强。后来采用的一个权宜的方法是，选择附着在实体表面的单元，然后用 plesol,ef,sum 来画出单元上的场强分布。可是，实体表面的单元参差不齐，画出的图很不好看，而且最大场强出现位置也可能被附着的单元所遮挡。当时甚至还想到，是否可以在实体表面做一个平整的空气薄层，然后再取出该薄层来画场强分布。而这无疑会使计算量大大增加，因此也没有采用。

实质上，最期待得到的结果是：选出实体表面的节点，然后对这些节点采用 plnsol,ef,sum 画场强图。这时需要注意的是，由于不同材料分界面处场强不连续，分界面上节点的场强实际上有两个值。为了画出实体的表面场强，亦即紧贴实体的空气层中的场强，必须在后处理时只读入空气层单元的计算结果。否则，将画出的是实体内部的场强。若该实体恰为导体，则其内部场强为零，根本无法得到想要的结果。然而，按这个方法画图，实体表面的场强将被外部空气域边界的场强分布所遮挡，尽管我们只选了导体表面的节点。

到昨天才无意中发现，若将 PlotCtrls => Hidden Line Options 中的 graphics display method 选为 full model，对应的宏命令为 /graphics,full，则可以仅显示实体表面的场强分布。这样，以后就无需麻烦地做截图，得到的最大场强值也更加准确。

用下面的命令流则可以作测试：

/prep7
et,1,solid123
mp,perx,1,1
mp,perx,2,1e5   

sph4,1,0,0.2
sph4,-1,0,0.2
sph4,0,0,10
vsel,all
vovlap,all
numcmp,all  

vsel,s,,,1,2
vatt,2,,1
smrtsize,1
esize,0.02
vmesh,all

vsel,s,,,3
vatt,1,,1
esize,0.5
vmesh,all

vsel,s,,,1
aslv,s
nsla,s,1
d,all,volt,10   

vsel,s,,,2
aslv,s
nsla,s,1
d,all,volt,-10   

asel,s,,,5,6
nsla,s,1
d,all,volt,0

/sol
antype,static,new
vsel,s,,,3
allsel,below,volu
solve
save

/post1
vsel,s,,,3
allsel,below,volu
set,1
/graphics,full
/show,win32c
/contour,1,128
vsel,s,,,1,2
aslv,s
nsla,s,1

plnsol,ef,sum

得到的小球表面场强如下：