圖案化介質(zhì)磁盤-磁頭氣體潤滑動(dòng)態(tài)特性分析

楊利花 楊 凱 徐郝亮 虞 烈

1.西安交通大學(xué)機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室,西安,7100492.陜西省先進(jìn)飛行器服役環(huán)境與控制重點(diǎn)實(shí)驗(yàn)室,西安,710049

圖案化介質(zhì)磁盤-磁頭氣體潤滑動(dòng)態(tài)特性分析

楊利花1,2楊 凱1徐郝亮1虞 烈1

1.西安交通大學(xué)機(jī)械結(jié)構(gòu)強(qiáng)度與振動(dòng)國家重點(diǎn)實(shí)驗(yàn)室,西安,7100492.陜西省先進(jìn)飛行器服役環(huán)境與控制重點(diǎn)實(shí)驗(yàn)室,西安,710049

在小擾動(dòng)條件下，提出應(yīng)用偏導(dǎo)數(shù)法在復(fù)數(shù)域內(nèi)求解可壓縮氣體潤滑的動(dòng)態(tài)Reynolds方程，以計(jì)算任意磁盤表面形狀的磁盤-磁頭氣膜潤滑的動(dòng)態(tài)剛度和阻尼系數(shù)。以圓柱形和立方形比特圖案化介質(zhì)(BPM)磁盤為研究對象，研究了磁單元尺寸參數(shù)對氣膜動(dòng)態(tài)系數(shù)的影響，并比較了兩種BPM對應(yīng)的動(dòng)態(tài)性能。數(shù)值結(jié)果表明，與光滑磁盤相比，圖案化介質(zhì)磁盤對應(yīng)的氣膜潤滑動(dòng)態(tài)剛度系數(shù)和阻尼系數(shù)大大減小,且隨磁單元高度增大，動(dòng)態(tài)系數(shù)逐漸減小。相比圓柱形BPM，立方形BPM對應(yīng)的氣膜具有較大動(dòng)態(tài)剛度和阻尼系數(shù)，計(jì)算結(jié)果可為圖案化介質(zhì)磁盤的設(shè)計(jì)提供重要的理論依據(jù)。

比特圖案化介質(zhì)；計(jì)算機(jī)磁盤；存儲密度；動(dòng)態(tài)性能

0 引言

隨著磁盤存儲密度的不斷增大，磁頭與磁盤間隙(即飛行高度)越來越小，可能發(fā)生的磁頭磁盤界面接觸和凸峰沖擊等都會(huì)影響磁頭的飛行穩(wěn)定性和讀寫可靠性,因此，精確估計(jì)磁頭超薄氣膜潤滑的動(dòng)態(tài)特性對提高磁盤存儲密度和讀寫可靠性至關(guān)重要。

