原先一直都是通過截面圖來間接顯示實體表面的場強分布，確定最大場強的位置。由于只可能做有限數(shù)量的截面，所以不能保證得到值就一定是最大場強。后來采用的一個權(quán)宜的方法是，選擇附著在實體表面的單元，然后用 plesol,ef,sum 來畫出單元上的場強分布。可是，實體表面的單元參差不齊，畫出的圖很不好看，而且最大場強出現(xiàn)位置也可能被附著的單元所遮擋。當(dāng)時甚至還想到，是否可以在實體表面做一個平整的空氣薄層，然后再取出該薄層來畫場強分布。而這無疑會使計算量大大增加，因此也沒有采用。

實質(zhì)上，最期待得到的結(jié)果是：選出實體表面的節(jié)點，然后對這些節(jié)點采用 plnsol,ef,sum 畫場強圖。這時需要注意的是，由于不同材料分界面處場強不連續(xù)，分界面上節(jié)點的場強實際上有兩個值。為了畫出實體的表面場強，亦即緊貼實體的空氣層中的場強，必須在后處理時只讀入空氣層單元的計算結(jié)果。否則，將畫出的是實體內(nèi)部的場強。若該實體恰為導(dǎo)體，則其內(nèi)部場強為零，根本無法得到想要的結(jié)果。然而，按這個方法畫圖，實體表面的場強將被外部空氣域邊界的場強分布所遮擋，盡管我們只選了導(dǎo)體表面的節(jié)點。

到昨天才無意中發(fā)現(xiàn)，若將 PlotCtrls => Hidden Line Options 中的 graphics display method 選為 full model，對應(yīng)的宏命令為 /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

得到的小球表面場強如下：