徐小農(nóng) 盧定偉 王智河 陳年林 王坤明

?

磁化強度的磁感應(yīng)強度依賴函數(shù)在Halbach型永磁體陣列研究中的應(yīng)用

徐小農(nóng)1盧定偉1王智河1陳年林1王坤明2

（1. 南京大學物理學院固體微結(jié)構(gòu)實驗室，南京 210093；2. 南京騰馳電子有限公司，南京 211151）

在Halbach型永磁體陣列磁場和磁化強度分布研究中，涉及永磁材料內(nèi)部磁化強度()、磁場強度()和磁感應(yīng)強度()三者的相互關(guān)系函數(shù)()、()和()。本文介紹了這些函數(shù)的典型實驗特征，說明為計算獲得永磁體陣列磁體內(nèi)部自洽的磁場和磁化強度矢量分布，必須使用()函數(shù)。強場永磁體陣列材料內(nèi)部的磁場和磁化強度矢量分布無法直接被進行實驗測量，不恰當?shù)厥褂?)和()函數(shù)易導致不符合實際情況的磁場和磁化強度分布狀態(tài)預測。

Halbach型永磁體陣列；磁感應(yīng)強度；磁化強度與磁場自洽分布

永磁體陣列是實現(xiàn)單側(cè)增強磁場的有效手段[1]，在電氣技術(shù)領(lǐng)域有極為廣泛的應(yīng)用[2-9]。在一定理想條件近似下，永磁體陣列在磁體內(nèi)外磁場分布的解析解[7-15]，在一定近似條件下更接近實際永磁體陣列磁體內(nèi)外磁場分布的數(shù)值解[16-18]，有助于理解永磁體陣列單側(cè)磁場的產(chǎn)生和增強原理，以及永磁體陣列構(gòu)成磁體材料內(nèi)部在強自場作用下的局部退磁甚至磁化強度方向翻轉(zhuǎn)的機制。至今材料內(nèi)部的實際磁化強度分布無法由實驗直接觀測，近年來不斷有針對各種用途永磁體陣列的研究新結(jié)果報道，然而涉及永磁體陣列材料內(nèi)部一些基本的物理現(xiàn)象，如退磁等問題仍然存在一些不完全一致的分析。

在涉及實際永磁體陣列材料的磁化狀態(tài)以及產(chǎn)生磁場空間分布的數(shù)值計算中，需要用到磁體陣列構(gòu)成材料（超高矯頑力硬磁，超高磁能積硬磁）各向異性的磁滯回線。在推導磁體內(nèi)外磁場分布的解析表達式過程中，也需要用到理想化的有關(guān)磁化強度（）、磁感應(yīng)強度（）和磁場強度（）三者的相互關(guān)系。電磁學中磁感應(yīng)強度矢量是代表磁場的基本物理量，而磁場強度矢量是為了描述磁場方便引入的輔助物理量，在電氣工程技術(shù)領(lǐng)域描述磁性材料特征時，常用()磁滯回線來表征永磁材料磁化強度對外加磁場強度的響應(yīng)。本文介紹永磁材料的不同方向的磁特征函數(shù)()，()和()是如何獲得的；說明在永磁體陣列材料內(nèi)部磁感應(yīng)強度和磁化強度分布分析中，使用()函數(shù)存在的局限性以及永磁體陣列磁化狀態(tài)和磁場分布若干基本物理問題存在爭議的原因；強調(diào)要用()函數(shù)才能獲得的材料內(nèi)部自洽的磁感應(yīng)強度和磁化強度分布，從而正確認識實際發(fā)生在永磁體陣列構(gòu)成材料內(nèi)部的、偏離材料初始磁化狀態(tài)分布的復雜磁化現(xiàn)象。

1 永磁體材料的磁化特征函數(shù)

1.1 磁化強度的磁場強度依賴函數(shù)M(H)

硬磁材料具有明顯的各向異性，圖1給出了釹鐵硼永磁材料樣品磁化強度對于磁場強度的實驗磁化曲線()。均勻磁化的橢球狀樣品內(nèi)部的磁場強度應(yīng)當修正為

=0（1）

