單動葉安裝角異常時軸流風(fēng)機(jī)壓力脈動特性的數(shù)值分析

葉學(xué)民，　丁學(xué)亮，　李春曦

(華北電力大學(xué) 電站設(shè)備狀態(tài)監(jiān)測與控制教育部重點(diǎn)實(shí)驗(yàn)室， 河北保定 071003)

?

單動葉安裝角異常時軸流風(fēng)機(jī)壓力脈動特性的數(shù)值分析

葉學(xué)民，丁學(xué)亮，李春曦

(華北電力大學(xué) 電站設(shè)備狀態(tài)監(jiān)測與控制教育部重點(diǎn)實(shí)驗(yàn)室， 河北保定 071003)

摘要：以O(shè)B-84型動葉可調(diào)軸流風(fēng)機(jī)為對象，對單動葉安裝角在不同偏離度下的風(fēng)機(jī)內(nèi)部流場特征和性能進(jìn)行了非定常計(jì)算，得到了安裝角異常下壓力脈動的時域和時頻特性、靜壓時均值及等效應(yīng)力分布.結(jié)果表明：動葉安裝角異常時，各監(jiān)測點(diǎn)的壓力信號及壓力脈動能量均出現(xiàn)異常特征，尤其是葉輪出口處的變化最為突出；動葉異常偏離對時域和時頻特性的影響程度隨偏離度的增大而提高，并導(dǎo)致壓力脈動強(qiáng)度增大；各監(jiān)測點(diǎn)的近似熵在大偏離度下變化顯著，對動葉異常偏離具有良好的表征能力；靜壓時均值的分布和數(shù)值隨偏離度改變均有較大變化，而葉片等效應(yīng)力僅在數(shù)值上有明顯改變.

關(guān)鍵詞：動葉可調(diào)軸流風(fēng)機(jī)； 偏離度； 安裝角異常； 壓力脈動

動葉可調(diào)軸流風(fēng)機(jī)因其結(jié)構(gòu)簡單緊湊、運(yùn)行范圍廣及變工況效率高等特點(diǎn)而廣泛用于航空航天、能源、動力和交通等領(lǐng)域.但其內(nèi)部結(jié)構(gòu)復(fù)雜，變工況運(yùn)行時，常因執(zhí)行機(jī)構(gòu)故障引發(fā)動葉安裝角異常偏離的現(xiàn)象，造成風(fēng)機(jī)不規(guī)則振動，從而影響風(fēng)機(jī)的安全運(yùn)行[1].研究表明，該不規(guī)則振動主要是因動葉發(fā)生異常偏離后、內(nèi)部非定常流場及壓力脈動的變化所致[2].因此，深入研究動葉發(fā)生異常偏離時風(fēng)機(jī)內(nèi)部的壓力脈動特性及其內(nèi)部流場分布至關(guān)重要.

目前，對于軸流風(fēng)機(jī)的研究多集中在內(nèi)部流場特性與結(jié)構(gòu)優(yōu)化[3-5]、失速與喘振[6-8]、葉頂間隙對風(fēng)機(jī)性能的影響[9-11]及噪聲[12-14]等方面，并取得了豐厚的研究成果.關(guān)于通風(fēng)機(jī)非定常流場和壓力脈動方面，陳慶光等[15]利用大渦模擬方法通過對旋軸流風(fēng)機(jī)的非定常計(jì)算，分析了各監(jiān)測點(diǎn)下的壓力脈動頻率、強(qiáng)度及成因.Ballesteros-Tajadura等[16]模擬了離心風(fēng)機(jī)葉輪及機(jī)殼處的壓力脈動，得到了壓力時域和頻譜分布特性，其計(jì)算結(jié)果與實(shí)驗(yàn)值吻合良好，但其研究僅分析了壓力脈動特性，未對內(nèi)部流場分布進(jìn)行研究.Jiro等[17]通過修正離心力、重力和其他主要誤差的影響，準(zhǔn)確測得某軸流風(fēng)機(jī)的瞬時壓力特征，但其風(fēng)機(jī)結(jié)構(gòu)較為簡單，所得結(jié)論不具有普遍性.

