顏強(qiáng) 魏靜 張世界 郭劍鵬

摘要：為精確對(duì)大型風(fēng)電機(jī)組進(jìn)行共振分析及優(yōu)化，通過(guò)研究不同縮聚點(diǎn)個(gè)數(shù)的柔性齒對(duì)系統(tǒng)模態(tài)的影響，指出了柔性輪齒縮聚的必要性和推薦數(shù)量，并建立了考慮主軸柔性和齒輪傳動(dòng)系統(tǒng)全柔性的傳動(dòng)鏈動(dòng)力學(xué)模型。根據(jù)求解的風(fēng)機(jī)傳動(dòng)系統(tǒng)固有頻率，提出共振甄別五大篩選原則，識(shí)別出系統(tǒng)潛在共振點(diǎn)，通過(guò)掃頻時(shí)/頻域分析，甄別出系統(tǒng)危險(xiǎn)共振點(diǎn)。根據(jù)共振頻率振型特點(diǎn)對(duì)齒輪箱進(jìn)行了模型優(yōu)化，依據(jù)共振點(diǎn)數(shù)量和風(fēng)電機(jī)組重量決定優(yōu)化方案，研究方法和結(jié)果可以為大型風(fēng)機(jī)多柔體動(dòng)力學(xué)高精度建模和共振甄別優(yōu)化設(shè)計(jì)提供參考。

關(guān)鍵詞：多柔體動(dòng)力學(xué)；大型風(fēng)機(jī)傳動(dòng)鏈；潛在共振點(diǎn)甄別；動(dòng)力學(xué)優(yōu)化設(shè)計(jì)

中圖分類(lèi)號(hào)： O313.7；TK83??? 文獻(xiàn)標(biāo)志碼： A??? 文章編號(hào)：1004-4523（2022）05-1157-08

DOI：10.16385/j .cnki .issn .1004-4523.2022.05.013

引言

隨著風(fēng)電機(jī)組的普及，降低其故障發(fā)生率已成為核心問(wèn)題。據(jù)統(tǒng)計(jì)，傳動(dòng)系統(tǒng)的主軸、齒輪箱等發(fā)生率較低的失效形式，維修時(shí)間往往比較長(zhǎng)，其中齒輪箱故障導(dǎo)致風(fēng)電機(jī)組停機(jī)的時(shí)間最長(zhǎng)[1]。風(fēng)電機(jī)組受到齒輪時(shí)變嚙合剛度、阻尼、嚙合誤差和軸承時(shí)變剛度等內(nèi)部激勵(lì)，以及隨機(jī)風(fēng)載和部件柔性變形等外部激勵(lì)的影響，可能會(huì)在某個(gè)轉(zhuǎn)速下發(fā)生劇烈的共振。一旦發(fā)生共振，將對(duì)傳動(dòng)系統(tǒng)中的共振部件產(chǎn)生不可逆轉(zhuǎn)的破壞，造成整個(gè)風(fēng)機(jī)的停機(jī)，因此，有必要在設(shè)計(jì)過(guò)程中進(jìn)行風(fēng)力發(fā)電機(jī)傳動(dòng)鏈危險(xiǎn)共振甄別，指導(dǎo)后期改進(jìn)設(shè)計(jì)。

目前風(fēng)機(jī)齒輪箱建模方法可概括為三類(lèi)：純扭轉(zhuǎn)模型、3，4或6自由度剛性多體模型和柔性多體模型[2]。純扭轉(zhuǎn)模型對(duì)系統(tǒng)做了大量的簡(jiǎn)化，而剛體模型忽略了部件的模態(tài)信息，不能完整地反映系統(tǒng)的動(dòng)態(tài)特征。近年來(lái)，柔性多體建模方法得到了學(xué)者的青睞。然而考慮的柔性體數(shù)量較少，常常忽略主軸、齒輪箱箱體或者齒輪的柔性[3?6]。由于風(fēng)力渦輪機(jī)尺寸的逐漸增大，更大的力和力矩集中在齒輪內(nèi)部，齒輪箱部件的柔性會(huì)影響整個(gè)傳動(dòng)鏈的響應(yīng)結(jié)果，導(dǎo)致計(jì)算結(jié)果不夠準(zhǔn)確，因此齒輪自身的柔性與齒輪副的嚙合柔性對(duì)風(fēng)電機(jī)組傳動(dòng)鏈振動(dòng)特性的準(zhǔn)確分析至關(guān)重要。此外，在傳動(dòng)鏈共振點(diǎn)方面的研究中，大多是利用系統(tǒng)的固有頻率求解坎貝爾圖，并結(jié)合系統(tǒng)能量分布來(lái)對(duì)傳動(dòng)鏈的共振點(diǎn)進(jìn)行甄別[7?8]，這種方法忽略了固有頻率的阻尼特性和速度特性，并不全面。