式中，為橢球樣品的形狀退磁因子；0為外加磁場強度。如圖1中的兩線[0||(0||),0^(0^)]是通過0(00)測量，經(jīng)過式（1）退磁修正后，等價于無退磁效應(yīng)（=0）樣品在平行與垂直磁化主軸兩個方向上的磁滯回線。

圖1 一種釹鐵硼永磁材料樣品在平行與垂直易磁化軸方向m0M||(m0H||), m0M^(m0H^), m0M||(B||), m0M^(B^)磁滯回線

式中，和都很大時（≥10，1），0≈1。這說明了有強退磁形狀的樣品，即使真實的磁化率達到102-4量級，在(0)圖中起始磁化曲線的斜率0也成為了1。

1.2 磁化強度的磁感應(yīng)依賴函數(shù)M(B)

獲得()函數(shù)后，利用在三維空間成立的矢量表達式

=0(+) （3）

再分別作平行與垂直易磁化軸兩個方向上磁化強度和磁感應(yīng)強度的疊加，進而獲得了0()函數(shù)如圖1中的0||(||)和0^(^)回線。

考慮()圖中磁化初始部分=，利用式（3）得到=0(1+)，0()圖中磁化初始部分的斜 率為

當很大時，=1，因此圖1中()回線中幾乎垂直向上的起始磁化線，變成()函數(shù)中約45°的傾角。

注意圖1中易磁化軸方向||(||)的回退曲線與橫軸相交位置處的斜率幾乎為無窮大，也就是說在=-c磁化強度矢量發(fā)生翻轉(zhuǎn)之處，磁化強度存在一定的不確定性，磁化強度并不是磁場強度的單值函數(shù)。易磁化軸方向()的回退曲線與橫軸相交的位置直線斜率為1，因此磁化強度是磁感應(yīng)強度的單值函數(shù)，在數(shù)值計算中使用()函數(shù)，就避免了使用()函數(shù)在-c處遇到的不確定性。

1.3 磁感應(yīng)強度的磁化強度依賴函數(shù)B(H)

圖2是由圖1演變過來，在平行與垂直易磁化軸的兩個方向上的硬磁材料樣品的||(||)，^(^)回線。其中第一、二象限回退曲線在零場附近的線性較好的部分，在易磁化軸方向可以近似表達為

||=0(r||||+R) （5）

式中，r||和0R是擬合常數(shù)，分別等于零場附近||(0||)回退曲線被看成直線時的斜率和在軸上的截距。

圖2 一種釹鐵硼永磁材料樣品樣品在平行與垂直易磁化軸方向B||(H||), B^(H^)回線

垂直于易磁化軸方向的^(^)可以近似表達為

^=0r^^（6）

Halbach最初提出類似于式（5）的()關(guān)系時[1]，是將r看成是擬合常數(shù)，后來有學者定義式（5）中的r為回退磁導率[8-11]。要注意這里的r不符合傳統(tǒng)磁導率的物理定義，因為回退磁化曲線M(H)上各點的磁化強度M，不是由對應(yīng)的磁場強度H激發(fā)并且維持的?；赝?)曲線上隱含的磁化強度，磁體制造工藝上是由脈沖強磁場飽和充磁獲得，磁化強度M能保留依靠了永磁材料內(nèi)部的磁通釘扎，與對應(yīng)的磁場強度H關(guān)聯(lián)度不是太大。

實際的硬磁材料擬合常數(shù)通常r||1r^，在其他與易磁化軸不同交角方向上r也有不同的值。有學者因為由式（5）、式（6）代表的()關(guān)系在平行與垂直易磁化軸的兩個不同方向上，回退曲線r||和r^值很接近，就認為材料是各向同性的[11-13]。實際上式（5）用于描述同一硬磁材料相對于磁化主軸不同方向的()關(guān)系，參數(shù)r即使在相對磁化主軸的不同外場方向很接近，參數(shù)R仍會變化很大。擬合參數(shù)r本來就不符合磁導率的嚴格物理定義，與材料真實的各向異性并無太大的關(guān)聯(lián)，即使r||=r^，也完全不能代表各向同性。圖2所示的釹鐵硼永磁材料樣品()回線，盡管r||與r^接近，對應(yīng)圖1中0(0)起始磁導率和剩余磁化強度的數(shù)值，沿著易化軸方向與沿著垂直易化軸方向之比，都超過了5倍，磁各向異性非常顯著。