分析磁頭氣膜的動(dòng)態(tài)特性需要求解超薄氣體潤滑Reynolds方程和磁頭系統(tǒng)的運(yùn)動(dòng)方程。HONCHI等[1]建立了包含簡化懸臂支架模型、氣膜潤滑模型和磁頭-磁盤接觸模型的系統(tǒng)動(dòng)力學(xué)分析模型，采用有限差分法同步求解氣體潤滑Reynolds方程和磁頭的運(yùn)動(dòng)方程，研究了磁頭-磁盤接觸過程中的磁頭振動(dòng)模態(tài)。ZENG等[2-3]應(yīng)用開發(fā)的動(dòng)態(tài)求解器求解小擾動(dòng)下磁頭浮動(dòng)塊的振動(dòng)響應(yīng)，后又通過模態(tài)分析法獲得磁頭-磁盤系統(tǒng)的模態(tài)參數(shù)以及剛度和阻尼系數(shù)。ONO[4]通過小擾動(dòng)法和有限差分法在頻域計(jì)算了氣膜的動(dòng)態(tài)壓力分布，分析了不同表面形貌浮動(dòng)塊的氣膜動(dòng)態(tài)系數(shù)及動(dòng)態(tài)響應(yīng)，重點(diǎn)研究了擾動(dòng)頻率對氣膜動(dòng)態(tài)剛度和阻尼的影響。MITSUYA等[5]推導(dǎo)了與粗糙度方向相關(guān)的氣膜動(dòng)態(tài)潤滑方程，并用有限元法求解了動(dòng)態(tài)潤滑方程,得到了氣膜動(dòng)態(tài)剛度和阻尼系數(shù)，研究了粗糙度方向和大小、擾動(dòng)頻率、磁盤速度和載荷對動(dòng)態(tài)系數(shù)的影響。LI等[6]在小擾動(dòng)假設(shè)下用頻率響應(yīng)法計(jì)算了三自由度磁頭系統(tǒng)的氣膜剛度和阻尼系數(shù)，分析了不同浮動(dòng)塊表面結(jié)構(gòu)對應(yīng)的氣膜動(dòng)態(tài)系數(shù)隨擾動(dòng)頻率的變化關(guān)系，提出改進(jìn)磁頭阻尼的設(shè)計(jì)方案。莊蘋等[7]應(yīng)用伽遼金法和有限元法求解了動(dòng)態(tài)Reynolds方程，并進(jìn)行了磁頭-磁盤系統(tǒng)的動(dòng)態(tài)響應(yīng)分析，研究了粗糙度方向?qū)Υ蓬^飛行穩(wěn)定性的影響。王玉娟等[8]利用攝動(dòng)法建立氣體潤滑的靜動(dòng)態(tài)Reynolds方程，采用交換方向法研究了氣膜動(dòng)態(tài)剛度和阻尼與最小膜厚、磁盤轉(zhuǎn)速和磁頭傾角的關(guān)系。 白少先等[9]采用攝動(dòng)法進(jìn)行氣體潤滑頻域響應(yīng)分析，討論了氣膜剛度和阻尼等隨飛行高度和擾動(dòng)頻率的變化規(guī)律，可用于指導(dǎo)磁頭/磁盤系統(tǒng)的設(shè)計(jì)。

以上關(guān)于磁頭/磁盤系統(tǒng)氣膜動(dòng)態(tài)性能的研究，均是針對飛行在連續(xù)或光滑表面磁盤上的磁頭開展的。隨著更高磁盤存儲密度的需求，西部數(shù)據(jù)公司(western digital company, WD)等磁盤制造商近年來一直致力于比特圖案化介質(zhì)(bit pattern media, BPM)磁盤的制備與研發(fā)[10]。這種磁盤表面不再是光滑的連續(xù)表面，而是由制備在其上的多個(gè)獨(dú)立的有序圖案化磁單元排列構(gòu)成的離散表面，磁盤表面與磁頭上的浮動(dòng)塊表面間形成動(dòng)壓氣膜，氣膜的靜態(tài)和動(dòng)態(tài)性能對磁盤的存儲性能有重要影響。LI等[11-12]分析了圖案化介質(zhì)磁盤形狀和參數(shù)對浮動(dòng)塊飛行高度、氣膜壓力等靜態(tài)性能的影響。張延瑞等[13]以三維凸臺氣膜潤滑模型模擬圖案化介質(zhì)磁盤表面，研究了氣膜壓力與凸臺高度和最小間隙的關(guān)系，以討論圖案化介質(zhì)磁盤-磁頭系統(tǒng)的氣體潤滑靜態(tài)性能。LI等[14]應(yīng)用有限元分析方法研究了浮動(dòng)塊和圖案化介質(zhì)磁盤界面間的接觸狀態(tài)、材料特性以及不同充填材料下圖案化介質(zhì)的平坦化對磁盤塑性變形和溫度的影響。此類磁盤因具有超高的存儲密度和較好的讀寫穩(wěn)定性而成為近年來的發(fā)展目標(biāo)，但由于其制備的困難以及高成本而限制了其應(yīng)用。