通過(guò)研究不同縮聚點(diǎn)個(gè)數(shù)的柔性齒對(duì)系統(tǒng)模態(tài)的影響，指出了柔性輪齒縮聚的必要性和推薦數(shù)量，并建立了考慮主軸柔性和齒輪傳動(dòng)系統(tǒng)全柔性的傳動(dòng)鏈動(dòng)力學(xué)模型；提出甄別系統(tǒng)潛在共振點(diǎn)的五大原則，更為準(zhǔn)確地甄別系統(tǒng)潛在共振點(diǎn)，根據(jù)時(shí)/頻域分析，得到系統(tǒng)的危險(xiǎn)共振點(diǎn)；根據(jù)共振頻率振型特點(diǎn)對(duì)齒輪箱進(jìn)行了模型優(yōu)化，依據(jù)共振點(diǎn)數(shù)量和風(fēng)電機(jī)組重量決定優(yōu)化方案，優(yōu)化后的傳動(dòng)系統(tǒng)不存在危險(xiǎn)共振點(diǎn)。研究方法和結(jié)果可以為大型風(fēng)機(jī)多柔體動(dòng)力學(xué)高精度建模和共振甄別優(yōu)化設(shè)計(jì)提供參考。

1 柔性多體動(dòng)力學(xué)理論

柔性多體系統(tǒng)動(dòng)力學(xué)理論主要以拉格朗日乘子法[9]和模態(tài)綜合法[10]為理論基礎(chǔ)。柔性體的運(yùn)動(dòng)包含剛體運(yùn)動(dòng)和自身的柔性變形兩部分，結(jié)合笛卡爾坐標(biāo)系、歐拉角以及模態(tài)坐標(biāo)可將其廣義坐標(biāo)表示為[11]：

式中ψ為旋進(jìn)角（又稱(chēng)進(jìn)動(dòng)角）：θ為章動(dòng)角；φ為自轉(zhuǎn)角（可以統(tǒng)稱(chēng)為歐拉角），q 為某個(gè)柔性體模態(tài)坐標(biāo)；為慣性參考系中的笛卡爾坐標(biāo)；為反映剛體方位的歐拉角；為模態(tài)坐標(biāo)。柔性體經(jīng)歷剛體平移運(yùn)動(dòng)和柔體變形后，其上任一位置的坐標(biāo)矢量為：

式中ro為動(dòng)坐標(biāo)系原點(diǎn)在慣性坐標(biāo)系中的位置向量；A 為將物理坐標(biāo)轉(zhuǎn)換為慣性坐標(biāo)的關(guān)系矩陣；si為節(jié)點(diǎn)i發(fā)生變形之前的坐標(biāo)位置；φi表示點(diǎn)i平移運(yùn)動(dòng)的模態(tài)矩陣。

對(duì)式（2）求導(dǎo)即得到該節(jié)點(diǎn)的速度矢量：

式中ω表示角速度向量；B 為歐拉角求導(dǎo)后與角速度矢量轉(zhuǎn)換的關(guān)系矩陣；E 為單位矩陣；～表示各矢量的對(duì)稱(chēng)矩陣。

結(jié)合拉格朗日方程建立柔性體的動(dòng)力學(xué)方程：

式中ψ為約束方程；λ為拉格朗日乘子；ξ為廣義坐標(biāo)；Q 為投影到ξ上的廣義力；L 為拉格朗日項(xiàng)；Γ為能量消耗系數(shù)。柔性體的運(yùn)動(dòng)微分方程為：