對于動葉安裝角異常時的軸流風(fēng)機(jī)，葉學(xué)民等[18-19]通過對兩級動葉可調(diào)軸流風(fēng)機(jī)的數(shù)值研究，得到了第一級動葉發(fā)生異常偏離及兩級動葉均異常偏離對風(fēng)機(jī)性能和內(nèi)部流場的影響，但其分析基于定常計(jì)算結(jié)果.在非定常計(jì)算方面，Posson等[20]利用寬頻噪聲模型分析了動葉偏移時的聲功率分布特征.李春曦等[21]通過非定常計(jì)算得到軸流風(fēng)機(jī)內(nèi)聲壓信號的時域和頻域特性，分析了軸流風(fēng)機(jī)內(nèi)氣動噪聲的分布，但對于風(fēng)機(jī)異常偏離時壓力的非定常特性和局部壓力脈動特性尚未展開研究.

鑒于目前動葉可調(diào)軸流風(fēng)機(jī)的非穩(wěn)態(tài)流動研究尚不完善，尤其是對動葉異常偏離時風(fēng)機(jī)內(nèi)部的壓力脈動特性尚未完全認(rèn)識.為此，筆者采用Fluent軟件對OB-84型動葉可調(diào)軸流風(fēng)機(jī)開展非定常計(jì)算，進(jìn)而分析動葉異常偏離對壓力時域和時頻特性、靜壓時均值及等效應(yīng)力分布的影響，為風(fēng)機(jī)安全穩(wěn)定運(yùn)行及優(yōu)化設(shè)計(jì)提供參考依據(jù).

1計(jì)算模型

以O(shè)B-84型單動葉可調(diào)軸流風(fēng)機(jī)為研究對象，其結(jié)構(gòu)包括集流區(qū)、動葉區(qū)、導(dǎo)葉區(qū)和擴(kuò)壓區(qū).該風(fēng)機(jī)有14片動葉、15片導(dǎo)葉，葉輪直徑為1 500 mm，葉頂間隙為4.5 mm，轉(zhuǎn)速為1 200 r/min，由此可得葉片通過頻率為280 Hz，風(fēng)機(jī)結(jié)構(gòu)和參數(shù)詳見文獻(xiàn)[22].實(shí)際運(yùn)行中，動葉安裝角正向偏離現(xiàn)象較常見，因此筆者對單動葉安裝角偏離度Δβ=0°～50°進(jìn)行研究，Δβ=0°表示動葉安裝角處于正常狀態(tài)(此時動葉安裝角為32°)，不同偏離度表示該動葉沿翼型軸線逆時針旋轉(zhuǎn)相應(yīng)的度數(shù)，如圖1所示.

圖1　動葉安裝角偏離度

基于動葉可調(diào)軸流風(fēng)機(jī)的結(jié)構(gòu)特點(diǎn)，采取分區(qū)和局部加密劃分方法，在動葉區(qū)因流動梯度及曲率變化大，故采用加密網(wǎng)格;而集流區(qū)、導(dǎo)葉區(qū)和擴(kuò)壓區(qū)采用稀疏網(wǎng)格，以減少計(jì)算量.模擬中，為驗(yàn)證網(wǎng)格無關(guān)性，選取163萬、197萬、235萬和260萬4組網(wǎng)格數(shù)進(jìn)行計(jì)算.模擬顯示，計(jì)算精度和計(jì)算時間均隨網(wǎng)格數(shù)的增加而提高，當(dāng)網(wǎng)格數(shù)為235萬時，可同時滿足計(jì)算時間和計(jì)算精度的要求，故選取該網(wǎng)格數(shù)對風(fēng)機(jī)性能進(jìn)行模擬，此時動葉區(qū)和導(dǎo)葉區(qū)的網(wǎng)格數(shù)分別為135萬和21萬.

定常計(jì)算采用Realizablek-ε湍流模型，該模型因引入旋轉(zhuǎn)和曲率相關(guān)項(xiàng)，而廣泛用于旋轉(zhuǎn)均勻剪切流、射流、混合流的自由流動和帶有分離的流動[18-19,21].同時，選用多參考系模型解決葉輪中葉片與機(jī)殼間動靜干涉面數(shù)據(jù)傳遞的問題.計(jì)算中，風(fēng)機(jī)入口邊界條件為入口速度，出口邊界條件為自由流出，不同區(qū)域的交界面上采用interface實(shí)現(xiàn)各區(qū)域交界面上流體參數(shù)的交換.模擬中，當(dāng)各參數(shù)殘差均小于10-4，且進(jìn)、出口截面的總壓均不隨時間改變時，則視計(jì)算已收斂.