現(xiàn)有文獻(xiàn)的研究方法多數(shù)是在頻域內(nèi)計(jì)算非BPM磁盤表面和磁頭浮動(dòng)塊表面間的氣膜動(dòng)態(tài)剛度和阻尼，以避免對與時(shí)間相關(guān)的氣膜潤滑動(dòng)態(tài)Reynolds方程的直接求解，而本文擬在小擾動(dòng)假設(shè)下，提出采用偏導(dǎo)數(shù)法在復(fù)數(shù)域內(nèi)直接求解氣體潤滑的動(dòng)態(tài)Reynolds方程，避免了求解中的時(shí)間相關(guān)性，以計(jì)算氣膜潤滑的動(dòng)態(tài)剛度和阻尼系數(shù)。偏導(dǎo)數(shù)法可適用于任何表面形貌的磁盤和磁頭間氣膜潤滑動(dòng)態(tài)剛度和阻尼的計(jì)算。本文分別以圓柱形和立方形圖案化磁盤為研究對象，給出了對應(yīng)計(jì)算結(jié)果,并研究了磁單元形狀和尺寸參數(shù)對圖案化介質(zhì)磁盤-磁頭的氣膜動(dòng)態(tài)特性的影響。

1 數(shù)學(xué)模型

1.1超薄氣體潤滑方程

超高密度磁記錄系統(tǒng)中，磁頭和磁盤間的氣膜間隙一般小于氣體分子的平均自由程的100倍，需要引入修正系數(shù)以計(jì)入氣體分子滑移效應(yīng)對氣膜壓力的影響，DUWENSEE[15]給出了三種修正系數(shù)計(jì)算公式,其中一階修正的修正系數(shù)為

Q1=1+6Kn

二階修正的修正系數(shù)為

Boltzmann修正的修正系數(shù)為

(1)

Kn=λ/h

式中，Kn為表征氣體稀薄程度的Knudsen數(shù)；λ、h分別為氣體分子的平均自由程和氣膜厚度。

本文計(jì)算采用 Boltzmann修正系數(shù)。建立圖1所示三自由度坐標(biāo)系，引入修正系數(shù)后的超薄氣膜潤滑量綱一Reynolds方程形式為

(2)

式中，p為有量綱氣膜壓力;h為有量綱氣膜厚度;P、H分別為量綱一氣膜壓力和厚度;hmin為最小氣膜間隙或飛行高度；pa為環(huán)境大氣壓力；X、Y分別為浮動(dòng)塊長度和寬度方向的量綱一坐標(biāo);L為磁頭浮動(dòng)塊長度；x、y為有量綱坐標(biāo)；Λx、Λy為軸承數(shù)；μ為氣體動(dòng)力黏度；U、V分別為浮動(dòng)塊相對于磁盤在長度和寬度方向的飛行速度(本文計(jì)算中U為對應(yīng)磁盤徑向位置的線速度，V=0)；T為量綱一時(shí)間;ω為磁盤角速度；t為時(shí)間；σ為壓縮數(shù)。

磁頭浮動(dòng)塊氣膜潤滑量綱一氣膜厚度為

圖1 磁頭浮動(dòng)塊坐標(biāo)系Fig.1 Coordinate system of slider

(3)

Z=z/hminXp=xp/LYp=yp/L

式中，Z為浮動(dòng)塊垂直方向的量綱一坐標(biāo)；z為有量綱縱向坐標(biāo)；α、β分別為俯仰角(pitch角)和翻轉(zhuǎn)角(roll角)；Xp、Yp為樞軸點(diǎn)的量綱一坐標(biāo)；xp、yp為樞軸點(diǎn)有量綱坐標(biāo)。

應(yīng)用有限差分法和牛頓-拉普森迭代法同時(shí)迭代求解超薄氣膜潤滑的靜態(tài)Reynolds方程和浮動(dòng)塊的力和力矩平衡方程，就可獲得浮動(dòng)塊在特定工況時(shí)的靜態(tài)平衡位置及氣膜壓力和氣膜厚度分布[16]。