2 磁化強度的磁感應(yīng)強度依賴函數(shù)M(B)在永磁體陣列磁場和磁化強度自洽計算的應(yīng)用

在永磁體陣列情形，磁場強度是由磁體提供的自場?？紤]如圖3所示硬磁材料理想磁滯回線0(0)和0()，在＜-c之前的=-c，磁化強度從+R變化到-R，=-c稱為磁極性轉(zhuǎn)變的臨界磁場強度，在臨界磁場上材料的磁化強度值需用0()函數(shù)自洽計算獲得。在圖3所示的0(0)回線中，設(shè)置了參數(shù)R＞c，以展示單一磁塊在自場作用下進入退磁臨界態(tài)=-c的情形。

圖3 永磁材料的理想m0M(m0H)和m0M(B)回線（為方便說明相關(guān)問題，此處設(shè)定參量MR＞Hc=MR/3）

下面以一個均勻磁化的旋轉(zhuǎn)橢球（半軸長度分別為,和,=，磁化主軸平行軸）為例，介紹永磁橢球進入轉(zhuǎn)變臨界態(tài)后，如何獲得自洽的磁化強度和磁感應(yīng)強度一組數(shù)值的思路。這是理解如何獲得永磁體陣列內(nèi)部全自洽的磁化強度和磁感應(yīng)強度矢量分布計算方法的基礎(chǔ)。

圖3所示的0(0)回線給出=0時，退磁效應(yīng)可以被忽略的磁橢球的磁化強度=R。如果磁橢球有退磁效應(yīng)，并且退磁場可以達到或者超過材料的磁矯頑力，即-c，粗略地結(jié)合0(0)回線判斷，就有可能得出磁化強度極性會翻轉(zhuǎn)的結(jié)論。然而，單一永磁橢球在自場作用下磁極性自動翻轉(zhuǎn)，這是違背物理學常識的，實際情形在自場作用下磁化強度會退化但是不會翻轉(zhuǎn)，在0()磁滯回線的回退分支上可以找出自洽的工作點（＞0，0＞0）。

現(xiàn)在假設(shè)剩余磁化強度為與矯頑力的關(guān)系為

R=c（7）

式中，為常數(shù)，＞1時R＞c。如橢球的退磁因子為，那么在外加磁場不存在（0=0）時球內(nèi)的磁場強度全部來源于自場，磁化強度為R的橢球，在球內(nèi)產(chǎn)生的均勻磁場強度可以由式（1）和式（7）得到為

==-R=-c（8）

如果＞1，那么＜-c，退磁必定會發(fā)生。

在圖3所示0()回線給出磁極性轉(zhuǎn)變過程的磁感應(yīng)強度變化范圍是-m0(R+c)＜＜0(R-c)；退磁發(fā)生的條件是=0(+R)=0(-R+R)＜0(R-c)，即R＞c。在0()回線上對應(yīng)磁化強度從+R到-R轉(zhuǎn)變的臨界態(tài)上

0=0c+（9）

同時，

由上列兩式求出一組自洽的(,0)值為

＞1情形在自場作用下，初始飽和磁化的橢球磁化強度值從R自發(fā)退化至=c/。利用式（7）可得c＜=R/＜R。假設(shè)相關(guān)參數(shù)如圖3所示為=3,=1/2,3/2＞1，則式（11）給出退磁因子=1/2的永磁扁橢球（=＜＝在自場作用下，自洽的磁感應(yīng)強度和磁化強度工作點位于圖3第一象限0()中回線上的小圓圈處。此處(,0)=(0c, 20c)=(0R/3, 20R/3)，對應(yīng)0(0)回線上的臨界磁場強度值=/0-=-(1/3)R=-c。