將上述計(jì)算所得流場作為非定常計(jì)算的初場，由于滑移網(wǎng)格計(jì)算時只需要計(jì)算交界面兩側(cè)的通量且使其相等，即可實(shí)現(xiàn)動靜干涉面間的信息傳遞，故非定常計(jì)算選用滑移網(wǎng)格模型進(jìn)行動靜干涉面的數(shù)據(jù)傳遞.因風(fēng)機(jī)內(nèi)的氣流在動葉區(qū)獲得機(jī)械能，在導(dǎo)葉區(qū)改變其流動方向，故動葉區(qū)和導(dǎo)葉區(qū)氣流的壓力脈動比集流區(qū)和擴(kuò)壓區(qū)氣流的壓力脈動更為顯著.為此，筆者主要研究動葉區(qū)和導(dǎo)葉區(qū)壓力信號的時域和時頻特性.經(jīng)驗(yàn)證同一個徑向截面監(jiān)測點(diǎn)的壓力特性相似，為此每個截面只取一個具代表性的監(jiān)測點(diǎn)進(jìn)行分析，將壓力信號監(jiān)測點(diǎn)分別布置在動葉入口(P1點(diǎn))、葉頂間隙(P2點(diǎn))、動葉出口(P3點(diǎn))和導(dǎo)葉出口(P4點(diǎn))，如圖2所示.

圖2　監(jiān)測點(diǎn)位置分布

2模擬結(jié)果與分析

2.1時域分析

在計(jì)算動葉異常情形前，首先對正常情形下的風(fēng)機(jī)進(jìn)行數(shù)值模擬，以驗(yàn)證數(shù)值模擬的準(zhǔn)確性.結(jié)果表明，在模擬流量范圍內(nèi)，所得全壓與實(shí)驗(yàn)數(shù)據(jù)相比，平均偏差為1.7%，保證了數(shù)值模擬的可靠性.

為研究不同偏離度下，壓力信號隨時間的變化特征，圖3給出了葉輪旋轉(zhuǎn)一個周期內(nèi)的靜壓時域分布.從圖3可以看出，正常情形下，不同監(jiān)測點(diǎn)的壓力信號隨時間變化均呈周期或類周期性波動.這是由于葉輪旋轉(zhuǎn)過程中動葉周期性擊打周圍空氣所致，每個周期具有14個明顯峰值點(diǎn)，其個數(shù)與動葉數(shù)一致.另外，沿氣流流動方向，P1～P4點(diǎn)靜壓時均值逐漸增大，脈動幅度卻呈減小趨勢，這是因?yàn)檠豍1～P4點(diǎn)所在截面，氣體不斷獲得機(jī)械能，進(jìn)而靜壓逐漸提高，由此使得靜壓時均值增大；同時，氣流逐漸遠(yuǎn)離提供能量的動葉，因而其脈動幅度逐漸減小.

從圖3還可以看出，異常情形下，動葉入口處的P1點(diǎn)峰值點(diǎn)個數(shù)未改變，在0.01～0.02 s間出現(xiàn)異常波動.在此區(qū)域內(nèi)，峰值處的靜壓值顯著增大、而谷值處的靜壓值則顯著降低，且Δβ越大，峰谷間的脈動幅度更加突出，當(dāng)Δβ=50°時，其脈動幅度約為正常情形下的11倍，而其他時間區(qū)的壓力脈動變化不明顯.葉頂間隙處的P2點(diǎn)在0.01 s附近的靜壓峰值異常增大，隨Δβ的增加異常波動區(qū)有所變化，其中，Δβ=10°時僅在0.01 s附近出現(xiàn)一異常峰值點(diǎn)，而Δβ=50°時其影響范圍則擴(kuò)大到0.015 s.動葉異常偏離對P3和P4點(diǎn)處靜壓值的影響整體呈現(xiàn)相似趨勢：正常情形下，靜壓值均為正值，而動葉異常偏離后，尤其是Δβ較大時，靜壓值為負(fù)值且脈動幅度較大，波動時間范圍也擴(kuò)至0.01～0.03 s，并呈“V”形分布.其原因是氣流通過異常動葉時，流動角改變造成沖角增大，使得動葉對通過該通道氣流的擊打強(qiáng)度增大，進(jìn)而造成峰谷間的壓力脈動幅度隨Δβ增大而增大，并對后續(xù)時間的壓力脈動產(chǎn)生明顯的影響.可見，動葉異常偏離后各監(jiān)測點(diǎn)的時域特性明顯不同，但壓力異常波動范圍相似，且其影響程度隨Δβ的增大而增強(qiáng).

2.2壓力脈動強(qiáng)度

為進(jìn)一步研究不同監(jiān)測點(diǎn)處的非定常流動特征，將瞬時壓力轉(zhuǎn)化成無量綱壓力系數(shù)Cp(t)=p(t)/0.5ρu22，其中ρ為流體密度，p(t)為某時刻t下的靜壓值，u2為葉輪外圓周速度.一個周期內(nèi)壓力系數(shù)的平均值和壓力脈動強(qiáng)度[23]由式(1)和式(2)確定：