2 風(fēng)機(jī)傳動(dòng)鏈動(dòng)力學(xué)建模

2.1? 系統(tǒng)整體建模

將完整的傳動(dòng)鏈結(jié)構(gòu)拆分為若干個(gè)子結(jié)構(gòu)，其中包括齒輪、軸系、箱體三大種類(lèi)的部件，其基本參數(shù)如表1所示。利用 MPC（Multi?Point Constraints）建立相應(yīng)的力元、約束和風(fēng)機(jī)傳動(dòng)鏈動(dòng)力學(xué)分析模型，如圖1所示。

軸承剛度是齒輪箱內(nèi)最重要的參數(shù)之一，在全柔體模型中，利用6×6的軸承剛度矩陣來(lái)建立齒輪與軸系及箱體之間的耦合關(guān)系。軸承的滾動(dòng)體與內(nèi)外圈之間的耦合關(guān)系為：

式中：

式中Fbody，i為外界載荷矩陣（i=1代表齒輪，i=2代表軸承）；Kbearing為軸承剛度矩陣；Cbearing為軸承阻尼矩陣；qbody，i為變形量；kaxial表示軸承軸向剛度；kradial表示軸承徑向剛度；ktilt1和 ktilt2表示沿徑向方向的轉(zhuǎn)動(dòng)剛度。

利用 MPC 連接點(diǎn)，通過(guò)施加齒輪力元、花鍵力元和軸承力元建立子部件的連接關(guān)系，組建系統(tǒng)的耦合動(dòng)力學(xué)模型，其三維模型如圖2所示。

2.2 輪齒建模研究

Helsen等[12]表明由于剛性 MPC 中縮聚點(diǎn)與主節(jié)點(diǎn)之間是剛性連接，這種連接方式額外增加了有限元模型的剛度，這將大幅度地削弱傳動(dòng)鏈零部件的柔性。而柔性 MPC 利用下式所示的位移加權(quán)算法將縮聚點(diǎn)的載荷分配到主節(jié)點(diǎn)上，允許主節(jié)點(diǎn)（外部/界面節(jié)點(diǎn)）之間存在相對(duì)變形，保留了子結(jié)構(gòu)連接界面的柔性，因此以下研究均采用柔性多點(diǎn)約束。

式中? Fi 為作用在縮聚點(diǎn)上的力；M 為縮聚點(diǎn)受到的扭矩；ωi為權(quán)重系數(shù)；ri為主節(jié)點(diǎn)i到縮聚點(diǎn)的距離；n 為連接界面上的主節(jié)點(diǎn)數(shù)量。

根據(jù)每個(gè)零部件之間的運(yùn)動(dòng)關(guān)系，分別在各部件創(chuàng)建柔性多點(diǎn)約束（MPC）。為了將齒輪輪齒部分建立為柔性體，在每個(gè)輪齒節(jié)圓位置沿齒寬方向創(chuàng)建均勻分布的柔性 MPC，如圖3所示。根據(jù) Guy? an[13]提出的模型自由度靜態(tài)縮減方法可知，由于模型忽略了從自由度的慣性力，使得計(jì)算精度隨著模態(tài)頻率的增大而降低，但這類(lèi)誤差可以通過(guò)增加縮聚點(diǎn)的數(shù)目來(lái)彌補(bǔ)。

以單級(jí)斜齒齒輪傳動(dòng)系統(tǒng)為研究對(duì)象，建立單級(jí)齒輪傳動(dòng)系統(tǒng)全柔性傳動(dòng)鏈動(dòng)力學(xué)模型，保持模型和加載方式不變，改變輪齒處縮聚點(diǎn)數(shù)，分別計(jì)算輪齒處不同縮聚點(diǎn)數(shù)對(duì)系統(tǒng)模態(tài)的影響，基本參數(shù)如表2所示。

2.3 結(jié)果與分析

