用Mason計算針板電極場強大小的適用性分析

張 凱，孫紅梅

（1.國網遼寧省電力有限公司，遼寧 沈陽 110006；2.國網遼寧省電力有限公司檢修公司，遼寧 沈陽 110006）

用Mason計算針板電極場強大小的適用性分析

張 凱1，孫紅梅2

（1.國網遼寧省電力有限公司，遼寧 沈陽 110006；2.國網遼寧省電力有限公司檢修公司，遼寧 沈陽 110006）

在進行XLPE電纜絕緣材料電樹枝化試驗時，廣泛采用針板電極系統，研究了針電極的針尖形狀及尖端曲率大小是否滿足要求，并采用Mason公式計算針板電極系統最大場強是否適用、是否能正確反映針尖場強實際大小，結果證明試驗用針電極尖端形態和尺寸均滿足試驗設計要求及Mason公式的適用條件，利用Mason計算試驗用針板電極最大場強是可行的。

Mason；ANSYS；針板電極；電樹枝化

隨著城市和現代工業的迅速發展，XLPE電纜在我國城市電網中得到越來越廣泛的應用［1－2］，但同其他高分子聚合物絕緣材料一樣，在電場的長期老化作用下，電樹枝化已成為電纜絕緣失效的重要因素。

為解決由于電樹枝化導致絕緣失效而影響電纜正常運行的問題，各國學者對絕緣材料電樹枝化進行了研究［3－12］，希望從不同的角度全面了解電纜絕緣材料電樹枝化的特性。其研究試驗手段之一是利用針板電極，通過施加不同形式的電壓或將試樣放置在不同溫度下，進行電樹枝的引發和生長特性的研究，目前已取得一定的成果［13－16］。

研究發現引發材料中電樹枝化的原因很多，包括導體的凸起、毛刺、絕緣材料中的氣隙和應力等，通常采用針板電極系統模擬上述因素進行試驗研究。為方便研究結果對比，國內高校普遍采用曲率半徑為5 μm、針板電極距離為2 mm的針板電極系統。在進行針板電極系統中針電極尖端電場強度計算時，均采用Mason公式進行估算。

在進行XLPE電纜絕緣直流電樹枝化試驗時，發現即使施加的直流電壓值達到70 kV，利用Mason公式估算的最大場強為2 844.22 kV／mm，遠高于材料的本征擊穿電壓800 kV／mm時，仍未能成功引發直流電樹枝。在已有的研究結果基礎上［17］，為找出造成XLPE電纜絕緣材料在直流電壓下電樹枝化難以引發的直接原因，研究是否由于Mason公式計算試驗條件下針尖場強的誤差增大及是否由于針板電極系統中針電極尖端形態不符合要求而導致試驗結果的出現。

目前已有文獻對用Mason公式計算尖電極場強進行誤差研究［18］。本文采用ANSYS有限元仿真軟件，對Mason公式是否能用來計算設定試驗條件下針板電極系統電場強度大小進行研究，研究Mason公式是否可以用于估算試驗用針板電極的最大場強。同時利用Matlab對實際針電極的形狀和尖端曲率半徑進行提取，研究試驗所用針電極形狀和曲率半徑是否滿足試驗設定條件及Mason公式使用條件的要求。

1 Mason和ANSYS計算結果對比

1.1 Mason公式

J.H.Mason在1951年提出針板電極系統最大場強的近似計算公式［19］。我國倪光正等從理論出發，詳細推導了標準雙曲面最大場強計算過程［20］，并最終得到了Mason公式在長旋轉橢球坐標系中的電位、尖端最大場強的表達式及直角坐標系中尖端最大場強的表達式。推導前提是為標準旋轉雙曲面的電位函數建立數學模型，并歸結為如下的邊值問題：

通過進一步推導，得出電位函數的解為

針尖最大場強表達式為

式中 φ——電位函數；

η——長旋轉橢球坐標；

φ0——電極間的外施電壓；

Emax——針尖最大場強；

h——針尖到平面的距離。

考慮到實際針電極表面不完全符合雙曲面要求，但在計算中可將針尖端部位近似成曲率為R的雙曲面，考慮到R?h，可將式（3）簡化得到Mason公式：

式中 R——曲率半徑；

H——針板電極間距；

U0——試驗所加電壓。

由上述推導過程可知，Mason公式可適用于尖端輪廓為標準雙曲面的針板電極系統。

1.2 用ANSYS計算及建模驗證

為實現ANSYS軟件進行后續計算仿真工作，首先驗證ANSYS建模、加載荷及計算的正確性，即用ANSYS對標準旋轉雙曲面針板電極系統場強分布進行計算和仿真，并將仿真所得結果同Mason公式計算結果進行比較。

為使計算更具有實用性，能和實際情況相結合，在ANSYS建模時，取針尖曲率半徑R＝5 μm，針板電極間距d＝2 mm。同時為了建模簡單，優化資源利用，將實際中的試樣進行簡化，如圖1所示。針電極模型為其尖端部位，板電極模型為相對針尖尺寸較大的平板，電極間區域填充的介質為XLPE，分別在針電極與板電極所在區域加電壓載荷后進行有限元計算。

