水輪發(fā)電機(jī)組推力軸承散熱系統(tǒng)研究綜述

張 楊

（中國電建集團(tuán)西北勘測設(shè)計研究院有限公司，陜西 西安 710065）

0 引言

隨著現(xiàn)代水力發(fā)電技術(shù)的不斷進(jìn)步，大型與巨型水輪發(fā)電機(jī)組的應(yīng)用不斷增多，對電網(wǎng)的穩(wěn)定、安全、可靠運(yùn)行的影響日益增大，要求機(jī)組的運(yùn)行穩(wěn)定性必須得到保證。水輪發(fā)電機(jī)組的推力軸承是保證機(jī)組安全穩(wěn)定運(yùn)行的重要部件，承擔(dān)的載荷往往比較大，在實際生產(chǎn)運(yùn)行過程中，因為機(jī)組偏心、軸承散熱系統(tǒng)油路循環(huán)不合理、潤滑油油質(zhì)劣化等因素引起軸承系統(tǒng)散熱不良、造成軸承瓦溫過高甚至燒損的現(xiàn)象時有發(fā)生，不僅直接影響了推力軸承的使用壽命，還給電力安全生產(chǎn)帶來了嚴(yán)重的威脅。

葛洲壩二江電廠1號機(jī)組在并網(wǎng)送電時曾發(fā)生推力瓦溫過高的情況，造成1號機(jī)開關(guān)自動斷開，機(jī)組與電網(wǎng)解列，停機(jī)后對推力軸瓦進(jìn)行檢查，發(fā)現(xiàn)瓦的徑向中部造成了嚴(yán)重?zé)龘p，根據(jù)巡檢記錄，有14塊瓦的溫度越限[1]；2015年云南梨園發(fā)電公司新投運(yùn)的型號為SF600-64/16880的水輪發(fā)電機(jī)推力軸承經(jīng)過長時間滿負(fù)荷運(yùn)行，推力軸承的瓦溫超出了技術(shù)合同允許值，油溫超出了報警值，嚴(yán)重影響了機(jī)組的供電可靠性[2]；三峽左岸電站6臺VGS機(jī)組推力軸承均為外循環(huán)冷卻方式，在夏季高溫時，即使推力軸承油冷卻器全部投入運(yùn)行，也會出現(xiàn)達(dá)83℃的高溫，通過仿真計算發(fā)現(xiàn)，冷卻器的容量未達(dá)到設(shè)計要求，2003年對推力軸承油冷卻器的結(jié)構(gòu)進(jìn)行了改造，采用針刺管型冷卻器并增加了冷卻水管個數(shù)，才使得平均油溫和瓦溫得到了降低[3]。實踐證明，大型水輪發(fā)電機(jī)推力軸承的可靠性最低，它的故障占水力機(jī)械總故障率的55%～70%[4]。國家電網(wǎng)2012年開始執(zhí)行的新的《國家電網(wǎng)公司安全事故調(diào)查規(guī)程》中，也把發(fā)電廠發(fā)電機(jī)組軸瓦燒損定為六級設(shè)備事件[5]。因此，如何提高推力軸承散熱系統(tǒng)的散熱性能，顯得日益重要。

1 推力軸承的潤滑冷卻方式

水輪發(fā)電機(jī)組推力軸承的潤滑冷卻可分為內(nèi)循環(huán)冷卻方式和外循環(huán)冷卻方式兩種。內(nèi)循環(huán)冷卻系統(tǒng)。推力軸承和油冷卻器安裝在同一油槽內(nèi)，潤滑油依靠軸承轉(zhuǎn)動部件的旋轉(zhuǎn)在軸承和冷卻器之間流動，由冷卻器中的冷卻水把軸承損耗在潤滑油中的熱量帶走，保證推力軸承在熱平衡的狀態(tài)下，油溫在規(guī)定的溫度下運(yùn)行[6]。利用油槽內(nèi)油的粘滯作用和對流換熱形成回路。

外循環(huán)方式中，推力軸承安放在油槽內(nèi)部，油冷卻器安放在油槽的外部，外循環(huán)動力的方式又分為自身泵與外加泵[7]。我國采用內(nèi)循環(huán)冷卻的機(jī)組有葛洲壩、巖灘、五強(qiáng)溪等水電站，采用外循環(huán)冷卻的機(jī)組有三峽、寶泉、白蓮河等水電站。