1.2氣膜動(dòng)態(tài)性能計(jì)算[17]

采用小擾動(dòng)法，假設(shè)浮動(dòng)塊繞其靜態(tài)平衡位置發(fā)生周期性小擾動(dòng)，任意時(shí)刻的浮動(dòng)塊位置可以表示為

(4)

Ω=ν/ω

同樣利用小擾動(dòng)假設(shè)，可將氣膜壓力P和氣膜厚度H表示為靜態(tài)項(xiàng)和動(dòng)態(tài)項(xiàng)，其形式為

(5)

(6)

將式(4)代入式(5)，并忽略無窮小項(xiàng)，得到任意時(shí)刻的氣膜厚度為

(7)

則動(dòng)態(tài)氣膜厚度擾動(dòng)幅值為

(8)

將式(5)和式(6)代入動(dòng)態(tài)Reynolds方程(式(2))，利用靜態(tài)Reynolds方程化簡并忽略高階無窮小量，得到小擾動(dòng)時(shí)超薄氣體潤滑的動(dòng)態(tài)Reynolds方程如下：

(9)

關(guān)于QZ和HZ的方程為

(10)

HZ=1

(11)

由于QZ和HZ為復(fù)變量，可令QZ=QZR+iQZI，HZ=HZR+iHZI，代入式(10)，則可得到關(guān)于QZR、QZI的實(shí)部和虛部方程，并將它們化成如下橢圓型偏微分方程-·(cu)+au=f形式：

(12)

(13)

應(yīng)用MATLAB的偏微分方程PDE工具箱同時(shí)迭代求解式(12)和式(13)， 就可得到QZ的解。

同理，可得到Qα和Qβ的偏微分方程并求得其結(jié)果，其中Hα和Hβ有如下形式：

(14)

(15)

求得QZ、Qα和Qβ后, 則浮動(dòng)塊表面氣膜的量綱一動(dòng)態(tài)剛度和阻尼系數(shù)為

(16)

m=z,α,β

對應(yīng)的動(dòng)態(tài)系數(shù)為

[8][32][48] 范子英、彭飛、劉沖：《政治關(guān)聯(lián)與經(jīng)濟(jì)增長：基于衛(wèi)星燈光數(shù)據(jù)的研究》，《經(jīng)濟(jì)研究》2016年第1期，第114-126頁。

對三自由度浮動(dòng)塊氣膜潤滑而言，應(yīng)用小擾動(dòng)法可得到其9個(gè)剛度系數(shù)kij(i,j=z，α，β)和9個(gè)阻尼系數(shù)dij，由于偏導(dǎo)數(shù)法消除了Reynolds方程中的時(shí)間項(xiàng)，將問題轉(zhuǎn)換到復(fù)數(shù)域進(jìn)行求解，從而大大降低了求解的復(fù)雜性和速度。

應(yīng)用這些動(dòng)態(tài)剛度和阻尼系數(shù)，可以得到氣膜潤滑的動(dòng)態(tài)氣膜力和力矩矩陣表達(dá)形式如下:

(17)

2 結(jié)果與分析

2.1氣膜動(dòng)態(tài)性能計(jì)算

圖2所示為磁頭上的浮動(dòng)塊幾何模型及表面形貌，可以看出浮動(dòng)塊工作表面并非光滑。浮動(dòng)塊長度L=850 μm，寬度B=700 μm，磁頭和浮動(dòng)塊連接的樞軸點(diǎn)位于浮動(dòng)塊幾何中心位置。圖3給出了磁盤上兩種不同圖案化磁單元形狀的結(jié)構(gòu)模型及參數(shù)，即圓柱形和立方形圖案化介質(zhì)BPM，以研究磁單元形狀對氣膜動(dòng)態(tài)特性的影響。其中，b為相鄰磁單元邊距，d為圓柱形磁單元直徑或立方形磁單元邊長，hb為磁單元高度，w為相鄰磁單元間距,w=b-d。磁盤工作過程中，圖2浮動(dòng)塊工作表面和圖3的BPM表面間會(huì)形成動(dòng)壓氣膜潤滑。