(1)

(2)

其中，N為一個周期內(nèi)的數(shù)據(jù)長度.

圖4給出了各監(jiān)測點(diǎn)在不同偏離度下的壓力脈動強(qiáng)度.從圖4可以看出，正常情形下，P2點(diǎn)的壓力脈動強(qiáng)度遠(yuǎn)高于其他監(jiān)測點(diǎn).動葉異常偏離后，P2點(diǎn)的壓力脈動強(qiáng)度隨Δβ的變化并不明顯，而P1、P3和P4點(diǎn)的壓力脈動強(qiáng)度則呈升高趨勢.這是因?yàn)樵谌~頂間隙處，異常動葉對靜壓值的影響有限；而在動葉入口、動葉出口和導(dǎo)葉出口處，隨Δβ的增大，異常動葉造成當(dāng)?shù)仂o壓值發(fā)生顯著變化，導(dǎo)致其壓力脈動強(qiáng)度提高，以及由此對上下游流動的干涉作用逐漸增強(qiáng).其中，P3點(diǎn)的增幅最為顯著，而P1和P4點(diǎn)的壓力脈動強(qiáng)度隨Δβ的增幅相差不大.可見，動葉異常偏離對葉頂間隙處壓力脈動強(qiáng)度的影響較小，而對動葉出口處壓力脈動強(qiáng)度的影響最大.

2.3時頻分析

時頻分析能同時描述信號在不同時間和頻率的能量密度和強(qiáng)度.選用Choi-Williams分布描述信號的時頻能量分布，該分布可消除典型Wigner-Vill分布中的交叉干擾項(xiàng)，但計(jì)算量無明顯變化[24].

圖5為不同偏離度下P2點(diǎn)的時頻分布.由圖5可知，正常情形下，壓力脈動信號在280 Hz附近出現(xiàn)最高點(diǎn)，這與葉片通過頻率相一致，同時壓力脈動能量主要集中于200～400 Hz頻率段，在500～600 Hz(包括2倍葉片通過頻率)間存在較低的壓力脈動能量.動葉異常偏離后，各偏離度對時頻分布的影

(a) P1

(b) P2

(c) P3

(d)P4

圖3時域特性

Fig.3Characteristics of static pressure in the time domain

圖4　不同偏離度下的壓力脈動強(qiáng)度

響有所不同.Δβ=10°時，在0.01～0.02 s間的壓力脈動能量出現(xiàn)異常分布，這與圖3(b)中的異常波動時間區(qū)間相符，除在1倍及2倍葉片通過頻率附近，其他頻率段也出現(xiàn)壓力脈動能量值，但這些能量處于較低水平；Δβ=20°時，t=0.013～0.025 s壓力脈動能量顯著減小，而小于此區(qū)間的能量降幅相對較小，其他區(qū)間的能量分布變化不明顯；Δβ=30°時，在0.012 s處產(chǎn)生最大壓力脈動能量，同時在整個頻率區(qū)間內(nèi)均存在壓力脈動能量，t<0.012 s時的壓力脈動能量減小；Δβ=40°時，t=0.01～0.02 s內(nèi)的壓力脈動能量分布變化不大，仍為能量集中區(qū)，而其他波動時間區(qū)間的能量顯著減小；Δβ=50°時，與Δβ=40°相比，壓力脈動能量僅在數(shù)值上存在差異，其分布規(guī)律相似.總體上，壓力脈動能量最大值隨Δβ的增大呈增大趨勢，即偏離度越大，動葉異常偏離對時頻分布的影響越明顯.

(a) Δβ=0°

(b) Δβ=10°

(c) Δβ=20°

(d) Δβ=30°

(e) Δβ=40°

(f) Δβ=50°

2.4近似熵

近似熵算法最早由Pincus提出，其值能夠表征信號的復(fù)雜度，因此越復(fù)雜的信號近似熵越大.近似熵能根據(jù)較短數(shù)據(jù)得出比較穩(wěn)健的估計(jì)值，同時具有良好的抗噪和抗干擾能力，其表達(dá)式為：

(3)

式中：r為相似容限；m為模式嵌入維數(shù)；Φm(r)為壓力信號的對數(shù)平均值[25].

選取m=2，r=0.2S，S為原始數(shù)據(jù)的標(biāo)準(zhǔn)差[25].