2 推力軸承國外研究概況

從19世紀(jì)末水輪發(fā)電機(jī)出現(xiàn)時起，推力軸承的承載能力曾一度制約著水電機(jī)組單機(jī)容量的增長，直到20世紀(jì)初期，分塊瓦的Kingsbury推力軸承的出現(xiàn)，才使得水輪發(fā)電機(jī)組的單機(jī)容量有了較快的增長[8]。

推力軸承的發(fā)展研究主要有物理和機(jī)械性能的研究、潤滑性能計算分析和試驗[9]。在推力軸承承載能力的研究方面，W.F.Hughes研究了油流慣性與靜壓推力軸承性能之間的關(guān)系，結(jié)果表明油流慣性對推力軸承的承載性能有較大的影響。印度學(xué)者V.K.Kapur和Kamlesh Vermal在對平行平板階梯軸承進(jìn)行研究時還考慮了溫度和慣性力對軸承壓力分布和承載能力的影響。N.Ganesan和T.Jayachandra Prabhu用線性三角單元有限元法對圓形腔圓錐靜壓推力軸承的承載能力進(jìn)行了深入研究，并考慮了油流慣性的影響，而且用試驗對結(jié)論進(jìn)行了驗證，試驗值與理論值吻合較好[10-16]。

對支撐較重軸向載荷的推力軸承，當(dāng)油膜厚度經(jīng)受高壓和高溫時，推力軸承容易發(fā)生過多的熱彈性變形，從而導(dǎo)致軸承故障。針對這種情況，Najar F A和Harmain G A通過在巴氏合金軸瓦瓦面下方埋設(shè)冷卻水管路來降低油膜溫度，然后采用有限差分法求解雷諾方程、能量方程和廣義傅里葉熱傳導(dǎo)方程對此結(jié)構(gòu)軸承進(jìn)行數(shù)值模擬，結(jié)果表明冷卻回路基本合理，與傳統(tǒng)冷卻系統(tǒng)相比，軸承整體溫度顯著降低[17]。Akbarzadeh P.開發(fā)了一種熱力學(xué)數(shù)值計算程序，用于對水電站中具有穩(wěn)定載荷的可傾瓦推力軸承的計算，對雷諾方程和油膜潤滑中的二維能量方程分別通過具有連續(xù)過松弛方案的中心差分技術(shù)和反向差分與迭代技術(shù)來求解。該計算程序考慮到了潤滑劑粘度隨溫度的變化和軸瓦的變形，計算結(jié)果可以揭示軸瓦數(shù)量、預(yù)載荷因子、轉(zhuǎn)速對軸承性能的影響[18]。LimingZhai、YongyaoLuo等人在考慮粘度-溫度效應(yīng)的前提下，把用于油膜的3D CFD模型和用于軸瓦和鏡板的3D FEA模型結(jié)合起來，對抽水蓄能機(jī)組中的雙向推力軸承進(jìn)行了數(shù)值模擬，預(yù)測了軸瓦和鏡板的變形，以及油膜厚度、壓力、溫度的分布，揭示了油膜中產(chǎn)生的高壓和高溫的原因以及鏡板的轉(zhuǎn)速對潤滑特性的影響，并在優(yōu)化設(shè)計中提出推力載荷應(yīng)該與油膜溫度和壓力平衡的建議[19]。E.Walicki和M.wasilczuk分別分析了軸承壁面多孔性對油流場壓力分布的影響以及推力軸承二維和三維流動，得到了速度、油膜厚度、溫度分布規(guī)律[20,21]。

伊泰普水電站水輪發(fā)電機(jī)推力軸承外徑520cm，總的推力負(fù)荷約3 600 t，在實際運(yùn)行過程中軸承損耗過大，Piotr Pajaczkowski和Michel Spiridon等人為了提高伊泰普水電站推力軸承的運(yùn)行效率，進(jìn)一步減少損耗，對推力軸承結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計，并在推力軸瓦表面添加PolypadTM涂層這種現(xiàn)代聚合物涂層材料，使得推力軸承可以在更高的溫度和更低的油膜厚度下運(yùn)行[22]。

