楊靜，廖暉

(東方電氣風(fēng)電有限公司，四川德陽，618000)

風(fēng)力發(fā)電機(jī)組有限元分析時(shí)偏航變槳軸承的非線性模擬

楊靜，廖暉

(東方電氣風(fēng)電有限公司，四川德陽，618000)

偏航變槳軸承是風(fēng)力發(fā)電機(jī)組構(gòu)成和傳遞載荷的重要部件，風(fēng)力發(fā)電機(jī)組主要部件有限元強(qiáng)度分析時(shí)，必須正確考慮偏航和變槳軸承滾子和滾道對載荷傳遞和應(yīng)力分布的影響，文章使用ANSYS軟件中的link10單元組模擬滾子，不僅可以正確考慮軸承的非線性特征，而且可以大大簡化計(jì)算的模型，使風(fēng)力機(jī)主要部件的有限元強(qiáng)度分析變得快捷和精確。

風(fēng)力發(fā)電機(jī)組，有限元分析，偏航軸承，變槳軸承，非線性

0 引言

偏航變槳軸承是風(fēng)力發(fā)電機(jī)組構(gòu)成和傳遞載荷的重要部件，風(fēng)力發(fā)電機(jī)組主要部件（比如輪轂、機(jī)架、塔頂法蘭等）有限元分析時(shí)，都需要考慮其臨近部件，如偏航或變槳軸承。如何模擬偏航變槳軸承是風(fēng)力發(fā)電機(jī)組部件有限元模型正確模擬的關(guān)鍵，偏航變槳軸承的非線性特性很大程度上影響著有限元分析的精確性，必須正確考慮，然而偏航變槳軸承的內(nèi)外圈和滾子全部實(shí)體建模將會帶來復(fù)雜的接觸關(guān)系，計(jì)算難以收斂。本文使用ANSYS軟件中的link10單元組模擬滾子，通過合理設(shè)置link10單元的剛度和初始應(yīng)變，使link10單元組和軸承滾道滾子接觸對具有相同的非線性接觸剛度，而且在載荷改變時(shí)也可以模擬出滾道滾子接觸角變化的特性。最后通過有限元計(jì)算和軸承理論計(jì)算比較，驗(yàn)證了其載荷傳遞的正確性。本文的方法不僅可以正確考慮風(fēng)力發(fā)電機(jī)組有限元模型中軸承的非線性特性，而且可以大大簡化模型，減少計(jì)算量，使風(fēng)力發(fā)電機(jī)組主要部件的有限元分析變得快捷和精確。

1 偏航變槳軸承結(jié)構(gòu)受載特點(diǎn)

偏航變槳軸承通常為四點(diǎn)接觸盤轉(zhuǎn)軸承或八點(diǎn)接觸盤轉(zhuǎn)軸承，圖1為典型的四點(diǎn)接觸盤轉(zhuǎn)軸承的截面示意圖。偏航變槳軸承通常采用負(fù)游隙或0游隙，軸承溝曲率半徑系數(shù)為0.52～0.53，滾道的半徑比滾子半徑略大，在軸承不受載荷的初始狀態(tài)時(shí)，滾道和滾子有4個(gè)接觸點(diǎn)，如圖1所示，這4個(gè)接觸點(diǎn)構(gòu)成2個(gè)接觸對，接觸對1處于45°接觸角，接觸對2處于135°接觸角。當(dāng)偏航變槳軸承受載后，隨著載荷的加大，滾子和滾道基本上只有一個(gè)接觸對接觸，在另一個(gè)接觸對出現(xiàn)間隙，并且接觸角也在變化，滾道和滾子的接觸剛度隨載荷的增大也會變大。軸承在受載后的變化過程是一個(gè)典型的非線性變化過程。

圖1 偏航變槳軸承截面示意圖

2 偏航變槳軸承基本方程

接觸對滾道之間的法向趨近量與接觸載荷的關(guān)系為非線性[1]，具體為：

接觸對1處接觸載荷Q1ψ為：

式中：δ1n和δ2n分別為接觸對1和接觸對2被滾動體隔開的2個(gè)滾道之間的法向趨近量，等于鋼球與每個(gè)滾道的趨近量之和；Kn為2個(gè)滾道間的載荷位移系數(shù)。

