• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    納米硬/軟磁復(fù)合永磁薄膜磁滯回線的微磁學(xué)理論綜述

    2016-05-06 03:05:03趙國平萬秀琳
    關(guān)鍵詞:磁滯回線

    趙國平, 萬秀琳

    (四川師范大學(xué) 物理與電子工程學(xué)院固體物理研究所, 四川 成都 610066)

    ?

    納米硬/軟磁復(fù)合永磁薄膜磁滯回線的微磁學(xué)理論綜述

    趙國平,萬秀琳

    (四川師范大學(xué) 物理與電子工程學(xué)院固體物理研究所, 四川 成都 610066)

    摘要:納米硬/軟磁交換彈簧復(fù)合材料是一種新型的高性能永磁體材料,在基礎(chǔ)研究和實(shí)際應(yīng)用2個(gè)方面都得到廣泛的重視.它有著豐富的磁學(xué)性質(zhì),特別是在磁能積方面潛力巨大,理論磁能積密度高達(dá)1 MJ/m3.微磁學(xué)是量子磁學(xué)和宏觀磁性的橋梁,是研究材料磁滯回線和宏觀磁性的重要理論,它能清晰地闡述材料的磁化反轉(zhuǎn)機(jī)制,預(yù)測(cè)磁滯回線和磁性能.以薄膜結(jié)構(gòu)為例,綜述3種重要的計(jì)算復(fù)合磁體磁滯回線的微磁學(xué)方法:一維全解析模型;一維半解析模型;三維數(shù)值計(jì)算模型,并分別討論3種方法的優(yōu)點(diǎn)與不足.將不同方法得到的計(jì)算結(jié)果進(jìn)行比較,重點(diǎn)分析硬/軟磁相厚度和易磁化軸取向?qū)?fù)合磁體宏觀磁滯回線以及矯頑力機(jī)制的影響,并與實(shí)驗(yàn)進(jìn)行比較,討論理論和實(shí)驗(yàn)產(chǎn)生差別的原因,提出一些提高材料磁性能特別是矯頑力和磁能積的的方法.

    關(guān)鍵詞:微磁學(xué)理論; 磁滯回線; 磁能積

    納米尺度內(nèi)的交換耦合復(fù)合材料[1-2]結(jié)合了軟磁相的高剩磁和硬磁相的高矯頑力,具有很高的磁能積.該材料一經(jīng)提出,就受到人們的廣泛關(guān)注[3-15],是目前磁學(xué)以及磁性材料領(lǐng)域中一個(gè)重要的研究方向[15-29].在制造磁電阻[30-40]、永磁材料[3-4,41]和巨磁致伸縮材料[42-44]等方面,交換彈簧都有著重要應(yīng)用.

    R. Skomski等[4]預(yù)言硬軟磁復(fù)合材料最大磁能積密度(BH)max可以達(dá)到1 MJ/m3,是目前最好的單相NdFeB永磁體理論值的2倍左右,展現(xiàn)出誘人的前景.許多磁學(xué)工作者試圖在實(shí)驗(yàn)中獲得這樣巨大的磁能積,但是實(shí)驗(yàn)所得到的結(jié)果僅有理論值的20%~40%[5,45-48],迄今為止,最好的實(shí)驗(yàn)結(jié)果也只能達(dá)到486 kJ/m3[23].由于磁性能的實(shí)驗(yàn)結(jié)果與理論值之間存在很大差距,研究者采用各種理論模型及數(shù)值模擬對(duì)硬磁/軟磁交換彈性薄膜進(jìn)行研究,試圖找到產(chǎn)生差距的原因,得到了一些定性甚至定量結(jié)果.

    T. Leineweber等[12,49]計(jì)算了Nd2Fe14B/α-Fe磁性三層膜的剩磁隨軟磁相厚度的變化曲線,以及成核場(chǎng)、最大磁能積隨軟磁相厚度變化的曲線,并將成核場(chǎng)隨軟磁相厚度的變化趨勢(shì)分為3個(gè)區(qū)間,隨著軟磁相厚度的增加,成核場(chǎng)先后經(jīng)歷平臺(tái)、急速降低和緩慢降低3個(gè)區(qū)間.

    G. P. Zhao等[9-11,29]經(jīng)過詳細(xì)分析發(fā)現(xiàn)以上的平臺(tái)區(qū)間并不存在,并進(jìn)一步指出復(fù)合磁體的磁化反轉(zhuǎn)一般要經(jīng)歷反向磁疇的成核、反向磁疇的擴(kuò)展和磁化矢量的不可逆翻轉(zhuǎn).計(jì)算表明:隨著軟磁層厚度的增加,矯頑力的機(jī)制從成核變?yōu)獒斣?并在此基礎(chǔ)上提出了自釘扎的矯頑力機(jī)制和成核釘扎二重性的反磁化機(jī)制[29],這一理論得到了S. Sawatzki等[50]實(shí)驗(yàn)的證明.文獻(xiàn)[51-54]還分析了Nd2Fe14B/α-Fe復(fù)合磁性材料實(shí)現(xiàn)最大磁能積的薄膜厚度,認(rèn)為在硬軟磁相厚度均為10 nm左右可以實(shí)現(xiàn)最大磁能積約為800 kJ/m3,并考慮了易軸取向?qū)?fù)合磁多層膜磁性能的破壞作用.

    J. S. Jiang等[13-14,55]通過實(shí)驗(yàn)及理論分析得到了在一定程度上提高磁能積的方法,比如利用制作工藝改變界面的結(jié)構(gòu),致使界面的耦合強(qiáng)度增大,從而提高材料的磁能積.G. Asti等[7-8,56-60]計(jì)算了磁晶易軸平行和垂直于膜面的交換耦合納米復(fù)合材料的磁相圖,將磁體分為剛性復(fù)合磁體(RM)、交換彈簧磁體(ES)和退耦合磁體(DM),得到了不同材料的軟磁相和硬磁相的最佳比例,推導(dǎo)出易軸垂直于膜面情況下的釘扎場(chǎng)公式,并分析了界面交換耦合強(qiáng)度對(duì)磁滯回線的影響.

    近年來,實(shí)驗(yàn)上也對(duì)交換彈簧進(jìn)行了很多研究[45,61-71],磁能積得到很大提高,其中,關(guān)于硬磁軟磁多層膜的文獻(xiàn)很多.研究發(fā)現(xiàn),相比于塊材,薄膜材料的晶體易軸和材料的均勻性都比較容易控制,可以得到較高的磁能積.此外,現(xiàn)代技術(shù)可以較容易地掌控薄膜厚度,其實(shí)驗(yàn)結(jié)果可以很方便地與理論進(jìn)行對(duì)比,據(jù)此可以分析薄膜厚度對(duì)磁性能特別是矯頑力機(jī)制的影響.2011年,S. Sawatzki[50]等通過在軟磁相Fe磁層中加入少量Co來提高SmCo5/Fe三層膜的磁性能,其最大磁能積可達(dá)到39 MGOe.2012年,W. B. Cui等[23]通過在界面加入很薄的非磁性Ta層對(duì)硬磁相和軟磁相的元素進(jìn)行隔離,合成的Nd2Fe14B/FeCo的磁能積達(dá)到了61 MGOe.A. M. Belemuk等[72-73]還研究了溫度對(duì)磁性能的影響.

    微磁學(xué)最重要的任務(wù)之一是計(jì)算磁滯回線,從磁滯回線上可以比較不同磁性材料的磁性能特別是矯頑力.矯頑力是表征材料在磁化以后保持磁化狀態(tài)能力的大小,是考察永磁材料的標(biāo)準(zhǔn)之一,也是劃分軟磁和永磁材料的依據(jù)[24].目前復(fù)合磁體微磁學(xué)計(jì)算理論和實(shí)驗(yàn)的主要差距就是在矯頑力方面,理論計(jì)算的矯頑力比實(shí)驗(yàn)大很多,導(dǎo)致理論計(jì)算的磁能積偏大.這一差距幾乎在所有材料中都存在,基于此W. F. Brown[74-75]提出了矯頑力悖論.這一悖論產(chǎn)生的根本原因是理論計(jì)算一般都忽略了某些細(xì)節(jié),比如說材料的缺陷、微結(jié)構(gòu)和溫度等,這些細(xì)節(jié)往往會(huì)對(duì)矯頑力造成致命的破壞.

    本文通過比較的方法主要討論3種微磁學(xué)方法及其優(yōu)缺點(diǎn):一維全解析模型、一維半解析模型以及三維數(shù)值模擬計(jì)算3種方法.通過細(xì)致的比較,對(duì)現(xiàn)有理論取得的成績(jī)和存在的問題作較為系統(tǒng)的綜述,并討論理論和實(shí)驗(yàn)結(jié)果存在差別的原因.

    一維全解析模型[18,76-77]的優(yōu)點(diǎn)在于公式和推導(dǎo)過程都很簡(jiǎn)單,具有清晰的物理意義,特別是可以清楚地給出材料的矯頑力機(jī)制,對(duì)磁化物理過程有合理的解釋.但是這種模型往往只考慮了磁晶各向異性能和外場(chǎng)能,忽略了磁矩之間的相互作用,只能粗略地描述磁化反轉(zhuǎn)過程,得到的矯頑力和磁能積與實(shí)驗(yàn)結(jié)果相差很大.

    與以上簡(jiǎn)單模型相反,許多數(shù)值模擬軟件可以實(shí)現(xiàn)較為復(fù)雜的微磁學(xué)模擬,如利用有限元法(FEM)的Nmag、Magpar和Ansys等軟件,以及利用有限差分法(FDM)的OOMMF和Matlab等軟件.有限元法適用于非立方體結(jié)構(gòu)磁性材料的模擬計(jì)算,H. Kronmüller等[78-79]利用這種方法模擬了Nd2Fe14B/α-Fe和Pr2Fe14B/Fe納米復(fù)合材料的磁滯回線;M. Yue等[80-81]也運(yùn)用有限元方法計(jì)算了MnBi/Fe復(fù)合磁體的磁滯回線.G. P. Zhao等[82-85]則運(yùn)用OOMMF軟件計(jì)算了Nd2Fe14B/Fe2Co、Sm-Co/Fe以及Nd2Fe14B/Fe等材料的磁滯回線和磁矩自旋分布.通過求解Landau-Lifshitz-Gilbert方程的三維數(shù)值模擬方法可以得到磁矩分布隨外場(chǎng)變化的三維動(dòng)態(tài)過程,并進(jìn)而求得磁滯回線.本文將重點(diǎn)討論微磁學(xué)模擬中最常使用的軟件OOMMF及其在復(fù)合材料微磁學(xué)計(jì)算中的應(yīng)用.除了有限元和有限差分方法,還有其他方法來研究磁學(xué)性質(zhì),比如Monte Carlo方法,利用概率和隨機(jī)性來計(jì)算磁矩的分布;原子鏈模型則是一維數(shù)值計(jì)算模型,其計(jì)算方法與以下討論的一維半解析模型中的數(shù)值計(jì)算部分類似.數(shù)值計(jì)算可以模擬多種形狀的磁性材料,但是也有一定局限性.由于磁疇壁里面的磁矩分布對(duì)磁滯回線有著重要影響,精確模擬磁疇壁是微磁學(xué)計(jì)算的一個(gè)重要任務(wù),其網(wǎng)格大小要遠(yuǎn)小于疇壁寬度.稀土永磁體的疇壁都比較窄,一般在2~5 nm,因此劃分的網(wǎng)格較多,計(jì)算量較大,需要高速度、大容量的計(jì)算機(jī).

    一維半解析模型[7-8,86-88]是介于三維數(shù)值計(jì)算和全解析模型之間的一種方法,它在簡(jiǎn)單和精確之間維持了一定的平衡,是硬磁/軟磁復(fù)合磁體計(jì)算中用得最多的一種方法.它考慮了相鄰磁矩之間的交換相互作用,通過公式推導(dǎo),可以反映出磁化反轉(zhuǎn)的物理過程.它計(jì)算量較小,適用于幾何形狀較為規(guī)則的體系如硬磁/軟磁多層膜,能夠比較準(zhǔn)確地描述出不同外場(chǎng)下的磁化反轉(zhuǎn)情況.與一維全解析模型相比,這種方法的解析推導(dǎo)要復(fù)雜得多,推導(dǎo)過程中有時(shí)需要進(jìn)行簡(jiǎn)化處理,導(dǎo)致計(jì)算結(jié)果的精確程度下降.

    以上理論研究可以從不同側(cè)面很好地揭示硬磁/軟磁多層膜的磁化反轉(zhuǎn)機(jī)制.然而大多數(shù)的理論工作都具有一定的局限性,比如數(shù)值計(jì)算模型往往過于復(fù)雜,而且計(jì)算結(jié)果只是針對(duì)某種材料而言的,不具有普遍性.目前文獻(xiàn)中對(duì)不同計(jì)算結(jié)果和計(jì)算模型進(jìn)行比較和分析的工作相對(duì)較為缺乏.本文力圖在這方面作一個(gè)綜述,起到拋磚引玉的效果.下面從最簡(jiǎn)單的模型,即著名的S-W一致轉(zhuǎn)動(dòng)模型,討論微磁學(xué)的計(jì)算方法,并分析磁滯回線產(chǎn)生的物理根源.

    1一維全解析模型

    1948年,E. C. Stoner和E. P. Wohlfarth[18]提出了一致轉(zhuǎn)動(dòng)模型(S-W模型).S-W模型假設(shè)磁性材料由微小晶粒構(gòu)成并且每個(gè)晶粒都具有單易化軸,不考慮晶粒之間的相互作用,內(nèi)部的磁化反轉(zhuǎn)是通過一致轉(zhuǎn)動(dòng)實(shí)現(xiàn)的.該模型適用于單疇粒子,其自由能E只由磁晶各向異性能和塞曼能構(gòu)成.對(duì)于較大的粒子或多晶粒子的磁化反轉(zhuǎn)主要是通過磁疇壁的移動(dòng)實(shí)現(xiàn)的,即疇壁移動(dòng)模型[7,20-23],后者所得到的理論值更加接近實(shí)驗(yàn)值,但是更為復(fù)雜.這里主要討論S-W模型.

    1.1S-W模型根據(jù)S-W模型,能量密度公式可以寫為

    (1)

    其中,K和Ms分別表示磁晶各向異性常數(shù)和飽和自發(fā)磁化強(qiáng)度,β表示磁晶易軸與外加磁場(chǎng)的夾角,通過對(duì)能量求極小值,即?E/?θ=0和?E2/?θ2>0就可以得到磁滯回線和矯頑力.

    圖1是根據(jù)S-W模型計(jì)算得到單軸磁性材料在不同易軸偏角下的磁滯回線.當(dāng)β=0°時(shí),對(duì)應(yīng)的磁化強(qiáng)度只有2個(gè)穩(wěn)定的態(tài),即當(dāng)外場(chǎng)等于臨界場(chǎng)Hcrit的時(shí)候,磁化矢量就從M=+Ms直接跳轉(zhuǎn)到M=-Ms,發(fā)生不可逆反轉(zhuǎn),得到的磁滯回線是方形的,此時(shí)矯頑力Hc=Hcrit=Hk=2K/Ms.在β=30°時(shí),磁滯回線的方形度變差,剩磁和矯頑力都降低.在β=90°時(shí),可以得到一條過原點(diǎn)的直線,沒有形成磁滯回線.在β小于45°時(shí),矯頑力大小等于臨界場(chǎng),矯矯力機(jī)制為釘扎,當(dāng)β大于45°之后,矯頑力不再等于臨界場(chǎng).在外場(chǎng)小于臨界場(chǎng)之前,磁化矢量都是可逆的,當(dāng)外場(chǎng)大于臨界場(chǎng),磁矩將發(fā)生不可逆反轉(zhuǎn).根據(jù)S-W模型,取向的磁性材料矯頑力等于其各向異性場(chǎng)Hk,而各向同性材料(如甩帶)的矯頑力為0.48Hk.對(duì)于磁王Nd2Fe14B來說,分別對(duì)應(yīng)于6.7 T和3.3 T,而實(shí)驗(yàn)上燒結(jié)Nd2Fe14B和各向同性甩帶Nd2Fe14B的矯頑力值一般在1~2 T,理論值比實(shí)驗(yàn)值大很多,跟W. F. Brown[74-75]提出的矯頑力悖論吻合.一般來說,S-W模型只能用于討論單晶單疇粒子的矯頑力和磁滯回線,它可以使物理問題大大簡(jiǎn)化,并清楚地反映出矯頑力機(jī)制,將Hc與Hk緊密地聯(lián)系起來.后來的研究發(fā)現(xiàn),雖然S-W模型給出的矯頑力數(shù)值偏大,但是永磁材料的矯頑力和各向異性場(chǎng)Hk確實(shí)存在一個(gè)線性關(guān)系[28-29,89-91].

    1.2一致轉(zhuǎn)動(dòng)模型的磁滯回線S-W模型一般適用于單疇粒子.實(shí)際上,熱退磁態(tài)的單疇粒子在剩磁態(tài)往往不再是單疇粒子,所以嚴(yán)格地說一致轉(zhuǎn)動(dòng)模型在多數(shù)情況下并不適用,即使對(duì)于所謂的(在熱退磁狀態(tài)下的)單疇粒子也不適用[28-29].但是如上所述,一致轉(zhuǎn)動(dòng)模型作為一個(gè)近似模型具有清晰的物理意義,對(duì)實(shí)際的磁化反轉(zhuǎn)過程也具有借鑒意義.基于此,G. P. Zhao等[76-77,92]將一致轉(zhuǎn)動(dòng)模型推廣到硬/軟磁復(fù)合材料并將一致轉(zhuǎn)動(dòng)模型的結(jié)果與更為精確的非一致轉(zhuǎn)動(dòng)模型進(jìn)行比較,發(fā)現(xiàn)一致轉(zhuǎn)動(dòng)模型的成核場(chǎng)較高,對(duì)應(yīng)一個(gè)不容易出現(xiàn)的成核模式.研究發(fā)現(xiàn),當(dāng)軟磁相厚度較小時(shí),一致轉(zhuǎn)動(dòng)模型給出的磁滯回線和非一致轉(zhuǎn)動(dòng)模型比較接近,可以粗略地反映磁化反轉(zhuǎn)機(jī)制.

    下面比較2種極端情況:一種是硬軟磁之間完全退耦合,磁矩只隨外場(chǎng)變化;第二種是硬軟磁之間完全耦合,整個(gè)系統(tǒng)的磁矩隨外場(chǎng)一致轉(zhuǎn)動(dòng),如圖2所示.

    對(duì)硬/軟磁復(fù)合材料,根據(jù)一致轉(zhuǎn)動(dòng)模型,總能量密度[24]可以寫成

    (2)

    其中,右邊的第一項(xiàng)為硬磁相的總能量密度,第二項(xiàng)表示軟磁相的總能量密度,A表示交換作用常數(shù),f是體積所占比例,上標(biāo)h和s分別表示硬磁層和軟磁層,θ是磁化矢量和外場(chǎng)方向的夾角.

    由于多層膜體系具有周期性和對(duì)稱性,因此多層膜可以看成雙層膜的疊加,這里考慮硬/軟磁雙層膜體系,其性能可以推廣到多層膜和三層膜體系.根據(jù)一致轉(zhuǎn)動(dòng)模型可以知道,在易軸和外場(chǎng)平行時(shí)矯頑力和臨界場(chǎng)是相等的.由方程(2)可以得到矯頑力和臨界場(chǎng)公式為

    (3)

    (4)

    假設(shè)硬磁相和軟磁相的磁晶各向異性差10倍,即λ等于0.1,代入(4)式可以計(jì)算出矯頑力并畫出磁滯回線如圖3所示.圖3(a)表示外場(chǎng)與易軸平行時(shí)根據(jù)一致轉(zhuǎn)動(dòng)模型得到的不同厚度下的硬/軟磁雙層膜的磁滯回線,圖3(b)表示外場(chǎng)與軟磁相易軸垂直時(shí)(同時(shí)與硬磁相易軸平行)根據(jù)一致轉(zhuǎn)動(dòng)模型得到的不同厚度下的硬/軟磁雙層膜的磁滯回線.2種情況下,磁滯回線都是方形的并且主要的矯頑機(jī)制是成核.矯頑力隨著軟磁相厚度的增加而減小,垂直時(shí)的矯頑力比平行情況下的矯頑力要小,也就是說磁矩在易軸平行取向時(shí)比垂直取向時(shí)更容易反轉(zhuǎn).

    (5)

    圖4表示的是根據(jù)完全退耦合模型模擬得到的雙層膜磁滯回線,其磁化狀態(tài)有3種,即硬磁相磁矩和軟磁相磁矩平行的狀態(tài)、反平行狀態(tài)和反向平行的狀態(tài).對(duì)應(yīng)的退磁過程可以分為兩部分:首先是軟磁相磁矩的反轉(zhuǎn),對(duì)應(yīng)磁化反轉(zhuǎn)的成核過程;然后是硬磁相的磁矩反轉(zhuǎn),對(duì)應(yīng)磁化反轉(zhuǎn)的釘扎過程.這2個(gè)過程是完全獨(dú)立的,在軟磁相磁矩完全反轉(zhuǎn)之后,硬磁相磁矩才開始偏轉(zhuǎn).軟磁相與硬磁相的反轉(zhuǎn)場(chǎng)分別等于各自的磁晶各向異性場(chǎng).

    利用以上簡(jiǎn)單模型,可以很容易推導(dǎo)出解析公式,得到矯頑力明確的物理意義及其隨厚度的變化,在分析實(shí)驗(yàn)數(shù)據(jù)的時(shí)候有一定的參考價(jià)值.但是這2種模型都有一定的局限性,一致轉(zhuǎn)動(dòng)模型要求晶粒之間完全耦合,磁矩在磁化反轉(zhuǎn)過程中磁矩矢量保持一致,而完全退耦合模型忽略了硬/軟磁材料磁晶之間的相互作用,軟磁相磁矩反轉(zhuǎn)和硬磁相磁矩反轉(zhuǎn)是獨(dú)立完成的,這2種模型并不能精確描述出磁矩在反轉(zhuǎn)過程中的分布,計(jì)算出的磁滯回線與實(shí)驗(yàn)差距也比較大.

    2一維半解析模型

    前面介紹了簡(jiǎn)單的一維全解析模型,接下來討論一維半解析模型.H. Kronmüller等[79,86-87]、G. Asti等[7-8,56-60]以及G. P. Zhao等[88,92]運(yùn)用這種模型研究了多種納米復(fù)合材料的磁滯回線,并推導(dǎo)出成核場(chǎng)公式以及厚度較大時(shí)的釘扎場(chǎng)公式.本文以外場(chǎng)和易軸取向平行,并且均在薄膜平面內(nèi)的情況為例,系統(tǒng)分析硬磁軟磁多層膜磁體在退磁過程中的特點(diǎn).薄膜的橫向尺寸比其厚度大很多,因此可以假設(shè)膜面無限延伸,磁矩保持在膜面內(nèi)偏轉(zhuǎn),靜磁相互作用可以忽略不計(jì).同時(shí)該模型假設(shè)磁矩在薄膜平面內(nèi)一致轉(zhuǎn)動(dòng),即磁矩只在厚度方向有一個(gè)分布,以便將退磁問題簡(jiǎn)化成為一維模型[7,9,12,93-94],通過解析分析來處理.根據(jù)W. F. Brown[95]的微磁學(xué)理論,易軸平行于膜面的三層膜體系的能量密度可以表示為

    (6)

    易軸垂直于膜面時(shí)的總能量密度則表示為

    (7)

    右邊的第一項(xiàng)能量表示交換能,后面兩項(xiàng)分別表示磁晶各向異性能和塞曼能,垂直時(shí)多了一項(xiàng)退磁能.這里的A、K和Ms分別表示交換作用常數(shù)、磁晶各向異性常數(shù)和飽和自發(fā)磁化強(qiáng)度,上標(biāo)h和s分別代表硬磁相和軟磁相,θ表示磁化矢量和外場(chǎng)方向的夾角.

    系統(tǒng)的邊界條件[8,88]為:

    (9)

    將(6)式代入歐拉方程,可以得到θ與z之間的關(guān)系:

    其中,Δ=π(A/K)1/2表示材料的布洛赫寬度,Hk=2K/Ms表示磁晶各向異性場(chǎng),h=H/Hk表示約化外場(chǎng).

    體系的邊界條件為

    將(10)和(11)式代入上面公式,邊界條件可以改寫成

    (12)

    在三層膜中,硬磁相所占的比例為2th/(ts+2th);而在多層膜中,硬磁相所占的比例為th/(ts+th).對(duì)三層膜來說,考慮到真空磁矩為零,體系的邊界條件為

    在多層膜中,考慮到對(duì)稱性,也可以得到同樣的體系邊界條件

    因此三層膜體系中硬磁相的厚度th可以對(duì)應(yīng)多層膜中的2th.一個(gè)硬磁相厚度為th的三層膜,其磁滯回線和一個(gè)硬磁相厚度為2th的多層膜相同.因此,多層膜最大磁能積所對(duì)應(yīng)的厚度是三層膜的2倍.而一個(gè)軟磁相厚度為ts的雙層膜則與一個(gè)軟磁相厚度為2ts的三層膜對(duì)應(yīng),因此,三層膜最大磁能積所對(duì)應(yīng)的厚度是雙層膜的2倍.

    圖5表示的是運(yùn)用一維半解析方法得到Fe/Sm-Co相同硬軟磁相厚度不同膜層數(shù)的退磁曲線[96](其總層數(shù)N=Nh+Ns=200,對(duì)于nd=3表示h/s/h三層膜;nd=5表示h/s/h/s/h五層膜).可以看到這種材料的退磁過程,從圖中可以分析出成核場(chǎng)、矯頑力和釘扎場(chǎng)的變化.對(duì)于nd=11,磁滯回線接近于矩形,此時(shí)磁體為剛性磁體;nd=7,對(duì)應(yīng)的磁體為交換彈簧磁體;nd=3,對(duì)應(yīng)的磁體為退耦合磁體.一維解析微磁學(xué)方法可以有效解決一些物理問題,比如給出一維疇壁.但是這種分析方法有一定的局限性,由于數(shù)學(xué)表達(dá)式的復(fù)雜性,在推導(dǎo)過程中需要采取一些近似的方法,這可能導(dǎo)致計(jì)算結(jié)果誤差較大.另外,這種方法適用范圍受到一定限制,對(duì)于圓柱形、線形和錐形等形狀計(jì)算起來較為困難.

    3三維模型

    三維數(shù)值模擬在過去十多年發(fā)展迅速,用這種方法可以模擬不同的磁體形狀,比如薄膜、球體、圓柱體等,甚至不規(guī)則的形狀,比如錐形體、線形等[97].對(duì)于較為復(fù)雜的不規(guī)則形狀磁體,用三維數(shù)值模擬方法處理更為方便.R.P.Boardman等[98-101]用這種方法研究了不同形狀的磁體磁化反轉(zhuǎn)過程,G.P.Zhao等[82-84]用三維有限差分方法研究了納米復(fù)合磁體退磁過程中,磁矩形成的形態(tài)以及磁滯回線和磁能積.這種三維微磁學(xué)計(jì)算方程是基于LLG[102]動(dòng)力學(xué)方程

    (13)

    (14)

    平均能量密度E可以用布朗方程[95]得到

    (15)

    A和K分別表示交換常數(shù)和磁晶各向異性常數(shù),H表示外加磁場(chǎng),Hd(r)表示退磁場(chǎng).

    G.P.Zhao等[82-84]成功地應(yīng)用以上方法計(jì)算了易軸平行和垂直于膜面的磁滯回線,并與一維模型計(jì)算結(jié)果進(jìn)行了細(xì)致比較.圖7表示Nd2Fe14B/α-Fe雙層膜易軸平行于膜面不同軟磁相厚度的磁滯回線[83],其中上半部分由三維模擬計(jì)算得出,下半部分則根據(jù)一維半解析方法得到.由圖7可見,二者符合得很好,證明了計(jì)算方法和結(jié)果的可靠性.在軟磁相厚度很小的時(shí)候,磁滯回線為矩形,對(duì)應(yīng)的復(fù)合磁體可以看成是剛性磁體,在厚度處于中間值時(shí),磁滯回線的矩形度比剛性磁體差,此時(shí)表現(xiàn)出交換彈簧行為,對(duì)應(yīng)的復(fù)合磁體磁相是交換彈簧.繼續(xù)增大ts,成核場(chǎng)與釘扎場(chǎng)明顯分離開來,對(duì)應(yīng)的復(fù)合磁體磁相從交換彈簧過渡到退耦合.G.Asti等[7-8]的一維計(jì)算也給出了類似的結(jié)果.

    其他材料的計(jì)算也得到了類似的結(jié)果.圖8[84]是易軸垂直于膜面時(shí),不同軟磁相厚度下的Sm-Co/α-Fe三層膜的磁滯回線.圖中實(shí)線是三維微磁學(xué)模擬得到的結(jié)果,虛線是一維模擬計(jì)算的結(jié)果.垂直取向與平行取向變化相似,隨著軟磁相厚度ts的增大,對(duì)應(yīng)的復(fù)合磁體磁相從剛性到交換彈簧最后過渡到退耦合.在圖中,可以發(fā)現(xiàn)三維數(shù)值模擬結(jié)果和一維解析計(jì)算結(jié)果很相近,只有在軟磁相厚度比較大的時(shí)候差距才比較大,成核場(chǎng)和矯頑力都隨著軟磁相厚度的增加而減小.

    4薄膜厚度對(duì)磁滯回線的影響

    4.1軟磁相厚度對(duì)磁滯回線的影響薄膜厚度特別是軟磁相厚度對(duì)磁滯回線有著比較大的影響,隨著軟磁相厚度的增加,使得磁滯回線的方形度變差,導(dǎo)致磁能積急劇地減小,T.Leineweber等[12]、E.E.Fullerton等[105-108]以及G.Asti等[7-8]對(duì)此進(jìn)行了深入研究.

    圖9中描述的是外場(chǎng)垂直于膜面時(shí)不同軟磁相厚度下Fe/FePt雙層膜的磁滯回線[8].單相永磁體(軟磁相厚度為零)的磁滯回線是一個(gè)矩形.隨著軟磁相厚度的增加,磁滯回線的方形度急劇下降,矯頑力也開始降低.在軟磁相厚度很小的時(shí)候,成核場(chǎng)和釘扎場(chǎng)之間的間隙很小,磁滯回線接近方形,成核和磁矩的完全偏轉(zhuǎn)幾乎在同一外場(chǎng)下完成,此時(shí)矯頑力機(jī)制主要以成核為主,對(duì)應(yīng)的磁體為剛性復(fù)合磁體.在軟磁相厚度較大的時(shí)候,矯頑力機(jī)制以釘扎為主.隨著軟磁相厚度的增加,磁滯回線的方形度逐漸變差,成核場(chǎng)和矯頑力逐漸分離.當(dāng)軟磁相厚度較大時(shí),疇壁擴(kuò)展過程變長(zhǎng),因此成核到釘扎的歷程變長(zhǎng),因此磁滯回線的方形度和磁能積都下降.同時(shí),軟磁相厚度的增加會(huì)導(dǎo)致矯頑力迅速降低;另一方面,如果軟磁相厚度太小,雜質(zhì)的影響會(huì)較大,導(dǎo)致矯頑力和磁能積降低.

    4.2硬磁相厚度的影響硬磁相厚度對(duì)磁滯回線的影響相對(duì)較小,僅在厚度比較小的時(shí)候影響比較明顯.圖10給出了不同硬磁相厚度下Nd2Fe14B/α-Fe三層膜的磁滯回線[109],計(jì)算中軟磁相厚度為80nm.

    由于軟磁層比較厚,圖中的幾個(gè)磁滯回線與圖4類似,表現(xiàn)出明顯的退耦合行為.在圖10(a)中選取的硬磁相厚度較小,由于軟磁相厚度遠(yuǎn)大于硬磁相厚度,磁滯回線的主要貢獻(xiàn)來自于軟磁相的磁矩偏轉(zhuǎn),此時(shí)磁滯回線的主導(dǎo)矯頑力機(jī)制是成核,隨著硬磁相厚度的增大,成核場(chǎng)沒有變化,釘扎場(chǎng)逐漸增大.在圖10(b)中硬磁相厚度較大,隨著硬磁相厚度的增加,成核場(chǎng)和釘扎場(chǎng)均沒有變化,但是磁滯回線的主要貢獻(xiàn)由軟磁相的偏轉(zhuǎn)逐漸過渡到硬磁相的偏轉(zhuǎn),矯頑力與成核場(chǎng)的差距逐漸增大,矯頑力機(jī)制從成核過渡到釘扎.

    5易軸取向的影響

    由R.Skomski和J.M.D.Coey兩位科學(xué)家提出的理論磁能積是基于易軸與外加磁場(chǎng)方向完全一致的一維模型,并且要求軟磁相材料結(jié)構(gòu)無缺陷.在實(shí)際的磁性材料中,硬磁相和軟磁相的易軸與外場(chǎng)都存在一定的夾角.在實(shí)驗(yàn)[45,110]中發(fā)現(xiàn),當(dāng)外場(chǎng)為零時(shí),磁性材料的磁矩并不完全垂直或者平行于膜面,而是出現(xiàn)了一定的分布,也就是說磁晶易軸與外加磁場(chǎng)不會(huì)完全平行,而是存在一定的夾角.很多的理論計(jì)算都沒有考慮磁晶易軸的有效偏角[60,88],因此這是磁能積實(shí)驗(yàn)值小于理論計(jì)算值的一個(gè)重要原因.

    圖11給出了th=10nm,ts=20nm易軸在膜面內(nèi)不同β值下的Nd2Fe14B/α-Fe/Nd2Fe14B三層膜的退磁曲線[51].從這些曲線中不難看出,β對(duì)于磁滯回線有著很大的影響,隨著β的增加磁滯回線變得越來越傾斜,剩磁也急劇減小.當(dāng)外場(chǎng)為零時(shí),由布朗方程可以知道,薄膜內(nèi)部磁矩的方向與易軸的方向保持一致,剩磁可以表示為Ms=Mr×cosβ,M隨著β的增大而減小.當(dāng)M為零時(shí),對(duì)應(yīng)外場(chǎng)的相反數(shù)是矯頑力Hc,表征磁體抵抗退磁能力的大小.從圖中可以看出,矯頑力隨著β的增大而減小,當(dāng)β等于0時(shí),矯頑力機(jī)制為釘扎,隨著β的增加,矯頑力小于釘扎場(chǎng)并且二者的差距逐漸增大.文獻(xiàn)[52-54]也給出了類似結(jié)果.由于磁晶易軸偏角的存在,導(dǎo)致剩磁和矯頑力比理想情況小,引起磁能積的減小,在一定程度上解釋了磁能積的理論值和實(shí)驗(yàn)值之間的巨大差距.

    6理論和實(shí)驗(yàn)的對(duì)比

    微磁學(xué)長(zhǎng)期以來一直為部分人所詬病,主要原因在于理論計(jì)算給出的矯頑力和磁能積比實(shí)驗(yàn)大很多,這其實(shí)不是微磁學(xué)本身的問題.一方面,理論計(jì)算采取的模型過于簡(jiǎn)化;另一方面實(shí)驗(yàn)制造出的材料還有很多缺陷,這一差距較好地反映在W.F.Brown[74-75]提出的矯頑力悖論和G.P.Zhao等[9]提出的磁能積悖論中.G.P.Zhao等[28-29,89-91]進(jìn)一步指出,隨著實(shí)驗(yàn)手段的不斷豐富和計(jì)算條件的不斷改善,微磁學(xué)理論和實(shí)驗(yàn)的差別越來越小.圖12給出了納米復(fù)合材料Nd2Fe14B/Fe65Co35多層膜磁滯回線[111]最新的理論和實(shí)驗(yàn)對(duì)比.理論計(jì)算所選取的參數(shù)根據(jù)實(shí)驗(yàn)條件所確定.從圖12可以看出理論得到的結(jié)果與實(shí)驗(yàn)結(jié)果比較接近,特別是理論得到的剩磁和矯頑力都和實(shí)驗(yàn)值差別不大.這也比較有力地證明了微磁學(xué)的可信度.

    7討論和結(jié)論

    微磁學(xué)可以用來研究硬/軟磁復(fù)合材料的磁滯回線,是比較有力且值得信賴的研究方法.本文中分別討論了一維全解析、一維半解析以及三維數(shù)值模擬得到的磁滯回線,簡(jiǎn)單地介紹了這3種方法的優(yōu)點(diǎn)以及缺點(diǎn).在用一維全解析和三維數(shù)值模擬這種模型時(shí),它可以明確給出磁化反轉(zhuǎn)過程,尤其是疇壁的形成、擴(kuò)展和釘扎.充分模擬疇壁的擴(kuò)展,可以發(fā)現(xiàn)材料的矯頑力機(jī)制為自釘扎.而材料本身也存在各種雜質(zhì)和缺陷,因此,修復(fù)材料本身缺陷和阻礙疇壁的擴(kuò)展是提高矯頑力的2種有效方法.對(duì)復(fù)合磁體三維數(shù)值計(jì)算來說,最重要的是網(wǎng)格的選取.一個(gè)永磁體的磁疇壁的寬度大約是3~5nm.因此,網(wǎng)格尺寸應(yīng)小于1nm×1nm×1nm,這樣在每個(gè)單元內(nèi)有著36°~60°的疇壁變化.這樣的精度對(duì)取向的復(fù)合磁多層膜來說基本可以滿足要求.由此,一個(gè)100nm×100nm×100nm的計(jì)算模型,需要106個(gè)網(wǎng)格,這個(gè)數(shù)量也是一臺(tái)計(jì)算機(jī)的計(jì)算極限.有些時(shí)候,三維計(jì)算并不是最適合的,而一維解析方法可以取得更好的結(jié)果,總的來說,在計(jì)算能力受限的情況下,將一維和三維計(jì)算結(jié)合起來并將計(jì)算結(jié)果細(xì)致比較,取長(zhǎng)補(bǔ)短,是一個(gè)較好的方法.

    由以上討論可知:矯頑力的降低是引起復(fù)合磁

    體磁能積減小的一個(gè)重要原因.因此,G.P.Zhao等[85]提出,在實(shí)驗(yàn)上可以通過提高矯頑力來提高磁能積,比如在硬/軟磁材料之間添加一層非鐵磁層材料來避免兩相材料在界面的擴(kuò)散,從而提高材料的矯頑力.日本的K.Hono小組[23]和中國燕山大學(xué)的張湘義小組[112]用不同方法提高了界面的釘扎,在這方面取得了成功.磁性材料的易軸取向也會(huì)導(dǎo)致剩磁和磁能積迅速減小[52-54,111],控制好薄膜的易軸取向和材料的均勻性都可以提高磁能積.此外,較好地掌控薄膜厚度,也是提高磁能積的途徑之一[7-8,12,105-108].

    本文討論了硬/軟磁相厚度對(duì)磁滯回線的影響,分析了矯頑力和成核場(chǎng)的變化.可以得到軟磁相厚度的增加會(huì)使得磁滯回線的方形度變差,矯頑力也逐漸降低,矯頑力機(jī)制逐漸從成核向釘扎過渡.由于軟磁相厚度較大時(shí),疇壁擴(kuò)展過程變長(zhǎng),成核到釘扎的歷程變長(zhǎng),導(dǎo)致磁滯回線的方形度和磁能積都下降.隨著硬磁相厚度的增加,矯頑力與成核場(chǎng)的差距逐漸增大,矯頑力機(jī)制從成核過渡到釘扎,此時(shí)釘扎場(chǎng)的大小主要由硬磁相決定.除此之外,也討論了磁晶易軸的取向?qū)Υ艤鼐€的影響,當(dāng)易軸偏離逐漸增大,磁滯回線的方形度嚴(yán)重受損.體系的剩磁單調(diào)減小,矯頑力也呈減小的趨勢(shì),這是磁能積實(shí)驗(yàn)值比理論值小很多的一個(gè)重要原因.

    參考文獻(xiàn)

    [1] COEHOORN R, DEMOOIJ D B, DEWAARD C. Meltspun permanent magnet materials containing Fe3B as the main phase[J]. J Magn Magn Mater,1989,80(1):101-104.

    [2] COEHOORN R, DEMOOIJ D B, DUCHATEAU J, et al. Novel permanent magnetic materials made by rapid quenching[J]. J de Physique,1988,49(C-8):669-670.

    [3] KNELLER E F, HAWING R. Exchange-spring magnet:a new material principle for permanent magnets[J]. IEEE Trans Magn,1991,27(4):3588-3591.

    [4] SKOMSKI R, COEY J. Giant energy product in nanostructured two-phase magnets[J]. Phys Rev,1993,B48(21):15812-15816.

    [5] LIU W, ZHANG Z D, LIU J P, et al. Exchange coupling and remanence enhancement in nanocomposite multilayer magnets[J]. Adv Mater,1993,14(24):1832-1834.

    [6] LIU W, ZHANG Z D, LIU J P, et al Nanocomposite (Nd,Dy)(Fe,Co,Nb,B)5.5/α-Fe multilayer magnets with high performance[J]. J Phys D:Appl Phys,2003,36(17):L63-L66.

    [7] ASTI G, SOLZI M, GHIDINI M, et al. Micromagnetic analysis of exchange-coupled hard-soft planar nanocomposites[J]. Phys Rev,2004,B69:174401.

    [8] ASTI G, GHIDINI M, PELLICELLI R, et al. Magnetic phase diagram and demagnetization processes in perpendicular exchange-spring multilayers[J]. Phys Rev,2006,B73:094406.

    [9] ZHAO G P, WANG X L. Nucleation, pinning, and coercivity in magnetic nanosystems:an analytical micromagnetic approach[J]. Phys Rev,2006,B74:012409.

    [10] ZHAO G P, CHEN L, HUANG C W, et al. Micromagnetic calculation of hysteresis loops in exchange-coupled nanolayers[J]. Solid State Commun,2010,150(31):1486-1488.

    [11] ZHAO G P, ZHOU G, ZHANG H W, et al. Angular dependence of magnetic reversal in two-phased nanolayers[J]. Comput Mater Sci,2008,44(1):117-121.

    [12] LEINEWEBER T, KRONMüLLER H. Micromagnetic examination of exchange coupled ferromagnetic nanolayers[J]. J Magn Magn Mater,1997,176(97):145-154.

    [13] JIANG J S, PEARSON J E, LIU Z Y, et al. Improving exchange-spring nanocomposite permanent magnets[J]. Appl Phys Lett,2004,85(22):5293-5295.

    [14] CHOI Y, JIANG J S, PEARSON J E, et al. Controlled interface profile in Sm-Co/Fe exchange-spring magnets[J]. Appl Phys Lett,2007,91(7):072509-3.

    [15] SAWATZKI S, HELLER R, MICKEL Ch, et al. Largely enhanced energy density in epitaxial SmCo5/Fe/SmCo5exchange spring trilayers[J]. J Appl Phys,2011,109(12):123922-7.

    [16] BELEMUK A M, CHUI S T. Comparative study of finite temperature demagnetization in Nd2Fe14B and SmCo5based hard-soft composites[J]. J Appl Phys,2011,110(7):073918-7.

    [17] BELEMUK A M, CHUI S T. Temperature-dependent demagnetization behaviour in perpendicular exchange-coupled SmCo5/FeCo multilayers[J]. J Phys D:Appl Phys,2012,45(12):125001-125014.

    [18] STONER E C, WOHLFARTH E P. A mechanism of magnetic hysteresis in heterogeneous alloys[J]. Magnetics,Phil Trans Roy Soc London,1948,A240:599-642.

    [19] NéEL L. Les lois de l’aimantation et de la subdivision en domaines élémentaires d'un monocristal de fer[J]. J Phys Radium,1944,5(11):265-276.

    [20] ZHAO G P, ZHAO M G, LIM H S, et al. From nucleation to coercivity[J]. Appl Phys Lett,2005,87(16):162513-3.

    [21] ZHAO G P, CHEN L, HUANG C W, et al. Evolution of domain walls and reversal mechanism in exchange-coupled nanolayers[J]. J Magn Magn Mater,2009,321(15):2322-2326.

    [22] ATKINSON D, ALLWOOD D A, XIONG G, et al. Magnetic domain-wall dynamics in a submicrometre ferromagnetic structure[J]. Nature Material,2003,2(2):85-87.

    [23] CUI W B, TAKAHASHI Y K, HONO K. Nd2Fe14B/FeCo anisotropic nanocomposite films with a large maximum energy product[J]. Adv Mater,2012,24:6530-6535.

    [24] 姜壽,李衛(wèi). 凝聚態(tài)磁性物理[M]. 北京:科學(xué)出版社,2003:344.

    [25] ARCAS J, HERNANDO A, BARANDIARn J, et al. Soft to hard magnetic anisotropy in nanostructured magnets[J]. Phys Rev,1998,B58(9):5193-5196.

    [26] FENG W C, GAO R W, YAN S S, et al. Effects of phase distribution and grain size on the effective anisotropy and coercivity of nanocomposite Nd2Fe14B/α-Fe magnets[J]. J Appl Phys,2005,98(4):044305-5.

    [27] HAN G B, GAO R W, FU S, et al. Effective anisotropy between magnetically soft and hard grains in nanocomposite magnets[J]. Appl Phys,2005,A81(3):579-582.

    [28] ZHAO G P, MORVAN F, WAN X L. Micromagnetic calculation for exchange-coupled nanocomposite permanent magnets[J]. Rev Nanosci Nanotechnol,2014,3(4):227-258.

    [29] ZHAO G P, ZHANG X F, MORVAN F. Theory for the coercivity and its mechanisms in nanostructured permanent magnetic materials[J]. Rev Nanosci Nanotechnol,2015,4:1-25.

    [30] JULLIERE M. Tunneling between ferromagnetic films[J]. Phys Lett,1975,A54(3):225-226.

    [31] SLONCZEWSKI J C. Conductance and exchange coupling of two ferromagnets separated by a tunneling barrier[J]. Phys Rev,1989,B39(10):6995-7002.

    [32] GIEBELER C, ADELERHOF D J, KUIPER A E T, et al. Robust GMR sensors for angle detection and rotation speed sensing[J]. Sensors and Actuators,2001,A91(1):16-20.

    [33] DAUGHTON J M. GMR and SDT sensor applications[J]. IEEE Trans Magn,2000,36(5):2773-2778.

    [34] CHRISTIDES C, STAVROYIANNIS S, KALLIAS G, et al. Giant magnetoresistance Co/Cu multilayer sensors for use in magnetic field mapping[J]. Sensors and Actuators,1999,A76(1):167-171.

    [35] SHENG L, XING D Y, SHENG D N, et al. Theory of colossal magnetoresistance in R1-xAxMnO3[J]. Phys Rev Lett,1997,79(9):1710-1713.

    [36] MILLIS A J, LITTLEWOOD P B, SHRAIMAN B I. Double exchange alone does not explain the resistivity of La1-xSrxMnO3[J]. Phys Rev Lett,1995,74(25):5144-5147.

    [37] WARREN E P, DAVID J S. Electronic structure and half-metallic transport in the La1-xCaxMnO3system[J]. Phys Rev,1996,B53(3):1146-1160.

    [38] LI X W, GUPTA A, XIAO G, et al. Low-field magnetoresistive properties of polycrystalline and epitaxial perovskite manganite films[J]. Appl Phys Lett,1997,71(8):1124-1126.

    [39] HWANG H Y, CHEONG S W, ONG N P, et al. Spin-polarized intergrain tunneling in La2/3Sr1/3MnO3[J]. Phys Rev Lett,1996,77(2):2041-2044.

    [40] SCHIFFER P, RAMIREZ A P, BAO W, et al. Low temperature magnetoresistance and the magnetic phase diagram of La1-xCaxMnO3[J]. Phys Rev Lett,1995,75(18):3336-3339.

    [41] GRIFFITHS M K, BISSHOP J E L, DAVIES H A. Computer simulation of single-phase nanocrystalline permanent magnets[J]. J Magn Magn Mater,1999,183(1):49-67.

    [42] 王博文,葛景巖. 巨磁致伸縮材料及其應(yīng)用[J]. 沈陽工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),1998,20(3):64-67.

    [43] CLARK A E, WUN-FOGLE M, RESTORFF J B, et al. Magnetostriction and magnetomechanical coupling of grain oriented Tb0.6Dy0.4sheet[J]. IEEE Trans Magn,1993,29(6):3511-3513.

    [44] DUC N H, MACKAY K, BETZ J, et al. Giant magnetostriction in amorphous (Tb1-xDyx)(Fe0.45Co0.55)yfilms[J]. J Appl Phys,1996,79(2):973-977.

    [45] ZHANG J, TAKAHASHI Y K, GOPALAN R, et al. Sm(Co,Cu)5/Fe exchange spring multilayer films with high energy product[J]. Appl Phys Lett,2005,86(12):122509-3.

    [46] YUE M, ZHANG J, TIAN M. Microstructure and magnetic properties of Isotropic bulk NdxFe94-xB6(x=6,8,10) nanocomposite magnets prepared by spark plasma sintering[J]. J Appl Phys,2006,99(8):08B502-3.

    [47] ZHANG J, LI Y X, WANG F, et al. Coercivity mechanism of nanocomposite Sm-Co/Fe multilayer films[J]. J Appl Phys,2010,107(4):043911.

    [48] LIU W Q, ZUO J H, YUE M, et al. Strucure and magnetic properties if bulk anisotropic SmCo5/α-Fe nanocomposite permanent magnets with different α-Fe content[J]. J Appl Phys,2011,109(7):07A741-3.

    [49] LEINEWEBER T, KRONMüLLER H. Magnetization reversal modes in inhomogeneous magnets[J]. Phys Stat Sol,1997,201(1):291-301.

    [50] SAWATZKI S, HELLER R, MICKEL C, et al. Largely enhanced energy density in epitaxial SmCo5/Fe/SmCo5exchange spring trilayers[J]. J Appl Phys,2011,109(12):123922-7.

    [51] XIA J, ZHAO G P, ZHANG H W, et al. Significant deterioration of energy products in exchange-coupled composite magnets[J]. J Appl Phys,2012,112(1):013918-5.

    [52] 夏靜,張溪超,趙國平. 易軸取向?qū)d2Fe14B/α-Fe雙層膜退磁過程影響的微磁學(xué)分析[J]. 物理學(xué)報(bào),2013,62(22):227502.

    [53] 夏靜,趙國平,張敏,等. 易軸取向?qū)τ?軟磁多層膜磁性能的影響[J]. 中國科學(xué):物理學(xué)·力學(xué)·天文學(xué),2013,43(1):39-47.

    [54] WAN X L, ZHAO G P, ZHANG X F, et al. Hysteresis of misaligned hard-soft grains[J]. J Magn Magn Mater,2016,397:181-187.

    [55] GUO Z J, JIANG J S, PEARSON J E, et al. Exchange-coupled Sm-Co/Nd-Co nanomagnets:correlation between soft phase anisotropy and exchange field[J]. Appl Phys Lett,2002,81(11):2029-2031.

    [56] ASTI G, CARBUCICCHIO M, RATEO M, et al. Exchange coupling in nanometric Fe/Co multilayers[J]. J Magn Magn Mater,1999,59:196-197.

    [57] GHIDINI M, ASTI G, PELLICELLI R, et al. Hard/soft composite magnets[J]. J Magn Magn Mater,2007,316(11):159-165.

    [58] ASTI G, CARBUCICCHIO M, GHIDINI M, et al. Microstructural and magnetic properties of exchange-coupled Co/Fe multilayers[J]. J Appl Phys,2000,87(9):6689-6691.

    [59] ASTI G, GHIDINI M, NERI F M, et al. Flux reversal in hard-soft composite magnets[J]. J Magn Magn Mater,2004,272(2):650-651.

    [60] ASTI G, SOLZI M, GHIDINI M. Reversal modes of the multilayer exchange-spring magnet[J]. J Magn Magn Mater,2001,226/230:1464-1466.

    [61] LIU Y G, XU L, WANG Q F, et al. Development of crystal texture in Nd-lean amorphous Nd9Fe85B6under hot deformation[J]. Appl Phys Lett,2009,94(17):172502-3.

    [62] CHENG Z H, ZHANG J X, KRONMüLLER H. Magnetically soft phase in magnetization reversal processes of nanocomposite Sm2Fe15Ga2Cx/α-Fe permanent magnetic materials[J]. Phys Rev,2003,B68(14):144417.

    [63] LI W, LI L L, NAN Y, et al. Controllable nanocrystallization in amorphous Nd9Fe85B6via combined application of severe plastic deformation and thermal annealing[J]. Appl Phys Lett,200,91(6):062509-3.

    [64] GABAY A M, ZHANG Y, HADJIPANAYIS G C. Die-upset hybrid Pr-Fe-B nanocomposite magnets[J]. Appl Phys Lett,2004,85(3):446-448.

    [65] LEE D, BAUSER S, HIGGINS A, et al. Bulk anisotropic composite rare earth magnets[J]. J Appl Phys,2006,99(8):08B516-3.

    [66] LIU J P, LUO C P, LIU Y, et al. High energy products in rapidly annealed nanoscale Fe/Pt multilayers[J]. Appl Phys Lett,1998,72(4):483-485.

    [67] LI B Z, ZHANG M, SHEN B G, et al. Non-uniform magnetization reversal in nanocomposite magnets[J]. Appl Phys Lett,2013,102(10):102405-4.

    [68] YAO Q, LIU W, CUI W B, et al. Growth mechanism and magnetic properties for the out-of-plane-oriented Nd-Fe-B films[J]. J Mater Res,2009,24(9):2802-2812.

    [69] GOLL D, BREITLING A, GU L, et al. Experimental realization of graded L10-FePt/Fe composite media with perpendicular magnetization[J]. J Appl Phys,2008,104(8):083903-4.

    [70] ZHOU T J, LIM B C, LIU B. Anisotropy graded FePt-TiO2nanocomposite thin films with small grain size[J]. Appl Phys Lett,2009,94(5):152505-3.

    [71] PANDEY K K M, CHEN J S, CHOW G M, et al. L10CoPt-Ta2O5exchange coupled multilayer media for magnetic recording[J]. Appl Phys Lett,2009,94(23):232502-3.

    [72] BELEMUK A M, CHUI S T. Comparative study of finite temperature demagnetization in Nd2Fe14B and SmCo5based hard-soft composites[J]. J Appl Phys,2011,110(7):073918.

    [73] BELEMUK A M, CHUI S T. Temperature-dependent demagnetization behaviour in perpendicular exchange-coupled SmCo5/FeCo multilayers[J]. J Phys D:Appl Phys,2012,45(12):125001.

    [74] BROWN W F. Theory of the approach to magnetic saturation[J]. Phys Rev,1940,58(8):736-743.

    [75] BROWN W F. The effect of dislocation on magnetization near saturation[J]. Phys Rev,1941,60(2):139-147.

    [76] ZHAO G P, WANG X L, FENG Y P, et al. Coherent rotation and effective anisotropy[J]. IEEE Trans Magn,2007,43(6):2908-2910.

    [77] 李建明,郭乃理,李貌,等. 硬/軟磁復(fù)合材料的一致轉(zhuǎn)動(dòng)模型及其磁滯回線[J]. 四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,36(2):263-268.

    [78] FISCHER R, KRONMüLLER H. Importance of ideal grain boundaries of high remanent composite permanent magnets[J]. J Appl Phys,1998,83(6):3271-3275.

    [79] KRONMüLLER H, GOLL D. Magnetic properties of exchange-spring composite films[J]. Solid State Physics,2011,248(10):2361-2367.

    [80] LI Y Q, YUE M, ZUO J H, et al. Investigation of magnetic properties of MnBi/α-Fe nanocomposite permanent magnets by micromagnetic simulation[J]. IEEE Trans Magn,2013,49(7):3391-3393.

    [81] LI Y Q, YUE M, ZUO J H, et al. Effect of α-Fe content on the magnetic properties of MnBi/α-Fe nanocomposite permanent magnets by micromagnetic calculation[J]. J Magn,2013,18(3):245-249.

    [82] YUAN X H, ZHAO G P, YUE M, et al. 3D and 1D calculation of hysteresis loops and energy products for anisotropic nanocomposite films with perpendicular anisotropy[J]. J Magn Magn Mater,2013,343(5):245-250.

    [83] ZHANG W, ZHAO G P, YUAN X H, et al. 3D and 1D micromagnetic calculation for hard/soft bilayers with in-plane easy axes[J]. J Magn Magn Mater,2015,324(24):4231-4236.

    [84] ZHANG X C, ZHAO G P, XIA J, et al. Micromagnetic simulation of Sm-Co/α-Fe/Sm-Co trilayers with various angles between easy axes and the film plane[J]. Chin Phys,2014,B23(9):097504-7.

    [85] SI W J, ZHAO G P, RAN N, et al. Deterioration of the coercivity due to the diffusion induced interface layer in hard/soft multilayers[J]. Scientific Reports,2015,5:16212.

    [86] KRONMüLLER H, GOLL D. Micromagnetic analysis of nucleation-hardened nanocrystalline PrFeB magnets[J]. Scripta Mater,2002,47(8):551-556.

    [87] FISCHER R, KRONMüLLER H. Static computational micromagnetism of demagnetization processes in nanoscaled permanent magnets[J]. Phys Rev,1996,B54(10):7284-7294.

    [88] ZHAO G P, DENG Y, ZHANG H W, et al. Accurate calculation of the nucleation field and hysteresis loops in hard-soft multilayers[J]. J Appl Phys,2011,109(7):07D340-3.

    [89] YUE M, LI Y Q, WU Q, et al. Bulk nanostructural permanent magnetic materials[J]. Rev Nanosci Nanotechnol,2014,3(4):276-288.

    [90] LIU W Q, YUE M, CUI B Z, et al. Permanent magnetic nanoparticles and nanoflakes prepared by surfactant-assisted high-energy ball milling[J]. Rev Nanosci Nanotechnol,2014,3(4):259-275.

    [91] ZHANG J, YAN A R, LIU J P. Permanent magnetic films[J]. Rev Nanosci Nanotechnol,2014,3(4):289-303.

    [92] ZHAO G P, YANG C, XIAN C W, et al. Analytical determination of nucleation field and magnetic reversal modes in exchange-coupled nanolayers[J]. Mod Phys Lett,2009,B23(25):2955-2961.

    [93] PELLICELLI R, SOLZI M, NEU V, et al. Characterization and modeling of the demagnetization processes in exchange coupled SmCo5/Fe/SmCo5trilayers[J]. Phys Rev,2010,B81(18):184430.

    [94] WILSON M, ZHU M, MYERS R, et al. Interlayer and interfacial exchange coupling in ferromagnetic metal/semiconductor heterostructures[J]. Phys Rev,2010,B81(4):045319.

    [95] BROWN W F. Virtues and weaknesses of the domain concept[J]. Rev Mod Phys,1945,17(1):15-19.

    [96] AMATO M, PINI M G, RETTORI A. Optimization study of the nanostructure of hard/soft magnetic multilayers[J]. Phys Rev,1999,B60(5):3414-3420.

    [97] SKOMSKI R, LIU J P, SELLMYER D J. Quasicoherent nucleation mode in two-phase nanomagnets[J]. Phys Rev,1999,B60(10):7359-7365.

    [98] ZZIMMERMANN J R P, BORDIGNON G, BOARDMAN R P, et al. Micromagnetic simulation of the magnetic exchange spring system DyFe2/YFe2[J]. J Appl Phys,2006,99(8):08B904-3.

    [99] BOARDMAN R P, ZIMMERMANN J R P, FANGOHR H, et al. Micromagnetic simulation studies of ferromagnetic part-spheres[J]. J Appl Phys,2005,97(10):10E305-3.

    [100] BOARDMAN R P, FANGOHR H, COX S J, et al. Micromagnetic simulation of ferromagnetic part-spherical particles[J]. J Appl Phys,2004,95(11):7037-7039.

    [101] BOARDMAN R P. Computer simulation studies of magnetic nanostructures[D]. Southampton:University of Southampton,2005.

    [102] DONAHUE M J, PORTER D G. OOMMF User’s Guide, Version 1.0[M/OL]. National Institute of Standards and Technology,1999. http://math.nist.gov/oommf/.

    [103] FULLERTON E E, JIANG J S, BADER S D. Hard/soft magnetic heterostructures:model exchange-spring magnets[J]. J Magn Magn Mater,1999,200(1):392-404.

    [104] CHEN B Z, ZHAO G P, ZHANG H W, et al. 3D calculation of hysteresis loops, magnetic orientations and reversal processes for exchange-spring bilayers with perpendicular anisotropy[J]. Physica,2012,B407:4574-4578.

    [105] FULLERTON E E, JIANG J S, GRIMSDITCH M, et al. Exchange-spring behavior in epitaxial hard/soft magnetic bilayers[J]. Phys Rev,1998,B58:12193-12200.

    [106] FULLERTON E E, JIANG J S, SOWERS C H. Structure and magnetic properties of exchange-spring Sm-Co/Co superlatices[J]. Appl Phys Lett,1998,72(3):380-382.

    [107] JIANG J S, FULLERTON E E, GRIMSDITCH M, et al. Exchange-spring behavior in epitaxial hard/soft magnetic bilayer films[J]. J Appl Phys,1998,83(11):6238-6240.

    [108] FULLERTON E E, SOWERS C H, PEARSON J E, et al. A general approach to the epitaxial growth of rare-earth-transition-metal films[J]. Appl Phys Lett,1996,69(16):2438-2440.

    [109] ZHAO G P, DENG Y, ZHANG H W, et al. Thickness dependent magnetic reversal process and hysteresis loops in exchange-coupled hard-soft trilayers[J]. J Appl Phys,2010,108(9):093928-8.

    [110] LIU Y, GEORGE T A, SKOMSKI R, et al. Aligned and exchange-coupled FePt-based films[J]. Appl Phys Lett,2011,99(17):172504-3.

    [111] 彭懿,趙國平,吳紹全,等. 不同易軸取向下對(duì)Nd2Fe14B/Fe65Co35磁性雙層膜的微磁學(xué)模擬[J]. 物理學(xué)報(bào),2014,63(16):167505.

    [112] LI H L, LOU L, HOU F, et al. Simulataneously increasing the magnetization and coercivity of bulk nanocomposite magnets via severe plastic deformation[J]. Appl Phys Lett,2013,103(14):142406.

    (編輯李德華)

    Review of Micromagnetic Theory for Hysteresis Loop of Nanocomposite Hard/Soft Permanent Magnets

    ZHAO Guoping,WAN Xiulin

    (CollegeofPhysicsandElectronicEngineering,SichuanNormalUniversity,Chengdu610066,Sichuan)

    Abstract:Exchange-spring nanocomposite hard/soft magnets, as a new type of high-performance permanent magnets, has fundamental researches and applications and attracte much attention. This new kind of magnets with numerous magnetic properties and a giant especial magnetic energy product’has a theoretical upper limit of 1 MJ/m3. Micromagnetics as a bridge between the quantum magnetism and macroscopic magnetism, is an important theory to investigate the properties of magnetic materials and can clearly demonstrate the reversal mechanism of a given material as well as predicting its hysteresis loop and magnetic properties. This paper mainly expounds the hysteresis loops of thin film composite magnets according to three important micromagnetic calculation models: full analytical, quasi-analytical and numerical models. The advantages and weaknesses of each model are discussed. The results obtained from the three different models are compared and the effect of the hard/soft phase thicknesses and easy axis orientation on the macro hysteresis loops and coercivity mechanism are analysed. The gap between the theoretical and the experimental values are disscussed to propose some methods to enhance the magnetic properties of a given material, especially its coercivity and energy product.

    Key words:micromagnetic; hysteresis loop; magnetic energy product

    doi:10.3969/j.issn.1001-8395.2016.01.024

    中圖分類號(hào):O482.5

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1001-8395(2016)01-0136-15

    作者簡(jiǎn)介:趙國平(1967—),男,教授,主要從事納米永磁材料、斯格明子、交換偏置和鳥類磁感應(yīng)機(jī)制的研究,E-mail:zhaogp@uestc.edu.cn

    基金項(xiàng)目:國家自然科學(xué)基金(10747007和11074179)

    收稿日期:2015-11-10

    猜你喜歡
    磁滯回線
    基于MATLAB處理大學(xué)物理實(shí)驗(yàn)數(shù)據(jù)探究
    變壓器激磁涌流優(yōu)化設(shè)計(jì)及應(yīng)用探討
    磁場(chǎng)強(qiáng)度波形畸變對(duì)交流磁滯回線形狀的影響
    虛實(shí)結(jié)合的磁滯回線實(shí)驗(yàn)教學(xué)設(shè)計(jì)與實(shí)踐
    磁性材料磁滯回線重要參數(shù)的應(yīng)用分析
    振動(dòng)樣品磁強(qiáng)計(jì)的磁性表征測(cè)量
    高頻脈沖激勵(lì)下磁滯回線動(dòng)態(tài)測(cè)量裝置的設(shè)計(jì)及分析
    磁性材料磁化曲線和磁滯回線的Matlab繪制與擬合
    科技視界(2016年26期)2016-12-17 18:02:13
    磁通門鐵芯渦流效應(yīng)磁場(chǎng)計(jì)算與HSPICE仿真*
    基于磁滯回線數(shù)學(xué)模型的變壓器勵(lì)磁電流分析
    黑龍江電力(2015年2期)2015-03-06 04:32:26
    国产欧美日韩精品亚洲av| 一进一出抽搐动态| 亚洲成av人片免费观看| 成年免费大片在线观看| 蜜桃久久精品国产亚洲av| 国产色婷婷99| 国产探花极品一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲专区国产一区二区| 国产精品一区二区性色av| 午夜影院日韩av| 草草在线视频免费看| 男女做爰动态图高潮gif福利片| 亚洲综合色惰| 国产蜜桃级精品一区二区三区| 亚洲av美国av| 欧美一区二区国产精品久久精品| 久久久久久大精品| 成人午夜高清在线视频| 色播亚洲综合网| 久9热在线精品视频| 成年免费大片在线观看| 精品人妻一区二区三区麻豆 | 中文资源天堂在线| 国产精品人妻久久久久久| 两个人的视频大全免费| 精品午夜福利在线看| 在线观看午夜福利视频| 国产色爽女视频免费观看| 免费av不卡在线播放| eeuss影院久久| 免费看美女性在线毛片视频| 国产v大片淫在线免费观看| 国产伦在线观看视频一区| 日韩精品青青久久久久久| 久久精品影院6| 青草久久国产| 99热精品在线国产| 久久久久久九九精品二区国产| 国产精品一区二区三区四区久久| 午夜精品久久久久久毛片777| 成人美女网站在线观看视频| 久久精品人妻少妇| 欧美一区二区亚洲| 欧美中文日本在线观看视频| 人妻制服诱惑在线中文字幕| 久久精品国产亚洲av香蕉五月| 亚洲av熟女| 蜜桃亚洲精品一区二区三区| 免费看光身美女| 国产真实伦视频高清在线观看 | 久久久久九九精品影院| 亚洲人成电影免费在线| 在线免费观看不下载黄p国产 | 午夜老司机福利剧场| 欧美xxxx性猛交bbbb| 精品日产1卡2卡| 国产色婷婷99| 少妇熟女aⅴ在线视频| 国产午夜福利久久久久久| 亚洲精品在线观看二区| 一进一出好大好爽视频| 亚洲性夜色夜夜综合| 午夜老司机福利剧场| 国产成人av教育| 老司机福利观看| 日韩 亚洲 欧美在线| 久久99热6这里只有精品| 97热精品久久久久久| 直男gayav资源| 麻豆成人av在线观看| 99久久精品热视频| 日韩亚洲欧美综合| 久久久久免费精品人妻一区二区| 欧美日本亚洲视频在线播放| 国产欧美日韩一区二区三| 男女床上黄色一级片免费看| 欧美潮喷喷水| 又紧又爽又黄一区二区| 亚洲专区中文字幕在线| 国产精品,欧美在线| a级毛片免费高清观看在线播放| 亚洲精品久久国产高清桃花| 亚洲人成伊人成综合网2020| 国产综合懂色| 日日摸夜夜添夜夜添av毛片 | 免费在线观看日本一区| 国产高潮美女av| 欧美一区二区精品小视频在线| 亚洲乱码一区二区免费版| 国产一区二区在线观看日韩| 国产精品久久久久久精品电影| 狠狠狠狠99中文字幕| 亚洲欧美日韩高清专用| 高清在线国产一区| 欧美高清成人免费视频www| 国产精品一区二区三区四区免费观看 | 亚洲,欧美,日韩| 小说图片视频综合网站| 人妻制服诱惑在线中文字幕| 色综合婷婷激情| 欧美精品啪啪一区二区三区| 日本成人三级电影网站| 亚洲精品在线美女| 日韩精品青青久久久久久| 美女免费视频网站| 九九久久精品国产亚洲av麻豆| 免费av毛片视频| 国产一区二区激情短视频| 一夜夜www| 亚洲欧美清纯卡通| 免费大片18禁| 美女cb高潮喷水在线观看| 热99re8久久精品国产| 性欧美人与动物交配| 淫妇啪啪啪对白视频| 可以在线观看毛片的网站| 99久久99久久久精品蜜桃| 老司机深夜福利视频在线观看| 亚洲电影在线观看av| 人人妻,人人澡人人爽秒播| 久久久成人免费电影| 草草在线视频免费看| 97超视频在线观看视频| 久久久久久大精品| 老鸭窝网址在线观看| 久9热在线精品视频| 女人被狂操c到高潮| 中文字幕免费在线视频6| 亚洲在线自拍视频| 欧美日韩瑟瑟在线播放| 三级国产精品欧美在线观看| 久久欧美精品欧美久久欧美| 热99在线观看视频| www.999成人在线观看| 琪琪午夜伦伦电影理论片6080| 欧美日韩乱码在线| 午夜老司机福利剧场| 精品人妻熟女av久视频| 亚洲不卡免费看| a级一级毛片免费在线观看| 人妻制服诱惑在线中文字幕| 日本黄大片高清| 亚洲 国产 在线| 日本五十路高清| 久久久成人免费电影| 丰满乱子伦码专区| 18美女黄网站色大片免费观看| 禁无遮挡网站| 久久久久九九精品影院| 在线免费观看不下载黄p国产 | 免费黄网站久久成人精品 | 精品一区二区免费观看| 亚洲18禁久久av| 日本免费a在线| 精品久久久久久久末码| 久久99热6这里只有精品| 真实男女啪啪啪动态图| 老司机深夜福利视频在线观看| 免费无遮挡裸体视频| av天堂中文字幕网| 国产av在哪里看| 少妇丰满av| 露出奶头的视频| 国产亚洲av嫩草精品影院| 国产久久久一区二区三区| 国产伦一二天堂av在线观看| 欧美日韩瑟瑟在线播放| 亚洲一区二区三区色噜噜| 亚洲片人在线观看| 亚洲无线在线观看| 中亚洲国语对白在线视频| 真人一进一出gif抽搐免费| 精品99又大又爽又粗少妇毛片 | 国产aⅴ精品一区二区三区波| 人人妻人人看人人澡| 永久网站在线| 久久国产乱子免费精品| 色哟哟·www| www.www免费av| 亚洲人成网站在线播| 日日摸夜夜添夜夜添av毛片 | 国产白丝娇喘喷水9色精品| 久久欧美精品欧美久久欧美| 久久久久久久久大av| 国产精品自产拍在线观看55亚洲| 一本一本综合久久| 老司机福利观看| 看黄色毛片网站| 亚洲中文字幕一区二区三区有码在线看| 国产蜜桃级精品一区二区三区| 欧美激情在线99| 欧美丝袜亚洲另类 | 国产精品野战在线观看| 99视频精品全部免费 在线| 精品久久久久久久久亚洲 | 老鸭窝网址在线观看| 精品日产1卡2卡| 欧美乱色亚洲激情| 欧美日本亚洲视频在线播放| 亚洲美女黄片视频| 女生性感内裤真人,穿戴方法视频| 欧美色视频一区免费| 美女高潮的动态| 国产亚洲精品综合一区在线观看| 一本综合久久免费| 成人国产一区最新在线观看| a级毛片a级免费在线| 少妇人妻精品综合一区二区 | 国产在线男女| 国产一级毛片七仙女欲春2| 色噜噜av男人的天堂激情| 亚洲av第一区精品v没综合| xxxwww97欧美| 麻豆一二三区av精品| 1000部很黄的大片| 亚洲av免费在线观看| 毛片女人毛片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产三级在线视频| 成年女人永久免费观看视频| 男人和女人高潮做爰伦理| 看免费av毛片| 欧美日韩黄片免| 国产国拍精品亚洲av在线观看| 久9热在线精品视频| 亚洲精品乱码久久久v下载方式| 天天一区二区日本电影三级| 午夜福利欧美成人| 亚洲黑人精品在线| 校园春色视频在线观看| 婷婷精品国产亚洲av| 久久久成人免费电影| 午夜福利视频1000在线观看| 国产精品一区二区免费欧美| 亚洲专区国产一区二区| 国产伦精品一区二区三区视频9| 日本一二三区视频观看| 淫妇啪啪啪对白视频| 久久久久久九九精品二区国产| 九九热线精品视视频播放| 熟女人妻精品中文字幕| 91在线观看av| 搡老岳熟女国产| 免费搜索国产男女视频| 中文在线观看免费www的网站| 欧美高清成人免费视频www| 午夜福利18| 最近中文字幕高清免费大全6 | 日韩av在线大香蕉| 成人亚洲精品av一区二区| 亚洲色图av天堂| 久久亚洲精品不卡| 国产亚洲精品久久久com| 亚洲在线观看片| 久久热精品热| 国产精品女同一区二区软件 | 97超级碰碰碰精品色视频在线观看| 内射极品少妇av片p| 一区福利在线观看| 在线观看午夜福利视频| 夜夜躁狠狠躁天天躁| 亚洲精品亚洲一区二区| 99国产综合亚洲精品| 国产成人啪精品午夜网站| 亚洲电影在线观看av| 亚洲国产欧美人成| 国产亚洲欧美98| 亚洲色图av天堂| 我的女老师完整版在线观看| 97热精品久久久久久| 国产亚洲欧美98| 国内精品一区二区在线观看| 久久久久亚洲av毛片大全| 久久久久久久精品吃奶| 成年免费大片在线观看| 又黄又爽又免费观看的视频| 国产精品免费一区二区三区在线| 久9热在线精品视频| 搡老熟女国产l中国老女人| 搡女人真爽免费视频火全软件 | 赤兔流量卡办理| 久久人人精品亚洲av| 欧美成人a在线观看| 欧美最黄视频在线播放免费| 99国产极品粉嫩在线观看| 日韩av在线大香蕉| 久久久久久久亚洲中文字幕 | 国内精品久久久久久久电影| 成人永久免费在线观看视频| 精品人妻视频免费看| 免费搜索国产男女视频| 午夜精品一区二区三区免费看| 99视频精品全部免费 在线| 十八禁网站免费在线| 一进一出抽搐动态| 久久精品国产亚洲av涩爱 | 18禁在线播放成人免费| 国产亚洲欧美98| 亚洲经典国产精华液单 | 极品教师在线视频| 三级男女做爰猛烈吃奶摸视频| 久久久久久久久大av| 禁无遮挡网站| 韩国av一区二区三区四区| 欧美日本亚洲视频在线播放| 亚洲国产精品成人综合色| 亚洲精品乱码久久久v下载方式| 白带黄色成豆腐渣| 精品午夜福利在线看| 欧美丝袜亚洲另类 | 久久精品91蜜桃| 久久久久久大精品| 精品久久久久久成人av| 午夜两性在线视频| 亚洲无线在线观看| 亚洲久久久久久中文字幕| 麻豆久久精品国产亚洲av| 88av欧美| 国产成人aa在线观看| 久久久久免费精品人妻一区二区| 亚洲最大成人av| 女同久久另类99精品国产91| 五月伊人婷婷丁香| 最近中文字幕高清免费大全6 | 九色成人免费人妻av| 亚洲第一欧美日韩一区二区三区| 亚洲av二区三区四区| 欧美午夜高清在线| av女优亚洲男人天堂| 成人精品一区二区免费| 国产男靠女视频免费网站| 一个人观看的视频www高清免费观看| 久久国产精品人妻蜜桃| 亚洲成人久久性| 亚洲最大成人av| 精品午夜福利视频在线观看一区| 国产黄a三级三级三级人| 成年女人毛片免费观看观看9| 日本精品一区二区三区蜜桃| 精品久久久久久久久亚洲 | 精品久久久久久久久亚洲 | 免费人成在线观看视频色| 久久精品国产99精品国产亚洲性色| 日本撒尿小便嘘嘘汇集6| 久久久精品大字幕| 青草久久国产| 国产黄a三级三级三级人| 亚洲国产色片| 欧美一区二区亚洲| 久久久精品欧美日韩精品| 欧美乱妇无乱码| 亚洲最大成人手机在线| 男人的好看免费观看在线视频| 亚洲熟妇熟女久久| 欧美乱妇无乱码| 婷婷色综合大香蕉| 国产精品久久久久久人妻精品电影| 欧美xxxx黑人xx丫x性爽| 少妇的逼好多水| 人人妻,人人澡人人爽秒播| 日本一本二区三区精品| 人妻制服诱惑在线中文字幕| 又紧又爽又黄一区二区| 国产一区二区三区视频了| 一a级毛片在线观看| 天美传媒精品一区二区| 日韩有码中文字幕| 99热6这里只有精品| 成熟少妇高潮喷水视频| 国产精品电影一区二区三区| 欧美日韩综合久久久久久 | 久久天躁狠狠躁夜夜2o2o| 成年免费大片在线观看| 我要搜黄色片| 18禁裸乳无遮挡免费网站照片| 九色国产91popny在线| 欧美一区二区精品小视频在线| 国产高清视频在线观看网站| 岛国在线免费视频观看| 亚洲欧美精品综合久久99| 日韩中文字幕欧美一区二区| 久久久成人免费电影| 久久6这里有精品| 久久久色成人| 热99re8久久精品国产| 人妻丰满熟妇av一区二区三区| 欧美激情久久久久久爽电影| 国产男靠女视频免费网站| 日韩欧美三级三区| 久久久久国内视频| 色av中文字幕| 国产乱人伦免费视频| 国产高清视频在线观看网站| 久久精品人妻少妇| 欧美一区二区亚洲| 欧美色欧美亚洲另类二区| 国产探花极品一区二区| 亚洲激情在线av| 很黄的视频免费| 网址你懂的国产日韩在线| 国产黄a三级三级三级人| 欧美成人免费av一区二区三区| 亚洲专区国产一区二区| 色哟哟·www| 国产中年淑女户外野战色| 一个人看的www免费观看视频| 免费高清视频大片| 成年女人毛片免费观看观看9| av在线蜜桃| 淫秽高清视频在线观看| 中文字幕免费在线视频6| 精品人妻偷拍中文字幕| 露出奶头的视频| 丰满人妻熟妇乱又伦精品不卡| 一区二区三区高清视频在线| 免费看美女性在线毛片视频| 国产大屁股一区二区在线视频| ponron亚洲| 亚洲av电影在线进入| 精品人妻一区二区三区麻豆 | 黄色一级大片看看| 久久精品国产清高在天天线| 国产人妻一区二区三区在| 18+在线观看网站| 久久精品久久久久久噜噜老黄 | 亚洲国产精品sss在线观看| 狠狠狠狠99中文字幕| 中文在线观看免费www的网站| 12—13女人毛片做爰片一| 淫秽高清视频在线观看| 亚洲久久久久久中文字幕| 人妻夜夜爽99麻豆av| 国产成人福利小说| 国产精品永久免费网站| 日韩免费av在线播放| 少妇人妻一区二区三区视频| 久久久色成人| 波野结衣二区三区在线| 免费高清视频大片| 特大巨黑吊av在线直播| 91久久精品电影网| 90打野战视频偷拍视频| 日韩有码中文字幕| 国产伦一二天堂av在线观看| 国产一区二区亚洲精品在线观看| 美女大奶头视频| 天堂动漫精品| 久久精品国产亚洲av天美| 人妻夜夜爽99麻豆av| 性插视频无遮挡在线免费观看| 9191精品国产免费久久| 男女床上黄色一级片免费看| 免费黄网站久久成人精品 | 夜夜爽天天搞| netflix在线观看网站| 97超视频在线观看视频| 亚洲不卡免费看| 亚洲欧美日韩无卡精品| 欧美成狂野欧美在线观看| 国产乱人伦免费视频| 51午夜福利影视在线观看| 国产av麻豆久久久久久久| 国产私拍福利视频在线观看| 国产精品久久久久久久久免 | 丰满人妻一区二区三区视频av| 女人被狂操c到高潮| 午夜视频国产福利| 一级a爱片免费观看的视频| 国产成人aa在线观看| 久久久久国产精品人妻aⅴ院| 综合色av麻豆| 国产精品一区二区免费欧美| 精品国产三级普通话版| 日本撒尿小便嘘嘘汇集6| 久久久久久九九精品二区国产| 午夜老司机福利剧场| 在线十欧美十亚洲十日本专区| 国产精品日韩av在线免费观看| 一本精品99久久精品77| 久久婷婷人人爽人人干人人爱| 女人被狂操c到高潮| 神马国产精品三级电影在线观看| 久久午夜亚洲精品久久| 国产乱人伦免费视频| 欧美激情在线99| 中文资源天堂在线| 国产三级黄色录像| 亚洲欧美日韩高清在线视频| 亚洲美女视频黄频| 国产精品亚洲av一区麻豆| 欧美+亚洲+日韩+国产| 免费观看的影片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 日日干狠狠操夜夜爽| 韩国av一区二区三区四区| 亚洲,欧美精品.| 老女人水多毛片| 在线观看一区二区三区| 国产69精品久久久久777片| 色哟哟哟哟哟哟| 一区二区三区免费毛片| 国产毛片a区久久久久| 少妇裸体淫交视频免费看高清| 亚洲中文日韩欧美视频| 哪里可以看免费的av片| 国产又黄又爽又无遮挡在线| 国产不卡一卡二| 91麻豆av在线| 丁香欧美五月| 国产精品伦人一区二区| 色在线成人网| 麻豆国产av国片精品| 成年女人看的毛片在线观看| 中国美女看黄片| 一个人看的www免费观看视频| 日本a在线网址| 亚洲成av人片在线播放无| 欧美一区二区国产精品久久精品| 亚洲在线观看片| 国产精品一区二区三区四区久久| 黄色丝袜av网址大全| 一级a爱片免费观看的视频| 成人鲁丝片一二三区免费| 欧美日本亚洲视频在线播放| 亚洲成a人片在线一区二区| 久久精品人妻少妇| 51国产日韩欧美| 我的老师免费观看完整版| 一个人观看的视频www高清免费观看| 亚洲欧美日韩东京热| 长腿黑丝高跟| 国产探花在线观看一区二区| 久久精品影院6| 脱女人内裤的视频| 国产国拍精品亚洲av在线观看| 真实男女啪啪啪动态图| 99热只有精品国产| 久久久久精品国产欧美久久久| 色综合亚洲欧美另类图片| 黄色丝袜av网址大全| 亚洲av.av天堂| 又黄又爽又刺激的免费视频.| 日韩欧美精品v在线| 在线免费观看不下载黄p国产 | 91av网一区二区| 精品午夜福利视频在线观看一区| 日韩中字成人| 中文字幕免费在线视频6| 热99在线观看视频| 中文字幕免费在线视频6| 久久精品夜夜夜夜夜久久蜜豆| 校园春色视频在线观看| 九九久久精品国产亚洲av麻豆| 精品日产1卡2卡| 草草在线视频免费看| 中文字幕人妻熟人妻熟丝袜美| 亚洲综合色惰| 国产色婷婷99| 国产av不卡久久| 成人精品一区二区免费| 成人特级av手机在线观看| 在线播放无遮挡| 男女之事视频高清在线观看| 激情在线观看视频在线高清| 久久亚洲真实| 国产视频一区二区在线看| 午夜福利视频1000在线观看| 亚洲国产高清在线一区二区三| 91在线精品国自产拍蜜月| 可以在线观看毛片的网站| 午夜福利免费观看在线| 国产伦精品一区二区三区四那| 欧美性猛交黑人性爽| 99久久精品一区二区三区| 国产成年人精品一区二区| 特大巨黑吊av在线直播| 亚洲成人中文字幕在线播放| 精品久久久久久成人av| 婷婷精品国产亚洲av在线| 亚洲欧美日韩高清专用| 老熟妇仑乱视频hdxx| 夜夜夜夜夜久久久久| 18禁黄网站禁片免费观看直播| 99热这里只有精品一区| 黄片小视频在线播放| 精品一区二区三区av网在线观看| 少妇的逼好多水| 国产熟女xx| 51午夜福利影视在线观看| 网址你懂的国产日韩在线| 给我免费播放毛片高清在线观看| 国产免费一级a男人的天堂| 精品无人区乱码1区二区| 日韩成人在线观看一区二区三区| 色av中文字幕| 久久精品国产清高在天天线| 免费观看精品视频网站| 露出奶头的视频| 熟妇人妻久久中文字幕3abv| 亚洲五月天丁香| 日本黄大片高清| 免费在线观看亚洲国产| 精品久久久久久久久亚洲 | 日本一本二区三区精品| 国产一区二区在线av高清观看| 国产精品野战在线观看| 精品免费久久久久久久清纯| 亚洲av二区三区四区| 欧美色视频一区免费| av中文乱码字幕在线| 日韩欧美在线二视频| 亚洲熟妇中文字幕五十中出| 国产亚洲精品av在线| 毛片一级片免费看久久久久 |