在對推力軸承潤滑性能的試驗研究方面，可用電阻法、電容法、光纖傳感器法、電渦流法等來測量彈性流體動力潤滑油膜厚度[23]。Glavatskih S B.，Reddy D S K等人采用多個熱電偶記錄了水潤滑推力軸承瞬態(tài)運(yùn)行和穩(wěn)態(tài)運(yùn)行時的水膜溫度，并進(jìn)行了分析比較[24]。UusitaloOS和Glavatskih SB用電渦流傳感器對傾斜瓦推力軸承的油膜厚度與溫度進(jìn)行了監(jiān)測，并對傳感器設(shè)計、校準(zhǔn)程序的準(zhǔn)確性和敏感性進(jìn)行了描述[25]。為了使推力軸承潤滑性能的研究與實際運(yùn)行情況更加符合，研究人員在考慮溫度對潤滑性能影響的前提下，利用熱彈流分析理論，對不同油膜溫度分布模型進(jìn)行了分析。Elrod和Brewe等人通過勒記德多項式溫度輪廓近似法模擬油膜溫度場的分布來研究流體動壓推力軸承的潤滑性能[26]。Sharma和Pandey評估了勒記德多項式溫度輪廓近似法和拋物線多項式溫度輪廓近似法所得結(jié)果的精確性，認(rèn)為勒記德多項式溫度輪廓近似法的數(shù)值解更加符合實際工況[27]。

3 推力軸承國內(nèi)研究概況

我國的哈爾濱電機(jī)廠和東方電機(jī)廠的推力軸承設(shè)計制造技術(shù)代表了國內(nèi)水輪發(fā)電機(jī)推力軸承的水平。以哈爾濱電機(jī)廠為代表，我國從20世紀(jì)60年代初期開始對推力軸承技術(shù)進(jìn)行研究與試驗，經(jīng)過多年的發(fā)展，形成了科研開發(fā)、技術(shù)設(shè)計、工藝、試驗研究等有機(jī)結(jié)合的推力軸承科研開發(fā)體系[28]。我國水輪發(fā)電機(jī)組推力軸承巴氏合金瓦已經(jīng)應(yīng)用了多年，隨著機(jī)組單機(jī)容量和推力負(fù)荷的增大，為解決鎢金瓦推力軸承在電站運(yùn)行中的欠缺，1989年我國進(jìn)行了彈性金屬塑料瓦的研制，1990年第一套100 t級彈性金屬塑料瓦投入運(yùn)行，1997年，我國成功研制三峽6 000 t級彈性金屬塑料瓦推力軸承，并完成了全尺寸模擬試驗，取得了大量可靠的數(shù)據(jù)[29,30]。在設(shè)計技術(shù)上，哈爾濱電機(jī)廠自主開發(fā)了水輪發(fā)電機(jī)組推力軸承潤滑性能計算分析軟件TBearing/ANSYS，可以準(zhǔn)確地計算鏡板推力頭的變形。

在對推力瓦型面的研究方面，李忠、王風(fēng)才、袁曉陽等人提出了適用于推力軸承的螺旋面扇形瓦，并對其熱動力潤滑性能進(jìn)行分析，最終發(fā)現(xiàn)螺旋面扇形瓦具有承載能力大、油膜厚度大、溫升較小、粘性耗散較小、制造加工容易等諸多優(yōu)點，可以很好地取代平面扇形瓦[31]；陳志瀾、王海林等人以三峽水輪發(fā)電機(jī)推力軸承為研究對象，利用瓦面二次曲面數(shù)學(xué)模型和三維熱彈流潤滑性能分析軟件，對橫彎形面、反橫彎形面、馬鞍形面、平面形面、圓柱形面、球形形面6種瓦面形面推力軸承的最大油膜壓力、最高油膜溫度、最小油膜厚度等進(jìn)行了數(shù)值計算，通過比較分析計算，對三峽水輪發(fā)電機(jī)推力軸承進(jìn)行瓦面形面多目標(biāo)優(yōu)化設(shè)計，提高了軸承的潤滑性能[32]。