各階頻率隨輪齒縮聚點(diǎn)數(shù)變化如表3所示。由表3可知，齒輪副的耦合作用會(huì)派生出新的頻率，例如第3階，其縮聚點(diǎn)數(shù)從1變化至5時(shí)，第3階頻率從397.43 Hz 減小到381.65 Hz，整體變化顯著。圖 4為第3階振型圖，由圖4可知，隨著縮聚節(jié)點(diǎn)數(shù)的增加，輪齒的柔性更強(qiáng)。由表3可以看出，系統(tǒng)頻率中各軸系單獨(dú)的頻率不受輪齒處縮聚節(jié)點(diǎn)的影響，如耦合系統(tǒng)第6階頻率為主動(dòng)輪軸第13階頻率。

結(jié)果表明，耦合系統(tǒng)中與單軸系頻率相近的頻率基本不受輪齒縮聚點(diǎn)數(shù)影響；耦合系統(tǒng)派生頻率中輪齒處變形較大的頻率對(duì)輪齒縮聚點(diǎn)數(shù)敏感性較大，隨著縮聚點(diǎn)數(shù)增加呈減小趨勢(shì)。在大型傳動(dòng)鏈模型中，由于誤差的累積與系統(tǒng)的耦合作用，精度會(huì)進(jìn)一步下降。因此為了準(zhǔn)確地表示輪齒變形，避免由簡(jiǎn)化輪齒柔性導(dǎo)致的振型缺失或不完整，建議輪齒處縮聚點(diǎn)數(shù)設(shè)置為3個(gè)，本文風(fēng)電機(jī)組模型中輪齒處均采用3個(gè)縮聚點(diǎn)建模。

3 傳動(dòng)鏈系統(tǒng)動(dòng)力學(xué)分析

3.1? 系統(tǒng)潛在共振點(diǎn)甄別

3.1.1 頻率原則

系統(tǒng)的內(nèi)、外部激擾頻率及其倍頻會(huì)與系統(tǒng)的固有頻率形成共振，在分析時(shí)，將激擾頻率的上限定為最高轉(zhuǎn)速軸轉(zhuǎn)頻的6倍頻或者最高轉(zhuǎn)速齒輪的嚙頻3倍頻[14]。

根據(jù)風(fēng)機(jī)傳動(dòng)鏈齒輪傳動(dòng)系統(tǒng)的設(shè)計(jì)參數(shù)，模型中最高轉(zhuǎn)速軸的軸頻六倍頻為60 Hz；最高轉(zhuǎn)速齒輪的嚙頻三倍頻為780 Hz 。為了將盡可能多的潛在共振點(diǎn)考慮在內(nèi)，本文選取嚙頻三倍頻780 Hz 為頻率篩選的上限頻率，即考慮到傳動(dòng)系統(tǒng)固有頻率的294階，如表4所示。

3.1.2 阻尼原則

阻尼原則主要是通過(guò)對(duì)比各階次固有頻率對(duì)應(yīng)的阻尼比范圍，達(dá)到縮減固有頻率排查范圍的目的。阻尼比ζ的數(shù)值范圍及代表的物理意義如表5所示。根據(jù)機(jī)械零部件的物理及機(jī)械性質(zhì)，無(wú)阻尼、臨界阻尼和過(guò)阻尼這三種振動(dòng)狀態(tài)實(shí)際并不存在，因此其對(duì)應(yīng)的固有頻率是無(wú)意義的，應(yīng)該排除。僅留下阻尼比滿(mǎn)足0＜ζ＜1的固有頻率階次，從而達(dá)到了縮小甄別范圍的目的。

根據(jù)表4中各階模態(tài)固有頻率的阻尼比進(jìn)行判斷，可以將第10階之前的頻率階次排除在外，系統(tǒng)還剩余285階固有頻率。

3.1.3 能量原則

第10階到18階頻率系統(tǒng)各部件分別在 x，y 和 z的平移運(yùn)動(dòng)方向及其對(duì)應(yīng)的轉(zhuǎn)動(dòng)方向α，β和γ的振動(dòng)能量占比，如圖5所示。各階次的模態(tài)總能量在六個(gè)方向上分布并不相同，不同階次模態(tài)能量的主要分布方向也并不相同。而傳動(dòng)鏈扭轉(zhuǎn)方向的阻尼比較小，吸振能力差，使傳動(dòng)鏈易于發(fā)生扭轉(zhuǎn)振動(dòng)，影響機(jī)組的安全穩(wěn)定運(yùn)行[11]，因此γ方向的模態(tài)能量需要重點(diǎn)關(guān)注。