上述具有圖3磁滯回線特征的永磁橢球，在自身磁場作用下磁化強度退化的例子說明，非自洽物理模型下給出的磁場強度值＜-c，并不是一定能導致磁化強度翻轉(zhuǎn)。磁化強度翻轉(zhuǎn)發(fā)生在0()回線的=-m0c處。從圖3可以看出，=-c表明磁化強度可能退化，＜-c則意味著磁化強度徹底翻轉(zhuǎn)達到反向飽和。達到反向磁飽和，需要的磁感應(yīng)強度值＜-m0(Rc)，也就是||＞0(Rc)，這是單個永磁橢球自身無法提供的反向磁場。

以上舉例說明了對于一個永磁體體元，如何在已知0()函數(shù)的基礎(chǔ)上獲得自洽的磁化強度與磁感應(yīng)強度。在永磁體陣列情形，要考慮到所有計算磁元之間的相互作用，磁元磁化強度對磁感應(yīng)強度場響應(yīng)存在的各向異性，磁化強度飽和，磁化強度轉(zhuǎn)向等因素，利用||(||)和^(^)函數(shù)，在采取一定近似的條件下，在永磁體陣列磁體材料內(nèi)部和外部獲得全自洽的磁化強度和磁感應(yīng)強度分布[16-18]。

2.1 乳腺葉狀腫瘤良惡性病理表現(xiàn) 納入研究107例患者的107個病灶其中良性乳腺葉狀腫瘤45個，交界性乳腺葉狀腫瘤33個，惡性乳腺葉狀腫瘤29個。其中，良性病灶出現(xiàn)膨脹性形狀，且間質(zhì)細胞較為疏松的排列，具有明顯的核分裂象(每10個高倍視野下0～4個核分裂象)，但并無出血壞死的病理表象；交界性乳腺葉狀腫瘤呈現(xiàn)明顯的浸潤性生長，且間質(zhì)細胞生長過度，具有明顯的核分裂象(每10個高倍視野下5～9個核分裂象)，可見少量出血壞死現(xiàn)象；惡性乳腺葉狀腫瘤間質(zhì)細胞排列緊密，但胞質(zhì)內(nèi)明確可見空泡，且導管稀少，每10個高倍視野下10個核分裂象，可見大片出血和壞死。

3 結(jié)果和討論

3.1 永磁體陣列材料內(nèi)部磁感應(yīng)強度，磁場強度和磁化強度三者的基本關(guān)系

一些有關(guān)永磁體陣列磁場分布的解析解的研究論文[10-11]，將在易磁化軸一維方向成立的式（5）推廣到三維矢量空間的等式

=0(r+R) （12）

式中，R代表剩余磁化強度矢量的三維空間分布，而后由已知的R空間分布，通過求出矢量磁勢或標量磁勢，導出磁感應(yīng)強度矢量分布或者磁化強度的矢量分布。我們認為將易磁化軸方向上磁場強度和磁感應(yīng)強度兩者關(guān)系的式（5），推廣成為二維或三維矢量形式的式（12），其合理性存疑。式（5）代表在易磁化軸方向上，磁場強度與磁感應(yīng)強度兩個變量，在一定磁化強度變化范圍內(nèi)的依賴關(guān)系，R的地位只是一個擬合常數(shù)，在數(shù)值上等于剩余磁化強度R=(=0)。永磁體材料內(nèi)部磁感應(yīng)強度，磁場強度和磁化強度三者基本矢量關(guān)系只能是式（3）。

式（5）中隱含了一個對場敏感的磁化強度模型，由式（5）和式（3）消去磁感應(yīng)強度得到磁化強度，磁場強度和剩余磁化強度三者矢量關(guān)系為

||=(r||-1)||+R（13）

對照圖1可見，在易磁化軸方向式（13）在＞-c大致成立，在＜-c實驗||(||)關(guān)系完全偏離式（8）表達的線性。文獻[11]在式（12）的基礎(chǔ)之上求解矢量磁勢的方法，給出了理想Halbach圓柱型陣列材料內(nèi)外磁場分布r≥1的解析解，并且否定了文獻[10]給出的磁場解析解。利用復變函數(shù)直接積分的方法，導出的理想Halbach圓柱型陣列磁體內(nèi)外磁感應(yīng)強度分布解析解[15]，與文獻[11]中的磁場解析表達式在（r=1）時完全一致，這也進一步證實了文獻[11]作者的判斷是正確的，文獻[10]給出的磁場分布解析解有誤。