如果軸承受軸向力Fa、徑向力Fr以及翻轉(zhuǎn)力矩M作用，根據(jù)平衡條件列出的方程組為：

同理，接觸對2處接觸載荷Q2ψ為：

式中：鋼球總數(shù)為z；第i個(gè)鋼球的方位角為ψ；接觸對1的接觸角為α1ψ；接觸對2的接觸角為α2ψ。

式（3）是一個(gè)非線性方程組，當(dāng)軸承的基本尺寸參數(shù)給定時(shí)，對應(yīng)一組外部載荷Fr、Fa和M，可以通過牛頓迭代法求解得到各個(gè)接觸對的法向趨近量[2]。該非線性方程的求解需要編譯程序求解，詳細(xì)的求解過程參考文獻(xiàn)[3]。

3 偏航變槳軸承在ANSYS中的模擬方式

link10單元具有獨(dú)一無二的雙線性剛度矩陣特性[4]，使其成為一個(gè)僅受拉或僅受壓的桿單元，這里設(shè)置成僅受壓的特性來模擬滾子和滾道的接觸屬性：滾子和滾道接觸時(shí)傳遞法向的接觸載荷，滾子和滾道分離時(shí)不傳遞載荷。通過一組link10單元來模擬一個(gè)滾子，設(shè)置合理的link10單元參數(shù)，就能很好地模擬該類軸承的非線性接觸特性，以及接觸角度變化和接觸剛度的變化特性。

本文采用7根桿單元模擬一個(gè)接觸對，如圖2所示，一個(gè)滾子有2個(gè)接觸對，共需要用14根link10單元模擬一個(gè)滾子，在45°和135°方向，滾子與滾道剛好接觸，在其他方向，滾子與滾道存在間隙，間隙值的大小可以通過滾子和滾道的幾何關(guān)系求解。

圖2 軸承滾子link10單元模擬示意圖

Link10單元可以通過賦予初始應(yīng)變來模擬滾子和滾道的間隙。不同角度上link10單元的初始應(yīng)變不一樣，在45°和135°方向上link10單元的初始應(yīng)變?yōu)?，見表1。

表1 link10單元參數(shù)取值參數(shù)表

對于模擬一個(gè)滾子的link10單元組，隨著接觸載荷的增加，滾道之間相互趨近，承載作用力的link10單元數(shù)量增加，接觸的剛度增大，符合滾球和滾道的接觸特點(diǎn)。通過合理設(shè)置link10單元的剛度（通過橫截面積大小實(shí)現(xiàn)）和初始應(yīng)變，使link10單元組和軸承滾道滾子接觸對具有相同的非線性接觸剛度。在風(fēng)力發(fā)電機(jī)組有限元模擬時(shí)通過這種方式模擬軸承來保證軸承接觸剛度和傳遞載荷的精度。

本文作者通過多個(gè)軸承有限元計(jì)算與理論的比較，給出單個(gè)link10單元橫截面積Alink10的近似公式為：

式中：Kn為2個(gè)滾道間的載荷位移系數(shù)；L為滾子直徑；E為軸承彈性模量；Alink10為link10的橫截面積；δn為偏離45°第一個(gè)link10所在角度球和滾道的間隙。

根據(jù)幾何關(guān)系，初始間隙的計(jì)算公式為：

式中：f為滾道曲率半徑系數(shù)；θ為link10所處的角度；A為滾道滾子初始接觸時(shí)的溝心距；R為滾子半徑；δθ為初始接觸時(shí)各角度上滾子和滾道的間隙。

各個(gè)link10單元初始應(yīng)變εθ的計(jì)算為：

4 有限元和理論計(jì)算對比驗(yàn)證

這里通過對某3 MW風(fēng)力發(fā)電機(jī)組的偏航軸承進(jìn)行有限元模擬，該軸承有2排滾子，共177× 2=354個(gè)滾珠；滾子直徑45 mm；溝曲率系數(shù)0.525；滾道直徑為47.25 mm。

