孫興偉，張繼偉，馬鐵強(qiáng)

(沈陽(yáng)工業(yè)大學(xué) 機(jī)械工程學(xué)院，遼寧 沈陽(yáng) 110870)

風(fēng)力發(fā)電機(jī)軸承接觸特性研究

孫興偉，張繼偉，馬鐵強(qiáng)

(沈陽(yáng)工業(yè)大學(xué) 機(jī)械工程學(xué)院，遼寧 沈陽(yáng) 110870)

針對(duì)風(fēng)力發(fā)電機(jī)內(nèi)部結(jié)構(gòu)和所選用的軸承的特殊性，選擇10 kW風(fēng)力發(fā)電機(jī)結(jié)構(gòu)，根據(jù)6311軸承的具體參數(shù)建立三維模型，合理設(shè)置模型邊界條件，利用有限元方法對(duì)接觸應(yīng)力與變形進(jìn)行仿真計(jì)算，得到了整體深溝球軸承與單個(gè)滾動(dòng)體與內(nèi)外滾道接觸的等效應(yīng)力、等效應(yīng)變、徑向變形和接觸應(yīng)力的變化情況，驗(yàn)證了利用有限元解決發(fā)電機(jī)軸承接觸問(wèn)題是可行的，同時(shí)也符合軸承受載荷作用下的接觸應(yīng)力與變形實(shí)際情況。對(duì)風(fēng)力發(fā)電機(jī)軸承設(shè)計(jì)、參數(shù)優(yōu)化以及裝配調(diào)試具有一定指導(dǎo)意義。

風(fēng)力發(fā)電機(jī)；深溝球軸承； 有限元分析； 接觸應(yīng)力

0 引言

風(fēng)力發(fā)電機(jī)是將風(fēng)能轉(zhuǎn)換為機(jī)械功，機(jī)械功帶動(dòng)轉(zhuǎn)子旋轉(zhuǎn)，最終輸出交流電的電力設(shè)備，即風(fēng)力發(fā)電機(jī)是一種以風(fēng)力為能源，以大氣為工作介質(zhì)的發(fā)電機(jī)。與普通的發(fā)電機(jī)不同在于其需要在波動(dòng)的機(jī)械能條件下運(yùn)轉(zhuǎn)，風(fēng)力發(fā)電機(jī)主要包括雙饋型風(fēng)力發(fā)電機(jī)、直驅(qū)式永磁風(fēng)力發(fā)電機(jī)、半直驅(qū)式永磁風(fēng)力發(fā)電機(jī)[1]。作為支撐風(fēng)力發(fā)電機(jī)的重要部件軸承，需要充分考慮其使用環(huán)境和高可靠性和高經(jīng)濟(jì)性并存的最佳設(shè)計(jì)。由于風(fēng)力發(fā)電機(jī)運(yùn)行環(huán)境惡烈，維護(hù)困難，對(duì)軸承的使用壽命和可靠性要求高；另外，由于發(fā)電機(jī)軸承連續(xù)高速工作，對(duì)軸承的精度提出了更高的要求。風(fēng)力發(fā)電機(jī)組中的軸承包括偏航系統(tǒng)軸承、變槳距系統(tǒng)軸承和傳動(dòng)系統(tǒng)軸承，傳動(dòng)系統(tǒng)軸承包括了主軸軸承、變速箱軸承和發(fā)電機(jī)軸承，發(fā)電機(jī)軸承一般采用3套深溝球軸承和圓柱滾子軸承，要求高精度、低振動(dòng)和低噪聲，發(fā)電機(jī)軸承需要進(jìn)行結(jié)構(gòu)和優(yōu)化設(shè)計(jì)[2-3]，因此對(duì)發(fā)電機(jī)軸承進(jìn)行接觸應(yīng)力分析是很有必要的，用有限元軟件模擬發(fā)電機(jī)軸承的復(fù)雜的實(shí)際工況可以更好的分析計(jì)算軸承接觸應(yīng)力問(wèn)題。由于發(fā)電機(jī)軸承的運(yùn)動(dòng)參數(shù)和受力是個(gè)復(fù)雜的問(wèn)題，這其中會(huì)涉及到很多變量，但各變量之間是相互聯(lián)系的[4]。國(guó)內(nèi)外學(xué)者對(duì)發(fā)電機(jī)軸承的研究甚少，因此工程上常采用擬靜力學(xué)方法分析發(fā)電機(jī)軸承接觸問(wèn)題。

