風(fēng)力發(fā)電機偏航減速器的運動設(shè)計與仿真

設(shè)計

風(fēng)力發(fā)電機偏航減速器的運動設(shè)計與仿真

李明倫王益軒劉瑋黃新武張永軍

(西安工程大學(xué)機電工程學(xué)院，陜西710048)

摘要：為了實現(xiàn)風(fēng)電偏航減速器的大減速比，借助虛擬樣機技術(shù)，用Pro/E軟件建立了四級行星齒輪減速器三維模型，然后導(dǎo)入到ADAMS中進行仿真分析，計算出所有齒輪的角速度，通過分析證實了各級傳動比分配是合理的。同時，為了避免減速器行星架在工作過程中發(fā)生共振，運用ANSYS/Workbench15.0計算了減速器行星架的固有頻率及振型，并進行了分析，最后得出了偏航驅(qū)動減速器的工作頻率范圍。

關(guān)鍵詞：行星齒輪減速器；仿真分析；偏航系統(tǒng)；傳動比；模態(tài)分析；風(fēng)力發(fā)電機

中圖分類號：TS 103.133文獻標(biāo)志碼：A

收稿日期：2015—04—16

基金項目：本文由陜西省教育廳產(chǎn)業(yè)化培育項目資助(項目編號：2013JC17)。

Kinematic Design and Simulation of Yawing Reducer

for Wind Generator

Li Minglun, Wang Yixuan, Liu Wei, Huang Xinwu, Zhang Yongjun

Abstract：To realize larger reduction ratio of yawing reducer of wind generator and adopt of virtual prototype technique, 3D model of four-stage planetary reducer has been established with Pro/E software, then simulative analysis has been conducted with the model leading to ADAMS to work out angular velocities of total gears, it has been proved that various drive ratio distribution was reasonable. Meanwhile, natural frequencies and vibration shapes of gear carrier have been calculated and analyzed by applying ANSYS/Workbench15.0 in order that sympathetic vibration of reducer gear carrier wouldn't occur during operation, finally working frequency scope of yawing drive reducer has been determined.

Key words：planetary reducer; analogue analysis; yawing system; drive ratio; mode analysis; wind generator

偏航系統(tǒng)是風(fēng)力發(fā)電機的重要部分，偏航系統(tǒng)的靈活可控性可以使風(fēng)機準(zhǔn)確穩(wěn)定對風(fēng)，使風(fēng)電機組能最大程度利用風(fēng)能，提升發(fā)電效率。減速器是偏航系統(tǒng)至關(guān)重要的一個組件， 行星減速器在風(fēng)力發(fā)電機上應(yīng)用較多，為了提高行星減速器的可靠性能及其使用年限，同時為了克服常規(guī)解析分析方法計算量大、耗時長、效率低等不足，運用有限元分析軟件對減速器的關(guān)鍵件行星架進行模態(tài)分析并計算其固有頻率及陣型是有必要的。

本文利用虛擬樣機技術(shù)，在Pro/E中建立了四級行星齒輪減速器各個零件的三維模型并裝配，通過Pro/E和ADAMS軟件之間的接口將模型導(dǎo)入ADAMS中，添加運動副和驅(qū)動之后進行仿真，得到各級行星輪和太陽輪、輸入軸和輸出軸的角速度曲線，通過分析數(shù)據(jù)和曲線，驗證了設(shè)計的合理性。運用ANSYS/Workbench15.0對其行星架進行了動態(tài)特性分析。

1風(fēng)力發(fā)電機的三維模型

風(fēng)力發(fā)電機組由風(fēng)輪和發(fā)電機兩大部分組成。風(fēng)輪的功能是將風(fēng)能轉(zhuǎn)換為機械能，發(fā)電機的功能是將機械能轉(zhuǎn)換為電能。風(fēng)力發(fā)電機一般由風(fēng)輪、發(fā)電機、偏航裝置、控制系統(tǒng)、塔架等主要部件所組成。如圖1所示。