偏航軸承的有限元模型如圖3～4所示。

圖3 偏航軸承有限元模型

圖4 偏航軸承有限元模型截面細(xì)節(jié)

這里通過2種方式來驗(yàn)證有限元模型中模擬的偏航軸承的非線性特性，第一是驗(yàn)證軸承徑向載荷位移關(guān)系，第二是驗(yàn)證軸承在翻轉(zhuǎn)力矩作用下各個(gè)滾子接觸對的接觸載荷和接觸角度。

（1）軸承徑向載荷位移關(guān)系驗(yàn)證

有限元計(jì)算時(shí)約束偏航軸承的外圈，在內(nèi)圈上施加徑向載荷，分多個(gè)載荷步逐漸增加徑向載荷，查看內(nèi)圈的徑向位移。圖5為有限元計(jì)算與理論計(jì)算的徑向載荷與徑向位移曲線比較，從圖中可以看出通過有限元計(jì)算與理論計(jì)算，2條曲線基本一致。

圖5 偏航軸承徑向載荷與徑向位移曲線比較

（2）接觸載荷和接觸角度驗(yàn)證

有限元計(jì)算時(shí)約束偏航軸承的外圈，在內(nèi)圈上施加6 000 kNm的翻轉(zhuǎn)力矩，計(jì)算完后通過后處理查看link10單組合成的接觸載荷以及合成的接觸角，并與理論計(jì)算比較。

通過圖6和圖7可以看出，有限元模型計(jì)算和理論計(jì)算的結(jié)果非常吻合，用link10單元來模擬軸承滾珠和滾道的非線性接觸非常合適。這樣模擬能夠保證軸承傳遞載荷的正確性，并且能夠大大簡化模型、減小計(jì)算量。

圖6 偏航軸承滾道滾子接觸載荷比較

圖7 偏航軸承滾道滾子接觸角度比較

5 結(jié)論

在風(fēng)力發(fā)電機(jī)組主要部件有限元強(qiáng)度分析時(shí)，用link10單元代替軸承滾子，能夠正確模擬載荷的傳遞，準(zhǔn)確地考慮軸承由于非線性對風(fēng)力發(fā)電機(jī)組有限元分析的影響，并能減小模型規(guī)模，縮短計(jì)算時(shí)間，提高機(jī)組設(shè)計(jì)的效率和可靠性。

[1]哈里斯,等著.滾動軸承分析[M].羅繼偉,等譯.北京:機(jī)械工業(yè)出版社,2009

[2]羅繼偉,羅天宇.滾動軸承分析計(jì)算與應(yīng)用[M].北京:機(jī)械工業(yè)出版社,2010

[3]楊靜,廖暉.風(fēng)力機(jī)偏航變槳軸承載荷分布分析及數(shù)值求解[J].東方汽輪機(jī),2012,(1)：35-38

[4]王新敏,李義強(qiáng),許宏偉,編著.ANSYS結(jié)構(gòu)分析單元與應(yīng)用[M].北京:人民交通出版社,2011

Nonlinear Modeling of Yaw and Pitch Bearings in Wind Turbine Finite Element Analysis

Yang Jing，Liao Hui
（Dongfang Electric Wind Power Co.,Ltd.,Deyang Sichuan,618000）

Yaw and pitch bearings are the important part of wind turbine for transmitting loads.Loads transfer and distribution between ball and race way of yaw and pitch bearings must be properly considered in finite element analysis of wind turbine main components.In this paper the bearing ball is modeled using link10 element group in ANSYS.We can not only consider nonlinear characteristics of bearing,but also greatly simplify the calculation model,which makes the finite element strength analysis of wind turbine main components become fast and accurate.

wind turbine,finite element analysis,yaw bearing,pitch bearing,nonlinear

TK83

B

1674-9987（2015）04-0064-04

10.13808/j.cnki.issn1674-9987.2015.04.014

楊靜（1984-），男，工學(xué)碩士，工程師，2009年畢業(yè)于重慶大學(xué)工程力學(xué)專業(yè)，現(xiàn)從事風(fēng)力發(fā)電機(jī)設(shè)計(jì)研發(fā)工作。