以10 kW的風(fēng)力發(fā)電機(jī)中用的6311深溝球軸承為研究對(duì)象，用有限元方法對(duì)其擬靜力學(xué)分析，建立6311整體和局部單個(gè)滾動(dòng)體與內(nèi)外滾道接觸的模型，研究軸承的接觸應(yīng)力和接觸變形規(guī)律。

1 發(fā)電機(jī)軸承內(nèi)部載荷分布

風(fēng)力發(fā)電機(jī)軸承的轉(zhuǎn)速通常不是很高的，滾動(dòng)體慣性力、陀螺力矩以及接觸角變化對(duì)載荷分布的影響可以忽略不計(jì)。研究發(fā)電機(jī)軸承內(nèi)部載荷分布是研究風(fēng)力發(fā)電機(jī)壽命和可靠性的前提條件，因此從理論上研究發(fā)電機(jī)軸承載荷變形關(guān)系[5]。如圖1所示的是深溝球軸承所受徑向載荷作用下的內(nèi)外圈在外力作用下發(fā)生的相對(duì)變形量是δr。根據(jù)已知第q個(gè)滾動(dòng)體與內(nèi)外圈之間的總的彈性變形為

δq=δrcosψq

(1)

式中，q為滾動(dòng)體的序號(hào)，位于徑向載荷作用下記為0，兩邊對(duì)稱分布，分別為1,2,3…；ψq為第q個(gè)滾動(dòng)體的位置角，徑向載荷作用下ψq=0。

從圖1中可以得知：當(dāng)ψq=0°時(shí)，δ0=δr=δmax，即0號(hào)滾動(dòng)體與內(nèi)外滾道接觸變形最大，其受徑向載荷最大；當(dāng)ψq=90°，δq=0，表明承載區(qū)的極限角度ψ1=±90°，受載荷為半圈滾動(dòng)體。

第q個(gè)滾動(dòng)體接觸應(yīng)力為

Qq=Qmax=cosnψq

(2)

內(nèi)圈的平衡方程為

(3)

(4)

式中，k為受載荷最小的滾動(dòng)體的序號(hào)；位置ψq角與軸承中滾動(dòng)體的數(shù)目Z和滾動(dòng)體的序號(hào)有關(guān)；指數(shù)n與接觸類(lèi)型有關(guān)。

考慮到軸承中的徑向游隙，因此對(duì)于深溝球軸承而言，其最大的接觸載荷公式為[6]

(5)

(6)

式中，ea、eb是赫茲接觸系數(shù)；Σρ為主曲率和；Z是滾動(dòng)體數(shù)目。

圖1 深溝球軸承內(nèi)部載荷分布

2 軸承有限元計(jì)算

2.1 模型建立和網(wǎng)格劃分

利用有限元軟件建立6311深溝球軸承整體模型和單一滾動(dòng)體與內(nèi)外圈滾道接觸局部模型，由于軸承倒角等細(xì)微結(jié)構(gòu)對(duì)接觸應(yīng)力分析結(jié)果有很微小的影響，再加上滾動(dòng)體與內(nèi)外圈滾道之間所產(chǎn)生的油膜厚度相對(duì)于軸承的整體尺寸很小，故綜上所述用有限元方法在分析接觸應(yīng)力與變形時(shí)可以對(duì)6311深溝球軸承模型簡(jiǎn)化，不考慮軸承油膜厚度對(duì)軸承接觸應(yīng)力的影響，忽略軸承保持架，所以有限元模型就僅考慮內(nèi)外圈、滾動(dòng)體等軸承的基本結(jié)構(gòu)。軸承的具體參數(shù)如表2所示，軸承材料選擇及其主要指標(biāo)如表1所示。

表1 軸承材料的各項(xiàng)指標(biāo)

表2 6311深溝球軸承主要參數(shù)

整體軸承有限元模型如圖2a所示，對(duì)單個(gè)滾動(dòng)體與內(nèi)外滾道接觸的有限元模型如圖2b所示。