圖6為各監(jiān)測點(diǎn)在不同流量系數(shù)下的近似熵分布.由圖6可知，各監(jiān)測點(diǎn)的近似熵呈現(xiàn)不同的特征.動葉入口處P1點(diǎn)的近似熵在正常情形下處于較低水平且隨流量系數(shù)的增大變化不明顯，而動葉異常偏離后，近似熵均明顯增大，同時隨Δβ的增大，近似熵的平均水平呈先增后減的趨勢.對于葉頂間隙處的P2點(diǎn)，其異常情形下的近似熵同樣顯著高于正常情形，但其近似熵總體隨Δβ的增大而增大.對于動葉出口處的P3點(diǎn)，其正常情形下的近似熵與Δβ=10°時相差不大，當(dāng)Δβ進(jìn)一步增大后，近似熵則顯著減小，其中，Δβ>20°時近似熵的降幅減小.對于導(dǎo)葉出口處的P4點(diǎn)，動葉異常偏離后的平均近似熵均小于正常情形，且隨Δβ的增大總體呈減小趨勢.

(a) P1

(b) P2

(c) P3

(d) P4

上述分析表明，P1、P2點(diǎn)處在正常情形下的近似熵小于異常情形，而P3、P4點(diǎn)則相反，這是由于近似熵大小與信號復(fù)雜度有關(guān).由圖3可知，P1、P2點(diǎn)在正常情形下的壓力信號呈周期性波動，其峰值、谷值隨時間變化不明顯，而動葉異常偏離后的壓力信號出現(xiàn)了異常波動，使壓力信號復(fù)雜度增加，因此近似熵變大；P3、P4點(diǎn)在正常情形下的峰值、谷值隨時間變化顯著，壓力信號復(fù)雜度較大，動葉異常偏離后的壓力信號呈較規(guī)律的“V”形分布，故動葉異常偏離后近似熵呈減小趨勢.可見，近似熵對動葉異常偏離有很好的辨識能力.

2.5靜壓時均值分布

靜壓時均值可反映靜壓隨時間的變化率及流道內(nèi)流動的非定常特征.圖7給出了葉輪中間截面上的靜壓時均值分布，葉輪旋轉(zhuǎn)方向如圖所示.由圖7可知，正常情形下，流道內(nèi)的靜壓時均值從壓力面到吸力面呈下降趨勢，在靠近葉頂處出現(xiàn)最小值.當(dāng)動葉異常偏離后(見圖7(b))，異常動葉壓力面?zhèn)鹊母邏簠^(qū)范圍擴(kuò)大，吸力面則在靠近輪轂處產(chǎn)生低壓區(qū).這是由于動葉異常偏離導(dǎo)致邊界層脫落的渦流與流道內(nèi)的主流相互作用增強(qiáng)，使得葉片壓力面和吸力面?zhèn)攘鞯纼?nèi)的靜壓時均值分布發(fā)生改變，而其他流道靜壓時均值變化不明顯.

當(dāng)Δβ進(jìn)一步增大(見圖7(c)～ 圖7(f))，異常動葉吸力面?zhèn)鹊牡蛪簠^(qū)范圍擴(kuò)大，且低壓區(qū)靜壓時均值也逐漸降低，在逆旋轉(zhuǎn)方向上與異常動葉相鄰葉片的流道也受到影響，尤其是Δβ=50°時異常動葉吸力面?zhèn)鹊恼麄€流道處于負(fù)壓水平，逆旋轉(zhuǎn)方向上與異常動葉相鄰葉片吸力面?zhèn)鹊牡蛪褐祻娜~片中上部移至中下部.可見，隨Δβ的增大，動葉異常偏離對靜壓時均值的影響程度提高.

(a)Δβ=0°(b)Δβ=10°

(c)Δβ=20°(d)Δβ=30°

(e)Δβ=40°(f)Δβ=50°

圖7葉輪中間截面上的靜壓時均值分布

Fig.7Contour of time-averaged static pressure on the middle surface

2.6等效應(yīng)力分布

以上分析表明，動葉發(fā)生異常偏離后改變了流道內(nèi)的壓力分布，從而導(dǎo)致葉片表面載荷也發(fā)生變化，因此，有必要分析動葉異常偏離對葉片等效應(yīng)力的影響.圖8給出了不同偏離度下的葉片等效應(yīng)力分布.由圖8(a)可知，正常情形下，最大等效應(yīng)力值集中于吸力面葉根靠近前緣處，且等效應(yīng)力沿葉高方向逐漸減小.這是由于葉根固定在輪轂上且截面積最大，葉輪轉(zhuǎn)動過程中承受整個葉片的離心力，因此應(yīng)力集中于葉片根部.Δβ=10°～30°時，如圖8(b)～圖8(d)所示，隨Δβ的增大，葉片等效應(yīng)力分布變化不明顯，但最大等效應(yīng)力值呈增大趨勢.Δβ=40°和50°時，如圖8(e)和圖8(f)所示，最大等效應(yīng)力值進(jìn)一步增大，且集中于異常葉片吸力面?zhèn)龋渌~片高等效應(yīng)力區(qū)范圍也有所擴(kuò)大.

