風(fēng)力發(fā)電機葉片振動的分析與處理

風(fēng)力發(fā)電機葉片振動的分析與處理

葉片作為風(fēng)機的重要部件，其短而長的柔性結(jié)構(gòu)與風(fēng)機的正常運行密切相關(guān)。而作用于葉片上的載荷包括離心力、氣動載荷和慣性載荷等都具有隨機性和交變性,必然會引起振動,以至于風(fēng)力發(fā)電機不能平穩(wěn)運行甚至?xí)l(fā)生斷裂破壞。ZhangJP等用流固耦合方法對葉片在氣動載荷作用下進行了穩(wěn)定性分析,ZhangH等通過理論計算得出了風(fēng)力機葉片模態(tài)頻率值大小,郭龍、康順等用模態(tài)分析法對風(fēng)力機葉片進行了動力學(xué)分析計算,單光坤等用振動試驗方法對風(fēng)力機葉片進行疲勞特性分析。由于風(fēng)力機葉片復(fù)雜的結(jié)構(gòu)形式與工作時環(huán)境因素的影響,葉片的研究還有諸多未解決的問題,其中有關(guān)來流風(fēng)速和轉(zhuǎn)速對葉片振動特性影響的研究相對較少。模擬風(fēng)力發(fā)電機的實際運行狀況,利用ANSYS有限元分析軟件對葉片進行模態(tài)和諧響應(yīng)分析,研究葉片的振動頻率和振型的變化規(guī)律;葉片在不同來流風(fēng)速與轉(zhuǎn)速的作用下所受的預(yù)應(yīng)力對振動特性的影響規(guī)律,及葉片在最大外載荷作用下的受迫振動規(guī)律。1固有振動特性對葉片振動特性進行有限元求解,首先要對其結(jié)構(gòu)模型進行有限元離散,再對離散了的每個單元節(jié)點求解振動方程。在此基礎(chǔ)上對單元組合進行整體結(jié)構(gòu)計算,加入適當?shù)倪吔鐥l件,采用直接解法或迭代解法進行求解,得到結(jié)構(gòu)振動方程的解?？紤]葉片在沒有預(yù)應(yīng)力作用下的固有振動特性,振動方程為式中:[M]為結(jié)構(gòu)質(zhì)量矩陣;[K]為結(jié)構(gòu)剛度矩陣;{U″}為節(jié)點加速度矢量;{U}為節(jié)點位移矢量。由于固有頻率和振型是系統(tǒng)的固有屬性,與外部激勵無關(guān),因此式(1)右端的載荷項可以忽略。對于線性多自由度系統(tǒng),自由振動為簡諧振動式中:{Φ}i為第i(i=1,2,…,n)階固有頻率對應(yīng)的特征向量,即振型;ωi為第i階固有頻率;t為時間。將式(2)代入式(1),得到其特征方程為由式(4)可解得特征值ωi2。將ωi2代入式(3),可求得對應(yīng)的特征向量{φ}i。當風(fēng)速引起的氣動力或者轉(zhuǎn)速引起的的離心力作為穩(wěn)態(tài)載荷作用于葉片上時,此時可以把外載荷對葉片的作用力作為預(yù)應(yīng)力加在葉片上。有預(yù)應(yīng)力的振動方程為2網(wǎng)格劃分模型建立風(fēng)力發(fā)電機葉片的振動特性分析模型包括:建立葉片三維實體模型,確定葉片所選用的材料和需要劃分網(wǎng)格的單元類型,進行網(wǎng)格劃分以及設(shè)置其邊界條件等。1工程軟件geomagic12對某小型風(fēng)力機葉片實物進行三坐標激光掃描儀掃描,將掃描獲得的葉片勻點圖利用逆向工程軟件GEOMAGIC12進行處理,使其實體化,轉(zhuǎn)成三維實體模型,如圖1所示。葉片內(nèi)部為實心結(jié)構(gòu),長度680mm,厚度10mm,寬度最大140mm,寬度最小51mm,攻角最大35°,攻角最小0.07°。2葉片有限元模型風(fēng)力機葉片材料為GFRP塑料,其材料特性如下:密度為1.95×103kg/m3、展向彈性模量為1.93×104MPa、剪切彈性模量為2.88×104MPa、泊松比為0.15、材料阻尼為0.48,疲勞強度為20MPa?？紤]到葉片復(fù)雜的扭曲形狀,網(wǎng)格劃分采用帶中間節(jié)點的四面體單元,單元類型選取SOLID186。網(wǎng)格劃分后,得到8203個網(wǎng)格單元,16647個節(jié)點。3風(fēng)輪齒輪中心平面位置依據(jù)風(fēng)力發(fā)電機葉片在工作中的安裝情況,在有限元分析、簡化時,依據(jù)右手笛卡爾直角坐標系定則,把坐標原點O設(shè)置在風(fēng)輪輪轂中心,X軸方向與風(fēng)輪旋轉(zhuǎn)周向一致,Y軸方向與風(fēng)輪半徑方向一致,Z軸方向與風(fēng)輪旋轉(zhuǎn)平面垂直。為了與風(fēng)力機實際運行情況接近,在邊界條件設(shè)置時,對葉片的葉根部分進行約束,葉根部分的迎風(fēng)面和背風(fēng)面限制其Y和Z方向的平動自由度,葉根部分的其余兩側(cè)面限制X方向的平動自由度。圖2為建立的葉片有限元分析模型。3葉片振動特性分析對風(fēng)力發(fā)電機葉片的模態(tài)分析包括靜頻特性與動頻特性分析,研究葉片的固有振動頻率和振型的變化規(guī)律,以及在預(yù)應(yīng)力的作用下葉片的振動頻率的變化規(guī)律。3.1振動方式與振型的關(guān)系靜頻特性分析是對振動系統(tǒng)在不受外力作用下的振動特性分析。在葉片的葉根部分施加約束后,通過ANSYS有限元仿真計算出葉片的固有頻率和振型,提取前6階模態(tài)響應(yīng),其各階頻率和相應(yīng)的葉片最大變形值如表1所示。仿真結(jié)果表明:葉片的振動形式主要有3種:揮舞、擺振和扭振。其中,揮舞是葉片在垂直于旋轉(zhuǎn)平面方向上的彎曲振動,擺振是葉片在旋轉(zhuǎn)平面內(nèi)的彎曲振動,扭振是繞葉片徑向的扭轉(zhuǎn)振動。圖3為葉片發(fā)生振動時前6階頻率下的振型變化情況。從圖3可以看出:葉片的第1階振型以揮舞振動為主,葉片的2階、5階振型以擺陣振動為主,葉片的3階振型為扭振與擺陣的耦合振動,葉片的第4階、6階振型主要是扭轉(zhuǎn)振動。依據(jù)振動理論,葉片振動過程的能量主要集中于1、2階頻率處,而1、2階振型為揮舞和擺陣振動,扭轉(zhuǎn)振動發(fā)生在高階振型處,因此揮舞和擺陣為葉片的主要振動形式,是引起葉片疲勞破壞的主要原因。3.2轉(zhuǎn)速對葉片振動頻率的影響葉片在工作過程中受到的外載荷主要是縱向載荷,包括來流風(fēng)速對葉片產(chǎn)生的氣動壓力和葉片在旋轉(zhuǎn)過程中產(chǎn)生的離心力等載荷。如果穿過風(fēng)力發(fā)電機的來流風(fēng)速均勻且恒定,則可認為葉片所受的風(fēng)壓為恒定載荷,因此,可把風(fēng)壓對葉片的作用力作為預(yù)應(yīng)力加在葉片上,來研究不同來流風(fēng)速時的葉片振動特性。如果風(fēng)力發(fā)電機葉片轉(zhuǎn)速恒定,則可認為葉片所受的離心力為恒定載荷,在模擬時可把葉片所受的離心力作為預(yù)應(yīng)力施加在葉片上,來研究不同轉(zhuǎn)速時的葉片振動特性。利用Fluent流體軟件在風(fēng)速v分別為4m/s、8m/s和12m/s時,對葉片所受到的氣動壓力進行數(shù)值模擬,獲得的葉片表面風(fēng)壓值p如表2所示,其中,葉片表面風(fēng)壓是指葉片在迎風(fēng)面和背風(fēng)面上的壓強差。將表2中的不同風(fēng)速v時的風(fēng)壓載荷p施加在葉片表面上,對葉片進行有預(yù)應(yīng)力的模態(tài)分析,獲得不同風(fēng)速下葉片前6階振動頻率如表3所示?？梢钥闯?在風(fēng)速v為4m/s時,葉片的各階振動頻率與固有振動頻率大小一致;在風(fēng)速v為8m/s時,葉片在第4階和第6階的振動頻率比固有頻率大0.01Hz;而在風(fēng)速v為12m/s時,葉片在第3階的振動頻率比固有頻率大0.001Hz,第4階頻率大0.01Hz,第6階頻率大0.02Hz。風(fēng)速在4m/s或者更小時,作用于葉片上的氣動壓力較小,即使把風(fēng)壓載荷加在葉片上,對葉片的固有頻率也沒有影響;而風(fēng)速在8~12m/s時,作用于葉片上的氣動壓力值有所增大,把此范圍內(nèi)的風(fēng)速對應(yīng)的風(fēng)壓載荷加在葉片上,對葉片的振動頻率影響也不大,只是在高階次的振動頻率有微小的增加;而葉片的振動主要發(fā)生在1階和2階頻率上,此時風(fēng)速對低階頻率并無影響,因此,來流風(fēng)速的變化對葉片振動頻率的影響很小。在轉(zhuǎn)速n分別為100r/min、150r/min和300r/min時,對葉片施加由于轉(zhuǎn)速引起的離心力載荷后,進行有預(yù)應(yīng)力模態(tài)分析,獲得不同轉(zhuǎn)速下葉片前6階振動頻率f,如表4所示?？梢钥闯?在轉(zhuǎn)速n為100r/min時,葉片的各階振動頻率比固有頻率都有所增大,而在轉(zhuǎn)速n為150r/min和300r/min時,葉片各階振動頻率增加的更多。這是因為,隨著葉片轉(zhuǎn)速的增加,葉片所受的離心力越大,從而引起葉片動力剛化效應(yīng)的增強,使得風(fēng)力機葉片的各階振動頻率隨著轉(zhuǎn)速的增加而增大。由于轉(zhuǎn)速對葉片振動頻率影響較大,引用模態(tài)頻率偏差系數(shù)ηm來分析轉(zhuǎn)速對葉片振動頻率的影響程度。模態(tài)頻率偏差系數(shù)ηm式中:fdm為考慮轉(zhuǎn)速的第m階動模態(tài)頻率;fjm為不考慮轉(zhuǎn)速的第m階靜模態(tài)頻率。圖4為葉片轉(zhuǎn)速n分別為100r/min、150r/min和300r/min時,不同的頻率階次m與模態(tài)頻率偏差系數(shù)ηm的關(guān)系曲線。由圖4知,葉片轉(zhuǎn)速n對葉片固有頻率的1階頻率影響最大,其次是3階頻率,而對第2、4、5、6階的頻率影響較小;隨著轉(zhuǎn)速n的增加,葉片的各階模態(tài)頻率偏差系數(shù)ηm也相應(yīng)的增大,如n為100r/min時,1階的ηm頻率偏差系數(shù)為0.56%,當n為300r/min時,1階頻率偏差系數(shù)ηm增加到4.9%。4振動頻率對葉片響應(yīng)的影響對葉片進行諧響應(yīng)分析是在模態(tài)分析的基礎(chǔ)上,計算出葉片的固有頻率后,確定其頻率范圍,再把旋轉(zhuǎn)葉片所受到的離心力和風(fēng)壓加在葉片上,通過計算葉片在周期載荷作用下,不同頻率時的頻率響應(yīng)值,找出葉片在受迫載荷作用下的最大位移和對應(yīng)的應(yīng)力。在進行葉片諧響應(yīng)分析時,葉片所施加的最大受迫載荷為來流風(fēng)速為12m/s的風(fēng)壓值和轉(zhuǎn)速為300r/min的慣性載荷。由于葉尖處的位移是葉片上位移響應(yīng)最大處,因此,在葉尖處取一個節(jié)點(節(jié)點號4371),研究其在振動頻率范圍內(nèi)因受迫載荷產(chǎn)生的位移響應(yīng)值。圖5給出了該節(jié)點處的X、Y、Z這3個方向的響應(yīng)位移s與頻率f的關(guān)系曲線,從圖中可以看出,葉片在最大受迫載荷作用下,振動頻率f從4Hz增加到20Hz時,該節(jié)點處的3個方向的響應(yīng)位移s都增加;頻率f再增加到60Hz時,該節(jié)點處的X與Z向的響應(yīng)位移s急劇下降,而Y向的響應(yīng)位移s逐步降低;因此可以進一步確定,葉片在最大受迫載荷作用下能夠激發(fā)起的振動頻率為20Hz,接近葉片固有頻率中的1階頻率(18.639Hz),此時該節(jié)點處的X、Y、Z這3向的響應(yīng)位移s分別為2.54mm、0.14mm、0.63mm。表明位移響應(yīng)s最大發(fā)生在X方向,其次發(fā)生在Z向,而Y向的位移響應(yīng)較小。在葉片振動頻率為20Hz時,對其施加最大受迫載荷,進行葉片整體位移和應(yīng)力的諧響應(yīng)分析,獲得的葉片振動響應(yīng)結(jié)果如圖6所示。從圖6a)的位移響應(yīng)云圖可知,葉片的響應(yīng)位移沿著葉根到葉尖的徑向方向逐漸增加,在葉尖處的位移最大,其值為2.605mm,在葉根處位移很小。從圖6b)的應(yīng)力響應(yīng)云圖可知,葉片在靠近葉根處的應(yīng)力最大,最大值為4.711MPa,而在葉片中間區(qū)域應(yīng)力值都在1.57MPa以下,遠低于葉片的疲勞強度(20MPa)。因此,葉片在風(fēng)速小于12m/s、轉(zhuǎn)速低于300r/min運行時,葉片的響應(yīng)位移和應(yīng)力都較小,不易發(fā)生疲勞破壞。5來流風(fēng)速對振動頻率的影響1)對葉片進行靜頻特性分析,得到了葉片的前6階固有頻率和振型,表明揮舞和擺陣振動為葉片的主要振動形式,是引起葉片疲勞破壞的主要原因。2)動頻特性分析結(jié)果表明:在0~12m/s范圍內(nèi)的來流風(fēng)速對葉片的振動頻率的影響很小;而葉片轉(zhuǎn)速在0~