文獻[10-11]采用物理模型本身是不自洽或者是有矛盾的。在推導磁化強度對磁場敏感的情形下永磁體陣列磁場分布解析解的過程中，作者一方面引入常數(shù)r＞1表明磁化強度的幅值是對外加磁場敏感的，另一方面在推演過程中使用了||=R的關(guān)系，用式（12）中=0時的剩余磁化強度矢量分布R導出同式中磁場強度分布，并未進一步考慮式（13）所示磁化強度隨外場的改變，也就是在實際推演過程中，最終又忽略了磁化強度對外加磁場敏感性。

文獻[11]在r＞1情形，即磁化強度對磁場強度敏感的假設(shè)下，以式（12）為基本函數(shù)導出了的雙極柱形永磁體陣列空腔內(nèi)部磁場的解析表達式，始終代表一個與r值有關(guān)，大小和方向都恒定的均勻磁場。而實驗測量或者依據(jù)實驗磁化回線()自洽計算出的雙極Halbach永磁體陣列空腔內(nèi)的磁感應(yīng)強度分布，是有一定變化規(guī)律的非均勻場[16]，所以在式（12）基礎(chǔ)之上導出的空腔內(nèi)磁感應(yīng)強度的解析表達式，在定性上并不符合Halbach永磁體陣列磁場分布的實驗變化趨勢。磁化強度對磁場敏感的永磁材料（r＞1，或者|1常數(shù)R）的永磁體陣列內(nèi)外磁感應(yīng)強度的分布解析表達式，至今依然是一個尚未解決的挑戰(zhàn)性問題。

3.2 在永磁體陣列磁場的分析和計算中該使用M(H)函數(shù)還是M(B)函數(shù)？

文獻[13]討論了磁塊局部區(qū)域的磁極性在永磁體陣列自身強自場作用下翻轉(zhuǎn)的可能性。該文依據(jù)永磁材料的磁化曲線M(H)（如圖1和圖3所示），在反向磁場強度H小于矯頑力場c時（H＜-c）磁化強度轉(zhuǎn)向，得出Halbach中空柱形永磁體陣列南北外側(cè)兩極處的局部材料磁極性，在材料磁參量值c＜R條件下會翻轉(zhuǎn)的結(jié)論。這里需要指出，得出這個翻轉(zhuǎn)結(jié)論所依據(jù)的數(shù)值結(jié)果，并不符合永磁體陣列磁體內(nèi)部磁化強度矢量和磁場強度矢量的自洽分布。

＜-m0c（14）

分別用復變函數(shù)積分導出的理想中空圓柱型永磁體陣列磁體內(nèi)外磁感應(yīng)強度的解析分布表達式，用求解泊松方程方法導出理想中空球形永磁體陣列磁體內(nèi)外磁感應(yīng)強度的解析分布表達式，進而依據(jù)理想永磁體陣列磁感應(yīng)強度在材料內(nèi)部解析分布，建立起解析的退磁函數(shù)[14-15]。用實驗()關(guān)系自洽數(shù)值計算獲得的實際永磁體陣列材料內(nèi)部的磁化狀態(tài)分布[16-18]，都不支持文獻[13]提出的Halbach中空雙極永磁體陣列外側(cè)南北兩極處的磁塊局部磁化強度方向、在材料磁參量值c＜R條件下會翻轉(zhuǎn)的結(jié)論。

3.3 軟鐵在永磁體陣列產(chǎn)生的磁場中是容易被飽和磁化的嗎？