在推力軸承潤滑性能的研究方面，武中德[33]開發(fā)了一套二維TEHD有限元程序，通過有限元法聯(lián)立求解油膜厚度方程、能量方程、溫粘關(guān)系方程和雷諾方程，并用結(jié)構(gòu)有限元程序ADINA對推力瓦和鏡板推力頭進(jìn)行熱彈變形分析；張艷芹、陳瑤、范立國等人利用數(shù)值模擬的方法，在轉(zhuǎn)速、油腔深度、油腔面積相同時針對重型靜壓推力軸承橢圓形、扇形、工字形、矩形4種形狀的油腔研究了間隙油膜的壓力場，分析了4種油腔的壓力分布規(guī)律，對油腔結(jié)構(gòu)進(jìn)行了優(yōu)化[34]；盧德平對內(nèi)循環(huán)抽屜式冷卻器推力軸承、內(nèi)循環(huán)立式冷卻器推力軸承、鏡板泵外循環(huán)推力軸承、導(dǎo)瓦泵外循環(huán)推力軸承和外加泵外循環(huán)推力軸承5種不同結(jié)構(gòu)推力軸承油槽內(nèi)部的流場特征進(jìn)行對比分析，得到推力軸承承載力與間隙形狀、間隙尺寸、轉(zhuǎn)動部件轉(zhuǎn)速之間的關(guān)系以及相同間隙不同結(jié)構(gòu)推力軸承承載力的變化，對推力軸承的結(jié)構(gòu)進(jìn)行了優(yōu)化設(shè)計[35]。

20世紀(jì)70年代以來，東方電機(jī)有限公司在推力軸承冷卻方式方面進(jìn)行了大量的科研工作。一方面在內(nèi)循環(huán)結(jié)構(gòu)中，轉(zhuǎn)動部件上開孔促進(jìn)油路循環(huán)改善，采用散熱性能好的翅片管冷卻器代替光管冷卻器；另一方面對推力軸承外循環(huán)冷卻系統(tǒng)進(jìn)行研究，完成了烏江渡外加泵外循環(huán)和龍羊峽自身泵外循環(huán)系統(tǒng)的研制，并將這一成果應(yīng)用于寶珠寺175 MW和李家峽400 MW水輪發(fā)電機(jī)組推力軸承設(shè)計中[36]。國內(nèi)一些水電廠的運(yùn)行實踐也為推力軸承油循環(huán)冷卻系統(tǒng)性能的研究積累了許多寶貴的經(jīng)驗。例如，水布埡電廠水輪發(fā)電機(jī)推力軸承采用自身鏡板泵外循環(huán)冷卻方式，機(jī)組投運(yùn)以來，3、4號機(jī)組推力瓦溫運(yùn)行溫度70℃左右，與全國同類型推力瓦溫相比偏高，研究表明，瓦溫偏高與進(jìn)、出油邊油膜厚度比有關(guān)，對油循環(huán)冷卻系統(tǒng)進(jìn)行改造時在瓦面上預(yù)加工出型面，增大了進(jìn)、出油邊油膜厚度比值，提高了瓦的周向傾角，提高了軸承的承載能力，改造后推力瓦溫和平均瓦溫持平，達(dá)到了預(yù)期效果[37]。三峽左岸VGS機(jī)組推力軸承在設(shè)計時采用彈簧束支承，油槽容積較小，油冷卻器換熱能力偏小，造成推力軸承瓦間油溫比其他部位高4℃左右，針對左岸推力軸承存在的不足，東方電機(jī)廠在設(shè)計右岸推力軸承時，增大了推力軸承油槽的尺寸，提高了外循環(huán)冷卻器的換熱能力，通過加大瓦厚度、瓦間距、瓦塊徑向?qū)挾鹊却胧┨岣咄屏S承的承載能力，改進(jìn)了推力軸承油循環(huán)冷卻路徑，最終使得在相近的運(yùn)行油溫下，右岸推力軸承平均瓦溫、最高瓦溫與左岸推力軸承相比均降低7 K以上[38]。

4 推力軸承散熱系統(tǒng)的數(shù)值模擬

水輪發(fā)電機(jī)組推力軸承散熱系統(tǒng)（油循環(huán)冷卻系統(tǒng)）的油循環(huán)是一個比較復(fù)雜的問題，目前尚沒有較完整的試驗數(shù)據(jù)，難以準(zhǔn)確定量地計算，通常根據(jù)經(jīng)驗估算確定[39,40]。按照經(jīng)驗與規(guī)范，在常規(guī)設(shè)計中，軸承系統(tǒng)安全系數(shù)裕度較大，往往引起不必要的浪費(fèi)。目前，針對火電廠空氣冷卻系統(tǒng)方面的散熱研究成果比較多，而針對水電廠軸承油循環(huán)冷卻系統(tǒng)方面的數(shù)值模擬研究成果較少，隨著計算機(jī)技術(shù)的發(fā)展，目前可以采用現(xiàn)代流場計算分析Fluent軟件對推力軸承油循環(huán)冷卻系統(tǒng)進(jìn)行數(shù)值模擬，對以往的生產(chǎn)實踐經(jīng)驗進(jìn)行校正。