圖2 浮動(dòng)塊幾何模型及表面形貌Fig.2 Geometrical model and surface topography of slider

(a)圓柱形BPM

(b)立方形BPM圖3 BPM磁盤模型Fig.3 BPM model

本文所研究的磁盤-磁頭浮動(dòng)塊表面間的潤滑氣體采用含有98%氦氣的空氣-氦氣混合工質(zhì)，動(dòng)力黏度μ=2.018×10-5kg/(m·s)，密度ρ=0.1801 kg/m3，環(huán)境壓力pa=65.86 kPa, 氣體分子平均自由程λ=190 nm。

2.2BPM尺寸對氣膜潤滑動(dòng)態(tài)性能的影響

本文所研究的磁盤-磁頭在浮動(dòng)塊長度方向線速度U=20 m/s，通過有限元分析獲得懸臂支架剛度，分別為kz=20 N/m，kα=kβ=3.5×10-5N·m/rad。樞軸點(diǎn)坐標(biāo)位置為xp=425 μm，yp=350 μm時(shí)，通過求解靜態(tài)Reynolds方程和浮動(dòng)塊的力平衡方程計(jì)算了空氣-氦氣混合工質(zhì)下圖案化介質(zhì)磁盤對應(yīng)的浮動(dòng)塊表面靜態(tài)壓力和厚度分布，其后在小擾動(dòng)假設(shè)下應(yīng)用上述偏導(dǎo)數(shù)法求解了不同磁單元尺寸對應(yīng)的氣膜潤滑動(dòng)態(tài)系數(shù)。由于動(dòng)態(tài)系數(shù)的交叉項(xiàng)相對于直接項(xiàng)很小，故本文只討論氣膜潤滑的直接動(dòng)態(tài)系數(shù)。

圖4是不同磁單元尺寸的立方形圖案化介質(zhì)BPM對應(yīng)的氣膜動(dòng)態(tài)剛度系數(shù)和阻尼系數(shù)隨磁單元高度的變化規(guī)律曲線。從圖4可以看出，由于氣膜厚度方向(Z向)為主要承載方向，故其對應(yīng)剛度系數(shù)和阻尼系數(shù)遠(yuǎn)大于其他兩個(gè)方向的剛度系數(shù)和阻尼系數(shù)。與磁單元高度為零的光滑磁盤相比，圖案化介質(zhì)BPM使氣膜潤滑的動(dòng)態(tài)系數(shù)大大降低，其中kzz降低了大約6倍。隨磁單元高度增大，所有的剛度系數(shù)和阻尼系數(shù)均呈微小減小趨勢。對于具有相同磁單元尺寸比b/w=2的兩種參數(shù)b=10 μm,w=5 μm和b=1.2 μm,w=0.6 μm而言，它們對應(yīng)的動(dòng)態(tài)系數(shù)也幾乎相同。而較大的b/w值則對應(yīng)較大的動(dòng)態(tài)剛度系數(shù)kzz、kαα、kββ和阻尼系數(shù)dββ，特別對較大b/w值的b=1.2 μm,w=0.3 μm而言，其剛度系數(shù)和阻尼系數(shù)明顯比其他參數(shù)時(shí)大。這些結(jié)果可為圖案化介質(zhì)磁盤的設(shè)計(jì)提供理論參考。而且結(jié)果表明，盡管圖案化介質(zhì)磁盤能夠提高磁盤的存儲密度，但它會(huì)降低氣膜潤滑的動(dòng)態(tài)性能。這是因?yàn)楦?dòng)塊表面的氣膜壓力主要分布在它與圖案化磁單元表面對應(yīng)區(qū)域，而與磁單元間隙對應(yīng)區(qū)域的氣膜壓力很小[16]。因此與連續(xù)光滑表面磁盤相比，圖案化介質(zhì)磁盤因?qū)?yīng)的有效氣膜壓力作用面積減小，使浮動(dòng)塊表面的氣膜合力減小，從而降低了對應(yīng)氣膜潤滑的動(dòng)態(tài)性能。且磁單元高度增大不會(huì)影響有效氣膜壓力作用面積，只會(huì)使與磁單元間隙對應(yīng)區(qū)域的氣膜壓力有微小減小，因此磁單元高度變化對氣膜動(dòng)態(tài)性能的影響很小。