圖1中的0||(0||)是包含了被測量硬磁材料樣品自身對磁場強度的貢獻以后的磁化曲線0||(0||)。注意起始磁化曲線起始部分幾乎是垂直向上的，在0||值達到約0.2T時就獲得磁飽和，這說明無退磁效應(yīng)的內(nèi)置磁化主軸的塊狀硬磁，在順著易磁化軸的方向很容易被飽和磁化；而軟磁材料的起始磁化曲線往往比硬磁材料更接近垂直，通常起始磁化率＞102。然而，若認為依靠永磁體或者永磁體陣列產(chǎn)生的靜磁場，大塊的軟磁輕易會被飽和磁化，則是不符合實際情況的。式（2）表明即使是極其容易被磁化的軟磁材料，若有比較強的樣品形狀退磁效應(yīng)（10），則在(0)圖中的起始磁化曲線的傾角也會接近45°而不是90°。有強退磁效應(yīng)的樣品，只有在外磁場00基本達到飽和磁化強度值0R的時候，才能被飽和磁化[17-18]。

永磁體陣列任一體元的磁化強度，應(yīng)該看作是當?shù)卮鸥袘?yīng)強度的函數(shù)，需計及磁元自身的貢獻。磁感應(yīng)強度矢量和磁化強度矢量自洽數(shù)值計算的結(jié)果表明：在強場永磁體陣列內(nèi)部磁化強度達到飽和的區(qū)域，適度地用軟鐵替代部分硬磁，圓柱體陣列空腔內(nèi)部的總體磁通量和磁場均勻度都可能得到改善；但是替代部分硬磁材料的大塊軟磁，并非輕易達到飽和磁化狀態(tài)。軟鐵難以飽和的原因可以從圖1所示0()回線的特征以及式（4）看出，軟磁材料的起始磁化線的斜率≤1，而軟磁材料的飽和磁化強度，通常會遠高于硬磁材料的飽和磁化強度。軟磁材料替代硬磁材料后主要起到匯聚磁通量到空腔的作用，流經(jīng)空腔磁隙內(nèi)的磁通量會增加的同時，永磁體陣列整體磁路中的總磁通量反而會減少[18]。

4 結(jié)論

在強場永磁體陣列磁體材料內(nèi)部，用非自洽方法獲得的磁化強度和磁場分布的分析或計算結(jié)果，往往與實際情形相差甚遠，也很可能導致物理預測情形上的誤判。磁性材料磁化強度矢量對外加磁場強度矢量的響應(yīng)與磁元的幾何形狀相關(guān)，磁化強度的磁場強度依賴函數(shù)()難以直接用于涉及強場永磁體陣列材料內(nèi)部磁化強度和磁場相互作用的自洽分析。相對于磁場強度矢量，磁感應(yīng)強度矢量為更加基本且能與材料單一磁化狀態(tài)相關(guān)聯(lián)的物理量。磁化強度的磁感應(yīng)強度依賴函數(shù)()，適用于強場Halbach永磁體陣列關(guān)于磁場和磁化強度的自洽分析與計算，研究永磁體陣列產(chǎn)生磁場對材料內(nèi)部原有磁化狀態(tài)的改變，如在強自場作用下的局部退磁、磁飽和、磁化強度矢量翻轉(zhuǎn)等復雜的磁化現(xiàn)象。

[1] K. Halbach. Design of permanent multipole mangnets with oriented rare earth cobalt material[J]. Nuclear Instruments and Methods, 1980, 169: 1-10.

[2] J. M. D. Coey. Permanent magnet applications[J]. Journal of Magnetism and Magnetic Materials, 2002, 248: 441-456.

[3] 羅陽. 全球NdFeB磁體產(chǎn)業(yè)的變化與發(fā)展[J]. 電氣技術(shù), 2008, 9(6): 5-11.

[4] 張一鳴, 喬德治, 高俊俠. Halbach陣列永磁體的研究現(xiàn)狀與應(yīng)用[J]. 分析儀器, 2010(2): 5-10.

[5] 宋玉晶, 張鳴, 朱煜, 等. 基于偽周期的Halbach永磁陣列三維磁場端部效應(yīng)建模研究[J]. 電工技術(shù)學報, 2015, 30(12): 162-170.

[6] 羅辭勇, 李竹田, 沈啟平. 異形永磁體圓柱型直線電機的優(yōu)化設(shè)計[J]. 電工技術(shù)學報, 2017, 32(17): 127-134.