2偏航減速器傳動系統(tǒng)的選擇及參數(shù)的擬定

2.1傳動系統(tǒng)的選擇

本文設(shè)計的風(fēng)力發(fā)電機組偏航減速器的傳動機構(gòu)如圖2所示，由四級NGW型行星齒輪傳動串聯(lián)而成。由于大功率直驅(qū)風(fēng)力發(fā)電機組整體結(jié)構(gòu)受力較多，所處場所的運行條件一般都很惡劣，故要求偏航減速器所用的材料要滿足機械強度的要求和極端條件下的材料特性。減速器各個傳動齒輪的材料及熱處理工藝見表1。

2.2各構(gòu)件參數(shù)的擬定及傳動比的計算

在多級行星齒輪傳動中，高速級傳動比的范圍一般為3.15～9；減速器低速級的傳動比一般為4～5.6；而中間級傳動比的范圍常取5～7.1。根據(jù)風(fēng)電機機型和功率，一般額定功率為500 kW～3 000 kW的風(fēng)電機組，偏航轉(zhuǎn)速為(0.1～0.5)r/min。已知本文設(shè)計的減速器輸入轉(zhuǎn)速為750 r/min，取偏航轉(zhuǎn)速為0.428 r/min，需要實現(xiàn)的總傳動比為1 752.3∶1，故采用四級減速。設(shè)計各級齒輪具體參數(shù)及尺寸見表2。

1—風(fēng)輪　2—發(fā)電機　3—偏航系統(tǒng)

采用反轉(zhuǎn)法，即：根據(jù)物體的相對運動原理，假使讓整個周轉(zhuǎn)輪系加上公共角速度，那么行星架可認為固定不動了，于是，周轉(zhuǎn)輪系便轉(zhuǎn)化成了定軸輪系。一級行星齒輪系統(tǒng)(見圖3)，各個構(gòu)件轉(zhuǎn)化之前和之后的角速度見表3。

本文中ω3為零，即行星輪內(nèi)齒圈是固定的，則傳動比i1H為

由公式(1)可得，第一級行星齒輪傳動比i14=第一級行星內(nèi)齒圈/第一級太陽輪+1=126/18+1=8(下標(biāo)14為圖2行星齒輪傳動系統(tǒng)簡圖中所對應(yīng)的構(gòu)件)，按照此方法推得第二級齒輪傳動比i47=117/21+1=6.57，第三級i7·10=120/24+1=6，第四級i10·14=123/27+1= 5.56，總傳動比i=i14i47i7·10i10·14=1 753.4。經(jīng)理論計算可得，行星齒輪減速器第一級(即高速級)傳動比i1=8，中間級的第二級傳動比和第三級傳動比i2=6.57、i3=6，第四級(即低速級)傳動比i4=5.56，均在要求的傳動比范圍內(nèi)，故該行星齒輪減速器各級減速比分配合理，理論上能實現(xiàn)總減速比為1 752.3∶1，接下來通過Pro/E軟件建立四級行星齒輪減速器三維模型，并運用ADAMS軟件平臺對減速器虛擬仿真，驗證傳動比的合理性。

3偏航四級行星齒輪減速器的建模及動力學(xué)仿真分析

3.1偏航四級行星齒輪減速器的三維模型

利用Pro/E參數(shù)化建模功能，按照表1中的數(shù)據(jù)建立偏航四級行星齒輪減速器中各個零件的三維模型，然后根據(jù)裝配原則進行裝配，見圖4。

表1　齒輪的材料及其熱處理工藝

表2　各級行星輪系參數(shù)

第一級 第二級 第三級 第四級

圖3　一級行星齒輪系統(tǒng)簡圖