在強(qiáng)度計(jì)算中，塑性材料的許用應(yīng)力[σ]=σs/ns，其中σs為材料的屈服極限，ns為塑性材料的安全系數(shù).對于彈性結(jié)構(gòu)加載靜力載荷的情況下，ns取值范圍為1.5～2.所研究風(fēng)機(jī)葉片的材料為鑄鋁ZL101，其屈服強(qiáng)度σs為180 MPa，ns取2，可得出葉輪的許用應(yīng)力值為90 MPa.圖8中最大等效應(yīng)力值為10.1 MPa，小于許用應(yīng)力值.因此，軸流風(fēng)機(jī)動葉異常偏離后仍滿足強(qiáng)度要求.

3結(jié)論

(1) 動葉異常偏離后各監(jiān)測點(diǎn)壓力信號的周期性或類周期性波動被破壞，均出現(xiàn)不同程度異常波動，且偏離度越大其影響程度越明顯；動葉出口處的壓力脈動呈“V”形分布，其影響時間也明顯擴(kuò)大；P1、P3和P4點(diǎn)壓力脈動強(qiáng)度隨偏離度的增大而提高，而P2點(diǎn)壓力脈動強(qiáng)度的變化較小.

(2) 正常情形下，P2點(diǎn)的最大壓力脈動能量值位于葉片通過頻率處；動葉異常偏離改變了P2點(diǎn)的時頻分布，最大壓力脈動能量值隨偏離度的增大而增大.動葉異常偏離后，各監(jiān)測點(diǎn)的近似熵均發(fā)生改變，尤其是在大偏離度下的變化更加顯著，表明近似熵對動葉發(fā)生異常偏離有良好的表征能力.

(3) 隨著偏離度的增大，葉輪中間截面上的靜壓時均值分布及數(shù)值均發(fā)生改變；當(dāng)Δβ>20°時，在逆旋轉(zhuǎn)方向上與異常動葉相鄰葉片的流道也受到影響；而葉片等效應(yīng)力分布隨偏離度的增大變化不大，但等效應(yīng)力值呈增大趨勢.

(a)Δβ=0°(b)Δβ=10°

(c)Δβ=20°(d)Δβ=30°

(e)Δβ=40°(f)Δβ=50°

圖8葉片等效應(yīng)力分布

Fig.8Equivalent stress distribution on the surface of blade

參考文獻(xiàn)：

[1]郭玉杰，任廣印，劉占輝，等．超臨界600 MW機(jī)組鍋爐一次風(fēng)機(jī)間歇性振動診斷[J]．熱力發(fā)電，2012，41(8)：56-58.

GUO Yujie, REN Guangyin, LIU Zhanhui,etal. Diagnosis of intermittent vibration on primary air fan of supercritical boiler[J]. Thermal Power Generation,2012,41(8):56-58.

[2]XU C, AMANO R S. Unsteady pressure field investigation of an axial fan-blade unsteady pressure field measurement[J]. International Journal of Rotating Machinery,2002,8(6): 385-395.

[3]LI Yang, LIU Jie, OUYANG Hua,etal. Internal flow mechanism and experimental research of low pressure axial fan with forward-skewed blades[J].Journal of Hydrodynamics,2008,20(3): 299-305.

[4]孟麗，林麗華，吳海英，等．軸流風(fēng)機(jī)后導(dǎo)葉三維數(shù)值優(yōu)化設(shè)計(jì)方法及其應(yīng)用[J]．風(fēng)機(jī)技術(shù)，2014(3)：48-52．

MENG Li, LIN Lihua, WU Haiying,etal. A 3-D numerical optimization design method and its application in axial fans with outlet guide vanes[J].Fan Technology,2014(3): 48-52.

[5]劉廣強(qiáng)，呂子強(qiáng)，郭宇，等．中低壓軸流風(fēng)機(jī)流型優(yōu)化的數(shù)值模擬與試驗(yàn)分析[J]．流體機(jī)械，2014，42(10)：1-5．

