圓柱形永磁體磁場建模及仿真研究

周恩權 鄭仲橋 張燕紅 王奇瑞

（1.海安交睿機器人科技有限公司，江蘇 南通 226000；2.常州工學院，江蘇 常州 213032；3.江蘇大學現代農業裝備與技術教育部重點實驗室，江蘇 鎮江 212013）

圓柱形永磁體磁場建模及仿真研究

周恩權1鄭仲橋2張燕紅2王奇瑞3

（1.海安交睿機器人科技有限公司，江蘇 南通 226000；2.常州工學院，江蘇 常州 213032；3.江蘇大學現代農業裝備與技術教育部重點實驗室，江蘇 鎮江 212013）

為了尋求圓柱形永磁體磁場的簡便算法，基于永磁體等效電流模型，把永磁體被磁化的效果視為永磁體表面存在面電流產生的磁場，得到了永磁體磁場的空間分布規律，推導出圓柱體型永磁體空間的磁感應強度解析公式，并通過MATLAB多重數值積分功能求出圓柱形永磁體的空間磁感應強度。結果表明：計算值和實驗測量值基本吻合，用永磁體等效模型及磁場數值積分方法計算永磁體的磁感應強度不僅簡單而且計算精度較高。

永磁體；磁場；數值積分；磁感應強度

1 永磁體等效電流模型的建立

1.1 永磁體電流等效模型

物質是由原子組成的，每個原子又由原子核和電子組成。電子繞原子核轉動形成電流，這些環流定向排列起來，在宏觀上顯示出N和S極。當磁介質均勻時，介質內部任何兩個分子環流中相鄰的那一對電流元方向總是彼此相反，它們的效果相互抵消。因此，永磁體的等效電流模型可以認為是，永磁體外部空間磁場是由永磁體側表面的束縛面電流產生的。其束縛面電流密度為：

其中，Jm表示束縛面電流密度；Br表示剩余磁化強度；μ0表示真空磁導率，μ0=4π×10-7H/m。

1.2 圓柱形永磁體三維磁場解析式推導

當半徑為r0、高度為h的圓柱形永磁體沿一個方向均勻磁化且達到飽和狀態后，外部空間中任意一點的磁場僅由永磁體表面閉合電流環路所激發。設環路內電流為I，則在與平面xoy平行的任意一平面上電流密度由于是圓柱形永磁體，選用柱坐標系，用()r0,θ,z0表示位于永磁體內部的點，建立如圖1所示坐標系。

圖1 圓柱永磁體的參數

考慮到由平面Z與z+dz構成的薄層，在其外部空間任一點p(x,y,z)處，由強度為Jmdz的電流環產生的磁場為dB→，從畢奧-薩伐爾定律和磁場疊加原理可知，圓柱體外任一點p(x,y,z)的磁感應強度為：

電流環產生磁感應強度B→在x,y,z的分量分別為：

圓柱永磁體空間磁感應強度為：

將（1）代入（6）得

1.3 MATLAB計算磁感應強度

已知函數的積分表達式時，理論上可以用牛頓萊布尼茨公式求解，但在磁場計算中并不實用，大多數函數找不到積分函數。MATLAB提供很多數值積分的命令函數，可以根據實際情況調用，不但方便而且精確。表1列出了部分數值積分命令。

表1 積分命令表

本文用到的積分命令是dblquad，dblquad實際上是重復調用quad實現的。quad的原理是自適應辛普生求積法。如果用二次插值多項式——拋物線y=g（x）所圍成的曲邊梯形的面積近似代y=f（x）所圍成的曲邊梯形的面積，這時所得的積分公式稱為辛普生公式：

當上述變步長辛普生積分法用于計算二重積分時，數值積分的處理方法是將二重積分（10）分解為式（11）和式（12）：

在MATLAB環境下編程，將（7）中的模型轉化為fun（θ，Z0）的M文件，通過B=dblquad（fun，a，b，c，d，tol）來實現模型的求積。其中，a和b為θ的上下限，c和d為Z0的上下限，tol為積分精度（默認精度10e-6）。