能量篩選原則概括為：在某階次頻率下傳動(dòng)系統(tǒng)中的所有部件繞旋轉(zhuǎn)方向的能量之和小于1，則該階次頻率不會(huì)激起系統(tǒng)的共振，故該頻率可以排除[15]，即：

式中Enγ為某階次頻率下傳動(dòng)系統(tǒng)中的第 n 個(gè)子結(jié)構(gòu)繞旋轉(zhuǎn)方向的能量。將能量總和低于1的頻譜階次排除后，得到結(jié)果如表6所示的16個(gè)階次。

3.1.4 同一速度級(jí)原則

同一速度級(jí)原則的內(nèi)涵是：在某一個(gè)潛在共振點(diǎn)的系統(tǒng)頻率上，如果與激擾部件在同一速度級(jí)別的部件的模態(tài)能量值均小于20%，則該頻率不足以激起系統(tǒng)的共振。若激擾頻率為第一級(jí)行星輪系嚙頻的3倍頻，其同一速度級(jí)別的部件應(yīng)包括：第一級(jí)行星輪、行星軸和太陽(yáng)輪，如圖6所示。排除不滿(mǎn)足同一速度級(jí)原則的頻率階次后，得到滿(mǎn)足條件的剩余階次，如表7所示。

3.1.5 系統(tǒng) Campbell 圖

坎貝爾圖可以將系統(tǒng)固有頻率與激擾頻率曲線(xiàn)繪制在同一個(gè)坐標(biāo)系內(nèi)，以工作轉(zhuǎn)速為橫坐標(biāo)，激擾頻率為縱坐標(biāo)，通過(guò)尋找固有頻率與激擾頻率曲線(xiàn)的交點(diǎn)來(lái)確定系統(tǒng)的潛在共振點(diǎn)。8 MW 風(fēng)力發(fā)電機(jī)傳動(dòng)系統(tǒng)的切入轉(zhuǎn)速為2.57 r/min，切出轉(zhuǎn)速為21.43 r/min，額定轉(zhuǎn)速為9 r/min，依據(jù)表7的固有頻率繪制出的系統(tǒng)坎貝爾圖如圖7所示。

坎貝爾圖中，并不是所有的交點(diǎn)都是潛在共振點(diǎn)，結(jié)合表7，若某個(gè)交點(diǎn)的固有頻率和對(duì)應(yīng)的激勵(lì)源與表6中相應(yīng)的固有頻率和激勵(lì)源不一致，則該點(diǎn)并不是共振點(diǎn)，可以排除。最終找到6個(gè)潛在共振點(diǎn)，如表8所示，對(duì)應(yīng)的振型如圖8所示。表8中： f2m 表示第二級(jí)的嚙頻。

3.2? 掃頻時(shí)/頻域分析

上述五大篩選原則給出的是系統(tǒng)的潛在共振點(diǎn)，若要得到系統(tǒng)的危險(xiǎn)共振點(diǎn)，需要對(duì)其進(jìn)行升速掃頻分析。限于篇幅僅以表8第一個(gè)潛在共振點(diǎn)為例進(jìn)行分析，此點(diǎn)是固有頻率134.048 Hz 與齒輪箱第2級(jí)嚙頻的交點(diǎn)，能量主要集中在第1級(jí)齒圈，此時(shí)輪轂的轉(zhuǎn)動(dòng)速度為14.14 r/min，該轉(zhuǎn)速附近的角加速度時(shí)域圖通過(guò) FFT 變換。如圖9所示為該轉(zhuǎn)速附近第一級(jí)齒圈角加速度時(shí)/頻域圖，在135.6 Hz 處存在峰值，與該潛在共振點(diǎn)對(duì)應(yīng)頻率134.048 Hz 重合，故此潛在共振點(diǎn)是危險(xiǎn)共振點(diǎn)，會(huì)在響應(yīng)部件處產(chǎn)生共振危險(xiǎn)，需要重點(diǎn)關(guān)注。通過(guò)對(duì)其余潛在共振點(diǎn)進(jìn)行時(shí)/頻域分析，系統(tǒng)共存在兩個(gè)危險(xiǎn)共振點(diǎn)，該模型有發(fā)生共振的危險(xiǎn)，需要對(duì)其進(jìn)行優(yōu)化。