對推力軸承油循環(huán)冷卻系統(tǒng)中潤滑油的數(shù)值模擬，既需要通過獲得潤滑油的速度分布和壓力分布來判斷潤滑油的流態(tài)情況和軸瓦受力情況，又需要通過獲得潤滑油的溫度分布來判斷散熱系統(tǒng)的散熱性能，而潤滑油的速度場、壓力場、溫度場又互相影響，這就涉及到三維實體、流場和溫度場的耦合分析，其難點在于計算很難收斂，而且需要占用大量的計算機(jī)資源，對計算機(jī)的配置要求很高，大量的時間都花費(fèi)在模型的簡化以及網(wǎng)格的劃分上，工作繁重，技術(shù)難度比較大。由于推力軸承油槽內(nèi)的油流場具有周期性變化的特點，在結(jié)合電站具體參數(shù)建立推力軸承的模型時，可以建立模型的1/Z（Z指推力軸瓦個數(shù)）進(jìn)行分析，從而減少網(wǎng)格數(shù)目，降低計算量。

采用Gambit或者ICEM建立推力軸承的物理模型，對建立好的模型進(jìn)行網(wǎng)格劃分之后，在Fluent中選擇合適的算法、設(shè)定初始條件和邊界條件、添加湍流模型對潤滑油的連續(xù)方程、動量方程、能量方程

進(jìn)行求解，獲得潤滑油的速度、壓力、溫度分布。（1）連續(xù)性方程

式中：ρ是油的密度；

t是時間；

ux、uy、uz分別是速度矢量在直角坐標(biāo)系

x、y、z方向上的分量。

（2）動量方程

動量守恒定律是任何流動系統(tǒng)必須滿足的基本定律，該定律實際上是牛頓第二定律，在流體力學(xué)中進(jìn)一步表示為Navier-Stokes方程：

式中：ρ是油體微元上的壓力；

Su、Sv、Sw是 動 量 方 程 的 廣 義 源 項 ，

Fx、Fy、Fz是微元體上的體力，若體力只有重力，且 z軸豎直向上，則有 Fx=0，F(xiàn)y=0，F(xiàn)z=-ρg。一般而言，sx、sy、sz是小量，其表達(dá)式如式（5）至（7）所示：

對于粘性為常數(shù)的不可壓縮流體，sx=sy=sz=0。

（3）能量方程

寫成展開形式為：

式中：cp是油的比熱；

T是油溫；

k為油的傳熱系數(shù)；

ST為油的內(nèi)熱源及由于粘性作用油的機(jī)

械能轉(zhuǎn)換為熱能的部分，簡稱為粘性耗散項。

推力軸承油循環(huán)冷卻系統(tǒng)中油的傳熱方式屬于對流傳熱。對流傳熱是指流體流過一個物體表面時流體與物體表面之間的熱量傳遞過程。冷卻器管內(nèi)通入冷卻水，在分析過程中將水看作低速不可壓縮流體，在假定油流沒有相變的基礎(chǔ)上，對控制方程進(jìn)行求解。

通過上述方法，可以對不同支承結(jié)構(gòu)的推力軸承油循環(huán)冷卻系統(tǒng)進(jìn)行數(shù)值模擬，并比較不同支承結(jié)構(gòu)推力軸承油循環(huán)冷卻系統(tǒng)的散熱性能；可以對不同推力瓦塊數(shù)的同一類型推力軸承油循環(huán)冷卻系統(tǒng)進(jìn)行數(shù)值模擬，從而確定最合適的推力瓦的個數(shù)；可以改變各冷卻器的布置間距對推力軸承油循環(huán)冷卻系統(tǒng)進(jìn)行數(shù)值模擬，來確定合適的冷卻器布置形式；可以對不同工況下某一確定的推力軸承油循環(huán)冷卻系統(tǒng)進(jìn)行數(shù)值模擬，來研究該推力軸承油循環(huán)冷卻系統(tǒng)的流態(tài)特征和散熱性能，為推力軸承油循環(huán)冷卻系統(tǒng)結(jié)構(gòu)改進(jìn)和安全穩(wěn)定運(yùn)行提供合理建議。