LIU Guangqiang, Lü Ziqiang，GUO Yu,etal. Numerical simulation of axial fan based on different circulation index of impeller[J]. Fluid Machinery,2014,42(10):1-5.

[6]桂良明，劉聰，秦明亮．動葉可調(diào)軸流雙級一次風(fēng)機(jī)喘振分析及對策[J]．流體機(jī)械，2012，40(8)：44-46．

GUI Liangming, LIU Cong, QIN Mingliang. Analysis and solution of surge of axial flow primary air fan with double stage adjustable movable blade[J].Fluid Machinery,2012,40(8):44-46.

[7]張磊，郎進(jìn)花，王松嶺．電站軸流風(fēng)機(jī)旋轉(zhuǎn)失速工況下的葉輪靜力特性研究[J]．動力工程學(xué)報(bào)，2015，35(5)：387-393．

ZHANG Lei, LANG Jinhua, WANG Songling. Static characteristics of the impeller of an axial flow fan under rotating stall conditions[J]. Journal of Chinese Society of Power Engineering,2015,35(5): 387-393.

[8]SALUNKHE P B, PRADEEP A M. Stall inception mechanism in an axial flow fan under clean and distorted inflows[J].Journal of Fluids Engineering,2010,132(12):121102.

[9]葉學(xué)民，李鵬敏，李春曦．葉頂開槽對軸流風(fēng)機(jī)性能影響的數(shù)值研究[J]．中國電機(jī)工程學(xué)報(bào)，2015，35(3)：652-659．

YE Xuemin, LI Pengmin, LI Chunxi. Numerical investigation of grooved blade tips effect on the performance of an axial flow fan[J]. Proceedings of the CSEE,2015,35(3):652-659.

[10]NHO Y C, PARK J S, LEE Y J,etal. Effects of turbine blade tip shape on total pressure loss and secondary flow of a linear turbine cascade[J].International Journal of Heat and Fluid Flow,2012, 33(1):92-100.

[11]蔡濤，周水清，王軍.加環(huán)控制軸流風(fēng)機(jī)葉頂間隙的分析與試驗(yàn)[J]．工程熱物理學(xué)報(bào)，2014，35(12)：2406-2411．

CAI Tao,ZHOU Shuiqing,WANG Jun. Analysis and test on the performance of the axial fan with a ring near the blade tip on the casing[J].Journal of Engineering Thermophysics,2014, 35(12): 2406-2411.

[12]HU Binbin, OUYANG Hua, WU Yadong,etal. Numerical prediction of the interaction noise radiated from an axial fan[J]. Applied Acoustics, 2013,74(4):544-552.

[13]LEE J W, KIM Y Y. Topology optimization of muffler internal partitions for improving acoustical attenuation performance[J]. International Journal for Numerical Methods in Engineering,2009,80(4):455-477.

[14]李業(yè)，周水清，王軍，等．轉(zhuǎn)速對彎掠軸流風(fēng)機(jī)氣動噪聲的影響分析[J]．工程熱物理學(xué)報(bào)，2014，35(1)：51-55．

LI Ye, ZHOU Shuiqing, WANG Jun,etal. Analysis of aerodynamic noise of a swept-curved axial flow fan with different rotate speed[J]. Journal of Engineering Thermophysics,2014,35(1):51-55.

[15]陳慶光，房菲，張振東，等．對旋風(fēng)機(jī)葉輪內(nèi)部流場的壓力脈動分析[J]．山東科技大學(xué)學(xué)報(bào)(自然科學(xué)版)，2011，30(2)：81-85．

CHEN Qingguang, FANG Fei, ZHANG Zhendong,etal. Analysis on pressure fluctuation of the interior flow field in impellers of a contra-rotating axial fan[J].Journal of Shandong University of Science and Technology(Natural Science),2011,30(2):81-85.

[16]BALLESTEROS-TAJADURA R,VELARDE-SUAREZ S, HURTADO-CRUZ J P,etal. Numerical calculation of pressure fluctuations in the volute of a centrifugal fan[J].Journal of Fluids Engineering,2006,128(2):359-369.

[17]JIRO F, TAKUYA F, YUSUKE O,etal. Measurements of minute unsteady pressure on three-dimensional fan[J].Journal of Environment and Engineering,2010,5(2): 298-310.

[18]葉學(xué)民，李新穎，李春曦．第一級葉輪單動葉安裝角異常對動葉可調(diào)軸流風(fēng)機(jī)性能的影響[J].中國電機(jī)工程學(xué)報(bào)，2014，34(14)：2297-2306．