圖2 軸承有限元模型

將軸承有限元模型進(jìn)行網(wǎng)格劃分，采用有限元下自動(dòng)劃分網(wǎng)格形式對(duì)軸承進(jìn)行劃分，設(shè)置內(nèi)外滾道的單元網(wǎng)格尺寸為2 mm，滾動(dòng)體的單元網(wǎng)格尺寸為1 mm，接觸區(qū)域進(jìn)行局部加密細(xì)化，深溝球軸承整體網(wǎng)格劃分共生成558 294個(gè)單元和838 308個(gè)節(jié)點(diǎn)；對(duì)單個(gè)滾動(dòng)體和內(nèi)外滾道間的接觸的網(wǎng)格劃分共生成23 735個(gè)單元和45 571個(gè)節(jié)點(diǎn)。

2.2 定義接觸設(shè)置

在有限元中提供了五種接觸的類(lèi)型，分析接觸應(yīng)力和接觸變形，本文采用摩擦接觸類(lèi)型；在分析軸承時(shí)建立了兩個(gè)接觸對(duì)：一是內(nèi)圈滾道與滾動(dòng)體的接觸對(duì)，二是外圈滾道與滾動(dòng)體的接觸對(duì)[7]。這兩個(gè)接觸對(duì)都是以滾動(dòng)體表面為接觸面。由于滾道之間的接觸區(qū)域是橢圓形的，長(zhǎng)度可以近視的認(rèn)為是長(zhǎng)方形，根據(jù)多面體部件面與面接觸時(shí)設(shè)置的原則及結(jié)合對(duì)6311球軸承的分析設(shè)置，本研究中設(shè)置接觸對(duì)的摩擦因數(shù)MU=0.003，法向接觸剛度FKN=0.1，初始接觸因子ICONT=0.1，接觸算法設(shè)置為增廣拉格朗日法。整個(gè)軸承滾動(dòng)體與單個(gè)滾動(dòng)體與內(nèi)外滾道間的接觸設(shè)置如圖3所示。

圖3 軸承網(wǎng)格設(shè)置

2.3 邊界條件設(shè)定

根據(jù)軸承工作的實(shí)際條件設(shè)定軸承接觸分析的邊界條件。

(1)根據(jù)工作條件設(shè)定在有限元分析中對(duì)外圈施加全約束來(lái)模擬軸承外圈固定和軸承座對(duì)軸承外圈在x、y、z方向的自由度；

(2)模擬軸承實(shí)際的工況條件對(duì)軸承的內(nèi)外圈的z方向平動(dòng)自由度約束；

(3)模擬軸承的軸向運(yùn)動(dòng)約束軸承的內(nèi)外圈側(cè)面；

(4)模擬軸承保持架的作用，限制每個(gè)滾動(dòng)體與內(nèi)外滾道間接觸的軸向與周向的自由度，允許其滾動(dòng)體在徑向載荷作用下發(fā)生彈性變形[8-9]。

由于深溝球軸承只承受徑向載荷，故對(duì)軸承的內(nèi)圈內(nèi)表面下半部分施加徑向載荷，在軸承內(nèi)圈施加轉(zhuǎn)速，模擬軸承與發(fā)電機(jī)軸一起旋轉(zhuǎn)。

3 有限元仿真結(jié)果分析

通過(guò)有限元分析對(duì)軸承的滾動(dòng)體與內(nèi)外滾道間的接觸應(yīng)力、接觸變形分析，通過(guò)分析可以得到軸承整體與單個(gè)滾動(dòng)體的等效應(yīng)力云圖、等效應(yīng)變圖。圖4為在Fr=2 kN時(shí)整體軸承的等效應(yīng)力、等效應(yīng)變與綜合位移圖，圖5為相同條件下的單個(gè)滾動(dòng)體與內(nèi)外滾道間的等效應(yīng)力、等效應(yīng)變與徑向變形云圖，圖6為整體軸承接觸應(yīng)力云圖。

圖4 整體軸承的等效應(yīng)力、等效應(yīng)變、綜合位移云圖

圖5 單個(gè)滾動(dòng)體的等效應(yīng)力、等效應(yīng)變、綜合位移云圖