圖1 ANSYS建模簡化

針電極建模時，采用軸對稱結構特有的簡化算法，即將所用單元的關鍵項設為軸對稱方式，然后建立1／4體積截面的針尖模型進行仿真計算。針尖是用標準雙曲方程繪制而成，如式（5）所示：

利用式（6）對標準雙曲形狀的針電極針尖進行繪制，然后完善模型，包括繪制平板電極、半圓形介質模型和圓環形的無限遠場，最終模型如圖2所示。

在進行計算時，選擇的分析類型為靜態電場分析，單元類型為PLANE121，開啟軸對稱選項，為了使計算結果更具有實用性，所加載荷為試驗所得的XLPE電纜50%絕緣雷電擊穿電壓為35 kV，選用ANSYS默認求解器，求解結果如圖3所示。

圖2 ANSYS建模

圖3 雙曲針尖仿真結果

由圖3可見，雙曲針板電極系統中，針尖曲率R＝5 μm，針板距離d＝2 mm，施加電壓U0＝35 kV時，仿真得出最大場強分布在針尖處，仿真計算結果為Emax＝1 971.96 kV／mm。在相同的條件下，即針板距離d＝2 mm，曲率半徑R＝5 μm，施加的電壓U0＝35 kV，代入Mason公式進行計算，Emax＝1 897.43 kV／mm，結果相對誤差為3.92%，可以認為Mason公式和ANSYS計算結果是相同的。由于Mason公式是通過純理論推導得出，完全適用于標準雙曲面，其計算結果能反映電壓值作用在針電極尖端處場強的大小，對比仿真結果可知，本文所采用的ANSYS建模、加載荷及計算方法切實可行，可用于計算針板電極下針尖的最大場強。

2 實際針尖模擬

2.1 提取實際針電極形態

為確定用Mason公式估算實際針板電極針尖場強的合理性，用ANSYS對XLPE電樹枝化試驗所用針電極尖端進行仿真計算。在ANSYS中建立的實際針尖模型數據是通過Matlab對針電極尖端圖像形狀提取得到的，Matlab中針電極尖端圖像形狀提取步驟如圖4所示。

得到針尖輪廓點的坐標后，在Matlab軟件中選取擬合結果較好的擬合函數進行擬合，經過嘗試發現用8級Fourier函數進行擬合較合理。擬合函數為

圖4 實際針尖擬合流程

將擬合所得結果與原數據點進行比較，標準方差為2.05，結果如圖5所示。根據曲率半徑的定義，求取擬合后的針電極曲線在針尖處的曲率半徑R＝6.988 μm，基于擬合過程中的誤差，由圖5（c）可以看出，擬合曲線與原輪廓線相比，曲率相對較大，故認為實際針尖處的曲率半徑R≈5 μm。

2.2 實際針電極和標準雙曲線形態對比

針尖輪廓曲線和標準雙曲線對比如圖6所示，在相同曲率半徑下，實際針尖形狀不能完全符合雙曲線的形狀，但在尖端局部有很好的重合，計算針電極尖端最大場強時，由于針尖電極尖端局部形狀滿足標準雙曲線要求，采用Mason公式進行針電極最大場強的估算是可行的。

2.3 用ANSYS計算實際針電極場強大小

用式（7）的擬合函數，繪制ANSYS模型中的針電極，進行仿真計算，所涉及的建模、加載荷和計算流程與上述標準雙曲模型仿真計算過程一致，結果如圖7所示。從仿真結果看，實際針電極的最大場強分布在針尖處，且針尖最大場強Emax＝1 486.05 kV／mm。該場強值是在只考慮施加電壓作用下，針電極尖端處的實際場強大小。

用Mason公式估算實際場強值，將R＝6.988 μm，d＝2 mm代入式（4），得Emax＝1 422.11 kV／mm。將該值與ANSYS仿真計算結果Emax＝1 486.05 kV／mm進行比較，兩者相差4.49%。從而可以看出，用Mason公式估算電樹枝化試驗用針板電極系統的場強時，其估算結果能切實反映所施加電壓作用在針電極尖端處產生場強的大小。

圖5 實際針尖提取過程

圖6 輪廓對比

圖7 實際針尖的ANSYS仿真結果

3 結論

利用ANSYS仿真分析了標準雙曲面和實際針尖形狀曲線的最大場強，并和Mason計算結果比較，發現ANSYS仿真實際針尖的結果同Mason計算結果相近，因此，可用Mason公式對實際使用的針電極尖端處的最大場強進行估算，當曲率半徑為5 μm，針板距離為2 mm時，用Mason計算針電極最大場強的大小所引起的誤差在可接受范圍內。試驗結果同時說明，實際使用的針電極尖端不滿足標準雙曲面的要求，但在尖端局部能與標準雙曲面很切合，從而滿足了Mason計算條件。結果還說明在實際電樹枝試驗中［14］，當在試樣的針電極上施加35 kV直流電壓時，針電極尖端處場強Emax可高達1 897.43 kV／mm，此值遠遠高于XLPE電纜絕緣材料的擊穿場強800 kV／mm，但實際中并沒有能夠成功引發電樹枝，這證明正是由于同極性空間電荷的作用削弱了針尖處場強，使試樣沒有引發電樹枝。