4? 系統(tǒng)模型優(yōu)化

傳動(dòng)鏈在經(jīng)過(guò)共振點(diǎn)甄別后，識(shí)別出2個(gè)危險(xiǎn)共振點(diǎn)，對(duì)應(yīng)系統(tǒng)的2階固有頻率，對(duì)應(yīng)的共振轉(zhuǎn)速分別為14.14和14.9 r/min，共振部件包括前箱體、一級(jí)齒圈和后箱體。由于各支撐剛度對(duì)固有頻率的影響并不明顯，而質(zhì)量對(duì)固有頻率的影響較大[16]，所以?xún)?yōu)化主要針對(duì)箱體結(jié)構(gòu)尺寸，在保證能夠?yàn)辇X輪系統(tǒng)提供穩(wěn)定支撐的基礎(chǔ)上，盡可能地降低傳動(dòng)鏈系統(tǒng)的共振風(fēng)險(xiǎn)。最簡(jiǎn)單的一種優(yōu)化方法是改變箱體壁厚；第二種優(yōu)化方法是根據(jù)共振對(duì)應(yīng)的振型，在適當(dāng)?shù)奈恢迷O(shè)置筋板，從而有針對(duì)性地阻止其發(fā)生共振變形。

基于已甄別出的共振點(diǎn)，第一種優(yōu)化方案為增加前箱體和后箱體在發(fā)生共振時(shí)變形最大的位置處的箱體壁厚以及齒圈外圈直徑，如圖10所示，將這三個(gè)部件的壁厚增加20%。第二種優(yōu)化方案考慮共振振型為一級(jí)齒圈變形和前后箱體整體擺動(dòng)，故在前箱體前端的上、下吊耳處一共設(shè)置8個(gè)筋板，在后箱體變形較大處共設(shè)置12個(gè)筋板，如圖11所示。

利用共振點(diǎn)甄別五大原則對(duì)優(yōu)化模型進(jìn)行潛在共振點(diǎn)甄別，通過(guò)掃頻時(shí)/頻域分析，方案一最終確定2個(gè)危險(xiǎn)共振點(diǎn)，方案二未發(fā)現(xiàn)危險(xiǎn)共振點(diǎn)，如表9所示。可見(jiàn)，在變形處增設(shè)筋板的優(yōu)化方法可以有效降低系統(tǒng)的共振風(fēng)險(xiǎn)。并且從表10中可以看出，采用方案一使總質(zhì)量增加了10.69%，遠(yuǎn)大于第二種方案的1.95%。

綜上所述，在變形處增設(shè)筋板的優(yōu)化方法在僅使齒輪箱增加了不到2%的重量的情況下，減少了危險(xiǎn)共振點(diǎn)的數(shù)量，是更優(yōu)的優(yōu)化方案。

5 結(jié)論

結(jié)合縮聚和多體動(dòng)力學(xué)理論，建立了某8 MW 大型風(fēng)力發(fā)電機(jī)的多柔體動(dòng)力學(xué)模型，考慮了主軸柔性和齒輪箱的全柔性。研究了輪齒處不同縮聚點(diǎn)數(shù)對(duì)系統(tǒng)模態(tài)的影響，得出耦合系統(tǒng)中與單軸系頻率相近的頻率基本不受輪齒縮聚點(diǎn)數(shù)影響；耦合系統(tǒng)派生頻率中輪齒處變形較大的頻率對(duì)輪齒縮聚點(diǎn)數(shù)敏感性較大，隨著縮聚點(diǎn)數(shù)增加呈減小趨勢(shì)。