圖6 滾動(dòng)軸承接觸應(yīng)力

從對(duì)整體深溝球軸承的有限元仿真結(jié)果可以得出，整體球軸承的最大等效應(yīng)力67.738 MPa，出現(xiàn)在徑向力所作用下的軸承內(nèi)圈所在的滾動(dòng)體接觸區(qū)域，即軸承的最大等效應(yīng)力的部位正好是軸承疲勞損壞的部位即軸承滾動(dòng)體與內(nèi)外滾道之間的接觸點(diǎn)處；從圖4b軸承的等效應(yīng)變圖可知，最大等效應(yīng)變是0.363 64 μm，承受徑向載荷的滾動(dòng)體處才可以產(chǎn)生等效應(yīng)變，從有限元仿真結(jié)果可以知道最大等效應(yīng)力與等效應(yīng)變發(fā)生的位置是相同的。由圖4c可知軸承的徑向載荷變形為4.4175 μm；由圖6軸承接觸應(yīng)力圖可知，最大接觸應(yīng)力是48.799 MPa，其最大接觸應(yīng)力的位置恰好是赫茲理論分析所得到的位置，也就是內(nèi)圈所受載荷最大的滾動(dòng)體上，對(duì)單個(gè)滾動(dòng)體與內(nèi)外滾道間接觸有限元仿真結(jié)果，最大等效應(yīng)力是21.383 MPa，最大等效應(yīng)變是1.0495 μm，最大徑向變形是1.3235 μm，其最大等效應(yīng)變和等效應(yīng)變力出現(xiàn)的位置和整體球軸承位置相同。從單個(gè)滾動(dòng)體與內(nèi)外滾道接觸的分析，可以更好的分析球軸承接觸應(yīng)力和接觸變形。

軸承的接觸應(yīng)力最大的部位正好是軸承所受載荷的滾動(dòng)體與內(nèi)外滾道間的接觸區(qū)域，滾動(dòng)體局部的接觸應(yīng)力很大，上半部分滾動(dòng)體與滾道的接觸應(yīng)力為零，每個(gè)滾動(dòng)體與內(nèi)、外圈滾道的接觸區(qū)域呈橢圓狀分布，接觸應(yīng)力值從中心向外逐漸變小。

4 不同工況下發(fā)電機(jī)軸承接觸應(yīng)力變化

利用有限元方法分析在不同工況條件下的發(fā)電機(jī)軸承的等效應(yīng)力、等效應(yīng)變、綜合位移、接觸應(yīng)力時(shí)，應(yīng)在不考慮其他條件下分別對(duì)軸承施加不同的載荷，施加的徑向載荷由2 kN到9 kN，用有限元軟件分析其等效應(yīng)力與應(yīng)變，徑向載荷變形和接觸應(yīng)力，所得到的結(jié)果如表3所示。隨著徑向負(fù)荷的增大，軸承滾動(dòng)體與內(nèi)外圈的接觸應(yīng)力由48.799 MPa增加到134.85 MPa,表明其接觸應(yīng)力在逐漸增大，但從表3中可以得知增大趨勢(shì)逐漸放慢。

表3 整體滾動(dòng)軸承不同載荷下的有限元仿真結(jié)果

5 結(jié)論

(1)采用有限元軟件分析計(jì)算風(fēng)力發(fā)電機(jī)6311軸承接觸狀態(tài)，可以較為直觀的反映出風(fēng)力發(fā)電機(jī)軸承工作運(yùn)行狀態(tài)，建立正確的軸承有限元模型可以詳細(xì)地描述軸承滾動(dòng)體與內(nèi)外滾道接觸力的傳遞和變化，可見(jiàn)利用有限元方法對(duì)解決風(fēng)力發(fā)電機(jī)軸承內(nèi)部結(jié)構(gòu)的接觸分析是行之有效的；

(2)利用有限元軟件較為方便的得出軸承的滾動(dòng)體與內(nèi)外滾道間的接觸應(yīng)力云圖、等效應(yīng)力與應(yīng)變?cè)茍D。從云圖中可以了解最大接觸應(yīng)力、最大等效應(yīng)力與徑向變形的位置具有一致性，其位置是最低端滾動(dòng)體與內(nèi)外滾道接觸的位置，從底端開(kāi)始兩側(cè)的接觸應(yīng)力逐漸減小，這與實(shí)際情況相同。