[7] 羅玲, 薛利昆, 吳先宇, 等. Halbach永磁陣列無刷直流電機轉(zhuǎn)矩的解析計算和分析[J]. 電工技術(shù)學報, 2017, 32(16): 124-135.

[8] 宋玉晶, 張鳴, 朱煜. Halbach永磁陣列磁場解析求解及推力建模[J]. 電工技術(shù)學報, 2014, 29(11): 61-67.

[9] 寇寶泉, 曹海川, 李偉力, 等. 新型雙層Halbach永磁陣列的解析分析[J]. 電工技術(shù)學報, 2015, 30(10): 68-76.

[10] Xia Z P, Zhu Z Q, Howe D. Analytical magnetic field analysis of Halbach magnetized permanent magnet machines[J]. IEEE Transactions on Magnetics, 2004, 40: 1864-1872.

[11] Bj?rk R, Smith A, Bahl C R H. Analysis of the magnetic field, force, and torque for two-dimensional Halbach cylinders[J]. Journal of Magnetism and Magnetic Materials, 2010, 322: 133-141.

[12] Bj?rk R, Smith A, Bahl C R H. The efficiency and the demagnetization field of a general Halbach cylinder[J]. Journal of Magnetism and Magnetic Materials, 2015, 384: 128-132.

[13] Insinga A R, Bahl C R H, Bj?rk R, et al. Permanence of Halbach magnet arrays with finite coercivity[J]. Journal of Magnetism and Magnetic Materials, 2016, 407: 369-376.

[14] Xu Xiaonong, Lu Dingwei, Xu Xibin, et al. Analytic solution of magnetic induction distribution of ideal hollow spherical field sources[J]. Journal of Mag- netism and Magnetic Materials, 2017, 443: 397-403.

[15] Xu Xiaonong, Lu Dingwei, Xu Xibin, et al. Analytic solution of field distribution and demagnetization function of ideal hollow cylindrical field source[J]. AIP Advances, 2017, 7: 095203.

[16] Xu X N, Lu D W, Yuan G Q, et al. Self-modification of the remanence distribution of a hollow cylindrical permanent magnet array[J]. Journal of Applied Physics, 2006, 100: 043906.

[17] Xu X N, Lu D W, Xu X B, et al. Simulations of remanence and field distributions of permanent magnet array with rectangular cavity[J]. Journal of Applied Physics, 2008, 104: 043919.

[18] Xu Xiaonong, Lu Dingwei, Xu Xibin, et al. Optimizing the field distribution of a Halbach type permanent magnet cylinder using the soft iron and superhard magnet[J]. Journal of Magnetism and Magnetic Materials, 2018, 445: 77-83.

[19] 劉亞丕. 高性能永磁材料磁性測量若干問題探討及測量技術(shù)進展[J]. 電氣技術(shù), 2006, 7(9): 9-16.

[20] 趙凱華, 陳熙謀. 電磁學[M]. 2版. 北京: 高等教育出版社, 1985.

Applications of the magnetic induction strength dependent magnetization function in the study of Halbach type permanent magnet array

Xu Xiaonong1Lu Dingwei1Wang Zhihe1Chen Nianlin1Wang Kunming2

(1. Physics College of Nanjing University, National Laboratory of Solid State Microstructures, Nanjing 210093;2. Nanjing Tengchi Electronics Corporation Limited, Nanjing 211151)

In the study of magnetic field and magnetization distribution generated by Halbach permanent magnet array, the relationship functions(),() and() among magnetization (), magnetic field strength () and magnetic induction () are involved. Here the experimental characteristics of these functions are introduced, and it is illustrated that() function should be used in the calculation for obtaining the self consistent magnetization and field distributions inside the permanent magnet array. At the material interior of permanent magnet array, the sophisticated vectorial distributions of field and magnetization may not be experimentally measured, inappropriate use of() and() functions possibly results in the untrue predication of magnetic field and magnetization distributions under the intense self field.

Halbach permanent magnet array; magnetic induction; self-consistent distribution of magnetic field and magnetization

2018-04-23

徐小農(nóng)（1962-），男，博士，教授，研究方向為磁性材料應(yīng)用等。

國家自然科學基金（50577030）