構(gòu)件原角速度轉(zhuǎn)化輪系齒輪1ω1ωH1=ω1-ωH齒輪2ω2ωH2=ω2-ωH齒輪3ω3ωH3=ω3-ωH支架4ω4=0ωH4=ω4-ωH=-ωH行星架HωHωHH=ωH-ωH=0

3.2ADAMS平臺上為模型添加約束條件

利用Pro/E和ADAMS之間良好的數(shù)據(jù)交換功能，將減速器的三維模型導(dǎo)入到 ADAMS/View 中，文中減速器各個零件之間的運動副主要包括齒輪副、旋轉(zhuǎn)副和固定副。減速器采用四級行星輪系減速，每一級所添加的約束都相同，各級之間用固定副連接。

第一級所加的約束(見圖5)有：三個行星輪與行星架之間添加旋轉(zhuǎn)副；分別在內(nèi)齒圈與行星輪、太陽輪與行星輪之間添加齒輪副，共6個；3個行星輪分別與行星輪旋轉(zhuǎn)軸之間添加旋轉(zhuǎn)副；在行星架與大齒圈、太陽輪與行星架之間施加旋轉(zhuǎn)副；行星齒輪內(nèi)齒圈與減速器固定外殼之間施加固定約束。其余各級均按此依次添加約束，所有約束添加結(jié)束所得減速器整機模型如圖6所示。

3.3為模型添加驅(qū)動并分析數(shù)據(jù)

在ADAMS為行星齒輪減速器選擇添加旋轉(zhuǎn)驅(qū)動，設(shè)定的電機轉(zhuǎn)速為750 r/min。即角速度ω=4500°/s，設(shè)置仿真時間為5 s，仿真步長為500，仿真后的結(jié)論如下：

(1)減速器能夠在ADAMS實現(xiàn)運動仿真，各部件都能夠協(xié)調(diào)勻速旋轉(zhuǎn)，裝配正確，完全不存在干涉現(xiàn)象，故說明該減速器結(jié)構(gòu)設(shè)計合理。

(2)通過ADAMS/Processor(后處理模塊)獲取各主要部件角速度曲線，即第一級太陽輪角速度為4 500°/s，第一級的行星架角速度為562.5°/s，第二級太陽輪角速度562.5°/s(因為第一級行星架和第二級太陽輪間設(shè)置了固定副，故角速度相同)，第二級行星架角速度為85.7°/s，第三級太陽輪角速度為85.7°/s(第二級行星架和第三級太陽輪設(shè)置固定副)，第三級行星架角速度為14.3°/s，第四級太陽輪角速度為14.3°/s(第三級行星架和第四級太陽輪設(shè)置固定副)，第四級行星架角速度為2.57°/s，即行星減速器輸出軸角速度為2.57°/s，即偏航角速度為2.57°/s(即0.428 r/min)。

圖4　Pro/E平臺上的四級行星齒輪減速器三維模型

1—輸出軸　2—行星架　3—內(nèi)齒圈

圖6　ADAMS平臺上加約束副之后的減速器

(3)整理計算出減速器各級傳動比i1=8,i2=6.56,i3=5.99,i4=5.57，總傳動比是i=1 751，與理論設(shè)計總傳動比1 752.3的總誤差小于1%，因此，該減速器傳動比分配合理。

4ANSYS/Workbench 15.0模態(tài)分析

4.1行星架模態(tài)分析理論計算

行星架模態(tài)計算就是將行星架三維實體分散為一個多自由度系統(tǒng)，并且計算該系統(tǒng)中每一個自由度的模態(tài)及陣型，并通過線性幾何的方法疊加這些參數(shù)，對于n自由度系統(tǒng)，其物理參數(shù)的模型需要要通過n個相互獨立的物理坐標(biāo)來表達，在該系統(tǒng)的物理坐標(biāo)系中n個主振動疊加為該系統(tǒng)自由振動的響應(yīng)，每一種主振動都會相對為一種特定形態(tài)的自由振動，主振動的振動形態(tài)即為該系統(tǒng)的模態(tài)，而該系統(tǒng)的主振動頻率即系統(tǒng)的固有頻率。