(3)利用有限元分析單一滾動(dòng)體與內(nèi)外滾道間的接觸應(yīng)力比整體分析更為直觀的表達(dá)了滾動(dòng)體的接觸應(yīng)力變化情況，可以方便地得出不同部位的接觸應(yīng)力變化規(guī)律，為風(fēng)力發(fā)電機(jī)軸承的選用設(shè)計(jì)和優(yōu)化提供良好的依據(jù)，對(duì)發(fā)電機(jī)軸承的疲勞壽命分析奠定基礎(chǔ)，有助于將在此基礎(chǔ)上研究在極限風(fēng)速狀態(tài)下深溝球軸承接觸應(yīng)力變化。

[1] 朱輝.風(fēng)力發(fā)電機(jī)組滾動(dòng)軸承故障的無(wú)傳感器診斷方法研究[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2015.

[2] 陳紅濤.兆瓦級(jí)風(fēng)力發(fā)電機(jī)組變槳軸承有限元分析[D].洛陽(yáng)：河南科技大學(xué)，2012.

[3] 王揚(yáng).風(fēng)力發(fā)電機(jī)組轉(zhuǎn)盤(pán)軸承疲勞壽命的研究[D].北京：華北電力大學(xué)，2012.

[4] 蔣立冬，應(yīng)麗霞.高速重載滾動(dòng)軸承接觸應(yīng)力和變形的有限元分析[J].機(jī)械設(shè)計(jì)與制造，2008(10): 62-64.

[5] 鄧四二，賈永川.推力球軸承摩擦力矩特性研究[J].兵工學(xué)報(bào)，2015，36(7)：1615-1623.

[6] 鄧四二，賈群義.滾動(dòng)軸承設(shè)計(jì)原理[M].北京:中國(guó)標(biāo)準(zhǔn)出版社，2008.

[7] 童寶宏，劉穎，蘇葒，等.滾動(dòng)軸承動(dòng)態(tài)接觸特性數(shù)值模擬中若干問(wèn)題[J].機(jī)械工程學(xué)報(bào)，2012，48(21)：116-123．

[8] 張福星，鄭源，汪清，等.基于ANSYS Workbench 的深溝球軸承接觸應(yīng)力有限元分析[J].機(jī)械設(shè)計(jì)與制造，2012(10):222-224.

[9] 朱福先，周金宇，孫奎洲，等.基于有限元的深溝球軸承應(yīng)力和剛度研究[J].機(jī)床與液壓，2014,42(17)：42-44.

Research on characteristics of wind turbine ball bearing

SUN Xing-wei，ZHANG Ji-wei，MA Tie-qiang

(School of Mechanical Engineering，Shenyang University of Technology, Shenyang 110870, China)

According to the internal structure of the wind turbine and the special characteristics of its bearing, the structure of 10 kW wind power generator is selected in this paper. The bearing 3D model was built according to the parameters of 6311 deep groove ball bearing with the reasonable setting of the boundary conditions, the contact stress and the deformation were simulation calculated by using the finite element method, the change of equivalent stress and strain, radial deformation and contact stress were obtained between overall deep groove ball bearings and single rolling body and the inner or outer raceways. Furthermore, the results show that solving the bearing contact problem with FEM is feasible and it coincides with the stress and deformation of the actual situation, and it can offer guidance for the design and debug of wind turbine bearing.

wind turbines; deep groove ball bearing；FEM；contact stress

2016-09-22；

2016-11-09

國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(51537007)，遼寧省科技創(chuàng)新重大專(zhuān)項(xiàng)(201303005)

孫興偉(1970-)，遼寧人，女，教授，博士，主要研究方向：復(fù)雜曲面制造技術(shù)及專(zhuān)用數(shù)控機(jī)床開(kāi)發(fā)。

張繼偉(1989-)，山西人，男，碩士生，主要研究方向：復(fù)雜曲面制造與大型能源裝備研發(fā)。

TH133.33，TK83

A

1001-196X(2017)03-0034-05