圖5是磁單元尺寸為b=1.2 μm,w=0.3 μm時(shí)，圓柱形和立方形BPM對應(yīng)的氣膜潤滑動(dòng)態(tài)系數(shù)隨磁單元高度的變化曲線。結(jié)果表明，在此磁單元尺寸下，隨磁單元高度增大，圓柱形和立方形BPM對應(yīng)的氣膜潤滑的6個(gè)直接剛度系數(shù)和阻尼系數(shù)均逐漸減小。與圓柱形BPM相比，立方形圖案化介質(zhì)對應(yīng)的氣膜具有較大的剛度和阻尼系數(shù)。也就是說，立方形BMP對應(yīng)的浮動(dòng)塊具有比圓柱形BMP更優(yōu)越的動(dòng)態(tài)性能。這是由于當(dāng)圓柱形磁單元直徑與立方形磁單元邊長相同時(shí)，立方形磁單元對應(yīng)的有效氣膜壓力作用面積大于圓柱形磁單元對應(yīng)的有效氣膜壓力作用面積，因而立方形磁單元對應(yīng)的浮動(dòng)塊表面氣膜合力較大,所以立方體磁單元的動(dòng)態(tài)性能更優(yōu)。

(a)剛度系數(shù)kzz (b)剛度系數(shù)kαα (c)剛度系數(shù)kββ

(d)阻尼系數(shù)dzz (e)阻尼系數(shù)dαα (f)阻尼系數(shù)dββ

圖4 立方形BPM對應(yīng)氣膜潤滑動(dòng)態(tài)性能隨磁單元高度變化
Fig.4 Changes of dynamic coefficients of air bearing over cuboidal bit pattern media with bit height

(a)剛度系數(shù)kzz (b)剛度系數(shù)kαα (c)剛度系數(shù)kββ

(d)阻尼系數(shù)dzz (e)阻尼系數(shù)dαα (f)阻尼系數(shù)dββ

圖5 立方形和圓柱形BPM對應(yīng)浮動(dòng)塊動(dòng)態(tài)性能比較
Fig.5 Comparison of dynamic coefficients for air bearing over cylindrical and cuboidal bit pattern media

3 結(jié)論

(1)在小擾動(dòng)假設(shè)下提出采用偏導(dǎo)數(shù)法在復(fù)數(shù)域內(nèi)求解浮動(dòng)塊的氣膜潤滑動(dòng)態(tài)Reynolds方程，避免了方程中時(shí)間項(xiàng)的求解。此方法可用于具有任意表面形狀的磁盤-磁頭氣膜潤滑動(dòng)特性分析，具有普遍適用性。應(yīng)用此方法分別計(jì)算了圓柱形和立方形圖案化磁盤介質(zhì)對應(yīng)的氣膜潤滑動(dòng)態(tài)系數(shù)，研究了磁單元結(jié)構(gòu)和尺寸參數(shù)對動(dòng)態(tài)系數(shù)的影響。

(2) 相比光滑磁盤表面，無論立方形還是圓柱形圖案化介質(zhì)，因其對應(yīng)的有效氣膜壓力作用面積減小使氣膜合力降低，從而造成氣膜潤滑剛度系數(shù)和阻尼系數(shù)均大大減小，且隨磁單元高度增大稍有減小。

