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

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

?

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

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

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

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

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

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

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

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

1.1 磁化強(qiáng)度的磁場(chǎng)強(qiáng)度依賴(lài)函數(shù)M(H)

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

=0（1）

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

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

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

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

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

=0(+) （3）

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

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

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

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

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

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

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

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

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

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

^=0r^^（6）

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

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

2 磁化強(qiáng)度的磁感應(yīng)強(qiáng)度依賴(lài)函數(shù)M(B)在永磁體陣列磁場(chǎng)和磁化強(qiáng)度自洽計(jì)算的應(yīng)用

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

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

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

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

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

R=c（7）

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

==-R=-c（8）

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

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

0=0c+（9）

同時(shí)，

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

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

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

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

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

3 結(jié)果和討論

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

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

=0(r+R) （12）

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

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

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

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

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

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

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

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

＜-m0c（14）

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

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

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

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

4 結(jié)論

在強(qiáng)場(chǎng)永磁體陣列磁體材料內(nèi)部，用非自洽方法獲得的磁化強(qiáng)度和磁場(chǎng)分布的分析或計(jì)算結(jié)果，往往與實(shí)際情形相差甚遠(yuǎn)，也很可能導(dǎo)致物理預(yù)測(cè)情形上的誤判。磁性材料磁化強(qiáng)度矢量對(duì)外加磁場(chǎng)強(qiáng)度矢量的響應(yīng)與磁元的幾何形狀相關(guān)，磁化強(qiáng)度的磁場(chǎng)強(qiáng)度依賴(lài)函數(shù)()難以直接用于涉及強(qiáng)場(chǎng)永磁體陣列材料內(nèi)部磁化強(qiáng)度和磁場(chǎng)相互作用的自洽分析。相對(duì)于磁場(chǎng)強(qiáng)度矢量，磁感應(yīng)強(qiáng)度矢量為更加基本且能與材料單一磁化狀態(tài)相關(guān)聯(lián)的物理量。磁化強(qiáng)度的磁感應(yīng)強(qiáng)度依賴(lài)函數(shù)()，適用于強(qiáng)場(chǎng)Halbach永磁體陣列關(guān)于磁場(chǎng)和磁化強(qiáng)度的自洽分析與計(jì)算，研究永磁體陣列產(chǎn)生磁場(chǎng)對(duì)材料內(nèi)部原有磁化狀態(tài)的改變，如在強(qiáng)自場(chǎng)作用下的局部退磁、磁飽和、磁化強(qiáng)度矢量翻轉(zhuǎn)等復(fù)雜的磁化現(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] 羅陽(yáng). 全球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永磁陣列三維磁場(chǎng)端部效應(yīng)建模研究[J]. 電工技術(shù)學(xué)報(bào), 2015, 30(12): 162-170.

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

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

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

[9] 寇寶泉, 曹海川, 李偉力, 等. 新型雙層Halbach永磁陣列的解析分析[J]. 電工技術(shù)學(xué)報(bào), 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] 劉亞丕. 高性能永磁材料磁性測(cè)量若干問(wèn)題探討及測(cè)量技術(shù)進(jìn)展[J]. 電氣技術(shù), 2006, 7(9): 9-16.

[20] 趙凱華, 陳熙謀. 電磁學(xué)[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àn)榇判圆牧蠎?yīng)用等。

國(guó)家自然科學(xué)基金（50577030）