一般結(jié)構(gòu)系統(tǒng)動力學(xué)分析的基本方程為：

式中，{F(t)}是力向量；[K]是剛度矩陣；[C]是阻尼矩陣；[M]是質(zhì)量矩陣；{x＂}是用物理坐標(biāo)描述的加速度列陣；{x′}是用物理坐標(biāo)描述的速度列陣；{x}是用物理坐標(biāo)描述的位移列陣。

當(dāng)系統(tǒng)的作用力為零時,該系統(tǒng)將處于自由振動狀態(tài),則：

當(dāng)行星架結(jié)構(gòu)發(fā)生自由振動的時候，在無阻尼時它們之間的頻率存在如下關(guān)系:

而實際中行星架阻尼比往往較小，即為λ<20%，對行星齒輪系統(tǒng)振動的影響可忽略不計,忽略后可得到該系統(tǒng)無阻尼自由振動方程：

設(shè)公式(5)的特解為：

將公式(6)代入方程(5)，得：

當(dāng)矩陣{?}為非零時, ([K]-μ2[M])為方程(7)的系數(shù)矩陣，{?}為特征向量，μ2為特征值。如果需要方程(7)有非零解，那么該方程系數(shù)矩陣的行列式必須等于零，方可滿足條件，即:

對方程(8)求解可得μ的n個根 μ0e(e=1，2，…n)，即該行星架的固有頻率。代入方程(7)得到矩陣{?}的n-1個獨立方程組,一共能求解得到n個線性無關(guān)的解,稱為行星架固有振型。其表達式為：

將特征矢量{?i}按列排成一個n階方陣，即為:

公式(10)即為系統(tǒng)的模態(tài)矢量，是最終要求解的模態(tài)矩陣。

因為矩陣[M]和矩陣[K]都為正定矩陣,故{?i}關(guān)于[M][K]正交，即表達為:

4.2行星架模態(tài)分析

模態(tài)分析用來確定目標(biāo)結(jié)構(gòu)的固有頻率和模態(tài)振型，同時要避免該結(jié)構(gòu)在外加載荷的作用下出現(xiàn)共振。本文運用ANSYS/Workbench15.0 對四級行星齒輪減速器行星架進行模態(tài)分析，確定其固有頻率及其振型，并計算出減速器各個部件的轉(zhuǎn)動頻率以及齒輪的嚙合頻率，然后和行星架固有頻率相比較，進而驗證結(jié)構(gòu)設(shè)計的合理性。

如圖7所示，在Pro/E軟件中建立行星架的三維模型，然后將在Pro/E中建立的實體模型導(dǎo)入ANSYS/Workbench15.0軟件中設(shè)置單位，為模型添加材料為structuralstell，采用自動劃分網(wǎng)格，生成六面體單元，并設(shè)置行星架在無約束下(自由狀態(tài)下)求解前二十階模態(tài)，通過ANSYS/Workbench15.0計算出行星架前二十階固有頻率，各階頻率見表4，其中4階的模態(tài)振型見圖8～圖11。

根據(jù)分析結(jié)果，可得知行星架的動態(tài)特性，在工作中盡量使工作頻率遠離求解的該系統(tǒng)的固有頻率，禁止機器在運行中發(fā)生共振現(xiàn)象，避免不必要的的人員傷亡及財產(chǎn)損失。由表4可以看出行星架前6階頻率幾乎都是0，這是由于前6階模態(tài)分別代表6個自由度的剛體模態(tài)，沿軸x、y、z方向的移動和轉(zhuǎn)動均不受約束，因此頻率等于0，第七階模態(tài)的固有頻率才是我們最終所得出的結(jié)果。經(jīng)過計算得出齒輪的嚙合振動頻率、各級齒輪的扭轉(zhuǎn)頻率，和該四級行星齒輪減速器各運動部件的工作頻率，可得知太陽輪與行星輪、行星輪與內(nèi)齒圈的嚙合頻率遠低于行星架的第七階固有頻率1 388.6Hz，故該結(jié)構(gòu)設(shè)計合理，能保證減速器正常運行，不至于在運行過程中發(fā)生共振現(xiàn)象。