(3)因磁單元尺寸相同時(shí)，立方形圖案化介質(zhì)對應(yīng)的有效氣膜壓力作用面積大于圓柱形圖案化介質(zhì)對應(yīng)的有效氣膜壓力作用面積，因此其對應(yīng)的氣膜具有比圓柱形 BPM更大的剛度系數(shù)和阻尼系數(shù)，表明立方形BPM比圓柱形BPM具有更好的動(dòng)態(tài)性能，可為圖案化介質(zhì)磁盤的設(shè)計(jì)提供重要的理論依據(jù)。

致謝本文研究工作是在美國加州大學(xué)圣地亞哥分校(University of California, San Diego)磁記錄研究中心(Center for Memory and Recording Research)Frank E. Talke教授的精心指導(dǎo)下和研究生的大力幫助下完成的，在此表示衷心感謝!

[1] HONCHI M, KOHIRA H, MATSUMOTO M. Numerical Simulation of Slider Dynamics during Slider-disk Contact[J]. Tribolology International,2003,36:235-240.

[2] ZENG Qinghua, CHEN Liangsheng, BOGY D B. A Modal Analysis Method for Slider Air Bearing in Hard Disk Drives[J]. IEEE Transactions on Magnetics,1997,33(5):3124-3126.

[3] ZENG Qinghua, BOGY D B. Stiffness and Damping Evaluation of Air Bearing Sliders and New Designs with High Damping[J]. Journal of Tribology,1999,121(2):341-347.

[4] ONO K. Dynamic Characteristics of Air-lubricated Slider Bearing for Noncontact Magnetic Recording[J]. Journal of Lubrication Technology,1975,97(2):250-260.

[5] MITSUYA Y, OTA H. Stiffness and Damping of Compressible Lubricating Films between Computer Flying Heads and Textured Media: Perturbation Analysis Using the Finite Element Method[J]. Journal of Tribology,1991,113(4):819-827.

[6] LI Jianhua, XU Junguo, SHIMIZU Y, et al. Design and Evaluation of Damped Air Bearings at Head-disk Interface[J]. Journal of Tribology,2010,32(3):031702-1-13.

[7] 莊蘋，王玉娟，陳云飛. 超薄氣體潤滑磁盤/磁頭系統(tǒng)動(dòng)力學(xué)分析[J]. 中國機(jī)械工程，2001,12(10): 1191-1194. ZHUANG Ping, WANG Yujuan, CHEN Yunfei. Ro-ughness Effects on Dynamic Response of High Density Disk Drives Under High Knudsen Number[J]. China Mechanical Engineering,2001,12(10):1191-1194.

[8] 王玉娟，陳云飛，莊蘋，等. 磁頭磁盤系統(tǒng)動(dòng)特性參數(shù)及系統(tǒng)穩(wěn)定性分析[J]. 中國機(jī)械工程，2002,13(3): 256-259. WANG Yujuan, CHEN Yunfei, ZHUANG Ping, et al. Dynamic Coefficients and Dynamic Stability in Magnetic Recording[J]. China Mechanical Engineering,2002,13(3):256-259.

[9] 白少先，彭旭東，孟永剛，等. 磁頭/磁盤超薄膜氣體潤滑動(dòng)態(tài)性能分析[J]. 應(yīng)用力學(xué)學(xué)報(bào)， 2009，26(4): 685-689. BAI Shaoxian, PENG Xudong, MENG Yonggang, et al. Dynamic Character of Ultra Thin Film Gas Lubrication for Magnetic Head/Disk[J]. Chinese Journal of Applied Mechanics,2009,26(4):685-689.

[10] ALBRECHT T R, ARORA H, AYANOOR-VITIKKATE V, et al. Bit Patterned Magnetic Recording: Theory, Media Fabrication, and Recording Performance[J]. IEEE Transactions on Magnetics,2015,51(5):1-42.

[11] LI Hui, ZHENG Hao, YOON Y, et al. Air Bearing Simulation for Bit Patterned Media[J]. Tribology Letters,2009,33:199-204.