利用模型求解了系統(tǒng)前294階固有頻率，利用五大篩選原則找到了系統(tǒng)的6個(gè)潛在共振點(diǎn)，通過(guò)時(shí)域分析和頻域分析確定2個(gè)危險(xiǎn)共振點(diǎn)。針對(duì)危險(xiǎn)共振點(diǎn)的振型和時(shí)頻域分析結(jié)果對(duì)系統(tǒng)模型進(jìn)行優(yōu)化，對(duì)比了改變箱體壁厚和在齒輪箱箱體變形位置增設(shè)筋板兩種優(yōu)化方案，增設(shè)筋板的方法優(yōu)化后，傳動(dòng)鏈沒(méi)有危險(xiǎn)共振點(diǎn)，模型優(yōu)化取得了優(yōu)良的效果。研究結(jié)果可以為大型風(fēng)機(jī)多柔體動(dòng)力學(xué)建模、共振點(diǎn)甄別和系統(tǒng)優(yōu)化等提供參考。

參考文獻(xiàn)：

[1] 金曉航，孫毅，單繼宏，等.風(fēng)力發(fā)電機(jī)組故障診斷與預(yù)測(cè)技術(shù)研究綜述[ J ].儀器儀表學(xué)報(bào)，2017，38（5）：1041-1053.

JinXiaohang，Sun Yi，Shan Jihong，et al . Fault diagno? sis? and? prognosis? for? wind? turbines ：an? overview [ J ].Chinese Journal of Scientific Instrument，2017，38（5）：1041-1053.

[2] 劉波，賀志佳，金昊.風(fēng)力發(fā)電現(xiàn)狀與發(fā)展趨勢(shì)[ J ].東北電力大學(xué)學(xué)報(bào)，2016，36（2）：7-13.

LIU? Bo，HE Zhijia，JIN? Hao . Wind power status and development trends[ J ]. Journal of Northeast Dianli Uni? versity，2016，36（2）：7-13.

[3] 何玉林，黃偉，李成武，等.大型風(fēng)力發(fā)電機(jī)傳動(dòng)鏈多柔體動(dòng)力學(xué)建模與仿真分析[ J ].機(jī)械工程學(xué)報(bào)，2014，50（1）：61-69.

HE Yulin，HUANG Wei，LI Chengwu，et al . Flexible multibody dynamics modeling and simulation analysis of large-scale? wind? turbine? drivetrain [ J ]. Journal? of? Me ? chanical Engineering，2014，50（1）：61-69.

[4] 張智偉，左扣成，周明明，等.基于RomaxWIND的大兆瓦級(jí)風(fēng)機(jī)傳動(dòng)鏈非線(xiàn)性動(dòng)力學(xué)分析[ J ].機(jī)械傳動(dòng)，2014，38（10）：58-61.

Zhang Zhiwei，ZuoKoucheng，Zhou Mingming，et al . Nonlinear dynamics analysis of large wind turbine drive? train based on RomaxWIND[ J ]. Journal of Mechanical Transmission，2014，38（10）：58-61.

[5] 李浪.大型風(fēng)力發(fā)電機(jī)組整機(jī)動(dòng)力學(xué)建模與性能預(yù)測(cè)分析[ D ].重慶：重慶大學(xué)，2016.

Li? Lang . Dynamics? modeling? and? performance? predic? tion analysis of large-scale wind turbine[D ]. Chongqing： Chongqing University，2016.

[6] 周林，鄭四發(fā)，連小珉.加速工況下傳動(dòng)系統(tǒng)扭轉(zhuǎn)振動(dòng)分析[ J ].振動(dòng)工程學(xué)報(bào)，2010，23（6）：601-605.

ZHOU? Lin ，ZHENG? Sifa，LIAN? Xiaomin . Modeling and research? on torsional vibration? of transmission? sys? tem? under? speeding-up? condition[ J ]. Journal? of Vibra ? tion Engineering，2010，23（6）：601-605.

[7] 杜靜，秦月，李成武.風(fēng)力發(fā)電機(jī)組傳動(dòng)鏈動(dòng)力學(xué)建模與仿真分析[ J ].太陽(yáng)能學(xué)報(bào)，2014，35（10）：1950-1957.

Du? Jing ，Qin? Yue ，Li? Chengwu . Dynamics? modeling and? simulation? analysis? of wind? turbine? drive? train[ J ].Acta? Energiae? Solaris? Sinica，2014，35（10）： 1950-1957.