YE Xuemin, LI Xinying, LI Chunxi. Effect of the first-stage impeller with abnormal blade on the performance of an axial fan[J]. Proceedings of the CSEE,2014,34(14):2297-2306.

[19]葉學(xué)民，王豐，李春曦．兩級動葉異常偏離下的軸流風(fēng)機(jī)性能和內(nèi)流特征[J]．動力工程學(xué)報(bào)，2015，35(5)：394-403．

YE Xuemin, WANG Feng, LI Chunxi. Performance and internal flow characteristic of an axial flow fan with installation deviation in two-stage blades[J].Journal of Chinese Society of Power Engineering,2015,35(5):394-403.

[20]POSSON H, MOREAU S, ROGER M. Broadband noise prediction of fan outlet guide vane using a cascade response function[J]. Journal of Sound and Vibration,2011,330(25):6153-6183.

[21]李春曦，林卿，葉學(xué)民．單動葉安裝角異常時軸流風(fēng)機(jī)的噪聲特性[J]．中國電機(jī)工程學(xué)報(bào)，2015，35(5):1183-1192．

LI Chunxi, LIN Qing, YE Xuemin. Acoustic characteristics of an axial flow fan with abnormal installation angle of single blade[J]. Proceedings of the CSEE,2015,35(5):1183-1192.

[22]索洛瑪霍娃 T C．通風(fēng)機(jī)氣動略圖和性能曲線[M]．北京：煤炭工業(yè)出版社，1986：336-348．

[23]王文杰，袁壽其，裴吉，等．時序效應(yīng)對導(dǎo)葉式離心泵內(nèi)部壓力脈動影響的數(shù)值分析[J]．機(jī)械工程學(xué)報(bào)，2015，51(4)：185-192．

WANG Wenjie, YUAN Shouqi, PEI Ji,etal. Numerical analysis of the clocking effect on the pressure fluctuation in the centrifugal pump with vaned diffuser[J]. Journal of Mechanical Engineering, 2015,51(4):185-192.

[24]胡勁松．面向旋轉(zhuǎn)機(jī)械故障診斷的經(jīng)驗(yàn)?zāi)B(tài)分解時頻分析方法及實(shí)驗(yàn)研究[D]．杭州：浙江大學(xué)，2003．

[25]胡紅英，馬孝江．局域波近似熵及其在機(jī)械故障診斷中的應(yīng)用[J]．振動與沖擊，2006，25(4)：38-41．

HU Hongying, MA Xiaojiang. Application of local wave approximate entropy in mechanical fault diagnosis[J]. Journal of Vibration and Shock,2006,25(4):38-41.

Numerical Analysis on Pressure Fluctuation of an Axial Flow Fan with Abnormal Installation Angles of a Single Blade

YEXuemin,DINGXueliang,LIChunxi

(Key Lab of Condition Monitoring and Control for Power Plant Equipment of Ministry of Education,North China Electric Power University, Baoding 071003, Hebei Province, China)

Abstract：Taking the OB-84 variable-pitch axial-flow fan as an object of study, an unsteady calculation was conducted on the internal dynamics and fan performance under different deviation degrees of a single blade, so as to obtain the time-domain and time-frequency characteristics, and the distribution of time-averaged static pressure and equivalent stress. Simulation results show that both the pressure signal and pressure fluctuation energy at each monitoring point exhibit unusual features with abnormal installation angles of blade, especially the most obvious changes at impeller outlet. The impacts of abnormal blade installation on time-domain and time-frequency characteristics are enhanced with rising deviation degrees, resulting in increased intensity of pressure fluctuations. The approximate entropy at each monitoring point changes obviously with large deviation degrees, indicating good characterization ability on the abnormal blade installation. Remarkable changes are found in the distribution and values of time-averaged static pressure, but the equivalent stress only varies obviously in values.

Key words：variable-pitch axial-flow fan; deviation degree; abnormal installation angle; pressure fluctuation

收稿日期：2015-07-13

基金項(xiàng)目：河北省自然科學(xué)基金資助項(xiàng)目(E2012502016)；中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)基金資助項(xiàng)目(13MS98)

作者簡介：葉學(xué)民(1973-)，男，河北邢臺人，教授，博士，主要從事流體機(jī)械、流體動力學(xué)理論及應(yīng)用和新能源技術(shù)等方面的研究.

文章編號：1674-7607(2016)05-0395-09中圖分類號：TK223.26

文獻(xiàn)標(biāo)志碼：A學(xué)科分類號：470.30

電話(Tel.)：13932203443;E-mail：yexuemin@163.com.