[12] LI Hui, TALKE F E. Numerical Simulation of the Head/Disk Interface for Bit Patterned Media[J]. IEEE Transactions on Magnetics,2009,45(11):4984-4989.

[13] 張延瑞，孟永剛. 基于連續(xù)模型和離散模型的微尺度三維凸臺軸承的氣體潤滑分析[J]. 摩擦學(xué)學(xué)報(bào)，2011，31(4): 397-407. ZHANG Yanrui, MENG Yonggang. Gas Lubrication Analysis of Three Dimensional Micro Scale Gas Slider Bearings with Bumps Based on Continuum and Discrete Models [J]. Tribology,2011,31(4):397-407.

[14] LI Longqiu, SONG Wenping, ZHANG Chunhui, et al. Investigation of Thermo-mechanical Contact between Slider and Bit Patterned Media[J]. Microsystem Technologies,2012,18(9/10):1567-1574.

[15] DUWENSEE M. Numerical and Experimental Investigation of the Head Disk Interface [D]. Q. San Diego: University of California, 2007.

[16] 楊利花，徐郝亮，虞烈. 比特圖案化介質(zhì)磁盤-磁頭的氣膜靜態(tài)特性研究 [J]. 西安交通大學(xué)學(xué)報(bào)，2016, 50(6): 110-115. YANG Lihua, XU Haoliang, Yu Lie. Study on Static Characteristics of Air Bearing Film over Bit Pattern Media in Hard Disk Drives[J]. Journal of Xi’an Jiaotong University,2016,50(6):110-115.

[17] YANG Lihua, QI Shemiao, YU Lie. Analysis on Dynamic Performance of Hydrodynamic Tilting-pad Gas Bearings Using Partial Derivative Method [J]. Transactions of the ASME Journal of Tribology, 2009,131(1):011703-1-8.

(編輯王艷麗)

AnalysisonDynamicCharacteristicsofAirBearingoverBitPatternMediainHardDiskDrives

YANG Lihua1,2YANG Kai1XU Haoliang1YU Lie1

1.State Key Laboratory for Strength and Vibration of Mechanical Structures,
Xi’an Jiaotong University, Xi’an, 710049 2.Shaanxi Key Laboratory of Environment and Control for Flight Vehicle, Xi’an, 710049

With the assumption of small perturbation, this paper applied the partial derivative method to solve the compressible gas-lubricated dynamic Reynolds equation in the complex domain. This method might be used to calculate air-film dynamic stiffnesses and dampings of hard disk drives with arbitrary surface topography of disks. The dynamic stiffness and damping coefficients of air bearing over bit pattern media were calculated. Moreover, the effects of bit height of cylindrical and cuboidal bit pattern media on the dynamic coefficients were studied and also the dynamic performances of two types of bit pattern media were compared. The results indicate that the dynamic stiffness and damping coefficients of air bearing over bit pattern media largely decrease compared with those over smooth disk. With the increasing of bit height, the dynamic coefficients gradually reduce. Compared with cylindrical bit pattern media, the gas bearing over cuboidal bit pattern media has larger dynamic stiffnesses and dampings. These results may provide important theoretical reference for the designs of bit pattern media.

bit pattern media; computer disk; storage density; dynamic characteristics

2017-01-23

國家自然科學(xué)基金資助項(xiàng)目(51575425)

TH117.2

10.3969/j.issn.1004-132X.2017.17.004

楊利花，女，1975年生。西安交通大學(xué)航天航空學(xué)院副教授。主要研究方向?yàn)槟Σ翆W(xué)、軸承-轉(zhuǎn)子動(dòng)力學(xué)。發(fā)表論文50余篇。E-mail:yanglihua_2@126.com。楊凱，男，1989年生。西安交通大學(xué)航天航空學(xué)院碩士研究生。徐郝亮，男，1988年生。西安交通大學(xué)航天航空學(xué)院博士研究生。虞烈，男，1948年生。西安交通大學(xué)機(jī)械工程學(xué)院教授。