[8] 秦月.風(fēng)力發(fā)電機(jī)組傳動(dòng)鏈動(dòng)力學(xué)建模與仿真分析研究[ D ].重慶：重慶大學(xué)，2014.

Qin? Yue . Analysis? and? research? on? the? wind? turbine drive? train’s? dynamic? modeling? and? simulation [ D ]. Chongqing：Chongqing University，2014.

[9]  SIMEON? B . On? lagrange? multipliers? in? flexible? multi?body dynamics[ J ]. Computer Methods in Applied Me? chanics? and? Engineering ，2006，195（50-51）： 6993-7005.

[10] CRAIG Jr R R . Review of time-domain and frequency-domain? component-mode? synthesis method[ J ]. Ameri? can Society of Mechanical Engineers，Applied Mechan? ics Division，AMD，1985，67：1-30.

[11]陳爽.兆瓦級(jí)風(fēng)力發(fā)電機(jī)組傳動(dòng)鏈動(dòng)態(tài)特性研究[ D ].重慶：重慶大學(xué)，2013.

Chen Shuang . Dynamic characteristics research of mega? watt? level? wind? turbine? drive? train [ D ]. Chongqing： Chongqing University，2013.

[12] Helsen? J ，Vanhollebeke? F ，Marrant? B ，et? al . Multi?body modelling of varying complexity for modal behav? iour analysis of wind turbine gearboxes[ J ]. RenewableEnergy，2011，36（11）：3098-3113.

[13] GUYAN R J . Reduction of stiffness and mass matrices[ J ]. AIAA Journal，1965，3（2）：380.

[14]顏誠(chéng).基于模態(tài)綜合法2 MW 風(fēng)電增速箱柔體動(dòng)力學(xué)分析及優(yōu)化設(shè)計(jì)[ D ].湘潭：湘潭大學(xué)，2015.

Yan? Cheng . Flexible body danamic simulation research for 2 MW wind turbine with component mode synthesis theory[ D ]. Xiangtan：Xiangtan University，2015.

[15]繆炳榮，羅仁，王哲，等. SIMPACK 動(dòng)力學(xué)分析高級(jí)教程[M].成都：西南交通大學(xué)出版社，2010.

Miao Bingrong，Luo Ren，Wang Zhe，et al . Advanced Tutorial on SIMPACK Dynamics Analysis[M]. Cheng? du：Southwest Jiaotong University Press，2010.

[16]魏靜，呂程，孫偉，等. NGW 型行星輪系模態(tài)特性及參數(shù)敏感性研究[ J ].振動(dòng)工程學(xué)報(bào)，2013，26（5）：654-664.

WEI Jing，Lü? Cheng，SUN Wei，et al . Study on the mode characteristics and parameter sensitivity for NGW planetary gear system[ J ]. Journal of Vibration Engineer? ing，2013，26（5）：654-664.

Modeling of large -scale wind turbine drive train and methods of resonance identification and model optimization

YAN Qiang，WEI Jing，ZHANG Shi-jie，GUO Jian-peng

（State Key Laboratory of Mechanical Transmission，Chongqing University，Chongqing 400044，China）

Abstract： To analyze and optimize a large-scale wind turbine accurately，this paper considers the flexibility of the main shaft and the full flexibility of the gear transmission system，and establishes a flexible multi-body dynamic model of the drive chain of an 8 MW wind turbine . The effects of number of condense points on the gear system modal are also studied . The natural frequencies of the system are obtained，and 6 dangerous resonance points are identified by applying five screening principles . Based on the charac ? teristics of the resonance points，the gearbox model is optimized by comparing two methods . The number of resonance points is re? duced to zero by optimization . The results can provide a reference for the dynamic modeling and optimal design of large-scale wind turbines

Key words ： flexible multi-body dynamics；large-scale wind turbine drive chain；potential resonance point identification；dynamic optimization and designing

作者簡(jiǎn)介：顏強(qiáng)（1997―），男。電話(huà)：18728193776；E-mail：slmt _yanqiang@foxmail .com。

通訊作者：魏靜（1978―），男，教授，博士生導(dǎo)師。電話(huà)：13629752837；E-mail：weijing _slmt＠163.com。