表4　前二十階固有頻率

圖7　 Pro/ E平臺上行星架三維模型

圖8　行星架第七階振型圖

圖9　行星架第十階振型圖

圖10　行星架第十五階振型圖

圖11　行星架第二十階振型圖

5結(jié)論

行星齒輪結(jié)構(gòu)比較復(fù)雜，特別是多級行星齒輪系統(tǒng)在設(shè)計制造時過程繁瑣。本文利用Pro/E軟件建立了四級行星齒輪減速器三維模型，在ADAMS中對行星齒輪減速器仿真分析，并針對仿真結(jié)果進行有效分析，驗證了理論傳動比分配是合理的，基本符合設(shè)計要求。并通過ANSYS/Workbench15.0對行星架進行分析計算，計算出結(jié)構(gòu)的各級固有頻率，可防止結(jié)構(gòu)在外加載荷影響下，發(fā)生共振，進一步改善減速器性能，為2MW大型風(fēng)力發(fā)電機偏航減速器的設(shè)計提供了參考依據(jù)。

參考文獻

[1]鄂加強，張彬，董江東，等．新型風(fēng)電偏航減速機動力學(xué)仿真分析[J].中南大學(xué)學(xué)報：自然科學(xué)版，2011，42(8):2324-2331.

[2]彭清東.多級行星齒輪減速器傳動比的優(yōu)化分配[J].中國高新技術(shù)企業(yè),2008(9):82-83.

[3]Chen,J.L.,Tseng.DesignConsiderationsforAssembleinPlanetaryGearTrains[J〗.ProceedingsoftheInternationalConferenceonGearing,Transmissions,andMechanicalSystems.NotinghamTrentUniversityUK.Conference.2000．

[4]李善飛.偏航行星減速機優(yōu)化設(shè)計研究[D].大連理工大學(xué)碩士學(xué)位論文,2008：52-60.

[5]郝則勝. 基于精確加載的大功率增速器行星架結(jié)構(gòu)拓撲優(yōu)化設(shè)計[J]. 機械設(shè)計與制造,2008(7)：61-62.

[6]趙麗娟,陳令國,劉紅梅.礦用減速器行星架的有限元分析[J]. 煤礦機械,2007,28(1):51-52.

[7]濮良貴,紀(jì)名剛.機械設(shè)計[M].西北工業(yè)大學(xué)機械原理及機械零件教研室,2005:186-200.

[8]饒振綱.行星齒輪傳動設(shè)計[M].北京：化學(xué)工業(yè)出版社，2003:54-78.

[9]朱孝錄.齒輪傳動設(shè)計手冊[S].北京：化學(xué)工業(yè)出版社,2005:190-210.

[10]劉銀虎,繆炳祺.多體動力學(xué)仿真軟件ADAMS理論基礎(chǔ)及其功能分析[J].電子與封裝,2005(4):25-28.

[11]行星齒輪裝配原理，施洪生，滕健.SolidWorks中行星齒輪傳動的快速裝配[J〗.中國制造業(yè)信息化，2006(6):71-73.

[12]李軍.2MW風(fēng)力發(fā)電機組塔架結(jié)構(gòu)分析研究[J].太原理工大學(xué)學(xué)報，2011(05):43-45.

[13]陳旭，王益軒，李端，等.大型龍門數(shù)控纖維纏繞/鋪放機動態(tài)分析[J].西安工程大學(xué)學(xué)報，2013(2):216.

編輯傅冬梅