5 結(jié)語

自20世紀(jì)初以來，在推力軸承的承載能力、潤滑性能的試驗、推力瓦型面的研究等方面取得了大量的研究成果，推動了推力軸承技術(shù)的發(fā)展，也促進(jìn)了水電機(jī)組單機(jī)容量的增長，但在推力軸承油循環(huán)冷卻系統(tǒng)內(nèi)部流場性能，特別是散熱性能的研究方面成果還比較少，隨著計算機(jī)技術(shù)的發(fā)展，這將成為一個重要的研究方向，不僅有助于我們進(jìn)一步了解推力軸承油循環(huán)冷卻系統(tǒng)的散熱性能，而且對冷卻器的結(jié)構(gòu)改進(jìn)、機(jī)組的穩(wěn)定運(yùn)行以及同類型電站的安裝檢修都具有重要參考價值。

[1]胡贛生，晏蓮清.對葛洲壩水利樞紐二江電廠1號機(jī)組推力軸瓦燒損的初步探討 [J].葛洲壩水電工程學(xué)院學(xué)報，1982（1）：1-11.

[2]張奇.600 MW大型水輪發(fā)電機(jī)組推力軸承瓦溫和油溫高原因分析及處理[J].華電技術(shù)，2015，37（10）：54-57.

[3]萬玉倩，田子勤.三峽左岸電站VGS機(jī)組推力瓦溫偏高及改善措施[J].人民長江，2009，40（2）：26-28.

[4]陳占發(fā).白山300 MW機(jī)組推力軸承燒瓦原因分析及改進(jìn)方法[J].水利水電技術(shù)，1988（8）：25-31.

[5]國家電網(wǎng)安監(jiān)〔2011〕2024號.關(guān)于印發(fā)《國家電網(wǎng)公司安全事故調(diào)查規(guī)程》的通知[Z]，2011：12-14.

[6]陳錫芳.水輪發(fā)電機(jī)結(jié)構(gòu)運(yùn)行監(jiān)測與維修[M].北京：中國水利水電出版社，2008：241～250.

[7]劉琪.水輪發(fā)電機(jī)組推力軸承技術(shù)的發(fā)展[J].黑龍江科學(xué)，2014，5(3)：67-68.

[8]Pradhan S，Hacioglu K，Krugler V，et al.Support Vector Learning for Semantic Argument Classification[J].Machine Learning，2005，60(1)：11-39.

[9]武中德，張宏.水輪發(fā)電機(jī)組推力軸承技術(shù)的發(fā)展[J].電器工業(yè)，2007(1)：32-36.

[10]張文光，王家序.玻纖增強(qiáng)改性酚醛樹脂復(fù)合材料在水潤滑軸承上的應(yīng)用研究[J].潤滑與密封，2007（9）：116-118.

[11]Sukumaran N V P，Prabhakaran N K.Finite element analysis of elastohydynamic circular journal bearing with micropolar lubricants[J].Finite element in Analysis and Design，2004，41(1)：59-75.

[12]Chambers W S，Mikua A M.Operational data for a large vertical thrust bearing of a pumped storage application[J].Transactionsof the Society of Tribologistsand eers，1988，31(1)：61-65.

[13]Raimondi A A.The influence of longitudinal and transverse profile on the load capacity of pivoted pad bearings[J].ASLE Trans，1960，3(2)：173-180.

[14]Tieu A K.A numerical simulation of finite-width thrust bearingstaking into account viscosity variation with temperature and pressure[J].JMech Eng Sic，1974，17(10)：1-10.

[15]Kurban A O.Analysisof a hydrodynamic thrust bearing with elastic deformation using a recurrent neural network[J].Tribology International，2003，36(12)：943-948.

[16]Reddy D S K，Swarnamani S，Prabhu B S.Thermo elastohydrodynamic analysis of tilting pad journal bearing-Theory and experiments[J].ASME Tribology Transactions，2000，43(1)：82-90.

[17]Najar F A，Harmain G A.Novel Approach Towards Thrust Bearing Pad Cooling[C]//ASME 2014 Gas Turbine India Conference.American Society of Mechanical Engineers，2014.

[18]Akbarzadeh P.Numerical Study of Thermohydrodynamic Characteristics of Oil Tilting-Pad Journal Bearings with a Self-Pumping Fluid Flow Circulation[J].Tribology Transactions，2015，58(1)：18-30.