［1］姜小兵.XLPE電纜線芯溫度計算方法研究［J］.東北電力技術，2014，35（3）：1－3.

［2］金 鑫，李書碩，王東燁.交聯電纜0.1 Hz試驗方法研究［J］.東北電力技術，2009，30（9）：5－7.

［3］FH Kreuger.Partial discharge detection in high－voltage equipment［M］.London：Butterworths，1989.

［4］LA Dissado，JC Fothergill.Electrical degradation and breakdown in polymers［M］.London：Peter Peregrinus，1992.

［5］J Densley.Ageing mechanisms and diagnostics for power cables－an overview［J］.IEEE Electrical Insulation Magazine，2001，17（1）：14－22.

［6］J Densley，T Kalickil，Z Nadolnf.Characteristics of PD pulses in electrical trees and interfaces inextruded cables［J］.IEEE Trans. Dielectrics and Electrical Insulation，2001，8（1）：48－57.

［7］Fukuda T.Technological progress in High voltage XLPE power cable in Japan［J］.Electrical Insulation Magazine，1998，4（6）：15－20.

［8］Muller K B，Tretow U，Dellby B，et al.XLPE cable Technology for operating voltage of 245 kV and above［J］.IEEE Transaction on Power Delivery，1990，5（4）：1 660－1 668.

［9］李盛濤，鄭曉泉.聚合物電樹枝化［M］.北京：機械工業出版社，2006.

［10］陳仕軍.不同條件下交聯聚乙烯電纜電樹枝生長特性研究［D］.重慶：重慶大學，2009.

［11］石 峰.交聯聚乙烯電纜耐壓試驗方法［J］.東北電力技術，2010，31（1）：47－49.

［12］張秀閣，賀景亮，王洪新，等.工頻電壓下XLPE絕緣的電樹枝老化實驗研究［J］.高電壓技術，1998，24（4）：57－58.

［13］鄭曉泉，G CHEN.XLPE電纜絕緣中的電樹枝種類及其影響因素［J］.電工電能新技術，2003，22（4）：21－24.

［14］劉 英，曹曉龍.直流電壓下交聯聚乙烯電纜絕緣中電樹枝的生長特性研究［J］.西安交通大學學報，2014，48（4）：41－46.

［15］鮑明暉，尹小根.高頻電壓下交聯聚乙烯中電樹枝的形態特性［J］.中國電機工程學報，2011，31（34）：184－191.

［16］陳向榮，徐 陽.高溫下110 kV交聯聚乙烯電纜電樹枝生長及局部放電特性［J］.高電壓技術，2012，38（3）：645－654.

［17］Liu ying，Cao xiaolong.Electrical tree initiation in XLPE cable insulation by application of DC and impulse voltage，IEEE Transactions on Dielectrics and Electrical Insulation，2013，20（5）：1 691－1 698.

［18］文康珍，文元芳，端木林楠，等.尖電極靜電場ANSYS解的誤差［J］.電瓷避雷器，2009，43（2）：14－16.

［19］J.H.Mason：The Deterioration and Breakdown of Dielectrics Resulting from Internal Discharges，PIEE，1951，Vol98.Pt I.

［20］倪光正，吳 炯.Mason公式的導出及其在樹枝放電研究中的應用［J］.西安交通大學學報，1984，18（4）：87－93.

Suitability Analysis on Needle Plate Electrode Field Strength Based on Mason Calculation

ZHANG Kai1，SUN Hong?mei2
（1.State Grid Liaoning Electric Power Co.，Ltd.，Shenyang，Liaoning 110006，China；2.Maintenance Company of State Grid Liaoning Electric Power Co.，Ltd.，Shenyang，Liaoning 110006，China）

Learners commonly use the needle plate electrode system when studying the properties of the electrical tree in XLPE cable insulation.This paper studies that whether the profile of the needle electrode tip meets the hyperbolic shape or the curvature of needle electrode tip meets the requirements.This paper also studies that whether the theoretical analysis of the field strength in the needle plate electrode system with Mason formula can correctly describe the field strength at needle electrode tip in the experiment.The results show that the shape and size of the needle electrode tip can meet the requirements of design of experiment.The Mason formula can be used to calculate maximum field strength of the needle plate electrode system in the experiment.

Mason；ANSYS；Needle plate electrode；Electrical tree

TM74

A

1004－7913（2016）06－0027－04

張 凱（1985—），男，碩士，工程師，從事電力系統運行工作。

2016－03－20）