2 實例求解及實驗驗證

2.1 模型驗證

根據以上推導的圓柱永磁體空間磁場強度的模型，在MATLAB環境下編程仿真。仿真和實驗用的永磁體各項參數如表2所示。

表2 圓柱永磁體材料和尺寸參數

以圓柱永磁體底面圓心為坐標原點，圓柱的中軸為Z軸，建立如圖2所示的坐標系。在永磁體上表面Z軸上取如表3所示的12個場點，通過計算、測量及和磁偶極子模型的結果進行比較。計算和測量結果如表3。

表3 測量值和計算值比較

當z=45mm時，在X軸上取如表4所示的12個點，計算和測量值如表4所示。

表4 測量和計算值比較表

由表2和表3可知，本文提出的模型的計算結果和測量值吻合。磁偶極子模型的計算結果和實際值相差較大，但能反映磁感應強度的變化趨勢。測量值略小于本文模型計算值，主要原因是理論剩磁大于實際剩磁［13］。本文提出的方法能滿足工程計算需求。

2.2 圓柱形永磁體空間磁場仿真研究

圖2給出了圓柱永磁體Bz隨離開上表面距離變化的曲線。由圖2可知，從永磁體上表面到z=60mm的范圍內，Bz值迅速下降；當z＞60mm時，Bz的值下降緩慢；圖3是當Z=25mm、Y=0mm時Bz、BX沿X軸的分布規律圖。Bz為軸對稱分布而BX為中心對稱分布。這與實際觀察情況相符。

圖2 BZ沿著Z軸的分布

圖3 Y=0mm，Z=25mm時，BX、BZ沿X軸分布

圖4、5是BZ、BX在xoy平面的分布情況，當x=r0或x=-r0時，BZ、BX會達到極值，BX在x=0時為0；BZ、BX在Z軸方向上開始衰減比較迅速。

圖4 BZ在xoz平面內的分布

圖5 BX在xoz平面內的分布

圖6、7是當z=25mm時，平行于xoy平面上的BZ、BX分布情況。

圖6 在z=25mm的平面內，BX分布

圖7 在z=25mm的平面內，Bz分布

By的分布和Bx相似，只是場強的方向不同。只要將圖4、5中的X軸換成Y軸就可以得到By的分布圖。

3 結論

①本文給出了圓柱形永磁體的建模方法，推導出求解圓柱形永磁體外空間任意一點的磁感應強度模型。通過和磁偶極子模型對比，本文提出的模型和實測值更吻合，完全滿足工程需求。

②通過MATTALAB的數值積分功能，省去了對積分公式進行離散化的復雜求導過程。因此，避免了積分公式遞推過程中的穩定性問題、精度問題和溢出問題等。

③MATLAB的可視化功能，方便地表示出感應強度的大小和分布規律。

④從宏觀上講，本文推導的模型中場點的值僅取決于所有源點對其的影響，它們之間是通過畢奧-薩伐定理聯系起來的。如果永磁體的參數是常數，那么計算值就不會因為方法本身而產生誤差。因此，本文采用的方法精度更高。

［1］ 張興，方亮，李國麗，等.無線內窺鏡中圓柱永磁體建模與仿真［J］.系統仿真學報，2007（3）：494-497.

［2］ 田錄林，李言，王山石，等.雙筒永磁向心軸承磁力工程化解析算法研究［J］.中國電機工程學報，2007（6）：57-61.

［3］ 林德華，蔡從中，策萬春.方型永磁體表面磁感應強度分布的研究［J］.工科物理，1999（2）：59.

［4］ 張冉，王秀和，楊玉波，等.基于等效剩磁法的永磁電動機轉子偏心磁場解析計算［J］.電工技術學報，2009（5）：7-10.

［5］ 樊明武，顏威利.電磁場積分方程法［M］.北京：機械工業出版社，1988.