[19]Zhai L，Luo Y，Wang Z，et al.3D Two-way coupled TEHD analysis on the lubricating characteristics of thrust bearings in pump-turbine units by combining CFD and FEA [J].Chinese Journal of Mechanical Engineering，2016，29(1)：112-123.

[20]Walicki E，Walicki A，Makhaniok A.Pressure distribution in a curvilinear thrust bearing with one porous wall lubricated by aBingham fluid[J].Meccanica，2001，36：709-716.

[21]Wasilczuk M，Rotta G.Modeling lubricant flow between thrust-bearing pads[J].Tribology International，2008，41：908-913.

[22]Pajaczkowski P，Spiridon M，Schubert A，et al.Itaipu Binacional hydro power plant thrust bearing design optimization for higher efficiency[C]//Proceedings of 2nd International Brazilian Conference on Tribology，F(xiàn)ozda Igacu，Parana，Brazil，2014：3-5.

[23]何春勇.潛水泵水潤滑推力軸承潤滑性能數(shù)值計算研究[D].武漢：武漢理工大學(xué)，2010.

[24]Glavatskih S B.A method of temperature monitoring in fluid filmbearings[J].TribologyInternational，2004，37(2)：143-148.

[25]Bathe K J，Zhang H，Ji S H.Finite element analysis fluid flows fully coupled with structural interactions[J].Computer sand Structures，1999，72(1-3)：1-16.

[26]Elord H G，Brewe D E.Thermohydrodynamic analysis for laminar lubricating films[R].NASA Technical Memorandu，1986.

[27]Sharma R K，Pandey R K.An investigation into the validity of the temperature profile approximations across the film thicknessin THD analysisof infinitely wide slider bearing[J].Tribology Online，2006，1(1)：19-24.

[28]吳軍令.3 000 t推力軸承試驗臺研制成功 [J].電氣應(yīng)用，1993（3）：39-40.

[29]武中德，張宏，梁廣泰，等.三峽水輪發(fā)電機(jī)推力軸承[J].中國三峽建設(shè)，2003（9）：7-8.

[30]武中德，張宏，王黎欽，等.大型水輪發(fā)電機(jī)彈性金屬塑料瓦推力軸承技術(shù)[J].潤滑與密封，2002（2）：59-62.

[31]李忠,王風(fēng)才,袁曉陽,等.螺旋面扇形瓦推力軸承熱動力潤滑性能分析[J].機(jī)械科學(xué)與技術(shù)，1999，18(3)：409-411.

[32]陳志瀾，袁小陽，王海林，等.推力軸承瓦面形面對潤滑性能影響的研究[J].摩擦學(xué)學(xué)報，2003，23（1）：56-59.

[33]武中德.大型水輪發(fā)電機(jī)推力軸承熱彈流潤滑性能研究[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2001.

[34]張艷芹，陳瑤，范立國，等.四種油腔形狀重型靜壓軸承承載性能理論分析[J].哈爾濱理工大學(xué)學(xué)報，2013，18（2）：68-71.

[35]盧德平.大型水輪發(fā)電機(jī)推力軸承油槽內(nèi)流場分析和結(jié)構(gòu)優(yōu)化[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2010.

[36]劉漢城，李康揚(yáng).三峽機(jī)組推力軸承油路循環(huán)及冷卻方式研究[J].東方電氣評論，1997，11（2）：123-136.

[37]唐衛(wèi)東，曹凱，趙公杰，等.水布埡電廠3、4號機(jī)組推力軸承運(yùn)行情況分析[J].水電與新能源，2014（2）：59-62.

[38]鐘海權(quán)，鐸林，張?zhí)禊i.三峽右岸水輪發(fā)電機(jī)組推力軸承設(shè)計與運(yùn)行[J].東方電機(jī)，2008（5）：22-29.

[39]張宏，劉琪，武中德，等.大型水輪發(fā)電機(jī)推力軸承油槽內(nèi)油溫度場[J].大電機(jī)技術(shù)，2014(7)：109-112.

[40]武中德，張宏.大型水輪發(fā)電機(jī)推力軸承設(shè)計[C]//水電站設(shè)備的研究與實踐—第十七次中國水電設(shè)備學(xué)術(shù)討論會論文集，2009：230-235.