［6］ 王淑紅，熊光煌.新型筒型永磁動圈式直線電動機氣隙磁場解析分析［J］.電工技術學報，2007（5）：41-45.

［7］ 鄭文鵬，施進浩，屠關鎮，等.電流域標量處理法在橫向磁場電機三維磁場分析中的應用［J］.電工技術學報，2006（6）：66-70.

［8］ 李永斌，袁瓊，江建中.一種新型聚磁式橫向磁場永磁電機研究［J］.電工技術學報，2003（5）：46-49.

［9］ 章躍進，江建中，屠關鎮.應用數值解析結合法計算旋轉電機磁場［J］.電工技術學報，2004（1）：7-11.

［10］ 曹輝.磁道釘磁場特性試驗研究［D］.武漢：武漢理工大學，2007.

［11］ 孔繁余，陳剛，曹衛東.磁力泵磁性聯軸器的磁場數值計算［J］.機械工程學報，2006（11）：213-217.

［12］ EdwardB.Magrab，ShapourAzarm，Balakumar Balachandran，et al.An Engineer’s Guide to MATLAB with Application from Mechanical,Aerospace,Electrical and Civil Engineering［M］.Beijing：Publishing House of Electronics Indus?try，2006.

［13］ 田錄林，賈嶸，楊國清，等.永磁鐵磁貼合體的磁場及磁力［J］.電工技術學報，2008（6）：7-14.

Modeling and Simulation for Cylinder Permanent Magnetic Field

Zhou Enquan1Zheng Zhongqiao2Zhang Yanhong2Wang Qirui3
（1.Haian Jaorui Talent Robot Technology Co.,Ltd.，Nantong Jiangsu 226000；2.Changzhou Institute of Technology，Changzhou Jiangsu 213032；3.Key Laboratory of Modern Agricultural Equipmentand Technology，Jiangsu University，Zhenjiang Jiangsu 212013）

To find the convenient and simple calculation method of magnetic field for cylinder permanent magnet,based on permanent magnet equivalent electric current model and,regarded magnetization magnetiz?er effect as surface current on magnetizer,this paper presents the law of magnetic field spatial distribu?tion,derives the analytical formulae of magnetism induction density in space for cylinder permanent,and get the answer of the magnetism induction density for cylinder permanent,through the numerical integration function of MATLAB.The result shows that calculated results are in agreement with those actually-mea?sured.It is convenient and simple and can derive precision calculation results by the magnetic field equiva?lent model of permanent magnet and numerical integration method.

permanent magnet；magnetic field；numerical integration；magnetic induction

永磁材料為重要的功能材料，目前應用最多的是圓柱形永磁體、方形永磁體、扇形永磁體和環形永磁體等［1］。隨著永磁體應用領域的擴大，準確得到永磁體磁場分布和大小成為關鍵問題。計算永磁體磁場有解析法［2-4］、電磁積分方程法［5］、有限元法［6-7］、等效磁網絡法［8］和數值解析結合法［9］。目前，國內研究圓柱永磁體磁場分布有磁偶極子［1］和經驗公式法［10］等。磁偶極子法要求場點半徑要遠大于圓柱永磁體的半徑，因此在圓柱體半徑不能忽略時誤差很大。經驗公式法，如曹輝［10］提到的圓柱形永磁體磁感應強度計算公式，只能估算圓柱軸線上的磁感應強度。經驗公式法中的經驗系數要通過大量試驗才能得到［11］，該方法計算精度較低。為了尋求圓柱永磁體磁場簡便而精確的算法，本文將根據永磁體等效電流模型，推導出用于計算圓柱形永磁體外空間任意場點磁感應強度數學模型，并通過MATLAB多重數值積分功能［12］，簡化求解過程，提高求解精度。通過實驗驗證，計算值和實驗值吻合。

TM144

A

1003-5168（2017）11-0139-05

2017-10-09

江蘇省自然科學基金資助項目（BK20151182）。

周恩權（1984-），男，碩士，工程師，研究方向：機械設計理論及高可靠性磁力系統。